Моделирует теплообмен между сырой воздушной сетью и сетью, которая может подвергаться изменению фазы
Simscape/Жидкости/Интерфейсы Гидравлической сети/Теплообменники
Блок Condenser Evaporator (2P-MA) моделирует теплообменник с одной сырой воздушной сетью, которая течет между портами A2 и B2, и одной двухфазной гидросистемой, которая течет между портами A1 и B1. Теплообменник может действовать как конденсатор или как испаритель. Потоки жидкости могут быть выровнены в параллельных встречных или поперечных строениях потока.
Пример теплообменника для холодильного оборудования
Можно смоделировать сторону сырого воздуха как поток в трубах, поток вокруг двухфазной гидравлической трубки или с помощью эмпирической, типовой параметризации. Сырой воздух содержит воздух, следящий газ и водяной пар, которые могут конденсироваться в течение всего цикла теплообмена. Модель блока учитывает передачу энергии от воздуха к слою конденсации жидкой воды. Этот жидкий слой не собирается на поверхности теплопередачи и, как принято, полностью удаляется из последующего потока сырого воздуха. Скорость конденсации влаги возвращается как физический сигнал в порт W.
Блок использует метод Effectiveness-NTU (E-NTU), чтобы смоделировать теплопередачу через общую стенку. Загрязнение на стенках теплообменника, которое увеличивает тепловое сопротивление и уменьшает теплообмен между двумя жидкостями, также моделируется. Вы также можете опционально смоделировать ребра как на стороне сырого воздуха, так и на стороне двухфазной жидкости. Падение давления из-за вязкого трения с обеих сторон теплообменника может быть смоделировано аналитически или путем типовой параметризации, которую можно использовать, чтобы настроить на собственные данные.
Можно смоделировать двухфазную сторону жидкости как поток внутри трубы или набора труб. Двухфазные гидравлические трубы используют следующую за ними модель, чтобы отслеживать подохлаждаемую жидкость (L), смесь пар-жидкость (M) и перегретый пар (V) в трех зонах. Относительное количество пространства, которое зона занимает в системе, называется дробью длины зоны в системе.
Фракции длины зоны в двухфазном гидравлическом трубопроводе
Сумма долей длины зоны в двухфазной гидравлической трубке равна 1
. Порт Z возвращает дроби длины зоны как вектор физических сигналов для каждой из трех фаз: [L, M, V].
Эффективность теплообменника основана на выбранных строениях теплообменника, свойствах жидкости в каждой фазе, конфигурации трубы и строения потока на каждой стороне теплообменника, а также использовании и размере ребер.
Параметр Flow arrangement присваивает относительные пути потока между двумя сторонами:
Parallel flow
указывает, что жидкости движутся в том же направлении.
Counter flow
указывает, что жидкости движутся параллельно, но в противоположных направлениях.
Cross flow
указывает, что жидкости движутся перпендикулярно друг другу.
Когда Flow arrangement установлено на Cross flow
Используйте параметр Cross flow arrangement, чтобы указать, разделены ли двухфазные потоки жидкости или сырого воздуха на несколько путей перегородками или стенками. Без этих разделений поток может свободно перемешиваться и рассматривается как смешанный. Обе жидкости, одна жидкость или ни одна жидкость не могут быть смешаны в поперечном потоке. Смешивание гомогенизирует температуру жидкости вдоль направления потока второй жидкости и изменяется перпендикулярно второму потоку жидкости.
Несмешанные потоки изменяются по температуре как вдоль, так и перпендикулярно пути потока второй жидкости.
Примеры строений кросс-потока
Обратите внимание, что направление потока во время симуляции не влияет на выбранную настройку расположения потока. Порты на блоке не отражают физические положения портов в системе физического теплообмена.
Все устройства для потока являются однопроходными, что означает, что жидкости не делают нескольких поворотов в теплообменнике для дополнительных точек теплопередачи. Чтобы смоделировать многоходовой теплообменник, можно расположить несколько блоков испарителя конденсатора (2P-MA) последовательно или параллельно.
Например, чтобы достичь двухпроходного строения на двухфазной стороне жидкости и однопроходного строения на влажной стороне воздуха, можно соединить двухфазные стороны жидкости последовательно и влажные стороны воздуха с тем же входом параллельно (такие как два блока Источника Массового Расхода с половиной общего массового расхода), как показано ниже.
Устройство Flow geometry наборов параметров потоком сырого воздуха либо внутри трубки, либо внутри набора трубок, либо перпендикулярно блоку трубок. Можно также задать эмпирическое типовое строение. Двухфазная жидкость всегда течет внутрь трубы или набора труб.
Когда Flow geometry установлено на Flow perpendicular to bank of circular tubes
, используйте параметр Tube bank grid arrangement, чтобы задать выравнивание двухфазной гидравлической трубки как Inline
или Staggered
. Красный, нисходящая стрела указывает направление потока сырого воздуха. Также на Inline рисунке указаны параметры Number of tube rows along flow direction и Number of tube segments in each tube row. Здесь направление потока относится к потоку сырого воздуха, а труба относится к двухфазной гидравлической трубе. Параметр Length of each tube segment in a tube row показан на Шаггерном рисунке.
Строение теплообменника без ребер, когда параметр Total fin surface area установлен в 0 m^2
. Ребра вводят дополнительную площадь поверхности для дополнительного теплопередачи. Каждая сторона жидкости имеет отдельную площадь плавника.
Скорость теплопередачи вычисляется для каждой жидкой фазы. В соответствии с тремя зонами жидкости, которые происходят на двухфазной стороне жидкости теплообменника, скорость теплопередачи вычисляется в трех секциях.
Теплопередача в зоне вычисляется как:
где:
C Min является меньшим из скоростей теплоемкости двух жидкостей в этой зоне. Скорость теплоемкости является продуктом удельного тепла жидкости, c p и массового расхода жидкости. C Min всегда положительный.
T In, 2P является температурой на входе в зону двухфазной жидкости.
T In, MA является температурой на входе в зону сырого воздуха.
ε - эффективность теплообменника.
Эффективность является функцией скорости теплоемкости и количества передаточных модулей, NTU, а также изменяется на основе расположения потока теплообменника, которое более подробно обсуждается в Effectiveness by Flow Organization. Значение NTU определяется как:
где:
z является отдельной долей длины зоны.
R - общее тепловое сопротивление между двумя потоками из-за конвекции, проводимости и любого загрязнения стенок трубы:
где:
U - коэффициент конвективной теплопередачи соответствующей жидкости. Этот коэффициент обсуждается более подробно в Двухфазной Корреляции Жидкости и Влажных Корреляциях Воздуха.
F - Fouling factor на двухфазной жидкой или сырой воздушной стороне, соответственно.
R W является Thermal resistance through heat transfer surface.
A Th - площадь поверхности теплопередачи соответствующей стороны теплообменника. A Th - это сумма площади поверхности стенки, A W и Total fin surface area, A F:
где η F является Fin efficiency.
Общая скорость теплопередачи между жидкостями является суммой тепла, передаваемого в трех зонах переохлаждаемой жидкостью (QL), смесью пар жидкости (QM) и перегретым паром (QV):
Эффективность теплообменника изменяется в зависимости от его строения потока и смешивания в каждой жидкости. Ниже приведены составы для эффективности, рассчитанные в жидкой и паровой зонах для каждого строения. Эффективность для всех строений в зоне смеси.
Когда Flow arrangement установлено на Parallel flow
:
Когда Flow arrangement установлено на Counter flow
:
Когда Flow arrangement установлено на Cross flow
и Cross flow arrangement установлено на Both fluids unmixed
:
Когда Flow arrangement установлено на Cross flow
и Cross flow arrangement установлено на Both fluids mixed
:
Когда одна жидкость смешивается, а другая не смешивается, уравнение эффективности зависит от относительных скоростей теплоемкости жидкостей. Когда Flow arrangement установлено на Cross flow
и Cross flow arrangement установлено на Two-Phase Fluid 1 mixed & Moist Air 2 unmixed
или Two-Phase Fluid 1 unmixed & Moist Air 2 mixed
:
Когда жидкость с Cmax смешивается и жидкость с Cmin не смешивается:
Когда жидкость с Cmin смешивается и жидкость с Cmax не смешивается:
C R обозначает отношение между скоростями теплоемкости двух жидкостей:
Со стороны сырого воздуха на поверхности теплопередачи может образоваться слой конденсации. Этот жидкий слой может влиять на количество тепла, передаваемого между сырым воздухом и двухфазной жидкостью. Уравнения для теплопередачи E-NTU выше приведены для сухой теплопередачи. Для коррекции влияния конденсации уравнения E-NTU дополнительно вычисляются с мокрыми параметрами, перечисленными ниже. Любое из двух расчетных тепловых скоростей потока жидкости, результатов в большем количестве охлаждения со стороны сырого воздуха, используется в тепловых расчетах для каждой зоны [1]. Чтобы использовать этот метод, число Льюиса принято близким к 1 [1], что верно для сырого воздуха.
Количества E-NTU, используемые для вычислений скорости теплопередачи
Расчет сухого раствора | Мокрый расчет | |
---|---|---|
Температура на входе в зону сырого воздуха | Tin,MA | Tin,wb,MA |
Скорость теплоемкости | ||
Коэффициент теплопередачи | UMA |
где:
Tin,MA - температура входного отверстия зоны сырого воздуха.
Tin,wb,MA - температура влажного воздуха, смоченного с луковицей, связанная с Tin,MA.
- сухой воздушный массовый расход жидкости.
- теплоемкость сырого воздуха на единицу массы сухого воздуха.
- эквивалентная теплоемкость. Эквивалентной теплоемкостью является изменение удельной энтальпии сырого воздуха (на модуль сухого воздуха). относительно температуры в условиях насыщенного сырого воздуха:
Массовый расход жидкости конденсированного водяного пара, выходящего из потока сырой воздушной массы, зависит от относительной влажности между входным отверстием сырого воздуха и стенкой канала и NTU теплообменника:
где:
W стенке МА - коэффициент влажности на поверхности теплопередачи.
W в, МА является отношением влажности на входном отверстии потока сырого воздуха.
NTU MA - количество передаточных модулей на стороне сырого воздуха, рассчитанное как:
Энергетический поток, связанный с конденсацией паров воды, основан на различии между специфической энтальпией паров, h водой, стенкой и специфической энтальпией испарения, h fg, для воды:
Конденсат не накапливается на поверхности теплопередачи и не влияет на геометрические параметры, такие как диаметр трубы. Сконденсированная вода полностью удаляется из последующего потока сырого воздуха.
Коэффициент конвективной теплопередачи изменяется в соответствии с числом Нуссельта жидкости:
где:
Nu - среднезонное число Нуссельта, которое зависит от режима течения.
k - теплопроводность жидкой фазы.
D H - гидравлический диаметр трубы.
Для турбулентных потоков в переохлажденных жидких или перегретых паровых зонах число Нуссельта вычисляется корреляцией Гнилинского:
где:
Re - жидкое число Рейнольдса.
Pr - жидкое число Прандтля.
Для турбулентных потоков в зоне смеси пар-жидкости число Нуссельта вычисляется корреляцией Каваллини-Цекхина:
где:
Re SL является числом Рейнольдса насыщенной жидкости.
Pr SL является числом Прандтля насыщенной жидкости.
ρ SL является плотностью насыщенной жидкости.
ρ SV является плотностью насыщенного пара.
a = 0,05, b = 0,8 и c = 0,33.
Для ламинарных течений число Нуссельта задается параметром Laminar flow Nusselt number.
Для переходных потоков число Нуссельта является смесью между ламинарным и турбулентным числами Нуссельта.
Когда параметр Heat transfer coefficient model установлен в Colburn equation
число Нуссельта для переохлаждаемой жидкости и перегретых паровых зон вычисляется эмпирическим уравнением Колбёрна:
где a, b и c заданы в параметрах Coefficients [a, b, c] for a*Re^b*Pr^c in liquid zone и Coefficients [a, b, c] for a*Re^b*Pr^c in vapor zone.
Число Нуссельта для зон смеси пар-жидкости вычисляется уравнением Каваллини-Цекчина с переменными, заданными в параметре Coefficients [a, b, c] for a*Re^b*Pr^c in mixture zone.
Падение давления из-за вязкого трения изменяется в зависимости от режима течения и строения. В вычислении используется общая плотность, которая является общей двухфазной массой жидкости, разделенной на общий двухфазный объем жидкости.
Для турбулентных потоков, когда число Рейнольдса выше Turbulent flow lower Reynolds number limit, падение давления на трение вычисляется в терминах коэффициента трения Дарси. Перепад давления между портом A1 и внутренним узлом I1:
где:
A1 - общая скорость потока жидкости через порт A1.
f D, A является коэффициентом трения Дарси, согласно корреляции Haaland:
где ε R - двухфазный гидравлический трубопровод Internal surface absolute roughness. Обратите внимание, что коэффициент трения зависит от числа Рейнольдса и вычисляется в обоих портах для каждой жидкости.
L является Total length of each tube на двухфазной стороне жидкости.
L Add является двухфазной жидкостной боковой Aggregate equivalent length of local resistances, которая является эквивалентной длиной трубки, которая вводит такое же количество потерь, как и сумма потерь из-за других локальных сопротивлений в трубке.
A CS является площадью поперечного сечения трубы.
Перепад давления между портом B1 и внутренним узлом I1:
где B1 - общая скорость потока жидкости через порт B1.
Коэффициент трения Дарси в B1 порта:
Для ламинарных течений, когда число Рейнольдса ниже Laminar flow upper Reynolds number limit, падение давления на трение вычисляется в терминах Laminar friction constant for Darcy friction factor, λ. λ является пользовательским параметром, когда Tube cross-section задано значение Generic
в противном случае значение вычисляется внутренне. Перепад давления между портом A1 и внутренним узлом I1:
where - двухфазная динамическая вязкость жидкости. Перепад давления между портом B1 и внутренним узлом I1:
Для переходных потоков перепадом давления из-за вязкого трения является сглаженная смесь между значениями для ламинарного и турбулентного падения давления.
Когда Pressure loss model установлено на Pressure loss coefficient
, падения давления из-за вязкого трения вычисляются эмпирическим коэффициентом падения давления, ξ.
Перепад давления между портом A1 и внутренним узлом I1:
Перепад давления между портом B1 и внутренним узлом I1:
Когда значение Flow geometry сырого воздуха установлено равным Flow inside one or more tubes
число Нуссельта вычисляют согласно корреляции Гнилинского таким же образом, как двухфазная переохлажденная жидкость или перегретый пар. Смотрите Коэффициент Теплопередачи для получения дополнительной информации.
Когда значение Flow geometry сырого воздуха установлено равным Flow perpendicular to bank of circular tubes
число Нуссельта вычисляется на основе числа Хагена, Hg, и зависит от Tube bank grid arrangement настройки:
где:
D является Tube outer diameter.
l L - Longitudinal tube pitch (along flow direction), расстояние между центрами труб вдоль направления потока. Направление потока относится к потоку сырого воздуха.
l T является Transverse tube pitch (perpendicular to flow direction), показанным на рисунке ниже. Поперечный тангаж является расстоянием между центрами двухфазной гидравлической трубки в одном ряду.
l D - диагональный интервал между трубами, вычисленный как
Для получения дополнительной информации о вычислении числа Хагена см. [6].
Продольное и поперечное расстояния тангажа одинаковы для обоих типов расположения решетки.
Поперечное сечение двухфазной гидравлической трубы с измерениями тангажа
Когда для Heat transfer coefficient model задано значение Colburn equation
или когда Flow geometry установлено на Generic
число Нуссельта вычисляется эмпирическим уравнением Колберна:
где a, b и c являются значениями, заданными в параметре Coefficients [a, b, c] for a*Re^b*Pr^c.
Когда значение Flow geometry сырого воздуха установлено равным Flow inside one or more tubes
Потеря давления вычисляется таким же образом, как и для двухфазных потоков, с соответствующим коэффициентом трения Дарси, плотностью, массовыми расходами жидкости и трубопроводом длинами стороны сырого воздуха. Для получения дополнительной информации см. раздел «Падение давления».
Когда значение Flow geometry сырого воздуха установлено равным Flow perpendicular to bank of circular tubes
Число Хагена используется для вычисления падения давления из-за вязкого трения. Перепад давления между портом A2 и внутренним узлом I2:
где:
μ MA является динамической вязкостью жидкости.
N R является Number of tube rows along flow direction. Это количество двухфазной гидравлической трубки, строк вдоль направления потока сырого воздуха.
Перепад давления между портом B2 и внутренним узлом I2:
Когда для Pressure loss model задано значение Euler number per tube row
или когда Flow geometry установлено на Generic
, потеря давления из-за вязкого трения вычисляется коэффициентом падения давления, в терминах числа Эйлера, Eu:
где ξ - эмпирический коэффициент падения давления.
Перепад давления между портом A2 и внутренним узлом I2:
Перепад давления между портом B2 и внутренним узлом I2:
Общая скорость накопления массы в двухфазной жидкости определяется как:
где:
M2P - общая масса двухфазной жидкости.
A1 - массовый расход жидкости в порту A1.
B1 - массовый расход жидкости в порту B1.
Поток положителен, когда течет в блок через порт.
Уравнение сохранения энергии связывает изменение удельной внутренней энергии с теплопередачей жидкостью:
где:
u 2P является двухфазной жидкостной удельной внутренней энергией.
φ A1 является энергетической скоростью потока жидкости в порту A1.
φ B1 является энергетической скоростью потока жидкости в порту B1.
Q - скорость теплопередачи, которая положительна при выходе из двухфазного объема жидкости.
Существует три уравнения для сохранения массы на стороне сырого воздуха: одно для смеси сырого воздуха, одно для конденсированного водяного пара и одно для следового газа.
Примечание
Если для Trace gas model задано значение None
в Moist Air Properties (MA) блоке следовой газ не моделируется блоками во влажной воздушной сети. В блоке Condenser Evaporator (2P-MA) это означает, что уравнение сохранения для следового газа установлено на 0.
Скорость накопления массы сырой воздушной смеси учитывает изменения всего потока сырой воздушной массы через порты теплообменника и массового расхода жидкости конденсации:
Уравнение сохранения массы для водяного пара учитывает транзит водяного пара через сторону сырого воздуха и образование конденсации:
где:
x w - массовая доля пара . - скорость изменения этой дроби.
- водяной пар, массовый расход жидкости в порту A2.
- водяной пар, массовый расход жидкости в порту B2.
- скорость конденсации.
Баланс массы следового газа:
где:
x g - массовая доля следового газа . - скорость изменения этой дроби.
- трассировочный газ, массовый расход жидкости в порту A2.
- трассировочный газ, массовый расход жидкости в порту B2.
Экономия энергии на стороне сырого воздуха учитывает изменение удельной внутренней энергии из-за теплопередачи и конденсации водяного пара из массы сырого воздуха:
где:
ϕ A2 является энергетической скоростью потока жидкости в порту A2.
ϕ B2 является энергетической скоростью потока жидкости в порту B2.
ϕ Cond - это расход энергии из-за конденсации.
Тепло, переданное к или от сырого воздуха, Q, равно теплу, переданному от или к двухфазной жидкости.
[1] Справочник ASHRAE 2013 - Основные принципы. Американское общество инженеров отопления, охлаждения и кондиционирования воздуха, Inc., 2013.
[2] Braun, J. E., S. A. Klein, and J. W. Mitchell. «Модели эффективности для градирен и охлаждающих катушек». Транзакции ASHRAE 95, № 2, (июнь 1989 г.): 164-174.
[3] Ченгель, Юнус А. Теплопередача: практический подход. 3rd ed, McGraw-Hill, 2007.
[4] Ding, X., Eppe J.P., Lebrun, J., Wasacz, M. "Модель охлаждающей катушки для использования в переходных и/или влажных режимах. Теоретический анализ и экспериментальная валидация ". Материалы третьей Международной конференции по системной симуляции в созданиях (1990 год): 405-411.
[5] Митчелл, Джон У. и Джеймс Э. Браун. Принципы отопления, вентиляции и кондиционирования воздуха в созданиях. Уайли, 2013.
[6] Шах, Р. К., и Душан П. Секулич. Основные принципы проекта теплообменника. John Wiley & Sons, 2003.
[7] Уайт, Фрэнк М. Механика Жидкости. 6-е изд, McGraw-Hill, 2009.
Condenser Evaporator (TL-2P) | E-NTU Heat Transfer | Heat Exchanger (G-TL) | Heat Exchanger (TL-MA) | Heat Exchanger (TL-TL) | Thermostatic Expansion Valve (2P)