![]() Главная страница Случайная страница КАТЕГОРИИ: АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника |
Уравнение Лапласа. ⇐ ПредыдущаяСтр 9 из 9
Многие стационарные физические задачи (исследования потенциальных течений жидкости, определение формы нагруженной мембраны, задачи теплопроводности и диффузии в стационарных случаях и др.) сводятся к решению уравнения Пуассонавида При F(x, у, z) = 0, уравнение Пуассона называют уравнением Лапласа. Для простоты будем рассматривать двумерное уравнение Лапласа Решение этого уравнения будем искать для некоторой ограниченной области G изменения независимых переменных х, у. Границей области G является замкнутая линия L. Для полной формулировки краевой задачи кроме уравнения Лапласа нужно задать граничное условие на границе L. Примем его в виде
Задача, состоящая в решении уравнения Лапласа (или Пуассона) при заданных значениях искомой функции на границе расчетной области, называется задачей Дирихле. Одним из способов решения стационарных эллиптических задач, в том числе и краевой задачи, является их сведение к решению некоторой фиктивной нестационарной задачи (гиперболической или параболической), найденное решение которой при достаточно больших значениях времени t близко к решению исходной задачи. Такой способ решения называется методом установления. Поскольку решение U(x, y) уравнения Лапласа не зависит от времени, то можно в это уравнение добавить равный нулю (при точном решении) член Это — известное нам уравнение теплопроводности, для которого мы уже строили разностные схемы. Остается только задать начальное условие. Его можно принять практически в произвольном виде, согласованном с граничными условиями. Положим Граничное условие при этом остается стационарным, т. е. не зависящим от времени.
Метод установления фактически представляет итерационный процесс, причем на каждой итерации значения искомой функции получаются путем численного решения некоторой вспомогательной задачи. В теории разностных схем показано, что этот итерационный процесс сходится к решению исходной задачи, если такое стационарное решение существует. Другой способ решения задачи Дирихле состоит в построении разностной схемы путем аппроксимации уравнения Лапласа. Введем в прямоугольной области G сетку с помощью координатных прямых х = const и у = const. Примем, для простоты значения шагов по переменным, х и у равными h (предполагается, что стороны области G соизмеримы). Значения функции U в узлах С помощью данного уравнения можно записать систему линейных алгебраических уравнений относительно значений сеточной функции в узлах в виде Значения сеточной функции в узлах, расположенных на границе расчетной области, могут быть найдены из граничного условия: Перейдем теперь к решению полученной системы. Каждое уравнение системы (за исключением тех, которые соответствуют узлам, расположенным вблизи границ) содержит пять неизвестных. Одним из наиболее распространенных методов решения этой системы линейных уравнений является итерационный метод. Каждое из уравнений записываем в виде, разрешенном относительно значения В ряде случаев уравнение с частными производными удобно привести к системе обыкновенных дифференциальных уравнений, в которых оставлены производные искомой функции лишь по одной переменной. Такой способ можно использовать и для решения уравнения Лапласа. Пусть требуется решить для него задачу Дирихле в прямоугольнике ABCD. Разобьем прямоугольник на полосы с помощью прямых, параллельных Построим разностную схему, для определения значений функций Таким образом, решение задачи Дирихле сводятся к решению краевой задачи для системы обыкновенных дифференциальных уравнений относительно значений искомой функции вдоль прямых Метод прямых широко, используется для решения нестационарных задач. Например, если имеются две независимые переменные х, t, а искомый параметр является гладкой функцией переменной х, то дискретизация вводится по этой переменной. Тогда исходная задача заменяется задачей Коши для системы обыкновенных дифференциальных уравнений вида
|