Смешанная конвекция в прямоугольной области с локальными источниками ввода и вывода массы в условиях неоднородного теплообмена Текст научной статьи по специальности «Физика»

УДК 536.2:532.5

СМЕШАННАЯ КОНВЕКЦИЯ В ПРЯМОУГОЛЬНОЙ ОБЛАСТИ С ЛОКАЛЬНЫМИ ИСТОЧНИКАМИ ВВОДА И ВЫВОДА МАССЫ В УСЛОВИЯХ НЕОДНОРОДНОГО ТЕПЛООБМЕНА

Г.В. Кузнецов, В.И. Максимов

Томский политехнический университет E-mail: [email protected]

Представлены результаты математического моделирования конвекции вязкой несжимаемой жидкости в прямоугольной полости с источниками ввода и отвода массы в сопряженной постановке в рамках приближения Буссинеска. Исследованы режимы смешанной конвекции в полости, имеющей две вертикальные и одну горизонтальную стенки конечной толщины, с двумя участками ввода и вывода жидкости, и одной свободной поверхностью. Рассмотрена плоская нестационарная задача в рамках модели Навье-Стокса для жидкой фазы и теплопроводности для твердой фазы. Получены распределения гидродинамических параметров и температур при использовании различных граничных условий на внешнем контуре рассматриваемой полости. Установлено, что даже при умеренном теплоотводе на внешних границах области решения приводит к масштабным изменениям в температурных полях жидкости.

1. Введение

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

2. Физическая модель

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

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

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

Рис. 1. Область решения задачи: 1) твердая фаза; 2,3) область ввода и вывода жидкости в полость; 4) полость с жидкостью; 5) свободная поверхность жидкости

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

Особое значение в задаче имеют граничные условия на внешних границах области решения (рис. 1). Эти условия определяют интенсивность теплоотвода с внешних границ и могут оказывать существенное влияние на режим течения жидкости. Соответственно, возможно изменение интенсивности теплопереноса в жидкости. Исходя из анализа различных вариантов реализации рассма-

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

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

При решении свободно-конвективных задач теплообмена чаще всего на внешних границах используются условия первого рода для уравнения энергии [7—10]. Такой подход удобен при численной реализации и в целом не выглядит противоречивым. Но анализ типовых примеров подобных задач показывает, что условия первого рода в целом почти адекватны условиям теплоизоляции. По этой причине они не соответствуют реальным режимам теплообмена на практике и могут быть использованы только в качестве условий для модельных задач. Но в этом качестве они являются элементами эталонных задач, решения которых играют роль своеобразных «измерительных приборов» при оценке тех или иных значимых физических эффектов на внешних границах областей решения. В данной работе оценка масштабности влияния теплоотвода с внешнего контура также проводилась путем сравнения с результатами решения такой «эталонной» задачи с граничными условиями пер-

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

Следует также отметить, что моделирование теплоотвода с трех внешних («твердых») границ области решения проводилось путем использования граничных условий второго рода. При этом величина отводимого теплового потока д являлась параметром задачи и выбиралась исходя из реального диапазона его возможного изменения. Такой прием позволяет не расширять область решения до размеров, при которых с одной стороны было бы правомерным применение условий четвертого рода на границах сопряжения разнородных сред. Но с другой стороны такое расширение области анализа делало практически невозможным сам численный анализ рассматриваемой нелинейной задачи. Выбор же значений д всегда может быть проведен с требуемой для практики или для обоснования достоверности результатов теоретического анализа точностью. Так, например, такой выбор может быть реализован по итерационной схеме путем расчета на каждой последующей итерации величины д, исходя из полученного на предыдущей итерации градиента температур вблизи границы. Кроме этого достаточно точно могут быть определены величины д, исходя из исходных данных каждой задачи. Последнее создает объективные условия для получения достоверных результатов.

3. Математическая модель и метод решения

Постановка задачи по пространственным координатам является двухмерной и описывается системой нестационарных уравнений Навье-Стокса для жидкой фазы и уравнением теплопроводности для твердой фазы [7, 8].

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

X=£, у = У, т= 1 и=±г V = А Ь Ь Г0 V¡n V,,

ДТ ’ ф0’

ДТ = Г, - То, / = V ь = р

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

Соответственно безразмерные уравнения Навье-Стокса в приближении Буссинеска в пере-

менных «вихрь скорости - функция тока - темпе ратура» для жидкой фазы (режим смешанной кон векции) и уравнение теплопроводности для твер дой фазы будут иметь вид:

дО ггдО ггдО 1 (д2О д2О ^ Ог (д0

дт

дХ дУ

д0 ттд0 тгд0>

----+ и----------+ V —

дт дХ дУ

+ и ™+ Яе I дХ2 + дУ2 Г Яе2 [дХ } ’(1)

1

Яе- Рг 15Х2

5 2Т 52 Т = О

д 20 д20 +

дУ/

ЭХ2 д0

дУ2

____=д0 де

д Бо = дХ2 дУ2 ’

(2)

(3)

(4)

здесь вг =

- число Грасгофа; в - темпера-

сТ

дХ

= 0,

д2Т дУ2

= 0,

д0 у (Х, У) д У

= Б1-0, (Х, У) + Бг

Т - Т

• 1П -* о

Тп - Т0

на внутренних границах раздела твердой и жидкой фаз, параллельных осям 0 7(0Х): дТ

Т = о, ---------= о,

дХ (У)

0 =0 д0

дХ (У)

/

= Я

Ч /

д0 / ; дХ (У)

- на участке ввода жидкости: 0 = 1, дТ

дУ

= 1;

- на участке оттока жидкости: д0( Х, У)

дУ

дТ

дУ

= 0,

= 1;

здесь К = -

qЬ

аЬ

турный коэффициент объемного расширения; g -ускорение, создаваемое массовыми силами;

2^ „V

Яе =----- - число Рейнольдса; Рг = — - число

V а

01

Прандтля; Бо = -Ь- - число Фурье; а - коэффициент температуропроводности.

Начальные условия:

Т( Х, У ,0) = 0,

О(Х, У ,0) = 0,

0( Х, У ,0) = 0.

Граничные условия:

- на внешнем контуре рассматриваемой области (кроме свободной поверхности жидкости) в одном варианте задаются граничные условия первого рода 0=0, а в другом - второго рода

д0 = Кц

дп

- на свободной поверхности жидкости:

- число Кирпичева; Б1 =

Я (Тп - Т) я

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

Задача (1-4) с соответствующими граничными и начальными условиями решена методом конечных разностей [7].

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

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

4. Анализ результатов численного моделирования

Численные исследования были проведены при следующих значениях:

- безразмерных величин: Рг=7,1, 100^е<1000, Gг<105;

- температур Т0=293 К, 300<Тп<343 К, 273<Т<293 К.

Диапазоны изменения безразмерных величин и

размерных характеристик выбирались исходя из следующих условий. Первое - соответствие диапазонам изменения параметров в реальных вариантах реализации таких систем, как «жидкость - твердое тело - газообразная внешняя среда». Второе - реальные масштабы существования как вынужденной, так и естественной конвекции без вырождения какого-либо из механизмов теплопереноса. Кроме того, контролировались условия геометрического подобия реальных систем и используемой при математическом описании процесса модели. Также рассматривались не самые благоприятные для реализации «сопряженного» теплопереноса [8] диапазоны изменения как размерных величин, так и безразмерных комплексов. Так, например, увеличение температуры вводимой в рассматриваемую полость жидкости, например, воды на 20...30° могло существенно изменить (при прочих равных условиях) соотношение между характерными временами реализации процесса вынужденной и естественной конвекции, но экстремальные варианты не рассматривались.

2

V

Рис. 2. Структура течения (а, в) и поле температур (б, г) для модели смешанной конвекции, когда на внешнем контуре, за исключением верхней свободной границы, заданы различные граничные условия: а, б) задано граничное условие первого рода Т=293 К, в, г) имеет место теплоотвод q=3 Вт/м2. Рг=7,1, Gr=105, Ре=1000. Размеры приведены в м, температуры - в К

На рис. 2 приведены типичные результаты решения системы (1-4) при различных граничных условиях на внешнем контуре.

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

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

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

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

Интересным является то, что в условиях достаточно интенсивного теплоотвода с нижней горизонтальной и с вертикальных границ значительно, по сравнению с первым вариантом (рис. 2, б), снижаются перепады температур во всех зонах области решения. Происходит выравнивание температуры в объеме. Более того, при наличии потерь тепловой энергии жидкость прогревается более равномерно при значительном снижении размеров области максимальных температур вблизи зоны ввода «горячей» жидкости. Полученный результат является следствием совместного протекания процессов естественной и вынужденной конвекции. «Горячая» жидкость фактически перетекает вдоль дна

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

5. Заключение

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

СПИСОК ЛИТЕРАТУРЫ

1. Водохранилища и водооградительные сооружения ГАЭС, ТЭС и АЭС / Под ред. Т.П. Доценко. - М.: Энергоатомиздат, 1989. -192 с.

2. Попов В.М., Рябцев В.И., Рябцев ГА. О тепловом загрязнении окружающей среды выбросами и сбросами ТЭС и ТЭЦ // Безопасность жизнедеятельности. - 2002. - № 4. - С. 26-28.

3. Rosengarten G., Graham L. Morrison Mixed convection in a narrow rectangular cavity with bottom inlet and outlet // International Journal of Heat and Fluid Flow. - 2001. - V. 22. - № 2. -P. 168-179.

4. Rossi di Schio B.E., Zanchini E. Combined forced and free flow in a vertical rectangular duct with prescribed wall heat flux // International Journal of Heat and Fluid Flow. - 2003. - V. 24. - № 6. -P. 874-887.

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

Результаты работы показывают возможности использования уравнений Навье-Стокса в переменных «вихрь скорости - функция тока» для моделирования достаточно сложных по своему характеру течений при умеренных значениях чисел Рейнольдса в режиме смешанной конвекции и при неоднородном теплообмене на внешних границах области решения.

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

Работа выполнена при финансовой поддержке РФФИ и Администрации Томской области (грант № 05-02-98006 конкурс робьа).

5. Иваненко С.А., Корявов П.П. Динамика вод и распространение загрязняющих веществ в водохранилище // Математическое моделирование. - 2002. - Т 14. - № 6. - С. 105-118.

6. Алексеев Л.П. Изучение и методы расчета крупномасштабной турбулентности глубоководного водоема: по материалам исследований на оз. Байкал. - Л.: Гидрометеоиздат, 1989. - 128 с.

7. Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена. - М.: Наука, 1984. - 288 с.

8. Лыков А.В., Алексашенко А.А., Алексашенко В.А. Сопряженные задачи конвективного теплообмена. - Минск: Наука и техника, 1971. - 346 с.

9. Джалурия Й. Естественная конвекция: Тепло- и массообмен. -М.: Мир, 1983. - 400 с.

10. Роуч П. Вычислительная гидродинамика. - М.: Мир, 1980. -616 с.