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

Конденсаторный испаритель (2P-MA) моделирует теплообменник с одной сетью влажного воздуха, которая протекает между портами A2 и B2, и одной двухфазной сетью текучей среды, которая протекает между портами A1 и B1. Теплообменник может действовать как конденсатор или как испаритель. Потоки текучей среды могут быть выровнены параллельно, встречно или поперечно.
Пример теплообменника для холодильных установок

Можно моделировать сторону влажного воздуха как поток внутри труб, поток вокруг двухфазной трубки текучей среды или с помощью эмпирической общей параметризации. Влажная воздушная сторона содержит воздух, следовый газ и водяной пар, которые могут конденсироваться в течение всего цикла теплообмена. Модель блока учитывает передачу энергии из воздуха в слой конденсации жидкой воды. Этот жидкий слой не собирается на поверхности теплопередачи и, как предполагается, полностью удаляется из потока влажного воздуха, расположенного ниже по потоку. Скорость конденсации влаги возвращается как физический сигнал в порту W.
Блок использует метод Effectivity-NTU (E-NTU) для моделирования теплопередачи через общую стенку. Также моделируется загрязнение стенок теплообменника, которое увеличивает тепловое сопротивление и уменьшает теплообмен между двумя текучими средами. Можно также дополнительно моделировать ребра как на стороне влажного воздуха, так и на стороне двухфазной текучей среды. Потеря давления из-за вязкого трения с обеих сторон теплообменника может быть смоделирована аналитически или с помощью общей параметризации, которую можно использовать для настройки на собственные данные.
Можно моделировать двухфазную сторону текучей среды как поток внутри трубы или набора труб. Двухфазные трубки для текучей среды используют следующую за границей модель для отслеживания субохлаждаемой жидкости (L), парожидкостной смеси (M) и перегретого пара (V) в трех зонах. Относительная величина пространства, занимаемого зоной в системе, называется дробью длины зоны в системе.
Доли длины зоны в двухфазном трубопроводе текучей среды

Сумма долей длины зоны в двухфазных трубах флюида равна 1. Порт Z возвращает доли длины зоны как вектор физических сигналов для каждой из трех фаз: [L, M, V].
Эффективность теплообменника основана на выбранной конфигурации теплообменника, свойствах жидкости в каждой фазе, геометрии трубы и конфигурации потока на каждой стороне теплообменника, а также использовании и размере ребер.
Параметр расположения потока задает относительные пути потока между двумя сторонами:
Parallel flow указывает, что текучие среды движутся в том же направлении.
Counter flow указывает, что жидкости движутся параллельно, но в противоположных направлениях.
Cross flow указывает, что жидкости движутся перпендикулярно друг другу.
Если для параметра «Расположение потока» установлено значение Cross flowс помощью параметра «Расположение поперечных потоков» укажите, разделяются ли двухфазные потоки текучей среды или влажного воздуха на несколько путей перегородками или стенками. Без этих разделений поток может свободно смешиваться и считается смешанным. Обе жидкости, одна жидкость или ни одна из них не могут быть смешаны в устройстве поперечного потока. Перемешивание гомогенизирует температуру текучей среды вдоль направления потока второй текучей среды и изменяется перпендикулярно второму потоку текучей среды.
Несмешанные потоки изменяются по температуре как вдоль, так и перпендикулярно пути потока второй текучей среды.
Примеры кросс-поточных конфигураций



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

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


Конфигурация теплообменника не имеет ребер, если для параметра Общая площадь поверхности ребер установлено значение 0 m^2. Ребра вводят дополнительную площадь поверхности для дополнительной теплопередачи. Каждая сторона текучей среды имеет отдельную область ребер.
Скорость теплопередачи рассчитывают для каждой фазы текучей среды. В соответствии с тремя жидкостными зонами, которые возникают на двухфазной жидкостной стороне теплообменника, скорость теплопередачи вычисляется в трех секциях.
Теплопередача в зоне рассчитывается как:
MA),
где:
CMin является меньшей из скоростей теплоемкости двух текучих сред в этой зоне. Скорость теплопроизводительности является произведением удельного тепла текучей среды, cp и массового расхода текучей среды. CMin всегда положительный.
TIn,2P - температура зоны на входе двухфазной текучей среды.
TIn, МА - температура зоны на входе влажного воздуха.
λ - эффективность теплообменника.
Эффективность является функцией скорости теплопроизводительности и количества узлов передачи, NTU, а также изменяется в зависимости от расположения потока теплообменника, что более подробно обсуждается в разделе Эффективность по устройству потока. NTU рассчитывается как:
zCMinR,
где:
z - доля длины отдельной зоны.
R - общее тепловое сопротивление между двумя потоками вследствие конвекции, проводимости и любого загрязнения стенок трубы:
+ 1UMAATh, MA,
где:
U - коэффициент конвективной теплопередачи соответствующей текучей среды. Этот коэффициент более подробно обсуждается в двухфазных корреляциях текучей среды и корреляциях влажного воздуха.
F - коэффициент обрастания на двухфазной текучей среде или влажной воздушной стороне соответственно.
RW - тепловое сопротивление через теплопередающую поверхность.
ATh - площадь поверхности теплопередачи соответствующей стороны теплообменника. ATh - сумма площади поверхности стенки, AW, и общей площади поверхности ребра, AF:
startFAF,
где «F» - эффективность Фина.
Общая скорость теплопередачи между текучими средами представляет собой сумму теплоты, передаваемой в трех зонах переохлажденной жидкостью (QL), смесью жидкость-пар (QM) и перегретым паром (QV):
Эффективность теплообменника изменяется в зависимости от его конфигурации потока и смешивания в каждой текучей среде. Ниже приведены составы для определения эффективности, рассчитанные в жидкостной и паровой зонах для каждой конфигурации. Эффективность - (−NTU) для всех конфигураций в зоне смеси.
Если для параметра «Расположение потока» установлено значение Parallel flow:
1 + CR
Если для параметра «Расположение потока» установлено значение Counter flow:
NTU (1 − CR)]
Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено в Both fluids unmixed:
1]}
Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено в Both fluids mixed:
− 1NTU] − 1
Когда одна текучая среда смешивается, а другая не смешивается, уравнение эффективности зависит от относительной теплоемкости текучих сред. Если для параметра «Расположение потока» установлено значение Cross flow и расположение поперечных потоков установлено на Two-Phase Fluid 1 mixed & Moist Air 2 unmixed или Two-Phase Fluid 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 - число Прандтля жидкости.
Для турбулентных потоков в зоне парожидкостной смеси вычисляют число Нуссельта с корреляцией Каваллини - Зекчина:
b) (αSLαSV − 1) (xOut − xIn).
где:
ReSL - число Рейнольдса насыщенной жидкости.
PrSL - число Прандтля насыщенной жидкости.
αSL - плотность насыщенной жидкости.
δSV - плотность насыщенного пара.
a = 0,05, b = 0,8 и c = 0,33.
Для ламинарных потоков число Нуссельта задается параметром числа Нуссельта потока Ламинара.
Для переходных потоков число Нуссельта представляет собой смешение между ламинарными и турбулентными числами Нуссельта.
Если для параметра модели коэффициента теплопередачи установлено значение Colburn equationчисло Нусельта для зон переохлажденной жидкости и перегретого пара вычисляется эмпирическим уравнением Колберна:
aRebPrc,
где a, b и c определены в коэффициентах [a, b, c] для a * Re ^ b * Pr ^ c в зоне жидкости и коэффициентах [a, b, c] для a * Re * b * Pr ^ c в параметрах паровой зоны.
Число Нуссельта для зон смеси жидкость-пар вычисляется с помощью уравнения Каваллини-Зекчина с переменными, указанными в коэффициентах [a, b, c] для параметра * Re ^ b * Pr ^ c в зоне смеси.
Потеря давления из-за вязкого трения изменяется в зависимости от режима потока и конфигурации. В расчете используется общая плотность, которая представляет собой общую двухфазную массу текучей среды, деленную на общий двухфазный объем текучей среды.
Для турбулентных потоков, когда число Рейнольдса выше предела числа турбулентного потока ниже Рейнольдса, потери давления из-за трения вычисляются в терминах коэффициента трения Дарси. Перепад давления между портом 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.
Для ламинарных потоков, когда число Рейнольдса ниже верхнего предела числа Рейнольдса потока Ламинара, потери давления из-за трения вычисляются в терминах постоянной трения Ламинара для коэффициента трения Дарси, λ. λ - определяемый пользователем параметр, если для параметра «Сечение трубы» установлено значение Genericв противном случае значение вычисляется внутренне. Перепад давления между портом A1 и внутренним узлом I1:
LAdd2),
где λ - двухфазная динамическая вязкость текучей среды. Перепад давления между портом B1 и внутренним узлом I1:
LAdd2).
Для переходных потоков перепад давления, обусловленный вязким трением, представляет собой сглаженную смесь между значениями ламинарных и турбулентных потерь давления.
Если для модели потери давления установлено значение Pressure loss coefficient, потери давления из-за вязкого трения вычисляются с эмпирическим коэффициентом потери давления,
Перепад давления между портом A1 и внутренним узлом I1:
Перепад давления между портом B1 и внутренним узлом I1:
Если для геометрии потока влажного воздуха задано значение Flow inside one or more tubes, число Нусельта вычисляется согласно корреляции Гниелинского таким же образом, как двухфазная переохлажденная жидкость или перегретый пар. Дополнительные сведения см. в разделе Коэффициент теплопередачи.
Если для геометрии потока влажного воздуха задано значение Flow perpendicular to bank of circular tubes, число Нуссельта рассчитывается на основе числа Хагена, Hg, и зависит от настройки расположения сетки банка Тюбик:
, в шахматном порядке
где:
, в шахматном порядке с lL < D
D - наружный диаметр трубы.
1L - продольный шаг трубы (вдоль направления потока), расстояние между центрами трубы вдоль направления потока. Направление потока относится к потоку влажного воздуха.
lT - поперечный шаг трубы (перпендикулярно направлению потока), показанный на рисунке ниже. Поперечный шаг представляет собой расстояние между центрами двухфазной жидкостной трубки в одном ряду.
lD - расстояние между диагональными трубами, рассчитанное как lL2.
Для получения дополнительной информации о вычислении числа Хагена см. [6].
Продольные и поперечные расстояния наклона одинаковы для обоих типов расположения групп решеток.
Поперечное сечение двухфазной трубки текучей среды с измерениями шага


Если для модели коэффициента теплопередачи установлено значение Colburn equation или если для геометрии потока задано значение Generic, число Нуссельта вычисляется эмпирическим уравнением Колберна:
aRebPrc,
где a, b и c - значения, определенные в коэффициентах [a, b, c] для параметра * Re ^ b * Pr ^ c.
Если для геометрии потока влажного воздуха задано значение Flow inside one or more tubesпотери давления вычисляются таким же образом, как и для двухфазных потоков, с соответствующим коэффициентом трения Дарси, плотностью, массовыми скоростями потока и длинами труб стороны влажного воздуха. Для получения дополнительной информации см. раздел Потеря давления.
Если для геометрии потока влажного воздуха задано значение 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:
Суммарная скорость накопления массы в двухфазной текучей среде определяется как:
где:
M2P - общая масса двухфазной текучей среды.
m˙A1 - массовый расход жидкости в канале A1.
m˙B1 - массовый расход жидкости в канале B1.
Поток является положительным при протекании в блок через порт.
Уравнение энергосбережения связывает изменение удельной внутренней энергии с теплопередачей текучей средой:
ϕB1 − Q,
где:
u2P - двухфазная внутренняя энергия, специфическая для текучей среды.
φA1 - расход энергии в порту A1.
φB1 - расход энергии в порту B1.
Q - скорость теплопередачи, которая положительна при выходе из двухфазного объема текучей среды.
Существует три уравнения сохранения массы на стороне влажного воздуха: одно для влажной воздушной смеси, одно для конденсированного водяного пара и одно для следового газа.
Примечание
Если для модели Trace gas установлено значение None в блоке свойств влажного воздуха (МА) следовой газ не моделируется блоками в сети влажного воздуха. В блоке конденсаторного испарителя (2P-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.
Испаритель конденсатора (TL-2P) | Теплопередача E-NTU | Теплообменник (G-TL) | Теплообменник (TL-MA) | Теплообменник (TL-TL) | Клапан термостатического расширения (2P)