О вычислении амплитуды рассеяния включения в сложном поле Текст научной статьи по специальности «Физика»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2009, том 19, № 2, c. 85-89

-МАТЕМАТИЧЕСКИЕ МОДЕЛИ =

УДК 534.29; 534.138 © Б. П. Шарфарец

О ВЫЧИСЛЕНИИ АМПЛИТУДЫ РАССЕЯНИЯ ВКЛЮЧЕНИЯ

В СЛОЖНОМ ПОЛЕ

В работе показано, что выражение для амплитуды рассеяния включения в поле, состоящем из совокупности плоских волн, удовлетворяет однородному уравнению Гельмгольца вне источника этих волн по переменным координатам включения. Это позволило получить простые асимптотические ряды по обратным степеням расстояния между источником волн и включением для вычисления амплитуды рассеяния. Приведены примеры применения этих рядов.

Кл. сл.: уравнение Гельмгольца, амплитуда рассеяния, асимптотический метод

ВВЕДЕНИЕ

В многочисленных приложениях теории волнового и квантовомеханического рассеяния необходимо уметь вычислять амплитуду рассеяния (ар) включения в сложном первичном поле, отличном от поля плоской бегущей волны. В линейной теории рассеяния решение этой задачи в конечном итоге сводится к вычислению в общем случае несобственного кратного интеграла Фурье (см., например, [1, 2]). Обычно для этого приходится пользоваться асимптотическими методами (перевала, стационарной фазы). В этом случае, как правило, обходятся только старшим членом разложения по обратным степеням большого параметра ввиду сложности вычисления последующих членов. В настоящей работе будет предложен эффективный метод вычисления амплитуды рассеяния произвольного включения в произвольном первичном поле, основанный на том, что она по параметрам местоположения включения удовлетворяет однородному уравнению Гельмгольца, что позволяет воспользоваться геометрооптическим рядом, все старшие члены которого легко вычисляются рекурсивно по известному первому члену. Проблема будет решена применительно к акустическому случаю ввиду его математической простоты и прозрачности.

ПОСТАНОВКА ПРОБЛЕМЫ

Пусть идеальное жидкое пространство, в котором находится включение, полагается однородным с плотностью р и скоростью с . В пространстве подразумевается наличие суммы падающей и (х) и рассеянной от включения и!, (х) гармонических волн с фактором е~,ш, который далее опускается:

u (x) = u (x) + us (x). (1)

Включение расположено в области D, внутри которой находится начало координат. Рассеянное поле в волновой зоне, как известно, описывается выражением

Us (x) = um (xx) —¡—j— + o(|x|-1), (2) lxl

где k = a/ c — волновое число; x = (x, y, z) = = (R, в, p) — точка наблюдения соответственно в

декартовых и сферических координатах; x

x = ^-т = (sinecosp,sinesinp,cose) — вектор на lxl

единичной сфере, направленный в точку наблюдения; uM (x) — амплитуда рассеяния включения, характеризующая угловое распределение поля рассеяния включения при его облучении полем u (x). Обычно выделяется случай плоской падающей волны u (x) = eM x. Здесь d = (sin a cos в, sin a sin в, cos a) — единичный вектор, совпадающий по направлению с волновым вектором плоской волны. В этом случае фиксируется зависимость рассеянной волны us (x) = = us (x, d) от направления распространения падающей плоской волны единичной амплитуды и нулевой фазы и, следовательно, аналогичная зависимость ар um (x) = um (x, d). Очевидно, что вектор x в функции uM (x, d) привязан к точке расположения включения. В силу линейности задачи очевидно, что в случае, когда падающая волна представляет собой совокупность плоских волн со спектром g (d)

(x) = J g (d)ekdx ds(d),

(3)

ар вычисляется следующим образом (см., например, [2]):

,(x) = J g (dK (X, d)d s(d).

ik 2n

П

2n ~~

J d^ J M(d)e

ikR(sinasinecos(p—p)+cosacose)

0 0 x sin a da.

начало координат. Пусть включение находится в точке х = х0 = (х0,Уо,г0) = . Тогда, со-

гласно (5), включение облучают плоские волны с амплитудой

(4)

g (d) = ^u(d)ea 2п

(6)

Здесь О — область определения спектра падающей волны; = 8тас1ас1в — элемент поверхности единичной сферы.

Известно [3], что если функция иппс (х) является решением однородного уравнения Гельмгольца во всем пространстве х е Я3 (так называемое целое решение), то в ее представлении по плоским волнам (3) участвуют только однородные плоские волны и интегрирование в (3), (4) осуществляется по поверхности единичной сферы О = О (область видимости). Если же иппс(х) является решением однородного уравнения Гельмгольца только в бесконечной подобласти Я3 и удовлетворяет условиям излучения Зоммерфельда, то интегрирование в (3), (4) происходит как в области видимости (однородные плоские волны), так и вне ее (неоднородные плоские волны).

РЕШЕНИЕ ЗАДАЧИ

Итак, для расчета результирующей амплитуды рассеяния необходимо вычислять интеграл (4). Однако вычисление этого интеграла может оказаться довольно сложной задачей. Ниже приводится достаточно простой алгоритм вычисления амплитуды рассеяния с заданной точностью.

При построении алгоритма будем исходить из самого сложного случая: первичное падающее поле представляет собой излученное решение.

Воспользуемся выражением для излученного решения в виде разложения по плоским волнам, приведенного в [3, 4, 5], для первичного падающего поля и ппс (х)

( ) i f f u(kx, ky) ikd-x d d

u,nc (x) = — J J—-e dkxdky =

—f —f z

С учетом этого запишем (4) в виде интеграла по сферическим углам типа (5):

имЯоА^о) =

2- J J u(a,ви (0,q>,a,P) x

П

Л--If

2п 2

2п

0 0

x e

ikR (sinasine0 cos(fi—p0 )+cosacose0 )

sin a da de.

(7)

Между выражением (7), где включение смещено в точку x = x0, и выражениями (2), (4), где включение находится в начале координат, противоречий нет, т. к. углы в и р в результирующей

амплитуде рассеяния ит привязаны, по определению функции ит (X, d) и согласно (4), и согласно (7), к системе координат с началом в точке расположения включения.

Функция ит(в,р,R0,e0,p0) в (7) зависит от R0, в0, р0 как от параметров, характеризующих координаты включения. Как видно из (7), функция ит(в,р,R0,e0,p0) отвечает всем признакам излученного решения вида (5) по переменным R0, в0 и р0 (излученное решение со спектром U(a, в)их (e,p,a,e)). Это означает, что при фиксированных значениях вир функция ит (в,р, R0,e0,p0) удовлетворяет однородному уравнению Гельмгольца везде, где справедливо представление (5) для первичного поля по плоским волнам. А это в свою очередь серьезно облегчает вычисление интеграла (7). А именно вместо того, чтобы производить прямые вычисления или пользоваться общим, но трудоемким методом перевала, можно воспользоваться изящной теоремой, полученной Аткинсоном (Atkinson F.V.) и Уилкоксом (Wilcox C.H.) [6, с. 84], которая в данном случае записывается так:

(5)

ikR,

Kf (в,р, R0,e0,(0)

<Fn (в,Ро)

R0 n=0 (kR0)n

(8)

Здесь и(кх,ку) — спектр плоских волн в первичном поле. В (5) подразумевается, что источник поля сосредоточен в области В, включающей в себя

где функции ¥п(в0,ф0) определяются из следующего рекуррентного соотношения:

□

f

Fn =

1 д в д )

--(sin0o-)-

sinö„ дв0 дв0

1

sin в0 др0

- n(n -1)

Fn

2in

а Fo(é>o,^o) равно

Fo(eo,Po) = û(90, <Po К (в,Р,во,Ро).

(9)

(1o)

При um(в,(р,60,ф0) = 1 равенство (8) сходится абсолютно и равномерно при R > R ', его можно дифференцировать по Ro, в0, р0 почленно любое число раз и полученные после этого ряды также сходятся абсолютно и равномерно [6]. Здесь R ' — радиус наименьшего шара, включающего в себя область D источника первичной волны. При функции um(в,р,в0,р0), отличной от константы, этот вопрос требует отдельного рассмотрения, и наверняка можно сказать только то, что асимптотический ряд (8) по обратным степеням волнового параметра kR0) определяется рекурсией (9), (10) единственным образом в силу единственности представления асимптотическими рядами и удовлетворения функции um(в,р,Ro,eo,po) однородному уравнению Гельмгольца в указанной области. А следовательно, и в этом случае ряд (8) позволяет вычислить совокупную амплитуду рассеяния um (в,р, Ro, в0, р0 ) с заданной точностью.

В случае, когда первичное поле является целым решением, вычисление интеграла (4), собственного в этом случае, является исходно существенно более простой задачей и может быть проведено точно или приближенно. Остановимся на асимптотике интеграла (4) в случае, когда падающие плоские волны являются только уходящими, т. е. имеют угол падения а е [0, п / 2]. Для этого случая интеграл (7) преобразуется к виду

u« (в,Р Ro,eo,Po) =

п

ik 2п 1

— f fu(a,ß)uM (в,р,а,Р) х 2п

0 0

ikRo (sinasineo cos(ß-p0 )+ cosacoseo )

ет асимптотику

^ ikRo

Гельмгольца, для совокупной амплитуды рассеяния um (в,р, Ro,eo,po) остается справедливым

п

представление (8)—(10). При в0 = 0 или в0 = —

асимптотика (12) для (11) нарушается, и пользоваться представлением (8)-(10) неправомерно.

В случае, когда u(a, ß) из (7) или (11) является суммой дельта-функций (например, случай бегущей или стоячей плоских волн), вычисление результирующей амплитуды рассеяния непосредственно из выражения (7) или (11) становится тривиальным.

Рассмотрим применение (8)—(10) на некоторых примерах.

Пример 1. Случай сферической падающей волны. Рассмотрим пример расчета радиационного давления на сферическую частицу, находящуюся в сферической волне. Как известно (см., например, [7]), амплитуда рассеяния сферической частицы при падении плоской волны в длинноволновом приближении равна сумме монопольного и ди-польного моментов

um (x, d) = Ao + Aj cos у =

= A + A1(sinesinacos(p-ß) + cosecosa). (13)

Здесь у — угол между векторами x и d ; A и A1 — постоянные коэффициенты, зависящие от частоты и свойств включения и среды. Пусть источник сферической волны расположен в начале координат. Центр рассеивателя находится в точке x0 = (0,0, z0), где z0 > 0. Для определения uM (X) воспользуемся выражениями (3), (4). Представление сферической волны при z > 0 по плоским волнам в этом случае имеет вид (см., например, [3,

4]):

(x) =■

x

- п

2п

2п 2

ik(sinacos ßx+sinasinßy + zcos a)

sinadadß. (14)

хe"-^"".................sinadadß. (11)

п

Согласно [4], при 0 <в0 < — интеграл (11) име-

В точке x0 = (0,0, z0) поле равно

u (0,0, zo) = ■

им(6,Ф,КоА,Фо) ~ Fo(0o,Фo)—, (12)

где Е0(в0,р0) определяется из (10). А в силу того, что (11) удовлетворяет однородному уравнению

2 —

— f f eikz°cosa sinadadß = 2п

, ikz0

Сравнивая последнее выражение с (3), находим

2

u

п

x

0

0

g (d) = l!^e*zocosa. ç учетом этого, подставляя (13) 2п

в (4), находим

"ю (X =

П

Л —гад ik 2П 2

2п

i J (

+A1 (sin в sin a cos(p -в) + cos в cos a) ) х

x eikz°cosa sin a da de =

f A.eikz° Л f

v zo

Л

-(kz. + i)

cos в.

Вычислим иад (в) по соотношениям (8)—(10). Имеем: R. = z., в. = 0 , г?(в.,%) = 1,

W(х, У ) =

1

2n°x°y

Дx,y,z) = J J^(kx,ky)e

i( kxX+kyy + kzZ )

dkxdky, z > 0.

ем:

U,ne (X, У, z) =

П

--гад

7 2 2n 2 1 > 2 2 '2

k Г Г —k a sin a

e 2 cosax

4П

2

0 0

0 0

x e

ik ( x sinacos в+y sinasin в+ z cos a)

sin a da de,

(15)

z > 0.

Пусть включение (13) находится в точке x0 = ^0,в0,р0). Тогда имеем:

иад (в,P, R0,в0,^0) =

П

Л —уад

k2 г 2г

= —j J J (A0 + A1(sinвsinacos(^-в) +

4n

+cosвcosa) )É

a sin a

x cosae

ikR0 (sin в) sinacos(e~^0 )+ cos a cos в))

sin a da de. (18)

(16)

Р0{в0,Р0>) = и(в0, Р0 К (в, РА-Ро) = = А0 + А1 СО80СО800,

^(в0,р0) = ,А1 со8всо8в0;

^(в0,р0) - 0, п = 2,3,... (17)

Подставляя (16), (17) в (8) и полагая в0 = 0 , получаем в точности выражение (15).

Пример 2. Случай гауссова пучка. В разобранном выше примере сферической волны тривиально вычислялся интеграл (4) и преимущества предлагаемого алгоритма были не столь очевидны. В рассматриваемом ниже случае гауссова пучка преимущества окажутся существеннее.

Под гауссовым пучком здесь понимается плоская волна в,ь, при г = 0 пропускаемая через маску с функцией прозрачности у/(х, у), совпадающей с гауссовой кривой

Последний интеграл, как видно, уже совсем не прост, поэтому прибегнем к описываемому алгоритму. В силу цилиндрической симметрии падающего поля ограничимся случаем р0 = 0. Тогда, согласно (7) и (10), имеем

Р0(в0,Р0) = к

2 ni

( A0 + A (sin в0 sin в cos p + cos в0 cos в)) x

x e

--k2 a2 sin2 в0 2 0

cos в0.

(19)

Из (9) получаем с точностью до третьей степе-

ни множителей cos в0 и sin в0

Примем далее, что ах = ау = а . Тогда нетрудно показать, что при г > 0 падающее поле определяется выражением (см., например, [8])

Здесь ^(kx, ky ) — фурье-образ функции у/(x, y). Переходя к сферической системе координат, име-

k -2k2 a2 sin2 в0

F^,^ e 2 x

4n

x(2cosв0A0 + A1(4cos2 в0 - 2sin2 в0)cosв +

+A1 (5 cos в0 sin в0 - ctg в0 cos2 в0 ) sin в cos p). (20)

Подставляя (19), (20) в (8), получаем с точностью до O(k-2 R0-3) значение интеграла (18).

Отметим, что в0 = 0 является особой точкой дифференциального оператора в (9), поэтому при вычислении Fn(в0 = 0), n = 2,3,... возможны особенности, которых можно избежать поворотом системы координат.

ВЫВОДЫ

Таким образом, благодаря выявлению того факта, что амплитуда рассеяния удовлетворяет однородному уравнению Гельмгольца по параметрам местоположения включения, предложен ори-

e

гинальный алгоритм расчета результирующей амплитуды рассеяния включения в поле произвольной волны, основанный на представлении решения в виде геометрооптического ряда.

СПИСОК ЛИТЕРАТУРЫ

1. Шарфарец Б.П., Курочкин В.Е., Князьков Н.Н. Радиационное давление в произвольном падающем поле. Связь с амплитудой рассеяния включения // ДАН. 2008. Т. 421, № 2. С. 186189.

2. Зацерковный А.В., Сергеев В.С., Шарфарец Б.П. Использование амплитуды рассеяния для решения задач дифракции волн в полупространстве // Акуст. журн. 2001. Т. 47, № 5. С. 650-656.

3. Devaney A.J., Wolf E. Multipole Expansions and Plane Wave Representations of the Electromagnetic Field // J. Math. Phys. 1974. V. 15, N 2. P. 235-244.

4. Hansen T.B., Yaghjian A.D. Plane-Wave Theory of Time-Domain Fields. New York: IEEE Press, 1999. 367 p.

5. Шарфарец Б.П. О некоторых свойствах амплитуды рассеяния // Научное приборостроение. 2007. Т. 17, № 4. С. 55-60.

6. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М.: Мир, 1987. 311 с. (Colton D., Kress R. Integral equation methods in scattering theory. New York: Wiley&Sons, 1983).

7. Морс Ф.М. Фешбах Г. Методы теоретической физики. Т. 2. М.: Изд-во иностр. лит-ры, 1960. 860 с. (P.M. Morse and H. Feshbach. Methods of Theoretical Physics. Vol. II. New York: McGraw-Hill, 1953).

8. Шарфарец Б.П., Курочкин В.Е., Князьков Н.Н. Радиационное давление на включение с заданной амплитудой рассеяния в произвольном внешнем поле // Акуст. журн. 2009. Т. 55, № 2. С. 147-152.

Институт аналитического приборостроения РАН,

Санкт-Петербург

Материал поступил в редакцию 8.02.2009.

ON CALCULATION OF A SCATTERING AMPLITUDE OF CLOSING IN A COMPLEX FIELD

B. P. Sharfarets

Institute for Analytical Instrumentation RAS, Saint-Petersburg

The works shows that the expression for the closing scattering amplitude in the field, consisting of a set of flat waves satisfies to homogeneous Helmholtz equation out of a source of these waves on variable co-ordinates of closing. It has allowed to receive simple asymptotic rows on return degrees of distance between a source of waves and closing for scattering amplitude calculation. Examples of these rows use are given.

Keywords: Helmholtz equation, scattering amplitude, asymptotic method