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



ГЛАВА XII

ЭЛЛИПТИЧЕСКИЕ УРАВНЕНИЯ

Глава XII посвящена методам решения краевых задач для эллиптических уравнений. В § 1 решение таких задач сводится к решению эволюционных задач для параболических уравнений до выхода на стационарный режим; последнее выполняется при помощи многомерных разностных схем, изложенных в гл. XI, § 2. Обсужден выбор оптимального шага по времени (или набора переменных шагов) в таких расчетах.

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

§ 1. Счет на установление

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

Примером полной математической постановки является задача с краевыми условиями первого рода, называемая задачей Дирихле; требуется найти непрерывное решение задачи

Ан(г) = -/(г), reG, иг (г) = 11 (г), (1)

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

div[fe(r)grad«(r)] = -/(г), rG, Ur{r) = ii{r), k{r)>0.



Wr (r, 0 = 0, W {/, 0) = йУо (r) = Уо {r) - и (r).

Поскольку начальные данные в (3) были выбраны произвольно, то без ограничения обшности можно считать, что начальные данные задачи (4) также выбраны произвольно.

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

w{r,t)==Zcge~hwg{r). (5)

Здесь Wg (г) и Хд - собственные функции и собственные значения многомерной задачи Штурма -Лиувилля:

dw[kir)gTauWg] + XgWg{r) = 0, г G, Wg(r)r = 0, (6)

Cgw {г, 0) Wg (г) dr

являются коэффициентами Фурье начальных данных (4) по системе функций Wg(r). Собственные значения задачи (6) положительны и образуют неубывающую последовательность

0<?ii=s;?i2<?i3< (7)

а собственные функции Wg (г) образуют полную ортонормированную систему в G (г).

Задачи с другими краевыми условиями мы не будем рассматривать.

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

= div [/г (г) grad у (л 0]+ /(/-), rG, 0t< + oo, vrir, t) = \iir), v(r, 0) = Uo(r).

Исследуем, насколько решение v (г, t) эволюционной задачи отличается от решения и (г) стационарной задачи. Вычитая (2) из (3) и учитывая, что du{r)/dt = 0, найдем, что разность w(r, t) = v{r, t) - u(r) удовлетворяет однородному параболическому уравнению с однородными краевыми условиями:

-w{r, t) = div[k{r)gradw], rG, 0<< + сх),



„. Л

/со \ 1/2

11 (Г, о 1U. = 2 cle- V е-V 2 4) = V н (г)

Оно означает, что разность ш(г, t) = v {г, t) - и{г) при /->-оо экспоненциально стремится к нулю по норме • Ц, так что решение v (г, t) эволюционной задачи (3) среднеквадратично сходится к решению и (г) стационарной задачи (2) при t-oo.

Замечание 1. Пусть граничные и начальные условия таковы, что решения задач (2) и (3) имеют в G (г) непрерывные производные, ограниченные равномерно по t. Тогда сходимость V (г, t) к и (г) будет равномерной.

Таким образом, вместо задачи (2) для эллиптического уравнения можно взять эволюционную задачу (3) для параболического уравнения с тем же пространственным оператором, произ; вольно выбрать начальные данные и вычислить решение v (г, t) при достаточно большом t. Стационарный (не зависящий от времени) предел и (г), к которому стремится v{r, t) при i->-oo, и будет решением стационарной задачи (2).

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

Установление стационарного решения происходит довольно быстро благодаря экспоненциальному характеру затухания начальных данных. Из (8) видно, что если нужна точность ~б, то надо вести вычисления до момента

- • 1п, (9)

где Ху есть наименьшее собственное значение соответствующей задачи Штурма - Лиувилля (6).

Замечание 2. На стационарное решение выходят не только решения параболических задач. То же происходит при других диссипативных процессах со стационарными граничными условиями, например при колебаниях с вязким трением, описываемых уравнением

o« + py, = div(gradf) + f, fr = x(/-), р>0. (10)

Можно формулировать эволюционную задачу для этого уравнения; однако это менее удобно.

Замечание 3. Можно составить разностную схему, непосредственно аппроксимирующую исходную задачу (2). Но в § 2

Из (5) и (7) нетрудно получить неравенство



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.0331