Главная Численные методы при исследовании физических задач



ГЛАВА IX

УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ

В главе IX рассмотрены методы численного решения задач для уравнений в частных производных. В § 1 обсуждены некоторые постановки задач и дан обзор методов, которыми решаются подобные задачи. Остальные параграфы содержат изложение основ наиболее широко применяемого и хорошо изученного метода -разностного. В § 2 рассмотрены способы построения разностных схем и введено понятие аппроксимации. В § 3 даны методы исследования устойчивости разностных схем. В § 4 доказаны основные теоремы о сходимости разностного решения к точному.

§ 1. Введение

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

Математические модели сплошной среды приводят к уравнениям в частных производных, которым удовлетворяют упомянутые средние величины. Например, изменение температуры в неподвижном теле описывается уравнением теплопроводности

с {и, r,t) = div [k («, /-, t) grad u] + q {u, r, t), (1)

где « - температура, с -теплоемкость, - коэффициент теплопроводности и -плотность источников тепла.

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

Независимыми переменными в физических задачах обычно являются время / и координаты г; бывают и другие переменные, например, скорости частиц v в задачах переноса. Решение требуется найти в некоторой области изменения независимых пере-



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

Если одной из переменных является t, то чаще всего рассматривают области вида

G(t, г, ...)=gir,...)x[to,n (2)

т. е. решение ищут в некоторой пространственной области g(r, ...) на отрезке времени ttT. В этом случае дополнительные условия, заданные при t = to, называют начальными, а дополнительные условия, заданные на границе Г (г) области g (г), ~ граничными или краевыми.

Задачу, у которой имеются только начальные условия, называют задачей Коши. Например, для уравнения теплопроводности (1) в неограниченном пространстве можно поставить задачу с начальными условиями

«(г, g = [x(r). (3)

Если [X (г) - кусочно-непрерывная ограниченная функция, то решение задачи (1), (3) единственно в классе ограниченных функций (при некоторых ограничениях на коэффициенты уравнения; см. [40]).

Задачу с начальными и граничными условиями называют смешанной краевой задачей или нестационарной краевой задачей. Для уравнения (1) дополнительные условия такой задачи могут иметь, например, вид

и {г, g = р(г), rg{г), и {г, t)r = [Xi (г, t), totT. (4)

Для этого уравнения допустимы и другие граничные условия, например, содержащие нормальную производную решения.

Зстречаются задачи, в которых область G (t, г) имеет другой вид. Примером является задача с условиями на характеристиках (см. [40]), возникающая при изучении процессов сушки, сорбции газов и многих других.

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

Мы ограничимся рассмотрением корректно поставленных задач, когда для некоторого класса начальных и граничных данных решение (в заданном классе функций) существует, единственно и непрерывно зависит от этих данных. Будем также предполагать, что решение непрерывно зависит от всех коэффициентов уравнения.



Для уравнений в частных производных существуют физически интересные задачи, являющиеся некорректно поставленными; обратные задачи тепло1;ро-водности, задачи на развитие неустойчивостей и другие. Так, рассмотренный в главе I пример Адамара связан с возникновением релей-тейлоровской неустойчивости, когда слой тяж-;лой жидкости налит поверх слоя легкой жидкости. Но здесь мы такие задачи не будем рассматривать (см. [39] и приведенную там библиографию).

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

Напомним классификацию таких уравнений. Они имеют следующий вид (для простоты мы ограничиваемся случаем двух переменных):

Аи + 2Виу + Cuyy -f D«., + Euy + F = Q. (5)

Коэффициенты уравнения (5), вообще говоря, зависят от и, х, у. Если коэффициенты не зависят от переменных, то это линейное уравнение с постоянными коэффициентами. Если F линейно зависит от и, а остальные коэффициенты от и не зависят, то это линейное уравнение с переменными коэффициентами. Если коэффициенты зависят от ы, то уравнение (5) называется квазилинейным.

Если А = B = C = Q, но ОфО и ЕфО,то уравнение (5) имеет первый порядок и называется уравнением переноса.

Уравнения второго порядка классифицируются по знаку дискриминанта В - АС: у гиперболических уравнений дискриминант положителен, у параболических - равен нулю, у эллиптических - отрицателен.

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

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

2. Точные методы решения. В курсах уравнений математической физики изложен ряд методов, позволяющих для некоторьщ классов задач найти точное решение (см. [40]). К таким методам



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 [96] 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170


0.0143