Задача пуассона для уравнения лапласа

Уравнения Лапласа и Пуассона

(1.110)

называется уравнением Пуассонав трехмерном пространстве. Если в этом уравнении , то оно называется уравнением Лапласа:

.

(1.111)

Если ввести оператор , называемый оператором Лапласа, то уравнения (1.110) и (1.111) запишутся соответственно

и .

К исследованию уравнений Лапласа и Пуассона приводит рассмотрение задач о стационарном процессе: это задачи гидродинамики, диффузии, фильтрации, распределения температуры, электростатики и др.

Эти уравнения относятся к уравнениям эллиптического типа.

Те задачи, которые приводят к уравнениям, содержащим время, называются динамическими или нестационарными задачами математической физики; задачи, приводящие к уравнениям, не содержащим время, называются стационарнымиили статическими.

О постановке задачи математической физики

И ее корректности

Как было показано, уравнения математической физики имеют бесчисленное множество решений, зависящее от двух произвольных функций (речь идет об уравнениях второго порядка для функции двух переменных). Для того, чтобы из множества решений выделить определенное, характеризующее процесс, необходимо на искомую функцию наложить дополнительные условия, которые диктуются физическими соображениями. Тут можно провести аналогию с обыкновенными дифференциальными уравнениями, когда для выделения из общего решения частного, удовлетворяющего некоторым дополнительным условиям, отыскивались по этим условиям произвольные постоянные. Таковыми условиями для уравнений в частных производных являются, чаще всего, начальные и граничные условия. Граничные условия – это условия, заданные на границе рассматриваемой среды; начальные условия – условия, относящиеся к какому-нибудь моменту времени, с которого начинается изучение данного физического явления. Дополнительные условия,

так же как и само дифференциальное уравнение, должны вводиться на основе физических соображений, связанных с самим процессом. Вместе с тем дополнительные условия должны быть такими, чтобы обеспечить выделение из всего множества решений единственного решения. Число граничных и начальных условий определяется типом уравнения, а их вид – заданным исходным состоянием на границе объекта и внешней среды. Для рассматриваемых нами уравнений число начальных условий равно порядку старшей производной по времени, входящей в уравнение, а число граничных условий – порядку старшей производной по координате.

Совокупность дифференциального уравнения и дополнительных условий представляет собой математическую формулировку физической задачи и называется задачей математической физики.

Физическая задача решается по схеме:

1) реальный физический процесс (явление, объект) заменяется некоторым идеальным процессом (явлением, объектом) так, что последний значительно проще первого и вместе с тем сохраняет его основные черты (идеализация процесса);

2) выбирается величина (функция), характеризующая процесс, и используются законы, по которым он происходит;

3) на основании выбранных законов выводится дифференциальное уравнение для величины, характеризующей процесс;

4) выводятся дополнительные условия – начальные и граничные – также в соответствии с выбранными законами.

Итак, задача математической физики состоит в отыскании решений уравнений в частных производных, удовлетворяющих некоторым дополнительным условиям, скажем, граничным и начальным.

Задача математической физики считается поставленной корректно, если решение задачи, удовлетворяющее всем ее условиям, существует, единственно и устойчиво; последнее означает, что малые изменения любого из данных задачи вызывают малое изменение решения. Требование устойчивости необходимо по следующей причине. В данных любой конкретной задачи, особенно если они получены из опыта, всегда содержится некоторая погрешность, и нужно, чтобы малая погрешность в исходных данных приводила к малой неточности в решении. Это требование выражает физическую определенность поставленной задачи.

Примеры

ПРИМЕР 2.36. Выяснить, являются ли приведенные ниже равенства дифференциальными уравнениями в частных производных:

а) ,

б) .

Решение. Преобразуем уравнение а)

.

Данное уравнение является уравнением в частных производных, так как в него входят частные производные второго порядка

и .

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

ПРИМЕР 2.37. Выяснить, какие из следующих уравнений являются линейными (однородными или неоднородными) и какие нелинейными:

а) ,

б) ,

в) .

Решение. Сравнивая данные уравнения с формой (1.4), заключаем, что

— уравнение а) есть неоднородное линейное уравнение второго порядка, для которого ;

— уравнение б) нелинейное, так как оно не является линейным относительно старших частных производных;

— уравнение в) является однородным линейным уравнением третьего порядка.

ПРИМЕР 2.38. Решить уравнение .

Решение. Ясно, что искомая функция не зависит от переменной , но может быть любой функцией от : , поскольку, дифференцируя по , получим ноль, а это значит, что данное равенство выполняется. Таким образом, решение уравнения содержит одну произвольную функцию .

ПРИМЕР 2.39. Решить уравнение , где заданная функция.

Решение. Интегрируя по , восстановим искомую функцию

, где произвольная функция.

Итак, решение уравнений в примерах 2.38 и 2.39 содержат одну произвольную функцию . Такое решение называется общим. В отличие от общего решения обыкновенного дифференциального уравнения первого порядка, которое содержит одну произвольную постоянную, решение уравнения в частных производных первого порядка содержит одну произвольную функцию.

ПРИМЕР 2.40. Решить уравнение .

Решение. Перепишем уравнение так: . Положим , после чего данное уравнение принимает вид . Как было установлено в примере 2.38, общее решение последнего уравнения имеет вид: , где произвольная функция. Исходное уравнение примет вид: . Проинтегрировав полученный результат по , получим

, иначе

,

где и произвольные дважды дифференцируемые функции.

Легко проверить, что найденная функция удовлетворяет данному уравнению.

Итак, решение уравнения в частных производных второго порядка содержит уже две произвольные функции. Такое решение называют общим.

Приведенные в качестве примеров уравнения дают основание сделать заключение: общее решение уравнения в частных производных первого порядка содержит одну произвольную функцию, а общее решение уравнения второго порядка – две произвольные функции. В этом заключается коренное отличие общего решения уравнения в частных производных от общего решения обыкновенного дифференциального уравнения, которое содержит одну и две произвольные постоянные.

В дальнейшем будет выяснено, какие дополнительные условия надо задать, чтобы с их помощью можно было выделить частное решение, т. е. функцию, удовлетворяющую как уравнению, так и дополнительным условиям.

2.17 КЛАССИФИКАЦИЯ И ПРИВЕДЕНИЕ К КАНОНИЧЕСКОМУ ВИДУ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ

ВТОРОГО ПОРЯДКА

Все многообразие линейных относительно старших производных (или просто линейных) уравнений может быть разделено на три класса (типа). В каждом классе есть простейшие уравнения, которые называются каноническими. Решения уравнения одного и того же типа (класса) имеют много общих свойств. Для изучения этих свойств достаточно рассмотреть канонические уравнения, так как другие уравнения данного класса могут быть приведены к каноническому виду.

Классификация уравнений вида (2.52) проводится в соответствии со знаком дискриминанта .

Говорят, что уравнение (1.3) в области принадлежит

а) гиперболическому типу, если ,

б) параболическому типу, если ,

в) эллиптическому типу, если .

(2.54)

называется каноническим уравнением гиперболического типа.

Второй канонический вид уравнения гиперболического типа таков:

(2.55)

(2.56)

называется каноническим уравнением параболического типа.

(2.57)

называется каноническим уравнением эллиптического типа.

Дата добавления: 2017-10-09 ; просмотров: 1430 ; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ

Уравнения Пуассона и Лапласа

Уравнения Пуассона и Лапласа являются основными дифференциальными уравнениями электростатики. Они вытекают из теоремы Гаусса в дифференциальной форме. Действительно, подставляя в уравнение

вместо величин Е_х; Е_у; Е_z их выражения через потенциал:

Это дифференциальное уравнение носит название уравнения Пуассона.

является решением уравнения Пуассона для случая, когда заряды распределены в конечной области пространства.

Если в рассматриваемой области пространства отсутствуют объемные электрические заряды, то уравнение Пуассона получает вид

и называется в этом частном случае уравнением Лапласа.

Отметим, что в цилиндрической и сферической системах координат уравнение Пуассона и Лапласа имеют другую форму записи. Поэтому данные уравнения часто записывают в виде, не зависящем от системы координат:

Оператор ? 2 часто обозначают и называют оператором Лапласа или лапласианом.

При интегрировании уравнения Лапласа (или Пуассона) в решение входят постоянные интегрирования. Их определяют из граничных условий.

Численное решение уравнений в частных производных эллиптического типа на примере уравнений Лапласа и Пуассона

Среди всех типов уравнений математической физики эллиптические уравнения с точки зрения вычислителей стоят особняком. С одной стороны, имеется хорошо развитая теория решения эллиптических уравнений и систем. Достаточно легко доказываются теоремы об устойчивости разностных схем для эллиптических уравнений. Во многих случаях получаются априорные оценки точности расчетов и числа итераций при решении возникающих систем сеточных уравнений . С другой стороны, системы сеточных уравнений , возникающие при решении уравнений методами сеток, имеют большую размерность и плохо обусловлены. Для решения таких систем разработаны специальные итерационные методы .

6.1. Постановка задачи. Простейшая разностная схема «крест». Устойчивость схемы «крест»

Будем рассматривать двухмерное уравнение Пуассона

в единичном квадрате с краевыми условиями первого рода на границе расчетной области

( — заданная на границе функция ).

В случае прямоугольной области граничные условия удобно записать в следующем виде:

Для простоты выкладок введем равномерную расчетную сетку с узлами m, y_l> , m, l = 0, 1, . , M с равным количеством шагов по каждому пространственному направлению, сеточную область D — совокупность всех узлов сетки, включая граничные, и сеточную функцию < u_ml >. В этом случае шаги по координатам предполагаются равными. В случае неравных шагов по каждому направлению полученные результаты не изменятся, а запись уравнений станет более громоздкой.

Выбираем простейший пятиточечный шаблон разностной схемы «крест» . На этом шаблоне аппроксимирующее разностное уравнение легко выписать. Для этого производные заменим вторыми разностями:

где h — шаг по координатам, или в операторной форме

Эту же разностную схему можно записать в каноническом виде для разностных схем для эллиптических уравнений:

Такую каноническую запись не следует путать с канонической формой записи итерационного метода, которая встретится ниже.

Такая схема обладает вторым порядком аппроксимации по обеим координатам. Это легко показать, применяя разложение в ряд Тейлора функции — проекции точного решения на сетку — вплоть до членов четвертого порядка включительно. Проведем такое разложение для одного из операторов, стоящих в данном разностном уравнении:

Здесь учтено разложение проекции точного решения в ряд Тейлора

и аналогичное разложение для u_{m — 1}.

Для рассматриваемого двухмерного уравнения получим выражение для главного члена невязки

Рассмотрим устойчивость полученной схемы. Отметим, что методы исследования на устойчивость , применяемые для эволюционных (зависящих от времени) уравнений, здесь не работают. Действовать приходится на основе определения устойчивости.

Сформулируем и докажем две леммы, которые облегчат процедуру доказательства устойчивости разностной схемы.