Математическое моделирование неустойчивости тонкого слоя вязкой жидкости

УДК 532.5

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НЕУСТОЙЧИВОСТИ ТОНКОГО СЛОЯ ВЯЗКОЙ ЖИДКОСТИ

Л.А. Прокудина, Е.А. Саламатов

MATHEMATICAL SIMULATION OF TREACLE SLICE UNSTEADINESS

LA. Prokudina, E.A. Salamatov

Осуществлено исследование неустойчивости тонкого слоя вязкой жидкости в линейном и нелинейном приближениях. Получены аналитические зависимости волновых характеристик от физико-химических факторов, оказывающих влияние на течение и его устойчивость. Представлены волновые характеристики и области неустойчивости для воды в диапазоне чисел Рейнольдса Re < 12. Исследована эволюция волновых возмущений в окрестности кривой нейтральной устойчивости.

Ключевые слова: тонкий слой, вязкая жидкость, неустойчивость, волновые характеристики.

The research of treacle slice unsteadiness in linear and nonlinear approximation was held. The analytic dependences of wave characteristics upon physical and chemical factors influencing upon its process and steadiness were obtained. The wave characteristics and unstable zones for water within the range of Reynolds number Re < 12 are presented. The evolution of wave disturbance near the neutral stability curve was explored.

Keywords: slice, treacle, unsteadiness, wave characteristics.

В основе многих технологических процессов различных отраслей промышленности лежат процессы течения тонких слоев вязких жидкостей (жидких пленок). Жидкие пленки находят применение в энергетической, металлургической, пищевой, атомной, химической и других отраслях промышленности.

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

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

Избежать нарушений технологического процесса и аварийных ситуаций можно, управляя течением жидкого слоя или совершенствуя методы

Прокудина Людмила Александровна - д.ф.-м.н., доцент, профессор кафедры прикладной математики ЮУр-ГУ; prima@prima.susu.ac.ru.

Саламатов Евгений Александрович - аспирант кафедры прикладной математики ЮУрГУ; зеа2009-1983@mail.ru.

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

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

В рамках математической модели неустойчивых режимов течения жидких пленок [1,2] иссле-

Prokudina Liudmila Alexandrovna - PhD, associate professor, professor of applied mathematics department of SUSU; prima@prima.susu.ac.ru.

Salamatov Evgeny Alexandrovich - post graduate student of applied mathematics department of SUSU; sea2009-1983@mail.ru.

д ЛЭш д4і|/ З3ш / \32ш 9ш Зш д2ці

а? _ 1 Ьг+аі +а4 тт+(ае+ дю )тт+ап ^+^ +«16^ —у+«і?

дх ) ді дх

дх'

дх*

дх

дх

дх

'д2Ч | дхді дх ді

+а

'22

Зш Э2ш Э3ш + <* & * ' ~

гду

2_+а28Ч'—у + ^ЗО ЙС

Эш 3 ш 3 ш

Т' + Ч'—г ,,. 3 д-4

дх дхъ дх4,

2 9ш 2 3 Ч/ 2 Э V

+ «37V^7 + «39У —у+«40У Г-т;

гас йг дхді

Эш Эш +ам^——+а^

Эу V дх у

Эш Э2ш Э3ш 2 34ш 2 ^Зш л

+ а^-^ + а^+0^ іГ *

О)

Ле2^ ЛеМ 3 „ Зг,, 5 „2г, ЛеМ „ „ „

«і =——, «4 =--2—’ °6 =-----2 + 40 р(х + р)’аі =й р’аю =----^~,аи = ~Кер~КеТ>

3 9

а14 = -2/?е/г - Лет, а16 = -ЯеМ + - Яе3Рх +—Ле3/*’2, а17 = 4а7, а22 = аХ6, а26 = а4, а28 = За4, а30 = 3^,

$ 20

/?<?м Зі 9^о Зі9'і/>

а37 = -ЯеР, а39 =- —н—/?е Fт+—Ле F , а40 = 6а7, а44 =120;, а45 =-ЛеМ +—Ле ^тч-—Ле ^ ,

«49 = 2«4> «л = 6а1; а55 = Заи а58 = За4

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

1. Постановка задачи

В прямоугольной системе координат ОХ рассмотрим гравитационно стекающий по вертикальной плоской поверхности с температурой Т„ тонкий слой вязкой жидкости.

Плоскость ОХг связана с поверхностью, по которой движется жидкая пленка; ось ОХ направлена по течению пленки, ось ОУ - по нормали в пленку жидкости.

Движение пленки сопровождается обтеканием ее свободной поверхности газовым потоком, который движется параллельно плоскости 0X2 и создает на свободной поверхности пленки постоянное касательное напряжение.

Учтем поверхностную вязкость, обусловленную наличием на свободной поверхности пленки нерастворимых поверхностно-активных веществ (НПАВ).

Безразмерная толщина пленки: 5 = б(дг, /) = = 1 + н/(х, /), где у(х, /) - отклонение свободной поверхности пленки от стационарного состояния 50 (рис. 1).

В таком случае эволюция свободной поверхности жидкой пленки описывается уравнением (1) [2]. Коэффициенты уравнения (1) учитывают физико-химические факторы, качественно влияющие на течение жидкой пленки и ее устойчивость: чис-

ло Рейнольдса Ле, число Фруда коэффициент поверхностной вязкости И, касательное напряжение на поверхности раздела газ-жидкость т, поверхностное натяжение а, параметр Марангони М.

В рамках уравнения (1) проведено исследование неустойчивости жидкой пленки в линейном и нелинейном приближениях для чисел Рейнольдса Яе < 12.

2. Неустойчивость жидкой пленки в линейном приближении

Пусть на свободной поверхности жидкой пленки развивается монохроматическое возмущение Ц1(х, 0 = Л-ехр(/(Ьг - (сог + /0»,)О), где к - волновое число, шг - частота, ш, - инкремент.

Подставив монохроматическое возмущение \)/

Дисперсионное уравнение (2) связывает волновые характеристики жидкой пленки с физикохимическими факторами ст, х, М и позволяет в явной форме выразить их, в частности, частоту сог:

{^-ахкА + а6кг^а1к+[аАкг-апк^

(а7Л) +1

инкремент О/ (рис. 2):

ахк,А - а6к2 + аук^к3 - апк)

<о, =-

(оуЛ) +1 найти фазовую скорость С (рис. 3):

С = аг/к,

групповую скорость Сгр (рис. 4):

Ор = Лог/й№.

Расчет волновых характеристик позволяет оценить эффективность протекающих на поверхности раздела газ-жидкость процессов [4].

Инкремент неустойчивого монохроматического возмущения с волновым числом к положителен, ш, > 0, (рис. 2), что позволяет в линейном приближении рассчитать область неустойчивости жидкой пленки (рис. 5) в плоскости параметров (Яе, к); также рассчитывается кривая максимального инкремента со,тах [5], соответствующая возмущениям, которые обладают наибольшей скоростью роста, с наибольшей вероятностью реализуются в эксперименте, оказывают наибольшее влияние на процесс течения жидкой пленки при работе пленочных аппаратов [6].

Рис. 1. Схема течения жидкой пленки:

1 - поверхность отекания;

2 - жидкая пленка;

3 - газовый поток

Рис. 2. График инкремента:

1 -/?е = 5;

2-Яе = 9;

3-Ле = 10

Рис. 3. График фазовых скоростей:

1 -Ле = 5;

2 - Ле = 9;

З-Ле = 10

у

0,05 0,1 к

Рис. 4. Групповые скорости

1 -Яе = 5;

2-Яе = 9;

3-Яе=10

Рис. 5. Область неустойчивости в линейном приближении:

1 - кривая нейтральной устойчивости;

2 - кривая максимального инкремента

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

Известно, что присутствие НПАВ на поверхности раздела газ-жидкость стабилизирует течение, способствует уменьшению скоростей роста возмущений, фазовых, групповых скоростей; газовый поток оказывает стабилизирующее влияние на течение жидкой пленки в режиме прямотока и дестабилизирующее - в режиме противотока [4]. В условиях дестабилизации наблюдается увеличение скоростей роста возмущений, фазовых, групповых скоростей. При этом можно найти такие значения физикохимических факторов, при которых течение станет устойчивым [4] или, наоборот, жидкая пленка разрушится [2].

Вычислительный эксперимент для воды в условиях свободного стекания (т = О, N = О, М = 0) показал, что течение жидкой пленки неустойчиво во всем диапазоне чисел Рейнольдса Яе < 12 (рис. 5).

Чем больше значение числа Рейнольдса Яе, тем зона неустойчивости шире по диапазону волновых чисел к, выше максимальная скорость роста возмущений (рис. 2, табл. 1), тем меньше фазовые (рис. 3) и групповые (рис. 4) скорости.

Фазовые скорости неустойчивых возмущений меньше фазовых скоростей устойчивых возмущений (рис. 3). Расчеты показали, что минимальное значение фазовых скоростей С достигается на кривой максимального инкремента.

Таблица 1

Расчет кривых нейтральной устойчивости и максимального инкремента

Re к Ю/тах

Ю/ = 0 ímax

1 0,02714 0,01919 0,00022

5 0,10378 0,07245 0,01536

9 0,16937 0,10974 0,05464

10 0,18491 0,11628 0,06419

Групповые скорости Сгр неустойчивых возмущений меньше, чем групповые скорости устойчивых возмущений (рис. 4).

3. Неустойчивость жидкой пленки в нелинейном приближении

Задача исследования неустойчивости течения тонкого слоя вязкой жидкости не исчерпывается линейным анализом. Учет нелинейностей позволяет уточнить результаты, полученные в линейной постановке, дает возможность выяснить влияние нелинейного взаимодействия возмущений на течение жидкой пленки.

Амплитуда возмущения, развивающегося на свободной поверхности жидкой пленки, представима в виде (7):

A(x,z,t) = % е'°^+ах е'(01«+М+

+а2 е'^'ЬМ+Яз е'

где а0 - амплитуда несущей волны; аъ а2 - амплитуды боковых волн; а3, а4 - наклонных волн.

= уоо + Р2 (2а0ахО2 sin Д + 2(^0$^ sin Д,) -

dt

-Pi(«o +2% ах + 2a0al +2а0а| +2a0aj +

+200^02 eos Д + 2о0а3а4 eos Aj);

у = (Мй-^2«1+у)аг

~P2 sin(-A)-2a2a3a4 sin(Aj - Д)|-

-$x[2c%ax + + 2axc% +2axa% +

+2alá^ + ^a2 eosД+ 2a2a3a4 cos

^ = -{kxdx+k2xax-y)a2--P2(a¿ajSÍn(-A)-2a1úr3a4sin(A1 - A)j--pj |2flQ«2 +2axa2 +a2 +2 a2a2 +

+2a2a4 + Qqüx cos

^- = (kzd2-k2a3+y)a3-

~P2 (aoa4 Sín + 2ala2a4 SÍ11 (^1 ~ ^)) -~Pj (2ajja3+2axa3 +2 a2a3 + a3 + 2a3a4 +

+aga4 cosAj +2a1o2íí4cos

^ = "(Мг + kza3 “y)«4 "

-P^aoí^sinA] +2^0303 sin (Aj -A)j--P, (2a¿a4 + 2axa4 + 2 a2a4 + 2a¡a4 + a4 +

+cfia3 eos A] +2aja2a3 cos

— = 4pj (axa2 sin A + a3a4 sin Aj) +

+4P2 (a{a2 eos Д+а3<з4 cosA,) +

+P2(-2 üq +af +a2^~ 2kxa2 +

+—(aoa2sin(—Д)—2a2a3a4sin(Ax - A))-

~—(ala2cosA+2a2a3a4cos(Ax -A)) + ax v >

+—^aoúfisin(-A)-2a1a3a4sin(A1 -Д))-

^aoa1cosA + 2a1o3a4cos

= 4pt (aj Аз sin A + a3a4 sin Aj) +

+4p2 (axa2 eos Д+a3a4 eos Aj) +

+P2 (-2ao + 2a\ + a3 + a4 j - 2k2a4 +

+—(aoa4sinAj + 2a1a2a4 sin^ -A))-

«3v '

~—(aoa4 eos A, + 2 axa2a4 cos(A-Aj)) +

+—{c¡la3 sin At +2a!a2a3 sin^ - A))-

-—

Из уравнения (1) выведена система дифференциальных уравнений для амплитуд a¡, i = 0,4 и сдвигов фаз волн Д, Ах (8)-{14), или в векторной форме:

f-FOO,

где X = [<?о «1 «2 Оз «4 А А[]г.

Коэффициенты системы (8)—(14) у, dx, d2, аь а2, а3, а4, Pi (рис. 6), р2 (рис. 7) зависят от физикохимических факторов (N, т, а, М).

Рассчитаем область неустойчивости жидкой пленки в рамках нелинейного приближения. Для этого исследуем неустойчивость системы (15) относительно амплитуд a¡, i = 0,4.

Перейдем к уравнениям первого приближения для системы (15):

^ = Об)

dt

где у - вектор вариаций неизвестных системы (15), J = |(^/9vl'/)[((=0|> '»7=0,6.

Запишем, х = [уг zT]r, где у = [а0 ах а2 а3 а4]г -вектор переменных, относительно которых исследуется неустойчивость жидкой пленки, вектор «контролируемых» переменных; z = [А AJr - вектор «неконтролируемых» переменных.

Тогда (16) представим в виде ¡ys=Ays+Bzs;

(А =су*+02„

где А, В, С, D - соответствующие подматрицы матрицы J; y« zs - векторы вариаций «контролируемых» и «неконтролируемых» переменных соответственно.

Для системы (17) можно показать, что:

1. В»0.

2. Для неравенства »

где R(ys, Zj) - остаточный член разложения вектор-функции F(x) в ряд Тейлора в окрестности точки

\|/ = 0, справедлива оценка % ~ |sup|a, |v=0| j .

В таком случае из обобщенной теоремы Ляпунова об устойчивости по первому приближению [8] следует, что для устойчивости возмущений, развивающихся на свободной поверхности тонкого слоя вязкой жидкости, относительно амплитуд

а, необходимо и достаточно, чтобы все собственные числа подматрицы А системы (17) имели отрицательные вещественные части. Если среди собственных чисел подматрицы А найдется хотя бы одно с положительной вещественной частью, то возмущения на свободной поверхности жидкой пленки будут неустойчивыми относительно амплитуд Я/.

Спектр подматрицы А исследуем с помощью критерия Зубова [9]: для того чтобы все собственные числа подматрицы А имели отрицательную вещественную часть, необходимо и достаточно

выполнения условия lim й" = 0, где п - целое

Л->СО

число, Q - функционально-преобразованная матрица, S2 = Е-2(Е-А)-1.

Для исследования неустойчивости возмущений на поверхности жидкой пленки в рамках критерия Зубова реализован алгоритм [11], в котором использованы достаточные формы критерия Зубова [10]. Чтобы вещественные части собственных чисел подматрицы А были отрицательными, достаточно, чтобы, начиная с некоторого выполнялось условие

|ос1<1.

С другой стороны, если, начиная с некоторого С,, справедливо неравенство

где п0 - порядок подматрицы А («0 = 5), то среди собственных чисел подматрицы А найдется хотя бы одно с положительной вещественной частью.

На каждом этапе алгоритма анализируется выполнение достаточных условий критерия Зубова

(18)—(19). Если выполняется какое-либо из условий (18>—(19), то возмущения в точке плоскости параметров (Яе, к) признаются, соответственно, устойчивыми или неустойчивыми.

Алгоритм реализован таким образом, чтобы осуществлялся анализ выполнения условий (18)—

(19) для последовательности матриц О, £22, О4, О16, О32 и т.д. Это дает возможность за сравнительно небольшое число шагов осуществить анализ высоких степеней матрицы О, поскольку на V этапе алгоритма анализируется 2у“|-я степень О.

Если за заранее заданное число этапов добиться выполнения условия (18) или (19) не удается, то в возмущения в точке плоскости параметров (Яе, к) считаются неустойчивыми. Такая ситуация возможна, если максимальное число этапов алгоритма задано неверно и необходимо увеличить число Уши, чтобы проанализировать выполнение условий (18>—(19) для еще более высоких степеней матрицы П. Если матрица Е-А оказывается вырожденной, т.е. осуществить расчет матрицы О невозможно, то возмущения также считаются неустойчивыми.

Алгоритм позволяет рассчитать кривую нейтральной устойчивости в плоскости параметров (Яе, к) в рамках системы дифференциальных уравнений (15).

Если в плоскости параметров (Яе, к) ввести сетку с шагом А по диапазону волновых чисел к, то алгоритм обеспечит расчет кривой нейтральной устойчивости с абсолютной погрешностью 8 < /г/2. Для расчета кривой нейтральной устойчивости с абсолютной погрешностью 80 в диапазоне волновых чисел к е [&„иП; Ащах] необходимо ввести сетку с ЧИСЛОМ узлов А'о = (¿шах ~ АтшУ(2е0).

Вычислительный эксперимент по исследованию неустойчивости жидкой пленки осуществлен для воды (т = 0, М= 0, N - 0), рассчитана область неустойчивости тонкого слоя вязкой жидкости (Яе < 12): Щ = 100, к е [0; 0,3], е < 0,0015, У™« = = 20.

Анализ результатов эксперимента (рис. 8 -вершина «треугольника» направлена в область неустойчивых возмущений) показал, что течение жидкой пленки неустойчиво во всем диапазоне чисел Рейнольдса Яе < 12. Чем больше значение числа Рейнольдса Яе, тем шире зона неустойчивости по диапазону волновых чисел к.

Области неустойчивости жидкой пленки, рассчитанные в линейном (рис. 5) и нелинейном (рис. 8) приближениях, согласуются между собой, рассогласование в вычислениях А к не превышает 6 % (табл. 2).

Таблица 2

Сравнение расчетов областей неустойчивости в линейном и нелинейном приближениях

Яе Кривая нейтральной устойчивости Ак, %

Линейное приближение Нелинейное приближение

1 0,02714 0,02559 5,71

5 0,10378 0,10357 0,20

7 0,13737 0,13655 0,60

9 0,16937 0,16954 0,10

10 0,18491 0,18454 0,20

12 0,21525 0,21453 0,33

Поскольку в результате волнового взаимодействия на свободной поверхности тонкого слоя вязкой жидкости происходит обмен и перераспределение энергии между возмущениями [7], интересным представляет исследование поведения составляющих амплитуды А(х, г, /) (7).

Для изучения эволюции во времени амплитуд я, волнового возмущения (7) система дифференциальных уравнений (8)—(14) решалась методом Рун-ге-Кутты 4-го порядка. Изучение осуществлялось в окрестности кривой нейтральной устойчивости.

Исследование показало, что по всему диапазону чисел Рейнольдса Яе в окрестности кривой нейтральной устойчивости наблюдается рост амплитуды несущей волны а0 (рис. 9) и амплитуд наклонных волн а3, а4. Амплитуды боковых волн аг (рис. 10) и а2 вблизи границы области неустойчивости затухают.

С увеличением числа Рейнольдса Яе увеличивается скорость роста амплитуды несущей волны а0 (рис. 9) и амплитуд наклонных волн а3, аЛ; увеличиваются скорости затухания амплитуд а3, а4 (рис. 10).

Рис. 6. График коэффициента Рі:

1 -Яе = 5;

2-Яе = 9\

Ъ-Яе= 10

Заключение

Исследование неустойчивости жидких пленок, подробное изучение процессов, протекающих на поверхности раздела газ-жидкость, актуально, так как имеет большое научное и прикладное значение.

Рис. 7. График коэффициента р2:

1 -Яе = 5\

2-Яе = 9\

3-Ле = 10

Рис. 8. Область неустойчивости в нелинейном приближении

Рис. 9. Поведение а0 в окрестности кривой нейтральной устойчивости:

1 -Яе = 5;

2-Яе-9;

2-Яе= 10

кривой нейтральной устойчивости:

1 -Яе = 5;

2-Яе = 9;

3-Ле= 10

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

В рамках математической модели неустойчивых режимов течения жидких пленок, представленной дифференциальным уравнением (1), осуществлено исследование неустойчивости тонкого слоя вязкой жидкости для диапазона чисел Рейнольдса Яе < 12 в линейном и нелинейном приближениях, проведен вычислительный эксперимент для воды (т = О, ЛГ= О, М= 0).

В линейном приближении получено дисперсионное уравнение (2), выведены аналитические зависимости волновых характеристик (оа„ ©/, С, Сгр) от физико-химических факторов (/V, т, а, М), рассчитана область неустойчивости (рис. 5) в плоскости параметров (Яе, к). Показано, что течение жидкой пленки (воды) неустойчиво в диапазоне чисел Рейнольдса Яе < 12, причем зона неустойчивости расширяется по диапазону волновых чисел к с увеличением значений числа Рейнольдса Яе. При этом увеличивается скорость роста волновых возмущений (рис. 2), уменьшаются фазовые (рис. 3) и групповые (рис. 4) скорости.

В нелинейном приближении в рамках системы дифференциальных уравнений (8—14) предложен алгоритм расчета области неустойчивости жидкой пленки. С помощью алгоритма вычислена кривая нейтральной устойчивости (рис. 8) для во-ды (т = 0, N = 0, М = 0). При этом рассогласование между результатами, полученными в линейном и нелинейном приближениях, не превышает 6% (табл. 2).

В рамках нелинейного приближения в окрестности кривой нейтральной устойчивости изучена эволюция во времени амплитуд составляющих возмущения (7). Показано, что по всему диапазону чисел Рейнольдса Яе наблюдается рост амплитуд ао, а3, а4 и затухание амплитуд аъ а2.

Литература

1. Прокудина, Л. А. Неустойчивость неизо-термической жидкой пленки / Л. А. Прокудина, Г. П. Вяткин // Доклады РАН, 1997. — №6. — С. 770-772.

2. Прокудина, Л. А. Волновое течение неизотермической жидкой пленки: препринт /

Л. А. Прокудина, Г. П. Вяткин. - Челябинск: Изд-воЮУрГУ, 1998.-42 с.

3. Трошенькш, Б. А. Циркуляционные и пленочные испарители и водородные реакторы / Б. А. Трошенькин. - Киев: Наукова Думка, 1985. -176 с.

4. Прокудина, Л. А. Моделирование неустойчивости жидкой пленки в тепломассообменных процессах при волнообразовании /Л. А. Прокудина, Е. А. Саламатов // Механика и процессы управления: труды XXXVI Уральского семинара. - Екатеринбург: УрО РАН, 2006. -Т. 1,- С. 183-191.

5. Свидетельство РФ об официальной регистрации разработки в ОФАП № 5392. Область неустойчивости жидкой пленки /Л. А. Прокудина, Е. А. Саламатов; зарегистрировано 18.11.2005.

6.Холпанов, Л. П. Гидродинамика и тепломассообмен с поверхностью раздела / Л. П. Холпа-нов, В. Я. Шкадов. - М.: Наука, 1990. - 271 с.

7. Прокудина, Л. А. Моделирование нелинейных явлений в жидкой пленке / Л. А. Прокудина // Известия Челябинского научного центра. - 2001. -№4. -С. 6-9.

8. Румянцев, В. В. Устойчивость и стабилизация движения по отношению к части переменных /В. В. Румянцев, А. С. Озиранер. - М.: Наука, 1987. -256 с.

9. Зубов, В. И. Об одном новом методе построения областей устойчивости в пространстве допустимых значений параметров системы автоматического регулирования / В. И. Зубов // Автоматика и телемеханика. - 1959. — № 3. -С. 331-334.

Ю.Дидук, Г. А. Методы теории матриц и га применение для исследования и проектирования систем управления / Г. А. Дидук. - СПб: Изд-во Санкт-Петербургского ун-та, 1993. -320 с.

11. Саламатов, Е. А. Математическое моделирование неустойчивости жидких пленок в тепломассообменных процессах с учетом нелинейных эффектов / Е. А. Саламатов, Л. А. Прокудина // Актуальные проблемы современной науки: труды 2-го международного форума. Естественные науки: ч. 1-3: Математика, математическое моделирование, механика. - Самара: СамГТУ, 2006. -С. 88-91.

Поступила в редакцию 1 февраля 2008 г.