Объект: HT01NN

Библиотека: HydroThermo2

Имя на уровне решателя: HT01NN

Аннотация: Расчёт мощности теплообменника. Ku = const * Qm^pow

Обозначение:

Порты (степени свободы) компонента:

№	Обозначение порта	Тип	Наименование порта
1	In1	THYDRAULIC	Теплогидравлический порт, 1 контур, вход
2	Out1	THYDRAULIC	Теплогидравлический порт, 1 контур, выход
3	In2	THYDRAULIC	Теплогидравлический порт, 2 контур, вход
4	Out2	THYDRAULIC	Теплогидравлический порт, 2 контур, выход

Пользовательские параметры модели:

№	Параметр	Тип	Описание	Единица	Значение по умолч.
1	A1	real	Площадь поверхности теплоотдачи 1 контура, м^2	m**2	1
2	A2	real	Площадь поверхности теплоотдачи 2 контура, м^2	m**2	1
3	Cd1	real	Коэффициент расхода		5
4	Cd2	real	Коэффициент расхода		5
5	F1	real	Суммарная проходная площадь сечения трубок 1 контура, м2	m**2	0.01
6	F2	real	Суммарная проходная площадь сечения трубок 2 контура, м2	m**2	0.01
7	Fluid1		Параметры среды 1 контура
8	Fluid2		Параметры среды 2 контура
9	K1	real	Коэффициент в формуле теплоотдачи, 1 контур, Вт/(м^2*К*(кг/с)^n1)		12.0
10	K2	real	Коэффициент в формуле теплоотдачи, 2 контур, Вт/(м^2*К*(кг/с)^n1)		12.0
11	L1	real	Длина трубок 1 контура, м	m	1
12	L2	real	Длина трубок 2 контура, м	m	1
13	PT_init1	list	Инициализирующие значения портов 1 контура: P1, T1, P2, T2, Па, К		101326,323.15,101325,323.15
14	PT_init2	list	Инициализирующие значения портов 2 контура: P1, T1, P2, T2, Па, К		101326,293.15,101325,293.15
15	Qm_eps	real	Масштаб Qm при сглаживании невязок знтальпии при Qm = 0, кг/с	kg/s	0.01
16	R12	real	Тепловое сопротивление стенки, К/Вт	K/W	0.01
17	Type	string	Тип теплообменника [Parallel_flow_single_pass,Counter_flow_single_pass,Cross_flow_both_mixed,Cross_flow_both_unmixed,Parallel_flow_dT_method,Counter_flow_dT_method,Cross_flow_dT_method]		Parallel_flow_single_pass
18	dT_mult	real	Корректирующий множитель для перепада температур при перекрёстном потоке	K	0.95
19	dp_min	real	Уровень давления для начала сглаживания, Па	Pa	0.01
20	n1	real	Степень в формуле теплоотдачи, 1 контур		1.5
21	n2	real	Степень в формуле теплоотдачи, 2 контур		1.5
22	ro1	real	Плотность 1 контура, кг/м3	kg/m**3	1
23	ro2	real	Плотность 2 контура, кг/м3	kg/m**3	1

Расчетные и измеряемые величины:

№	Параметр	Описание	Единица
1	V_Pressure1	Давление на порте	Pa
2	V_SpecificEnthalpy1	Удельная энтальпия;	J/kg
3	V_Density1	Плотность рабочей среды;	kg/m**3
4	V_Temperature1	Температура рабочей среды на порте;	K
5	S_MassFlow1	Массовый расход на порте ;	kg/s
6	S_EnthalpyFlow1	Поток энтальпии на порте ;	W
7	V_Pressure2	Давление на порте	Pa
8	V_SpecificEnthalpy2	Удельная энтальпия;	J/kg
9	V_Density2	Плотность рабочей среды;	kg/m**3
10	V_Temperature2	Температура рабочей среды на порте;	K
11	S_MassFlow2	Массовый расход на порте ;	kg/s
12	S_EnthalpyFlow2	Поток энтальпии на порте ;	W
13	V_Pressure3	Давление на порте	Pa
14	V_SpecificEnthalpy3	Удельная энтальпия;	J/kg
15	V_Density3	Плотность рабочей среды;	kg/m**3
16	V_Temperature3	Температура рабочей среды на порте;	K
17	S_MassFlow3	Массовый расход на порте ;	kg/s
18	S_EnthalpyFlow3	Поток энтальпии на порте ;	W
19	V_Pressure4	Давление на порте	Pa
20	V_SpecificEnthalpy4	Удельная энтальпия;	J/kg
21	V_Density4	Плотность рабочей среды;	kg/m**3
22	V_Temperature4	Температура рабочей среды на порте;	K
23	S_MassFlow4	Массовый расход на порте ;	kg/s
24	S_EnthalpyFlow4	Поток энтальпии на порте ;	W

Пример использования компонента

Модель теплообменника (HT01NN)

Зависимость коэффициента теплоотдачи для каждого контура принята следующей: Alpha = C * Qm^n. К нулю зависимость не продлевалась. Элемент имеет два контура. Гидравлика — Qm = sqrt(dp / Rh). Используется в термодинамических и гидравлических расчётах для моделирования теплообменных аппаратов, радиаторов, охладителей и нагревателей.

Теплообмен между водой и воздухом при постоянных перепадах давлений и расходах. Движение жидкости происходит в направлении от 1 Out к 1 In. Воздух охлаждает воду и нагревается. Температура воды уменьшается и выходит на постоянное значение. Знак расхода на первом контуре отрицательный.

Начальные данные:

Параметры модели (HT01NN2)

type = Parallel_flow_single_pass - Тип теплообменника

F1 = S - Площадь сечения 1 контура, м²

L1 = 0.1 - Длина 1 контура, м

F2 = S2 - Площадь сечения 2 контура, м²

L2 = 0.1 - Длина 2 контура, м

ro1 = 1000 - Плотность жидкости на входе в 1 контур (для расчёта Rh), кг/м^3

Cd1 = 10 - Приведённый коэффициент расхода, 1 контур

ro2 = 1.17 - Плотность воздуха на входе во 2 контур (для расчёта Rh), кг/м^3

Cd2 = 10 - Приведённый коэффициент расхода, 2 контур

F3 = 10 - Площадь поверхности теплоотдачи 1 контура, м²

F4 = 10 - Площадь поверхности теплоотдачи 2 контура, м²

K1 = 320.0 - Амплитудный коэффициент 1 контура

n1 = 1.5 - Степенной коэффициент 1 контура

K2 = 320.0 - Амплитудный коэффициент 2 контура

n2 = 1.01 - Степенной коэффициент 2 контура

R12 = 0.001 - Тепловое сопротивление стенки, К/Вт

dT_mult = 0.05 - Корректирующий множитель для случая перекрёстных потоков

PT_init1 = 101450, 323.15, 101326, 323.15 - Начальные значения давлений и температур на 1 и 2 портах

PT_init2 = 101450, 293.15, 101325, 293.15 - Начальные значения давлений и температур на 3 и 4 портах

Fluid1 = Water - Вода, 1 контур

Fluid2 = DryAir - Воздух, 2 контур

Параметры дросселей и площади

D = 0.05 - Гидравлический диаметр трубы, первый контур

D2 = 0.5 - Гидравлический диаметр трубы, второй контур

S = 0.25 * π * D * D - Площадь сечения первого контура

S2 = 0.25 * π * D2 * D2 - Площадь сечения второго контура

MDT = 0.72 - Коэффициент расхода

Параметры интегрирования

end = 15 с - Время интегрирования

control = 1e-9 - Точность контроля

dabs = 0.1 - Абсолютная погрешность

Схема тестирования:

Схема содержит теплообменник HT01NN2, включённый между двумя контурами. Первый контур (вода) подключён к граничным условиям HT1SPT1 (P = 101450 Па, T = 323.15 К) и HT1SPT3 (P = 101325 Па, T = 323.15 К) через дроссели HT1RSP1 и HT1RSP3. Второй контур (воздух) подключён к граничным условиям HT1SPT2 (P = 101450 Па, T = 293.15 К) и HT1SPT4 (P = 101325 Па, T = 293.15 К) через дроссели HT1RSP2 и HT1RSP4. Начальные параметры: PT_init1 = 101450, 323.15, 101326, 323.15; PT_init2 = 101450, 293.15, 101325, 293.15. Расходы через контуры измеряются датчиками Qm1 и Qm2. Температуры на входе и выходе каждого контура регистрируются как Tin_1, Tout_1, Tin_2, Tout_2. Результаты выводятся на осциллографы DISP1 (температуры) и RUN (расходы и температуры).

Результат:

На графике DISP1 показано изменение температур во времени: Tin_1 (вход воды) остаётся постоянной и равной 323.15 К; Tout_1 (выход воды) резко падает от 323.15 К и выходит на постоянное значение примерно 317.5 К (вода охлаждается); Tin_2 (вход воздуха) остаётся постоянной и равной 293.15 К; Tout_2 (выход воздуха) резко возрастает от 293.15 К и выходит на постоянное значение примерно 299.5 К (воздух нагревается). На графике RUN показано изменение расходов и температур: Qm1 (расход воды) и Qm2 (расход воздуха) практически равны нулю, что связано с малыми перепадами давлений; температуры Tin_1, Tout_1, Tin_2, Tout_2 соответствуют графикам DISP1. Это подтверждает корректность работы модели теплообменника: воздух охлаждает воду, вода нагревает воздух, температуры выходят на постоянные значения.

Формулы для расчёта:

\[\begin{split}\begin{aligned} \alpha &= C \cdot Q_m^{\,n} \\ Q_m &= \sqrt{\frac{\Delta P}{R_h}} \end{aligned}\end{split}\]

где: \(\alpha\) — коэффициент теплоотдачи, \(C, n\) — эмпирические коэффициенты (\(K_1, n_1\) для 1-го контура; \(K_2, n_2\) для 2-го контура), \(Q_m\) — массовый расход, \(R_h\) — гидравлическое сопротивление, \(\Delta P\) — перепад давления.

Тепловой поток:

\[Q = \frac{\Delta T}{R_{12}}\]

Схема теста

Рисунок 1 - Схема теста

Графики результатов теста

Рисунок 2 - Температуры HT01NN.Tin_1, HT01NN.Tout_1, HT01NN.Tout_2, HT01NN.Tin_2

Рисунок 3 - Совмещённый график RUN (HT01NN.Tin_1, HT01NN.Tout_2, HT01NN.Tin_2, HT01NN.Qm2, HT01NN.Qm1)