Урок 6. Разделение переменных в декартовых координатах

advertisement
1. Электростатика
1.
1
Электростатика
Урок 6
Разделение переменных в декартовых координатах
1.1. (Задача 1.49) Плоскость z = 0 заряжена с плотностью σ (x, y) = σ0 sin (αx) · sin (βy),
где σ0 , α, β – постоянные. Найти потенциал этой системы зарядов.
Решение Потенциал ϕ удовлетворяет уравнению Лапласа
∂ 2ϕ ∂ 2ϕ ∂ 2ϕ
+ 2 + 2 = 0.
∂x2
∂y
∂z
На заряженной плоскости нормальная составляющая электрического поля терпит разрыв
¯
¯
E2z ¯z=0 − E1z ¯z=0 = 4πσ.
(1)
(2)
Поскольку поверхностная плотность σ есть произведение двух функций, одна из которых зависит
только от координаты x, а другая от y, будем искать решение уравнения (1) в виде
ϕ(x, y, z) = f1 (x) · f2 (y) · f3 (z).
Подставим функцию (3) в уравнение (1) и поделим каждое слагаемое на ϕ, получим
µ 00
¶
f300 (z)
f1 (x) f200 (y)
=−
+
.
f3 (z)
f1 (x)
f2 (y)
(3)
(4)
Чтобы уравнение (4) удовлетворялось при всех x, y, z, нужно приравнять каждое слагаемое константе. Пусть
f300 (z)
= −γ 2 .
f3 (z)
Общим решением этого уравнения является функция
f3 (z) = Ae−γz + B eγz .
Понятно, что при стремлении z → ±∞ , f3 (z) не должна стремиться к бесконечности, поскольку поле создается знакопеременными зарядами. Этому условию можно удовлетворить, если
записать решение, например, так:
f3 (z) = A e−γ|z| .
Пусть
f100 (x)
= −γ12 ,
f1 (x)
f200 (y))
= −γ22 , γ 2 = γ12 + γ22 .
f2 (y)
2
Решениями этих уравнений являются гармонические функции. Чтобы удовлетворить гармоническому условию (2), нужно выбрать:
значит, γ =
p
f1 (x) = sin αx ,
α2 = γ12 ;
f2 (y) = sin βy ,
γ22 = β 2 ;
α2 + β 2 . Итак,
ϕ = A e−γ|z| sin αx sin βy.
Нормальная составляющая электрического поля Ez = −∂ϕ/∂z. Поэтому
¯
E2z ¯z=0 = γA sin αx sin βy,
¯
E1z ¯z=0 = −γA sin αx sin βy.
Подставляя эти выражения в уравнение (2), найдем, что A = 2σ0 /γ. Окончательно распределение потенциала будет иметь вид
ϕ(x, y, z) =
где γ =
p
1.2.
Y
b
Ψ=0
Ψ=0
Ψ=0
2σ0 −γ|z|
e
sin αx sin βy,
γ
α2 + β 2 .
Металлическая труба, бесконечная вдоль оси z, имеет прямоугольное
сечение a × b. Одна из сторон имеет потенциал U , а остальные
- потенциал 0. Найти распределение потенциала внутри труΨ=U
бы.
a X
Решение Расположим сечение z = 0 по отношению к осям как показано на рисунке. Поскольку труба бесконечна вдоль z, то зависимость от z отсутствует и мы имеем задачу с двумя
переменными. Поскольку внутри трубы зарядов нет, потенциал внутри удовлетворяет уравнению
Лапласа
∂2ψ ∂2ψ
4ψ =
+ 2 =0
(1)
∂x2
∂y
и граничным условиям ψ(0, y) = ψ(x, 0) = ψ(x, b), ψ(a, y) = U. Предположим, что решение
уравнения (1) (или, хотя бы частное его решение) можно представить в виде
ψ(x, y) = X(x) · Y (y).
(2)
Тогда, разделив уравнение (1) на X(x)Y (y), его можно переписать в виде
1 d2 X
1 d2 Y
+
= 0.
X dx2
Y dy 2
(3)
1. Электростатика
3
Перенеся одно из слагаемых в правую часть уравнения, получим, что функция от x равна во всей
области определения функции от y, что возможно только в случае когда обе эти функции – константы, причем
1 d2 X
1 d2 Y
=
−
= λ2 .
(4)
X dx2
Y dy 2
Можно было бы просто в правой части записать константу и не определять ее знак. Впрочем приведенная запись не гарантирует положительность правой части (например, если λ – мнимая). На
самом деле далее будет показано, что только такой выбор знака для нашей задачи позволяет получить нетривиальное решение. Очевидно, что вместо 1 уравнения с частными производными мы
получили систему двух уравнений, каждое из которых - обыкновенное дифференциальное уравнение.
Y 00 + λ2 Y = 0,
X 00 − λ2 X = 0.
(5)
Общее решение первого уравнения имеет вид
Y (y) = A sin(λy) + B cos(λy).
Граничные уравнения для этого уравнения получаются из общих граничных условий
ψ(x, 0) = ψ(x, b) = 0, ⇒ X(x)Y (0) = 0, X(x)Y (b) = 0,
откуда получаем требуемые два условия
Y (0) = 0, Y (b) = 0.
(6)
Из уравнения для Y и приведенных выше граничных условий следует, что если бы в первом уравнении системы (5) знак у константы был противоположным, то полученные в этом случае экспоненциальные решения не смогли бы одновременно удовлетворить граничным условиям (6). Используя граничные условия можно получить следующий набор (бесконечный) частных решений,
удовлетворяющих и уравнению и граничным условиям
³ nπ ´
nπ
y , λn =
, n = 1, 2, ...
Yn (y) = An sin
b
b
Подставив значение λn во второе уравнение системы (5), получим решение этого уравнения
X(x) = A sh
nπ
nπ
x + B ch
x.
b
b
4
Из граничного условия на границе x = 0 следует, что X(0) = 0, откуда получается, что частные
решения этого уравнения, удовлетворяющее условию на левой границе имеют вид
Xn = Bn sh
nπ
x.
b
(7)
Используя полученные частные решения для X и Y , а также принцип суперпозиции, можно записать общее решение всей задачи (удовлетворяющее только 3 граничным условиям) в виде
ϕ(x, y) =
∞
X
An sh
n=1
nπ
nπ
x sin
y.
b
b
(8)
Теперь можно использовать четвертое граничное условие (ϕ(a, y) = U ) для определения коэффициентов An . Это граничное условие переписывается в виде
ϕ(a, y) =
∞
X
An sh
n=1
nπa
nπ
sin
y = U.
b
b
(9)
Для определения коэффициентов этого ряда можно воспользоваться формулами для определения
коэффициентов ряда Фурье, поскольку именно ряд Фурье по y записан в выражении для граничного условия. А можно явно домножить правую и левую часть уравнения (9) на sin mπ
y ипроb
интегрировать обе стороны полученного равенства по y от 0 до b/ Поскольку на этом интервале
функции sin mπ
y для разных n и m ортогональны, то в сумме останется только член с n = m.
b
mπa
Am sh
b
Zb
mπy
sin2
dy = U
b
0
Zb
sin
mπy
dy.
b
0
Интеграл слева вычисляется следующим образом
Zb
mπy
b
sin2
dy =
b
mπ
Zmπ
b mπ
b
sin2 tdt =
= .
mπ 2
2
0
0
Интеграл справа вычисляется элементарно
Zb
0
mπy
b
sin
dy = −
b
mπ
Zmπ
b
d(cos t) = −
(cos mπ − 1)
mπ
0
(10)
1. Электростатика
5
Из этого выражения ясно, что отличны от нуля коэффициенты только при нечетном m, т.е. при
m = 2k + 1. Окончательно получаем
A2k+1 =
4U
(2k + 1)π sh (2k+1)πa
b
,
или, окончательно получаем
∞
(2k+1)πx
4U X sh
(2k + 1)πy
b
ϕ(x, y) =
sin
.
(2k+1)πa
π k=0 sh
b
b
таким образом, использованный метод разделения переменных позволил нам получить общее решение о потенциале внутри прямоугольника в общем виде. При этом следует сделать два замечания
1. Если необходимо решить задачу с разными значениями потенциала на всех четырех границах, то это делается на основании принципа суперпозиции, используя полученное выше
решение.
2. Ряд сходится довольно плохо, особенно вблизи правой границы.
Для демонстрации этого была написана программа для расчета численного примера. Результаты
расчета показаны на графиках. Видно (на третьем графике) что даже при использовании 20 членов
ряда значение потенциала на границе колеблется, т.е. граничное условие выполняется приближенно.
6
Изолинии для V0=1,a=5, b=2, N=10
2
y
1.5
1
0.5
0
2.5
3
3.5
1
V(x,y=1)
x
4
4.5
5
0.8
V
0.6
0.4
0.2
0
2
2.5
3
3.5
4
4.5
5
x
V(x=a,y=1)
1.2
1.1
1
0.9
0.8
0
2
4
6
8
10
12
n - число членов ряда
14
16
18
20
1.3. Найти потенциал между двумя бесконечными плоскостями, если потенциал плоскости x =
0 равен A sin βy, а потенциал плоскости x = d равен ϕ0 + B sin βy.
Решение Сформулируем математическую задачу: найти решение уравнения Лапласа в области
0≤x≤d
4ϕ(x, y) = 0,
с граничными условиями ϕ(x = 0, y) = A sin βy, ϕ(x = d, y) = ϕ0 + B sin βy. Эта задача
не сводится очевидным образом к методу разделения переменных, поэтому, используя принцип
суперпозиции, разобъем эту задачу на две, точнее, представим решение задачи в виде суммы двух
решений ϕ(x, y) = ϕ(x, y)1 + ϕ(x, y)2 ..
Система I
4ϕ1 = 0,
ϕ1 (x = 0, y) = 0,
ϕ1 (x = d, y) = ϕ0 ,
Система II
4ϕ2 = 0,
ϕ2 (x = 0, y) = A sin βy,
ϕ2 (x = d, y) = B sin βy,
Решаем сначала уравнение I. Поскольку граничные условия не зависят от y, можно предположить,
что решение этого уравнения будет зависеть только от x. Тогда уравнение сводится к обыкновен-
1. Электростатика
7
ному дифференциальному уравнению
d2 ϕ1
= 0, при условиях ϕ1 (0) = 0, ϕ1 (d) = ϕ0 .
dx2
Откуда получаем
ϕ0
x.
d
Теперь решим второе уравнение методом разделения переменных. Как обычно, предположим, что
искомое решение можно представить в виде ϕ2 (x, y) = X(x)Y (y) и подставим это в исходное
уравнение для ϕ2 (x, y). Тогда, разделив обе части уравнения на X(x)Y (y), получим
ϕ1 (x) =
X 00 (x)
Y 00 (y)
=−
= γ 2.
X(x)
Y (y)
Каждая из дробей равна константе, поскольку если функция только от x равна во всей области
определения функции от y? то такая функция может быть только константой. Знак константы
выбран из соображений удобства и из тех (эвристических) соображений, что функция от y должна быть скорее всего ∼ sin(y). Используя сделанные предположения, получим частное решение
Y (y) = sin γy, откуда делаем вывод, что γ = β. Тогда уравнение для X(x) имеет решение
X(x) = C1 exp (−γx) + C2 exp (γx).
Используя приведенные выше граничные условия для решения задачи II, запишем систему уравнений для констант C1 и C2
¡
(C1 + C2 ) sin βy =A sin βy,
¢
C1 e−γd + C2 eγd sin βy =B sin βy.
Сокращая левую и правую части этих уравнений на sin βy, получим систему уравнений
(C1 + C2 ) =A,
¡
¢
C1 e−γd + C2 eγd =B,
решение которой имеет вид
Aeγd − B
C1 = γd
,
e − e−γd
Ae−γd + B
C2 = γd
e − e−γd .
Тогда окончательное решение можно записать в виде
µ γd
¶
ϕ0
Ae − B −γx Ae−γd + B γx
ϕ(x, y) = ϕ(x, y)1 + ϕ(x, y)2 =
x+
e
+ γd
e
sin βy.
d
eγd − e−γd
e − e−γd
8
1.4. Рассмотрим плоскую задачу, т.е. потенциал не зависит от z. Область определения представляет собой полубесконечную полосу 0 ≤ x ≤ ∞, 0 ≤ y ≤ b. Потенциал Ψ(x, y) удовлетворяет уравнению Лапласа, Ψ(0, y) = V, Ψ(x, 0) = Ψ(x, b) = 0.
Решение Уравнение Лапласа в декартовой системе координат записывается в виде
∂ 2Ψ ∂ 2Ψ
+
= 0.
∂x2
∂y 2
Предположим, как обычно, что решение представимо в виде
Ψ(x, y) = X(x) · Y (y).
Подставив это в исходное уравнение и разделив уравнение на функцию Ψ(x, y), получим
1 d2 X
1 d2 Y
=
−
= γ 2.
2
2
X dx
Y dy
Поскольку функция x (в левой части уравнения) равна во всей области определения функции y (в
правой части уравнения) эти комбинации функций должны равняться константе. Рассмотрим эти
получившиеся уравнения по отдельности.
d2 Y
+ γ 2 Y = 0.
dy 2
Общее решение этого уравнения имеет вид
Y (y) = A sin γy + B sin γy.
Из граничных условий для потенциала можно получить условия для функции Y
Ψ(x, 0) = Ψ(x, b) = 0, ⇒ Y (0) = Y (b) = 0.
Из этих условий следует, что B = 0 и γk = kπ
, k = 1, 2, ... Подставляя полученные допустиb
мые значения для γ в уравнение для X, получим
d2 X
+ γk2 X = 0.
dx2
Общее решение этого уравнения, удовлетворяющее условию ограниченности решения, имеет вид
Xk = Ce−γk x ,
1. Электростатика
9
а общее решение исходной задачи записывается в виде
Ψ(x, y) =
∞
X
ak e−
kπx
b
sin
k=1
kπy
.
b
Для определения коэффициентов ak используем граничное условие на левой границе
Ψ(0, y) =
∞
X
ak sin
k=1
kπy
= V.
b
Для нахождения коэффициентов ak можно воспользоваться общими формулами рядов Фурье,
а можно вычислить их явно, для чего домножим левую и правую части граничного условия на
sin mπy
и проинтегрируем левую и правую части по y. В силу ортогональности функций sin mπy
b
b
и sin kπy
при
m
=
6
k
в
сумме
останется
только
один
член
с
m
=
k,
а
уравнение
для
вычисления
b
коэффициента будет иметь вид
Zb
mπy
dy = V
sin
b
Zb
2
am
sin
mπy
dy,
b
0
0
или, преобразуя переменные под интегралом
mπ
R
am = V
0
mπ
R
sin tdt
=
sin2 tdt
4V
, при m = 2n + 1, n = 0, 1, ...
mπ
0
окончательное выражение приобретает вид
¤
£
∞
πy
4V X sin (2n + 1) b − (2n+1)πx
b
Ψ(x, y) =
e
.
π n=0
2n + 1
Самое удивительное, что этот ряд можно свернуть в конечное выражение. Для этого перепишем
полученное выражение в виде
( ∞ £ πx
)
¤
X e− b · ei πyb 2n+1
4V
Ψ(x, y) =
Im
.
π
2n
+
1
n=0
Введя обозначение Z = e−
πx
b
· ei
πy
b
и используя соотношение
∞
X
1 1+Z
Z 2n+1
= ln
,
2n + 1
2 1−Z
n=0
10
получаем в результате
·
¸
sin πy
2V
b
Ψ(x, y) =
arctg
.
π
sh πx
b
Download