Уравнение Лапласа.

Поскольку решение U(x, y) уравнения Лапласа не зависит от времени, то можно в это уравнение добавить равный нулю (при точном решении) член . Тогда уравнение примет вид

Это — известное нам уравнение теплопроводности, для которого мы уже строили разностные схемы. Остается только задать на­чальное условие. Его можно принять практически в произвольном виде, согласованном с граничными условиями. Положим

Граничное условие при этом остается стационарным, т. е. не зави­сящим от времени.

Процесс численного решения такого уравнения состоит в переходе при от произвольного значения к искомому стационарному решению. Счет ведется до выхода решения на стационарный режим. Естественно, ограничиваются решением при некотором достаточно большом t, если искомые значения на двух после­довательных слоях совпадают с заданной степенью точности.

Метод установления фактически представляет итерационный процесс, причем на каждой итера­ции значения искомой функции получаются путем численного решения

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

Другой способ решения задачи Дирихле состоит в построении разностной схемы путем аппроксимации уравнения Лапласа. Вве­дем в прямоугольной области G сетку с по­мощью координатных прямых х = const и у = const. Примем, для простоты значения шагов по переменным, х и у равными h (пред­полагается, что стороны области G соизме­римы). Значения функции U в узлах заменим значениями сеточной функции Тогда, аппроксимируя в урав­нении Лапласа вторые производные с помощью отношений конечных раз­ностей, получим разностное уравнение.

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

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

Перейдем теперь к решению полученной системы. Каждое уравнение системы (за исключе­нием тех, которые соответствуют узлам, расположенным вблизи границ) содержит пять неизвестных. Одним из наиболее распространенных мето­дов решения этой системы линейных уравнений является итерационный метод. Каждое из уравнений записываем в виде, разрешенном относи­тельно значения в центральном узле:

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

Такой способ можно использовать и для решения уравнения Лапласа. Пусть требуется решить для него задачу Дирихле в прямоугольнике ABCD. Разо­бьем прямоугольник на полосы с помощью прямых, параллельных оси х. Для определенности проведем три отрезка , которые разделят прямоугольник на четыре полосы постоянной ширины h. Решение U задачи Дирихле приближенно заменим набором функций , каждая из которых определена на отрезке U и зависит только от одной переменной х, т. е. = для =1, 2, 3. На отрезках значения заданы граничными условиями.

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

Таким образом, решение задачи Дирихле сводятся к ре­шению краевой задачи для системы обыкновенных дифференциальных уравнений относительно значений искомой функции вдоль пря­мых . В этом состоит метод прямых. Граничные условия при х=а, х = b можно получить из уравнений

Метод прямых широко, используется для решения нестационарных задач. Например, если имеются две независимые переменные х, t, а ис­комый параметр является гладкой функцией переменной х, то дискре­тизация вводится по этой переменной. Тогда исходная задача заменяется задачей Коши для системы обыкновенных дифференциальных уравнений вида