ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2007, том 17, № 1, c. 62-72
ИССЛЕДОВАНИЯ, ПРИБОРЫ, МОДЕЛИ-И МЕТОДЫ АНАЛИЗА
УДК 534.23 © Б. П. Шарфарец
К ВОПРОСУ О ВЫЧИСЛЕНИИ АМПЛИТУДЫ РАССЕЯНИЯ ОБЪЕМНЫХ И ПОВЕРХНОСТНЫХ РАССЕИВАТЕЛЕЙ
В работе приводятся некоторые подходы к вычислению амплитуды рассеяния акустических рассеивателей, применяемые в различных областях теории рассеяния волн и частиц. Рассматриваются как объемные, так и поверхностные рассеиватели. Приводятся примеры конкретных расчетов, проведенных различными методами.
ВВЕДЕНИЕ
Одной из важнейших характеристик использующихся в задачах рассеяния волн и частиц является амплитуда рассеяния (ар). Ар является как бы "паспортом" рассеивателя, характеризующим его свойства. С помощью этой характеристики решаются всевозможные прямые и обратные задачи теории и практики рассеяния. В частности, знание ар отдельных рассеивателей является необходимым условием для проведения вычислений по приведенному в работе [1] и других работах автора методу учета многократного рассеяния в акустике и электродинамике. В настоящей работе некоторые подходы к вычислению ар, применяемые в различных областях теории рассеяния, адаптируются к нуждам теории рассеяния в акустике. Рассматриваются как объемные, так и поверхностные рассеиватели.
АМПЛИТУДА РАССЕЯНИЯ ОБЪЕМНЫХ РАССЕИВАТЕЛЕЙ.
ОБЩИЕ ВЫРАЖЕНИЯ
Рассмотрим случай объемного рассеивателя, вызванного наличием в однородном пространстве Я3 некоторой ограниченной области Б, включающей в себя начало координат, с возмущенными плотностью среды р1(х) и скоростью звука сДх), х е Б:
c = const, x е R3 \ D; Ci(x), xе D;
p = const, x е R3 \ D; p(x), x е D,
(1)
где си р — невозмущенные значения скорости звука и плотности среды пространства Я3.
Запишем уравнение для звукового давления р(х) в этом случае [2]:
2 1
Ap (x) + к2 n(x) p(x)--Vp(x) • Vp (x) = 0,
P(x)
(2)
x e R3
где х = (х, у, z); к = ю/ с — волновое число; п(х) = с2/ с 2(х) — показатель преломления, равный единице в области Я3 \ Б . Заметим, что в (2) предполагается дифференцируемость плотности р. После замены и = р / р получаем:
Аи + К2 (х)и = 0, (3)
где
1 3 f 1 ^
K2 (x) = к2 n(x) +-Ap---Vp
2p 4 p
(4)
Перепишем (3) с учетом (1) и (4) в виде неоднородного уравнения Гельмгольца:
Au + к 2ы =
к 2(1 - n( x)) + -
f
± Vp
p
- fAp
2p
u. (5)
Решение (5) при условии наличия первичной падающей волны ui ищется формально из интегрального уравнения (ИУ) Липпмана—Швингера, которое при стандартном допущении о том, что ui представляет собой плоскую волну ui = вЛк,х, где к {— соответствующий ей волновой вектор, имеет вид [3]
u (x) = e'k,x - к2 J
ехрОк x - y)
——г1-r~m(y)u(y)d 3
4п x - y\ (6)
xe R3, где с учетом (5)
m(y) = 1 - n + —
к2
f
3 1
1
--Vp--Ap
4 p 2p
ехр(/£| х - y|) G(x, y) =-1-l.
фундаментальное ре-
4п |х - y|
шение уравнения Гельмгольца.
Отметим, что p(x) = u (х) в области х е R3 \ D . Наконец, полагая р(х) = р = const, х е R3, приходим к обычной форме уравнения Липпмана— Швингера [3] для m(х):
m(x) = 1 - п(х) = 1 -
сДх)2
(7а)
В работе [3, теорема 8.7] показано, что уравнение Липпмана—Швингера (6) при всяком к > 0 имеет единственное решение, непрерывно зависящее от амплитуды падающей волны ui.
Из (6) следует выражение для ар [3], получаемое предельным переходом |х| ^ ^ под знаком интеграла и выделением множителя при сфериче-
Л\х|
ской волне
u (х) = ui + us =
ik тх
= e ' +•
2п m, п = 1,2.
J
><rr J
exp(ik П х)
a
-г: (к,. (9)
Здесь k1 = (kx, ky ,a), k2 = (kx, ky, -a), a =
= ^к2 - кх2 - ку2 , а Тт (кi.,к^) определяется из (8).
При т = 1 плоская падающая волна распространяется снизу вверх в сторону увеличения координаты г, при т = 2 — в сторону уменьшения координаты г . При п = 1 рассеянная волна распространяется снизу вверх в сторону увеличения координаты г, при п = 2 — рассеянная волна рас-
пространяется сверху вниз в сторону уменьшения координаты г .
Приближение Борна для ар имеет вид [3, 6]:
Тт (к., к,)=Тт (к т - к п)= к2
= -— / ехр(.(кт -кП)у)т(у)а3у. (10) 4п и3
Как видно из (10), ар в приближении Борна за-
V.» V.» ■ т 1 п
висит только от одной переменной р = к { - к ^, равной вектору, соединяющему две точки на сфере радиусом к .
Следуя идее, изложенной в работе [6, с. 76],
к 2
умножим обе части в (9) на--ехр(-.к1 х)т(х),
4п
1 = 1,2 и проинтегрируем по х в области В . В результате получаем уравнения, связывающие ар (8) и ар в приближении Борна (10):
к2
Т(к.,к,) = -—] ехр(-к,у)т(у)и(у)а3у . (8)
4п В
Здесь учтено, что функция и (х) из (6) и (8) зависит от направления падения плоской волны и (х) = и(х, к ) . Из (8) очевидно, что нахождение ар в общем случае является нетривиальной задачей, т. к. для этого предварительно необходимо решить интегральное уравнение (6) и вычислить функцию и(х, к {) для всех необходимых значений вектора к..
Выражение для ар также можно получить, воспользовавшись разложением сферической волны в (6) по плоским волнам [2] и техникой, изложенной в работах [4, 5]. В этом случае выражение (6) перепишется в виде
т: (к, к)=т: (к, к)+
0
, i fTn' (к,, к)
nV - '-T2 (к, к, )dkxdky; (11)
2п ¿2 ал 1, т, п = 1,2.
Для получения интегрального уравнения необходимо в (11) потребовать выполнения условия п = 1. В результате получаем интегральное уравнение
0
тт (к., к)=тт (к,, к)+
0
+
i J Т/ (к,, к T
R
2п r a
l,m = 1,2.
Tlm (кЛ^М^; (12)
Отметим, что ИУ в виде (12) отличаются от полученных в [6] и других работах тем, что в (12) рассматривается не в 3-, а в 2-мерном импульсном пространстве.
Теперь, предварительно вычислив ар в приближении Борна (10), можно из (12) вычислить искомую амплитуду рассеяния. В работе [6] приведен пример решения ИУ (12) для потенциала Юкава.
Как видно из (12), ар Тт (к i, к) представляет
собой разрешающее ядро (резольвенту) ядра
0
Тт (к i, к), совпадающего с ар в приближении Бор-на, что может быть использовано для изучения свойств ар Тт (к i, к).
Таким образом, задача нахождения ар в общем случае может быть сведена либо к решению инте-
2
c
х
грального уравнения (6) для и(х), х е Б с последующим вычислением ар из выражения (8), либо напрямую из интегрального уравнения (12).
РЕШЕНИЕ ИСХОДНОГО ИНТЕГРАЛЬНОГО УРАВНЕНИЯ (6)
Приведем некоторые методы решения исходного интегрального уравнения (6) с учетом того, что достаточно найти решение и(х) только в области Б. После этого можно воспользоваться выражением (8) для вычисления ар, с помощью которой можно затем найти поле рассеяния [1] (такой подход является альтернативным методу нахождения поля рассеяния путем решения интегрального уравнения (8) в области Я3 \ Б (см. библиографию в [3])). Отметим, что ядро интегрального уравнения (6) имеет слабую особенность, поэтому его необходимо заменять близким по норме непрерывным ядром.
Аппроксимация ядра интегрального уравнения близким по норме ядром
Ядро К (х, у) ИУ (6) является ядром со слабой особенностью, оставаясь при этом фредгольмов-ским [7, 8]:
) exp(/£|x -y|) ( ) K (x, y) = —¡-^ m(y).
4п x - y
(13)
L(x, y) =
ex
p(ik |x - y|)
4n |x - y| exp(ik |x - y|)
m(y), |x - y\ >e;
(14)
4ne
m(y), |x - y|
<e,
где £ — некоторая малая константа, норма разностного ядра И(х, у) = К (х, у) - Ь(х, у) равна
||*(x, y )|| „ = max J| R(x, y)\ly
<-
.Me'
(15)
где M = ||m(x)||^ = max|m(x)|, что при e ^ 0 и при
достаточно малых величинах M дает хорошее приближение.
Пусть далее ядро L из каких-либо соображений в свою очередь необходимо аппроксимировать другим (к примеру, вырожденным) ядром Ln (x, y). Тогда норма результирующего разностного ядра удовлетворяет неравенству
= |\K(x, y) - Ln (x, y)||^ =
= |\К(х, у) - Ь(х,у) + Ь(х,у) - Ьы(х, у<
<1 И -+| К| L=a, (16)
где 1К|[ =11Ь - I, и при малом £ Ц^Ц_ практически не превосходит нормы ядра = Ь - Ьы . После определения нормы ЦИ можно воспользоваться различными оценками [7, 9, 10] для получения нормы разности решений интегрального уравнения соответственно с исходным ядром и аппроксимирующим ядром. Из всего многообразия оценок приведем наиболее простую [11, с. 359], которая применительно к ядру (13) имеет вид
<(x) - Un (x)| xeD = O(a).
(17)
В работе [8, с. 299] показано, что при замене ядра со слабой особенностью (13) непрерывным ядром
Здесь и(х) — решение инт. уравнения (6); иы (х) — решение инт. уравнения вида (6), где вместо ядра К фигурирует Ьы. Отметим, что в указанных работах приведены более точные оценки невязки решения инт. уравнения при аппроксимации исходного ядра. Впрочем, в рамках настоящей работы ограничимся оценкой (17).
Далее рассмотрим некоторые методы решения интегрального уравнения (6).
Метод последовательных приближений
Одним из часто применяющихся является метод последовательных приближений. Пусть В — шар минимального радиуса а с центром в начале координат, описывающий включающую в себя все неоднородности область Б . Тогда (6) для х, у е Б может быть переписано в виде (в силу того, что т(у) = 0 в области В \ Б)
( ) k 2fexp(ik |x - y|) ( ) ( з u(x) = ut - k I , '-r~m(y)u(y)d y,
в 4n x - y\
(18)
x е D.
Как показано в [3], метод последовательных приближений сходится к точному решению при условии
(19)
ka <
В этом случае справедлива оценка [9]
(Ma2 k2/2)r+1 1 - Ma2 k2/2 .
u (x) - ur (x) I <■
Здесь
,(x) = J(-1) k2nKn[u,](x),
(20)
(21)
n
где
Кп[щ ](х) =
= ИIв (х, г^тЫв (11, 2)...
в [_ в в
...в(1 п-1, у)т(1 п |х
п-1
X т(у )и. (у)ф>. В случае если условие (19) не выполнено, к > 1 [2
то для к >—л— резольвенту интегрального опе-а\Ы
ратора (18) Я(х, у, к2 +Я) можно вычислить с помощью аналитического продолжения [10]
Я( х, у, к02 +Я) =
= Я(х, у, к02) + ЛЯ2(х, у, к02) + Я2 ЯДх, у, к02) + ...,
где Яп — повторные ядра резольвентного ядра Я; к0 — некоторое значение к , для которого выполняется условие (19).
Ниже на рис. 1 приведена зависимость величи-
V2/M
12
Но
10
20
SO 40
С1/200, м/с
Рис. 1. Зависимость л от скорости звука в неоднородном включении
при c1 е[1225,2122] до ka при c1 ^ 1500±0. Вне этой окрестности ka < 2 .
Редукция исходного трехмерного интегрального уравнения к системе одномерных
В задачах с радиальной симметрией, когда :(х) = m(r), возможна редукция трехмерного ИУ (6) к системе одномерных [12, 13 и др.]. В этом случае в ИУ (6) используется разложение [3, с. 29]
= f £ (2n + D/n (r, r)Pn (cos 0). (22) 4п|х - y| 4п n=0
где
/n (r, r) =
Ijn (kr )h n (1)(kr)
j (kr )hn(1)(kr)
r < r , r > r ;
ны . — от значения c1 для случая, когда в (7а) \M
с1(х) = c1 = const. При этом c = 1500м/с. В точке c1 = 1500 на графике сингулярность. Очевидно, что в окрестностях этой точки ka из (18) может варьироваться в широких пределах: от ka = 2
г = х ; г = у ; 0 — угол между х и у ; Рп — полиномы Лежандра; ]п, Ьп(1)— сферические функции Бесселя соответственно первого и третьего рода. В реальных вычислениях необходимо в (22) оперировать конечной суммой. В качестве аппроксимирующего ядра Ьы возьмем конечную сумму
^ (х у) =
ik 4п
m(r )£(2n + 1)/, (r, r ) Pn(cos 0). (23)
В качестве L примем ядро (14). Для оценки N необходимо вычислить норму разностного ядра |L - Ln|| . Анализ выражения (22) показывает, что модуль остаточного члена ряда как функции y при r = |х| = const и при R = |х - y| = const возрастает с ростом r. При r = 0 ряд (22) вырождается в единственный член ряда с индексом n = 0, и тем самым ^Дх, y) точно описывает ряд (22) (остаточный член тождественно равен нулю). На рис. 2 приведены зависимости 50х|Lj0| и |L500| в центре
сферической волны |х - y| = 0 от величины r удаления этого центра от начала координат. Из рисунка видно, что, во-первых, точность падает с ростом r , а во-вторых, достаточно точно выполняется зависимость N2 х |ln1 (r, r)| = N1 х |ln^ (r, r)|,
N * 0.
На рис. 3-5 представлены зависимости модуля разностного ядра RLN в зависимости от N
и r1, где r1 = |y| , а точка y лежит на луче, соединяющем начало координат и точку х .
т, Чг
500
400 300 200 100
Цъ, 44
Рис. 2. Зависимости п1 = 50 х |Ь10(г, г ) (и П = Ь500(Г, г ')| (-) от г
45, Ц6
50
Рис. 4. Зависимость п5 = Я, ( ■
' 5 Ь50
и П = Я^ (-) пРи г = 5
10 г, м
51 52
^/10, м
)
50
Рис. 3. Зависимости п3 = ЯЬ (и П 4 = ИАоо (-) пРи г = 1
Пъ Цъ
89 90
Рис. 5. Зависимость п7 = Я^ ( -и П = Я, (-) при г = 9
51 52
Г:/10, м
- )
91
92
^/10, м
)
)
На рис. 3-5 утолщенные и обычные кривые соответственно рассчитаны с одинаковой точностью. Видно, что по мере роста г = |х| (г соответственно равно 1, 5 и 9) число членов суммы N в (23) монотонно растет для утолщенных и обычных кривых соответственно в пропорции 9 : 50 : 95 и 100 : 180 : 310. Очевидно, что при достаточно больших N разностное ядро существенно отлично от нуля только в ближайшей окрестности центра сферической волны г = г1 и чем дальше этот центр от начала координат, тем большим должно быть N для обеспечения заданной точности.
Во всех приведенных выше расчетах было принято к = 1. Отметим, что оценку нормы остаточного члена (15), (16) достаточно уметь вычислять для к = 1 . Действительно, норма разностного ядра (15) не зависит от к , а для нормы (16) легко показать, что имеет место следующее равенство:
Л ЯЬя (к, х, у фу = к-Г IК (к = 1, X, ¥)|аг.
В к кВ
Здесь X = кх ; ¥ = ку ; В = кВ — шар радиусом, в к раз превышающим радиус шара В . От-
метим, что величина £, фигурирующая в непрерывном ядре (14), при переходе к координатам (k = 1, X, Y) трансформируется в величину £ = k£ , а шар B£ ^ Bk£.
После определения величины N можно записать ИУ с аппроксимирующим ядром LN :
un(х) = e^ - k2 JLn(х,yК(y)dy, (24)
в
где
ln (х у ) =
k N
m(r ')£ (2n +1)/, (r, r ) Pn (cos 0). (25)
4n n=0
В рассматриваемом случае область D представляет собой шар B радиусом a с центром в начале координат (радиус a может быть и бесконечным) и :(х) = m(r). Положим в силу симметрии задачи, что плоская первичная волна распространяется вдоль положительного направления оси z . Тогда (6) преобразуется к виду
u(r) = eikz - k2 J
exi
p(ik | r - r '|)
4n r - r'
m(r ')u (r ')d3r'
(26)
г е И3.
Как показано в [6, 13], решение интегрального уравнения (26) сводится к решению счетной системы одномерных интегральных уравнений вида
X (r) = j„ (kr) + k2 J gn (r, r ')m( r ')%„ (r ')dr';
(27)
n = 0,1,2,...; rе [0,-).
Здесь gn (r, r') =
1
k
Jn (kr)nn (kr '), r < r',
ция Грина задачи;
i!+k 2+nin^l'
dr2 r2
V
n = 0,1,2,...;
—j„ (kr ')nn (kr), r > r' k
„ (r, r') = 5(r - r'),
(28)
Jn(kr)=
nkr
J 1(kr) = kr jn (kr);
2 n+—
2
Пп (кг) = (-1)п+\гкг 3 1 (кг) = кг Пп (кг);
V 2 -(п+2)
I — цилиндрическая функции Бесселя первого рода; п п (кг)— сферические функции Бесселя вто-
рого рода; jn (kr), n n (kr) — линейно независимые решения однородного уравнения (28).
Тогда решение (26) выражается через решения (27) в виде:
u(r) = k- £inA (2n + 1)Xn (r) Pn (cos 0), (29) kr n=0 (29)
r е R3
где
An =-
1 + ik J m( r') j-„ (kr ')x„ (r ')dr
Асимптотика поля (29) имеет вид [17]
1
u(r) = — £in(2n + 1)eiSn sin(kr - — + 8n)Pn(cos0),
kr n=0 2
откуда следует выражение для амплитуды рассеяния [4, 17]
/(0) = 2k£(2n + 1)(e2'5" - 1)Pn(cos0). (30)
2ik n=0
Здесь 8n определяется выражением
a
tg(5n) = -kJm(r)j-„ (kr)x„ (r)dr . (31)
0
Таким образом, после решения счетной системы (27) и нахождения из последних выражений постоянных S„ может быть рассчитана амплитуда рассеяния из (30). Очевидно, что замена ряда (30) конечной суммой первых N +1 членов эквивалентна решению интегрального уравнения (24) с усеченным аппроксимирующим ядром (25).
Еще больше упрощается задача расчета ар в случае радиально-симметричной неоднородности, сосредоточенной в ограниченном шаре B .
СЛУЧАЙ РАДИАЛЬНО-СИММЕТРИЧНОЙ ОГРАНИЧЕННОЙ НЕОДНОРОДНОСТИ
Пусть задан сферически симметричный ограниченный рассеиватель, заключенный в шаре B радиусом a с центром в начале координат. Решение задачи рассеяния u = ut + us внутри области
B, где ui = eikz — первичная падающая плоская волна, us — рассеянная волна, имеет применительно к акустике вид [13]:
1 "
u(r, 0) =—£il (2l + 1)Xi (k, r) Pi (cos 0), (32)
kr l =0
где
— функ-
1
Xt" + [k 2(1 - m(r)) - Щ^-Х = 0,
r
x (0) = 0, r е [0, a] или с учетом (7а):
c 2 l (l + 1)
(33)
X + [ k 2 ) 2 C (x) r
Xi(0) = 0, r е [0,a].
]Xi = 0,
(34)
вид
u = — Jil (2l + 1)x
kr l=0
X
1
j (kr) + (kr)
Pl (cos0), (35)
где h1 (1) (kr ) = kr h/1 (kr); hl^(kr) — сферическая функция Ханкеля; a — неизвестные постоянные. При всех a = 0 выражение (35) представляет собой падающую плоскую волну.
Из (35) следует выражение для ар [13]:
f (0)=2k J(2l+1)a Pl (cos 0). (36)
2ik l=0
Коэффициенты al определяются из соотношения [13]:
a, = -2
Llj (ka)- kajl (ka) L1fi()(ka)- kahf )(ka)
(37)
где Ll =
Г
dln Xi dlnr
— логарифмическая производ-
ная на границе г = а; производные в (37) берутся по аргументу кг . Логарифмические производные Ц легко рассчитываются при помощи решений X . Например, решением (34) при постоянном показателе преломления п(г) = п = с2/ с12 является функция
СЛУЧАЙ ОДНОРОДНОГО ШАРА
Наиболее простым для анализа является случай, когда в качестве объемного рассеивателя выбран шар В, радиусом а с центром в начале координат. Плотность шара совпадает с плотностью окружающей среды, скорость звука постоянна:
m(x) =
Вне шара в решение задачи рассеяния имеет
|1 -c2/c^| = M = const, xе B;
|0, x g в.
I {к,г) = Л^ ПТ '+ЛЬг),
где к1 = к—; Л1 — некоторые постоянные коэф-
с1
фициенты. При этом для нахождения величин Ц отпадает необходимость вычисления постоянных Л1, которые в этом случае просто сократятся. Отметим, что имеет место соотношение [13]
а = (въ5п -1).
Рассмотрим ряд (36). Учитывая, что |Р(х)| < 1, х е [-1,1], а также то, что при I >> 1 справедливо 21 +1 ,
отношение —3 ~ 1, приходим к выводу, что
скорость сходимости ряда (36) определяется поведением коэффициентов а . Для случая однородной сферы при различных радиусах а поведение коэффициентов а представлено в таблице. Скорость звука в окружающей среде с1 = 1500 м / с , частота / = 1000 Гц . Как видно, при а = 1м последовательность а мажорируется геометрической прогрессией со знаменателем 10-1 уже начиная с 10 = 5 ; при а = 2м — начиная с 10 = 13 и т. д. Поэтому погрешность оценки ряда (36) конечной суммой 10 слагаемых (I = 0,1,..., 10 -1), не превышает 1°а0 . Реально с ростом а число слагаемых
в конечной сумме ряда можно принимать значительно меньшим соответствующего значения 10 (см. таблицу).
Ниже на рис. 6-10 представлены кривые модуля ар в зависимости от угла падения при а, равных соответственно 1м, 2м,..., 5м. На рис. 11 они представлены в одних координатах. Как видно, величина первого минимума кривых уменьшается от ~ ±70° при а = 1 м до ~ ±10° при а = 5 м. Амплитуда максимума при этом возрастает соответственно от 0.679 до 68.8.
Отметим, что близкие результаты получены автором при вычислении ар из системы одномерных интегральных уравнений (27) при п = 4 для случая, относящегося к рис. 6 (/ = 1000 Гц , а = 1,
с = 1500м/с , с1 = 1500м/с ). Уравнения решались методом квадратур с использованием квадратурной формулы Симпсона [15]. Отметим, однако, что скорость вычислений при это уменьшилась практически на порядок по сравнению со случаем использования формул (36), (37).
Коэффициенты а1 в случае объемного рассеивателя в виде однородной сферы при различных радиусах а
1 а = 1 м а = 2 м а = 3 м а = 4 м а = 5 м
0 0.455 1 1.459 1.732 1.914
1 0.549 0.971 1.360 1.732 1.946
2 0.407 0.998 1.433 1.705 1.909
3 0.163 0.842 1.341 1.730 1.936
4 0.039 0.895 1.327 1.653 1.903
5 0.006 0.821 1.340 1.704 1.905
6 0.0006 0.533 1.172 1.620 1.903
13 1.1-10"13 110-5 0.047 1.07 1.595
14 2.28-10-15 9.19-10-7 0.012 0.723 1.571
20 1.2-10"26 3.28-10-14 1.52-10-7 1.510-3 0.249
21 1.1-10"28 1.27-10-15 1.48-10-8 3.2-10-4 0.109
27 1.2-10"41 7-10-25 1.62-10-15 2.3-10-9 3.7-10-5
28 6.3-10-44 1.5-10-26 8.4-10-17 2.3-10-10 6.56-10-6
35 1.9-10-60 5.5-10-39 1.2-10-26 2.8-10-18 3.5-10-12
36 3.8-10-63 212-10-41 3.6-10-22 1.61019 3.3-10-13
ар(0)
0.6
0.5
р. 4
/ 0.3
/ 0.2
/ 0.1
-150 -100 -50
ар(0)
Рис. 6. Зависимость модуля ар от угла 0 при
ка = 4/3 п, а = 1
50 100 150
0, град
-150 -100 -50
50 100 150
0, град
Рис. 7. Зависимость модуля ар от угла 0 при ка = 8/3 п, а = 2
ар(0)
31
....... , .. „ —j
-150 -100 -50
'и I',
■ V
0, град
V-...............
50 100 150
0, град
Рис. 8. Зависимость модуля ар от угла 0 при
ka = 4 п, a = 3
Рис. 9. Зависимость модуля ар от угла 0 при ka = 16/3 п, a = 4
ар(0)
70
60
50
40
30
20
J
-150 -100 -50 50 ICO 150
0, град
Рис. 10. Зависимость модуля ар от угла 0 при ka = 20/3 п, a = 5
ар(0)
70 ь
-150 -100 -50 50 100 150
0, град
Рис. 11. Зависимость модуля ар от угла 0 для всех случаев, показанных на рис. 6-10
АМПЛИТУДА РАССЕЯНИЯ ПОВЕРХНОСТНЫХ РАССЕИВАТЕЛЕЙ
Пусть в И3 в поле плоской падающей волны и. = вЯк,х находится поверхностный рассеиватель, занимающий область Б с гладкой замкнутой границей дБ. Рассмотрим две задачи: на границе рассеивателя для суммарного поля падающей
и рассеянной волн u(x) = ui + us выполняется 1) однородное условие Дирихле
u(x) = 0, xе dD (38)
или
2) однородное условие Неймана du(x)
dn
= 0, x edD.
(39)
где п — вектор внешней нормали к поверхности дБ.
Тогда суммарное поле можно определить соответственно из выражений [12]:
u(x) = u(x) - JV(y)G(x - y)dsy
dD
(40)
для условия Дирихле и
и(х) = и, (х) + Ыу)-тв(х - у)^ (41)
дБ дп
для условия Неймана. В (40) и (41) х е Я3 \ Б ; п' в (41) — вектор внешней нормали к дБ в точке у ;
1 exp(ik|x - y I)
G(x -У) = ---1|
4n |x - y|
(42)
Для определения функций ^(x) и (p(x) из (40), (41) следуют интегральные уравнения [12]
^(x) + J^(y G( x - y )ds
2 dD dn
= — ut (x), x e dD; on
1 Г d
-P(x) -J p(y) G(x - y )ds
(43)
dD
dn
= ui (x), x e dD.
(44)
Из (40), (41) следуют обычным порядком выражения для амплитуд рассеяния для обоих рассматриваемых случаев
Га(к,,к,) = -4-^(у)ехрНк, (45)
T
-L ЛТ
(k,,k,) = -!- J (k,n')p(y) exp(-ik,y)dsy . (46) 4n J
Здесь по-прежнему учтено, что имеют место зависимости
W(x) = ^(x, k,), p(x) = p(x, k,).
(47)
В работе [12] приведено решение интегрального уравнения для шара с однородными условиями Дирихле и Неймана. В общем случае интегральные уравнения (43), (44) рассматриваются в работе [14].
плитуды рассеяния объемных и поверхностных рассеивателей, принятые в теории рассеяния волн и частиц. Приведены точные и приближенные методы вычислений. Указаны погрешности использования приближенных методов. Приведены примеры расчетов и некоторые вычислительные тонкости.
СПИСОК ЛИТЕРАТУРЫ
1. Шарфарец Б.П. Метод решения задач множественного рассеяния на нескольких телах в однородной безграничной среде // Акуст. журн. 2005. Т. 51, № 5. С. 672-681.
2. Бреховских Л.М. Волны в слоистых средах. М.: Наука, 1973. 344 с.
3. Colton D., Kress R. Inverse acoustic and electromagnetic scattering theory. Springer, New York,
1998.
4. Шарфарец Б. П. Поле направленного излучателя в слоисто-неоднородном волноводе // Акуст. журн. 1985. Т. 31, № 1. С. 119-125.
5. Hansen T.B. and Yaghjian A.D. Plane-wave theory of time-domain fields. IEEE Press, New York,
1999.
6. Морс Ф.М., Фешбах Г. Методы теоретической физики. Т. 2. М.: Изд-во иностр. лит-ры, 1960. 860 с.
7. Забрейко П.П., Кошелев А.И. и др. Интегральные уравнения. М.: Наука, 1968. 448 c.
8. Владимиров В. С. Уравнения математической физики. М.: Наука, 1981. 512 c.
9. Иванов В.В. Методы вычислений на ЭВМ. Справочное пособие. Киев: Наукова думка, 1986. 584 с.
10. Краснов М.Л. Интегральные уравнения. М.: Наука, 1975. 304 с.
11. Михлин С.Г., Смолицкий Х.Л. Приближенные методы решения дифференциальных и интегральных уравнений. М.: Наука, 1965. 384 с.
12. Хёнл Х., Мауэ А., Вестпфаль К. Теория дифракции. М.: Мир, 1964. 428 с.
13. Флюгге З. Задачи по квантовой механике. Т. 1. М.: Мир, 1974. 341 с.
14. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М.: Мир, 1987. 311 с.
15. Полянин А.Д., Манжиров А.В. Справочник по интегральным уравнениям. М.: Физматлит, 2003. 608 с.
ВЫВОДЫ
В настоящей работе применительно к акустике рассмотрены различные методы вычисления ам-
Санкт-Петербург
Материал поступил в редакцию 26.06.2006.
ON CALCULATION OF SCATTERING AMPLITUDES OF BULK AND SURFACE SCATTERERS
B. P. Sharfarets
Saint-Petersburg
The paper describes some approaches to calculation of acoustic scatterers' scattering amplitudes, which are used in various fields of wave and particle scattering theory. Both bulk and surface scatterers are considered. Specific examples of calculations obtained by various methods are presented.