ПУАССОНА УРАВНЕНИЕ


ПУАССОНА УРАВНЕНИЕ

; численные методы решения - методы, заменяющие исходную краевую задачу для уравнения Пуассона

(1)

системой из Nлинейных алгебраич. уравнений

LN(uN)=fN,(2)

решение к-рой позволяет построить нек-рую аппроксимацию pNuN для решения исходной задачи,.

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

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

Проекционные методы включают в себя ряд методов: вариационные, наименьших квадратов, Галеркина, проекционно-разностные, проекционно-сеточные, конечных элементов. Для всех них характерно сведение исходной краевой задачи к операторному уравнению

L(u)=f (3)

(оператор Lдействует, напр., из гильбертова пространства Нв Н).с последующим выбором конечномерных подпространств Н N и ; сама задача (3) в этих методах заменяется задачей нахождения такой, что для любого


Тогда при заданных базисах в Н N и FN система (2) является системой относительно коэффициентов разложения по базису Н N и за pNuN можно принять саму функцию ; погрешность метода естественно определить как . В наиболее важных случаях Нявляется нек-рым подпространством пространства Соболева , и если y1(x),y2(x), . . .,yN(x).- базис HN, то система (2) принимает вид

(4)

Погрешность метода при этом определяется расстоянием в H от решения исходной задачи до подпространства HN (см. [1], [5]-[9]). В современных вариантах проекционных методов подпространства Н N стремятся выбирать так, чтобы функции yi(x).имели локальные носители и в каждом уравнении (4) лишь конечное число коэффициентов было отлично от нуля. Методы такого типа и наз. проекционно-сеточными методами (проекционно-разностными, вариационно-разностными, конечных элементов) (см. [1], [4], [7] - [9], [11]). Наибольшим достоинством этих методов является их применимость при достаточно сложной геометрии области W, в к-рой рассматривается краевая задача. К проекционным методам примыкает и относительно редко применяемый коллокаций метод.

Разностные (конечноразностные) методы используют аппроксимацию исходной области W нек-рой сеточной областью WN, содержащей Nузлов сетки, и обычно приводят к системе (2) на основе аппроксимации П. у. и соответствующих граничных условий их разностными (сеточными) аналогами, использующими лишь значения функции в выбранных узлах (см. [1]). Погрешность метода обычно получается сравнением вектора uN и вектора, получаемого сужением искомого решения на множество рассматриваемых узлов. Корректность и аппроксимация могут изучаться при различных выборах норм, в частности возможно использование принципа максимума; сходимость получается как следствие корректности и аппроксимации (см. [1]-[4]).

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

Методы решения систем сеточных уравнений (2) наиболее интенсивно изучались для простейших разностных аналогов на параллелепипедной сетке (см. [1], [9], [11]-[19]). В случае двух переменных и области W на плоскости, являющейся прямоугольником, часто применяются для ряда граничных условий прямые методы, позволяющие найти решение (2) при затрате O(Nln N).арифметич. действий. Это - метод разделения переменных, использующий дискретное преобразование Фурье, и метод редукции (см. [12], [13]); известны и методы с оценкой O(N).действий (см. [14]). При и наличии разделения переменных (W в атом случае - параллелепипед) решение (2) с точностью е>0 можно найти при затрате O(Nln N|lne|) действий с помощью итерационного метода переменных направлений (см. [1], [12]); итерационные методы с факторизуемыми операторами (метод последовательной сверхрелаксации с симметризацией, неполной матричной факторизации, попеременно-треугольный) позволяют найти решение (2) с точностью e при затрате 0(N1+1/2d |ln e|) действий (см. [12], [15]) для довольно общих ситуаций.

В случае однесвязных и многосвязных областей W на плоскости, составленных из конечного числа прямоугольников, решение системы (2) с точностью е может быть найдено при затрате 0(Nln N+N3/4 |lne|In N) действий на основе разрезов W на прямоугольники (см. [15], [16]). Похожая асимптотика для дискретных аналогов задачи Дирихле в случае некоторых областей получена с помощью метода емкостей (см. [17]) и фиктивных неизвестных (см. [18]). Для ряда систем (2), являющихся проекционно-разностными аналогами исходных задач, затраты типа O(Nln N |lne|), a иногда и О(Nln N) (при ) достигаются с помощью итерационных методов, использующих эквивалентность операторов по спектру (см. [11], [17], [19]). Использование последовательностей сеток в ряде случаев позволяет получить итерационные методы, дающие решение (2) с точностью , при асимптотически минимальных затратах вычислительной работы (число действий есть 0(N)).(см., напр., [1], [9], [11], [19]).

Лит.:[l] Годунов С. К., Рябенький В. С., Разностные схемы, 2 изд., М., 1977; [2] Ладыженская О. А.., Краевые задачи математической физики, М., 1973; [3] Марчук Г. И., Шайдуров В. В., Повышение точности решений разностных схем, М., 1979; [4] СамарскийА. А.,Андреев В. Б., Разностные методы для эллиптических уравнений, М., 1976; [5] Михлин С. Г., Численная реализация вариационных методов, М., 1966; [6] Красносельский М. А. [и др.], Приближенное решение операторных уравнений, М., 1969; [7] Обэн Ж. П., Приближенное решение эллиптических краевых задач, пер. с англ., М., 1977; [8] Стренг Г., Фикс Дж., Теория метода конечных элементов, пер. с англ., М., 1977; [9] Оганесян Л. А., Ривкинд В. Я., Руховец Л. А., Вариационно-разностные методы решения эллиптических уравнений, ч. 1-2, Вильнюс, 1973-74 (Дифференциальные уравнения и их применение, в. 5, 8); [10] Самарский А. А., Фрязинов И. В., "Успехи матем. наук", 1976, т. 31, в. 6, с. 167-97; [11] ДьяконовЕ. <Г., в кн.; Вариационно-разностные методы в математической физике, Новосиб., 1978, с. 149-64; [12] Самарский А. А., Николаев Е. С., Методы решения сеточных уравнений, М., 1978; [13] Ваrkеr R. J., в кн.: Mathematical models and numerical methods, Warsz., 1978, p. 255-68 (Banach center publications, v. 3); [14] Bank R. E., Rоse D. J., "SIAM J. Numer. Anal.", 1975, v. 12, № 4, p. 529-40; [15] Кузнецов Ю. А., в кн.: Вариационно-разностные методы в математической физике, Новосиб., 1978, с. 178-212; [16] Дьяконов Е. Г., в кн.: Численные методы в математической физике, Новосиб., 1979, с. 45-68; [17] Астраханцев Г. П., в кн.: Разностные и вариационно-разностные методы, Новосиб., 1977, с. 63- 72; [18] Капорин И. Е., Николаев Е. С., "Дифференциальные уравнения", 1980, т. 16, № 7, с. 1211-25; [19] Корнеев В. Г., Схемы метода конечных элементов высоких порядков точности, Л., 1977. Е. Г. Дьяконов.


Математическая энциклопедия. — М.: Советская энциклопедия. . 1977—1985.

Смотреть что такое "ПУАССОНА УРАВНЕНИЕ" в других словарях:

  • ПУАССОНА УРАВНЕНИЕ — дифференциальное уравнение д2u/дx2+д2u/дy2+д2u/дz2= 4pr(x, y, z) одно из осн. ур ний теории потенциала. Так, П. у. определяет потенциал и в точке с координатами х, у, z в электростатич. поле, создаваемом электрич. зарядами с объёмной плотностью… …   Физическая энциклопедия

  • ПУАССОНА УРАВНЕНИЕ — уравнение с частными производными вида ?u= f, где ? Лапласа оператор. Изучено С. Пуассоном …   Большой Энциклопедический словарь

  • Пуассона уравнение — уравнение с частными производными вида Δu = f, где Δ  Лапласа оператор. Впервые изучено С. Пуассоном. * * * ПУАССОНА УРАВНЕНИЕ ПУАССОНА УРАВНЕНИЕ, уравнение с частными производными вида Du= f, где D Лапласа оператор (см. ЛАПЛАСА ОПЕРАТОР).… …   Энциклопедический словарь

  • Пуассона уравнение —         уравнение с частными производными вида Δu = f, где Δ оператор Лапласа:                   При n = 3 этому уравнению удовлетворяет Потенциал u (х, у, z) объёмных масс, распределённых с плотностью f (x, у, z)/4π (в областях, где f = 0… …   Большая советская энциклопедия

  • ПУАССОНА УРАВНЕНИЕ — уравнение с частными производными вида Дельта u = f, где Дельта Лапласа оператор …   Большой энциклопедический политехнический словарь

  • Пуассона уравнение — Уравнение Пуассона эллиптическое дифференциальное уравнение в частных производных, которое, среди прочего, описывает электростатическое поле, стационарное поле температуры, поле давления, поле потенциала скорости в гидродинамике. Оно названо в… …   Википедия

  • ПУАССОНА УРАВНЕНИЕ — дифференциальное уравнение с частными производными, к рому удовлетворяет объемный потенциал внутри областей, занятых создающими этот потенциал массами. Для ньютонова потенциала в пространстве , и логарифмического потенциала в П. у. имеет вид где… …   Математическая энциклопедия

  • ПУАССОНА УРАВНЕНИЕ — ур ние с частными производными вида дельта u = f, где Д оператор Лапласа. Впервые изучено С. Пуассоном …   Естествознание. Энциклопедический словарь

  • Уравнение Власова — Уравнение Власова  система уравнений, описывающих динамику плазмы заряженных частиц с учётом дальнодействующих кулоновских сил посредством самосогласованного поля. Впервые предложена А. А. Власовым в статье[1] и позднее излагается… …   Википедия

  • Уравнение Лапласа — Уравнение Лапласа  дифференциальное уравнение в частных производных. В трёхмерном пространстве уравнение Лапласа записывается так: и является частным случаем уравнения Гельмгольца. Уравнение рассматривают также в двумерном и одномерном… …   Википедия

Книги