
Решение уравнения Пуассона
В достаточно большом количестве случаев наиболее удобным методом поиска напряженности поля является решение дифференциального уравнения для потенциала. Получим его, используя в качестве основы теорему Остроградского - Гаусса в дифференциальной форме:
где ρ -- плотность распределения заряда, ε0 -- электрическая постоянная, div→E=→∇→E=∂Ex∂x+∂Ey∂y+∂Ez∂z) -- дивергенция вектора напряженности и выражение связывающее напряженность поля и потенциал:
Подставим (2) в (1), получим:
Учитываем, что divgradφ=∇2φ=∂2φ∂x2+∂2φ∂y2+∂2φ∂z2, где △=∇2- оператор Лапласа, тогда равенство (3) запишем как:
Уравнение (4) называется уравнением Пуассона (для вакуума) в системе СИ. Если заряды отсутствуют, то уравнение (4) преобразуется в уравнение Лапласа:
После того, как найден потенциал из уравнения Пуассона, обычно вычисляется напряженность по формуле (2). Решения уравнения Пуассона должны удовлетворять таким требованиям:
- Потенциал должен быть непрерывной функцией.
- Потенциал должен быть конечной функцией.
- Производные от потенциала как функции по координатам должны быть конечными.
Если заряды сосредоточены в объеме V, то решением уравнения (4) будет выражение для потенциала вида:
Итак, общая задача электростатики сводится к нахождению решения дифференциального уравнения (уравнения Пуассона), которое удовлетворяет выше перечисленным требованиям. Нахождение решения -- задача весьма сложная. Теоретические решения известны для небольшого количества частных случаев. Если удалось подобрать функцию φ, которая удовлетворяет всем условиям задачи, то она единственная.
Не всегда есть необходимости задавать заряды или потенциалы во всем пространстве. Например, если необходимо найти электрическое поле в полости, которая окружена проводящей оболочкой, то можно найти поле только для тел внутри самой полости.
Каждое решение уравнения Пуассона в ограниченной области однозначно определяется краевыми условиями, которые накладывают на поведение решения. На границе перехода из одной среды в другую выполняются граничные условия:
E2n−E1n=4πσ, или ∂φ1∂n−∂φ2∂n=0.где σ- поверхностная плотность свободных зарядов, n- единичный вектор нормали к границе раздела, проведенный из среды 1 в 2, τ− единичный вектор, касательный к границе.
Данные уравнения выражают скачок нормальных составляющих вектора напряженности и непрерывность касательной составляющей вектора напряженностей электрического поля при переходе через любую заряженную поверхность не зависимо от формы этой поверхности и наличия и отсутствия зарядов вне ее.
Уравнение Пуассона в сферических, полярных и цилиндрических координатах
Уравнение Пуассона может быть записано не только в декартовых координатах, но также в сферических и цилиндрических, полярных.
В сферических координатах (r,θ,ϑ) уравнение Пуассона имеет следующий вид:
1r2⋅∂∂r(r2∂φ∂r)+1r2sinθ∂θ(sinθ⋅∂φ∂θ)+∂2φr2sin2θ∂φ2=−1ε0ρ(7).В полярных координатах (r,θ) система координат уравнение имеет вид:
1r⋅∂∂r(r∂φ∂r)+∂2φr2∂θ2=−1ε0ρ(8).В цилиндрических координатах (r,ϑ,z) уравнение имеет вид:
1r⋅∂∂r(r∂φ∂r)+∂2φ∂z2+∂2φr2∂ϑ2=−1ε0ρ(9).Задание: Найдите поле между двумя коаксиальными цилиндрами с радиусами r1 и r2, разность потенциалов между которыми равна △U=φ1−φ2.
Рис. 1
Решение:
Запишем уравнение Лапласа в цилиндрических координатах с учетом аксиальной симметрии:
1r⋅∂∂r(r∂φ∂r)=0(1.1)
Оно имеет решение φ=−Aln(r)+B. Выберем нулевой потенциал на наружном цилиндре, найдем, получим:
φ(r2)=0=−Alnr2+B, следовательно
B=Alnr2.В результате имеем: φ(r)=−△Uln(r2r1) ln(r)+△Uln(r2r1) lnr2
Ответ: Поле между двумя коаксиальными цилиндрами задается функцией φ(r)=−△Uln(r2r1) ln(r)+△Uln(r2r1) lnr2
Задание: Найти, используя уравнение Пуассона потенциал поля, которое создает бесконечно длинный круглый цилиндр радиуса R с объемной плотностью заряда ρ.
Решение:
Ось Z направим по оси цилиндра. Так как цилиндрическое распределение заряда аксиально симметрично, то потенциал обладает той же симметрией, то есть он является функцией φ(r), где r -- расстояние от оси цилиндра. Поэтому используем цилиндрическую систему координат. Запишем уравнение Пуассона в ней с учетом симметрии:
\frac{1}{r} \cdot \frac{\partial}{\partial r} \left(r\frac{\partial {\varphi }_1}{\partial r}\right)=- \frac{1}{ \varepsilon_0} \rho \ \left(0 Общие решения уравнений (2.1) (2.2) имеют вид:\[{\varphi }_1=-\frac{1}{4}\frac{\rho }{{\varepsilon }_0}r^2+C_1lnr+{C'}_1\left(2.3\right),
где C1, C′1,C2, C′2 -- постоянные интегрирования. Так как потенциал во всех точках должен быть конечным, а limr→0lnr =∞, следовательно, C1=0. Пронормируем потенциал условием:φ1(0)=0, тогда: C′1=0.
Так как поверхностных зарядов нет, то напряженность электрического поля на поверхности шара непрерывна, то есть непрерывна производная от потенциала при r=R. И непрерывен сам потенциал. Эти условия дают два алгебраических уравнения для того, чтобы найти постоянные C2, C′2:
C2lnR+C′2=−14ρε0R2(2.5)Следовательно, получаем выражения для потенциалов:
\[{\varphi }_1\left(r\right)=-\frac{1}{4}\frac{\rho }{{\varepsilon }_0} r ^ 2 \left(0 Ответ: Потенциал поля равен: $$ \varphi_1 \left(r \right)=-\frac{1}{4}\frac{\rho }{\varepsilon_0} r^2 \left (0