Линейные ДУ в частных производных второго порядка. Основные уравнения математической физики. Уравнения Максвелла, страница 12

Пусть дана эллиптическая краевая задача вида LU = f в области Д (1), [¡1¶U/¶n + ¡2U]|S = j(м) (2). Для разрешимости этой задачи необходимо: ¡1, ¡2 ³ 0, ¡12 + ¡22 ¹ 0. Для того чтобы найти решение (1), (2) нужно решить задачу со специальной правой частью. LG(M,M0) = – d(M,M0) (3) в Д. [¡1¶G/¶n + ¡2G]|S = 0 (4). Потребуем чтобы G(М,М0) была непрерывна всюду вместе с частными производными 1го порядка в замкнутой области Д. Если функция G найдена, то тогда можно решить задачу (1), (2). Применим вторую формулу Грина: ò{D}[G(M,M0)LU – ULG(M,M0)]dt = ò{S}[G(M,M0)¶U/¶n – U¶G(M,M0)/¶n]dS (5). LU = f(M) по условию, LG = – d(M,M0), тогда ò{S}G(M,M0)f(M)dtм + ò{D}U(M)d(M,M0)dtм = ò{S}[G(M,M0)¶U/¶n – U¶G(M,M0)/¶n]dS. На основании свойств d-функции ò{D}U(M)d(M,M0)dtм = U(M0). Тогда U(M0) = –ò{D}G(M,M0)f(M)dtм + ò{S}[G(M,M0)¶U/¶n – U¶G(M,M0)/¶n]dS (6).

1. Решение с помощью формулы (6) задачи Дирихле.

В (2) ¡1 = 0, ¡2 ¹ 0, U|S = j(M), G(M,M0)|S = 0. U(M0) = –ò{D}G(M,M0)f(M)dtм – ò{S}j(M)¶G(M,M0)/¶n dS.

2. Решение с помощью формулы (6) задачи Неймана.

В (2) ¡1 ¹ 0, ¡2 = 0, ¶U/¶n|S = j1(M), ¶G(M,M0)/¶n|S = 0. U(M0) = –ò{D}G(M,M0)f(M)dtм + ò{S}G(M,M0)j1(M)dS.

3. Третья краевая задача

¡1 ¹ 0, ¡2 ¹ 0. [¡1¶G/¶n + ¡2G]|S = 0, [¡1¶U/¶n + ¡2U]|S = j(M). ¶G/¶n = – (¡2/¡1)G|S, ¶U/¶n = – (¡2/¡1)U|S + j(M)/¡1; U(M0) = –ò{D}G(M,M0)f(M)dtм + ò{S}[G(M,M0)[ – (¡2/¡1)U + j(M)/¡1] + (¡2/¡1)GU]dS => U(M0) = – ò{D}G(M,M0)f(M)dtм + ò{S}j(M)G(M,M0)/¡1 dS.

Опр.: Точки М и М1 называются инверсными относительно плоскости в пространстве или прямой в плоскости, если они симметричны относительно этой плоскости или прямой.

М(х0, у0, z0) -> М(х0, у0, –z0) {3D}, М(х0, у0) -> М(х0, –у0) {2D}.

Опр.: Точки М и М1 называются сопряжёнными относительно сферы, если они лежат на одном плече, исходящем из центра, а произведение их расстояний от центра равно квадрату радиуса.

Рис.

R2 = r*r1

Лекция №10 за 3.11.2000

Теория потенциала.

Рис.

r = [(x-a)2 – (y-b)2 – (z-c)2]1/2. Пусть в точку А помещаем электрический заряд, тогда заряд создаёт электростатическое поле. Напряжённость поля в любой точке пространства, не совпадающей с точкой А, определяется следующим образом:`Е = q`r/r3. Это поле имеет проекции: Ех = q(x-a)/r3, Еy = q(y-b)/r3, Еz = q(z-c)/r3 (1). Правые части соотношения (1) с противоположным знаком равны частным производным функции U(м) = q/r + const (2) – потенциал электростатического поля. U(м) –> 0 при r –> 0, поэтому const = 0. U(м) = q/r = q/[(x-a)2 – (y-b)2 – (z-c)2]1/2. (3)  Предположим, что у нас несколько точечных источников, тогда потенциал U = åUi. Если заряд распределён по некоторому объёму Д, т.е. задана плотность заряда, тогда U(м) = ò{D}r(N)/r dtN (4). Правая часть (4) – объёмный потенциал. Заряд распределён по поверхности, его плотность r1(N), тогда его потенциал U(м) = ò{S}r1(N)/r dS (5). r – расстояние от точки М до переменной точки поверхности S. (5) – потенциал простого слоя.

Рис.

А – центр, h << расстояния до точки М. Рассмотрим величину qh. Обозначим qh = p = const. U(м) = lim{h->0}q(1/r’ – 1/r”) = lim{h->0}p(1/r’ – 1/r”)/h = p*¶(1/r)/¶L = p*cos(`АМ,`L)/r2 (6). При помощи двух таких зарядов диполь может быть представлен очень приближённо.

Рис.

S – строго ориентированная поверхность (это поверхность, на которой строго определены внешняя и внутренняя нормали), ni – внутренняя нормаль, nL – внешняя нормаль. На поверхности S распределён диполь с плотность m(N). В каждой точке N направление оси диполя совпадает с направлением внутренней нормали поверхности S. Тогда потенциал, который создаёт этот диполь будет равен: W(M) = ò{S}m(N)*cos(`NM,`ni)/r2 dS (7), где вектор r направлен от N к М. (7) – интеграл двойного слоя, т.к. распределение диполя может быть приближённо представлено как два, наложенных на поверхность S распределения зарядов с плотностью m(N)/n и –m(N)/n. В дальнейшем будем полагать, что r направлено от точки М к N, а нормаль будем брать внешнюю. W(M) = –ò{S}m(N)cosj/r2 dS.