Новый тип бифуркаций в модифицированной задаче о конвекции Рэлея-Бенара

Бифуркации ^^^^^^^^^^^^^^^^^

в динамических системах

Изв. вузов «ПНД», т. 21, № 2, 2013 УДК 532.516, 517.957

НОВЫЙ ТИП БИФУРКАЦИЙ В МОДИФИЦИРОВАННОЙ ЗАДАЧЕ О КОНВЕКЦИИ РЭЛЕЯ-БЕНАРА

И. А. Судаков, С. А. Вакуленко, Т. Г. Сукачева

Классическая конвекция Рэлея-Бенара является стандартным примером системы, в которой при изменении управляющего параметра возникают бифуркации. В данной статье рассматривается модифицированная задача о конвекции Рэлея-Бенара, включающая радиационные эффекты и источник газа на поверхности. В такой формулировке в указанной задаче возникает новый тип бифуркаций, помимо хорошо известных ячеек Бенара. Данная задача имеет большое значение для математической теории климата, так как обосновывает возникновение «критической точки» климатической системы, связанной с эмиссией парниковых газов в атмосферу Земли.

Ключевые слова: Конвекция Рэлея-Бенара, бифуркации, радиация, эмиссия газа, критическая точка, климат.

Введение

Явлению термогравитационной конвекции обязаны многие атмосферные и климатические процессы, поэтому развитие теории конвекции является важной задачей, широко востребованной в математической теории климата. Современная точка зрения на климатическую систему как на динамическую систему [1], а на климат -как на статистический ансамбль состояний, проходимых климатической системой за промежутки времени в несколько десятилетий [2], способствует применению различных методов математического анализа и математического моделирования для описания эволюции климатической системы.

Применение идей теории динамических систем к климатическим задачам началось после работ Э. Лоренца (1963) [3]. Исследуя конвективное движение газа в атмосфере, и применив метод Галеркина к системе типа Буссинеска, он ввел систему трех обыкновенных дифференциальных уравнений, известную теперь как система Лоренца. Далее им был изучен так называемый странный аттрактор Лоренца - компактное инвариантное множество Ь в трехмерном фазовом пространстве гладкого потока, которое имеет определённую сложную топологическую структуру и является асимптотически устойчивым [4]. Система Лоренца привлекла внимание крупных

специалистов по теории динамических систем, что способствовало её детальному изучению (например, Дж. Гукен-хеймер, 1976), однако строгое доказательство существования хаоса в этой модели было получено лишь в 1999 году В. Такером [5].

С 50-х годов XX века интерес к изучению климата методами теории динамических систем несколько упал вследствие перехода от аналитических моделей климата к достаточно сложным комплексным численным моделям всей земной системы. Из работ этого периода можно выделить работу Р. Гуди (1956), предложившего радиационно-конвективную модель атмосферы [6], устойчивость которой в последующем изучал В. Ларсон (1999) и подтвердил существование в ней ячеек Бенара [7] (см. рисунок).

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

В основополагающих работах [6, 8, 9] исследовался слой газа, который заключен между параллельными пластинами при заданных температурах, конвекция в котором происходит по типу Рэлея-Бенара. Кроме того, учтен тепловой перенос излучения, однако температура верхней пластины фиксирована. Отметим, что данная модель имеет некоторые преимущества над комплексными моделями атмосферы. Во-первых, работоспособность данной модели была проверена путем лабораторных экспериментов, о чем свидетельствует работа [8]. Во-вторых, результаты моделирования можно сравнить непосредственно с подобными моделями конвекции Рэлея-Бенара, выделяя таким образом эффект влияния теплового переноса излучения на конвективную устойчивость. В работах [10-12] отмечено, что введение в модель радиации оказывает два стабилизирующих воздействия на возникновение конвекции Рэлея-Бенара. Во-первых, радиационное затухание ведет к уменьшению малых температурных возмущений. Во-вторых, радиация приводит к более устойчивому температурному градиенту для основных состояний профиля температуры в исследуемой области.

Гуди [6, 8, 9] использовал вариационные методы и двухпотоковую радиационную модель, чтобы найти критические условия линейной устойчивости. Он рассматривал пределы оптически тонких и толстых газов и свободное скольжение, оптиче-

100001

Рис. Ячейки Бенара в гравитационном поле -возникновение упорядоченности в виде конвективных ячеек в форме правильных шестигранных структур в слое вязкой жидкости с вертикальным градиентом температуры

ски черные границы. В следующих работах рассматриваются более реалистичные системы, и делаются дальнейшие расчеты. В [13,14] был рассчитан полный спектр оптических толщин, доказана смена устойчивости для линейного состояния основного профиля температуры и представлен приблизительный критерий устойчивости в форме радиационного числа Рэлея. Некоторыми авторами [7] были изучены отдельные вклады основного состояния профиля температуры и радиационного затухания в конвективную устойчивость и сделаны расчеты для вертикального потока тепла при возникновении конвекции. Например, некоторые авторы рассчитывали граничные коэффициенты излучения. В [7,15] была изучена модель Гуди в пределе почти прозрачного газа. В этой работе получены два основных результата: во-первых, модель Гуди упрощена в пределе, когда коэффициент тепловой диффузии стремится к нулю, и было получено полное аналитическое решение задачи линейной стабильности. Линейный порог стабильности в точности определяется единственным параметром (радиационным числом Рэлея), который напоминает число Рэлея, используемое для характеристики конвекции Рэлея-Бенара [16-18]. Во-вторых, изучена нелинейная устойчивость решений для модели Гуди.

В настоящем исследовании предлагается распространить радиационно-кон-вективную модель атмосферы на случай существования постоянного динамического источника на земной поверхности, из которого эмитирует парниковый газ, и провести её бифуркационный анализ. В такой постановке это, прежде всего, позволит выявить новый тип бифуркаций в модифицированной задаче о конвекции Рэлея-Бенара, помимо ячеек Бенара. Кроме того, данное исследование имеет важное прикладное значение: модифицированная модель может быть применена для изучения эмиссии метана из вечной мерзлоты в атмосферу [19]. Метан - опасный парниковый газ; его влияние на дополнительный разогрев атмосферы в двадцать один раз больше, чем у углекислого газа. Огромный метановый пул захоронен в виде газогидратов в вечной мерзлоте. В связи с изменением климата вечная мерзлота оттаивает, и метан поступает в атмосферу, где влияет на ее разогрев. Это еще больше усиливает оттаивание мерзлоты, благодаря явлению положительной обратной связи. Таким образом, это явление ведет климатическую систему к прохождению точки бифуркации, идентификация которой является на сегодняшней день достаточно актуальной задачей.

1. Классическая модель Гуди и её обобщение

Модель, которая исследована в настоящей работе, рассматривает атмосферу в области типа бесконечной пластины, ограниченной верхней и нижней твердыми, свободными от скольжения границами. Ненулевой адиабатический градиент скорости Г* = д*/ер^ (здесь обозначаем звездочками размерные величины) входит в уравнения, так как вертикальный градиент имеет большое значение для моделирования процессов в атмосфере. Тепловое излучение поглощается, но не рассеивается в газе. Предположим, что рассматриваемая система - радиационно-серое тело, то есть ее оптические свойства считаются независимыми от длины волны излучения. Рассмотрим только прозрачные слои газа, то есть слои газа с малой оптической глубиной.

Следуя [7,15], запишем в аппроксимации Обербека-Буссинеска уравнение движения:

v*í + =--V*/* + #*ае,9*г + ^У^*,

р*

где V* - скорость, 9** - отклонения температуры от некоторого основного уровня 9*, Р* - давление, Ь* - время, р* - постоянная, относящаяся к плотности, ае* - постоянный коэффициент объемного расширения, г - компонента вертикального вектора, V* - кинематическая вязкость и д* - гравитационная постоянная. Уравнение теплопроводности:

9*4 +

Р*ср*

Здесь 9 = 9* + 9* - общая температура, полученная путем сложения температуры в равновесном состоянии и температуры после перехода из этого состояния, ю* - вертикальная компонента скорости V*, Е* - радиационный поток, ср* - удельная теплоемкость при постоянном давлении, адиабатический градиент скорости Г* = = д*/ср*, к* - коэффициент теплопроводности.

В приближении Обербека-Буссинеска уравнение неразрывности имеет вид

У^* = 0.

Следуя [15], введем новые безразмерные переменные и перепишем уравнения (1) (3) для модели Гуди в безразмерной стандартной форме в приближении Обербека-Буссинеска с учетом действия внешней периодической силы и записи некоторых физических констант для удобства проведения последующих вычислений

Vt +

е4 +

Vv = 0.

Здесь V = (и(х, у, г, Ь)^(х, у, г, Ь), ю(х, у, г, Ь)) - скорость потока газа с вертикальной компонентой ю, 9(х, у, г, Ь) - температура потока, Р - давление, С - концентрация эмитирующего газа (фактически, примесей), Г - безразмерный адиабатический градиент, о - число Прандтля, 01 - параметр плавучести, z - единичный вектор в вертикальном направлении, ©о = 9т* - основное температурное состояние 9, вычисленное в средних точках слоя, f - внешняя сила, у - коэффициент, подобный числу Рэлея, который умножается на член, описывающий плавучесть в уравнении движения, Е - радиационный поток и ^ - тепловой источник. Для простоты предположим, что движение происходит в тонком слое определенном как 0 < г < Н, (х,у) € П, здесь П - плоский прямоугольник с размерами 11, 12.

Уравнение радиационного потока Е может быть получено путем аппроксимации Эддингтона для уравнения Шварцшильда [20-23]. В результате

У

где предполагается, что а является достаточно малой постоянной, пропорциональной коэффициенту поглощения излучения в единицу объема. Тогда (7) может быть упрощено

-УЕ = 3а6.

В самом деле, подстановка (8) в (7) дает, что (7) удовлетворяется с точностью до слагаемых порядка о(а2). Тогда (5) примет вид

е4 + ^У)е + юГ = де - 3а6 + д.

Это уравнение уже не содержит радиационный поток. Таким образом, основные уравнения для модели даются формулами

С +

где С(х, у, г, -) является концентрацией эмитирующего газа, Ь0 - положительная постоянная (член —Ь0С описывает деградацию молекул газа в атмосфере в результате химических реакций), с! - коэффициент диффузии эмитирующего газа.

Уравнение (10) может быть решено, если известна скорость потока в сплошной среде V и уравнения в, V и С независимы [24,25]. Для описания влияния газа на циркуляцию в рассматриваемой системе, предполагается, что коэффициент а в (9) зависит от С: а = а(С). Так как концентрация газа С очень мала (исходя из физических соображений), естественно предположить, что эта зависимость линейная. Тогда, используя ряд Тейлора для а, получим

а & а0 + а1С,

где ао, а1 - постоянные; а1 - феноменологический коэффициент, который может быть определен путем анализа экспериментальных данных. Тогда уравнения (4), (6), (9) и (10) становятся взаимосвязанными и, следовательно, эта система может описывать обратную связь между сплошной средой (скорость потока, температурные поля) и концентрацией газа.

Рассмотрим граничные условия. Условия для V - это стандартные условия отсутствия скольжения

v

Граничные условия для температуры е имеют вид

и означают отсутствие потока на границе.

Граничные условия для С описывают отсутствие потока газа при г = Н и задают поток газа с поверхности г = 0:

Сг (х,У,г,Щг^ = 0

С г (х, у, г, -)|г=0 = -ц(х,у, е(х,у, 0,-)),

где ^-функция, которая описывает интенсивность потока газа с поверхности. Предполагается, что этот поток, вообще говоря, неоднороден в пространстве и зависит от температуры. Для упрощения математического анализа используем периодические граничные условия по ж, у. Если П = {(ж, у) :0 <х <11, 0 <у < 12}, следовательно, имеем:

v(ж,y, г, г) = v(ж + 11, у, г, г) = v(x, у + 12,г,г),

Р

0(х, у, г, г) = 0(ж + 11 ,у, г, г) = 0(ж + 12 ,у, г, г).

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

=д

где У0 является потенциальным барьером для химической реакции образования газа, кв - постоянная Больцмана. Здесь будем считать, что газ с поверхности поступает не сплошным потоком, а из динамических отверстий (подобно тому, как метан поступает в атмосферу из мерзлотных озер [26]), тогда можно предположить, что

N

9(ж,у,г) = сдXVОi(t),

г=1

где сд - константа, Ог - двумерная область, занятая г-м отверстием, и XV является характеристической функцией множества V: XV = 1, если (ж, у) € V и XV = 0 в противном случае. Обратим внимание, что области Ог могут медленно меняться, тогда необходимо рассматривать случай, когда д зависит от времени. Основываясь на физических соображениях, будем считать, что изменение формы отверстий (по аналогии с мерзлотными озерами) идет медленно [26]: Ог = Ог(т), где т является медленным временем т = е^, £1 ^ 1, и возможно предположить, что в первом приближении Ог - фиксированные области.

2. Существование и единственность решения краевой задачи

Опишем общий подход, который позволяет доказать существование и единственность решения для системы дифференциальных уравнений (4), (6), (9)-(11) с краевыми условиями (12)-(16). Этот подход хорошо известен - см. [27,28].

Функциональные пространства. Определим некоторые функциональные пространства. Обозначим через Ва,р стандартные интерполяционные пространства

где Ад - оператор Лапласа при нулевых граничных условиях Дирихле и символ ||...|| обозначает норму функции. Интерполяционные пространства Ва,р тесно связаны со стандартными пространствами Соболева. Их теория подробно изложена в [27,29].

Краевая задача как абстрактное эволюционное уравнение. Обозначим через u тройку (v,6,C) и Va,p - пространство векторных полей (vi,v2,v3). Предполагается, что каждый компонент лежит в соответствующем функциональном пространстве: v € Va,p и 6, C лежат в Ba,p. Для реализации стандартной схемы доказательства существования решений уравнения (4) скорости потока будем применять определенный для данной задачи оператор Лере P [27, 30]. Обратим внимание, что PVP = 0, поэтому этот метод исключает давление [27,31]. В результате можно переписать всю систему как эволюционное уравнение:

щ = Au + F

где A - некоторый самосопряженный отрицательно определенный оператор A = (oPA, A, dA)tr, F определяет другие нелинейные и линейные члены, которые не входят в A, tr(u) - след функции на нижней границе s, подчиняющийся следующему вложению [31]: ||tr(u)||ip(s) < c ||u||a,p, a > 1/2p.

Локальное существование. Чтобы показать локальное существование решений во времени, применим стандартный подход [27]. Выберем такие p, a, что

Cs(Q) С Bap, 0 < s < 2(a - 1/p).

Тогда sup |V6| < C ||6|| , если 9 € Ba,p. Это показывает, что отображение 6 ^ (vV)9 является C1-отображением из Ba,p в Lp. Аналогично, для отображения C ^ (vV)C. Таким образом, F является C1-отображением из Va,p х в Vo,p х Lp. Это доказывает существование локальных (по времени) решений (17).

Априорные оценки. Описанная выше конструкция позволяет доказать существование решений краевой задачи на ограниченном интервале времени. Чтобы показать, что на самом деле решения краевой задачи существуют для всех t, необходимо получить априорную оценку решений. Допустим, что имеет место

Лемма. Пусть функция ^ > 0, ^ € C1 ограничена в 6:

sup x,y, 6) < c0 < то

x,yeQ,6e

и внешние источники и силы Q, f € C0определены так, что

sup Q(x, y, z, t) = q0 < то, Q > 0,

x,y€Q,t>0

sup |f(x,y,z,t)l = f0 < то.

x,yeQ,t>0

Тогда поля концентрации эмитирующего газа C и температуры 6 равномерно ограничены.

Доказательство. Учитывая, что supx y,z 6(x, y, z, 0) = M0, при помощи верхних и нижних решений [32] можно легко доказать, что

0 < 6(x,y,z,t) < max {M0,q0(3a0)-1} = Mi и 0 <C(x,y,z,t) < c0z,

если начальные поля в, С(х, у, г, -) положительны при всех (х, у, г) € Если температура в равномерно ограничена, тогда уравнение для скорости потока дает стандартную оценку

2

—У—<Г —Л

что в соответствии с рассуждениями [32] доказывает, что норма |^(-)|| ограничена по

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

Теорема. Пусть выполнены условия Леммы. Тогда решение краевой задачи (4)-(16) существует и единственно при всех - € (0, то) и лежит в пространстве В1рБа>рБа>р (то есть V- > 0, в € Ба,р, С € Ба,р, у^Ба,р, ] = 1, 2, 3), а > 1/р.

Соотношение для средней температуры. Сначала получим некоторые интегральные соотношения для средней температуры ((в)), где((/))обозначает значения функции /, усредненной по всей области ((/)) = |^|-1 и |П| = - объем элементарной периодической ячейки. Проинтегрируем уравнение теплопроводности (9) по Предполагая, что нет источников на границах, проинтегрируем по частям, и, учитывая граничные условия, получаем

^Ш = _з(((ао + а1с)е>>+ «<?>>.

Тогда вклад эмитирующего газа определяется как

Яшеь = -3а1

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

3. Бифуркационный анализ обобщенной модели Гуди

Здесь возможны несколько различных ситуаций. Для их описания введем безразмерный параметр

квбо'

который определяет типичные скорости эмиссии газа, когда выполняется соотношение типа Аррениуса. Здесь 9о является температурой на поверхности. Будем считать, что ai = е - малый параметр в теории возмущений, что позволяет оценить Rmet. Существует три основных возможных случая.

I. Параметр k ограничен, и система находится далеко от точки бифуркации.

II. Параметр k ограничен, и система близка к точке бифуркации.

III. Параметр k имеет достаточно большие значения, и система близка к точке бифуркации.

Случай I (система вдали от точки бифуркации). Фактически, цель этого раздела состоит в построении теории возмущений, что позволяет получить систему уравнений для некоторого X и вычислить X, оценивая влияние эмитирующего газа на динамику системы. Это подход из теории возмущений, основан на том, что концентрация этого газа мала (как и в реальных физических ситуациях).

Пусть отверстия, из которых эмитирует газ, имеют относительно малый диаметр VI по отношению к размерам 1\, 12 области П. Таким образом, можно предположить, что температура 0(ж, у, 0, Ь) является постоянной внутри всей области отверстия. Обозначим усредненное значение через 0(ж,у, 0, Ь) для г-го отверстия (Ь), тогда формально получим

Wi = area(Qi) 1 Ц 0(x,y, 0,t)%Qi (x,y)dxdy.

Q

Кроме того, представим неизвестные функции в виде суммы основного вклада (переменные с индексом «0») и малых поправок, тогда основные условия для vo могут быть найдены из (4), (6) и (9) и из уравнения переноса тепла при а = ао:

004 + (^У)0о + юоГ = Д0о - 3ао0о-

Следующим шагом является решение краевой задачи для С:

С + КУ)С = йДС - ьоС,

С

( V \ N

Cz

\z=о

\ Л L> W Л' I I I /

i=l

Явные решения (19)-(21) можно найти только в некоторых частных случаях. Рассмотрим случай, когда невозмущенный поток газа стационарен во времени и внешняя сила отсутствует. Тогда вспомогательная краевая задача запишется так:

Ои + (v0V)Gi = dAGi - b20Gi,

GiZ (x,y,z,t)\z=h =

Giz (x,y,z,t)\z=0 = S(t)XQi (x,y)-Здесь решение для C может быть выражено через Gj:

N 1

С(х, y,z,t) = J2 J y,z,t- t') ехр(-—(22) i=l 0

Подставляя соотношение (22) в (18), получаем

N t

N ^ r Vo

^met(i) = -3ai Е J {{Gi[x,y,z,t -t')Qo(x,y,z,t)))exp(-—

i=l 0

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

Случай II. Система около точки бифуркации. Фундаментальные вариационные соотношения, полученные В.С. Сорокиным в 1954 году [33], а также в работах [3436], могут быть применены в данном исследовании, что позволяет оценить влияние эмиссии газа с поверхности на бифуркацию в рассматриваемой системе. Положим Е, ^ = 0 и удалим нелинейные члены из основных уравнений Навье-Стокса (4), (6) и (9). Постулируем, что основное состояние системы дано соотношениями vo = 0, температура во = Т(г) = Ьг + а, радиационный поток Е = Г(г^, и радиационный коэффициент поглощения а является постоянным, а = ао. Получим уравнения, линеаризованные вблизи основного решения

vt = о^ - УР + о^,

еt = Дв + ш

Vv = 0.

При этих предположениях скорость потока и тепловых полей можно определить с точностью до малых значений при помощи нормальных мод [15, 34], при этом V(x, у, г), в = Т(х, у, г)являются собственными функциями некоторого линейного оператора Ь [7, 34], действующего на малые поправки v, в. Если рассмотреть линеаризованные уравнения (предполагаем, что внешняя сила равна нулю, тепловой источник тоже ноль и члены а1С пренебрежимо малы), то X(-) = X(0)ехр(-А,-), где X (-) является неизвестной функцией времени, определяющей временную эволюцию, А - в общем случае, комплексное собственное значение А = в + ш, где г = \/— I. Точки бифуркации определяются условием |3 = 0.

Теперь можно получить некоторые важные вариационные отношения для А, повторяя известные расчеты В.С. Сорокина [33, 34] применительно к (23)-(25), и тогда для скорости потока V = (V!, V2, Vз) и температуры Т имеют место соотношения:

-XV = оДV -УР + 0]_Тг,

-АТ = ДТ + Уэ

Vv = 0.

Умножим (26) на комплексно-сопряженную V* и (27) на Т*. Тогда получим два соотношения, обобщающие известные соотношения для конвекции Рэлея-Бенара [33]:

-А / |V|2= -о / ^ V|2dv + 01 / TzV*dv, Уа Уа Уа

-А I T2dv = - / (|УТ|2 + 3аT2)dv + у' / Т*zVdv, а а а

где = (1х(1ус1г и у' = у-Г. Комбинируя эти соотношения, при некоторых условиях получаем 1т А = 0 (следовательно, здесь нет колебаний по времени). Далее, следуя

В.С. Сорокину [33], можно получить, что

Но = 3у'-1о^1 / С |Т|2

./а

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

Величина Но может быть вычислена в ряде случаев, когда возможно получить некоторые выражения для С. При условии, что конвективный член мал, можно это выражение найти из уравнения

dЛC -Ъ0С = 0

при граничном условии

С (х,у,г,г)\г=0 = ц(х,у, во),

где функция (х определена в (16).

Если параметр Ьо мал по отношению к \Л/, то концентрация газа достаточно однородна и возможно получить новое соотношение для

Но = 3у/-101а1 «С»/ |Т|2 dv,

а

где С является усредненной концентрацией примесей газа в системе. Для малых возмущений 6т^Х параметра релаксации X, связанных с воздействием эмитирующего газа, получаем

-8те1Х = НоК-1,К = / (|У|2 - У'-Ч |Т|2^-1.

а

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

Случай III. Система около точки бифуркации: взрывообразные выбросы. Третий случай описывает ситуацию, когда переходные процессы связаны с эмиссией газа, которая происходит взрывообразно с очень большой скоростью. Например, в реальности, для некоторых температур

р(х,у,г) = д(х,у)Б (во в),

где и0 = в0в, Б - монотонно возрастающая функция, такая что Б(и0) = 1

и Б(ио) = 0.

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

Б(и0) = 0, и0 < 0, Б = и0, 0 <и0 < 1, Б(и0) = 1 и0 > 1.

Получим здесь уравнения для амплитуд мод X^(¿), г = 1, ...Ы, где N представляет собой число бифуркационных мод, участвующих в медленной динамике, в точке бифуркации. Например, для случая II, такие уравнения достаточно просты и напоминают известные модели Лоренца [3]. Кроме того, для случая III и N = 1 уравнение эволюции является относительно простым. Наиболее интересный случай - это случай III с N > 1, где уравнения могут быть более сложными. Остановимся на нем более детально. Рассмотрим случай, когда N > 1 - некоторые основные режимы могут сосуществовать. Представим здесь температуры и скорость потока как

N

в = в0 + Е Хп(1)Тп (х,у,г),

п=1

N

V = Е Хп (¿)Уп (х,у,г),

п=1

где \п, Тп - заданные функции (которые обычно выбирают в качестве собственных функций линейной части), Хп являются неизвестными функциями времени. Рассмотрим сначала обычный случай для а^ Тогда, применяя стандартную процедуру Галеркина, получаем

^ = — ЬпХп + /га

где М - матрица, А,п > 0 - коэффициенты, ¡п - определяет внешний эффект силы, -\пХп является диссипативным членом, а первый член нелинейно связан с конвекцией. Если Мпч антисимметрична для п, г, то (30) является случаем известной системы Колмогорова. К тому же, известную систему Лоренца можно рассматривать как частный случай (30).

В рассматриваем случае, получаем (учитывая, что 3а 1С6 достаточно мало):

<1Хп и и

Х3 + ¡п + а1 Р„(Х),

г,3=1 3 = 1

где Л является положительной матрицей с элементом А,п,, который отвечает за диссипацию в системе,

N

Рп = £ ВпзХзБ(Рз(X)). 3=1

N

Здесь Рз(X) = 60(Хг,уз,Ь) + ^ ^азтХт), а а,т,Бп, - матрицы, (х,,у,) - ко-

т=1 т

ординаты, которые определяют центр j-го отверстия. Теперь опишем, как получить вклад Р. Применяя стандартную процедуру Галеркина, имеем

¥п = 3 (С (Х) 6,Тп) ,

где (/, д) обозначает стандартное скалярное произведения в ¿2^). Чтобы найти С, необходимо использовать приближения, где можно пренебречь эффектами, связанными с конвекцией и распадом примесей газа в уравнении (19) для С. Тогда можно искать для плотности квазистационарных потоков газа решение краевой задачи

дС = 0,

С

С

Тогда решение (31)-(33) для С может быть выражено через функцию Грина, и здесь получается (с учетом того, что отверстия достаточно малы по размерам) следующее соотношение:

Макез N

С (х, у, X, Х) ~ ^ Ог(х, у, X, Хг ,уг)Б(6о + ^ Т, (хг ,уг)Х, ).

г=1 3

В результате,

а,т — Т3

N¡^68

Бпт ^ У ] (С%(х,у,х,х% ,уг)Тт,Тп).

г=1

Рассмотрим теперь возможный механизм хаоса, связанного с взрывообразными выбросами газа с поверхности. Даже для а1 = 0, то есть без эмитирующего газа, влияние динамики может быть сложным (частный случай этих систем - система Лоренца, может порождать динамику со сложными бифуркациями и хаосом). Рассмотрим некоторые дополнительные хаотические эффекты, связанные с учетом эмиссии газа

в модели Гуди. Эти эффекты возможны даже в случае II. Удалим все нелинейные члены из (30), кроме членов, отвечающих за эмиссию газа, а также удалим член с ajmXm. Получим систему

dXn N

J, = ~ ^njXj + fn + a\EnjXjS

j=i

Исследуем теперь динамику системы (34). Будем рассматривать nj = = S(0o(xj,yj, t)) как сезонный длительный фактор времени (так как в реальном случае озера, из которых эмитирует метан, то тают, то замерзают), включающийся или выключающийся в случайные моменты времени (из-за сигмоидального характера S). Динамика системы (34) при ai = 0 вполне тривиальна, поскольку все траектории сходятся к состоянию равновесия, а динамика системы при ai =0 может генерировать фрактальные аттракторы.

Чтобы показать это, рассмотрим отображение Пуанкаре P : X (t) ^ X (t+AT) для (34), где AT является длительным промежутком времени (например, год). С учетом рассуждений, сделанных выше, ясно, что за год интервал (nAT, (n + 1)AT) этого отображения P определяется линейным отображением

PnX =

где A есть фиксированная матрица (A = exp(-ЛАТ)), En- случайная матрица в зависимости от n. Эта случайность возникает из-за сезонного фактора nj = = S(60(xj,yj, t)). Отображение, подобное (35), детально изучено в [31]. Оно определяет так называемые системы итерируемых функций (IFS) [37]. К тому же, хорошо известны аттракторы этих случайных отображений, в общем случае это - фрактальные множества.

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

Заключение

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

распаду (фактически, такое поведение в атмосфере наблюдается у любого парникового газа).

Далее была сформулирована краевая задача для такой системы и показано существование и единственность решений этой задачи на основе построения оператора Лере и метода верхних и нижних решений.

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

Таким образом, применение данной задачи в математической теории климата позволит изучить «критические точки» (точки бифуркации) в климатической системе, что является, несомненно, важным для прогнозирования будущих изменений климата.

Работа выполнена в рамках государственного задания НовГУ (проект «Исследование моделей, методов и алгоритмов многопараметрических систем», № 438/АиГ-с (2011-2013)), программы стратегического развития НовГУ (проект «Исследование обобщенных решений нелинейных интегро-дифференциальных и разностных уравнений», № 456/АиГ-с (2012-2016)) и при поддержке РФФИ (грант 12-01-90803-мол _рф_нр).

Библиографический список

1. Thompson J.M. T., Sieber J. Predicting climate tipping as a noisy bifurcation: A review // Int J. Bif. Chaos. 2011. Vol. 21, № 2. P. 399.

2. Монин А.С., Шишков Ю.А. Климат как проблема физики // Успехи физических наук. 2000. № 4. С. 419.

3. Lorenz E.N. Deterministic nonperiodic flow // Journal of the Atmospheric Sciences. 1963. Vol. 20, № 2. P. 130.

4. Чуличков А.И. Математические модели нелинейной динамики. М.: Физматлит, 2000.

5. Tucker W. The Lorenz attractor exists // C. R. Acad. Sci. Paris. 1999. Vol. 328. P. 1197.

6. Goody R.M. The influence of radiative transfer on cellular convection // J. Fluid Mech. 1956. V. 1. P. 424.

7. Larson V.E. The effects of thermal radiation on dry convective instability // Dynamics of Atmospheres and Oceans. 2001. Vol. 34. P. 45.

8. Goody R.M. Atmospheric Radiation. I. Theoretical Basis. Oxford University Press, NewYork, 1964.

9. Goody R.M. Corrigendum // J. Fluid Mech. 1956. Vol. 1. P. 670.

10. Goody R.M., Yung Y.L. Atmospheric Radiation: Theoretical Basis, 2nd Edition. Oxford University Press, New York, 1989.

11. Gille J., Goody R.M. Convection in a radiating gas // J. Fluid Mech. 1964. Vol. 20. P. 47.

12. Goody R.M. Principles of Atmospheric Physics and Chemistry. Oxford University Press, New York, 1995.

13. Spiegel E.A. The convective instability of a radiating fluid layer // Astrophys. J.V. 1960. Vol. 132. P. 716.

14. Spiegel E.A., Veronis G. On the Boussinesq approximation for a compressible fluid //Astrophys. J. 1960. Vol. 131. P. 442.

15. Larson V.E. Stability properties of and scaling laws for a dry radiative-convective atmosphere // Q. J. R. Meteorol. Soc. 2000. Vol. 126. P. 145.

16. Murgai M.P., Khosla P.K. A study of the combined effect of thermal radiative transfer and a magnetic field on the gravitational convection of an ionized fluid // J. Fluid Mech. 1962. Vol. 14. P. 433.

17. Narasimha R., Vasudeva Murthy A.S. The energy balance in the Ramdas layer // Bound. Layer Meteorol. 1995. Vol. 76. P. 307.

18. Vasudeva Murthy A.S, Srinivasan J., Narasimha R. A theory of the lifted temperature minimum on calm clear nights // Phil. Trans. R. Soc. London A. 1993. Vol. 344. P. 183.

19. Sudakov I., Vakulenko S. Bifurcations of the climate system and greenhouse gas emissions // Philos. Trans. A Math. Phys. Eng. Sci. 2013: 371(1991):20110473.

20. Bledoui F., Soufani A. The onset of Rayleigh-Benard instability in molecular radiating gases // Phys. Fluids A. 1997. Vol. 9. P. 3858.

21. Getling A. V. On the scales of convection flows in a horizontal layer with radiative energy transfer // Atmos. Oceanic Phys. 1980. Vol. 16. P. 63.

22. Veronis G. Penetrative convection //Astrophys. J. 1963. Vol. 137. P. 641.

23. Vincenti W.G., Traugott S.C. The coupling of radiative transfer and gas motion // Annu. Rev. Fluid Mech. 1971. Vol. 3. P. 89.

24. Полежаев В.И., Яремчук В.П. Численное моделирование двумерной нестационарной конвекции в горизонтальной слое конечной длины, подогреваемом снизу // Механика жидкости и газа. 2001. № 4. C. 34.

25. Joseph D.D. Nonlinear stability of the Boussinesq equations by the method of energy // Arch. Ration. Mech. Anal. 1965. Vol. 22. P.163.

26. Судаков И.А. Динамика протаивания мерзлотных озер и изменения климата // Научно-технические ведомости Санкт-Петербургского государственного поли-

технического университета. Физико-математические науки. 2011. № 2. С. 74.

27. Генри Д. Геометрическая теория полулинейных параболических уравнений. М.: Мир, 1984.

28. Drazin P.G., Reid W.H. Hydrodynamic Stability. Cambridge University Press, New York, 1981.

29. Бесов О.В., Ильин В.П., Никольский С.М. Интегральные представления функций и теоремы вложения. М.: Наука, 1975.

30. Ладыженская О.А. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Наука, 1970.

31. Стейн И.М. Сингулярные интегралы и дифференциальные свойства функций. М.: Мир, 1973.

32. Фридман А. Уравнения с частными производными параболического типа. М.: Мир, 1968.

33. Сорокин В.С. О стационарных движениях жидкости, подогреваемой снизу // ПММ. 1954. № 18, вып. 2. C. 197.

34. Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972.

35. Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость // Итоги науки и техники. Серия «Механика жидкости и газа». 1978. Т. 11. С. 66.

36. Straughan B. The Energy Method, Stability, and Nonlinear Convection. Springer, New York, 1992.

37. Barnsley M.F., Demko S. Iterated function systems and the global construction of fractals// Proc. Roy. Soc. London Ser. A. 1978. Vol. 399. 1817. P. 243.

Университет Юты, США Поступила в редакцию 16.10.2012

Институт проблем машиноведения РАН После доработки 1.05.2013

Новгородский госуниверситет

NEW TYPE OF BIFURCATIONS IN THE MODIFIED RAYLEIGH-BENARD

CONVECTION PROBLEM

I. A. Sudakov, S. A. Vakulenko, T. G. Sukacheva

The original Rayleigh-Benard convection is a standard example of the system where bifurcations occur with changing of a control parameter. In this paper we consider the modified Rayleigh-Benard convection problem including radiative effects as well as gas sources on a surface. Such formulation leads to the identification of new type of bifurcations in the problem besides the well-known Benard cells. This problem is very important for mathematics of climate because it proves the occurrence of the climate system tipping point related to greenhouse gas emission into the atmosphere.

Keywords: Rayleigh-Benard convection, bifurcations, radiation, gas emission, tipping point, climate.

Судаков Иван Алексеевич - родился в Свердловской области (1985), окончил физический факультет Уральского государственного университета им. А.М. Горького (2008). После окончания университета обучался в аспирантуре Санкт-Петербургского государственного университета по специальности «Климатология», а также Новгородского государственного университета им. Ярослава Мудрого по специальности «Математическое моделирование». Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в НовГУ (2012) в области прикладной математики. В настоящее время занимает должность профессора, ассистирующего в исследованиях, математического департамента Университета Юты (США). Награжден именными стипендиями: имени Фритьофа Нансена (2008), имени Отто Шмидта (2011), имени Эдварда Лоренца (2012). Лауреат премии правительства СПб в области научно-педагогической деятельности за 2010 год.

University of Utah, 155 S 1400 E, Room 233 Salt Lake City, UT 84112-0090 USA E-mail: sudakov@math.utah.edu

Вакуленко Сергей Августович - родился в Ленинграде (1953), окончил математико-механический факультет Ленинградского государственного университета (1976), защитил кандидатскую диссертацию в 1981, докторскую диссертацию (по теоретической и математической физике) в 1993 в Санкт-Петербургском государственном университете. Являлся приглашенным профессором в Лионском Университете (1998-2003), университетах г. Ренна и г. Монпелье (2004-2012). С 2006-2011 заведовал кафедрой высшей математики и информатики Санкт-Петербургского государственного университета технологии и дизайна. Соавтор двух книг и автор свыше 70 работ по математической и теоретической физике, математической биологии и механике, опубликованных в ведущих журналах. Был одним из ведущих организаторов международной конференции «Patterns and Waves» (Пушкин, 2002). В настоящее время работает в должности ведущего научного сотрудника Института проблем машиноведения РАН.

199178 Санкт-Петербург, Большой проспект В.О., д. 61 Институт проблем машиноведения РАН E-mail: vakulenfr@mail.ru

Сукачева Тамара Геннадьевна - родилась в Новгороде (1955), окончила физико-математический факультет Новгородского государственного педагогического института (1977). Защитила диссертацию на соискание ученой степени кандидата физико-математических наук по специальности «дифференциальные уравнения» в Воронежском государственном университете (1990) и доктора физико-математических наук по специальности «математическое моделирование, численные методы и комплексы программ» в Новгородском государственном университете (2004). Автор монографии «Расширенные фазовые пространства моделей Осколкова». Имеет более 100 различных публикаций. В настоящее время работает в должности профессора кафедры алгебры и геометрии Института Электронных и Информационных Систем Новгородского государственного университета имени Ярослава Мудрого.

173003 Великий Новгород, ул. Большая Санкт-Петербургская, д. 41 Новгородский государственный университет им. Ярослава Мудрого, Институт электронных и информационных систем E-mail: tamara.sukacheva@novsu.ru