Численные методы решения уравнений эллиптического типа
Введение
Наиболее распространённым уравнением эллиптического типа является уравнение Пуассона.
К решению этого уравнения сводятся многие задачи математической физики, например задачи о стационарном распределении температуры в твердом теле, задачи диффузии, задачи о распределении электростатического поля в непроводящей среде при наличии электрических зарядов и многие другие.
Для решения эллиптических уравнений в случае нескольких измерений используют численные методы, позволяющие преобразовать дифференциальные уравнения или их системы в системы алгебраических уравнений. Точность решения определяется шагом координатной сетки, количеством итераций и разрядной сеткой компьютера [1]
Цель публикации получить решение уравнения Пуассона для граничных условий Дирихле и Неймана, исследовать сходимость релаксационного метода решения на примерах.
Уравнение Пуассона относится к уравнениям эллиптического типа и в одномерном случае имеет вид [1]:
(1)
где x – координата; u(x) – искомая функция; A(x), f(x) – некоторые непрерывные функции координаты.
Решим одномерное уравнение Пуассона для случая А = 1, которое при этом принимает вид:
(2)
Зададим на отрезке [xmin, xmax] равномерную координатную сетку с шагом ∆х:
(3)
Граничные условия первого рода (условия Дирихле) для рассматриваемой задачи могут быть представлены в виде:
(4)
где х1, xn – координаты граничных точек области [xmin, xmax]; g1, g2 – некоторые
константы.
Граничные условия второго рода (условия Неймана) для рассматриваемой задачи могут быть представлены в виде:
(5)
Проводя дискретизацию граничных условий Дирихле на равномерной координатной сетке (3) с использованием метода конечных разностей, получим:
(6)
где u1, un – значения функции u(x) в точках x1, xn соответственно.
Проводя дискретизацию граничных условий Неймана на сетке (3), получим:
(7)
Проводя дискретизацию уравнения (2) для внутренних точек сетки, получим:
(8)
где ui, fi – значения функций u(x), f(x) в точке сетки с координатой xi.
Таким образом, в результате дискретизации получим систему линейных алгебраических уравнений размерностью n, содержащую n – 2 уравнения вида (8) для внутренних точек области и уравнения (6) и (7) для двух граничных точек [1].
Ниже приведен листинг на Python численного решения уравнения (2) с граничными условиями (4) – (5) на координатной сетке (3).
Разработанная мною на Python программа удобна для анализа граничных условий.Приведенный алгоритм решения на Python использует функцию Numpy — u=linalg.solve(a,b.T).T для решения системы алгебраических уравнений, что повышает быстродействие при квадратной матрице . Однако при росте числа измерений необходимо переходить к использованию трех диагональной матрицы решение для которой усложняется даже для очень простой задачи, вот нашёл на форуме такой пример:
Программа численного решения на равномерной по каждому направлению сетки задачи Дирихле для уравнения конвекции-диффузии
(9)
Используем аппроксимации центральными разностями для конвективного слагаемого и итерационный метод релаксации.для зависимость скорости сходимости от параметра релаксации при численном решении задачи с /(х) = 1 и 6(х) = 0,10. В сеточной задаче:
(10)
Представим матрицу А в виде суммы диагональной, нижней треугольной и верхней треугольных матриц:
(10)
Метод релаксации соответствует использованию итерационного метода:
(11)
При \ говорят о верхней релаксации, при — о нижней релаксации.
На графике показана зависимость числа итераций от параметра релаксации для уравнения Пуассона (b(х) = 0) и уравнения конвекции-диффузии (b(х) = 10). Для сеточного уравнения Пуассона оптимальное значении параметра релаксации находится аналитически, а итерационный метод сходиться при .
- Приведено решение эллиптической задачи на Python с гибкой системой установки граничных условий
- Показано что метод релаксации имеет оптимальный диапазон () параметра релаксации.
Ссылки:
- Рындин Е.А. Методы решения задач математической физики. – Таганрог:
Изд-во ТРТУ, 2003. – 120 с. - Вабищевич П.Н.Численные методы: Вычислительный практикум. — М.: Книжный дом
«ЛИБРОКОМ», 2010. — 320 с.
Уравнение Пуассона и математическая постановка задач электростатики
Существует большое количество случаев, когда самым удобным методом нахождения напряженности поля считается решение дифференциального уравнения для потенциала. После его получения применим в качестве основы теорему Остроградского-Гаусса в дифференциальной форме:
где ρ является плотностью распределения заряда, ε 0 — электрической постоянной, d i v E → = ∇ → E → = ∂ E x ∂ x + ∂ E y ∂ y + ∂ E z ∂ z — дивергенцией вектора напряженности и выражением, связывающим напряженность поля и потенциал.
Произведем подстановку ( 2 ) в ( 1 ) :
Учитывая, что d i v g r a d φ = ∇ 2 φ = ∂ 2 φ ∂ x 2 + ∂ 2 φ ∂ y 2 + ∂ 2 φ ∂ z 2 , где ∆ = ∇ 2 — это оператор Лапласа, равенство ( 3 ) принимает вид:
Выражение ( 4 ) получило название уравнения Пуассона для вакуума. При отсутствующих зарядах запишется как уравнение Лапласа:
После нахождения потенциала переходим к вычислению напряженности, используя ( 2 ) . Решения уравнения Пуассона должны удовлетворять требованиям:
- значение потенциала как непрерывная функция;
- потенциал должен быть конечной функцией;
- производные потенциала как функции по координатам должны быть конечными.
При наличии сосредоточенных зарядов в объеме V , решение уравнения ( 4 ) будет выражаться для потенциала вида:
Общая задача электростатики сводится к нахождению решения дифференциального уравнения, то есть уравнения Пуассона, удовлетворяющего вышеперечисленным требованиям. Теоретические вычисления известны для небольшого количества частных случаев. Если возможно подобрать функцию φ , удовлетворяющую условиям, то она является единственным решением.
В таких задачах не всегда необходимо задавать заряды или потенциалы во всем пространстве. Для нахождения электрического поля в полости, окруженной проводящей оболочкой, достаточно вычислить поле тел, находящихся внутри нее.
Любое решение уравнения Пуассона ограниченной области может быть определено краевыми условиями, накладывающимися на поведение решения. Границы перехода из одной среды в другую имеют условия, которые должны быть выполнены:
E 2 n — E 1 n = 4 π σ , или ∂ φ 1 ∂ n — ∂ φ 2 ∂ n = 0 .
где σ — это поверхностная полость свободных зарядов, n – единичный вектор нормали к границе раздела, проведенный из среды 1 в 2 , τ — единичный вектор, касательный к границе.
Эти уравнения выражают скачок нормальных составляющих вектора напряженности и непрерывность касательной вектора напряженностей электрического поля при переходе через любую заряженную поверхность независимо от ее формы и наличия или отсутствия зарядов вне ее.
Уравнение Пуассона в сферических, полярных и цилиндрических координатах
Запись уравнения может быть как при помощи декартовых координат, также и сферических, цилиндрических, полярных.
При наличии сферических r , θ , υ уравнение Пуассона запишется как:
1 r 2 · ∂ ∂ r r 2 ∂ φ ∂ r + 1 r 2 sin θ ∂ θ sin θ · ∂ φ ∂ θ + ∂ 2 φ r 2 sin 2 θ ∂ φ 2 = — 1 ε 0 ρ .
В полярных r , θ :
1 r · ∂ ∂ r r ∂ φ ∂ r + ∂ 2 φ r 2 ∂ θ 2 = — 1 ε 0 ρ .
В цилиндрических r , υ , z :
1 r · ∂ ∂ r r ∂ φ ∂ r + ∂ 2 φ ∂ z 2 + ∂ 2 φ r 2 ∂ υ 2 = — 1 ε 0 ρ .
Примеры решения задач
Найти поле между коаксиальными цилиндрами с радиусами r 1 и r 2 и с имеющейся разностью потенциалов ∆ U = φ 1 — φ 2 .
Решение
Необходимо зафиксировать уравнение Лапласа с цилиндрическими координатами, учитывая аксиальную симметрию:
1 r · ∂ ∂ r r ∂ φ ∂ r = 0 .
Решение имеет вид φ = — A ln ( r ) + B . Для этого следует выбрать нулевой потенциал на нужном цилиндре, тогда:
φ ( r 2 ) = 0 = — A ln r 2 + B , следовательно
φ ( r 1 ) = ∆ U = — A ln r 1 + B , получим:
A = ∆ U ln r 2 r 1 .
φ ( r ) = — ∆ U ln r 2 r 1 ln ( r ) + ∆ U ln r 2 r 1 ln r 2 .
Ответ: поле с двумя коаксиальными цилиндрами может быть задано при помощи функции φ ( r ) = — ∆ U ln r 2 r 1 ln ( r ) + ∆ U ln r 2 r 1 ln r 2 .
Найти потенциал поля, которое создает бесконечно круглый цилиндр с радиусом R и объемной плотностью заряда ρ . Использовать уравнение Пуассона.
Решение
Необходимо направить ось Z по оси цилиндра. Видно, что цилиндрическое распределение заряда аксиально симметрично, потенциал имеет такую же симметрию, иначе говоря, считается функцией φ ( r ) с r , являющимся расстоянием от оси цилиндра. Для решения используется цилиндрическая система координат. Уравнение Пуассона в ней запишется как:
φ 2 = C 2 ln r + C ‘ 2 .
C 1 , C ‘ 1 , C 2 , C ‘ 2 — это постоянные интегрирования. Имеем, что потенциал во всех точках должен быть конечным, а l i m r → 0 ln r = ∞ . Отсюда следует, что C 1 = 0 . Далее необходимо пронормировать потенциал, задействовав условие φ 1 ( 0 ) = 0 . Получим C ‘ 1 = 0 .
Поверхностные заряды отсутствуют, поэтому напряженность электрического поля на поверхности шара является непрерывной. Следовательно, что и производная от потенциала также непрерывна при r = R , как и сам потенциал. Исходя из условий, можно найти C 2 , C ‘ 2 :
C 2 ln R + C ‘ 2 = — 1 4 ρ ε 0 R 2 .
C 2 R = — 1 2 ρ ε 0 R .
Значит, полученные выражения записываются как:
Ответ: потенциал поля равняется:
Уравнение Пуассона
Определение и формула уравнения Пуассона
Уравнение Пуассона описывает адиабатический процесс, происходящий в идеальном газе. Адиабатический процесс — это процесс, в котором нет теплообмена между рассматриваемой системой и окружающей средой:
Уравнение Пуассона имеет вид:
Здесь V — объем, занимаемый газом, P — его давление, а значение k называется адиабатическим индексом.
Адиабатический индекс в уравнении Пуассона
Адиабатический индекс можно рассчитать как отношение изобарной теплоемкости газа к его изохорной теплоемкости:
В практических расчетах удобно помнить, что для идеального газа адиабатический индекс равен для двухатомного и для трехатомного .
Что относительно реальных газов, когда силы взаимодействия между молекулами начинают играть важную роль? В этом случае адиабатический индекс для каждого испытательного газа может быть получен экспериментально. Один из таких методов был предложен в 1819 году Климентом и Дезормом. Мы наполняем цилиндр холодным газом, пока давление в нем не достигнет Р1. Затем мы открываем клапан, газ начинает адиабатически расширяться, а давление в цилиндре падает до атмосферного ПА. После того, как изохорный газ нагрелся до температуры окружающей среды, давление в цилиндре повысится до P2. Тогда адиабатический индекс можно вычислить по формуле:
Адиабатический индекс всегда больше 1, поэтому при адиабатическом сжатии газа — как идеального, так и реального — температура газа всегда поднимается до меньшего объема, а при расширении газ охлаждается. Это свойство адиабатического процесса, называемого пневматическим кремнем, используется в дизельных двигателях, где горючая смесь сжимается в цилиндре и воспламеняется теплом. Напомним первый закон термодинамики: , где — внутренняя энергия системы, а А — выполненная на ней работа. Поскольку работа, выполняемая газом, идет только для изменения ее внутренней энергии — и, следовательно, температуры. Из уравнения Пуассона можно получить формулу для расчета газовой операции в адиабатическом процессе:
Здесь n — количество газа в молях, R — универсальная газовая постоянная, T — абсолютная температура газа.
Уравнение Пуассона для адиабатического процесса используется не только при расчетах двигателей внутреннего сгорания, но и при проектировании холодильных машин.
Стоит вспомнить, что уравнение Пуассона точно описывает только равновесный адиабатический процесс, состоящий из непрерывно меняющихся состояний равновесия. Если на самом деле мы открываем клапан в цилиндре так, чтобы газ расширялся адиабатически, то возникнет нестационарный переходный процесс с газовой турбулентностью, который будет испаряться из-за макроскопического трения.
Примеры решения проблем
Одноатомный идеальный газ был адиабатически сжат, так что его объем увеличился в 2 раза. Как изменится давление газа?
Адиабатический индекс для одноатомного газа равен . Однако его можно вычислить по формуле:
где R — универсальная газовая постоянная, а і — степень свободы молекулы газа. Для одноатомного газа степень свободы равна 3: это означает, что центр молекулы может выполнять поступательное движение вдоль трех координатных осей.
Поэтому адиабатический индекс:
Представьте себе состояние газа в начале и конце адиабатического процесса через уравнение Пуассона:
Давление уменьшится в 3.175 раз.
100 молей двухатомного идеального газа было адиабатически сжато при 300 К. В то же время давление газа увеличилось в 3 раза. Как изменился газ?
Степень свободы двухатомной молекулы равна i = 5, так как молекула может двигаться постепенно вдоль трех координатных осей и вращаться вокруг двух осей.
Рассчитайте диатомический адиабатический индекс:
Определите, как изменяется объем газа при адиабатическом сжатии, из уравнения Пуассона:
Это означает, что объем газа уменьшился в 2,19 раза.
Вычислите работу газа, используя следующую формулу:
http://zaochnik.com/spravochnik/fizika/elektricheskoe-pole/uravnenie-puassona/
http://www.homework.ru/spravochnik/uravnenie-puassona/