CC BY
7ЧИСЛЕННАЯ МОДЕЛЬ ЧАСТИЦЫ ПРИМЕНИТЕЛЬНО К ИССЛЕДОВАНИЮ ЭНЕРГЕТИЧЕСКИХ И ПОЛЯРИЗАЦИОННЫХ ХАРАКТЕРИСТИК ИЗЛУЧЕНИЯ, ПРОШЕДШЕГО ЧЕРЕЗ АНИЗОТРОПНУЮ АТМОСФЕРНУЮ СРЕДУ
О. В. Шефер
Институт кибернетики Национального исследовательского Томского политехнического университета, 634050, Томск, Россия
УДК 535.36
Разработана оптическая модель частицы для исследования энергетических и поляризационных характеристик излучения, прошедшего через полупрозрачную пластинку. В рамках метода физической оптики получены соотношения, определяющие элементы матриц рассеяния и экстинкции. Представленные выражения связывают физические параметры частицы (геометрические размеры, показатель преломления) и параметры падающего излучения (длину волны, состояние поляризации) с характеристиками излучения, прошедшего через кристалл, при различных положениях источника, приемника и кристалла. Полученные выражения позволяют рассчитать характеристики трансформированного оптического излучения с учетом возможных колебаний частиц относительно их устойчивого положения в пространстве.
Ключевые слова: физическая оптика, ослабление, рассеяние, поляризация, ориентированные кристаллы.
An optical model of particle is developed for study of the energy and polarization characteristics of radiation, transmitted through an individual, semitransparent plate. In the frameworks of the physical optics method, we obtained the relations determining the elements of the scattering and extinction matrices. The presented expressions relate the particle physical parameters (geometric sizes and refractive index) and the parameters of the incident radiation (wavelength and polarization state) with the characteristics of radiation, transmitted through the crystal for different positions of the source, receiver, and crystal. The obtained relations make it possible to calculate the characteristics of transformed optical radiation, taking into account the possible oscillations of particles relative to their stable position in space.
Key words: physical optics, extinction, scattering, polarization, oriented crystals.
Введение. Кристаллические облака оказывают существенное влияние на радиационные процессы, происходящие в атмосфере. Поэтому в настоящее время проблеме изучения оптических характеристик кристаллических облаков уделяется большое внимание [1—3]. Эта проблема решается в рамках различных международных научно-исследовательских проектов, связанных
Работа выполнена в рамках государственного задания на выполнение НИОКР и при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 12-05-00169а).
с изучением атмосферы Земли. Результатом таких исследований является разработка ряда моделей атмосферы (например, ЕСНАМ, ECMWF, иКМО, LMD). Для дальнейшего повышения точности создаваемых моделей необходимо решить различные проблемы, возникающие, в частности, при определении оптических характеристик для несферических частиц.
Для изучения характеристик излучения, прошедшего через атмосферные образования, прежде всего следует разработать модель отдельной частицы, позволяющую адекватно установить зависимость параметров трансформированного излучения от параметров рассеивателя и падающего излучения. Проблема рассеяния и ослабления света хаотически ориентированными кристаллами рассматривалась в большом количестве работ. Однако в случае ориентированных кристаллов, особенно для ИК-области спектра, она изучена недостаточно.
Одной из важнейших оптических характеристик является ослабление. Для определения ослабления используются два способа. Первый способ основан на интегрировании вектора Пойн-тинга и законе сохранения энергии: экстинкция равна сумме рассеяния и поглощения. В частности, для отдельной частицы справедливо соотношение
Бвх! = Б аЬ$ + ,
где 8х - сечение ослабления; 8аы - сечение поглощения; - сечение рассеяния. В основе второго способа лежит "гипотетический эксперимент", в результате которого определяется ослабление. Этот способ основан на интерференционном аспекте экстинкции: экстинкция равна интерференции между падающим и рассеянным в направлении вперед излучением [4].
При расчете оптических характеристик для мягких частиц, показатель преломления которых близок к единице, как правило, используется метод аномальной дифракции [1—4]. Полное поле, рассеянное в направлении вперед, формируется в результате когерентного сложения дифракционного и рефракционного полей. При этом фактор ослабления осциллирует в окрестности значения, равного двум. При построении моделей кристаллического облака следует учитывать, что показатель преломления ледяных частиц существенно отличается от единицы, особенно при рассмотрении оптических характеристик в ИК-диапазоне длин волн.
Для крупных полупрозрачных кристаллов поле, рассеянное в направлении вперед, определяется, главным образом, дифракционным полем. В этом случае сечение ослабления равно удвоенной площади геометрической тени частицы в направлении распространения излучения [4, 5]. Это возможно в том случае, если рефракционные пучки затухают внутри кристалла либо на выходе существенно отклоняются от направления вперед. Такой подход обеспечивает удовлетворительную точность, когда расчеты экстинкции проводятся для больших хаотически ориентированных кристаллов. При определении ослабления света ориентированными кристаллами, имеющими плоскопараллельные грани, следует учитывать соизмеримость дифракционного поля и поля рефракционных лучей, проходящих через частицу [6]. В сечении пучка параллельных лучей, вышедших из кристалла, фаза волны постоянна. При этом метод физической оптики позволяет естественным путем получить конечные амплитуды полей в направлении выхода пучков из кристалла. Единый подход к определению дифракционного и рефракционного полей позволяет сложить их когерентно. Следует отметить, что метод геометрической оптики ограничивает возможности корректно сложить эти поля. Численное исследование, проведенное в данной
О""]
- 1
- 2
О
п
20
40
60
80
Рис. 1. Зависимость фактора эффективности ослабления от толщины кристалла при различных значениях показателя поглощения (I2/I1=1, у = 0o, а = 0o, l = 500 мкм, ß = 60o, X = 10,6 мкм, n = 1,31):
работе в рамках метода физической оптики, показало, что для ориентированных кристаллов с плоскопараллельными гранями фактор эффективности ослабления может быть отличен от двух [6, 7].
На рис. 1 представлены зависимости фактора ослабления для прямоугольного (кривые 1-3) и гексагонального (кривые 4-6) столбиков от их толщины d (угол у задает ориентацию плоскости поляризации, угол Р - ориентацию частицы относительно направления падения излучения, угол а - поворот кристалла вокруг своей оси; I - длина столбика; п, % - реальная и мнимая части показателя преломления частицы соответственно). На кристалл падает линейно поляризованное излучение (отношение параметров вектора Стокса падающего излучения 71/72=1). Из рис. 1 следует, что для прямоугольного столбика с малой толщиной (так же, как для пластинки) фактор ослабления может изменяться в диапазоне от 0 до 4. При увеличении количества пар плоскопараллельных граней частицы интервал изменения этой характеристики сужается (для гексагонального столбика - от 1 до 3). Заметим, что чем выше значение мнимой части показателя преломления частицы, тем ближе величины ослабления к их предельному значению. В свою очередь нарушение "параллельности" граней также обеспечивает стремление фактора ослабления к асимптотическому значению, равному двум.
На рис. 2 показаны значения фактора эффективности ослабления, вычисленные для пластинчатого кристалла со скошенной гранью. Угол В определяет наклон верхней грани пластинки относительно ее нижней грани. С увеличением угла В между гранями амплитуды осцилля-ций быстро уменьшаются. Нарушение параллельности граней более чем на 10о можно считать достаточным для того, чтобы не учитывать влияние толщи кристалла на характер ослабления оптического излучения. Если атмосферный кристалл не имеет плоскопараллельных граней, то,
1-3 - для прямоугольного столбика (1 - %= 10 4, 2 - %= 10 2, 3 - %=10 1), 4-6 - для гексагонального столбика (4 - %= 10"4, 5 - %= 10"2, 6 - %=10-1)
Рис. 2. Зависимость фактора эффективности ослабления от толщины пластинки со скошенной гранью (/2//1 = 1, у = 0о, а = 250 мкм, р = 0о, Х= 10,6 мкм, п = 1,26): 1 - 0 = 0о, 2 - 0 = 2о, 3 - 0 = 4о, 4 - 0 = 7о, 5 - 0 = 10о как правило, эта "непараллельность" превышает 10°. Таким образом, среди всех крупных атмосферных частиц полупрозрачные ориентированные пластинки выделяются тем, что для кристаллов такой формы фактор эффективности ослабления может принимать значения в наибольшем интервале от 0 до 4 [6, 7].
На рис. 3 показаны зависимости (0) при различных значениях длины волны и показателя преломления. Видно, что при учете рефракционного поля, когда ослабление определяется с заданной точностью, большей длине волны соответствует большее предельное значение 0. На ход кривой ОСХ (0) существенно влияет значение показателя преломления.
а вх1 \
3- \ \ \ \ \ \
V Л/ \ ^ —., _______ < * _ _ - -_ *
Z 1- ■•'А ' N ' / ^ / / / ^^-- - 1 ---2 3 • - 4
0 1 | 1 | 1 2 4 6 8 10 12 1 | 1 14 0,град
Рис. 3. Влияние на зависимость Оесхг* длины волны падающего излучения
и показателя преломления кристалла со скошенной гранью (с1 = 20 мкм, а = 250 мкм): 1 - X = 0,55 мкм, п = 1,31, % = 10-4, 2 - X = 10,6 мкм, п = 1,31, % = 10-4 , 3 - X = 10,6 мкм, п = 1,11, % = 10-4 , 4 - X = 10,6 мкм, п = 1,51, % = 10-4
После прохождения через ориентированную частицу могут измениться не только энергетические, но и поляризационные характеристики падающего излучения. Для полного представления этих характеристик, как правило, рассматриваются матрица рассеяния и матрица экстинк-ции. В данной работе с использованием метода физической оптики строится оптическая модель пластинчатого кристалла для расчета указанных матриц.
Постановка задачи. В качестве рассеивателя рассмотрим гексагональную пластинку толщиной d с длиной стороны ее основания а. Полупрозрачная ледяная пластинка имеет комплексный показатель преломления П - п +1% . Под углом р к основанию пластинки падает плоская волна эллиптической поляризации.
Определим оптические характеристики поля, рассеянного в передней полусфере. Передней будем называть полусферу, содержащую основание пластинки и лучи, рассеянные в направлении вперед. Для этого предлагается следующая схема решения поставленной задачи. Сначала получим полное рассеянное поле как когерентное сложение рефракционного и дифракционного полей. Далее для определения трансформированного поля в дальней зоне используем метод физической оптики и учтем все необходимые максимумы в расчетах интеграла Фраунгофера. Для характеристик излучения, рассеянного в передней полусфере, получим соотношения для сечений рассеяния, приведенных к соответствующим параметрам вектора Стокса, а также представим матрицу экстинкции.
Взаимосвязь положений источника, приемника и частицы для определения характеристик рассеяния излучения в передней полусфере. Применительно к определению характеристик рассеяния оптического излучения в передней полусфере представим необходимые угловые характеристики, связывающие положения источника, приемника и частицы (рис. 4). Пусть в точке О1 находится источник излучения, в точке О2 - приемник, в точке О3 - частица; Оху2 -абсолютная система координат, относительно которой вводятся еще три системы координат О1х1у121. О2х2_у2^2, О3х3у3£3, связанные соответственно с источником, приемником и рассеивате-лем. Координатная плоскость Оху параллельна поверхности земли, а вдоль оси 02 направлена нормаль к ней. Направление распространения падающего излучения совпадает с положительным направлением оси О^. Волновой вектор Пойнтинга к указывает направление падения света на основание пластины. Очевидно, что к || О^. Электрические составляющие падающей волны эллиптической поляризации (Е1, Е2) располагаются соответственно вдоль осей О1х1 и О1_у1, п± - нормаль к плоскости падения. В качестве рассеивающей частицы представлена пластинка.
Плоскость, проходящая через основание кристалла, есть координатная плоскость О3х3у3 ортогональной системы координат О3х3у323, Р - угол между направлением падения излучения к (или осью О^) и нормалью п к основанию пластинки (или осью О323). Электромагнитное поле, представленное лучами, вышедшими после ряда внутренних переотражений, формируется в направлении О323. Вектор к/ указывает направление выхода из пластинки рефракционных лучей. Очевидно, что к || к/1| Угол высоты у и угол широты задают возможные положения пластинки, обусловленные ее колебаниями относительно оси О323. При фиксированном значении у и при непрерывном изменении от 0 до 2п нормаль к основанию пластинки опишет ко-
Рис. 4. Совокупность систем координат для описания рассеяния света в передней полусфере
нус с осью О3^3. Прием рассеянного излучения осуществляется в направлении к^ (или оси О2^2), а ось О2у2 параллельна горизонтальной плоскости (или плоскости земли). Величины Ег, Еп -
составляющие поля, регистрируемого на приемнике, расположенные вдоль осей О2.х2 и О2у2 соответственно. Отклонение направления приема от линии выхода лучей из пластинки (т. е. угол между направлениями Ои к^) обозначим через О.
Для представления нормированных характеристик светорассеяния достаточно определить угловое положение единичных векторов, задающих составляющие рассеянного поля. В связи с этим центры всех четырех систем координат совместим в одной точке О и определим угловые зависимости ортов (х,у,г) абсолютной системы координат с ортами (х{,у{,г{) соотношением
V z J
= S
^хЛ
y»
v z. J
1=1,2,3,
где
S =
cos ( cos fy - sin ( sin ( cos fy cos ( - sin fy 0
cos ( sin fy л sin ( sin fy cosfy
Очевидно, что углы определяют положение базисных векторов х г, уг, г г ( г = 1, 2, 3) каждой
из трех соответствующих систем координат 0ху121 (г = 1, 2, 3) относительно абсолютной системы 0ху2.
Вращение тела в пространстве относительно декартовой системы координат принято описывать матрицей Эйлера. В связи с этим определим положение составляющих падающего поля Е1, Е2 в системе координат, связанной с рассеивателем. Орты систем координат 0х1у1г1 и 0х3у3г3 можно представить с помощью линейного преобразования
(х 1 А3 ( X11
У 3 = A У1
V z 3 J V z1 J
где
A =
cosacos fícosy-sinasiny - cosacos físiny-sinacosy cosasin fí
sinacos fícosy + cosasiny - sinacos físin y + cosacosy sinasin fí - sin fícosy sin físiny cos fí
Углы Эйлера a, fí, y определяются как некоторые комбинации углов 9Х(Х и fy3,(3.
Элементы матрицы A определяют положение векторов k, Ei, E2 в системе координат Ox3y3z3. Следует отметить, что при произвольном угле ориентации плоскости поляризации y составляющие E1, E2 не лежат в плоскости падения волны. Поэтому для дальнейших расчетов характеристик светорассеяния и применения к ним формул Френеля необходимо преобразовать составляющие E1, E2, так чтобы одна из них была перпендикулярна плоскости падения, а другая находилась в ней. Для этой цели используем линейное преобразование
( V ( E11
Ex = F e2
Vk j V k J
где
F =
- cosy sin y 0Л sin y cosy 0 0 0 1
Векторы Е± и Е,, в системе координат, связанной с пластинкой, определяются элементами
первого и второго столбцов матрицы ЛР = АР . Определим углы, задающие положение г2 отно-
' cos^ sin^ 0 - sin^ cosф 0 0 0 1
сительно векторов E^ E± ,k f . Если матрицу В определить как S3S, то искомые проекции вектора z2 можно записать в виде
cos Wx = Af11 B13 + AF21 B23 + AF31 B33 , cos Wy = Af12 B13 + Af22 B23 + Af32B33 , cos Wz = Af13 B13 + AF23 B 23 + AF33 B33 •
Учитывая, что
coswx — sin^cos^, coswy = sin^sin^, coswz — cos0, нетрудно определить 0 и ф. Углы 0 и ф отсчитываются от направления kf. Введем новую систему координат Oxsyszs, связанную с Ox3y3z3:
Aj — Ap S P .
Здесь
r cos ф cos 0 - sinф cosфsin0л St — sinфcos0 cosф sinфsin0 , P — - sin0 0 cos0 y
Электрические составляющие рассеянного поля ESj и E определяются элементами первого и второго столбцов матрицы AT, причем эти векторы лежат в плоскости, перпендикулярной направлению приема Oz2. Напомним, что Ер Еп - составляющие поля, регистрируемого на приемнике. Тогда взаимосвязь между Es, E и Et, Enопределяется с помощью линейных соотношений
Et — - cos u ES + sin u E^ , EII — sin u ES + cos u Es ,
где
cos u — Aj Bj2 ^ Aj B22 ^ Aj B32, sin u — Aj B12 ^ Aj B22 ^ Aj B32.
После падения излучения на пластину часть его отражается в направлении л—Р, а часть проходит внутрь частицы, претерпевая преломление и поглощение. После ряда внутренних переотражений вышедшие из пластинки рефракционные лучи формируются главным образом в направлении р.
Характеристики поля, рассеянного в передней полусфере. Электромагнитное поле падающей плоской волны в системе координат Ox1 y1 z1 можно представить в виде электрической составляющей E — x1E1 + y2E2, магнитной составляющей H — -x1H2 + y1H1 и волнового вектора k. Заметим, что амплитуды электрических и магнитных составляющих связаны между собой, поэтому выражения для магнитных составляющих в данной работе не приводятся.
Полное рассеянное поле образуется при прохождении первичной волны через частицу. В направлении распространения падающего излучения оно является суперпозицией первичной и рассеянной волн. Для ориентированных кристаллов следует учитывать, что поле, рассеянное в направлении вперед, определяется как когерентное сложение двух соизмеримых по величине полей: дифракционного и рефракционного. Это объясняется тем, что, с одной стороны, частица
ослабляет волну за счет возмущения ее фронта, обусловленного образованием в нем геометрической тени, а с другой - изъятая часть фронта после прохождения через частицу при определенных фазовых соотношениях добавляется к возмущенному фронту волны в противофазе, еще более ослабляя его.
Вектор электрической составляющей полного рассеянного поля зададим в виде
Е = Е - Е
где Е. и Е, - векторы падающего и рассеянного полей, которые представляются известными соотношениями [4]
е<кг
Ег - Е , Е, - А-. (1)
гкг
В свою очередь комплексные амплитуды соответствующих полей А - А1 + А2, А1 ± А2, ( Е - Е1 + Е2, Е1 ± Е2) могут быть представлены в виде
Ai 1 (Sii Si
Г К Л f Q Ci21( E^
V A2 J VS21 S22 JVE2 J
(2)
где S¡j - элементы амплитудной матрицы размером 2x2.
Для определения рассеянного поля используем метод физической оптики, который позволяет когерентно сложить дифракционное поле и поле рефракционных пучков с учетом фазовых набегов [6]:
A1 = (AD2 - AR2 ) (fy0 cos ( - Ф0 sin ( A2 = (AD2 - AR2 ) (fy0 sin ( + Ф0 cos (). (3)
Отметим, что орты сферической системы координат Ф0, ф0, r0 представлены элементами матрицы AT, а углы д и <р определены в системе координат Oxiy^, связанной с падающей волной. Согласно принципу Бабине амплитуды дифракционного поля AD¡ и A^ определяются по формулам
Ad = (1 + cos fy) Ei F (fy, () e¥, Ad2 = |- (1 + cosfy)E2F(fy, ()ey, (4)
1 4k 4k
где k = 2k / Л - волновое число; Л - длина волны; у - фазовый набег невозмущенной части
фронта волны после ее прохождения через объем, содержащий пластинку. Угловая функция
F (fy, () является интегралом Фраунгофера
F(fy, () = jjjexp(-ikx cos ( sin fy - iky sin ( sin fy)dxdy, (5)
S
где S - площадь геометрической тени частицы.
Амплитуды рассеянного поля AR¡ и A^ , прошедшего через пластинку, определяются как
сумма рассеянных пучков, образующихся при многократных прохождениях через пластину части фронта первичной волны. Для расчета AR¡ и Ar используем следующие соотношения:
k2 k2
Ar =-(1 + cosfy)(E1a11 + E2a12), AR =-(1 + cosfy)(E2a21 + E2 a22),
1 4k 2 4K
a11 = cos ucosyBy + sin usinyBL, a12 — cos ucosyBL- cos usinyB^, (6)
a21 — cosu siny Bl - sinu cosy By, a22 = sinu siny By + cosu cosy Bl. Угловые функции By и Bl определяются как комбинации коэффициентов Френеля TIT,R^Rl и угловых функций рассеяния j-x пучков F.(3,ф) eWi:
By(&,p) = Т^R^Fj(&,ф) e¥i, B±(&,ф) _ TITlXR^Fj(ЗД e*' . (7)
j=i м
Коэффициенты прохождения TpTl ,T^,T± и отражения Ry, R± определяются известными соотношениями
Т _ 2П cos & _ 2ncos3 ~ _ IcosP
Т _ ~ , _ ~ , Т _ ■
cosí!? + ñcosP ñ cos3 + cosP ñcosP + cos3
(8)
~ 2cosP cos3 - ñcos P ñ cos 3 - cos P —-- , R —---, —---
cos P + ñ cos 3 cos 3 + ncos P ñ cos 3 + cos P
Комплексный угол преломления 3 определяется из закона Снеллиуса sin 3 — sin 3/ñ , где ñ -комплексный показатель преломления частицы. Сомножитель F. (3,ф) в (7) представляет собой интеграл Фраунгофера
F.(3,ф) — j"jexp(-i'kxcos^sin3-ikysinфsin3)dxdy , (9)
где Sj - сечение пучка, выходящего из пластины после 2j -1 прохождений через нее в направлении |3. Фазовые набеги пучков различной кратности отражения обозначим через у... В случае
Ф=0 выражения (5) и (9) существенно упрощаются, при этом F(3, ф) — S, F . (3, ф) — S..
Определим энергетические и поляризационные свойства излучения, рассеянного в передней полусфере. Для этого рассмотрим сечения рассеяния, пропорциональные соответствующим параметрам вектора Стокса:
а — ^ I. (10)
1i
(I1 - интенсивность электромагнитного поля падающей волны). Параметры вектора Стокса If
выражаются через амплитуды трансформированного поля:
If1 — |Ej|2 + |En|2, If2 — |E,|2 -|En|2, f — 2Re(EIEII), Iu — 2Im(EIEII). (11)
В данной работе сечениями рассеяния будем называть все четыре величины af , несмотря
на то что классическим сечением рассеяния является только . Учитывая соотношения
(3)—(11) и проводя необходимые алгебраические преобразования, для сечений рассеяния получаем выражения
^ _ Му Ь-, у=1, 2, 3, 4,
1=\ А
к2 (1 + 008^) 2
IV _— --- ,
л 2
I (г = 1, 2, 3, 4) - параметры вектора Стокса падающего излучения,
м _ |Ьп| + |^121 + Ы + |Ь221 м _ |Ьп| ~ |Ь121 + |Ь211 ~\Ь.
11 2 , 12 2
М13 _ + Ь21Ь2*2 К М14 _- 1т(Ь11Ь1*2 + Ь21Ь22 ) ,
2 |2 |2 |2 |2
__/22 |
2 ' ~~22 2
М _ Ы + |Ь121 ~ Ы ~ |Ь221 м _ |Ьп| ~ Ы ~ |Ь211 + |Ь: М 21 _ ^ , М 22 _
М23 _ МЬ„Ь1*2 - Ь21Ь22 ) , М24 _-1т(Ь11Ь1*2 ^^ ,
М31 _ ^11^1 + Ь21Ь2*2) , М32 _ ^11^1 - ¿Л) ,
М33 _ Яе^Ь* + Ь^Ь*), М34 _- МЬцЬ; - ¿12^2*1) ,
М41 _ 1т(ЬпЬ2*1 + Ь12Ь22 ) , М42 _ 1т(Ь11Ь2*1 - ^^ ,
М43 _ 1т(ЬпЬ2*2 + Ь12Ь2*1), М44 _ ММм - Ь12Ь*) •
Элементы матрицы рассеяния Му представляют собой комбинации элементов Ьк/ (к=1, 2; /=1, 2) (или элементов амплитудной матрицы (2)). Следует отметить, что величины Ькг и Ski различаются только коэффициентом. Элементы Ьк/ определяются как разность соответствующих составляющих амплитуд дифракционного и рефракционного полей (3).
После прохождения падающего излучения через частицу с фиксированной ориентацией в пространстве могут измениться не только его энергетические, но и поляризационные характеристики. Для полного описания процесса ослабления введем так называемую (4х4)-матрицу экстинкции К. Физический смысл матрицы экстинкции состоит в следующем. Рассмотрим пучок параллельных лучей, распространяющихся в направлении к через среду, содержащую частицу. Изменение вектор-параметра светового пучка 1(к) вдоль части пути ds описывается уравнением [1]
^^ _- К (к )1(к ). ds
Матрица К может быть представлена в виде суммы двух когерентных слагаемых. Одно из них содержит характеристики трансформированного излучения, для которого различные компоненты поляризованного падающего поля после прохождения через среду по-разному поглощаются, другое характеризует ситуацию, когда указанные компоненты по-разному преломляются, обеспечивая анизотропность среды [8].
В работах [1-3, 8] выражения, описывающие матрицу экстинкции (МЭ), получены на основе оптической теоремы с использованием вектор-параметрического описания. В этих работах изложен вывод МЭ (а в [2] и ее разновидности - амплитудная матрица экстинкции, матрица
К = ^ к2
(12)
экстинкции Стокса) в операторной форме. Такой подход определяет связь МЭ с параметрами вектора Стокса и амплитудной матрицей. На заключительном этапе вывода каждый элемент матрицы экстинкции представляет собой линейную комбинацию элементов амплитудной матрицы. Матрица экстинкции, полученная на основе оптической теоремы, может быть представлена в следующей форме:
'Яе(Яц + ¿22 ) Яе(¿11 - ¿21 ) Ке(¿12 + ¿21) -1т(¿12 - ¿21 ) Яе(Яц - ¿22) Яе(¿ц + ¿22) Яе(¿12 - ¿21) - 1т(¿12 + ¿21) Яе(¿12 + ¿21) -Яе(¿12 - ¿21) Яе(( + ¿22) 1т(( - ¿22) ^ 1т(¿12 - ¿21) 1т(¿12 + ¿21) - 1т(( - ¿12) Яе(¿п + ¿22) Из (12) следует, что для элементов матрицы К справедливы соотношения
К11 _ К22 _ К33 _ К44 , К12 _ К21 , К13 _ К31, К14 _ -К41 , К32 _ -К23 , К24 _ -К42 , К34 _ -К43 .
Для отдельной частицы сечение экстинкции определяется через элементы первой строки матрицы (12) [2]:
¿ехг = К11 + К12 ~Т + К13 "Г + К14
11 А 11
Для представления сечения ослабления в случае поляризованных полей используется известная формула [4]
4 п Яе {ЕА| ¿.с} к2 Е*Е
Один из основных выводов, который следует из анализа (13), состоит в том, что излучение различной поляризации ослабляется по-разному [6]. Заметим, что выражение (13) представляет собой одну из частных форм оптической теоремы. Общность выражений (12), (13) становится более очевидной, если в первом из них в качестве множителя ввести величину 4п/к2 . В ряде представлений матрицы экстинкции этот коэффициент опускается, но он, очевидно, учитывается при определении ¿¡^. В представлении (12) элементы матрицы экстинкции имеют размерность площади (так же как и величина ¿ех1 (13)). В данном случае в выражение для элементов амплитудной матрицы ¿1] в качестве сомножителей входят квадрат волнового числа и дифракционный интеграл, имеющий размерность площади [4]. Следует отметить, что для пластинок в случае рассеяния в направлении вперед (т. е. при 9=0°) для элементов амплитудной матрицы справедливо соотношение ¿12= ¿21 (так же как для а12= а21 и Ь12= Ь21), а значит, элементы матрицы экстинкции, определяемые как разность ¿12 и ¿21, обращаются в нуль.
Термин "сечение экстинкции" является классическим. Этим и объясняется выбор названия матрицы (12) как "матрицы экстинкции". Следует отметить, что существуют и другие названия матрицы (12). Различаются также формы ее представления [1, 2, 8] , главным образом, в зависимости от исходного описания рассеянного поля (см. соотношения (1)). Однако при расчетах значений элементов МЭ получается один и тот же результат, поскольку одна форма представления рассеянного поля сводится к другой. Также следует отметить, что элементы амплитудной
^^ . (13)
матрицы и дифракционный интеграл могут рассматриваться как безразмерные величины либо иметь размерность длины или площади. В этом случае указанные соответствующие характеристики различаются, как правило, только множителем.
Заключение. С использованием метода физической оптики для исследования характеристик электромагнитного излучения, прошедшего через атмосферные ледяные образования, разработана численная модель. В качестве отдельной частицы рассмотрена полупрозрачная пластинка, для которой свойствен наибольший интервал возможных значений фактора ослабления 0^4. С помощью метода физической оптики получены соотношения для элементов матриц рассеяния и экстинкции с учетом возможных колебаний частиц относительно их устойчивого положения. Эти выражения позволяют рассчитать как энергетические, так и поляризационные характеристики светорассеяния в зависимости от размеров частиц и показателя преломления при различных положениях источника, приемника и рассеивателя, а также при любом значении длины волны из оптического диапазона и любом состоянии поляризации падающего излучения. Представленная численная модель позволяет исследовать сложную структуру зависимости рассеяния вблизи направления вперед.
Список литературы
1. Light scattering by nonspherical particles. Theory, measurements, and application / Ed. by M. I. Mishchenko, J. W. Hovenier, L. D. Travis. California: Acad. Press, 2000.
2. Scattering, absorption, and emission of light by small particles / Ed. by M. I. Mishchenko, L. D. Travis, A. A. Lacis. Cambridge: NASA, 2002.
3. Light scattering reviews: single and multiple light scattering / Ed. by A. Kokhanovsky. Berlin: Springer, 2006.
4. Борен К. Поглощение и рассеяние света малыми частицами / К. Борен, Д. Хафмен. M.: Мир, 1986.
5. Волковицкий О. А. Оптические свойства кристаллических облаков / О. А. Волковицкий, Л. Н. Павлова, А. Г. Петрушин. Л.: Гидрометеоиздат, 1984.
6. Popov A. A. Characteristics of extinction for a hexagonal crystal // SPIE. 1998. V. 3220. P. 380-389.
7. Popov A. A., Shefer O. V. Estimate of the extinction of optical radiation by crystals lacking planeparallel faces // Atm. Opt. 1990. V. 3, N 5. P. 416-421.
8. Розенберг Г. В. Вектор-параметр Стокса // Успехи физ. наук. 1955. Т. 56, вып. 1. C. 78-110.
Шефер Ольга Владимировна - канд. физ.-мат. наук, доц. Института кибернетики Томского политехнического университета; тел.: (382-2) 42-07-27; e-mail: [email protected]
Дата поступления - 02.10.12 г.
