urmat
.pdf
41
|
∫α |
(M ) |
∫ |
|
||
В этом случае u(P) = |
|
ϕ(M ) |
G(M , P)dσM − |
|
G(M , P) f (M )dvM . |
|
S |
2 |
|
V |
|||
|
|
|
|
|||
Таким образом, исходная задача сводится к задаче о нахождении функции Грина.
16. Нахождение функции Грина методом электростатических изображений.
Функция Грина находится явно лишь для областей частного вида. При построении функции Грина полезно воспользоваться следующей ее физической интерпретацией. Из курса физики известно, что электрический заряд величины q, помещенный в точку Р, создает в свободном неограниченном пространстве электростатическое поле, потенциал которого (при определенном выборе системы единиц)
равен: |
u0 (M ) = |
|
q |
. |
(16.1) |
|
|
|
|||||
|
|
4πρMP |
|
|
||
Поскольку в замкнутой области V функция Грина G(M,P) имеет вид |
||||||
|
G(M , P) = |
1 |
+ g(M , P), |
(16.2) |
||
|
|
|
||||
|
|
|
4πρMP |
|
||
где g(M,P), как функция точки М, гармоническая в области V и непрерывна вместе с первыми производными в замыкании [V] области V, то первое слагаемое в правой части (16.2) является потенциалом точечного единичного заряда, помещенного в точке Р области V. Второе слагаемое g(M,P) можно также интерпретировать, как потенциал электростатического поля, созданного одним или несколькими зарядами, но расположенными обязательно вне области V. Это возможно сделать потому, что потенциал электростатического поля, т.е. функция вида (16.1), является гармонической функцией в любой области, свободной от зарядов (т.е. при М≠Р). Заряды вне V надо выбрать так, чтобы они уничтожили на поверхности S действие заряда в точке P, т.е. чтобы выполнялось соотношение G S = 0;эти заряды называют электроста-
тическими изображениями единичного заряда в Р, а сам метод нахождения функции Грина − методом электростатических изображений.
Пример 1. Функция Грина для полупространства. Пусть S есть плоскость z=0, а V − полупространство z>0. Если в точке P V поместить единичный положительный заряд, то его действие на S уничтожится, очевидно, единичным отрицательным
зарядом, помещенным в точке P1,которая является зеркальным изображением точки P относительно S. Потенциал поля, созданного зарядом в P1 , есть
g(M , P) = − |
|
|
1 |
|
и, следовательно, функция Грина − |
|||||
4πρMP |
||||||||||
|
|
|
|
|
||||||
|
|
|
|
1 |
|
|
|
|
||
|
1 |
|
1 |
|
|
1 |
|
|||
G(M , P) = |
|
|
− |
|
||||||
|
|
|
|
|
|
|||||
|
4π |
|
|
ρMP |
|
ρMP1 |
|
|||
|
|
|
|
|
||||||
42
Пример 2. Функция Грина для шара. Пусть V − шар, ограниченный сферой S:
x2 + y2 + z2 = R2 с центром в начале координат. Поместим единичный заряд в точку P, расположенную внутри сферы S. Покажем, что действие этого заряда на S может быть уничтожено некоторым зарядом, помещенным в точке P1, являющейся инверсией точки P относительно сферы S; точка P1,лежит на прямой ОР вне шара, причем
ρOP ρOP = R2 . |
(16.3) |
1 |
|
Пусть P − произвольно зафиксированная точка сферы S. Рассмотрим два треугольника OPP и OP1 P. Эти треугольники подобны, так как они имеют общий угол при вершине О и стороны, образующие этот угол, пропорциональны в силу (16.3). Из
|
|
|
|
|
ρ |
|
|
|
|
= |
|
ρ0 |
|
|
ρ0 = ρOP ; откуда |
||||
подобия треугольников следует |
PP |
, где |
|||||||||||||||||
|
ρ |
|
|
|
|
R |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
PP |
|
|
|
|
|
|||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||
1 |
|
|
− |
|
R |
|
|
|
1 |
|
|
|
= 0 |
(16.4) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
4πρ |
|
|
|
|
ρ0 |
4πρ |
|
|
|
|
||||||||
PP |
PP1 |
|
|||||||||||||||||
при любом положении точки P на сфере.
Рисунок 7.
Из (16.4) следует, что действие заряда q =1 в Р уничтожается на S зарядом
q = − |
R |
, |
помещенном в P . Следовательно, функция Грина − потенциал поля, |
|
|||
|
ρ0 |
1 |
|
|
|
||
|
1 |
|
1 |
|
R |
|
1 |
|
|
созданного этими зарядами − есть G(M , P) = |
|
− |
|
|
|||||
|
|
|
|
||||||
|
|
|
ρMP |
|
|
|
|
|
|
|
4π |
|
ρ0 ρMP1 |
||||||
тростатических изображений можно применить и в плоском случае, однако физическая интерпретация здесь будет несколько иной. Именно, если на прямой, проходя-
щей через точку Р ортогонально плоскости (x,y), разместить положительные электрические заряды с единичной плотностью, то они создадут плоское поле (т.е. поле, не зависящее от координаты z), потенциал которого (при соответствующем выборе
системы единиц) будет равен u0 (M ) = |
1 |
ln |
1 |
. |
2π |
|
|||
|
|
ρMP |
||
43
17. Решение задачи Дирихле для шара.
Зная функцию Грина, можно построить решение задачи Дирихле для уравнения Лапласа. Поскольку в этом случае f (M ) = 0, то искомое решение примет вид
|
∫ |
|
|
|
|
|
|
|
|
∂n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
u(P) = − |
|
|
ϕ(M ) |
|
∂G |
dσM . Поскольку в этом случае |
|
|||||||||||||||||||
S |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
1 |
|
|
1 |
|
|
R |
|
|
|
1 |
|
|
|
|
|
|
|||||
G(M , P) = |
|
|
|
|
|
|
|
|
|
|
|
, то |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
ρMP |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
4π |
ρ0 ρMP1 |
|
|
|
|
|||||||||||||||||
|
|
|
1 |
|
∫ |
|
∂ |
|
|
|
|
1 |
|
|
R |
|
|
1 |
|
|
|
|
||||
u(P) = − |
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
. Таким образом, |
||||||||||
|
4π |
|
|
|
|
|
|
ρMP |
ρ0 |
|
|
|
|
ϕ(M )dSM , где ϕ(M ) = u |
S |
|||||||||||
|
|
S |
∂n |
|
|
ρMP1 |
|
|||||||||||||||||||
остается произвести дифференцирование в подынтегральном выражении. Пусть
Т− переменная точка, расположенная внутри шара.
Рисунок 8.
1
Из треугольников ОТР и OTP1 получим ρPT = (ρ2 + ρ02 − 2ρρ0 cosγ )2 ,
1
ρTP1 = (ρ2 + ρ12 − 2ρρ1 cosγ )2 ,
где ρ = ρOT , ρ0 = ρOP , ρ1 = ρOP1 , γ − угол при вершине О. Учитывая, что направление внешней нормали к сфере совпадает с направлением радиуса (т.е. с направлением роста ρ ), получим
|
|
∂ |
|
|
|
1 |
|
|
|
R |
|
|
1 |
|
|
∂ |
|
|
1 |
|
|
R |
|
1 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
. Учитывая выражения для |
|
|
∂n |
|
ρ |
|
|
ρ |
|
ρ |
|
∂P |
ρ |
|
|
ρ |
|
ρ |
|
|||||||||||||||||||
|
|
|
P M |
|
|
|
0 |
|
|
|
|
|
TP |
|
0 |
|
|
|
|
|
|
|||||||||||||||
|
|
|
M |
|
|
|
|
|
|
|
|
MP1 |
|
|
|
|
|
|
|
|
|
|
|
|
TP1 |
|
|
ρ=R |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R2 |
|
|
|
|
|
|
|
|
|
|||||||||||
|
ρ |
TP |
и для ρ |
TP |
, а также, что ρ |
|
= |
подсчитаем правую часть последнего выра- |
||||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
ρ0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
жения и получим, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
u(P) = |
|
1 |
|
∫ |
|
|
|
|
R2 − ρ02 |
|
|
|
|
|
|
|
|
ϕ(M )dSM . Эта формула называется фор- |
||||||||||||||||||
4πR |
|
|
|
|
|
− 2Rρ |
|
|
|
|
|
|
|
3 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
cosγ )2 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
S (R2 + ρ2 |
0 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мулой Пуассона, а интеграл, стоящий справа − интегралом Пуассона.
44
Можно проверить, что функция u(P) действительно является решением задачи Дирихле для шара при любой непрерывной функции ϕ(M ). Для этого следует записать подынтегральное выражение в декартовых координатах
|
R2 − x |
2 − y2 |
− z2 |
|
|
|
ϕ(M ), |
|||
|
|
0 |
0 |
0 |
|
|
|
|
||
|
|
|
|
|
|
|
3 |
|
||
((x |
−ξ)2 + ( y |
0 |
|
−η)2 |
+ (z |
0 |
−ς)2 ) |
2 |
|
|
0 |
|
|
|
|
|
|
|
|
||
где (x0 , y0 , z0 ), (ξ,η,ς) −соответственно координаты точек P(x0 , y0 , z0 ) и
M (ξ,η,ς).
Непосредственным дифференцированием убеждаемся, что это выражение как функция точки P(x0 , y0 , z0 ) удовлетворяет уравнению Лапласа.
45
Литература
1.Государственный образовательный стандарт высшего профессионального образования. − М.: Госкомвуз России, 2000.
2.Петровский И.Г. Лекции об уравнениях с частными производными. М.: Физмат-
гиз, 1961.
3.Смирнов В.И. Курс высшей математики, тт. 2 и 4. − М.: Физматгиз, 1958.
4.Арсенин В.Я. Математическая физика. Основные уравнения и специальные функции. − М.: Наука, 1961.
5.Владимиров В.С. Уравнения математической физики. − М.: Наука, 1971.
6.Бицадзе А.В. Уравнения математической физики. − М.: Наука, 1976.
7.Михлин С.Г. Линейные уравнения в частных производных. − М.: Высшая школа, 1977.
8.Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Уравнения в частных производных математической физики. − М.: Высшая школа, 1970.
9.Араманович И.Г., Левин В.И. Уравнения математической физики. − М.: Наука, 1969.
10.Кальницкий Л.А., Добротин Д.А., Жевержеев В.Ф. Специальный курс высшей
математики для втузов. − М.: Высшая школа, 1976.
11.Смирнов М.М. Дифференциальные уравнения в частных производных второго порядка. − М.: Наука, 1964.
12.Лебедев Н.Н., Скольская И.П., Уфлянд Я.С. Сборник задач по математической
физике. − М.: Гостехиздат, 1955.
13.Смирнов М.М. Задачи по уравнениям математической физики. − М.: Наука, 1968.
14.Чудесенко В.Ф. Сборник заданий по специальным курсам высшей математики. − М.: Высшая школа, 1983.
15.Данко П.Е., Попов А.Г. Высшая математика в упражнениях и задачах. Ч.3. − М.: Высшая школа, 1971.
