CC BY
69
2013 Механика № 1
УДК 533.15:536.25
В.Н. Косов1, О.В. Федоренко2, Ю.И. Жаврин2, А.Т. Нысанбаева2, М.К. Асембаева2
Казахский национальный педагогический университет им. Абая, Алматы, Казахстан Институт экспериментальной и теоретической физики при Казахском национальном университете им. аль-Фараби, Алматы, Казахстан
ОБ УСТОЙЧИВОСТИ КОНВЕКТИВНЫХ ТЕЧЕНИЙ ПРИ ИЗОТЕРМИЧЕСКОМ ДИФФУЗИОННОМ СМЕШЕНИИ ТРЕХКОМПОНЕНТНЫХ ГАЗОВЫХ СМЕСЕЙ
В рамках линейной теории устойчивости проведен анализ экспериментальных данных по изучению неустойчивости механического равновесия при диффузионном смешении трехкомпонентных газовых смесей в изотермических условиях. Рассматривается канал в виде вертикального цилиндра конечной высоты при граничных условиях, предполагающих отсутствие переноса компонентов через стенки диффузионного канала. Результаты теории сравниваются с полученным в опытах положением границ устойчивости для системы 0,4722 Не + 0,5278 Аг - N2 при варьировании диаметра канала. Показано, что в условиях развитой конвекции существование максимумов интенсивности переноса компонентов в зависимости от термодинамических параметров связано с модой возмущений, определяющей различные виды конвективных течений. Структурные образования, движущиеся навстречу друг другу, число которых в поперечном сечении диффузионного канала определяется модой возмущения, взаимодействуют и тем самым определяют волнообразное изменение интенсивности неустойчивого процесса. Для рассматриваемой системы первый максимум интенсивности неустойчивого процесса в зависимости от диаметра, возникающий при б = 6 мм, характеризуется модой возмущений п = 3 и критическим числом Рэлея К3 = 972,7. При этих условиях происходит стабилизация конвективного переноса, характеризующаяся определенным типом движения конвективных формирований, т.е. движением по шести каналам (три тока и три противотока). Также результаты показали, что увеличение моды возмущений и диаметра диффузионного канала приводит к нелинейному увеличению критического числа Рэлея для всех компонентов, участвующих в переносе. Приведенные результаты хорошо согласуются с экспериментальными данными для системы 0,4722 Не + 0,5278 Аг - N2.
Ключевые слова: диффузия, конвекция, мода возмущений, число Рэлея, линейная теория устойчивости.
V.N. Kossov1, O.V. Fedorenko2, Yu.I. Zhavrin2,
A.T. Nysanbaeva2, M.K. Asembaeva2
Abay Kazakh National Pedagogical University, Almaty, Kazakhstan Institute of Experimental and Theoretical Physics at al-Farabi Kazakh National University, Almaty, Kazakhstan
CONVECTIVE FLOWS STABILITY AT THE ISOTHERMAL DIFFUSION MIXING OF TERNARY GAS MIXTURES
Analysis of the experimental data for the study of mechanical equilibrium instability by the diffusion mixing of ternary gas mixtures in isothermal conditions is carried out within the linear theory of stability. A channel is considered in the form a vertical cylindrical channel of a final height at the boundary conditions implying the absence of components’ transfer through the diffusion channel walls. Predictions of the theory are compared with the experimentally determined location of stability borders for the mixture 0.4722 He + 0.5278 Ar - N2 by a variation for different diameter of the channel. The existence of the intensity maximum of components’ transfer in dependence on the thermodynamic parameters under conditions of the developed convection connects with the mode of disturbances determining the various types of convective flows. Structural formations moving opposite to each other which numbers are defined by the mode of disturbances in the cross-section of the diffusion channel interacts and thereby defines the wave-like change of unstable process intensity. For the examined system first maximum of the unstable process intensity subject to the diameter appearing at d = 6 mm is characterized by the mode of disturbances n = 3 and the critical Rayleigh number R3 = 972.7. Stabilization of the convective transfer defining by the certain motion type of convective forming i.e. the motion along six channels (three flows and three counterflows) occurs under such conditions. Also, results reveal that increase of the mode of disturbances and the diameter of diffusion channel result in the non-linear increase of the critical Rayleigh number for all components taking part in the transfer. Further, the present results are in good agreement with the experimental data for the system 0.4722 He + 0.5278 Ar - N2.
Keywords: diffusion, convection, mode of disturbance, Rayleigh number, linear theory of stability.
Введение
Экспериментальное изучение диффузионного массопереноса в изотермических трехкомпонентных газовых смесях показало, что за счет различия в коэффициентах диффузии компонентов, а также при определенных значениях давления, температуры, исходных составов смеси, диаметра диффузионного канала и т.д. возникают условия для возникновения концентрационной конвекции, обусловленной действием сил Архимеда [1-4]. Визуальные наблюдения за конвекцией с помощью теневого метода, метода катарометрических датчиков позволили выявить многообразие конвективных течений [1, 5], которые условно разделяют на следующие группы: 1) хаотические течения с множеством противотоков; 2) «капельный» режим конвекции с отчетливо
регистрируемыми конвективными формированиями; 3) «неподвижные страты», определяющие ламинарную конвекцию. Определение закономерности их возникновения является актуальной задачей. В данной работе в рамках линейной теории устойчивости рассматриваются условия, приводящие к различным типам конвективных течений.
1. Математическая модель устойчивости
1.1. Основные уравнения
Макроскопическое движение изотермической тройной газовой смеси описывается общей системой уравнений гидродинамики, которая включает в себя уравнения Навье-Стокса, сохранения числа частиц
смеси и компонентов. При условии независимой диффузии, когда для
3 ^ 3
изотермической газовой смеси ^ = 0 ; ^ сг = 1, эта система уравне-
I-1 I-1
ний имеет следующий вид [6;7]:
Эи
■+(uVu)
2- Гл ^ - -
= -Vp + ^V u + — + £, Vdivu + pg, V3 J
Hn r/c ^
----v div (rn) = 0, —L + vVc. = -divj,., (1)
dt v ' dt 1
* / A _ sjc _ \
j1 =-(А^С1 + D12VC2 ) ,
* / sk ______ sk _ \
j2 =-(D21VC + D22VC2 ) ,
где D* - «практические» коэффициенты тройной диффузии.
Уравнения (1) дополняются уравнением состояния среды
р = р(с1,с2р), T = const,
которое связывает термодинамические параметры в (1).
При решении системы уравнений (1) применялся метод малых возмущений [6, 7]. Учитывая, что при L^r различия между возмущениями среднечисловой v и среднемассовой U скоростей в уравнении Навье-Стокса будут несущественны [7], окончательная система уравнений гравитационной концентрационной конвекции для возмущенных значений в безразмерных величинах примет вид [7]
(2)
= 0,
где Рг, = у/В* - диффузионное число Прандтля; Я* = gfiiAid4 /уВ* -парциальное число Рэлея; = В*./В22 - параметры, определяющие
соотношение между практическими коэффициентами диффузии.
Для решения системы уравнений (2) необходимо конкретизировать граничные условия. Для этого нами была рассмотрена задача
о неустойчивости диффузионного смешения в цилиндрическом канале конечной высоты.
В диффузионном канале конечной длины (рис. 1) существенны трехмерные движения, поэтому при аппроксимации скорости следует считать все компоненты вектора и отличными от нуля.
Рассматривая в цилиндрической системе координат периодические по ф движения и удовлетворяя условиям на твердых границах ^ = ±И (И = Ь/й - геометрический параметр, характеризующий устойчивость), можно записать аппроксимацию скорости в виде [6]
1.2. Решение задачи для цилиндрического канала конечной длины
2
Рис. 1. Ограниченный цилиндрический канал
и,,
= 4 (И2 - ) и (Г ) СОБ пф,
= 2(И2 - 22 )и(Г) СОБ пф,
= 2(И2 -22 )ю(Г)б1и пф, (п = 0,1,2...).
и
и
(3)
Радиальные функции и, и, ю должны обращаться в нуль на твердой боковой поверхности цилиндра (при г=1). Из уравнения непрерывности следует соотношение, которое связывает эти функции,
1 ё , ч п
----( г Ы +—ю- и = 0,
г ёг г
причем
и
= Ып ( кг) Ы (к)
Г
и
ю =
кЫп ( к )
п
[ Л ( кг )- Ып ( к ) ‘
п+1
(4)
Г
( кг )“ ( к )‘
.п+1
кЫп ( к )
где Ып (кг) - функция Бесселя п-го порядка, а параметр к находится из
уравнения
Ып(к) = (п +1) Ып(к).
(5)
В зависимости от значения п формулы (3)-(5) описывают движения разного вида. Значение п = 0 соответствует осесимметричным движениям; при п = 1 получается антисимметричное движение с границей раздела по вертикальной плоскости, проходящей через ось цилиндра, и т.д. Каждому значению п соответствует спектр значений параметра к из уравнения (5). Эти различные значения отвечают движениям разной радиальной структуры: с ростом к увеличивается число радиальных узлов скорости. Численное решение уравнения (5) позволило определить наименьшие корни:
п = 0, к = 5,1356; п = 1, к = 2,8064; п = 2, к = 4,0863.
1
1
до
Полагая для первых двух уравнений (2), что —- = 0 , найдем кон-
Э/
центрации компонентов из уравнений
V2о. =-и К ,; = 1,2,
* 7. . ” 7 7
где К1 =
12
А
\
К2 =
^11 А ^21
V А2 1
(тп ^12^21) (^11 ^12 ^21)
Будем считать, что о1 = (г, 7)соб пф, тогда
дг г дг г
Ып (кг) _Гп
Ып (к)
(И2 - 72) К
(6)
Полагая, что на торцевых поверхностях исчезают возмущения концентраций, из (6) вытекают дополнительные условия - обращение в
д2 /! /
нуль второй производной •)у 2 на торцах диффузионного канала.
/ 07
Таким образом, имеем следующие условия:
д2
при 7 = ±И = 0.
07
(7)
Эти условия позволяют выбрать следующую аппроксимацию:
(Г, 7) = (И2 - 7 2)(5И2 - 7 2)С; (Г), (8)
где С1 (г) - радиальная функция концентрации.
Для определения С;(г) применим метод Канторовича. Подставляя (8) в (6), умножая на зависящую от 7 часть функции /;(г, 7) и интегрируя в пределах от -И до И, получим уравнения
1
С; +-С;
Г
2
п 2 Т +а
VГ У
С =
11
248
Ып(кг)
Л (к)
Г
(9)
2 153 . т _
где а =тт7у, * = 1,2.
62И
При нахождении концентраций компонентов уравнения (9) реша-
5С- А
= 0, тогда конечное в центре реше-
лись с граничным условием
ние имеет вид [8]
дг
Г=1
11К (И2 -725И2 -72)СОБпф 248а2 (к2 +а2)
х
а
+
( кг) Ы ( к)
1п (аГ) <
о1п (а)
(к2 +а2)
п
(а2 + к2 )-а2к
(10)
Полученные выражения для и и о позволяют вычислить критическое число Рэлея с помощью интегрального соотношения Галеркина [6]:
Я =
|и7 -Аи7 &У
|и2 ■ о ёУ
(11)
Формула (11) определяет критические числа Рэлея в зависимости от параметра И и моды возмущения п. Критические числа Рэлея для рассматриваемого канала имеют следующие значения: Я1 = 67, 95,
Я2 = 329,1, Я3 = 942,5, Я4 = 2102, Я5 = 4033... Приведенные значения критических чисел Рэлея соответствуют наиболее вероятным течениям газовой смеси, образующимся в канале, и соответствуют случаю кратного деления площади поперечного сечения.
Для определения границы монотонной устойчивости рассматриваемой задачи умножим скалярно третье уравнение системы (2) на вектор и и проинтегрируем по всему объему У диффузионного канала
Эй
при условиях, что Ур = 0,— = 0 :
Э/
|иУ2иёУ + Я1х111и7о1ёУ + Я21и7о2ёУ = 0 .
(12)
Это уравнение дает на плоскости (Я1, Я2) прямую ММ, разделяющую области затухающих (диффузия) и нарастающих (конвекция) монотонных возмущений.
о
Описанный выше подход позволяет для каждого компонента получить свое критическое число Рэлея, которое зависит от термодинамических величин и моды возмущений, что позволяет говорить о структурированных свойствах гидродинамических течений каждого компонента.
2. Результаты численного эксперимента
2.1. Анализ экспериментальных данных
Возникновение конвективных течений при многокомпонентной диффузии определяется геометрическими размерами диффузионного канала и конкретными значениями параметров состояния: давления, исходного состава газовой смеси и т.п. Каждый из перечисленных параметров при постоянстве остальных характеристик может иметь граничное значение, определяющее переход системы из устойчивого состояния в неустойчивое и обратно. Однако экспериментальное изучение неустойчивости механического равновесия при диффузионном смешении трехкомпонентных газовых смесей в изотермических условиях показало, что в условиях развитой конвекции возможно существование максимумов интенсивности в зависимости от давления [9]. Дальнейшее увеличение давления приводит к уменьшению интенсивности переноса компонента, т.е. устойчивости конвективных течений относительно предыдущей моды возмущения. При этом данный эффект наблюдается не только когда смесь газов состоит из компонентов, коэффициенты диффузии которых отличаются в несколько раз, но и когда компоненты имеют близкие по значению коэффициенты диффузии [10]. Аналогичные нелинейные зависимости наблюдаются при различных диаметрах диффузионного канала [11] и угла его ориентации относительно вертикали [12].
Рассмотрим влияние диаметра канала на возникновение неустойчивого режима при диффузии в тройных газовых смесях, а также на устойчивость конвективных течений.
Полученные данные для системы 0,4722 Не + 0,5278 Аг - N показали [11], что в области давлений Р = 0,1-18,0 МПа при комнатной температуре для данной системы критический диаметр канала = 2,65 мм.
При увеличении диаметра выше наблюдается возникновение конвективных течений. Результаты измерения перешедшего количества гелия и аргона в зависимости от диаметра диффузионного канала пока-
зали, что со значения, превышающего d = 6 мм, фиксируется волнообразное изменение концентрации этих газов и интенсивности неустойчивого процесса [11]. На (рис. 2) представлены зависимости продиффундиро-вавших концентраций компонентов от диаметра канала. Экстремумы плотности потока при d = 6 мм и d = 18 мм можно объяснить тем, что существует некоторый характерный размер dc «ячейки конвекции».
0,2 и
0,18 ■
0,16 ■
0,14£
§ 0,12 ■
1 ОД •
о
0,08 ■
0,06 ■
0,04 ■
0,02 -----1-----1-----1-----1-----1---1----1-------1-
2 4 6 8 10 12 14 16 18
й, 10"3 м
Рис. 2. Концентрации продиффундировавших при неустойчивом переносе гелия и аргона в системе 0,4722 Не + 0,5278 Аг - N для различных диаметров диффузионного канала; Р = 2,54 МПа, Т = 298,0 К; точки •,▲ - экспериментальные данные для Аг
и Не соответственно
При некоторых отношениях d/dc будут создаваться наиболее удобные структурированные течения, целое число раз размещающиеся в капилляре с диаметром d. Другими словами, появление максимумов (по концентрации) объясняется взаимодействием структурных образований, движущихся навстречу друг другу.
2.2. Анализ результатов расчета
Как было показано выше, каждому п соответствует свое масс-штабное возмущение, определяющее переход от одного типа течения к другому. Уменьшение масштаба возмущений, т.е. рост числа п, может происходить за счет увеличения одного из параметров при соблюдении условия неизменности других параметров системы. Как показали опытные данные (см. рис. 2), изменение перешедшей концентрации компонентов в зависимости от диаметра носит нелинейный характер. Если предположить, что областям, связанным с нелинейными измене-
ниями интенсивности процесса, соответствуют новые масштабы возмущений, то экспериментальные парциальные числа Рэлея Яг должны
находиться вблизи граничных линий, отвечающих разным модам возмущений п.
Для сравнения предлагаемого в п. 1 подхода с опытными данными, приведенными в [11], представим их в виде парциальных чисел Рэлея. Парциальные числа Рэлея в соответствии с (2) можно записать следующим образом:
_ §пг4А^1 дсх
Я] —----------------;-, Я2
руВ11 &
§пг4 Дт2 дс2 руД*2 дх
(13)
т.
т3. Если
где т1 - масса молекулы г-го сорта; Ат1 = т1 - т3, Ат2
известны условия проведения опыта (давление, температура, состав смесей в каждой из колб, размеры диффузионного канала), то по формулам (13) можно найти Я1 и Я2 и тем самым определить изображающую данный опыт точку на плоскости (Я1, Я2). Из опыта известно, какой режим (диффузия или конвекция) имеет место при заданных условиях.
На (рис. 3) приведены опытные числа Рэлея и критические линии для различных мод возмущений трехкомпонентной системы 0,4722 Не + + 0,5278 Аг - N2.
Рис. 3. Числа Рэлея при различных диаметрах диффузионного канала для системы 0,4722 Не + 0,5278 Аг - ^; точки - экспериментальные данные, соответствующие диаметру: 1 - 3,4; 2 - 4,0; 3 - 6,0; 4 - 6,4; 5 - 7,0; 6 - 7,6; 7 - 8,0; 8 - 9,6; 9 - 10,6;
10 - 12,0 мм
Сравнительный анализ (см. рис. 2, 3) показывает, что для смеси
0,4722 Не + 0,5278 Аг - N первая смена режима происходит при диаметре dif » 2,65 мм, т.е. наблюдается переход от диффузионного процесса смешения к конвективному, соответствующему антисимметричному критическому движению. Для моды п = 1 критическое число Рэлея составило Я1 = 73,05 . При диаметре d„ » 6 мм происходит
изменение масштаба возмущений уже к установившемуся конвективному режиму течения и возникает следующий конвективный режим, которому соответствует критическое число Рэлея Я3 = 972,7, вычисленное для п = 3, и т.д.
Таким образом, проведенные исследования показали, что различные виды конвективных течений, возникающих в цилиндрическом диффузионном канале при изотермической трехкомпонентной диффузии, связаны с модой возмущений и критическим числом Рэлея, зависящим от рассматриваемой системы и термодинамических параметров.
Заключение
Приведенная модель описания неустойчивого массопереноса в изотермических трехкомпонентных газовых смесях позволяет для различных мод возмущения определить спектр критических чисел Рэлея, характеризующих области устойчивой и неустойчивой диффузии, и определить параметры, при которых реализуются различные типы конвективных движений.
Работа выполнена при финансовой поддержке Фонда науки Министерства образования и науки Республики Казахстан (грант № 1134/ГФ).
Библиографический список
1. Влияние давления на устойчивость диффузии в некоторых трехкомпонентных газовых смесях / Ю.И. Жаврин [и др.] // ЖТФ. -1984. - Т. 54, № 5. - С. 943-947.
2. Жаврин Ю.И., Косов В.Н. Образование структур и концентрационная конвекция при изотермической диффузии трехкомпонентных газовых смесей через переменное число каналов равной площади // Письма в ЖТФ. - 1993. - Т. 19. - Вып. 10. - С. 18-21.
3. Жаврин Ю.И., Косов В.Н. Влияние температуры на процесс диффузионной неустойчивости // ИФЖ. - 1988. - Т. 55, № 1. - С. 92-97.
4. Влияние частоты вращения диффузионного аппарата на процесс смешения в трехкомпонентной газовой смеси / Ю.И. Жаврин [и др.] // Письма в ЖТФ. - 2003. - Т. 29. - Вып. 3. - С. 53-57.
5. Применение теневого метода для визуализации конвективных течений, образующихся при диффузии в многокомпонентных газовых смесях / Н.Д. Косов [и др.] // Теплофизика и аэромеханика. - 1994. -Т. 1, № 1. - С. 87-90.
6. Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. - М.: Наука, 1972. - 392 с.
7. Косов В.Н., Селезнев В.Д. Аномальное возникновение свободной гравитационной конвекции в изотермических тройных газовых смесях; УрО РАН. - Екатеринбург, 2004. - 151 с.
8. Аномальная гравитационная конвекция трехкомпонентной газовой смеси в вертикальном цилиндре конечной высоты / Ю.И. Жаврин [и др.] // Вестник ТГТУ. - Тамбов, 2005. - Т. 11, № 1А. - С. 94-102.
9. Жаврин Ю.И., Косов Н.Д., Белов С.М. Исследование неустойчивости при диффузии смеси гелия с аргоном в азот в области давлений 1,5-15 МПа // Молекулярный массоперенос и струйные течения -Алматы, 1984. - С. 3-7.
10. Experimental study of diffusion instability in three-component gas mixture without density gradient / Yu.I. Zhavrin [and others] // Technical Physics Letters. - 2011. - Vol. 37. - No. 8. - P. 721-723.
11. Косов B.H., Жаврин Ю.И. Экспериментальное исследование на диффузионную устойчивость некоторых изотермических трехкомпонентных газовых систем // Изв. АН КазССР. Сер. Физ.-мат. науки. -1990. - № 2. - С. 66-69.
12. Болотов И.В., Жаврин Ю.И., Косов В.Н. Исследование диффузионной неустойчивости в наклонном канале // Вопросы тепломассообмена - Алматы, 1989. - С. 7-11.
References
1. Zhavrin Yu.I., Kosov N.D., Belov S.M., Tarasov S.B. Vliyanie dav-leniya na ustoichivost diffuzii v nekotorykh trekhcomponentnykh gasovykh smesyakh [Effect of pressure on the diffusion stability in some three-component gas mixtures]. Zhurnal Tekhnicheskoy Fiziki, 1984, vol. 54, no. 5, pp. 943-947.
2. Zhavrin Yu.I., Kosov V.N. Obrazovanie structur i concentracion-naya konvectsiya pri isotermicheskoi diffusii trekhcomponentnykh gas-ovykh smesei cherez peremennoe chislo kanalov ravnoi ploschadi [Structures formation and concentration gravitation convection under isothermal diffusion of three-component gas mixtures through the variable number equal area channels]. Pis’ma v Zhurnal Tekhnicheskoy Fiziki, 1993, vol. 19, no. 10, pp. 18-21.
3. Zhavrin Yu.I., Kosov V.N. Vliyanie temperatury na process diffu-sionnoi neustoichivosti [Effect of temperature on the process of diffusion instability], Inzhenerno-Fizicheskii Zhurnal, 1988, vol. 55, no. 1, p. 92-97.
4. Zhavrin Yu.I., Kosov V.N., Kulzhanov D.U., Poyarkov I.V., Anku-sheva N.B. Vliyanie chastoty vrascheniya diffusionnogo apparata na process smesheniya v trekhcomponentnoi gasovoi smesi [Effect of rotation frequency of the diffusion apparatus on the mixing process in three-component gas mixture]. Pis’ma v Zhurnal Tekhnicheskoy Fiziki, 2003, vol. 29, no. 3, pp. 53-57.
5. Kosov N.D., Bychkov A.G., Zhavrin Yu.I., Kosov V.N. Primenenie tenevogo metoda dlya vizualizatsii konvektivnykh techenii, obrazuyuscikh-sya pri diffuzii v mnogokomponentnykh gasovykh smesyakh [Application of the schlieren method for the visualization of convective flows organizing under diffusion in multicomponent gas mixtures]. Teplofizika i aeromek-hanika, 1994, vol. 1, no. 1, pp. 87-90.
6. Gershuni G.Z., Zhukhovitskii E.M. Konvectivnaya ustoichivost’ neszhimaemoi zhidkosty [Convective stability of incompressible fluids]. Moscow: Nauka, 1972. 392 p.
7. Kosov V.N., Seleznev V.D. Anomal’noe vozniknovenie svobodnoi gravitatsionnoi konvektsii v isotermicheskikh troinykh gasovykh smesyakh [Anomalous occurrence of the free gravitational convection in isothermal ternary gaseous mixtures]. Yekaterinburg: Uralskoje otdelenie Rossijskoy akademii nauk, 2004. 151 p.
8. Zhavrin Yu.I., Seleznev V.D., Poyarkov I.V., Fedorenko O.V. Anomal’naya gravitatsionnaya konvetsiya trekhkomponentnoi gasovoi smesi v vertikal’nom tsilindre konechnoi vysoty [Abnormal gravitational convection of three-component gaseous mixture in vertical cylinder of finite height]. Vestnik tambovskogo gosudarstvennogo tekhnicheskogo univer-siteta, 2005, vol. 11, no. 1A, pp. 94-102.
9. Zhavrin Yu.I., Kosov N.D., Belov S.M. Issledovanie neustoichi-vosti pri diffuzii smesi geliya s argonom v oblasti davlenii 1,5-5 MPa [Study of instability by the diffusion of helium and argon mixture into the nitrogen in the pressure range 1,5-5 MPa]. Molekulyarnyi massoperenos i struinye techeniya, Almaty, 1984, pp. 3-7.
10. Zhavrin Yu.I., Moldabekova M.S., Poyarkov I.V., Mukameden-kyzy V. Experimental’noe issledovanie diffusionnoi neustoichivosti v trekhkomponentnoi gasovoi smesi pri nulevom gradient plotnosty [Experimental study of diffusion instability in three-component gas mixture without density gradient]. Pis’ma v Zhurnal Tekhnicheskoy Fiziki, 2011, vol. 37, no. 8, pp. 721-723.
11. Kosov V.N., Zhavrin Yu.I. Experimentalinoe issledovanie na dif-fusionnyu ustoichivost’ nekotorykh isotermicheskikh trkhkomponentnykh gasovykh system [Experimental study on the diffusional stability of some isothermal tree-component gas systems]. Izvestia Akademii nauk KazSSR. Fiziko-matematicheskie nauki, 1990, no. 2, pp. 66-69.
12. Bolotov I.V., Zhavrin Yu.I., Kosov V.N. Issledovanie diffuzionnoi neustoichivosti v naklonnom kanale [Study of diffusion instability in the inclined channel]. Voprosy teplomassoobmena, Almaty, 1989, pp. 7-11.
Об авторах
Косов Владимир Николаевич (Алматы, Казахстан) - член-корреспондент АН РК, доктор физико-математических наук, профессор, проректор по научной работе Казахского национального педагогического университета им. Абая (050100, г. Алматы, пр. Достык, 13, e-mail: kosov_vlad_nik@list.ru).
Федоренко Ольга Владимировна (Алматы, Казахстан) - кандидат физико-математических наук, доцент кафедры теплофизики и технической физики Казахского национального университета им. аль-Фараби (050038, г. Алматы, пр. аль-Фараби, 71, e-mail: fedor23.04@mail.ru).
Жаврин Юрий Иванович (Алматы, Казахстан) - доктор физико-математических наук, профессор, профессор кафедры теплофизики и технической физики Казахского национального университета им. аль-Фараби (050038, г. Алматы, пр. аль-Фараби, 71, e-mail: zhavrin@physics.kz).
Нысанбаева Асель Талиповна (Алматы, Казахстан) - магистрант кафедры теплофизики и технической физики Казахского национального университета им. аль-Фараби (050038, г. Алматы, пр. аль-Фараби, 71, e-mail: a.nysanbaeva@gmail.com).
Асембаева Мансия Кабыловна (Алматы, Казахстан) - кандидат физико-математических наук, доцент кафедры теплофизики и технической физики Казахского национального университета им. аль-Фараби (050038, г. Алматы, пр. аль-Фараби, 71, e-mail: zhavrin@bk.ru).
About the authors
Kosov Vladimir Nikolaevich (Almaty, Kazakhstan) - corresponding member of the Academy of Sciences of RK, Doctor of Physical and Mathematical Sciences, Professor, Vice-rector on the scientists research, Kazakh National Pedagogical University (13, Dostyk, 050100, Almaty, Kazakhstan, e-mail: kosov_vlad_nik@list.ru).
Fedorenko Olga Vladimirovna (Almaty, Kazakhstan) - Ph.D. in Physical and Mathematical Sciences, Ass. Professor, Department of Thermophysics and Technical Physics, al-Farabi Kazakh State University (71, al-Farabi, 050038, Almaty, Kazakhstan, e-mail: fedor23.04@mail.ru).
Zhavrin Yuri Ivanovich (Almaty, Kazakhstan) - Doctor of Physical and Mathematical Sciences, Professor, Department of Thermophysics and Technical Physics, al-Farabi Kazakh State University (71, al-Farabi, 050038, Almaty, Kazakhstan, e-mail: zhavrin@physics.kz).
Nysanbaeva Asel Talipovna (Almaty, Kazakhstan) - undergraduate, Department of Thermophysics and Technical Physics, al-Farabi Kazakh State University (71, al-Farabi, 050038, Almaty, Kazakhstan, e-mail: a.nysanbaeva@gmail.com).
Asembaeva Mansiya Kabylovna (Almaty, Kazakhstan) - Ph.D. in Physical and Mathematical Sciences, Ass. Professor, Department of Thermophysics and Technical Physics, al-Farabi Kazakh State University (71, al-Farabi, 050038, Almaty, Kazakhstan, e-mail: zhavrin@bk.ru).
Получено 15.02.2013
