close

Вход

Забыли?

вход по аккаунту

?

Исследование линейной устойчивости адвективного течения методом сеток.

код для вставкиСкачать
Вычислительные технологии
Том 6, № 6, 2001
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ
АДВЕКТИВНОГО ТЕЧЕНИЯ МЕТОДОМ СЕТОК
Е. Л. Тарунин, К. Г. Шварц
Пермский государственный университет, Россия
e-mail: tarunin@psu.ru
The problem of the linear stability theory for advective flow in horizontal rotating
liquid layer is solved by reduction to the boundary problem for the system of linear onedimensional partial differential equations. The distinctions of grid method realization are
described for the given problem. The calculated neutral curves for normal perturbations
define the dependence on the critical Grashoff number by wave number for the different
Taylor number and Pr = 6.7 (water).
Исследование спектра малых нормальных возмущений конвективных течений и их линейной устойчивости сводится, как правило, к задаче на собственные значения для системы обыкновенных дифференциальных уравнений высокого порядка с переменными
коэффициентами и малыми параметрами при старших производных [1]. Наиболее распространенными численными методами решения таких задач являются метод Галеркина,
метод пошагового интегрирования с ортогонализацией и метод дифференциальной прогонки [2]. В данной работе задача линейной теории устойчивости решается путем сведения
ее к начально-краевой задаче для системы линейных одномерных уравнений в частных
производных, которая решается затем методом сеток. Аналогичная методика поиска критических чисел Рэлея для замкнутых областей использована в [3]. Возможности метода
демонстрируются на примере построения нейтральных кривых для адвективного течения,
возникающего во вращающемся горизонтальном слое жидкости [4]. Метод легко реализуется и, кроме того, дает полезную информацию о характере установления надкритических
движений. Тому, кто освоил метод сеток для решения нелинейных задач, для построения
нейтральных кривых достаточно убрать из уравнений нелинейность и зависимость от “излишних” координат.
Адвективное течение, возникающее во вращающемся плоском горизонтальном слое
жидкости, в котором ось вращения совпадает с вертикальной осью координат, для случая,
когда температура на границах линейно изменяется с продольной координатой, впервые
описано аналитически С. Н. Аристовым [5]. В дальнейшем было построено целое семейство
точных решений при различных граничных условиях. Особенностью этих течений является отсутствие вертикальной компоненты скорости. Вектор скорости в потоке ориентирован
перпендикулярно к силе плавучести (имеются обе его горизонтальные компоненты). Адвективное течение во вращающемся слое имеет экмановскую спиралевидную форму [6] и
наряду с изотермическим течением Экмана является основным пограничным течением в
c Е. Л. Тарунин, К. Г. Шварц, 2001.
°
108
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ АДВЕКТИВНОГО ТЕЧЕНИЯ 109
атмосфере и океане. Работ по исследованию устойчивости адвективного течения во вращающемся горизонтальном слое жидкости до последнего времени не было. Первая работа
в этой области [4], к сожалению, не лишена неточностей, которые здесь исправлены.
1. Постановка задачи
Рассмотрим плоский горизонтальный слой несжимаемой жидкости, ограниченный твердыми плоскостями Z = ±h и вращающийся с постоянной угловой скоростью Ω0 . Направление оси вращения совпадает с вертикальной осью координат. На обеих границах заданы
температура, линейно меняющаяся с координатой x, условие прилипания и замкнутости
потока
Zh
Z = ±h : T = Ax, A = const, v = 0,
vx dz = 0.
(1)
−h
Возникающее адвективное течение описывается аналитически [3]:
1
u0 (z) = √ Imf1 (z),
Ta
T0 = x + GrPr τ0 (z),
1
τ0 (z) = √ v0 (z),
Ta
1
v0 (z) = √ [z − Ref1 (z)],
Ta
GrPr
p0 = x(z + 1) + √
Ta
f1 (z) =
sinh λz
,
sinh λ
λ=
√
4
Zz
v0 (ζ)dζ,
−1
1+i
Ta √ ,
2
i=
√
−1,
(2)
где z = Z/h; u0 , v0 — горизонтальные компоненты вектора скорости, T0 — температура;
µ
¶2
ν
gβAh4
2Ω0 h2
p0 — давление; Pr = — число Прандтля; Gr =
— число Грасгофа; Ta =
χ
ν2
ν
— число Тейлора; Re, Im — действительная и мнимая части комплексного значения; ν —
кинематическая вязкость жидкости; χ — коэффициент теплопроводности; g — ускорение
свободного падения; β — коэффициент теплового расширения [4].
Линеаризованные уравнения для возмущений скорости V, температуры θ, давления P
в безразмерном виде имеют следующее представление во вращающейся системе координат:
√
∂V
+ Gr [(V∇)v0 + (v0 ∇)V] + Ta(iz × V) = −∇P + ∆V + θiz ,
∂t
divV = 0,
1
∂θ
+ Gr [(V∇)T0 + (v0 ∇)θ] =
∆θ,
∂t
Pr
z = ±1 :
V = 0,
θ = 0;
v0 = (u0 , v0 , 0),
V = (u, v, w),
(3)
(4)
iz = (0, 0, 1).
(5)
В силу большой сложности исследования устойчивости течения в трехмерной постановке будем рассматривать предельные случаи: пространственные винтовые периодические
по x возмущения в виде валов с осью, перпендикулярной к оси X, и пространственные
спиральные периодические по y возмущения в виде валов с осью, параллельной оси X.
Заметим, что при отсутствии вращения для Ta = 0 винтовые возмущения вырождаются
в плоские.
Е. Л. Тарунин, К. Г. Шварц
110
2. Винтовые возмущения
Уравнения винтовых возмущений выводятся из системы (3) – (5) в предположении, что
производные по y от всех функций равны нулю. Имеются все три компоненты вектора
возмущения скорости, которые так же, как и возмущения температуры, являются функциями времени t и двух пространственных переменных x, z. Учитывая дивергентность
возмущений скорости, введем функцию тока возмущений
u=−
∂ψ
,
∂z
w=
∂ψ
∂x
(6)
и вихря возмущения скорости
∂u ∂w
ϕ=
−
= −∆ψ = −
∂z
∂x
µ
∂ 2ψ ∂ 2ψ
+ 2
∂x2
∂z
¶
.
(7)
Представим неизвестные функции в виде
ϕ = [ϕ1 (t, z) + iϕ2 (t, z)] exp ikx x,
ψ = [ψ1 (t, z) + iψ2 (t, z)] exp ikx x,
v = [v1 (t, z) + iv2 (t, z)] exp ikx x,
θ = [θ1 (t, z) + iθ2 (t, z)] exp ikx x,
(8)
где kx — волновое число. В результате задача сведется к системе линейных уравнений в
частных производных по времени t и переменной z:
√ ∂v1
∂ 2 ϕ1
∂ϕ1
− kx Gr[u0 (z)ϕ2 + u000 ψ2 ] − Ta
=
− kx2 ϕ1 + kx θ2 ,
∂t
∂z
∂z 2
∂ 2 ψ1
− kx2 ψ1 + ϕ1 = 0,
∂z 2
√ ∂v2
∂ 2 ϕ2
∂ϕ2
+ kx Gr[u0 (z)ϕ1 + u000 ψ1 ] − Ta
=
− kx2 ϕ2 − kx θ1 ,
∂t
∂z
∂z 2
∂ 2 ψ2
− kx2 ψ2 + ϕ2 = 0,
2
∂z
√ ∂ψ1
∂v1
∂ 2 v1
− kx Gr[u0 (z)v2 + v00 ψ2 ] − Ta
=
− kx2 v1 ,
∂t
∂z
∂z 2
√ ∂ψ2
∂ 2 v2
∂v2
+ kx Gr[u0 (z)v1 + v00 ψ1 ] − Ta
=
− kx2 v2 ,
∂t
∂z
∂z 2
·
¸
∂ψ1
1 ∂ 2 θ1
∂θ1
0
2
− kx Gr[u0 (z)θ2 + GrPrτ0 ψ2 ] − Gr
=
− kx θ1 ,
∂t
∂z
Pr ∂z 2
·
¸
∂θ2
∂ψ2
1 ∂ 2 θ2
0
2
+ kx Gr[u0 (z)θ1 + GrPrτ0 ψ1 ] − Gr
=
− kx θ2
∂t
∂z
Pr ∂z 2
с граничными условиями
¯
¯
¯
¯
¯
¯
¯
∂ψ1,2 ¯¯
¯
¯
¯
ψ1,2 ¯
= v1,2 ¯
=
= θ1,2 ¯
= 0.
¯
¯
¯
¯
∂z ¯
z=±1
z=±1
z=±1
(9)
(10)
z=±1
Одномерная система уравнений (9) с граничными условиями (10) решается по численной методике, аналогичной схемам двухполевого метода [7], который используется для
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ АДВЕКТИВНОГО ТЕЧЕНИЯ 111
решения двумерных задач в переменных функции тока и вихря скорости. В нашем случае
функция тока и вихрь скорости зависят от одной пространственной переменной поперек слоя. Вторая переменная исчезла в предположении периодичности решения от этой
координаты. Уравнения для возмущения вихря и функции тока решаются с помощью
классической неявной схемы:
¶
µ
1 n+1
1 n+1
1
2
2
ϕ1,j+1 −
+ + kx ϕn+1
1,j + 2 ϕ1,j−1 =
2
2
hz
hz τ
hz
n
n
√ v1,j+1
− v1,j−1
1
n
n
] − Ta
− kx θ2,j
= − ϕn1,j − kx Gr[u0 (zj )ϕn2,j + u000 (zj )ψ2,j
τ
2hz
(0 < j < N ),
n
n
ϕn+1
ϕn+1
1,0 = ϕ1,0 ,
1,N = ϕ1,N ,
¶
µ
1 n+1
2
1 n+1
n+1
2
+ 2 ψ1,j−1
ψ1,j+1 −
+ kx ψ1,j
= −ϕn+1
1,j
2
2
hz
hz
hz
(0 < j < N ),
n+1
n+1
= 0, ψ1,N
= 0,
ψ1,0
¶
µ
1 n+1
1 n+1
1
2
ϕ2,j+1 −
+ + kx2 ϕn+1
2,j + 2 ϕ2,j−1 =
2
2
hz
hz τ
hz
n
n
√ v2,j+1
− v2,j−1
n+1
n+1
00
n
+ kx Gr[u0 (zj )ϕ1,j + u0 (zj )ψ1,j ] − Ta
+ kx θ1,j
2hz
1
= − ϕn2,j
τ
(0 < j < N ),
n
n
ϕn+1
ϕn+1
2,0 = ϕ2,0 ,
2,N = ϕ2,N ,
¶
µ
1 n+1
1 n+1
2
n+1
2
ψ
+
ψ
−
+
k
ψ
= −ϕn+1
x
2,j
2,j+1
2,j
h2z
h2z
h2z 2,j−1
n+1
ψ1,0
= 0,
(0 < j < N ),
n+1
ψ1,N
= 0,
(11)
где τ — шаг по времени; n = 0, 1, ... — номер шага по времени; hz = 2/N — шаг по
вертикальной координате; N — число интервалов разбиения. Вихри ϕ1 , ϕ2 на границах
аппроксимируются по формуле Вудса [4]:
ϕn+1
k,0 = −
3 n+1 1 n+1
ψ
− ϕk,1 ,
h2z k,1
2
ϕn+1
k,N = −
1
3 n+1
ψk,N −1 − ϕn+1
2
hz
2 k,N −1
(k = 1, 2).
(12)
Остальные уравнения системы (9), (10) решались с помощью схемы Кранка — Николсона
[8], погрешность аппроксимации порядка O(τ 2 + h2z ):
¶
µ
n
n
n
v1,j+1
− 2v1,j
+ v1,j−1
1 n+1
1 n+1
1 kx2
1 n
kx2 n
1
n+1
v1,j + 2 v1,j−1 = − v1,j −
v
−
+ +
+ v1,j −
2h2z 1,j+1
h2z τ
2
2hz
τ
2h2z
2
n
−kx Gr[u0 (zj )v2,j
+
n+1
v00 (zj )ψ2,j
]
−
√
n+1
n−1
ψ1,j+1
− ψ1,j−1
Ta
2hz
(0 < j < N ),
n+1
v1,0
= 0,
n+1
v1,N
= 0,
Е. Л. Тарунин, К. Г. Шварц
112
1 n+1
v
−
2h2z 2,j+1
µ
1 kx2
1
+
+
h2z τ
2
¶
n+1
v2,j
+
n
+kx Gr[u0 (zj )v1,j
n
n
n
v2,j+1
− 2v2,j
+ v2,j−1
1 n
kx2 n
1 n+1
v
=
−
+
v
−
v +
2h2z 2,j−1
τ 2,j
2h2z
2 2,j
+
n+1
v00 (zj )ψ1,j
]
−
√
n+1
n−1
ψ2,j+1
− ψ2,j−1
Ta
2hz
(0 < j < N ),
n+1
n+1
= 0,
v2,0
= 0, v2,N
¶
µ
n
n
n
θ1,j+1
− 2θ1,j
+ θ1,j−1
1
1
1 n
kx2
1
1
n+1
n+1
n+1
θ
+
θ
−
+ +
θ
= − θ1,j −
+
2Prh2z 1,j+1
Prh2z τ
2Pr 1,j
2Prh2z 1,j−1
τ
2Prh2z
+
n−1
n+1
ψ1,j+1
− ψ1,j−1
kx2 n
n+1
n
θ1,j − kx Gr[u0 (zj )θ2,j
+ PrGrτ00 (zj )ψ2,j
] − Gr
2Pr
2hz
(0 < j < N ),
n+1
n+1
θ1,0
= 0, θ1,N
= 0,
¶
µ
n
n
n
θ2,j+1
− 2θ2,j
+ θ2,j−1
1
1
1
1 n
kx2
1
n+1
n+1
n+1
θ
+
θ
−
+
θ
=
−
+
+
θ
−
2Prh2z 2,j+1
Prh2z τ
2Pr 2,j
2Prh2z 2,j−1
τ 2,j
2Prh2z
n+1
n−1
ψ2,j+1
− ψ2,j−1
kx2 n
n+1
n+1
0
θ + kx Gr[u0 (zj )θ1,j + PrGrτ0 (zj )ψ1,j ] − Gr
+
2Pr 2,j
2hz
(0 < j < N ),
n+1
θ2,0
= 0,
n+1
= 0.
θ2,N
(13)
В качестве начальных возмущений для неизвестных ψ1 , ψ2 , v1 , v2 , θ1 , θ2 берутся функции, пропорциональные sin2 πz, удовлетворяющие граничным условиям (10). Число узлов
сетки, определяемое шириной экмановского пограничного слоя и точностью расчетов, изменялось от 101 до 201. Расчеты проводились при фиксированном Pr = 6.7 (вода). Число
Тейлора менялось от 0 до 106 .
Для построения нейтральной кривой, описывающей зависимость критического числа
Грасгофа от волнового числа kx при фиксированном числе Тейлора, для каждого выбранного значения kx требуется найти такое число Грасгофа, при котором действительная
часть декремента возмущений λ = λ1 + iλ2 равна нулю. Иными словами, решается задача о поиске корня λ1 = 0 для неявной функции λ1 (kx , Gr, Ta). Мы строим эту функцию
по точкам (дискретно) с помощью многократного решения эволюционной задачи (9), (10)
методом сеток.
Для нахождения действительной части декремента возмущений λ1 прослеживалась
эволюция во времени максимумов по модулю неизвестных. В качестве аппроксимации
зависимости амплитуд от времени использовалась экспоненциальная формула C exp λ1 t.
Неизвестные C и λ1 определялись методом наименьших квадратов [9] по ходу вычислений
уравнений (11) – (13). Для определения неизвестных необходимо производить расчеты до
20 – 25 единиц безразмерного времени при шаге по времени порядка h2z , где hz — шаг по
вертикальной координате. Использованный алгоритм экономичен и позволяет рассчитывать коэффициенты аппроксимации по ходу основных расчетов. Нулевое значение декремента возмущений уточняется методом хорд. Характер поведения возмущений во времени
существенно зависит от всех параметров задачи; в области неустойчивости (λ1 > 0) все
возмущения нарастают, а в области устойчивости (λ1 < 0) затухают. Судить о характере
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ АДВЕКТИВНОГО ТЕЧЕНИЯ 113
устойчивости достаточно по какой-нибудь одной амплитуде. Опыт решения этой задачи
показал, что наиболее информативной является амплитуда возмущения температуры.
Расчеты показали, что для винтовых возмущений результаты работы [4] были ошибочны. Причиной ошибки явился неверный знак, стоящий перед одним из слагаемых в
разностной схеме. При заданном числе Прандтля наиболее опасной является колебательная неустойчивость, а не монотонная. Декремент возмущений λ в этом случае является
комплексным числом. Тестовые расчеты, проведенные для случая отсутствия вращения
(Ta = 0), показали, что критическое число Грасгофа Grk ≈ 126.7(1 − 31.6h2z ). Из этой
оценки видно, что относительная погрешность вычисления критического числа Грасгофа
составляет менее 0.32 % при hz ≤ 0.01.
Типичные зависимости амплитуды возмущений температуры θ1 от времени для двух
значений числа Грасгофа при фиксированных значениях числа Тейлора (Ta = 1000) и
волнового числа kx = 5.5 представлены на рис. 1. Первая кривая (Gr = 700) соответствует области колебательной неустойчивости с λ ≈ 1.644, вторая (Gr = 630) — области
устойчивости с λ1 ≈ −0.379. Граница устойчивости (λ1 ≈ 0) соответствует Gr = 643.5.
Рис. 1. Эволюция температуры θ1 во времени при Ta = 1000, k = 5.5: 1 — Gr = 700, 2 — Gr = 630.
Рис. 2. Нейтральные кривые для плоских возмущений при различных числах Тейлора (а); зависимость мнимой части декремента λ2 от Ta (б).
Мнимая часть декремента возмущения λ2 определялась после построения нейтральной кривой. Для этого выбиралось критическое число Грасгофа и соответствующее ему
волновое число и прослеживалась эволюция во времени какой-либо неизвестной (напри-
Е. Л. Тарунин, К. Г. Шварц
114
мер, u1 ) при фиксированном значении аргумента z (например, вблизи начала пограничного слоя). В качестве аппроксимации зависимости от времени использовалась формула
C sin λ2 t. Неизвестная C определялась по максимуму модуля прослеживаемой функции, а
период — методом наименьших квадратов.
На рис. 2, а представлены нейтральные кривые зависимости критического числа Грасгофа от волнового числа при различных значениях числа Тейлора. Как видно, вращение повышает устойчивость адвективного течения. С увеличением Ta растет критическое волновое число kx , соответствующее минимуму критического числа Грасгофа. На рис. 2, б
представлена зависимость мнимой части декремента возмущений от числа Тейлора при
Gr = Grk . Как видно, λ2 убывает при 0 < Ta ≤ 200. Для значений числа Тейлора больше
200 частота колебаний возмущений возрастает с увеличением Ta.
3. Спиральные возмущения
Уравнения спиральных возмущений выводятся из системы (3) – (5) в предположении, что
производные по x от всех функций равны нулю. Имеются все три компоненты вектора
возмущения скорости, которые так же, как и возмущения температуры, являются функциями времени t и двух пространственных переменных y, z. Учитывая дивергентность
возмущений скорости, введем функцию тока возмущений
v=−
∂Ψ
,
∂z
w=
∂Ψ
∂y
(14)
и вихря возмущения скорости
∂v ∂w
Φ=
−
= −∆Ψ = −
∂z
∂y
µ
∂2Ψ ∂2Ψ
+
∂y 2
∂z 2
¶
.
(15)
Представим все неизвестные функции в виде
Φ = [Φ1 (t, z) + iΦ2 (t, z)] exp iky y,
u = [u1 (t, z) + iu2 (t, z)] exp iky y,
Ψ = [Ψ1 (t, z) + iΨ2 (t, z)] exp iky y,
θ = [θ1 (t, z) + iθ2 (t, z)] exp iky y,
(16)
где ky — волновое число. В результате задача сведется к системе линейных уравнений в
частных производных по времени t и переменной z:
√ ∂u1
∂ 2 Φ1
∂Φ1
− ky Gr[v0 (z)Φ2 + v000 Ψ2 ] − Ta
=
− ky2 Φ1 + ky θ2 ,
∂t
∂z
∂z 2
∂ 2 Ψ1
− ky2 Ψ1 + Φ1 = 0,
∂z 2
√ ∂u2
∂ 2 Φ2
∂Φ2
+ ky Gr[v0 (z)Φ1 + v000 Ψ1 ] − Ta
=
− ky2 Φ2 − ky θ1 ,
∂t
∂z
∂z 2
∂ 2 Ψ2
− ky2 Ψ2 + Φ2 = 0,
2
∂z
√ ∂Ψ1
∂ 2 u1
∂u1
− ky Gr[v0 (z)u2 + u00 Ψ2 ] + Ta
=
− ky2 u1 ,
∂t
∂z
∂z 2
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ АДВЕКТИВНОГО ТЕЧЕНИЯ 115
√ ∂Ψ2
∂u2
∂ 2 u2
+ ky Gr[v0 (z)u1 + u00 Ψ1 ] + Ta
=
− ky2 u2 ,
∂t
∂z
∂z 2
·
¸
∂θ1
1 ∂ 2 θ1
0
2
− ky Gr[v0 (z)θ2 + GrPrτ0 Ψ2 ] − Gru1 =
− ky θ1 ,
∂t
Pr ∂z 2
·
¸
∂θ2
1 ∂ 2 θ2
0
2
+ ky Gr[v0 (z)θ1 + GrPrτ0 Ψ1 ] + Gru2 =
− ky θ2
∂t
Pr ∂z 2
с граничными условиями
¯
¯
¯
¯
¯
¯
¯
¯
∂Ψ
¯
¯
¯
1,2 ¯
Ψ1,2 ¯
= u1,2 ¯
=
= θ1,2 ¯
= 0.
¯
¯
¯
¯
∂z ¯
z=±1
z=±1
z=±1
(17)
(18)
z=±1
Система уравнений (17) с граничными условиями (18) решается с помощью неявных
разностных схем аналогично случаю винтовых возмущений. В качестве начальных возмущений для неизвестных Ψ1 , Ψ2 , u1 , u2 , θ1 , θ2 также берутся функции, пропорциональные
sin2 πz, удовлетворяющие граничным условиям (18). Число Тейлора менялось в пределах
от 0 до 40 000.
Расчеты подтвердили результаты работы [4]: при спиральных возмущениях возможна
колебательная неустойчивость.
На рис. 3, а представлены нейтральные кривые зависимости критического числа Грасгофа от волнового числа ky при различных значениях числа Тейлора. Заметно, что вращение
существенно повышает устойчивость адвективного течения. На рис. 3, б изображена зависимость “критической” мнимой части декремента λ2 от числа Тейлора (λ2 при этом соответствует минимуму нейтральной кривой Gr(ky )). Из представленной зависимости видно,
что с ростом числа Тейлора λ2 растет и, следовательно, увеличивается частота колебаний возмущений. С ростом Ta критическое волновое число увеличивается. Этот факт
необходимо учитывать при выборе шага по времени при расчетах методом сеток задачи
устойчивости.
Рис. 3. Нейтральные кривые для спиральных возмущений при различных числах Ta (а), зависимость мнимой части декремента от числа Тейлора (б).
Анализ нейтральных кривых для обоих случаев показал, что при малых значениях
числа Тейлора наиболее опасными являются спиральные возмущения, а при Ta > 60 —
винтовые пространственные возмущения (рис. 4).
116
Е. Л. Тарунин, К. Г. Шварц
Рис. 4. Зависимость критического числа Грасгофа от числа Тейлора: 1 — для случая спиральных
возмущений, 2 — для случая винтовых возмущений.
Выводы
1. Показана работоспособность представленного алгоритма для исследования линейной устойчивости адвективных течений при наличии вращения. Использовать метод особенно просто тем, кто освоил метод сеток.
2. Выявлена погрешность метода.
3. Исправлены результаты расчетов [4] и получены новые данные о границах устойчивости винтовых возмущений.
4. Построены нейтральные кривые для винтовых и спиральных возмущений при различных значениях числа Тейлора.
Список литературы
[1] Гершуни Г. З., Жуховицкий Е. М., Непомнящий А. А. Устойчивость конвективных течений. М.: Наука, 1989.
[2] Гольдштик М. А., Штерн В. Н. Гидродинамическая устойчивость и турбулентность. Новосибирск: Наука, 1977.
[3] Тарунин Е. Л. Определение границы конвективной устойчивости равновесия методом сеток // Алгоритмы и программы для ЭВМ. М., 1978. Деп. в ВИНИТИ 31 мая
1978, №1803–78(78).
[4] Аристов С. Н., Шварц К. Г. Об устойчивости адвективного течения во вращающемся горизонтальном слое жидкости // Изв. РАН. МЖГ. 1999. Т. 4. C. 3–11.
[5] Аристов С. Н., Зимин В. Д. Адвективные волны во вращающемся шаровом слое.
Свердловск, 1986 (Препр. / АН СССР. Уральский научный центр. Ин-т механики
сплошных сред).
[6] Аристов С. Н., Фрик П. Г. Динамика крупномасштабных течений в тонких слоях
жидкости. Свердловск, 1987 (Препр. / УрО АН СССР. ИМСС).
ИССЛЕДОВАНИЕ ЛИНЕЙНОЙ УСТОЙЧИВОСТИ АДВЕКТИВНОГО ТЕЧЕНИЯ 117
[7] Тарунин Е. Л. Вычислительный эксперимент в задачах свободной конвекции. Иркутск: Изд-во Иркутского ун-та, 1990.
[8] Самарский А. А. Теория разностных схем. М.: Наука, 1989.
[9] Демидович В. П., Марон И. А., Шувалова Э. З. Численные методы анализа.
Приближение функций, дифференциальные и интегральные уравнения. М.: Наука,
1967.
[10] Гершуни Г. З., Жуховицкий Е. М. Конвективная устойчивость несжимаемой
жидкости. М.: Наука, 1972.
Поступила в редакцию 28 декабря 2000 г.
в переработанном виде — 14 сентября 2001 г.
Документ
Категория
Без категории
Просмотров
3
Размер файла
193 Кб
Теги
методов, линейной, устойчивость, адвективного, сетон, исследование, течение
1/--страниц
Пожаловаться на содержимое документа