Главная Обратная связь

Дисциплины:






Принцип замороженных коэффициентов. Критерий Гельфанда – Бабенко



Обозначим du(ut)= v(ht)u(ht)сеточную функцию вычислительной погрешности. В силу линейности операторов Lh и lh вычислительная погрешность является решением разностной задачи

+ Lht(du(ht)) = e(ht)в Dht ,

lht(du(ht)) = e(ht)на дDh´Dt ,

du(ht) = e(ht), =e(ht) при j=Dh´{t=0}.

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

В задаче, определяющей эволюцию вычислительной погрешности du(ht), возмущения присутствуют в правой части разностных уравнений и в начальных и граничных условиях. Задание произвольного возмущения e(h) на начальном временном слое приведет к появлению возмущений на следующих слоях. С другой стороны, при исследовании задачи можно проследить за развитием возмущения, начиная с любого временного слоя, приняв его за начальный. Поэтому, достаточно рассмотреть задачу для однородного уравнения с однородными граничными условиями.

+ Lht(du(ht)) = 0 в Dht , lht(du(ht)) = 0 на дDh´Dt , du(ht) = e(h), = e(h) при j=Dh´{t=0} (3.5)

Основной метод исследования устойчивости разностной задачи называется принципом замороженных коэффициентов. Рассмотрим одномерную эволюционную задачу (2.17) с коэффициентами, зависящими непрерывно от x и t. Тогда каждый узел (m,k) сетки Dh принадлежит области существования решения дифференциальной задачи с некоторой окрестностью W(m,k), в которой коэффициенты уравнения можно считать приблизительно постоянными. При неограниченном измельчении сетки в W(m,k) окажется множество узлов. Поскольку для любого внутреннего узла расстояние до границы, измеряемое шагами сетки, велико, распространение возмущений в W(m,k) будет описываться задачей Коши с постоянными коэффициентами, причем для каждого узла (m,k) свои условия устойчивости. Из всей совокупности условий необходимо оставить наиболее ограничительные. Тогда ни в одном узле возмущение не будет возрастать, разностная схема будет устойчивой.

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

Принцип замороженных коэффициентов дает асимптотически необходимое условие устойчивости, т.е. условие, которое выполняется при предельном переходе t ® ¥ и h ® ¥. Реализация принципа замороженных коэффициентов приводит к исследованию в бесконечной и полубесконечной области задач с постоянными коэффициентами, повторяющими все значения коэффициентов реальной задачи в «замороженном» виде.



Принцип замороженных коэффициентов реализован в критерии Гельфанда – Бабенко, состоящем в следующем.

1. В окрестности каждого внутреннего узла исследуется задача Коши с постоянными коэффициентами в бесконечной области. За коэффициенты задачи принимаются их значения в узле. Возмущение ограничено при m ® ¥ и m ® -¥ на каждом временном слое.

+ Lht(du(ht)) = 0, -¥ < m < +¥, kt £ T, du(ht) = e(h), -¥ < m < +¥, k = 0, |du(ht)| < ¥ при m ® ±¥, kt £ T. (3.6)

2. В окрестности каждого граничного узла (0,k) исследуется краевая задача с постоянными коэффициентами в полубесконечной области. За коэффициенты задачи принимаются их значения в узле. Возмущение вызвано краевым условием и при m ® ¥ убывает на каждом временном слое.

+ Lht(du(ht)) = 0, 0 < m < +¥, kt £ T, lht(du(ht)) = 0, … m = 0, kt £ T, du(ht) = e(h), 0 £ m < +¥, k = 0, |du(ht)| ®0при m ® ¥. (3.7)

3. В окрестности каждого граничного узла (М,n) исследуется краевая задача с постоянными коэффициентами в полубесконечной области. За коэффициенты задачи принимаются их значения в узле. Возмущение вызвано краевым условием и при m ® -¥ убывает на каждом временном слое.

+ Lht(du(ht)) = 0, -¥ < m < M, kt £ T, lht(du(ht)) = 0, m = M, kt £ T, du(ht) = e(h), -¥ < m < M, k = 0, |du(ht)| ®0при m ® -¥. (3.8)

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

| | £ d.

Определимся с выбором начального возмущения. Задача Коши в бесконечной области для дифференциального уравнения с постоянными коэффициентами имеет точное решение

u(x,t) = ,

где i = – мнимая единица. В задаче переноса c = -iaw, в задаче теплопроводности c = - aw2, в волновой задаче c = . Решение представляет суперпозицию волн (гармоник) вида , функция C(w) определяется из разложения начальной функции в ряд Фурье. В соответствии с аппроксимацией разностную гармонику можно представить в виде =Cw lkeiwhm, а все решение – как суперпозицию гармоник. Множитель l может быть комплексным числом и является показателем роста или затухания соответствующей гармоники. Следовательно, в задаче 1 критерия Гельфанда – Бабенко за начальное возмущение можно взять em = e eiwhm, где ||e(h)|| = e,а решение задачи искать в виде

= e lkeiwhm

Поскольку |eiwhm| =1, все условия задачи (3.6) будут выполнены.

Подставим в однородное разностное уравнение и получим уравнение относительно l, называемое характеристическим. Чтобы вычислительная погрешность в расчетной области kt £ T не искажала результат, необходимо, чтобы при любом w во всех узлах выполнялось неравенство

| | = e×|l|k£ d,

откуда

ln(|l|ln(d/|e|)/k £ ln(d/|e|)t/T.

Положим

c = ln(B/|e|)/T.

Тогда

|l| £ ect= 1 + ct + ½ ecq(ct) 2, 0 < q < t.

Усилим неравенство, отбросив в формуле остаточный член. Получим необходимое условие устойчивости Неймана

|l| £ 1 + ct. (3.9)

Константу c можно заранее оценить. Пусть, например, расчет необходимо провести для 0£ t £10с погрешностью 5знаков после запятой. Пусть погрешность, порождаемая вычислениями с округлением, не превосходит 10–12в каждом узле. Приняв T = 10, e = 10–12, d = 10–5, получим

c = ln(10–5/10–12)/10 = 1,612.

Если будет выполнено условие Неймана с константой c £ 1,612, начальная погрешность возрастет до величины, не превосходящей 10–5.

Если характеристическое уравнение не содержит t,условие Неймана

|l| £ 1 (3.10)

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

При исследовании разностной задачи Коши с постоянными коэффициентами решается однородное разностное уравнение, граничные условия не используются. Разностная схема задачи Коши с постоянными коэффициентами при h ® 0и t ®0 устойчива,если решение характеристического уравнения при любом w удовлетворяет условию (3.9) или (3.10). Таким образом, для задачи Коши с постоянными коэффициентами условие Неймана является необходимым и достаточным.

Выпишем вспомогательные формулы:

= e lkeiwhm, = e lk+1eiwhm, = e lk-1eiwhm, = e lkeiwh(m+1), = e lkeiwh(m-1), eiwh = cos(wh) + i sin(wh), cos(wh) = (eiwh + e-iwh)/2, sin(wh) = (eiwh + e-iwh)/(2i), 1 - cos(wh) = 2 sin(½ wh), sin(wh) = 2 sin(½ wh) cos(½ wh). (3.11)

 

 





sdamzavas.net - 2020 год. Все права принадлежат их авторам! В случае нарушение авторского права, обращайтесь по форме обратной связи...