Научная статья на тему 'О соотношениях, связывающих поле рассеяния с амплитудой рассеяния'

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

CC BY
91
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Научное приборостроение
ВАК
RSCI
Область наук
Ключевые слова
АМПЛИТУДА РАССЕЯНИЯ / РАЗЛОЖЕНИЕ АТКИНСОНА-УИЛКОКСА / ПРЕДСТАВЛЕНИЕ УИТТЕКЕРА / ПРЕДСТАВЛЕНИЕ ДЕВАНИ-ВОЛЬФА / ПАДАЮЩАЯ ПЛОСКАЯ ВОЛНА / КАНОНИЧЕСКАЯ АМПЛИТУДА РАССЕЯНИЯ / РЕЗУЛЬТИРУЮЩАЯ АМПЛИТУДА РАССЕЯНИЯ / ВОЛНОВАЯ ФУНКЦИЯ ГЕРГЛОТЦА / SCATTERING AMPLITUDE / ATKINSON-WILCOX DECOMPOSITION / WHITTAKER''S REPRESENTATION / DEVANEY-WOLF REPRESENTATION / INCIDENT PLANE WAVE / THE CANONICAL SCATTERING AMPLITUDE / RESULTING SCATTERING AMPLITUDE / HERGLOTZ WAVE FUNCTION

Аннотация научной статьи по физике, автор научной работы — Шарфарец Борис Пинкусович

В работе рассмотрены различные функциональные соотношения между основными параметрами процесса рассеяния: полем рассеяния и амплитудой рассеяния. Это соотношения либо геометрооптического типа (представление в виде ряда Аткинсона-Уилкокса), либо связь через представление обеих функций в ряды по сферическим функциям (мультипольные представления), либо в виде интегрального представления (представление Девани-Вольфа). Такое многообразие представлений возможно вследствие аналитических свойств обеих функций: поле рассеяния является излученным решением уравнения Гельмгольца, любое дважды дифференцируемое решение которого есть аналитическая функция своих аргументов, а также тем, что амплитуда рассеяния является целой аналитической функцией своих аргументов. Приведен также аналог разложения типа Аткинсона-Уилкокса для обобщенной амплитуды рассеяния, что возможно только в случае, когда первичное падающее сложное поле является излученным решением уравнения Гельмгольца. Результирующая амплитуда рассеяния в этом случае также подчиняется уравнению Гельмгольца. Показано, что волновая функция Герглотца с точностью до постоянного множителя совпадает с разложением Уиттекера. Приведенные результаты весьма полезны для приложений, и, в частности, в задачах научного приборостроения.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ON THE RELATIONS CONNECTING THE SCATTERING FIELD WITH THE SCATTERING AMPLITUDE

Various functional relationships between the main parameters of the scattering process are considered in this paper: the scattering field and the scattering amplitude. These relations are either geometrically-optic type (decomposition in the form of the Atkinson-Wilcox series), or a connection through the representation of both functions in series in spherical functions(multipole representations), or as an integral representation (Devaney-Wolf representation). Such a variety of representations is also possible due to the analytic properties of both functions: the scattering field is an radiation solution of the Helmholtz equation, any two differentiable solutions of which are analytic functions of their arguments, and also because the scattering amplitude is an entire analytic function of its arguments. An analog of an Atkinson-Wilcox representation for the generalized scattering amplitude is also given, which is possible only in the case when the primary incident complex field is a radiating solution of the Helmholtz equation. The resulting scattering amplitude in this case also obeys the Helmholtz equation. It is shown that the Herglotz wave function coincides, up to a constant factor, with the Whittaker representation. The above results are very useful for applications, and in particular, in the problems of scientific instrument making.

Текст научной работы на тему «О соотношениях, связывающих поле рассеяния с амплитудой рассеяния»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2017, том 27, № 3, c. 128-136

МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ ^ =

В ПРИБОРОСТРОЕНИИ

УДК 534.874+534.26 © Б. П. Шарфарец

О СООТНОШЕНИЯХ, СВЯЗЫВАЮЩИХ ПОЛЕ РАССЕЯНИЯ С АМПЛИТУДОЙ РАССЕЯНИЯ

В работе рассмотрены различные функциональные соотношения между основными параметрами процесса рассеяния: полем рассеяния и амплитудой рассеяния. Это соотношения либо геометрооптического типа (представление в виде ряда Аткинсона—Уилкокса), либо связь через представление обеих функций в ряды по сферическим функциям (мультипольные представления), либо в виде интегрального представления (представление Девани—Вольфа). Такое многообразие представлений возможно вследствие аналитических свойств обеих функций: поле рассеяния является излученным решением уравнения Гельмгольца, любое дважды дифференцируемое решение которого есть аналитическая функция своих аргументов, а также тем, что амплитуда рассеяния является целой аналитической функцией своих аргументов. Приведен также аналог разложения типа Аткинсона—Уилкокса для обобщенной амплитуды рассеяния, что возможно только в случае, когда первичное падающее сложное поле является излученным решением уравнения Гельмгольца. Результирующая амплитуда рассеяния в этом случае также подчиняется уравнению Гельмгольца. Показано, что волновая функция Герглотца с точностью до постоянного множителя совпадает с разложением Уиттеке-ра. Приведенные результаты весьма полезны для приложений, и, в частности, в задачах научного приборостроения.

Кл. сл.: амплитуда рассеяния, разложение Аткинсона—Уилкокса, представление Уиттекера, представление Девани—Вольфа, падающая плоская волна, каноническая амплитуда рассеяния, результирующая амплитуда рассеяния, волновая функция Герглотца

ВВЕДЕНИЕ

В прямых и обратных задачах классического рассеяния волновых полей, а также квантовомеха-нического рассеяния частиц большое внимание уделяется двум функциям: самому полю рассеяния и его характеристике дальнего поля — амплитуде рассеяния (ар) (scattering amplitude, или far field pattern в английской транскрипции). Несмотря на то что ар является асимптотической характеристикой поля рассеяния в дальней зоне, существует достаточно много соотношений, связывающих строго функционально поле рассеяния вне рассеи-вателя с его ар. Более того, оказалось, что подобные соотношения верны и для различных связанных между собой амплитуд рассеяния. Дело в том, что возникает необходимость вычисления ар включения как в поле плоской волны, так и в сложном первичном поле, отличном от поля плоской бегущей волны (результирующей ар). В работе [1] был предложен алгоритм вычисления результирующей ар на включении при падении на него сложного поля, вытекающий из установленного в ней факта о том, что результирующая ар удовлетворяет однородному уравнению Гельм-гольца по переменным, совпадающим с координатами включения. Это позволило получить простые

асимптотические ряды для вычисления результирующей ар через ар в поле плоских волн. В работах [6-10] полученная в [1] теория была апробирована на различных достаточно нетривиальных примерах рассеяния волн как в теории радиационного давления [6], так и в различных случаях теории многократного рассеяния [7-10] и показала свою эффективность. Решению этой задачи способствовала достаточно полная изученность аналитических свойств амплитуды рассеяния [5], [11, § 3.2.4], [12-15 и др.].

Ряды, полученные в [1], по структуре оказались аналогичными рядам, использующимся при асимптотических оценках собственно волновых полей, удовлетворяющих вне некоторой конечной области однородному уравнению Гельмгольца и условию излучения Зоммерфельда на бесконечности. Эти разложения определяются теоремой Аткинсона—Уилкокса и связывают между собой волновые поля и отвечающие им характеристики дальнего поля, т. е. диаграммы направленности или амплитуды рассеяния (см. [2-4], [5, с. 84]). Теорема Ат-кинсона—Уилкокса является, по-видимому, первым прецедентом, когда обе характеристики (поле рассеяния и его ар) были связаны функционально.

К настоящему времени появился целый цикл работ, посвященных решению обратных задач

рассеяния (см. работу [13] и обширную библиографию в ней). При их решении эффективно используется некоторая функция, называемая волновой функцией Герглотца (ВФГ), которая оказывается тесно связанной с техникой, применявшейся в работе [1]. Однако значительно ранее было предложено разложение Уиттекера [16], [17, с. 247-250], которое фактически совпадает с ВФГ, а в работе [18] было предложено разложение авторов Девани—Вольфа, расширившее область преобразования Уиттекера с поверхности единичной сферы на комплексную область углов сферических координат, связывающее волновое решение уравнения Гельмгольца с его амплитудой рассеяния. Разложение Девани—Вольфа существенно отличается от разложения Аткинсона—Уилкокса тем, что первое связывает волновое поле с амплитудой рассеяния (диаграммой направленности) посредством интегрального преобразования, а второе посредством ряда.

Изучение взаимосвязей всех упомянутых подходов и учет сильных сторон каждого из них позволяет качественно решать важнейшую задачу теории рассеяния — нахождение всевозможных связей между полем рассеяния и амплитудой рассеяния, что весьма актуально в теории рассеяния в целом и для применения к задачам научного приборостроения по примеру работ [6-10] в частности.

Постановка проблемы

Настоящая работа посвящена изучению связей разложений Уиттекера, Девани—Вольфа, аппарата ВФГ и изложенного в [1]. Кроме того, приведен оригинальный вывод связи двух важнейших характеристик процесса рассеяния: поля рассеяния и амплитуды рассеяния через разложение Дева-ни—Вольфа для более общего случая рассеяния, чем это проделано в [18].

В работе для удобства читателя дается краткая информация о свойствах ВФГ, разложений Аткин-сона—Уилкокса, Уиттекера и Девани—Вольфа.

РЕШЕНИЕ ПРОБЛЕМЫ

Общие положения

Рассмотрим акустический случай рассеяния на примере идеального жидкого однородного безграничного пространства ж е ^ \ D , в конечной области D которого присутствует либо проницаемое неоднородное возмущение акустических свойств среды, либо некоторое непроницаемое неоднородное включение. Поэтому далее в соответствии с этим будем различать соответственно поверхностное и объемное рассеяние. Постановка

задачи акустического рассеяния для обоих этих случаев приведена, например, в работе [13, с. 2].

Далее приведем коротко математическую постановку задачи, которая принимается в данной работе за основу. Пусть идеальное жидкое пространство, в котором находится включение, полагается однородным с плотностью р и скоростью c . Рассматривается скалярная функция u(x), характеризующая акустическое поле. В пространстве подразумевается наличие суммы произвольной падающей и (х) и рассеянной от включения

М (х) гармонических волн с фактором е~, который далее опускается

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

(1)

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

us(x) = u™ (x)

i k\x|

e л i-Ч ^ + °( x ),

x

(2)

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

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

x = ^-г = (sm0cos^,sm0sin^,cos0) — вектор

x

на единичной сфере, направленный в точку наблюдения; um (xx) — ар включения, характеризующая угловое распределение поля рассеяния на включении при его облучении произвольным полем u (x). Обычно выделяется канонический

случай плоской падающей волны u (x) = e1 kd x. Здесь d = (sina cos ^,sinasin f3,cosa) — единичный вектор, совпадающий по направлению с волновым вектором плоской волны. В этом случае фиксируется зависимость рассеянной волны us (x) = us (x, d) от направления распространения падающей плоской волны единичной амплитуды и нулевой фазы, и, следовательно, аналогичная зависимость ар um (xx) = um (xx, d), которую назовем для удобства канонической ар. Очевидно, что вектор x в функции um (x, d) привязан к точке расположения включения.

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

uinc (x) = J g (d)e1 kdxds(d),

(3)

ар вычисляется следующим образом (см., например, [16] для излученных решений, [13, с. 56] для целых решений, определение которых приведено ниже):

представимо в виде разложения

uM (x) = { g (dK (x, d)ds(d).

(4)

8u ( x )

8| x|

- iku (x ) = (

,x,

VI I/

u (x )=£ ^m

r n=0 '

Здесь Q — область определения спектра падающей волны; ds(d) = sin a da d^ — элемент поверхности единичной сферы.

Определение: Если u (x) является решением уравнения Гельмгольца

Au (x ) + k 2u ( x ) = 0

в области x e R3 \ D , где D e R3 — конечная область и удовлетворяет условиям излучения Зом-мерфельда

которое справедливо при г > Л0 и сходится абсолютно и равномерно по г , в , р. Этот ряд можно почленно дифференцировать по г, в, р любое число раз, и полученные дифференцированием ряды также сходятся абсолютно и равномерно.

Элементы последнего ряда находятся из рекурсии [5, с. 86]

2iknFn (в,р) =

= В(Fn_l (в,р)) + п(п -1)Fn_l (в,р), п = 1,2,3,...,

D 1 8Г .8) 1 82

где B = ——-—| srn6 — | +

86 / sin2 6 8q>2

опера-

то такое решение называется излученным решением [13, с. 19]. Если u (x) является решением уравнения Гельмгольца во всем пространстве x е R3, то такое решение называется целым решением [13, с. 20].

Известно [16-18], что если функция u(x) является целым решением однородного уравнения Гельмгольца, то в ее представлении по плоским волнам (3) участвуют только однородные плоские волны и интегрирование в (3), (4) осуществляется по поверхности единичной сферы Q = Q (область видимости). Если же u(x) является излученным решением уравнения Гельмгольца, то интегрирование в (3), (4) происходит как в области видимости (однородные плоские волны), так и вне ее (неоднородные плоские волны evanescent wave).

Далее приведем необходимые сведения по всем рассматриваемым разложениям.

Теорема разложимости Аткинсона—Уилкокса

[2-4], [5, с. 84], [15]

Приведем теорему разложимости Аткинсона— Уилкокса в редакции работы [5, с. 84].

Теорема 3.6. Пусть u е C2 (R3 \ D) является

решением однородного уравнения Гельмгольца, удовлетворяющим условию излучения Зоммер-фельда. Пусть R0 таково, что Qr =

= {x е R3,|x| = R0}c R3 \ D, и пусть (r,6,q>) —

сферические координаты точки x. Тогда u (x)

sin в дв\

тор Бельтрами на сфере; F0 (в,^) — ар или диаграмма направленности поля F0 (в,^) = = u„ (x) = u„ (в,<?) .

Перепишем последние ряд и рекурсию в виде геометрооптического анзаца

u (x )=V £

Fn (6,ф)

(kr )n ,

Fn (М = ( 1 8

sin 6 86

• л 8 ) 1

srn6 — | + -

8

2 Л

86/ sin2 6 8ц>

+ n(n -1)

^ (в,р),

п = 1,2,3,...

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

Разложение результирующей ар в сложном падающем поле

Для эффективного вычисления суммарной ар в [1] предложен аналог теоремы разложимости Аткинсона—Уилкокса, представленный следующими выражениями, совершенно идентичными по структуре ряду и рекуррентной формуле из теоремы разложимости Аткинсона—Уилкокса, приведенным выше,

и. (*.* =£ ± ,

К0 п=0 (к«0)

где функции ¥п(60,ф0) с необходимыми параметрами (6,ф) определяются из следующего рекуррентного соотношения

F =

( 1 а (

sin0o дв0

sine„

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

_д_ да

+

0 J

+

1

sin в0 дде

д 2 Л ^

+ n(n -1) (2in) Fn-1,

a F0 равна

(x ) = J elkxd g (d) ds (d), x e R3

n

(5)

где g е Ь (О) ; О — поверхность единичной сферы; d = (smacos Д8та8т /3,со$а) — единичный вектор, совпадающий по направлению с волновым

вектором плоской волны, лежащий на поверхности единичной сферы О; g (d) — ядро Герглотца

функции V(х). В [13, с. 55] отмечается, что для

заданного g е Ь (О) функция

v(x) = Je-kxdg (d)ds (d), x e R3

n

(6)

также определяет волновую функцию Герглотца.

Для заданной функции g е Ь2 (О) решение проблемы рассеяния при падающей волне

(x) = J elkxdg (d ) ds (d), x e R3

Fo(e,P,eo,%) = ^o^oK(e,V,eo,Vo) .

Здесь um (в,ф, Ro,eo,yo) — результирующая ар; xo =( Ro,eo,^o) — координаты включения; u(eo,g>o) значение диаграммы направленности источника первичного поля, находящегося в начале координат по направлению (eo,g>o) на включение; um (e,V,eo,tyo) — каноническая ар включения при падающей плоской волне с единичным вектором

do =(1,eo,^o) .

Такое разложение справедливо, когда первичное падающее поле является излученным решением обычного уравнения Гельмгольца, и тогда сама обобщенная ар оказывается излученным решением уравнения Гельмгольца по переменным, совпадающим с координатами включения. Этот вывод следует из необходимых условий теоремы разложимости Аткинсона—Уилкокса, приведенной выше.

Таким образом, приведенное разложение связывает результирующую ар (аналог поля рассеяния в разложении Аткинсона—Уилкокса) с ар падающего поля и канонической ар включения (аналог ар в разложении Аткинсона—Уилкокса).

Функция Герглотца

Далее изложим основные определения и свойства волновой функции Герглотца (Herglotzwave function) [13, 2o, 21].

Волновая функция Герглотца (ВФГ) v (x) имеет следующее определение [13, c. 55]

дается выражением

,(x) = Jus (x,d)g (d)ds (d), x e R3 \ D

(7)

(8)

и характеризуется следующей амплитудой рассеяния:

,(x) = {uM(x,d)g(d)ds(d), x en .

(9)

Очевидно, что выражения (3) и (4) эквивалентны соответственно выражениям (7) и (9) при рассмотрении области видимости, т. е. на единичной сфере О . Отметим, что выражения (5) и (9) приведены также в [1, 19] независимо и для всей области определения, а не только для области видимости.

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

Для случая акустических волн ВФГ (5), (6) V (х) является целым решением уравнения Гельмгольца [13, с. 55], [20]

Av + k 2v = o.

(1o)

Кроме того, для функций и. (х, d) из (9) справедлив принцип взаимности для рассеяния, который выглядит как при объемном рассеянии, так и при многообразных краевых условиях на границе рассеивателя при поверхностном рассеянии следующим образом [11, с. 261], [13, с. 54, 55, 223, 285], [21], [22, с. 127]

v

v

v

,(x d) = М-^ -x) •

(11)

Это свойство тесно связано со следующей симметрией функции Грина для оператора Гельм-гольца (называемой также фундаментальным решением оператора Гельмгольца) [23, с. 746]:

G (х У ) = G (У, х ),

где G ( х, у ) и G (у, х ) определяются уравнениями

Ах G(х,у) + к2G(х,у) = 8(х - у),

Ау G(у,х) + к2G (у, х) = 8(у - х) .

Здесь Ах и Ау — операторы Лапласа, действующие по переменным х и у соответственно. Фундаментальное решение оператора Гельмгольца при принятом выше условии излучения Зоммер-фельда равно [24, с. 119, 204]

G ( x y ) = G ( y, x ) = G (Ix - y| ) = -

1 e

k|x-y|

4n x - y

(12)

и% (х, d) = к21G (х - у )т (у )и (у, d)dy . (15)

D

Получение асимптотики (2) из (14), (15) приводит к следующему виду канонической ар [13, с. 10]:

к 2 г

ит (х, d) = -—|в-кхут (у)и (у, d)¿>>. (16)

4П D

В случае поверхностного рассеяния при падении плоской волны в1 х каноническая ар поля и% (х, d) в общем случае определяется выражением [22, с. 499], [13, с. 20]

им (х, d) =

1 Л , ..ав-1кх'у 5и(у,^ -1кх/

Jju(x, d)

3n(y) 3n(y)

>ds( y).

(см. подробнее [23, c. 746-752]).

Отметим, что аналоги (5), (6), полученные в [1, 19], являются излученными решениями уравнения Гельмгольца.

О свойствах амплитуды рассеяния

[5], [11, § 3.2.4], [12-15 и др.]

В случае простейшей задачи объемного акустического рассеяния поле в области D подчиняется возмущенному уравнению Гельмгольца (см., например, [13, с. 2])

Au ( x ) + k2 n ( x ) u ( x ) = 0, (13)

где n (x) = c2 / c2 (x) — показатель преломления; c (x) — скорость звука, x e D . В случае, когда

(\ 1 kd -x

x) = e , решение

задачи рассеяния ищется из интегрального уравнения Липпмана—Швингера [13, с. 5]

u ( x, d ) = ekd x + k2 J G (x - y )m (y ) u (y ) dy , (14)

D

где u ( x, d ) определяется в (1) (в поле u ( x, d ) появился параметр d, отражающий неявную зависимость поля u, а значит, и поля рассеяния us и канонической ар um от направления распространения плоской волны ekd ' x ); m (y ) = 1 - n (y ) ; G(x - y) определяется из (12). Из (1) и (14) видно, что поле рассеяния us (x ) равно

Если u (x, d) обозначает давление p(x, d), то для um (xx) из (2) справедливо равенство [11, с. 109]

u„ d) = --1- х 4ж

х J |ikx • n(y)p(y,d) + 1®pn(y)•v(y,d)}e lkx-yds(y).

3D

Здесь v(y) — вектор колебательной скорости на поверхности 3D.

Соотношения в виде рядов по сферическим функциям (мультипольные представления)

Справедливо следующее утверждение [13, с. 34].

Пусть u(x) — излученное решение уравнения Гельмгольца вне шара |x| < R, R > 0. Тогда для u(x) справедливо разложение по сферическим функциям

да n

u (x) = Ц anm h n(1)(k| x| )C (x). (18)

n=0 m=-n

Полю (18) отвечает ар [13, с. 35]

1 да 1 n

u»(x) =1Iтп+TI OV"(x). (19)

k n=0 1 m=-n

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Здесь h n(1)(k |x| ) — сферическая функция Ханкеля первого рода с асимптотикой при большом аргу-

менте

h «'(kr) = I (-1 f" e- {1 + О ( ± )} ; anm —

u

искомые коэффициенты; Ynm (x) — сферические функции на единичной сфере, равные

Ynm (ад=А I2"+,1(n 7 'm! V1 e)em; p}m e —

]¡ 4л(п + |m|)! присоединенные полиномы Лежандра.

Выражение (19) следует из (18) при подстановке в нем асимптотик функций h n(1)(kr) при r ^ да и сравнении с выражением (2), являющимся фактически математическим определением ар.

В силу ортонормальности системы функций Ynm (x) на единичной сфере из (19) получаем

anm = k in+1 { u„ (x) (Y"m (x) )*ds(d) =

Q

2л л

= kin+1 Цu„ (в, y) (Y"m (в, y)) sin в de dy.

0 0

Звездочка в верхнем индексе в последнем выражении означает комплексное сопряжение.

В случае, когда падающая волна является плоской u (x) = e1 x, в выражениях (18) и (19) вместо u(x) и ида (x) должны фигурировать соответственно функции u(x, d) и um (x, d).

Согласно [13, с. 30], элементарными излученными решениями уравнения Гельмгольца являются следующие формы:

П Г (x) = hn(1) (k|x|)Ynm (в,ф), (20)

а элементарными целыми решениями уравнения Гельмгольца — элементарные формы

Л" (x) = J, (kl x| )Ynm (x),

u(x) = — J dßjM(d)e'kdx sina da ,

4n

л: (х)=]п к х| ^пт (в,ф)=

1 2я я

= (Ч)п — | Ynm(а,руы'х sшa da. (23) 4я 0 0

Тогда, если разложить образ й(6) по сферическим гармоникам

да п

= й(а, Р) = ХЕ (" (а, Р), (24)

п=0 т=-п

то окажется, что прообраз равен

(x) = к ZZ bnm Л nm (x).

u (x,

n=0 m=-n

Для bnm из (24) вытекает

(25)

(21)

являющиеся действительными составляющими форм (20) для излученного решения.

Разложение Уиттекера и Девани—Вольфа

[18]

РазложениеУиттекера. Уиттекером установлена следующая формула разложения по плоским волнам для целых решений однородного уравнения Гельмгольца, являющаяся преобразованием Фурье на сфере единичного радиуса:

k '

(22)

0 0

где d = (sinacosР,sinаsinР,cosа). Причем для элементарных целых решений Лпт (х) = = ]п (к х ^пт (х) имеет место равенство

bnm = in j dp\(Ynm(a,р))*«(a,P)sinada . (26)

о о

Как видно, целому решению (25) соответствует образ (24) с коэффициентами, определяемыми из (26). Разложение осуществляется по однородным плоским волнам с углами распространения из круга видимости ве[0,я], е[0,2я], что отвечает компонентам волнового числа kx = k sin^cos^ и ky = k sin^sin^, не превосходящим k, т. е. лежащим внутри круга К2 + ky 2 < k 2.

Равенство (23) означает, что с точностью до постоянного коэффициента функция Y]m (a, Р) является фурье-образом функции j; (kr)Yjm(в,ф) в круге видимости.

Сравнение функции Герглотца (5) и преобразования Уиттекера (22) с учетом того, что ds в (5) равно ds = sin adadp , показывает их совпадение с точностью до постоянного множителя.

Разложение Девани—Вольфа [18]. Преобразование Фурье от функции Y'¡" (a, Р) на всей плоскости ( kx, ky )е R2 равно [18]

П m (x) = h й(+) (k| x| )Ynm (в,ф) =

я .

Л --i да

• 2я 2

= (-i)n — f dp f Ynm(a,p)eMx sinada. (27) 2я 0 о

Показатель экспоненты в (27) равен d • x = |x|(sinesinacos(^-Р) ±cosecosa) ; знак + соответствует положительной ординате точки наблюдения x, а отрицательной — знак -.

В [18] используется версия h (к |х|) сферической функции Ханкеля, равная h (кг) = = ^n(1)(kr). Следовательно Ппт (х) = 1Ппт (х). Перепишем (27) в терминах функций h п(Г)(кг), внеся (- 1)п+1 под интеграл:

П пт (х) = h п(1) (к| х| )7пт (в,р) =

Л

--1 да

• 2л 2

— Гдр Г (- 1)п+17пт(а,р)вМхsmаdа. (27а) 2л

Сформируем в левой части (27а) поле (18), проведя соответствующее суммирование слева и справа в (27а), а также умножив и поделив справа на к :

XI anm h„(1)(kx|)Ynm(X) =

ki 2n

Ж

Л --1 да

2n 2

j dß j \т II (" i)n+1 anmYnm (a, ß)

k

0 k n=0 m=-n

x e1 kd'x sin a da.

Или иначе, учитывая (18) и (19):

(x) = — f dß f ида (d)e1 kd x sin a da . 9*7Г * *

2n 2

u(x) =

2n

(28)

Последнее выражение, связывающее ар ида с полем рассеяния и(х), хорошо известно (см., например, [22, с. 499, 500], [11, с. 122], [12]). Оно обычно записывается в двух формах: в сферических и в декартовых координатах:

(x) = — f dß f u(d)elkd x sin a da 9*7Г * *

Ж

2n ~

u(x) =

2n

2n

U(kx, ky) e'(kxx+kyy ±kzz) dkxdky. k „

—да —да z

n jj

Здесь u(d) = uM(d), а знак сверху формально отражает факт того, что функция u(d) является прообразом функции u(x); kz = ^jk2 - kx2 - ky2 =

= k cosa — вертикальная составляющая волнового вектора; в экспоненте знак + соответствует положительной ординате точки наблюдения x , а отрицательной знак -.

Замечание. Отметим, что в настоящей работе основное выражение (28) разложения Девани—

Вольфа получено на основе выражений (18) и (19), одинаково справедливых как для объемных, так и поверхностных рассеивателей в области x е R \ D, в то время как в [18] оно получено только для случая объемного рассеяния.

ВЫВОДЫ

Выше приведены и рассмотрены различные функциональные соотношения между ключевыми параметрами процесса рассеяния: полем рассеяния и амплитудой рассеяния. Таковыми соотношениями может быть либо соотношение геометроопти-ческого типа (разложение в виде ряда Аткинсона—Уилкокса), либо связь через представление обеих функций в ряды по сферическим функциям (мультипольные представления), либо в виде интегрального представления (разложение Девани— Вольфа). В первых двух случаях область определения обеих функций ограничивается т. н. областью видимости, в третьем случае область определения полярного угла распространяется на всю комплексную плоскость. Такое многообразие представлений возможно еще и вследствие аналитических свойств обеих функций: поле рассеяния является излученным решением уравнения Гельм-гольца, любое дважды дифференцируемое решение которого есть аналитическая функция своих аргументов [5, с. 83], а также тем, что амплитуда рассеяния является целой аналитической

функцией своих аргументов [11, §3.2.4, с. 126128]. Приведен также аналог разложения типа Ат-кинсона—Уилкокса для обобщенной амплитуды рассеяния, что возможно только в случае, когда падающее сложное поле является излученным решением уравнения Гельмгольца. В этом случае результирующая амплитуда рассеяния также подчиняется уравнению Гельмгольца. Показано, что волновая функция Герглотца с точностью до постоянного множителя совпадает с разложением Уиттекера. Приведенные результаты весьма полезны для приложений, и в частности, в задачах научного приборостроения.

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

1. Шарфарец Б.П. О возможности эффективного вычисления амплитуды рассеяния на включении в сложном поле // Акуст. журн. 2010. Т. 56, № 2. С. 166-171.

2. Atkinson F.V. On Sommerfeld's "Radiation conditions" // Philos. Mag. 1949. Vol. 40. P. 645-651. Doi: 10.1080/14786444908561291.

3. Wilcox C.H. Ageneralization of the orems of Rellichand Atkinson // Proc. Amer. Math. Soc. 1956. Vol. 7. P. 271276. Doi: 10.1090/S0002-9939-1956-0078912-4.

4. Wilcox C.H. An expansion theorem for electromagnetic

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

n=0 m=— n

X

n

да

fields // Comm. Pure Appl. Math. 1956. Vol. 9. P. 115— 134. Doi: 10.1002/cpa.3160090202.

5. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М.: Мир, 1987. 311 с.

6. Шарфарец Б.П. Радиационное давление при рассеянии произвольного поля на включении сложной формы // Акуст. журн. 2010. Т. 56, № 6. С. 767-772.

7. Шарфарец Б.П. Приближенный метод решения задач множественного рассеяния // Пятая Всероссийская научно-техническая конференция "Технические проблемы освоения Мирового океана". Владивосток, 30 сентября-4 октября 2013 г. С. 461-465.

8. Шарфарец Б.П. Приближенный метод решения задач многократного рассеяния в слоистых волноводах // XXVII сессия Российского акустического общества, Санкт-Петербург, 16-18 апреля 2014 г. С. 1309-1319. URL:

http://library.akin.ru/Rao/sess27/%D1%88%D0%B0%D1 %80%D1%84%D0%B0%D1%80%D0%B5%D1%86.pdf.

9. Шарфарец Б.П. Приближенный метод решения задач множественного рассеяния в полупространстве // Научное приборостроение. 2014. Т. 24, № 3. С. 75-79. URL: http://213. 170.69.26/mag/2014/abst3.php#abst9.

10. Шарфарец Б.П. К вопросу о приближенном методе решения задач множественного рассеяния. Решение на примере идеального волновода // Научное приборостроение. 2014. Т. 24, № 3. С. 80-86. URL: http://213.170.69.26/mag/2014/abst3.php#abst10.

11. Hansen T.B., Yaghjian A.D. Plane-wave theory of timedomain fields. N.Y.: IEEE Press, 1999. 367 p.

12. Шарфарец Б.П. Уточнение понятия "диаграмма направленности" // Акустические исследования жидкости с фазовыми включениями. Владивосток: ТОИ ДВНЦ АН СССР, 1984. С. 64-72.

13. Colton D., Kress R. Inverse Acoustic and Electromagnetic Scattering Theory. N.Y.: Springer, 1998. 331 p.

14. Алексеев Г.В., Бурштейн А.Б., Шарфарец Б.П. О некоторых свойствах диаграммных функций направленных излучателей // Электромагнитные и акустические процессы в океане. Владивосток: Изд-во ДВГУ, 1987. С. 130-141.

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

URL: http://213.17o.69.26/mag/2oo7/abst4.php#abst9.

16. Whittaker E.T. Math. Ann. 19o3. Vol. 57, no. 3. P. 333355. Doi: 1o.1oo7/BFo144429o.

17. Уиттекер Э.Т., Ватсон Дж.Н. Курс современного анализа. Т. II. Трансцендентные функции. М.: ГИФМЛ, 1963. 516 с.

18. Devaney A.J., Wolf E. Multipole expansions and plane wave representations of the electromagnetic field // J. Math. Phys. 1974. Vol. 15, no. 2. P. 235-244. Doi: Ю.Ю63/1.1666629.

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

20. Colton D., Kress R. On the denseness of Herglotz wave functions and electromagnetic Herglotz pairs in Sobolev spaces // Math. Meth. Appl. Sci. 2oo1. Vol. 24. P. 1289-13o3. Doi: 1o.1oo2/mma.277.

21. Colton D., Paivarinta L., Sylvester J. The interior transmission problem // Inverse Problems and Imaging. 2oo7. Vol. 1, no. 1. P. 13-28. Doi: 1o.3934/ipi.2oo7.1.13.

22. Морс Ф.М., Фешбах Г. Методы теоретической физики. Т. II. М.: Изд-во иностр. лит-ры, 196o. 886 с.

23. Морс Ф.М., Фешбах Г. Методы теоретической физики. Т. I. М.: Изд-во иностр. лит-ры, 1958. 93o с.

24. Владимиров В.С. Уравнения математической физики. М.: Наука, 1981. 512 с.

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

Контакты: Шарфарец Борис Пинкусович, sharb@mail.ru

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

ISSN 0868-5886

NAUCHNOE PRIBOROSTROENIE, 2017, Vol. 27, No. 3, pp. 128-136

ON THE RELATIONS CONNECTING THE SCATTERING FIELD WITH THE SCATTERING AMPLITUDE

B. P. Sharfarets

Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia

Various functional relationships between the main parameters of the scattering process are considered in this paper: the scattering field and the scattering amplitude. These relations are either geometrically-optic type (decomposition in the form of the Atkinson—Wilcox series), or a connection through the representation of both functions in series in spherical functions(multipole representations), or as an integral representation (Devaney— Wolf representation). Such a variety of representations is also possible due to the analytic properties of both functions: the scattering field is an radiation solution of the Helmholtz equation, any two differentiable solutions of which are analytic functions of their arguments, and also because the scattering amplitude is an entire analytic function of its arguments. An analog of an Atkinson—Wilcox representation for the generalized scattering amplitude is also given, which is possible only in the case when the primary incident complex field is a radiating solution of the Helmholtz equation. The resulting scattering amplitude in this case also obeys the Helmholtz equation. It is shown that the Herglotz wave function coincides, up to a constant factor, with the Whittaker representation. The above results are very useful for applications, and in particular, in the problems of scientific instrument making.

Keywords: scattering amplitude, Atkinson—Wilcox decomposition, Whittaker's representation, Devaney—Wolf representation, incident plane wave, the canonical scattering amplitude, resulting scattering amplitude, Herglotz wave function

REFERENCES

1. Sharfarets B.P. [The Possibility of Efficient Calculation of the Scattering Amplitude Caused by an Inclusion in a Complicated Incident Field]. Akusticheskij Zhurnal [Acoustic journal], 2010, vol. 56, no. 2, pp. 166-171. (In Russ.).

2. Atkinson F.V. On Sommerfeld's "Radiation conditions". Philos. Mag, 1949, vol. 40, pp. 645-651.

Doi: 10.1080/14786444908561291.

3. Wilcox C.H. Ageneralization of the orems of Rellichand Atkinson. Proc. Amer. Math. Soc., 1956, vol. 7, pp. 271-276. Doi: 10.1090/S0002-9939-1956-0078912-4.

4. Wilcox C.H. An expansion theorem for electromagnetic fields. Comm. Pure Appl. Math., 1956, vol. 9, pp. 115-134. Doi: 10.1002/cpa.3160090202.

5. Colton D., Kress R. Metody integral'nyh uravnenij v teorii rasseyaniya [Methods of the integrated equations in the theory of dispersion]. Moscow, Mir Publ., 1987. 311 p. (In Russ.).

6. Sharfarets B.P. [Radiation pressure in the case of the scattering of an arbitrary field by a complex-shaped inclusion]. Akusticheskij Zhurnal [Acoustic journal], 2010, vol. 56, no. 6, pp. 767-772. (In Russ.).

7. Sharfarets B.P. [Approximate method of solution of problems of multiple scattering]. Pyataya Vserossijskaya nauchno-tekhnicheskaya konferenciya "Tekhnicheskie problemy osvoeniya mirovogo okeana" [Fifth All-Russian scientific and technical conference "Technical Problems of Development of the World Ocean"]. Vladivostok, 30 September-4 October 2013, pp. 461-465. (In Russ.).

8. Sharfarets B.P. [Approximate method of solution of prob-

lems of multiple scattering in stratified wave guides]. XXVII sessiya Rossijskogo akusticheskogo obshchestva [XXVII session of the Russian acoustic society], Saint-Petersburg, 16-18 April 2014, pp. 1309-1319. URL: http://library.akin.ru/Rao/sess27/%D1%88%D0%B0%D1 %80%D1%84%D0%B0%D1%80%D0%B5%D1%86.pdf. (In Russ.).

9. Sharfarets B.P. [Approximate method of the solution of tasks multiple dispersion in the half-space]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2014, vol. 24, no. 3, pp. 75-79.

URL: http://213.170.69.26/mag/2014/abst3.php#abst9. (In Russ.).

10. Sharfarets B.P. [To the question of the approximate method of the decision problems of multiple dispersion. The decision on the example of the ideal wave-guide]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2014, vol. 24, no. 3, pp. 80-86.

URL: http://213.170.69.26/mag/2014/abst3.php#abst10. (In Russ.).

11. Hansen T.B., Yaghjian A.D. Plane-wave theory of timedomain fields. N.Y., IEEE Press, 1999. 367 p. (In Russ.).

12. Sharfarets B.P. [Specification of the concept "directional diagram"]. Akusticheskie issledovaniya zhidkosti s fazo-vymi vklyucheniyami [Acoustic researches of liquid with phase inclusions]. Vladivostok, TOI DVNTS OF ACADEMY OF SCIENCES OF THE USSR, 1984, pp. 64-72. (In Russ.).

13. Colton D., Kress R. Inverse Acoustic and Electromagnetic Scattering Theory. N.Y., Springer, 1998. 331 p.

14. Alekseev G.V., Burshtein A.B., Sharfarets B.P. [About some properties of chart functions of directional radia-

tors]. Elektromagnitnye i akusticheskie processy v okeane [Electromagnetic and acoustic processes in the ocean].Vladivostok, FEFU Publ., 1987, pp. 130-141. (In Russ.).

15. Sharfarets B.P. [On some properties of the scattering amplitude]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2007, vol. 17, no. 4, pp. 55-60.

URL: http://213.170.69.26/mag/2007/abst4.php#abst9. (In Russ.).

16. Whittaker E.T. Math. Ann,. 1903. vol. 57, no. 3, pp. 333355. Doi: 10.1007/BF01444290.

17. Whittaker, E.T., Watson, G.N. A course of modern analysis; an introduction to the general theory of infinite processes and of analytic functions; with an account of the principal transcendental functions. Cambridge, University Press, 1920, 632 p. (Russ. ed.: Whittaker E.T., Vatson Dzh.N. Kurs sovremennogo analiza. T. II. Tran-scendentnye funkcii. Translate G.M. Goluzin. Moscow, GIFML Publ., 1963. 516 p.).

18. Devaney A.J., Wolf E. Multipole expansions and plane wave representations of the electromagnetic field. J. Math. Phys, 1974, vol. 15, no. 2, pp. 235-244.

Doi: 10.1063/1.1666629.

19. Zatserkovnyi A.V., Sergeev V.A., Sharfarets B.P. [Em-

ployment of the Scattering Amplitude in Solving the Diffraction Problems for Waves in a Halfspace]. Akusti-cheskij Zhurnal [Acoustic journal], 2001, vol. 47, no. 5, pp. 650-656. (In Russ.).

20. Colton D., Kress R. On the denseness of Herglotz wave functions and electromagnetic Herglotz pairs in Sobolev spaces. Math. Meth. Appl. Sci., 2001, vol. 24, pp. 12891303. Doi: 10.1002/mma.277.

21. Colton D., Paivarinta L., Sylvester J. The interior transmission problem. Inverse Problems and Imaging, 2007, vol. 1, no. 1, pp. 13-28. Doi: 10.3934/ipi.2007.1.13.

22. Morse P.M., Feshbach H. Methods of Theoretical Physics. Part 2. McGraw-Hill Science/Engineering/Math, 1953. 1000 p. (Russ. ed.: Mors F.M., Feshbah G. Metody teore-ticheskojfiziki. T. II. Moscow, IIL Publ., 1960. 886 p.).

23. Morse P.M., Feshbach H. Methods of Theoretical Physics. Part 1. McGraw-Hill Science/Engineering/Math, 1953. 1997 p. (Russ. ed.: Mors F.M., Feshbah G. Metody teore-ticheskoj fiziki. T. I. Moscow, IIL Publ., 1958. 930 p.).

24. Vladimirov V.S. Uravneniya matematicheskoj fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1981. 512 p. (In Russ.).

Contacts: Sharfarets Boris Pinkusovich, sharb@mail.ru

Article received in edition: 24.04.2017

i Надоели баннеры? Вы всегда можете отключить рекламу.