Научная статья на тему 'Об одном классе аналитических решений стационарной осесимметричной конвекции Бенара–Марангони вязкой несжимаемой жидкости'

Об одном классе аналитических решений стационарной осесимметричной конвекции Бенара–Марангони вязкой несжимаемой жидкости Текст научной статьи по специальности «Математика»

CC BY
211
80
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОСЕСИММЕТРИЧНАЯ ТЕРМОКАПИЛЛЯРНАЯ КОНВЕКЦИЯ (КОНВЕКЦИЯ БЕНАРА–МАРАНГОНИ) / ЛОКАЛИЗОВАННЫЙ ПАРАБОЛИЧЕСКИЙ НАГРЕВ / ТОЧНОЕ РЕШЕНИЕ / ИЗОЛИНИИ / МАТРИЦА ГЕССЕ / СОБСТВЕННЫЕ ЧИСЛА / ЛОКАЛИЗАЦИЯ КОРНЕЙ ПОЛИНОМОВ / ЛОКАЛИЗАЦИЯ СОБСТВЕННЫХ ЧИСЕЛ МАТРИЦЫ / AXISYMMETRIC THERMOCAPILLARY CONVECTION (CONVECTION BéNARD–MARANGONI) / LOCALIZED PARABOLIC HEATERS / EXACT SOLUTION / ISOLINES / HESSIAN MATRIX / EIGENVALUES / LOCALIZATION OF POLYNOMIALS ROOTS / LOCALIZATION OF EIGENVALUES OF THE MATRIX

Аннотация научной статьи по математике, автор научной работы — Аристов Сергей Николаевич, Просвиряков Евгений Юрьевич

Целью настоящей работы является нахождение решений системы уравнений Обербека–Буссинеска плоской конвекции Бенара–Марангони вязкой несжимаемой жидкости, в которых радиальная составляющая градиента температуры может обратиться в нуль. Показано, что исходная система может быть сведена к системе обыкновенных дифференциальных уравнений одиннадцатого порядка. Получено точное решение в точке экстремума температуры (при нулевом числе Грасгофа). Интегрирование уравнений осуществлено в безразмерных переменных, которые введены неклассическим способом: введен характерный масштаб по каждой переменной, а не по линейному характерному размеру слоя. Найденное решение служит начальным приближением для построения решения конвекции Бенара–Марангони при числах Грасгофа, больших, чем нуль.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Аристов Сергей Николаевич, Просвиряков Евгений Юрьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

On one class of analytic solutions of the stationary axisymmetric convection Bénard–Maragoni viscous incompreeible fluid

The purpose of this work is to find solutions for the system of equations Oberbeck– Boussinesq flat convection Bénard–Marangoni a viscous incompressible fluid. In this viscous incompressible fluid the radial component of the temperature gradient may become zero. It is shown that the initial system may be reduced to the system of equations of ordinary differential equations of the eleventh order. We obtain the exact solution at the point of the extremum of the temperature (at zero including Grasgof’s). Integration of equations is carried out in dimensionless variables, which are non-classical way: put the scale factor for each variable, and not by linear characteristic size of the layer. The solution is the initial approximation to the solution of convection Bénard–Marangoni in numbers Grasgof’s, the big zero.

Текст научной работы на тему «Об одном классе аналитических решений стационарной осесимметричной конвекции Бенара–Марангони вязкой несжимаемой жидкости»

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2013. № 3 (32). С. 110—118

УДК 517.958:532.51

ОБ ОДНОМ КЛАССЕ АНАЛИТИЧЕСКИХ РЕШЕНИЙ СТАЦИОНАРНОЙ ОСЕСИММЕТРИЧНОЙ КОНВЕКЦИИ БЕНАРА-МАРАНГОНИ ВЯЗКОЙ НЕСЖИМАЕМОЙ ЖИДКОСТИ

С. Н. Аристов1, Е. Ю. Просвиряков2

1 Институт механики сплошных сред УрО РАН, Россия, 614013, Пермь, ул. Акад. Королёва, 1.

2 Казанский государственный национальный исследовательский университет им. А. Н. Туполева, Россия, 420111, Казань, ул. Карла Маркса, 10.

E-mails: asn@icmm.ru, evgen_pros@mail.ru

Целью настоящей 'работы является нахождение 'решений системы уравнений Обербека—Буссинеска плоской конвекции Бенара—Марангони вязкой несжимаемой жидкости, в которых радиальная составляющая градиента температуры может обратиться в нуль. Показано, что исходная система может быть сведена к системе обыкновенных дифференциальных уравнений одиннадцатого порядка. Получено точное решение в точке экстремума температуры (при нулевом числе Грасгофа). Интегрирование уравнений осуществлено в безразмерных переменных, которые введены неклассическим способом: введен характерный масштаб по каждой переменной, а не по линейному характерному размеру слоя. Найденное решение служит начальным приближением для построения решения конвекции Бенара—Марангони при числах Грасгофа, больших, чем нуль.

Ключевые слова: осесимметричная термокапиллярная конвекция (конвекция Бенара—Марангони), локализованный параболический нагрев, точное решение, изолинии, матрица Гессе, собственные числа, локализация корней полиномов, локализация собственных чисел матрицы.

Введение. Не вызывает сомнений, что конвективное течение жидкости, связанное с неоднородным нагревом, является самым распространенным во Вселенной [1]. Систематические исследования конвективных течений вязкой несжимаемой жидкости датируются тысяча восемьсот восемьдесят восьмым годом, когда Томсон наблюдал образование сотовых структур в мыльной воде и мясном бульоне [2]. Это эмпирическое наблюдение легло в основу экспериментальной гидродинамики и регулярно используется для подтверждения или опровержения новых теоретических концепций, позволяющих количественно и качественно оценить основные характеристики течений.

В работах Бенара [3-5] отражено продолжение экспериментальных исследований конвективных течений и описана структура течения, представляющая собой пчелиные соты. При интерпретации опыта Бенаром в качестве механизма образования конвекции при подогревании плоского горизонтального слоя указывалась связь вязкости жидкости с термокапиллярными силами, возникающими на свободной поверхности [3-5], а не подъемная сила, как предполагали другие исследватели [1, 2]. Несмотря на то, что изучение конвекции жидкости началось с плоских течений, не меньший интерес вызывают осесимметричные течения. Осесимметричная конвекция Бенара— Марангони очень важна при изучении ряда технологических процессов, например, при описании выращивания кристаллов методом Чохральского [6], в

Сергей Николаевич Аристов (д.ф.-м.н.), главный научный сотрудник, лаб. гидродинамической устойчивости. Просвиряков Евгений Юрьевич (к.ф.-м.н.), докторант, доцент, каф. основ конструирования.

астрофизике [7], геофизической гидродинамике [7-10], метеорологии и климатологии [7, 9, 10].

Точное решение уравнений, описывающих термокапиллярную конвекцию, впервые было найдено Бирихом [11]. Развитие и обобщение решения Бириха, интерпретацию физического смысла механизмов термокапиллярной конвекции можно проследить по работам [11-16]. Отметим, что до сих пор исчерпывающе не решена задача об описании осесимметричных термокапиллярных течений вблизи экстремальных значений температуры, в которых производная, вычисленная по радиальной координате, обращается в нуль. Построенные решения [11-16] не позволяют изучить анонсируемый класс течений, поэтому целью настоящей работы является построение точного решения системы Обербека—Буссинеска термокапиллярной конвекции, описывающего экстремальное температурное поле при числе Грасгофа, стремящемся к нулю. Найденное решение при необременительных ограничениях может быть использовано в качестве начального приближения для поиска устойчивых и неустойчивых решений осесимметричной конвекции Бенара—Марангони.

1. Постановка задачи. Запишем систему уравнений в цилиндрической системе координат, описывающих термокапиллярную конвекцию при установившемся осесимметричном движении вязкой несжимаемой жидкости (система уравнений Обербека—Буссинеска [17, 18]):

+ ---= + V Д V--2

дг дг г дг V г2

V-д^ + V*^ + — = VК - ^),

дг дг г V ^ г2 / (1)

дК дК дР V- ^ + V* —* = - — + + V ДР*, дг дг дг

т,дТ дТ 1 д (гК) д^* V.--+ V*— = хДТ,---—— +--* = 0.

г I г\ I г * ГЛ /V 5 ГЛ I ГЛ

дг дг г дг дг

В системе уравнений (1) введены следующие обозначения: V-, V* — радиальная, азимутальная и вертикальная скорости соответственно; Р — отклонение давления от гидростатического, делённое на постоянную среднюю плотность жидкости р; Т — отклонение от средней температуры; V, % — коэффициенты кинематической вязкости и температуропроводности жидкости соответственно; Д = ЦТ (гЛ) + —оператор Лапласа, определённый в цилиндрической системе координат.

При достаточно малом нагреве решения системы (1) могут быть разложены в ряды Тейлора с различным числом удерживаемых членов для анализа

течений [19, 20]. Решение системы (1) будем искать в следующем классе [19]:

г2 г2

VI = ги, ^ = ги, V* = ад, Т = То + -Тш Р = Ро + -Р1Ь (2)

Класс решений (2) обобщает решения, полученные Бирихом [11], Карманом [17], Сидоровым [20], Линем [21]. Подставим класс (2) в (1), получим полиномиальные выражения, которые в силу линейной независимости степенного базиса формируют систему уравнений, распадающуюся на две под-

системы:

dw -Pii ^ 2 du 2 ^ d2u

2u + jr=°> ^Т1 = 9PTii, u2 + w — - v2 = -Pu + v—^,

dz dz dz dz2

0 , dv d2 v dTii д 2Тц

2uv + w— = v—,, 2uTn + w—- = X ô 2 ;

dz dz2 dz dz2

dw dPo d2w dTo / d2To

(3)

w

дг + + ^^ = х{та + . (4)

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

нице задан локализованный источник тепла:

t = 9 2

i - ©2]

(5)

Здесь 9 — характерная разность температур, R — максимальное радиальное значение. Для основной системы (3) восьмого порядка граничные условия имеют следующий вид:

dw

при z = -h : v = w = —— = 0;

dz

при z = 0: w = ^ = 0, Pii =0, Tii = -R2, ^ = -Tii, dz R2 dz n

где a = da/dT температурный коэффициент поверхностного натяжения а, П — коэффициент динамической вязкости [18]. Для вспомогательной системы третьего порядка граничные условия имеют следующий вид:

r2 9

при z = -h : To + Tii — = 0; при z = 0: To = ^, Po = A. (6)

Здесь A — атмосферное давление [13, 17, 18]. Граничное условие на твёрдой поверхности задается из условия, согласно которому имеет место нулевая температура на радиальной границе. Раскладывая в ряды Тейлора граничные условия, убеждаемся, что класс решений (2) удовлетворяет граничным условиям.

2. Решение редуцированной системы уравнений Обербека—Буссинеска.

Для решения краевой задачи (3)—(6) приведем уравнения (3), (4) и граничные условия (5)-(6) к безразмерному виду. В качестве базиса выберем следующие размерные комплексы: h — вертикальный (поперечный) характерный размер; R — радиальный характерный размер; 9 — температура; gfî9h4/(vR2) — вертикальная скорость; gfî9h — давление, делённое на постоянную плотность. Используя выбранный базис, получим краевую задачу, записанную в безразмерной форме для систем (3) и (4):

dw dPii ^ ( 2 du 2\ л i d2u

2u + —= 0, —- = Tii; Gr u2 + w— - v2 = -Pn + —, ; dz dz V dz ) dz2

( dv \ d2v ( ^ dTii \ 1 d2Ti

Gr(^2uv + w—J = -q¿¿ , Gr ^2uTii + w-^^^ =

il

dz2 ' V dz ) Pr dz2

Сг62ад дад дР° + 62 д2ад

— = дг2

дг дг

-1 : V = дад

при г = = ад = дг

0 : ад = дv

при г = дг 0,

,гр г дТо 1/2 , д2Т°ч

+ То, с™— = -(6Ти + ;

1

Рг

д 2Т°

ди

Р11 =0, Т11 = -1, — = ^Т

дг

Т 1

при г = -1 : Т° + ^ = 0; при г = 0 : Т° = ^, Р° = 1.

(9)

Здесь Сг = дв^Л5/^2Я2) — осесимметричное число Грасгофа, записанное для вертикального и радиального характерных масштабов; Рг = и Mg = = ст/(рЛ2дв) — число Прандтля и Марангони соответственно; 6 = Л/Я — безразмерный геометрический параметр. Отметим, что если осуществлять обез-размеривание с одним характерным масштабом, то есть Л = Я, то получим классическое число Грасгофа: Сг = дв#Л3/V2 [1, 2, 6-18, 20].

При числе Грасгофа, равном нулю (Сг = 0), краевая задача (7)-(9) имеет аналитическое решение:

Т11 = 15(Mg - г - 1, Р11 = 15^ - 4) У - г,

и г4 г3 3 1

и = 151 Mg -7-77-7г + Mgz + - Mg - —,

3!

5

8

96'

1\ г5 г4 г2 /3 1

ад = 15 ( -2Mg + 2] 5!+ 2- - 2Mgу + (-3Mg + 48

Т° = 1562( - - Mg) ^т- + 62^т- + ( — - — Mg + -Mg62 - — )г + -, ° и ъ) 3! 2! V 8 2 6 26 8 / 2

1 ч г4 ог2г3 (15 15 5 г2 62ч г2 Р° = 4562 - - ^ — + 362— + ( — -— Mg + - Mg62 - — —+

1

3

2

2

15 15

5

2

62

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4

4!

3! 8

2

2

2!

+ - 26%) г + 1.

(10)

3. Исследование полученных решений. Проанализируем знаки значений гидродинамических полей (10). Рассмотрим для начала вертикальную скорость Уг = ад. Найденное решение принадлежит классу многочленов пятой степени. Отметим, что в силу граничных условий (8), (9) многочлен ад можно факторизовать следующим образом [22]:

ад = г (г + 1)2 д = г (г + 1)2

16 4) г + V 2 2^ ^ +48 4

Поскольку координата г определена на отрезке [—1; 0], то + 1)2 ^ 0, следовательно, знак скорости будет определяться знаком квадратичного полинома д, коэффициенты которого непрерывно зависят от числа Марангони. При Mg £ [0; 1/36] скорость ад принимает только отрицательные значения на области определения координаты (рис. 1), при Mg £ (1/36; 1/12) — как положительные, так и отрицательные значения, а при Mg ^ 1/12 — положительные значения (рис. 1). Таким образом, в зависимости от числа Марангони на

свободной (верхней) и твёрдой (нижней) границах скорость может одновременно иметь разные знаки (рис. 1). Полученный эффект смены знака вертикальной скорости на свободной границе при изменении числа Марангони наблюдался в экспериментальных и теоретических исследованиях [1, 11, 13, 18].

Знак радиальной скорости Уг = ги, очевидно, зависит от знака полинома и, который факторизуется следующим образом:

и = (г + 1)8 = (г + 1)

(5Mg 5 \ 3 (5Mg 1 \ 2

V 8 32;г3 "I — ++96^г2+

/5Ме 1 Л 3Ми 1

+ ( — + ^ +

8 96 8 96

На рис. 2 изображено множество значений функции и, что в силу класса (2) и области определения переменной г € [0; 1] с точностью до положительного множителя определяет множество значений радиальной скорости. При Mg € [0; 1/36) и (1/12;+го) полином в имеет один корень в области определения, а при Mg € [1/36; 1/12] —два корня. Данная локализация корней поясняет геометрию множества значений радиальной скорости.

Используя решение (10), приведем (рис. 3) распределение температуры на твёрдой границе (г = "1), зависящей, как нетрудно видеть, только от числа Марангони. Отметим, что на твёрдой границе при М^ € [0; 11/60] жидкость имеет отрицательную температуру, при Mg > 0 — положительную. Если рассматривать г € ("1;0), то знак и распределение температуры существенно зависят от величины 5, причем существуют такие значения числа Марангони, что параметр 5 принимает значения, большие шести. Это обстоятельство подчеркивает важность учета и введения различных масштабных параметров, а не одного [19].

Отметим, что для обобщённого давления можно провести анализ, аналогичный температурному, и уже на твёрдой границе (г = "1) проявляется зависимость распределения знаков приведенного давления от числа Маран-гони Mg и параметра 5.

На рис. 4-8 представлены изолинии обобщённой (модифицированной) функции тока [17], давления, отнесенного к плотности, и температуры при числе Марангони Mg = 0 и 5 = 1/2.

Для уточнения типа изолиний, изображенных на рис. 4-8, приведем алгоритм, позволяющий определить характер изолинии. Не ограничивая общности, рассмотрим алгоритм для температурного поля. Известно, что тип функции нескольких переменных (геометрически — поверхность) можно определить, анализируя знак гауссовой кривизны, который определяется знаками собственных значений ее матрицы Гессе [23-26]. Матрица Гессе, записанная для температуры, имеет следующий вид:

(дт / Ти гдГи.\

н(Т) = & дЩ = Д. ГЛ .

\dzdr дх2 ) V дх дх2 У

Изолинии температуры будут принадлежать к параболическому типу тогда, когда матрица Гессе имеет собственные значения одного знака, гиперболического— собственные числа разных знаков [23-26]. Воспользовавшись методами локализации собственных значений [22, 27, 28], получим, что изолинии температуры и других гидродинамических полей могут принадлежать только к гиперболическому типу.

г

-0,2 -0,4 -0,6 -0,8 - 1

0 0,05 0,1 0,15 0,2 0,25 Мё Рис. 1. Множество значений вертикальной скорости ш при различных числах Марангони Mg (закрашенная область соответствует положительным значениям, незакрашенная — отрицательным)

-0,2

-0,4 -0,6 -0,8

~ \) 0,05 0,1 0,15 Мё Рис. 2. Множество значений функции и при различных числах Марангони Mg (закрашенная область соответствует положительным значениям, незакрашенная — отрицательным)

0 0,2 0,4 0,6 0,8 Мё Рис. 3. Множество значений температуры на твёрдой границе при различных числах Марангони Mg (закрашенная область соответствует положительным значениям, незакрашенная — отрицательным)

-0,5 0 Рис. 4. Изолинии модифицированной функции тока при Mg = 10, 5 = 1/2 (изолинии построены для значений, принадлежащих отрезку [0; 1,4], с шагом, равным 0,1)

Р

-0,2

-0,4 -0,6 -0,8 - 1

-1 -0,5 0 0,5 г Рис. 5. Изолинии обобщённого давления при Mg = 0,01, 5 = 1/2 (изолинии построены для значений, принадлежащих отрезку [0,25; 0,9], с шагом, равным 0,05; для значений, принадлежащих отрезку [0,9; 1] с шагом, равным 0,01)

Рис. 6. Изолинии обобщённого давления при Mg = 10, 5 = 1/2 (изолинии построены для значений, принадлежащих отрезку [-30; 40], с шагом, равным 5)

-0,6

-0,8

-0,2

-0,4

- 1

Р

— 1

-1 -0,5 0 0,5

-1 - 0,5 0 0,5 2

Рис. 7. Изолинии температуры при Mg = = 0,01, 5 = 1/2 (изолинии построены для значений, принадлежащих отрезку

Рис. 8. Изолинии температуры при Mg = 10, 5 = 1/2 (изолинии построены для значений, принадлежащих отрезку [0; 60],

[—1; 0,4], с шагом, равным 0,05)

с шагом, равным 5)

На рис. 5, 6 и рис. 7, 8 показана зависимость направления теплового потока от величины числа Марангони. Пограничное значение числа Марангони вычисляется методами локализации корней степенных решений (10), которые были использованы выше [22, 27, 28].

Заключение. Построено новое точное решение системы уравнений Обер-бека—Буссинеска, описывающее осесимметричное термокапиллярное течение вязкой несжимаемой жидкости при локальном параболическом нагреве свободной поверхности. Рассмотрен случай, когда число Грасгофа стремится к нулю. Определены области постоянств знаков у множеств значений решений в зависимости от числа Марангони и безразмерного параметра 5. Найдено граничное значение числа Марангони, при котором происходит смена знака вертикальной и радиальной скоростей. Показано, что изолинии гидродинамических и температурного полей могут принадлежать только к гиперболическому типу.

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

Работа выполнена при поддержке РФФИ (проект № 12-01-00023-a) и ФСР МФП НТС (программа СТАРТ).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. А. В. Гетлинг, "Формирование пространственных структур конвекции Рэлея-Бенара"// УФН, 1991. Т.161, №9. С. 1-80; англ. пер.: A. V. Getting, "Formation of spatial structures in Rayleigh-Bénard convection" // Sov. Phys. Usp., 1991. Vol.34, no. 9. Pp. 737-776.

2. Ф. А. Гарифуллин, "Возникновение конвекции в горизонтальных слоях жидкости" // Соросовск. образоват. ж., 2000. Т. 6, №8. С. 108-114. [F. A. Garifullin, "Free convection in horizontal liquid layers" // Sorosovsk. Obrazovat. Zh., 2000. Vol.6, no. 8. Pp. 108-114].

3. H. Bénard, Les tourbillons cellulaires dans une nappe liquide propageant de la chaleur par convection, en regime permanent. These. Paris: Gauthier-Villars, 1901. 88 pp.; Ann. de Chim. et Phys., 1901. Vol. 23. Pp. 62-144.

4. H. Bénard, "Etude experimentale des courants de convection dans une nappe liquide. —

Régime permanent: tourbillons cellulaires"// J. Phys. Theor. Appl., 1900. Vol. 9, no. 1. Pp. 513-524.

5. H. Bénard, "Les tourbillons cellulaires dans une nappe liquide. - Methodes optiques d'observation et d'enregistrement" // J. Phys. Theor. Appl., 1901. Vol. 10, no. 1. Pp. 254-266.

6. Н. В. Никитин, С. А. Никитин, В. И. Полежаев, "Конвективные неустойчивости в гидродинамической модели роста кристаллов методом Чохральского" // Усп. мех., 2003. Т. 2, №4. С. 3-45. [N. V. Nikitin, S. A. Nikitin, V. I. Polezhaev, "Convective Instabilities in the Hydrodynamic Model of the Czochralski Crystal Growth" // Usp. Mekh., 2003. Vol.2, no. 4. Pp. 3-45].

7. Г. С. Голицин, Природные процессы и явления: волны, планеты, конвекция, климат, статистика. М.: Физматлит, 2004. 344 с. [G. S. Golitsin, Natural processes and phenomena: waves, planets, convection, climate, statistics. Moscow: Fizmatlit, 2004. 344 pp.]

8. В. В. Алексеев, А. М. Гусев, "Свободная конвекция в геофизических процессах" // УФН, 1983. Т. 141, №2. С. 311-343; англ. пер.: V. V. Alekseev, A. M. Gusev, "Free convection in geophysical processes" // Sov. Phys. Usp., 1983. Vol. 26, no. 2. Pp. 906-922.

9. А. С. Монин, Теоретические основы геофизической гидродинамики. Л.: Гидрометеоиз-дат, 1988. 424 с. [A. S. Monin, Fundamentals of Geophysical Fluid Dynamics. Leningrad: Gidrometeoizdat, 1988. 424 pp.]

10. А. Гилл, Динамика атмосферы и океана. Т. 1. М.: Мир, 1986. 397 с. [A. Gill, Atmosphere-ocean Dynamics. Vol. 1. Moscow: Mir, 1986. 397 pp.]

11. Р. В. Бирих, "О термокапиллярной конвекции в горизонтальном слое жидкости"// ПМТФ, 1966. №3. С. 69-72; англ. пер.: R. V. Birikh, "Thermocapillary convection in a horizontal layer of liquid" // J. Appl. Mech. Tech. Phys., 1966. Vol. 7, no. 3. Pp. 43-44.

12. L. G. Napolitano, "Plane Marangoni-Poiseuille flow of two immissible fluids" // Acta Astronaut., 1980. Vol.7, no. 4-5. Pp. 461-478.

13. В. К. Андреев, Решения Бириха уравнений конвекции и некоторые его обобщения: Препринт № 1-10. Красноярск: ИВМ СО РАН, 2010. 68 с. [V. K. Andreev, Birikh solution of convection equations and some of its generalizations: Preprint No. 1-10. Krasnoyarsk: Inst. Computer Modeling, Sib. Branch, Russian Acad. of Sci., 2010. 68 pp.]

14. O. Goncharova, O. Kabov, "Gas flow and thermocapillary effects of fluid flow dynamics in a horizontal layer" // Microgravity Sci. Technol., 2009. Vol. 21, no. 1. Pp. 129-137.

15. Л. Х. Ингель, М. В. Калашник, "Нетривиальные особенности гидротермодинамики морской воды и других стратифицированных растворов"// УФН, 2012. Т. 182, №4. С. 379-406; англ. пер.: L. Kh. Ingel, M. V. Kalashnik, "Nontrivial features in the hydrodynamics of seawater and other stratified solutions" // Phys. Usp., 2012. Vol. 55, no. 4. Pp. 356-381.

16. В. С. Бескин, "Осесимметричные стационарные течения в астрофизике" // УФН, 2003. Т. 173, №11. С. 1247-1253; англ. пер.: V. S. Beskin, "Axisymmetric steady flows in astrophysics" // Phys. Usp., 2003. Vol.46, no. 11. Pp. 1209-1214.

17. Л. Д. Ландау, Е. М. Лифшиц, Теоретическая физика. Т. VI: Гидродинамика. М.: Наука, 2006. 736 с. [L. D. Landau, E. M. Lifshits, Theoretical physics. Vol. VI: Fluid dynamics. Moscow: Nauka, 2006. 736 pp.]

18. Г. З. Гершуни, Е. М. Жуховицкий, Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972. 392 с. [G. Z. Gershuni, E. M. Zhukhovitskii, Convective stability of incompressible fluids. Moscow: Nauka, 1972. 392 pp.]

19. С. Н. Аристов, Вихревые течения в тонких слоях жидкости: Автореф. .. . дис. докт. физ.-мат. наук.. Владивосток, 1990. 32 с. [S. N. Aristov, Eddy currents in thin liquid layers: Dr. Phys. & Math. Sci. Thesis. Vladivostok, 1990. 32 pp.]

20. А. Ф. Сидоров, "Об одном классе решений уравнений газовой динамики и естественной конвекции" / В сб.: Численные и аналитические методы решения .задач механики сплошной среды ; ред. А. Ф. Сидоров, Ю. Н. Кондюрин. Свердловск: УНЦ АН СССР, 1981. С. 101-117. [A. F. Sidorov, "On a class of solutions of the equations of gas dynamics and natural convection" / In: Numerical and analytical methods for solving problems of

continuum mechanics; ed. A. F. Sidorov, Yu. N. Kondyurin. Sverdlovsk: Akad. Nauk SSSR, Ural. Nauchn. Tsentr, 1981. Pp. 101-117].

21. C. C. Lin, "Note on a class of exact solutions in magneto-hydrodynamics" // Arch. Rational Mech. Anal., 1958. Vol. 1, no. 1. Pp. 391-395.

22. Е. Е. Тыртышников, Матричный анализ и линейная алгебра. М.: Физматлит, 2007. 480 с. [E. E. Tyrtyshnikov, Matrix analysis and linear algebra. Moscow: Fizmatlit, 2007. 480 pp.]

23. С. П. Новиков, А. Т. Фоменко, Элементы дифференциальной геометрии и топологии. М.: Наука, 1987. 432 с. [S. P. Novikov, A. T. Fomenko, Elements of differential geometry and topology. Moscow: Nauka, 1987. 432 pp.]

24. В. И. Арнольд, Теория катастроф. М.: Наука, 1990. 128 с. [V. I. Arnol'd, Catastrophe theory. Moscow: Nauka, 1990. 128 pp.]

25. Т. Постон, И. Стюарт, Теория катастроф и её приложения. М.: Мир, 1980. 608 с. [T. Poston, I. Stewart, Catastrophe theory and its applications. Moscow: Mir, 1980. 608 pp.]

26. Р. Гилмор, Прикладная теория катастроф. Т. 1. М.: Мир, 1984. 350 с. [R. Gilmore, Catastrophe theory for scientists and engineers. Vol. 1. Moscow: Mir, 1984. 350 pp.]

27. Р. Хорн, Ч. Джонсон, Матричный анализ. М.: Мир, 1989. 655 с. [R. Khorn, Ch. Dzhonson, Matrix analysis. Moscow: Mir, 1989. 655 pp.]

28. М. М. Постников, Устойчивые многочлены. М.: Наука, 1981. 176 с. [M. M. Postnikov, Stable polynomials. Moscow: Nauka, 1981. 176 pp.]

Поступила в редакцию 22/1/2013; в окончательном варианте — 21/111/2013.

MSC: 76F02, 76M45; 76F45, 76R05, 76U05

ON ONE CLASS OF ANALYTIC SOLUTIONS OF THE STATIONARY AXISYMMETRIC CONVECTION BENARD-MARAGONI VISCOUS INCOMPREEIBLE FLUID

S. N. Aristov, E. Yu. Prosviryakov

1 Institute of Continuous Media Mechanics UB RAS, 1, Academician Koroleva st., Perm, 614013, Russia.

2 Kazan National Research Technical University named after A. N. Tupolev,

10, Karl Marx st., Kazan, 420111, Russia.

E-mails: asn@icmm.ru, evgen_pros@mail.ru

The purpose of this work is to find solutions for the system of equations Oberbeck-Boussinesq flat convection Benard-Marangoni a viscous incompressible fluid. In this viscous incompressible fluid the radial component of the temperature gradient may become zero. It is shown that the initial system may be reduced to the system of equations of ordinary differential equations of the eleventh order. We obtain the exact solution at the point of the extremum of the temperature (at zero including Grasgof's). Integration of equations is carried out in dimensionless variables, which are non-classical way: put the scale factor for each variable, and not by linear characteristic size of the layer. The solution is the initial approximation to the solution of convection Benard-Marangoni in numbers Grasgof's, the big zero.

Key words: axisymmetric thermocapillary convection (convection Benard-Marangoni), localized parabolic heaters, exact solution, isolines, Hessian matrix, eigenvalues, localization of polynomials roots, localization of eigenvalues of the matrix.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Original article submitted 22/I/2013; revision submitted 21/III/2013.

Sergey N. Aristov (Dr. Phys. & Math. Sci.), Chief Researcher, Lab. Hydrodynamical Instability. Evgeii Yu. Prosviryakov (Ph.D. Phys. & Math.), Doctoral Candidate, Dept. of Bases of Construction.

i Надоели баннеры? Вы всегда можете отключить рекламу.