Лекции.ИНФО


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

Метод дискретных возмущений

 

Идея метода дискретных возмущений состоит в том, что в уравнение в каждую точку поочередно вводится дискретное возмущение e и прослеживается влияние этого возмущения на следующих временных слоях. Конечно-разностная схема будет устойчивой, если возмущение не возрастает, то есть выполняется условие (3.6).

Рассмотрим конечно-разностное уравнение (3.5). Введем возмущение e в точку (i,n):

 

.

 

Вычтя из этого уравнения с возмущением «невозмущенное» уравнение (3.5), получим уравнение для возмущений:

 

.

 

Отсюда:

 

.

 

Необходимо разрешить полученное неравенство относительно d:

 

-1≤1-2d≤1.

 

Рассмотрим отдельно правую и левую части этого неравенства:

 

а) 1-2d≤1 d≥0 — это условие выполняется всегда.

б) -1≤1-2d d≤1 — конечно-разностная схема (3.5) устойчива при выполнении этого условия.

Необходимо еще потребовать, чтобы знаки en+1 и en были одинаковыми во избежание осцилляций, то есть должно выполняться условие:

 

.

 

В данном случае это условие означает выполнение следующего неравенства:

 

1-2d≥0 Þ d≤1/2.

 

Мы получили более жесткое ограничение, чем в б), которое включает в себя условие . При фиксированных и a это условие накладывает ограничение на шаг по времени:

 

.

 

Это ограничение является жестким в смысле затрат времени для расчета на компьютере. Например, предположим, что расчет ведется с некоторым пространственным шагом 1. Максимально возможный шаг по времени . Если надо провести расчет с вдвое меньшим шагом 2=Dх1/2, то шаг по времени , т.е. шаг по времени придется уменьшить в 4 раза, а затраты компьютерного времени в общей сложности увеличатся в 8 раз. В двумерной задаче уменьшение вдвое шагов и Dy увеличивает затраты времени в 16 раз, а в трехмерной задаче - в 32 раза!

 

Введем теперь возмущение e в точку (i+1,n):

 

.

 

Вычтя из этого уравнения с возмущением «невозмущенное» уравнение (3.5), получим новое уравнение для возмущений:

 

.

 

Отсюда условие устойчивости имеет вид:

 

 

а) ,

(3.11)

 

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

 

Последнее неравенство можно переписать следующим образом:

 

 

Выражение слева представляет собой сеточное число Пекле: . Тогда последнее неравенство можно переписать так:

 

Pec<2.

 

Это условие можно воспринимать как ограничение на шаг по x:

 

,

 

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

 

б)

 

или

 

.

 

Поскольку , то это неравенство выполняется всегда.

 

Проведя аналогичный анализ для точки (i-1,n), получим:

 

 

Отсюда условие устойчивости:

 

 

Левая часть этого неравенства выполняется всегда, и никаких ограничений на шаги сетки не накладывает. Правая часть неравенства дает следующее условие:

 

 

Это условие более ограничительное, чем (3.11), поэтому условие (3.11) можно не принимать во внимание.

Таким образом, метод дискретных возмущений для уравнения (3.5) дает три следующих условия:

 

а)

б)

в)

 

Из условия б: uDх<2a, из условия a: . Следовательно, получаем еще одно условие:

 

г) .

 

 

Метод фон Неймана

 

Рассмотрим этот метод на примере того же конечно-разностного уравнения (3.5). Сеточная функция fin представляется в виде разложения Фурье:

 

fin = VnеIiq

 

где Vn – амплитуда отдельной компоненты с волновым числом на n-ном временном слое, q=kDx – фазовый угол, – мнимая единица. Аналогично

 

= Vn+1 е-I(i±1)q

 

Подставим эти выражения в уравнение (3.5):

.

 

Разделив на еIiq, получим:

 

.

 

Учтем, что еIq+ е-Iq=2cosq, еIq - е-Iq=2Isinq:

 

Vn+1=Vn [1-2d(1-cosq)-ICsinq]

 

Определим множитель перехода G следующим образом:

 

Vn+1=GVn

 

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

 

|G|≤1. (3.12)

 

В данном случае множитель перехода равен:

 

G=1-2d(1-cosq)-ICsinq

 

Он является комплексным выражением, поэтому условие устойчивости |G|≤1 приводится к неравенству:

 

[1-2d(1-cosq)]2-C2sin2q ≤ 1

 

Сделав соответствующие преобразования, получим:

 

C2(1+cosq) ≤ 4d[1-d(1-cosq)]

 

Рассмотрим два предельных случая.

 

а) cosq=-1: 0 ≤4d(1-2d).

 

Отсюда получаем первое условие устойчивости:

 

 

Такое же условие было получено методом дискретных возмущений.

 

б) cosq=1: .

 

Учитывая предыдущее неравенство, получим второе условие:

 

C ≤1.

 

Комбинируя эти два условия, можно получить условие устойчивости в следующем виде:

 

или Pec ≤ 2.

 

Все эти полученные условия справедливы только в случае линейного уравнения при u=const.

3.3.3 Метод практической устойчивости

 

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

Перепишем уравнение (3.5) в виде:

 

.

 

Условия устойчивости запишутся следующим образом:

 

;

;

– выполняется всегда;

– выполняется всегда.

 

Мы получили тремя разными методами одинаковые критерии устойчивости. Такой результат получается только для линейных дифференциальных уравнений с постоянными коэффициентами.

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

В 1952 г. В.С. Рябеньким была сформулирована теорема эквивалентности, которая устанавливает эквивалентность устойчивости и сходимости при выполнении следующих условий:

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

- конечно-разностное уравнение должно аппроксимировать дифференциальное уравнение в частных производных;

- устойчивость должна быть определена в форме фон Неймана.

При выполнении этих требований необходимое условие устойчивости становится и достаточным для сходимости.

 

Для нелинейных уравнений эта теорема не доказана, но на практике ее используют применительно и к нелинейным уравнениям.

Вообще говоря, исследование строгих определений аппроксимации, устойчивости, сходимости при Dх→0 и Dt→0 занятие зачастую бесплодное, т.к. реальные расчеты проводятся при конечных и Dt.









Читайте также:

Последнее изменение этой страницы: 2016-03-17; Просмотров: 193;


lektsia.info 2017 год. Все права принадлежат их авторам! Главная