Моделирование теплообмена между сетью влажного воздуха и тепловой жидкостной сетью
Интерфейсы сети Simscape/Fluids/Fluid/Теплообменники

Блок теплообменника (TL-MA) моделирует теплообменник с одной сетью влажного воздуха, которая протекает между портами A2 и B2, и одной сетью термической жидкости, которая протекает между портами A1 и B1. Потоки текучей среды могут быть выровнены параллельно, встречно или поперечно.
Термический жидкостно-влажный воздушный теплообменник не подходит для холодильных систем охлаждения. См. Конденсаторный испаритель (2P-MA) или конденсаторный испаритель (TL-2P) для теплообменников, которые могут использоваться в холодильных установках.
Можно смоделировать сторону влажного воздуха как поток внутри труб, поток вокруг труб для термической жидкости или с помощью эмпирической общей параметризации. Влажная воздушная сторона содержит воздух, следовый газ и водяной пар, которые могут конденсироваться в течение всего цикла теплообмена. Модель блока учитывает скрытую теплоту, которая выделяется, когда вода конденсируется на поверхности теплопередачи. Этот жидкий слой не собирается на поверхности и, как предполагается, полностью удаляется из потока влажного воздуха, расположенного ниже по потоку. Скорость конденсации влаги возвращается как физический сигнал в порту W.
Блок использует метод Effectivity-NTU (E-NTU) для моделирования теплопередачи через общую стенку. Также моделируется загрязнение стенок теплообменника, которое увеличивает тепловое сопротивление и уменьшает теплообмен между двумя текучими средами. Можно также дополнительно моделировать ребра как на стороне влажного воздуха, так и на стороне термической жидкости. Потеря давления из-за вязкого трения с обеих сторон теплообменника может быть смоделирована аналитически или с помощью общей параметризации, которую можно использовать для настройки на собственные данные.
Можно моделировать сторону термической жидкости как поток внутри труб, поток вокруг труб влажного воздуха или с помощью эмпирической общей параметризации.
Эффективность теплообменника основана на выбранной конфигурации теплообменника, свойствах жидкости, геометрии трубы и конфигурации потока на каждой стороне теплообменника, а также использовании и размере ребер.
Параметр расположения потока задает относительные пути потока между двумя сторонами:
Parallel flow указывает, что текучие среды движутся в том же направлении.
Counter flow указывает, что жидкости движутся параллельно, но в противоположных направлениях.
Cross flow указывает, что жидкости движутся перпендикулярно друг другу.
Если для параметра «Расположение потока» установлено значение Cross flowиспользуйте параметр «Расположение поперечного потока», чтобы указать, разделены ли тепловые потоки жидкости или влажного воздуха на несколько путей перегородками или стенками. Без этих разделений поток может свободно смешиваться и считается смешанным. Обе жидкости, одна жидкость или ни одна из них не могут быть смешаны в устройстве поперечного потока. Перемешивание гомогенизирует температуру текучей среды вдоль направления потока второй текучей среды и изменяется перпендикулярно второму потоку текучей среды.
Несмешанные потоки изменяются по температуре как вдоль, так и перпендикулярно пути потока второй текучей среды.
Примеры кросс-поточных конфигураций



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

Параметр Геометрия потока (Flow geometry) задает расположение потока жидкости соответствующей диалоговой вкладки либо внутри трубы или набора труб, либо перпендикулярно пучку труб. Можно также указать эмпирическую общую конфигурацию.
Если для геометрии потока задано значение Flow perpendicular to bank of circular tubes, используйте параметр расположения сетки банка труб, чтобы определить выравнивание банка труб другой жидкости как Inline или Staggered. Красная стрелка, направленная вниз, на рисунке ниже указывает направление текучей среды, протекающей снаружи от пучка труб. На рисунке Inline также показано количество рядов труб вдоль направления потока и количество сегментов труб в параметрах каждого ряда труб. Здесь направление потока относится к текучей среде соответствующей диалоговой вкладки, а труба относится к трубке другой текучей среды. Параметр Длина каждого сегмента трубы в ряду труб показан на рисунке в шахматном порядке.


Только для одной жидкости может быть задано значение Геометрия потока (Flow geometry) Flow perpendicular to bank of circular tubes за раз. Другая жидкость должна быть сконфигурирована так, чтобы Flow inside one or more tubes или Generic. Если для геометрии потока для обеих жидкостей задано значение Flow perpendicular to bank of circular tubes, вы получите ошибку.
Конфигурация теплообменника не имеет ребер, если для параметра Общая площадь поверхности ребер установлено значение 0 m^2. Ребра вводят дополнительную площадь поверхности для дополнительной теплопередачи. Каждая сторона текучей среды имеет отдельную область ребер.
Скорость теплопередачи рассчитывают по усредненным свойствам обеих текучих сред.
Теплопередача рассчитывается как:
MA),
где:
CMin является меньшей из скоростей теплоемкости двух жидкостей. Скорость теплопроизводительности является произведением удельного тепла текучей среды, cp и массового расхода текучей среды. CMin всегда положительный.
TIn, TL - температура на входе термической жидкости.
TIn, МА - температура на входе влажного воздуха.
λ - эффективность теплообменника.
Эффективность является функцией скорости теплопроизводительности и количества узлов передачи, NTU, а также изменяется в зависимости от расположения потока теплообменника, что более подробно обсуждается в разделе Эффективность по устройству потока. NTU рассчитывается как:
1CMinR,
где R - общее тепловое сопротивление между двумя потоками вследствие конвекции, проводимости и любого загрязнения стенок трубы:
+ 1UMAATh, MA,
и где:
U - коэффициент конвективной теплопередачи соответствующей текучей среды. Этот коэффициент более подробно обсуждается в разделе Коэффициенты теплопередачи.
F - коэффициент обрастания на термической жидкой или влажной воздушной стороне, соответственно.
RW - тепловое сопротивление через теплопередающую поверхность.
ATh - площадь поверхности теплопередачи соответствующей стороны теплообменника. ATh - сумма площади поверхности стенки, AW, и общей площади поверхности ребра, AF:
startFAF,
где «F» - эффективность Фина.
Эффективность теплообменника изменяется в зависимости от его конфигурации потока и смешивания в каждой текучей среде.
Если для параметра «Расположение потока» установлено значение Parallel flow:
1 + CR
Если для параметра «Расположение потока» установлено значение Counter flow:
NTU (1 − CR)]
Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено в Both fluids unmixed:
1]}
Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено в Both fluids mixed:
− 1NTU] − 1
Когда одна текучая среда смешивается, а другая не смешивается, уравнение эффективности зависит от относительной теплоемкости текучих сред. Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено на Thermal Liquid 1 mixed & Moist Air 2 unmixed или Thermal Liquid 1 unmixed & Moist Air 2 mixed:
Когда жидкость с Cmax смешивается и жидкость с Cmin размешивается:
(− NTU)}})
Когда жидкость с Cmin смешивается и жидкость с Cmax не смешивается:
)]}
CR обозначает отношение между скоростями теплоемкости двух текучих сред:
CMinCMax.
На стороне влажного воздуха на поверхности теплопередачи может образовываться слой конденсации. Этот жидкий слой может влиять на количество тепла, передаваемого между влажным воздухом и термической жидкостью. Уравнения теплопередачи E-NTU приведены выше для сухой теплопередачи. Для коррекции влияния конденсации дополнительно рассчитываются уравнения E-NTU с приведенными ниже влажными параметрами. Какой из двух расчетных расходов тепла приводит к большему количеству охлаждения стороны влажного воздуха, используется в тепловых расчетах для каждой зоны [1]. Для использования этого метода предполагается, что число Льюиса близко к 1 [1], что верно для влажного воздуха.
Количество E-NTU, используемое для расчета скорости теплопередачи
| Сухой расчет | Влажный расчет | |
|---|---|---|
| Температура на входе в зону влажного воздуха | Олово, Массачусетс | Олово, wb, MA |
| Расход тепловой мощности | ||
| Коэффициент теплопередачи | УМА | p, MA |
где:
Олово, МА - температура влажного воздуха на входе.
Олово, wb, MA - влажная температура влажного воздуха, связанная с оловом, MA.
- массовый расход сухого воздуха.
MA - теплоемкость влажного воздуха на единицу массы сухого воздуха.
MA - эквивалентная теплоемкость. Эквивалентная теплоемкость представляет собой изменение удельной энтальпии влажного воздуха (на единицу сухого воздухаMA, по отношению к температуре в условиях насыщенного влажного воздуха:
∂h¯MA∂TMA) s.
Массовый расход сконденсированного водяного пара, выходящего из массового потока влажного воздуха, зависит от относительной влажности между впускным отверстием для влажного воздуха и стенкой канала и NTU теплообменника:
− NTUMA),
где:
Wwall, MA - отношение влажности на поверхности теплопередачи.
Win, MA - отношение влажности на входе потока влажного воздуха.
NTUMA - количество передающих устройств на стороне влажного воздуха, рассчитанное как:
Энергетический поток, связанный с конденсацией водяного пара, основан на разнице между удельной энтальпией пара, hwater, wall и удельной энтальпией испарения, hfg, для воды:
hfg).
Предполагается, что конденсат не накапливается на поверхности теплопередачи и не влияет на геометрические параметры, такие как диаметр трубы.
Приведенные ниже уравнения применимы как к термической жидкости, так и к влажному воздуху и используют соответствующие свойства жидкости.
Коэффициент конвективной теплопередачи изменяется в соответствии с числом Нусельта жидкости:
NukDH,
где:
Nu - среднее число Нуссельта, которое зависит от режима потока.
k - теплопроводность жидкости.
DH - гидравлический диаметр трубы.
Для турбулентных потоков число Нуссельта рассчитывается с корреляцией Гниелинского:
(Pr2/3 − 1),
где:
Re - жидкостное число Рейнольдса.
Pr - число Прандтля жидкости.
Для ламинарных потоков число Нуссельта задается параметром числа Нуссельта потока Ламинара.
Для переходных потоков число Нуссельта представляет собой смешение между ламинарными и турбулентными числами Нуссельта.
Если для геометрии потока задано значение Flow perpendicular to bank of circular tubes, число Нуссельта рассчитывается на основе числа Хагена, Hg, и зависит от настройки расположения сетки банка Тюбик:
, в шахматном порядке
где:
, в шахматном порядке с lL < D
D - наружный диаметр трубы.
1L - продольный шаг трубы (вдоль направления потока), расстояние между центрами трубы вдоль направления потока. Направление потока - направление потока внешней текучей среды.
lT - поперечный шаг трубы (перпендикулярный направлению потока), расстояние между центрами трубы в одном ряду другой жидкости.
lD - расстояние между диагональными трубами, рассчитанное как lL2.
Для получения дополнительной информации о вычислении числа Хагена см. [6].
Измерения lL и lT показаны в поперечном сечении пучка труб ниже. Эти расстояния одинаковы для обоих типов расположения банка сетки.
Поперечное сечение труб с измерениями шага


Если для параметра модели коэффициента теплопередачи установлено значение Colburn equation или если для геометрии потока задано значение Generic, число Нуссельта вычисляется эмпирическим уравнением Колберна:
aRebPrc,
где a, b и c определены в коэффициентах [a, b, c] для параметра * Re ^ b * Pr ^ c.
Приведенные ниже уравнения применимы как к термической жидкости, так и к влажному воздуху и используют соответствующие свойства жидкости.
Потеря давления из-за вязкого трения изменяется в зависимости от режима потока и конфигурации.
Для турбулентных потоков, когда число Рейнольдса выше нижнего предела числа Рейнольдса турбулентного потока, и когда для модели потери давления установлено значение Correlations for flow inside tubes, потеря давления из-за трения рассчитывается в терминах коэффициента трения Дарси.
Для стороны термической жидкости разность давлений между портом A1 и внутренним узлом I1 составляет:
LAdd2),
где:
m˙A1 - общий расход через порт A1.
fD, A - коэффициент трения Дарси, согласно корреляции Хааланда :
1,11]} -2,
где αR - абсолютная шероховатость внутренней поверхности трубы для термической жидкости. Следует отметить, что коэффициент трения зависит от числа Рейнольдса и рассчитывается в обоих окнах для каждой жидкости.
L - общая длина каждой трубки на стороне термической жидкости.
LAdd - тепловая жидкостная сторона Агрегатная эквивалентная длина локальных сопротивлений, которая представляет собой эквивалентную длину трубки, которая вносит ту же величину потерь, что и сумма потерь из-за других локальных сопротивлений в трубке.
ACS - общая площадь поперечного сечения трубы.
Перепад давления между портом B1 и внутренним узлом I1:
LAdd2),
где m˙B1 - общий расход через порт B1.
Коэффициент трения Дарси в порту B1:
1,11]} -2.
Для ламинарных потоков, когда число Рейнольдса ниже верхнего предела числа Рейнольдса потока Ламинара, и когда для модели потери давления установлено значение Correlations for flow inside tubes, потеря давления из-за трения рассчитывается в терминах постоянной трения Ламинара для коэффициента трения Дарси, λ. λ - определяемый пользователем параметр, если для параметра «Сечение трубы» установлено значение Genericв противном случае значение вычисляется внутренне.
Перепад давления между портом A1 и внутренним узлом I1:
LAdd2),
где λ - динамическая вязкость жидкости. Перепад давления между портом B1 и внутренним узлом I1:
LAdd2).
Для переходных потоков, если для модели потери давления установлено значение Correlations for flow inside tubesперепад давления, обусловленный вязким трением, представляет собой сглаженную смесь между значениями ламинарных и турбулентных потерь давления.
Если для модели потери давления установлено значение Pressure loss coefficient или если для геометрии потока задано значение Generic, потери давления из-за вязкого трения вычисляются с эмпирическим коэффициентом потери давления, Одни и те же уравнения применяются к сторонам влажного воздуха и термической жидкости и используют соответствующие свойства жидкости.
Для стороны термической жидкости разность давлений между портом A1 и внутренним узлом I1 составляет:
Перепад давления между портом B1 и внутренним узлом I1:
Если для геометрии потока задано значение Flow perpendicular to bank of circular tubes, число Хагена используется для расчета потери давления из-за вязкого трения. Одни и те же уравнения применяются к сторонам влажного воздуха и термической жидкости и используют соответствующие свойства жидкости.
Для стороны влажного воздуха перепад давления между портом A2 и внутренним узлом I2 составляет:
(Re),
где:
λ - динамическая вязкость влажной воздушной жидкости.
NR - количество рядов труб вдоль направления потока. Когда влажный воздух течет наружу к пучку труб, это количество рядов труб для термической жидкости вдоль направления потока влажного воздуха.
Перепад давления между портом B2 и внутренним узлом I2:
(Re).
Если для модели потери давления задано значение Euler number per tube row или если для геометрии потока задано значение Generic, потеря давления из-за вязкого трения рассчитывается с коэффициентом потери давления, в терминах числа Эйлера, Eu:
где λ - эмпирический коэффициент потери давления.
Перепад давления между портом A2 и внутренним узлом I2:
Перепад давления между портом B2 и внутренним узлом I2:
Суммарная скорость накопления массы в термической жидкости определяется как:
где:
MTL - общая масса термической жидкости.
m˙A1 - массовый расход жидкости в канале A1.
m˙B1 - массовый расход жидкости в канале B1.
Поток является положительным при протекании в блок через порт.
Уравнение энергосбережения связывает изменение удельной внутренней энергии с теплопередачей текучей средой:
ϕB1 − Q,
где:
uTL представляет собой тепловую жидкую специфическую внутреннюю энергию.
φA1 - расход энергии в порту A1.
φB1 - расход энергии в порту B1.
Q - скорость теплопередачи, которая положительна при выходе из термического объема жидкости.
Существует три уравнения сохранения массы на стороне влажного воздуха: одно для влажной воздушной смеси, одно для конденсированного водяного пара и одно для следового газа.
Примечание
Если для модели Trace gas установлено значение None в блоке свойств влажного воздуха (МА) следовой газ не моделируется блоками в сети влажного воздуха. В блоке теплообменника (TL-MA) это означает, что уравнение сохранения для следового газа установлено в 0.
Массовая скорость накопления влажной воздушной смеси учитывает изменения массового расхода всего влажного воздуха через отверстия теплообменника и массового расхода конденсации:
Уравнение сохранения массы для водяного пара учитывает прохождение водяного пара через сторону влажного воздуха и образование конденсации:
=m˙w,A2+m˙w,B2−m˙Cond,
где:
xw - массовая доля пара. - скорость изменения этой фракции.
- массовый расход водяного пара в порту A2.
- массовый расход водяного пара в порту B2.
- скорость конденсации.
Массовый баланс следового газа составляет:
=m˙g,A2+m˙g,B2,
где:
xg - массовая доля следового газа. - скорость изменения этой фракции.
- массовый расход следового газа в порту A2.
- массовый расход следового газа в порту B2.
Экономия энергии на стороне влажного воздуха объясняет изменение удельной внутренней энергии вследствие теплопередачи и конденсации водяного пара из массы влажного воздуха:
Q −
где:
ϕA2 - расход энергии в порту A2.
ϕB2 - расход энергии в порту B2.
β Cond - расход энергии, обусловленный конденсацией.
Теплота, передаваемая влажному воздуху или от него, Q, равна теплоте, передаваемой от термической жидкости или к ней.
[1] Справочник ASHRAE 2013 - Основы. Американское общество инженеров отопления, охлаждения и кондиционирования воздуха, Inc., 2013.
[2] Браун, Дж. Э., С. А. Кляйн и Дж. У. Митчелл. «Модели эффективности для градирен и охлаждающих змеевиков». Операции ASHRAE 95, № 2, (июнь 1989 года): 164-174.
[3] Ченгель, Юнус А. Тепло и массообмен: практический подход. 3-е изд., Макгро-Хилл, 2007.
[4] Дин, Х., Эппе Дж. П., Lebrun, J., Wasacz, M. "Модель охлаждающего змеевика для использования в переходных и/или влажных режимах. Теоретический анализ и экспериментальная валидация ". Материалы третьей Международной конференции по моделированию систем в зданиях (1990 год): 405-411.
[5] Митчелл, Джон У. и Джеймс Э. Браун. Принципы отопления, вентиляции и кондиционирования воздуха в зданиях. Уайли, 2013.
[6] Шах, Р. К., и Душан П. Секулич. Основы проектирования теплообменников. Джон Уайли и сыновья, 2003.
[7] Уайт, Фрэнк М. Механика жидкости. 6-е изд., Макгро-Хилл, 2009.
Испаритель конденсатора (2P-MA) | Испаритель конденсатора (TL-2P) | Теплопередача E-NTU | Теплообменник (G-TL) | Теплообменник (TL-TL)