К вопросу о приближенном методе решения задач множественного рассеяния. Решение на примере идеального волновода Текст научной статьи по специальности «Физика»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2014, том 24, № 3, c. 80-86

- МАТЕМАТИЧЕСКИЕ АНАЛИЗ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ

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

К ВОПРОСУ О ПРИБЛИЖЕННОМ МЕТОДЕ РЕШЕНИЯ

ЗАДАЧ МНОЖЕСТВЕННОГО РАССЕЯНИЯ. РЕШЕНИЕ НА ПРИМЕРЕ ИДЕАЛЬНОГО ВОЛНОВОДА

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

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

ВВЕДЕНИЕ

Задача учета множественного рассеяния в условиях влияния либо границ волновода, либо наличия множества рассеивателей, когда нельзя пренебречь многократным переотражением волн между ними, не является тривиальной. Эта тематика активно рассматривалась в последнее время различными авторами ([1-15] и др.). При этом для решения проблемы приходилось применять довольно громоздкие математические методы.

В настоящей работе предлагается приближенный достаточно простой метод решения задачи множественного рассеяния в плоскослоистом волноводе с границами. Метод основывается на полученном ранее автором результате [16], позволяющем достаточно просто рассчитывать амплитуду рассеяния (АР) произвольного включения при падении на него сложного первичного поля. При этом полагается известной каноническая АР, соответствующая условиям падения на включение плоской волны с текущего направления, имеющей единичную амплитуду и нулевую фазу в центре включения, где обычно помещается начало координат привязанной к включению сферической системы координат.

ПОСТАНОВКА ЗАДАЧИ

Алгоритм решения означенной задачи вкратце выглядит так. В волноводе присутствует источник

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

В настоящей работе описанный алгоритм, носящий универсальный характер, для простоты и наглядности применяется к простейшему случаю плоскослоистого волновода — однородному водному слою, ограниченному абсолютно отражаю-

щими плоскопараллельными границами с однородными условиями Дирихле и Неймана соответственно на верхней и нижней границах.

ОСНОВНЫЕ ПОЛОЖЕНИЯ РАБОТЫ [16]

Вначале кратко изложим основные, используемые ниже, положения работы [16]. Пусть включение находится в свободном пространстве, и начало координат привязано к включению (является либо его геометрическим центром, либо просто находится внутри области, занятой включением). Пусть далее йт (х, d) — каноническая АР, т. е. АР данного включения в случае, когда на него падает плоская волна, имеющая в начале координат единичную амплитуду и нулевую фазу с волновым вектором kd , где d = (smacos^этаэт fЗ,cosa) — единичный вектор, задающий направление распространения плоской волны, k — волновое число. Когда падающая волна итс(х) в общем случае представляет собой совокупность плоских волн со спектром g

c(x) =J g(d)e'Ad XdS(d):

(1)

то результирующая АР ида (X) вычисляется следующим образом:

,(X) = J g (d)ü„ (X, d)ds(d).

(2)

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

x

верхности единичной сферы;

X = т-г =

= (sm0cos^,sm0sm^,cos0) — вектор на единичной сфере, направленный в точку наблюдения; ит (х) — АР включения, характеризующая угловое распределение поля рассеяния на включении при его облучении полем ишс (х).

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

Гельмгольца только в бесконечной подобласти Я3 и удовлетворяет условиям излучения Зоммер-фельда на бесконечности (т. н. излученное решение), то интегрирование в (1), (2) происходит как

в области видимости (однородные плоские волны), так и вне ее (неоднородные плоские волны).

Отметим, что в случае, когда область Q совпадает с областью видимости ае[0,ж], ^е[0,2я], функция (1) носит название волновой функции Герглотца (Herglotz wave function) [17]. В этой же работе приведены некоторые свойства функции Герглотца.

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

U(kx,k) ...

v x' у/ ikd-x_

(x)=¿ íí^-dkA, =

—да —да z

ik_ 2ж

Ж

2ж 2

í dp í ü(d)emsmasmecos(p—')+cosacos 6) sin a da, (3) 0 0

выражение (2) записывается в виде

(в,', RoA'o) = ~~ í í ü(a,p)üx (6,',a,p) х 2ж J J

Ж

2ж ~2 —

(4)

о о

xe^R (smasm0o coe(p-% )+cosacos0o) sin a da dp Сравнение (4) с (2) дает ik

g(d) = {¡(a p)e1kRo(sinasin0oCos(P-^o)+cosacos0o) (5)

2n

В (3-5) u(kx,ky), {(d), {(a,P) — спектр плоских

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

x = xo = (xo. Уo, zo) = (RoA*o).

Далее, согласно [16], функция (в,*, Ro,eo,*o) из (4) является решением уравнения Гельмгольца относительно переменных (Ro,do,*o) и может быть представлена в виде ряда

M^RoA.*)=IT É ^Алх*^' (6)

Ro n=o (kRo)

где функции Fn(A,*,Ao,*o) определяются из следующего рекуррентного соотношения:

u

u

ü

□

X

x

F =-

n ^ •

2ш

а F0 равна

sin0o дв0 у

sine„

1 ^ да

+

o у

+-

1

д2

Л

sin в0 дф0

■ + n(n —1)

F

(7)

сти

O

i)

и выше для рассматриваемого

V /

здесь случая идеального волновода

Pinc = ^tJ— Z D (in ) sin (anz0 ) sin (anz) > H \nr n=0

x e

ж

"r-- 4 Je -1/2

C1/2 + O ((inr )-3/2 )■

да 1

•Z—

n=0 H\

2ж -i7 • i -e 4 sin (с

i r

n

x (e ' - e

c( e*

да .

Z An (r, Z0) ( el(<

,) D (in У

(4nr-anz)\ _

4r+anz) ei(inr—anz)

= (д,^,до,^о). (8)

О разложении (6) известно, что при ит (6,ф,во,фо) = 1 оно сходится абсолютно и равномерно при Яо > Я'; его можно дифференцировать по Яо, во, <ро почленно любое число раз и полученные после этого ряды также сходятся абсолютно и равномерно [18]. Здесь Я' — радиус наименьшего шара, включающего в себя область D источника первичной волны. При функции ит(9,ф,9о,фо), отличной от константы, этот вопрос требует отдельного рассмотрения, и наверняка можно сказать только то, что асимптотический ряд (6) по обратным степеням волнового параметра kЯо определяется рекурсией (7), (8) единственным образом в силу единственности представления асимптотическими рядами и удовлетворения функции иш(д,ф,Яо,до,фо) однородному уравнению Гельмгольца в указанной области. А следовательно, и в этом случае ряд (6) позволяет вычислить совокупную АР ит (д,ф, Яо,до,фо) с заданной точностью.

ПРИБЛИЖЕННЫИ МЕТОД РЕШЕНИЯ ЗАДАЧ РАССЕЯНИЯ

Рассмотрим идеальный волновод с абсолютно отражающими границами. Выражение для поля направленного излучателя в таком волноводе можно получить из выражений для общего случая слоистого волновода [19, 2о]. Опуская достаточно простые выкладки, выпишем искомое выражение для поля первичных нормальных волн с точностью до отбрасываемых слагаемых порядка мало-( з Л

где Н — глубина волновода; г = zо, г = о — точка локализации источника первичной волны;

(п +1/2)я

(r, z) — текущая точка;

С = -

H

п = о,1,2,... , — собственные числа поперечной

краевой задачи; Еп = у]к2 - ап2 , D(Еп) — коэффициенты возбуждения нормальных волн, обусловленные направленностью источника первич-, а

ных волн; к =--волновое число;

с

An Z0) =■

H V

2ж

— sin (ад)D (in) . (10)

i r

n

Как известно, приближение (9) справедливо в волновой зоне на больших по сравнению с длиной волны расстояниях. Тогда можно считать, что на включение, согласно (9), падает счетная система парных плоских волн с углами падения Е

дп = аг^—- и п-дп, п = о,1,2,.... Ограничимся " а

п

только числом N однородных нормальных волн,

2Н

определяемых условием N +1/2 <-< N + 3 / 2 .

Л

Углы дп, п < N, при этом остаются действительными. Выражение (9) для однородных нормальных волн перепишется в виде

(9а)

Рпс ~ IЛ (г, го) ( г**+-)- а)) .

п=о

Рассмотрим для простоты в качестве включения произвольное сферически симметричное тело (это допущение не является принципиальным и при необходимости легко распространяется на общий случай канонической АР). Каноническая АР сферически симметричного включения мм (ж, d) представляется в виде

,(©) = Z bmPm (cos ©).

(11)

Здесь ® — угол между единичными векторами d = (sin в cos gj,, sin а sin q>1, cos в) (направление дви-

д

1

)

e

жения плоской падающей волны) и x = т—¡- =

x

= (sin у cos р2, sin у sin р2, cos у) (направление из центра включения в точку наблюдения); bm — известные коэффициенты; Pm — полиномы Лежанд-ра; M — число, зависящее от волнового параметра ka и, как правило, выбираемое конечным; a — характерный размер включения.

Как известно, косинус угла ® между двумя единичными векторами x и d равен

cos © = cosycos0 + sin у sin в cos(p2 -р). (12)

Полагая, что азимут из точки источника первичного поля на точку включения р = 0, окончательно получаем

cos © = cosycose + sin у sin в cos р2. (13)

Первичная падающая волна приближенно определяется рядом (9а), (10), т. е. представляет собой сумму взвешенных плоских волн (в данном случае функция g (d) из (1) представляет собой набор взвешенных дельта-функций). Тогда, если на включение падают плоские волны из набора (9а) с векторами нормали d± = (sin вп ,0, ± cos вп) (знак (+) для углов падения вп и знак (-) для п -вп), то соответствующая ей результирующая АР с учетом (2), (9а) и (11) будет иметь вид

M /

Ду,Р2) = An (Г \ Z0) Z(b™ Х

m=0 V

х( el('

Pm (cos ©+) - e11

P

(cos ©- ) )).

Здесь (r', z') — точка локализации включения; cos ©+ = cos у cos вп + sin у sin вп cos р2; cos ©- = = - cos у cos вп + sin у sin вп cos р2; An (r', z0) определяется из (10), а также принято во внимание, что Р= 0.

Теперь, используя последнее выражение, легко получить результирующую АР включения

uM(1) (у,р2), вызванную первичным полем стороннего источника (9а)

N M

N N M /

ю(1) (у,Р2) =Е п (у,Р2) =ЕЕ(У

п=0 п=0 m=0 *

Ап (r z0)bm Х

х^^'+^Ч (cos ©;) - (cos ©-))). (14)

Включение с АР (14), рассеивая первичные волны, является в свою очередь источником вторичного поля, которое также может рассеиваться

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

Для облегчения задачи решим проблему создания третичного поля только самим включением за счет однократного рассеяния вторичного поля от границ волновода, пренебрегая при этом рассеянием на источнике. Таким образом, речь идет

о нахождении добавки и^2^(у,р2) к ар им^(у,р2), вызванной рассеянием на включении однократно отразившегося от границ вторичного поля.

Итак, поиск АР и^2\у,р2) формально сводится к следующей задаче: задан источник вторичного поля с АР иг}1\у,р2), совмещенный с включением, задана каноническая АР включения мм (©) из (11). Необходимо найти добавку иг}2\у,р2), состоящую из двух слагаемых, каждое из которых выражает вклад верхней и нижней границы волновода по отдельности при условии однократного отражения от каждой из границ вторичного поля с АР и/>(у,р2).

Эта задача редуцируется к решению двух самостоятельных задач. В одной из них необходимо решить задачу однократного рассеяния в жидком полупространстве 2 > 0 (ось аппликат направлена вниз) поля источника с АР им(1)(у,р2) на включении, совмещенном с источником. Коэффициент отражения от границы 2 = 0 обозначим для общности через V1 (в настоящей работе принято VI =-1). Акустические свойства жидкого полупространства совпадают со свойствами жидкого слоя. Во второй задаче необходимо решить то же только в полупространстве г < Н с коэффициентом отражения V2 на границе г = Н (в настоящей работе принято V2 = 1). Включение с источником совмещены в пространстве и локализованы в точке (г = 0, г') . Эти задачи легко решаются с применением техники, изложенной в работе [16] и коротко приведенной в начале настоящей работы. Пример подобного решения, только для случая, когда им(1)(у,р2) соответствует сферической волне, приведен в [21]. Здесь приведем коротко решение для случая полупространства г < Н , а для полупространства 2 > 0 выпишем конечный результат.

Итак, источник, он же рассеиватель, локализован в точке (г = 0, г. Граница проходит по плоскости 2 = Н . Задача эквивалентна следующей.

u

r -a.. z

u

Оставляя источник поля с АР им(1)(у,ф2) в исходной точке (г = о, г'), мысленно перемещаем рас-сеиватель зеркально относительно границы г = Н в точку (г = о, 2Н - г') . Таким образом, расстояние между источником и рассеивателем равно Яо2 = 2 (Н - г') . Угол падения с источника на рас-сеиватель до2 = ж . Для слагаемого результирующей АР, учитывающего влияние только нижней границы, им(2)2(у,ф2), при использовании выражениями (6)-(8) учтем только нулевое приближение:

е1кЯо2

и«( 2(Уф F0(r,^2,д0,^1) =

>ikR02

R02

R0.

(15)

Здесь и(до = ж,фо) определяется через им(1)(у,ф2) из (14) и коэффициент отражения V2 (см. [21]) следующим образом:

й(в0 = Ж,ф0) = йда(1) (У = ж,ф2)¥2(ж) = = uM(1) (у = ж)У (0).

(16)

йда (0)1 в0 =ж = йда (У,Ф2,в0,ф)\в0 =ж

м

= Z bmPm (— СоУ

(17)

Учитывая (16) и (17), получаем для йда(2)2(у,ф2):

йда(2)2(У,ф2) йда1(У = Л)У2(0) X

R02

йда(2)1(У,Ф2) ^^ йда(1) (у = Ж^(0) X

к

Z bmPm (cosy) = йД(у). (19)

Здесь Яо1 = 2г'.

Теперь можно выписать окончательное выражение для АР с учетом первичного и вторичного полей:

им (г, Ф2) = им(1) (г, Ф2) + им(2)1 (г) + им(2)2 (г), (2о)

где соответствующие функции определяются выражениями (14), (18) и (19).

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

в рассматриваемом волноводе создает поле (9). Включение, имеющее результирующую АР (2о) и находящееся в точке (г ',ф = о, го) , создает поле рассеяния р., которое может быть рассчитано по известной АР обычным порядком [19, 2о, 22]. Суммарное поле р = ртс + р. в произвольной точке волновода нужно рассчитывать с учетом отсутствия цилиндрической симметрии геометрии задачи.

Несмотря на упрощения рассмотренного частного случая волновода, предложенный подход может быть легко использован в слоисто-неоднородном волноводе при условии, что слой жидкости, содержащий включение, является однородным. В этом случае первичное поле также будет представлять совокупность парных плоских волн, а коэффициенты отражения на границах слоя могут быть легко определены [23]. Кроме того, задача легко решается и тогда, когда излучатель и включение совмещены [21]. При отсутствии сферической симметрии включения никаких сложностей не возникает. В этом случае каноническая АР будет зависеть не от одного угла, как в (11), а от четырех в общем виде йх (у,ф2,до,ф1) (см., например, (15).

ВЫВОДЫ

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

(2)

(2)

В (16), как видно из (12), (13), cos ® перестает зависеть от ф2, а функция V(ж) = V(0). Далее из

(11)—(13) имеем cos©I =— cosy, а значит, из

(11) имеем

ЬР (-СО.у) = им(2)2(Х). (18)

т=о

Из (18) видно, что им(2\(у,ф2) зависит только от угла падения у.

Отсюда же очевидно также, что при отражении от границы г = о функция им(2)1(у,ф2) должна выглядеть следующим образом:

e^ (1) (у _.

m=0

kR02

ния, но легко обобщается на общий вид включения.

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

1. Кравцов Ю.А., Кузькин В.М., Петников В.Г. Дифракция волн на регулярных рассеивателях в мно-гомодовых волноводах // Акуст. журн. 1984. Т. 30, № 3. С. 339-343.

2. Hackman R.H., Sammelmann G.S. Acoustic scattering in an inhomogeneous waveguide: Theory // J. Acoust. Soc. Am. 1986. Vol. 80, N 5. P. 1447-1458.

3. Sammelmann G.S., Hackman R.H. Acoustic scattering in a homogeneous waveguide: Theory // J. Acoust. Soc. Am. 1987. Vol. 82, N 1. P. 324-336.

4. Hackman R.H., Sammelmann G.S. Multiple-scattering analysis for a target in an oceanic waveguide // J. Acoust. Soc. Am. 1988. Vol. 84, N 5. P. 1813-1825.

5. Белов В.Е., Горский С.М., Зиновьев А.Ю., Хиль-ко А.И. Применение метода интегральных уравнений к задаче о дифракции акустических волн на упругих телах в слое жидкости // Акуст. журн. 1994. Т. 40, № 4. С. 548-560.

6. Gaunaurd J.C., Huang H. Acoustic scattering by a Spherical body a plane boundary // J. Acoust. Soc. Am. 1994. Vol. 96, N 4. P. 2526-2536.

7. Gaunaurd J.C., Huang H., Strifors H.C. Acoustic scattering by a pair of spheres // J. Acoust. Soc. Am. 1995. Vol. 96, N 1. P. 495-507.

8. Лебедев A.B., Хилько А.И. Рассеяние плоской волны на двух упругих шарах и сферических оболочках // Акуст. журн. 1997. Т. 43, № 4. С. 521-530.

9. Кузькин В.М. Об излучении и рассеянии звуковых волн в океанических волноводах // Акуст. журн. 2001. Т. 47, № 5. С. 678-684.

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

11. Кузькин В.М. Рассеяние звуковых волн на теле в плоскослоистом волноводе // Акуст. журн. 2003. Т. 49, № 1. С. 77-84.

12. Шарфарец Б.П. Метод расчета поля непрозрачного источника и поля рассеяния неоднородного включения в плоскослоистой среде // Акуст. журн. 2004. Т. 50, № 1. С. 123-128.

13. Шарфарец Б.П. Метод решения задачи множественного рассеяния электромагнитных волн в однородной безграничной среде // Научное приборо-

строение. 2004. Т. 14, №. 4. С. 57-69.

14. Шарфарец Б.П. Метод решения задач множественного рассеяния на нескольких телах в однородной безграничной среде // Акуст. журн. 2005. Т. 51, № 5. С. 672-681.

15. Григорьева Н.С., Фридман Г.М. Рассеяние звука сферической упругой оболочкой, помещенной в волновод с жидким дном // Акуст. журн. 2013. Т. 59, № 4. С. 424-432.

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

17. Colton D., Kress R. Inverse acoustic and electromagnetic scattering theory. Springer, 1998. 334 p.

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

19. Шарфарец Б.П. Поле направленного излучателя в слоисто-неоднородном волноводе // Акуст. журн. 1985. Т. 31, № 1. С. 119-125.

20. Шарфарец Б.П. Поле протяженного источника в регулярном океаническом волноводе // Акуст. журн. 1991. Т. 37, № 4. С. 794-799.

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

22. Шарфарец Б.П. Поле протяженного направленного излучателя в регулярном океаническом волноводе // Акуст. журн. 1989. Т. 35, № 1. С. 132-137.

23. Бреховских Л.М. Волны в слоистых средах. М.: Наука. 1973. 343 с.

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

Контакты: Шарфарец Борис Пинкусович, [email protected]

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

86

E. n. fflАP®АPЕЦ

UDK 534.874+534.26

TO THE QUESTION OF THE APPROXIMATE METHOD OF THE DECISION PROBLEMS OF MULTIPLE DISPERSION. THE DECISION ON THE EXAMPLE OF THE IDEAL WAVE GUIDE

B. P. Sharfarets

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

In work on the example of the elementary liquid layer with flat borders the general rather simple approximate algorithm of calculation of the field formed by dispersion of primary field radiated by a third-party source on inclusion and dispersion on inclusion of a secondary field after its single reflection from borders of a wave guide is offered. The secondary field is understood as the field of dispersion caused by falling on inclusion of primary field. The task is considered for inclusion in the form of any spherically a symmetric body, however this assumption is easily generalized on any inclusion. The solution of this task was possible thanks to the algorithm of stay of total amplitude of dispersion received earlier by the author when falling on inclusion of a wave of a difficult form. The algorithm can be used in layered and non-uniform wave guides on condition of uniformity of the layer containing inclusion. The problem is easily solved and when the radiator and inclusion are combined.

Keywords: ideal wave guide, layered wave guide, primary field, spherically symmetric inclusion, dispersion field, dispersion amplitude, uniform half-space, flat limits of the section

REFERENCES

1. Hackman R.H., Sammelmann G.S. Acoustic scattering in an inhomogeneous waveguide: Theory. J. Acoust. Soc. Am., 1986, vol. 80, no. 5, pp. 1447-1458.

2. Sammelmann G.S., Hackman R.H. Acoustic scattering in a homogeneous waveguide: Theory. J. Acoust. Soc. Am., 1987, vol. 82, no. 1, pp. 324-336.

3. Hackman R.H., Sammelmann G.S. Multiple-scattering analysis for a target in an oceanic waveguide. J. Acoust. Soc. Am., 1988, vol. 84, no. 5, pp. 1813-1825.

4. Gaunaurd J.C., Huang H. Acoustic scattering by a

Spherical body a plane boundary. J. Acoust. Soc. Am., 1994, vol. 96, no. 4, pp. 2526-2536.

5. Gaunaurd J.C., Huang H., Strifors H.C. Acoustic scattering by a pair of spheres. J. Acoust. Soc. Am., 1995, vol. 96, no. 1, pp. 495-507.

6. Colton D., Kress R. Inverse acoustic and electromagnetic scattering theory. Springer, 1998, 334 p.

Contacts: Sharfarets Boris Pinkusovich, [email protected]

Article arrived in edition: 21.05.2014