Алгоритм поиска приближенных решений уравнения Пуассона | Статья в журнале «Молодой ученый»

Автор:

Рубрика: Математика

Опубликовано в Молодой учёный №3 (83) февраль-1 2015 г.

Дата публикации: 22.01.2015

Статья просмотрена: 257 раз

Библиографическое описание:

Шевченко А. С. Алгоритм поиска приближенных решений уравнения Пуассона // Молодой ученый. — 2015. — №3. — С. 18-23. — URL https://moluch.ru/archive/83/15123/ (дата обращения: 24.10.2018).

Численное решение уравнения Пуассона является важным элементом многих задач вычислительной физики.

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

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

 

Рассмотрим в прямоугольной области  (см. рис. 2) уравнение Пуассона для функции :

                                                                       (1)

со смешанными краевыми условиями:

                                                                                                         (2)

Здесь  ‒ известная правая часть, ,  ‒ часть границы  ,  ‒ единичный вектор внешней нормали к границе ,  ‒ некоторые постоянные.

Рис. 2. Область и граница  для модельной задачи

 

Замечание. Сделав в задаче (1), (2) замену зависимой переменной :

,

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

Рис. 3. График вспомогательной достаточно гладкой функции

 

Следовательно, не нарушая общности, будем полагать далее, что в (2) функция  удовлетворяет на  однородному условию .

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

1.параболическая регуляризация

,                                                                                                          (3)

2. регуляризация Соболева

,                                                                                                 (4)

3.гиперболическая регуляризация

.                                                                                                (5)

Здесь  ‒ новая неизвестная функция,  ‒ постоянная, ‒ временная переменная. Краевые условия (2) для уравнений (3 ‒ 5) перепишутся так

.                                                                             (6)

Добавляя начальные условия

,                                                                                            (7)

стационарное решение исходной смешанной краевой задачи ищется как предел  при .

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

Для нахождения приближенного решения смешанной задачи (3), (6), (7) или (4), (6), (7) или (5), (6), (7) предложена дифференциально-разностная модель, основанная на методе прямых и принципах схем без насыщения. При этом производную  в уравнениях (3) — (5) мы заменим не разностным отношением, а используем для аппроксимации ее интерполяционный многочлен [1]

с узлами в нулях многочлена Чебышева

В итоге для каждой из регуляризаций будем иметь:

,                                                                                                       (8)

,                                                                                  (9)

,                                                                                            (10)

где , , , , , , ,   .

Далее в системах (8) ‒ (10) проведем дискретизацию по переменной . С этой целью вводим обозначения: , , , , ‒ шаг разностной сетки по времени.

Заменяя ,  получаем систему ОДУ второго порядка:

                                                                                                               (11)

Выражения для элементов матрицы  и компонент вектора  зависят от вида примененной регуляризации:

,

,

 ‒ единичная матрица порядка N.

К системе (11) добавляем граничные краевые условия (см. (6)):

                                                                             (12)

-диагональные матрицы порядка N(диагональные элементы равны либо 0, либо 1), причем:

Учитывая определенные трудности, возникающие при практическом использовании формулы для нахождения приближенного решения краевой задачи (11), (12), мы предложим следующую технологию построения приближенного решения этой краевой задачи. Будем искать приближенное решение в виде интерполяционного кубического сплайна класса (см. [2]):

                                   (13)

     

С учетом (13) из (11), (12) получаем трехточечную разностную схему:

                                        (11)

Здесь   ‒ диагональная матрица порядка N (диагональные элементы равны либо 0, либо 1). Систему алгебраических уравнений (11) можно решить методом матричной прогонки. Из неравенства  следуют условия хорошей обусловленности разностной краевой задачи (11).

Данный алгоритм был запрограммирован на языке Object Pascal (среда Delphi 6). Были проведены тестовые расчеты, которые показали хорошую работоспособность, приведенного выше алгоритма.

Данный алгоритм хорошо зарекомендовал себя при поиске решения 1D задачи о баллистическом диоде и 2D задачах о переносе заряда в транзисторе MESFET и MOSFET.

 

Литература:

 

1.                  Бабенко К. И. Основы численного анализа. Москва-Ижевск: НИЦ " Регулярная и хаотическая динамика», 2002.

2.                  Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л. Методы сплайн ‒ функций. М.: Наука, 1980.

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


Похожие статьи

Априорная оценка для решения первой краевой задачи для...

Для решения краевой задачи (1) — (2), мы используем приближенной (численный) метод (метод разностных схем). Для доказательства устойчивости разностной модели, мы используем априорную оценку (3).

Численная реализация разностного метода решения одной...

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

Разностная краевая задача для уравнения смешанного типа

Краевая задача: Найти в области решение уравнение (1) условие: (2). Численное решение краевой задачи (1)-(2) является непростой задачей ввиду

Рахмонов Х. О. О первой краевой задаче для одного уравнения смешанного типа в пространстве. — Новосибирск, 1985. -22с.

Решение методом продолжения задач математической физики...

Для дифференциальных уравнений различают три типа задач: задача Коши; краевая задача; смешанная задача.

Используя это решение, мы можем решить задачу Коши. Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия

Организация вычислений решения краевой задачи для...

1. Краевая задача для ОДУ. Сведения о приближённых методах.

В проекционных методах [1,2] приближённое решение отыскивается в виде конечной суммы с неопределёнными коэффициентами

Аппроксимация первой краевой задачи разностной моделью для...

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

Модульный анализ сеточных методов решения...

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

Решение трехмерного уравнения Пуассона с использованием...

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

Априорная оценка для решения первой краевой задачи для...

Для решения краевой задачи (1) — (2), мы используем приближенной (численный) метод (метод разностных схем). Для доказательства устойчивости разностной модели, мы используем априорную оценку (3).

Численная реализация разностного метода решения одной...

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

Разностная краевая задача для уравнения смешанного типа

Краевая задача: Найти в области решение уравнение (1) условие: (2). Численное решение краевой задачи (1)-(2) является непростой задачей ввиду

Рахмонов Х. О. О первой краевой задаче для одного уравнения смешанного типа в пространстве. — Новосибирск, 1985. -22с.

Решение методом продолжения задач математической физики...

Для дифференциальных уравнений различают три типа задач: задача Коши; краевая задача; смешанная задача.

Используя это решение, мы можем решить задачу Коши. Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия

Организация вычислений решения краевой задачи для...

1. Краевая задача для ОДУ. Сведения о приближённых методах.

В проекционных методах [1,2] приближённое решение отыскивается в виде конечной суммы с неопределёнными коэффициентами

Аппроксимация первой краевой задачи разностной моделью для...

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

Модульный анализ сеточных методов решения...

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

Решение трехмерного уравнения Пуассона с использованием...

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

Обсуждение

Социальные комментарии Cackle

Похожие статьи

Априорная оценка для решения первой краевой задачи для...

Для решения краевой задачи (1) — (2), мы используем приближенной (численный) метод (метод разностных схем). Для доказательства устойчивости разностной модели, мы используем априорную оценку (3).

Численная реализация разностного метода решения одной...

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

Разностная краевая задача для уравнения смешанного типа

Краевая задача: Найти в области решение уравнение (1) условие: (2). Численное решение краевой задачи (1)-(2) является непростой задачей ввиду

Рахмонов Х. О. О первой краевой задаче для одного уравнения смешанного типа в пространстве. — Новосибирск, 1985. -22с.

Решение методом продолжения задач математической физики...

Для дифференциальных уравнений различают три типа задач: задача Коши; краевая задача; смешанная задача.

Используя это решение, мы можем решить задачу Коши. Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия

Организация вычислений решения краевой задачи для...

1. Краевая задача для ОДУ. Сведения о приближённых методах.

В проекционных методах [1,2] приближённое решение отыскивается в виде конечной суммы с неопределёнными коэффициентами

Аппроксимация первой краевой задачи разностной моделью для...

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

Модульный анализ сеточных методов решения...

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

Решение трехмерного уравнения Пуассона с использованием...

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

Априорная оценка для решения первой краевой задачи для...

Для решения краевой задачи (1) — (2), мы используем приближенной (численный) метод (метод разностных схем). Для доказательства устойчивости разностной модели, мы используем априорную оценку (3).

Численная реализация разностного метода решения одной...

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

Разностная краевая задача для уравнения смешанного типа

Краевая задача: Найти в области решение уравнение (1) условие: (2). Численное решение краевой задачи (1)-(2) является непростой задачей ввиду

Рахмонов Х. О. О первой краевой задаче для одного уравнения смешанного типа в пространстве. — Новосибирск, 1985. -22с.

Решение методом продолжения задач математической физики...

Для дифференциальных уравнений различают три типа задач: задача Коши; краевая задача; смешанная задача.

Используя это решение, мы можем решить задачу Коши. Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия

Организация вычислений решения краевой задачи для...

1. Краевая задача для ОДУ. Сведения о приближённых методах.

В проекционных методах [1,2] приближённое решение отыскивается в виде конечной суммы с неопределёнными коэффициентами

Аппроксимация первой краевой задачи разностной моделью для...

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

Модульный анализ сеточных методов решения...

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

Решение трехмерного уравнения Пуассона с использованием...

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

Задать вопрос