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

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

В.И. Лалетин ФГБОУ ВПО «ЮРГУЭС», г.Шахты

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

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

Рассмотрим общие закономерности, присущие рассматриваемой физической модели, описанные в [2]. В основу изучения движения вязкого газа положим следующие допущения.

1). Газ совершенен, т. е. давление «р», плотность «р » и абсолютная температура «Т» удовлетворяют уравнению состояния — закону Менделеева - Клапейрона:

р-Я  Т

Р = ^ ,

т

где: Я - универсальная газовая постоянная;

т - молекулярная масса газа.

2). Коэффициент удельной теплоемкости «с» не зависит от абсолютной температуры газа и является физической константой газа.

3). При истечении газа учитывается только вязкость первого рода (сопротивление окружающего теплого воздуха потоку холодного воздуха). Коэффициент теплопроводности газа «X» пропорционален коэффициенту динамической вязкости « ц », то есть выполняется критерий Прандтля:

а-с

—— = а = const.

X

4). К вышеописанным уравнениям присоединяется уравнение неразрывности движения:

др | д(ри) | 8(pv) | d(pw) 0

8t 8x 8y 8z

5). Газ представляет собой «ньютоновскую» среду, подчиненную известному обобщенному закону Ньютона о линейной связи между тензором напряжений и тензором скоростей деформаций.

Основные уравнения Навье — Стокса динамики вязкого газа, отнесенные к единице массы:

du ^ др

Р— = ----------+ д

dt х дх

dv ^ др

Р— = PFV------------------+ Д

dt ду

dw ^ др

Р-Г = Рр, -^ + Д

dt дг

д и д и д u

7 ^-------7 ^----7

удх ду дг у

д v д v д v

—7 ^-------7 ^----7

удх ду дг у

^ д2 w д2 w д2 w^

д — 2 д —

+ Д — (divV(д-divV);

дху / 3 дху /

д — 2 д —

+дду

д — 2 д —

+ д—(d„v у - (д. sw);

У дх ду дг у

здесь: и,V, - проекции вектора скорости V на оси ОХ, ОУ, 07;

^ = ( Рх; Ру; ^)- вектор внешних объемных сил, действующих на газ в

каждой точке пространства;

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

7. — ди дv

=-----1----1---

дх ду дг

<

du dv dw

Полные производные ~~Г,~Г,~Г можно представить в развернутой

dt dt dt

форме:

du ди ди ди ди

= ь — ил v-\ w;

dt дt дх ду дг

dv дv дv дv дv

= = ь — и v + — w;

dt дt дх ду дг

dw дw дw дw дw

= 1 и v w.

dt дt дх ду дг

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

^ =0; р'у =0; ^ ^ ,

где g - ускорение силы тяжести.

При сформулированных выше допущениях система уравнений, описывающая истечение воздуха из холодильной камеры, содержит семь неизвестных функций: и, v, w, р, p, T, д. В вышеописанной постановке задачи таких уравнений предложено шесть, в качестве седьмого уравнения для замыкания системы может быть предложено уравнение диффузии, либо уравнение теплового баланса.

С формулируем граничные условия задачи.

1). Равенство нулю скорости на неподвижной твердой границе внутри холодильного шкафа.

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

3). Задана температура воздуха внутри шкафа в начальный момент времени

истечения холодного воздуха из шкафа и равна Тхол = 50 С.

4). В начальный момент времени температура снаружи холодильного шкафа постоянна и равна Ткомн = 250 С .

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

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

При таких допущениях имеют место равенства:

Т = сот^, р = сот^, д = сот^.

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

При сделанных допущениях уравнение неразрывности выполняется.

В силу условия (8) два последних слагаемых в каждом уравнении системы (4) равны нулю, учитывая формулы (1), (4) и допущения (6), (7), систему (5) преобразуем к более простому виду:

ґ ди ди ди ди ' д и д2 и д2 Л д и

р И и + V + — w = Д + +

ч дї дх ду д2 у удх2 ду2 д2 у

ґ дv дv дv дv л > гм ГО д \ > гм ГО

р И и + — v + — w = Д + 9 И" 9

ч дї дх ду д2 у удх2 ду д2 у

ґ дw дw дw дw гм ГО д2 w гм го

р — — — и + V + — w = Ртё + Д — — 9 +

ч дї дх ду д2 у ч дх ду2 д22 у

После сделанных допущений система (9) есть система трех уравнений относительно трех неизвестных функций и, V, w, т. е. система (9) - замкнута.

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

1 - верхнетемпературная камера;

2 - низкотемпературная камера.

Рис. 1 - Поведение свободно истекающего потока воздуха в плоскости Х07 Предположение о неизменности геометрии потока в направлении оси 0У влечет равенство нулю значений функции V = у(х, у, г) = 0 и функции и = и(х, г); w = ^(х, г), то есть они не зависят от переменной « у ». Таким образом, реальный трехмерный поток заменяется моделью двухмерного.

В этих условиях система (9) упростится:

Р

Р

ди ди ди

\

------1-----ил---------w

V д1 дх дг у

= Д

^ д2 и д2 и^ дх2 дг2

дw дw дw

\

+

д1 дх

и

— w

дг у

{ \2 ^2 *4

д w д w

—Г + —Т

дх дг

<

Будем дополнительно предполагать, что процесс истечения воздуха стационарный по времени, тогда:

* - 0, ^ - 0

ді ді

и система (10) упростится к виду:

ґдп дп Л

----П Л-------

V дх дг у

ґ дw

р

2„ Л

д п д п

дх2

дг2

у

дw

— п Л-----w

дх дг у

И

= Ш£Л---

р

д w д w

—т ^--------т

у дх дг у

Так как система (12) содержит два уравнения и две неизвестные функции, то она является замкнутой.

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

И

Выражения — Р

2. \

д п д п , Ь--------:

дх

Р

2_Л

д w д w

дхА дг2

, входящие в оба уравнения

V ^ ^ у р ч ^ ^ у

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

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

Р

Р

2. \

д п д п

дх2 дzz

дzw д2w'''s

----Г +------г

дх дг

п2;

И к 2 ---w ,

<

здесь к - коэффициент пропорциональности, определяемый экспериментально.

Так как коэффициент динамической вязкости и имеет единицы Кг с

измерения------, то из уравнений (13) получаем, что к измеряется в —-.

м  с м

Знак минус в правых частях равенств (13) указывает на то, что вторые частные производные от проекций вектора скорости отрицательны из -за тормозящего влияния вязкой среды.

В силу формул (13) система (12) преобразуется к виду:

ґ дп дп Л

— п Л------w

V дх дг у

ґ дм дм Л

п Л-----w

дх

дг

Ик и2

И к 2 = mg-----------w .

У

Система (14) является базовой для нахождения скоростей ламинарного потока при его свободном истечении из холодильной камеры в стационарном режиме.

<

Литература.

1. Левкин, В.В. Тепловые расчеты сборочных единиц бытовых холодильников / В.В. Левкин. Шахты.: АП «Полиграфист», 1994. - 229 с.

2. Лойцянский, Л.Г. Механика жидкости и газа / Л.Г. Лойцянский. - 5-е изд. -М.: Наука, 1978. - 736 с.

3. Корн, Г. Справочник по математике для научных работников и инженеров/ Г. Корн, Т. Корн. - М.: Наука, 1970. - 720 с.