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

Моисеев, Константин Валерьевич

Влияние функциональной зависимости вязкости от температуры на свободную конвекцию жидкости тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Моисеев, Константин Валерьевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Тюмень МЕСТО ЗАЩИТЫ

2009 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

Диссертация по механике на тему «Влияние функциональной зависимости вязкости от температуры на свободную конвекцию жидкости»

Автореферат диссертации на тему "Влияние функциональной зависимости вязкости от температуры на свободную конвекцию жидкости"

На правах рукописи

МОИСЕЕВ КОНСТАНТИН ВАЛЕРЬЕВИЧ

ВЛИЯНИЕ ФУНКЦИОНАЛЬНОЙ ЗАВИСИМОСТИ ВЯЗКОСТИ ОТ ТЕМПЕРАТУРЫ НА СВОБОДНУЮ КОНВЕКЦИЮ ЖИДКОСТИ

01.02.05 — Механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Тюмень — 2009

003473259

Работа выполнена в Институте механики Уфимского научного центра Российской академии наук

Научный руководитель: доктор физико-математических наук

Урманчеев Сайд Федорович

Научный консультант: кандидат физико-математических наук

Ильясов Айдар Мартисович Официальные оппоненты: доктор физико-математических наук

Родионов Сергей Павлович

кандидат физико-математических наук, доцент Гореликов Андрей Вячеславович Ведущая организация: Институт теплофизики СО РАН

Защита диссертации состоится « 19 » июня 2009 года в 15°° часов на заседании диссертационного совета ДМ 212.274.09 в Тюменском государственном университете по адресу: 625003, г. Тюмень, ул. Перекопская, 15А.

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета.

Автореферат разослан « 18 » мая 2009 года.

Ученый секретарь диссертационного совета, к.ф,-м.н., доцент ' Мусакаев Н. Г.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

тельных средств.

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

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

Научная новизна.

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

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

• Исследовано влияние угла наклона ячейки относительно горизонта на свободную конвекцию жидкости с немонотонной зависимостью вязкости.

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

Работа выполнена в рамках программы исследований по гранту Президента Российской Федерации для государственной поддержки ведущих научных школ РФ (грант НШ-3483.2008.1) и при финансовой поддержке Программы фундаментальных исследований

ОЭММПУ РАН «Динамика многофазных и неоднородных жидкостей» (2006-2008 гг.).

Результаты работы использовались при составлении отчета о НИР Института механики УНЦ РАН № 01200614458 инв. № 02.2.00 951105 за 2009 г.

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

Апробации работы. Основные положения и результаты, представленные в диссертации, докладывались на следующих конференциях:

• Российская научно-техническая конференция «Мавлютовские чтения», посвященная 80-летию со дня рождения чл.-корр. РАН, проф. P.P. Мавлютова, Уфа, 2006 г.

• Четвертая Российская национальная конференция по теплообмену (РНКТ-4), Москва, 2006 г.

• Российская конференция «Механика и химическая физика сплошных сред», Бирск, 2007 г.

• Международная конференция «Потоки и структуры в жидкостях», Санкт-Петербург, 2007 г.

• IV Всероссийская конференция «Актуальные проблемы прикладной математики и механики», посвященная памяти академика А.Ф. Сидорова, Абрау-Дюрсо, 2008 г.

Кроме того, результаты работы докладывались на научных семинарах в Институте механики Уфимского научного центра РАН под руководством профессора С. В. Хабирова, профессора В. Ш. Шагалова и в Тюменском филиале Института теоретической и прикладной механики им. С. А. Христиановича СО РАН под руководством профессора А. А. Губайдуллина.

В 2009 г. работа была удостоена гранта Республики Башкортостан на конкурсе работ молодых ученых и молодежных научных коллективов.

Публикации. Основные результаты исследований опубликованы в 6 печатных работах, список которых представлен в конце автореферата.

Объем и структура работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа содержит 113 страниц, 7 таблиц и 80 рисунков. Список литературы включает 82 наименования работ российских и зарубежных авторов.

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

Во второй главе описан численный метод и представлены результаты тестовых расчетов. Дано подробное описание метода контрольного объема и алгоритма SIMPLE (Semi - Implicit Method for Pressure - Linked Equations). Проверка используемого в работе численного метода была проведена на ряде тестовых расчетов: течение изотермической жидкости в канале, свободная и вынужденная конвекции в квадратной ячейке. Оценена погрешность численного метода на вложенных сетках по правилу Рунге.

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

зонтальные — изотермическими.

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

ди ду

ди ди ди дЬ дх ду

_3_

дх дх

2 Мб)

( дь

ди\ дх)

(2)

дv дь дь др д _| ^^ | | ^ ^^^^ — ___ ^^__ |, ^

дЬ дх ду ду дх

д_ ду

+ Сг • б,

(3)

— +и— V— — — — (—] —

дЬ дх ду Рг дх \дх) Рг ду \ду

(4)

где вил - горизонтальная и вертикальная составляющие вектора скорости; р — отклонение от гидростатического давления; в — температура и ц(в) — динамическая вязкость как функция температуры.

¡/0

Параметры подобия: число Прандтля Рг = — и число Грасгофа

Хо

д/3 (Гг — I? Мо

ьг = -5-, где 1Уо — — и Уо — характерные значения

^о Ро

кинематической вязкости и коэффициента температуропроводности;

до и ро — характерные динамическая вязкость и плотность; 7\ и Тг —

размерные температуры верхней и нижней границ ячейки (22 > Тг);

(3 — коэффициент объемного расширения; д — ускорение свободного

падения; — характерный размер ячейки. Кроме того, при изучении

процессов свободно-конвективных течений используем число Рэлея

11а = Рг ■ Сг.

Полагалось, что функция вязкости в размерном виде д(Т) определена на интервале температур Т\ < Т < г1\, а область значений разбивалась на область больших цо < ц < /хкир и меньших

Иш < /л < /К) по отношению к Цо вязкостей. Неоднородность значений вязкости задавалась параметром г, который в области больших вязкостей полагался равным г = а в области меньших

ВЯЗКОСТеЙ — Г = Цо/.

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

Монотонные зависимости вязкости от температуры моделировались показательными функциями.

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

(а) (б)

Рис. 1. Графики исследуемых зависимостей в области больших вязкостей (а) и в области меньших вязкостей (б)

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

А»2(б) = 4-(1 -г)-в2 +Г,

Мб(б)=г(в+0,5))

ß3{ß) =4-(г-1)т"1 -в2+г-\ щ(в) = 4-(1-г)т-1-02 + 1,

/i7(ö) = г

М&)

_ „-(0+0,5)

Дв-0,5)

Полагалось, что в начальный момент времени жидкость находилась в состоянии покоя:

<=0 = 0, «|(=0 = 0, в\,=0 = 0.

Для компонент вектора скорости на всех границах задавались условия прилипания (и = 0 и и = 0). Граничные условия для температуры в соответствии с постановкой задачи записывались в виде:

дв

дх

= 0,

<?U = -0,5, ö|w=0 =0,5.

Для исследования интенсивности теплообмена на изотермических границах вычислялись числа Нуссельта:

NuH

где Nuh и Nue — числа Нуссельта на подогреваемой и охлаждаемой границах соответственно.

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

№)1 =

оо

J Nu(i)<

>iu,t dt

где t vi ш — безразмерное время и безразмерная круговая частота соответственно.

Все расчеты проводились на сетке 50 х 50 контрольных объемов для фиксированных параметров Рг = 118 и г = 10. Числа Грасгофа варьировались в пределах 1 < Gr < 40000.

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

зависимостью вязкости от температуры вида Mi (б) обнаружено существование трех типов стационарных течений. При последовательном увеличении числа Грасгофа сначала появляется одновихревой (рис. 2(а)), затем двухвихревой симметричный (рис. 2(6)), и наконец, двухвихревой асимметричный (рис. 2(в)) режимы, в то время, как для жидкости с постоянной вязкостью при заданном числе Прандтля (Рг = 118) выявлен только стационарный одновихревой режим.

(а) (б) (в)

Рис. 2. Линии тока в установившихся режимах конвекции: одновихревом (а), двухвихревом симметричном (б) и двухвихревом асимметричном (в). Зависимость вязкости от температуры вида /и (б)

Для жидкости с немонотонной зависимостью вязкости от температуры вида Цз(в), как и для жидкостей с монотонными температурными зависимостями вязкости — ц${в), Цо(0), ¡17(в) и /^(б) обнаружено существование двух установившихся режимов. Свободная конвекции начинается с одновихревого течения, а при увеличении числа Грасгофа происходит переход к двухвихревому симметричному режиму.

Немонотонные зависимости вязкости от температуры с отрицательной второй производной и ц.\{0) качественных изменений в картины установившихся течений свободной конвекции не вносят. Для жидкостей с этими функциональными зависимостями, как и для жидкости с постоянной вязкостью цо, обнаружены стационарные одновихревые режимы свободной конвекции.

Для жидкостей с монотонными зависимостями вязкости от тем-

пературы ^(О), Це{0), и /хв(^) установлено, что в одновихревом режиме конвекции для монотонно возрастающих функций вязкости На(в) и [н' > 0) циркуляция происходит против часовой стрелки, а для монотонно убывающих /л^(в) и ц?(в) (ц! < 0) — по часовой. Также обнаружено влияние знака первой производной функции вязкости на направления вихрей в двухвихревом симметричном режиме. Для жидкостей с монотонно возрастающими зависимостями выявлены течения с нисходящим движением в центре ячейки, а для жидкостей с монотонно убывающими зависимостями — с восходящим движением в центре.

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

(а)

(б)

Рис. 3. Изменение чисел Нуссельта со временем (а) и спектр мощности числа Нуссельта на подогреваемой границе (б) в автоколебательном режиме конвекции. Зависимость вязкости от температуры вида 1^2 (0)

Увеличение числа Грасгофа приводит к потере устойчивости стационарных течений и возникновению автоколебательных режимов свободной конвекции (рис. 3). Колебания тепловых потоков вызва-

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

0.6

_0.4

о.:

Рис. 4. Спектр мощности числа Нуссельта на подогреваемой границе (а) и проекция фазового пространства (и, v, du/dt, dv/dt) в фиксированной точке ячейки (ж = 0,5; у = 0,04) на двумерную плоскость {v, dv/dt) (б) в квазииериодическом режиме конвекции. Зависимость вязкости от температуры вида /Х4 (в)

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

(а) (б)

Рис. 5. Линии тока (а) и изменение чисел Нуссельта со временем (б) в хаотическом режиме конвекции. Зависимость вязкости от температуры вида д5(б)

(а) (б)

Рис. 6. Спектр мощности числа Нуссельта на подогреваемой границе (а) и проекция фазового пространства (и, V, ¿и/сИ, (IV/(1,1) в фиксированной точке ячейки (х = 0,5; у = 0,04) на пространство (и, йи/сй, йг//ай) (б) в хаотическом режиме конвекции. Зависимость вязкости от температуры вида

МО)

На рис. 7 представлено изменение числа Нуссельта на подогреваемой границе ячейки в зависимости от числа Рэлея для жидкостей с немонотонными (рис. 7(а)) и монотонными (рис. 7(6)) температурными зависимостями вязкости, сплошной линией обозначена кривая Ми(11а) для жидкости с постоянной вязкостью до- Видно, что интенсивность теплообмена для жидкостей с зависимостями вязкости от температуры Цг{6), /12(6), и Мв(в) ниже, а для жидкостей с функциями вязкости /из(@), 14{0), 14(0) и ¡¿&(в) выше, чем интенсивность теплообмена для жидкости с постоянной вязкостью цо. Также нетрудно заметить, что кривые Ми(11а) для жидкостей с монотонными функциями вязкости /15(6!) {р!ь < 0) и (/4 > 0), так же, как и для жидкостей с зависимостями (¡л'7 < 0) и (х&(в) (¡л'& > 0), практически совпадают. То есть для жидкостей с монотонными зависимостями вязкости от температуры интенсивность теплообмена не зависит от знака первой производной функции вязкости по температуре. Для жидкостей с немонотонными зависимостями картина представляется более сложной.

I'5

: ......— ~~ ц4(в)

. ------ — Ц,(в)

М-о

--- ИЛ8)

- --- ......

, , . . 1 , , , , 1

1.5x10'" „ 3.0x10"' 4.5х!0,°"

Ка (б)

Рис. 7. Изменение числа Нуссельта на подогреваемой границе ячейки в зависимости от числа Рэлея для жидкостей с немонотонными (а) и монотонными (б) зависимостями вязкости от температуры

Для лучшего понимания влияния зависимости вязкости на теп-

1.5x10 „ 3.0x10* 4.5*10"

Ка (а)

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

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

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

Четвертая глава посвящена численному исследованию влияния угла наклона квадратной ячейки относительно горизонта на свободную конвекцию жидкости с квадратичной зависимостью вязкости от температуры с минимумом 1Л\(в). Верхняя и нижняя границы ячейки считались изотермическими, боковые — адиабатическими. Полагалось, что жидкость подогревается снизу.

В качестве математической модели использовалась безразмерная система уравнений свободной конвекции (1)-(4), в которой отклонение ячейки на угол а учитывалось в проекциях уравнений движения Навье - Стокса на оси координат:

тур

Т1

ди ди

ди др д ду дх дх

+ Сг • в ■ вта.

дм ду ду

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

Расчеты проводились на сетке 50 х 50 контрольных объемов для значений параметров Рг = 118, г = 10 и фиксированных углов а = 15°, 30°, 45°, 60°, 75°, 90°. Числа Грасгофа изменялись в интервале 1 < вг < 40000.

В отличие от угла наклона а — 0°, при а = 15° было обнаружено существование двух типов стационарных режимов свободной конвекции. При последовательном увеличении числа Грасгофа сначала появляется одновихревой режим, затем его сменяет одновихревой с вложенными внутренними вихрями. Дальнейшее увеличение чисел Грасгофа приводит к возникновению автоколебательных режимов. Установлено существование режима одночастотных колебаний, в котором числа Нуссельта на подогреваемой и охлаждаемой границе колеблются в противофазе (рис. 9). Также выявлен режим двухча-стотных синфазных колебаний (рис. 8).

При а = 30°, 45°, 60°, 75°, 90° обнаружены только стационарные одновихревые течения и стационарные течения с вложенными внутренними вихрями. С ростом угла наклона колебаний интегральных тепловых потоков на изотермических стенках не обнаружено.

На рис. 10(а) представлено изменение числа Нуссельта на подогреваемой границе в зависимости от числа Рэлея при различных углах наклона ячейки. Видно, что с ростом угла наклона интенсивность теплообмена повышается, однако для угла а = 90° становится меньше, чем для углов наклона а = 60°, 75°. На рис. 10(6) показано изменение числа Нуссельта на подогреваемой границе в зависимости от угла наклона ячейки при фиксированных числах Рэлея. Видно, что интенсивность теплообмена изменяется немонотонно и при угле наклона а яа 75° достигает своего максимального значения.

Рис. 9. Изменение числа Нуссельта со временем (а) и спектр мощности числа Нуссельта на подогреваемой границе ячейки (б), 11а = 1888 • 103

1«

- Иа=23600

- Яа^НВООО

- Яа=295000 -- Яа-1062000 -- Ча=23 60000

(б)

Рис. 10. Зависимость числа Нуссельта на подогреваемой границе от числа Рэлея при фиксированных углах наклона ячейки (а) и зависимость чисел Нуссельта от угла наклона ячейки при фиксированных числах Рэлея (б)

Обнаружено, что с ростом угла наклона минимальные критические числа Рэлея убывают.

В заключении представлены основные результаты, полученные в работе.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Проведено численное исследование свободной конвекции жидкостей с монотонными и немонотонными зависимостями вязкости от температуры. Установлено, что вид зависимости вязкости от температуры влияет на стационарные режимы конвекции. Для жидкостей с монотонными зависимостями и немонотонными зависимостями с положительной второй производной были обнаружены установившиеся одновихревые и двухвихре-вые режимы конвекции. Для жидкостей с немонотонными зависимостями с отрицательной второй производной, как и для жидкости с постоянной вязкостью, — только установившиеся одновихревые режимы.

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

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

4. Установлено, что увеличение угла наклона ячейки относительно горизонта оказывает стабилизирующее воздействие на режимы свободной конвекции для жидкости с температурной зависимостью вязкости, как и для жидкости с постоянной вязкостью. При углах наклона ячейки а > 30° в исследуемом диапазоне чисел Грасгофа обнаружены только стационарные режимы свободной конвекции.

5. Обнаружено, что с ростом угла наклона ячейки минимальные критические числа Рэлея убывают, а интенсивность теплообмена изменяется немонотонно и имеет локальный максимум.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Численное моделирование термоконвекции жидкости с квадратичной зависимостью вязкости от температуры / А. М. Ильясов, К. В. Моисеев, С. Ф. Урманчеев // Сибирский журнал индустриальной математики. — 2005. — Т. VIII, № 4(24). — С. 51-59.

2. Численное исследование термогравитационной конвекции жидкости с квадратичной зависимостью вязкости от температуры / К. В. Моисеев, А. М. Ильясов // Мавлютовские чтения: Российская научно - техническая конференция: сб. трудов. — Уфа: УГАТУ, 2006. - Т. 4. - С. 74-79.

3. Численное моделирование конвекции термовязкой жидкости в квадратной полости / А. М. Ильясов, К. В. Моисеев // Труды четвертой российской национальной конференции по теплообмену: в 8 томах. Т. 3. Свободная конвекция. Тепломассообмен при химических превращениях. — М.: Издательский дом МЭИ,

2006. - С. 113-116.

4. Влияние угла наклона полости на теплообмен при свободной конвекции аномально термовязкой жидкости / К. В. Моисеев // Труды Института механики Уфимского научного центра РАН. Вып. 4 / Под ред. С. Ф. Урманчеева, С. В. Хабирова. — Уфа: Изд-во «Гилем», 2007. — С. 166-173.

5. Численное исследование свободноконвективных течений термовязких жидкостей / А. М. Ильясов, К. В. Моисеев, С. Ф. Урманчеев // Сборник докладов международной конференции «Потоки и структуры в жидкостях». — Санкт-Петербург. —

2007. - С. 233-236.

6. Численное исследование конвекции термовязкой жидкости с границей раздела / К. В. Моисеев, С. Ф. Урманчеев // Тезисы докладов IV Всероссийской конференции «Актуальные проблемы прикладной математики и механики», посвященной памяти академика А.Ф. Сидорова. — Екатеринбург: УрО РАН, 2008. — С. 45.

Моисеев Константин Валерьевич

ВЛИЯНИЕ ФУНКЦИОНАЛЬНОЙ ЗАВИСИМОСТИ ВЯЗКОСТИ ОТ ТЕМПЕРАТУРЫ НА СВОБОДНУЮ КОНВЕКЦИЮ ЖИДКОСТИ

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Лицензия на издательскую деятельность ЛР№ 021319 от 05.01.99 г.

Подписано в печать 18.05.2009 г. Формат 60x84/16. Усл. печ. л. 1,38. Уч.-изд. л. 1,44. Тираж 100 экз. Заказ 283.

Редакционно-издательский центр Башкирского государственного университета 450074, РБ, г. Уфа, ул. Заки Валиди, 32.

Отпечатано на множительном участке Башкирского государственного университета 450074, РБ, г. Уфа, ул. Заки Встиди, 32.

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Моисеев, Константин Валерьевич

Введение

Глава 1. Аналитический обзор литературы

1.1. Обзор литературы.

1.2. Математическая модель свободной конвекции жидкости с температурной зависимостью вязкости.

Глава 2. Численный метод и тестовые расчеты

2.1. Введение

2.2. Метод контрольного объема и процедура SIMPLE.

2.3. Течение изотермической жидкости в плоском канале.

2.4. Вынужденная конвекция.

2.5. Свободная конвекция.

2.6. Выводы.

Глава 3. Влияние типа температурной зависимости вязкости на свободную конвекцию несжимаемой жидкости

3.1. Введение.

3.2. Математическая модель и постановка задачи.

3.3. Параметры модельной задачи

3.4. Результаты численного исследования.

3.5. Выводы.

Глава 4. Влияние угла наклона па свободную конвекцию жидкости с температурной зависимостью вязкости

4.1. Введение

4.2. Математическая модель и постановка задачи.

4.3. Параметры модельной задачи

4.4. Результаты численного исследования.

4.5. Выводы.

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

Актуальность работы

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

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

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

Цель работы

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

Научная новизна

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

Достоверность результатов

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

Практическая значимость

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

• Четвертая Российская национальная конференция по теплообмену (РНКТ-4), Москва, 2006 г.

• Российская конференция «Механика и химическая физика сплошных сред», Бирск, 2007 г.

• Международная конференция «Потоки и структуры в жидкостях», Санкт-Петербург, 2007 г.

Кроме того, результаты работы докладывались на научных семинарах в Институте механики Уфимского научного центра РАН под руководством профессора С. В. Хабирова, профессора В. Ш. Шагапова и в Тюменском филиале Института теоретической п прикладной механики им. С. А. Хри-стиановича СО РАН под руководством профессора А. А. Губайдуллина.

Публикации

Основные результаты по теме исследований опубликованы в 6 печатных работах.

Объем и структура работы

Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа содержит 113 страниц, 7 таблиц и 80 рисунков. Список литературы включает 82 наименования работ российских и зарубежных авторов.

Заключение диссертации по теме "Механика жидкости, газа и плазмы"

4.5. Выводы

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

Установлено, что угол наклона ячейки оказывает стабилизирующее воздействие на свободную конвекцию жидкости с температурной зависимостью вязкости, как и для жидкости с постоянной вязкостью. При а > 30° в исследуемом диапазоне чисел Грасгофа наблюдались только установившиеся режимы.

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

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

Выявлено, что минимальные критические числа Рэлея убывают с ростом угла наклона.

Заключение

1. Проведено численное исследование свободной конвекции жидкостей с монотонными и немонотонными зависимостями вязкости от температуры. Установлено, что вид зависимости вязкости от температуры влияет па стационарные режимы конвекции. Для жидкостей с монотонными зависимостями и немонотонными зависимостями с положительной второй производной были обнаружены установившиеся одновихревые и двухвихревые режимы конвекции. Для жидкостей с немонотонными зависимостями с отрицательной второй производной, как п для жидкости с постоянной вязкостью, — только установившиеся одновихревые режимы.

Список источников диссертации и автореферата по механике, кандидата физико-математических наук, Моисеев, Константин Валерьевич, Тюмень

1. О.Н. Аликина, E.J1. Тарунип Подкрнтнческие движения жидкости в случае вязкости, зависящей от температуры// МЖГ. 2001. № 4. С. 55 - 62.

2. Андерсон Д., Таннехплл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. Т. 2. М.: Мир, 1990. 728 с.

3. Ахатов И.Ш., Вайков В.А., Хасанов М.М. Неустойчивость и хаос в гидродинамике. Учебное пособие. Уфа, 1991. 92 с.

4. Базаров И.П. Термодинамика. М.: Высшая школа. 1991. 376 с.

5. Бармин А.А., Мельник О.Э. Течение загазованной магмы в канале вулкана // Известия РАН, сер. МЖГ. 1990. № 5. С. 35 40.

6. Бармин А.А., Мельник О.Э.Об особенностях динамики извержения спльновязких газонасыщенных магм // Известия РАН, сер. МЖГ. 1993. № 2. С. 49 60.

7. Бармин А.А., Мельник О.Э.Гидродинамика вулканических извержений // Успехи механики. 2002. № 1. С. 32 60.

8. Бердников В. С., Гришков В. А. Ламинарно-турбулентный переход в свободноконвективном слое и теплоотдача вертикальных стенок//

9. Труды четвертой российской национальной конференции по теплообмену. Т.З. Свободная конвекция. Тепломассообмен при химических превращениях. М: Изд. МЭИ. 2006. С. 67 70.

10. Буссе Ф.Г. Переход к турбулентности в конвекции Рэлея Бена-ра // Гидродинамические неустойчивости и переход к турбулентности. М.: Мир, 1984. С. 124 - 168.

11. Ф. Буссе, Д.В. Любимов, Т.П. Любимова, Г.А. Седельников Трехмерные режимы конвекции в кубической полости // МЖГ. 2008. N° 1. С. 3-11.

12. Ван- ДайкМ. Альбом течений жидкости и газа. М.: Мир, 1986. 184 с.

13. Вельтищев Н.Ф. Конвекция в горизонтальном слое жидкости с постоянным внутренним источником тепла // МЖГ. 2004. № 2. С. 24 33.

14. Галиев И.М., Зубков П.Т. Влияние инверсии плотности воды на плоскопараллельное течение и теплоперенос в канале постоянной ширины // МЖГ. 2000. № 1. С. 72 78.

15. Гебхарт В., Джалурия И., Махаджан Р., Саммакия Б. Свободнокон-вективные течения, тепло и массообмен. М.: Мир, 1991. Т. 1. 700 с.

16. Гебхарт В., Джалурия И., Махаджан Р., Саммакия Б. Свободнокон-вективные течения, тепло и массообмен. М.: Мир, 1991. Т. 2. 528 с.

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

18. Гершуни Г.З., Жуховицкий Е.М., Шихов В.М. Об устойчивости конвективного течения жидкости, с вязкостью зависящей от температуры. // Теплофизика высоких температур. 1975. Т. 13, № 4. С. 771 -778.

19. Гершуни Г.З., Жуховицкий Е.М., Непомнящий А.А. Устойчивость конвективных течений. М.: Наука, 1989. 320 с.

20. Гетлинг А.В. Формирование пространственных структур конвекции Рэлея Бенара // УФН. 1991. Т. 161, № 9. С. 1 - 80.

21. Грязнов В.Л., Полежаев В.И. Численное моделирование турбулентного режима конвекции в вертикальном слое // МЖГ. 1977. № 5. С. 8 24.

22. Дайковский А.Г., Полежаев В.И., Федосеев А.И. Исследование структуры переходного и турбулентного режимов конвекции в вертикальном слое // МЖГ. 1978. № 6. С. 66 75.

23. Зубков П.Т., Калабин Е.В. Численное исследование естественной конвекции воды вблизи точки инверсии плотности при числах Грасгофа до 106 // МЖГ. 2001. № 6. С. 103 110.

24. Зубков П.Т., Калабин Е.В. Яковлев А.В. Исследование естественной конвекции пресной воды вблизи 4 °С в кубической полости // МЖГ. 2002. № 6. С. 3 10.'

25. Зубков П.Т., Климин В.Г. Численное исследование естественной конвекции чистой воды вблизи точки инверсии плотности // МЖГ. 1999. № 4. С. 171 176.

26. Идельчик И.Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1992. 464 с.

27. Ильясов А. М., Моисеев К. В., Урманчеев С.Ф. Численное моделирование термоконвекции жидкости с квадратичной зависимостью вязкости от температуры// Сибирский журнал индустриальной математики, 2005. Т.VIII, № 4(24). С. 51 59.

28. Ильясов А. М., Моисеев К. В., Урманчеев С.Ф.Численное исследование свободнокопвективных течений термовязких жидкостей// Сборник докладов международной конференции потоки и структуры в жидкостях. Санкт Петербург. 2007. С. 233 236.

29. В.Н. Коровкин, А.П. Андриевский Свободная конвекция при нелинейной зависимости плотности от температуры: плоские задачи// ИФЖ. 2000. Т. 73. № 2. С. 381 386.

30. Ладыженская О.А. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Физматлит, 1961. 204 с.

31. Ландау Л.Д., Лифшпц Е.М. Механика сплошных сред. М.: Гостех-издат, 1953. 788 с.

32. Ландау Л.Д., Лифшпц Е.М. Статистическая физика. М.: Наука, 1964. Т. 5. 567 с.'

33. Ланда П.С. Нелинейные колебания и волны. М.: Наука, 1997. 496 с.

34. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1987. 840 с.

35. Лоренц Э. Детерминированное непериодическое течение.// Странные аттракторы. М.: Мир, 1981. С. 88 116.

36. С.Л. Марпл мл. Цифровой спектральный анализ и его приложения. Пер. с англ. М.: Мир, 1990. 584 с.

37. Моисеев К. В., Ильясов А. М.Численное исследование термогравитационной конвекции жидкости с квадратичной зависимостью-вязкости от температуры// Сборник трудов российской научно технической конференции Мавлютовские Чтения. 2006. Т 4. С. 74 - 79.

38. Моисеев К. В. Влияние угла наклона полости на теплообмен при свободной конвекции аномально термовязкой жидкости// Труды Института механики Уфимского научного центра РАН. Вып. 4. Уфа: Гилем, 2007. С. 166 173.

39. Моргун Д.А. Численное исследование влияния инверсии плотности на конвекцию холодной воды в квадратной полости // Диссер. на соискание уч. степени к.ф. м.н. Уфа, 2002. 101 с.

40. Мун Ф. Хаотические колебания. М.: Мир, 1990. 312 с.

41. Пасконов В.М., Полежаев В.И., Чудов В.А. Численное моделирование процессов тепло и массообмена. М.: Наука, 1984. 288 с.48,49,50,51.