Проблема восприимчивости погранслойных течений и бифуркации периодических режимов
Автор: Дармаев Тумэн Гомбоцыренович
Журнал: Вестник Бурятского государственного университета. Философия @vestnik-bsu
Рубрика: Функциональные уравнения и их приложения
Статья в выпуске: 9, 2011 года.
Бесплатный доступ
В работе метод инвариантной конечномерной проекции уравнений Навье-Стокса применяется для исследования плоскопараллельного течения вязкой несжимаемой жидкости над плоской полубесконечной пластиной. При этом нелинейная начально-краевая задача для периодических возмущений основного течения сводится к конечномерной рекуррентной системе линейных краевых задач для обыкновенных дифференциальных уравнений. Предложенный универсальный алгоритм позволяет численно находить амплитудные поверхности устойчивых и неустойчивых режимов и точки тангенциальных бифуркаций периодических режимов для произвольных частот и чисел Рейнольдса. Проведено сравнение численных результатов с проблемой восприимчивости и экспериментами по ламинарно-турбулентному переходу.
Ламинарно-турбулентный переход, метод инвариантной конечномерной проекции
Короткий адрес: https://sciup.org/148180507
IDR: 148180507
Текст научной статьи Проблема восприимчивости погранслойных течений и бифуркации периодических режимов
Основная задача исследования восприимчивости заключается в поиске, выборе, экспериментальном и теоретическом моделировании таких механизмов возбуждения колебаний, которые играют основную роль в формировании возмущенного течения и последующей его турбулизации [1]. К основным внешним факторам, которые вызывают возмущения, относятся шероховатость или геометрия поверхности, завихренность и турбулентность свободного потока (рис.1).
основной поток

Рис.1. Механизмы восприимчивости
В зависимости от геометрии обтекаемого тела, размера шероховатостей, степени турбулентности или завихренности основного течения появляются возмущения и в дальнейшем может произойти их рост и переход в турбулентность. Кроме этого в экспериментах для генерации вынужденных возмущений применяются акустические методы, вибрирующая лента, вдув или отсос в поток, которые в зависимости от частоты или амплитуды могут перейти в турбулентное состояние.
В настоящее время развиваются следующие теоретические направления исследования механизма восприимчивости – параболизованные уравнения [2], прямое численное моделирование [3], возмущения в профиле параллельного пограничного слоя [4-5], взаимодействие возмущений [6-7]. В них рассчитываются коэффициенты восприимчивости возмущений от начальных данных, задающих внешние факторы, амплитуды и частоты возмущений.
В данной работе метод инвариантной конечномерной проекции уравнений Навье-Стокса, разработанный Б.Ю.Скобелевым [8], применяется для исследования периодических трехмерных возмущений в плоскопараллельном течении вязкой несжимаемой жидкости над плоской полубесконечной пластиной. При этом нелинейная начально-краевая задача для возмущений основного течения сводится к конечномерной рекуррентной системе линейных краевых задач для обыкновенных дифференциальных уравнений [9]. Предложенный универсальный алгоритм позволяет численно находить амплитудные поверхности устойчивых и неустойчивых режимов и точки тангенциальных бифуркаций периодических режимов для произвольных частот и чисел Рейнольдса.
Метод инвариантной проекции
Суть метода инвариантной конечномерной проекции состоит в том, что решение уравнений Навье-Стокса для возмущений:
u + ( u ∇ ) u + Lu = 0,
∂ t
Lu = (u ∇) + (u∇)u - ∇q - ∆u,
0 0 Re div u = 0
ищется в виде:
u = v2n + u⊥, где v2n – линейная функция m собственных функций оператора L, а ортогональное дополнение u⊥ ищется в виде функции от v2n u⊥ =g(v2n).
Условие того, что решение принадлежит инвариантному многообразию уравнений Навье-Стокса однозначно определяет вид функции g ( v 2 n ) . Существенным достоинством метода инвариантной проекции является то, что он гарантирует правильное описание асимптотического поведения решений (т.е. при t → ∞ ) и учитывает дискретный и непрерывные спектры возмущений.
В работе рассматривались моногармонические решения, периодические по продольной координате с волновым числом α и по поперечной координате с волновым числом β .
В соответствии с теорией инвариантной конечномерной проекции решение уравнения для функции тока представлялось в следующем виде:
∞s ф (x, y, z, t) = 2 A(t) Re [ f (y) ё(ax+et)) ] + £ As(t) £ gsk (y) e(az+e-e(t», s=1 k=-s где A(t) и θ(t) – амплитуда и фаза возмущения, δ – символ Кронекера, f(у) – собственная функция задачи Орра-Зоммерфельда для течения Блазиуса, а gsk удовлетворяют следующей рекуррентной системе уравнений:
ik a _ ( №. - U ) A k + DU ] g sk - ^( A k )2 g sk = -5 k !( b s - i + ic s — 1 ) A 1 f - 5 k , — j *
R (1)
-
* ( b s - 1 — ic s - 1) A 1 f — ik £ c q A k g pk - i a £ £ [ lgqiD A j g pj - jDg ql A l g pj ] , q + p = s q + p = s l + j = k
где D = d , ∆ = D 2 - ( k α ɶ ) 2 , α ɶ 2 = α 2 + β 2 , U(y) – профиль невозмущенного ламинарного потока, f dy k
– сопряженная к собственной функции задачи Орра-Зоммерфельда, σ = σ i + i σ r – собственное число задачи Орра-Зоммерфельда.
Из условий прилипания на пластине и затухания возмущений на бесконечности ставятся следующие граничные условия:
g sk (0) = Dg sk (0) = 0, g sk ( ∞ ) = Dg sk ( ∞ ).
Коэффициенты bs - 1, cs - 1 определяются из условий ортогональности и нормировки:
∞∞
∫ f ( y ) ∆ 1 qs 1( y ) dy = 0, ∆ 1 = D 2 - α ɶ 2 , bs - 1 + ics - 1 = ∫ f * Fs 1 dy , 00
где f*(y) – решение уравнения, сопряженного к уравнению Орра-Зоммерфельда.
Результаты расчетов и их анализ
В работе проводился численный расчет приближенных уравнений для амплитуд периодических режимов
∞
γ + ∑ b 2 n Α 2 n = 0. n = 1
при N =2,3,4,5, где γ = ασ i – линейный коэффициент нарастания.
При численном интегрировании системы (1) граничные условия на бесконечности заменялись граничными условиями на интервале (0,L), L>>1. Затем преобразованием y = L y , 0 ≤ y ≤ L , 0 ≤ y ≤ 1 задача сводится к интегрированию на интервале (0,1). Вычисление присутствующих в правой части собственных функций задачи Орра-Зоммерфельда, а также решение системы (1) проводилось методом ортогональной прогонки [10] на неравномерной разностной сетке.
В результате вычислений выявилась следующая картина (рис.2-3). При малых значениях волнового числа α, соответствующих нижней ветви линейной нейтральной кривой, от течения Блазиуса ответвляется устойчивый периодический режим. При некотором α=α 2 передняя складка амплитудной поверхности из нефизической области отрицательных значений квадратов амплитуд выходит в область положительных значений и происходит смена закритической бифуркации докритической (α=α 3) . С увеличением числа Рейнольдса амплитуда этого режима нарастает, затем уменьшается и далее этот режим исчезает в результате слияния с неустойчивым режимом, ответвляющимся от верхней ветви линейной нейтральной кривой. Точки слияния этих режимов называются точками тангенциальной бифуркации [11], соответствующие точкам складки в теории катастроф [12]. Далее с увеличением α амплитудная поверхность периодических решений отрывается от линейной нейтральной кривой (α=α 4) и при некотором значении α периодические решения исчезают.

Рис.2 Рис.3
На рисунке 4 приведен срез амплитудной поверхности при α=0.206906, из которого видно, что при определенных начальных амплитудах периодические возмущения могут перейти либо в устойчивый периодический режим, либо в неустойчивый. На рисунке 5 приведены линейная нейтральная кривая (сплошная линия) и точки складок, или, что то же самое, точки тангенциальной бифуркации трехмерных режимов. Внутри нейтральной кривой возмущения нарастают, а вне – затухают.
В настоящее время экспериментально обнаружены три типа перехода. Первый был обнаружен в классической работе Клебанова и др. [13]. Он характеризуется быстрым нарастанием возмущений, высокочастотными всплесками и последующим образованием турбулентных пятен. Другой тип перехода был обнаружен в экспериментах Качанова, Козлова, Левченко [14]. Он характерен нарастанием высших гармоник двумерной волны Толлмина-Шлихтинга с последующим ростом трехмерной субгармоники. В экспериментах Козлова, Левченко, Сарика [15] наблюдался третий тип перехода – промежуточный между первыми двумя типами.
График амплитуды при «=0.206906

Рис.4 Рис.5
Для объяснения этих явлений в настоящее время существуют две основные теоретические модели. Первая была разработана Крейком [7] и по сути является обобщением метода Стюарта-Ватсона на случай взаимодействия трех резонансных волн. Другая модель – модель вторичной неустойчивости [6] – исследует устойчивость двумерных нелинейных волн относительно трехмерных линейных возмущений. Полученные с помощью этих моделей области быстрого роста трехмерных возмущений достаточно хорошо совпадают с экспериментальными результатами по наблюдению первых стохастических пульсаций. Но эти работы не объясняют причину появления стохастичности.
Как известно, в экспериментах по генерации ламинарно-турбулентного перехода с помощью вибрирующей ленточки клебановский тип перехода наблюдается только при относительно больших начальных амплитудах возмущения и характеризуется резким ростом возмущений. Таким образом клебановский тип перехода можно связать с передней верхней частью амплитудной поверхности, причем амплитуда соответствующего возмущения должна быть не меньше амплитуды периодического решения в передней точке тангенциальной бифуркации.
Если начальная амплитуда меньше критической клебановской, то возможен другой тип перехода, который можно соотнести с промежуточным типом, наблюдавшимся в эксперименте. Этот тип перехода связан с тем, что в результате нарастания возмущения становятся трехмерными и внутри нелинейной нейтральной поверхности для плоских возмущений их амплитуда может стать равной амплитуде точки тангенциальной бифуркации для трехмерных возмущений (β≠0).
И, наконец, при субгармоническом типе перехода в эксперименте наблюдается рост высших гармоник начальной двумерной волны Толлмина-Шлихтинга. Причем начальные амплитуды возмущений, вызывающих этот тип перехода, достаточно малы, и начальные возмущения действительно являются двумерными. Поэтому естественно связать субгармонический тип перехода с задней точкой тангенциальной бифуркации в окрестности верхней ветви линейной нейтральной кривой.
Таким образом, каждому типу перехода соответствует точка тангенциальной бифуркации периодических режимов (трехмерные или двумерные). Из теории динамических систем [16] следует, что в окрестности точки тангенциальной бифуркации возникает явление перемежаемости (т.е. чередование во времени ламинарного и турбулентного режимов). Причем длительность турбулентных режимов зависит от "расстояния" до точки тангенциальной бифуркации. Наблюдающаяся в экспериментах картина дает качественное подтверждение высказанным на основе проведенных расчетов предположениям.