Моделирует теплообмен между тепловой гидравлической сетью и сетью, которая может подвергаться изменению фазы
Simscape/Жидкости/Интерфейсы Гидравлической сети/Теплообменники
Блок Condenser Evaporator (TL-MA) моделирует теплообменник с одной тепловой гидравлической сетью, которая течет между портами A1 и B1, и одной двухфазной гидросистемой, которая течет между портами A2 и B2. Теплообменник может действовать как конденсатор или как испаритель. Потоки жидкости могут быть выровнены в параллельных встречных или поперечных строениях потока.
Пример теплообменника для холодильного оборудования
Можно смоделировать сторону тепловой жидкости как поток в трубах, поток вокруг двухфазной гидравлической трубки или путем эмпирической, типовой параметризации.
Блок использует метод 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, чтобы указать, разделены ли двухфазные потоки жидкости или тепловой жидкости на несколько путей перегородками или стенками. Без этих разделений поток может свободно перемешиваться и рассматривается как смешанный. Обе жидкости, одна жидкость или ни одна жидкость не могут быть смешаны в поперечном потоке. Смешивание гомогенизирует температуру жидкости вдоль направления потока второй жидкости и изменяется перпендикулярно второму потоку жидкости.
Несмешанные потоки изменяются по температуре как вдоль, так и перпендикулярно пути потока второй жидкости.
Примеры строений кросс-потока
Обратите внимание, что направление потока во время симуляции не влияет на выбранную настройку расположения потока. Порты на блоке не отражают физические положения портов в системе физического теплообмена.
Все устройства для потока являются однопроходными, что означает, что жидкости не делают нескольких поворотов в теплообменнике для дополнительных точек теплопередачи. Чтобы смоделировать многоходовой теплообменник, можно расположить несколько блоков испарителя конденсатора (TL-2P) последовательно или параллельно.
Для примера, чтобы достичь двухпроходного строения со стороны двухфазной жидкости и однопроходного строения со стороны тепловой жидкости, можно соединить стороны двухфазной жидкости последовательно и стороны тепловой жидкости с тем же входом параллельно (такие как два блока Источника Массового Расхода с половиной общего массового расхода), как показано ниже.
Устройство потока тепловой жидкости 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, TL является температурой входного отверстия зоны тепловой жидкости.
ε - эффективность теплообменника.
Эффективность является функцией скорости теплоемкости и количества передаточных модулей, 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 установлено на Thermal Liquid 1 mixed & Two-Phase Fluid 2 unmixed
или Thermal Liquid 1 unmixed & Two-Phase Fluid 2 mixed
:
Когда жидкость с Cmax смешивается и жидкость с Cmin не смешивается:
Когда жидкость с Cmin смешивается и жидкость с Cmax не смешивается:
C R обозначает отношение между скоростями теплоемкости двух жидкостей:
Коэффициент конвективной теплопередачи изменяется в соответствии с числом Нуссельта жидкости:
где:
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, падение давления на трение вычисляется в терминах коэффициента трения Дарси. Перепад давления между портом A2 и внутренним узлом I2:
где:
A2 - общая скорость потока жидкости через порт A2.
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 является общей площадью поперечного сечения трубы.
Перепад давления между портом B2 и внутренним узлом I2:
где B2 - общая скорость потока жидкости через порт B2.
Коэффициент трения Дарси в B2 порта:
Для ламинарных течений, когда число Рейнольдса ниже Laminar flow upper Reynolds number limit, падение давления на трение вычисляется в терминах Laminar friction constant for Darcy friction factor, λ. λ является пользовательским параметром, когда Tube cross-section задано значение Generic
в противном случае значение вычисляется внутренне. Перепад давления между портом A2 и внутренним узлом I2:
where - двухфазная динамическая вязкость жидкости. Перепад давления между портом B2 и внутренним узлом I2:
Для переходных потоков перепадом давления из-за вязкого трения является сглаженная смесь между значениями для ламинарного и турбулентного падения давления.
Когда Pressure loss model установлено на Pressure loss coefficient
, падения давления из-за вязкого трения вычисляются эмпирическим коэффициентом падения давления, ξ.
Перепад давления между портом A2 и внутренним узлом I2:
Перепад давления между портом B2 и внутренним узлом I2:
Когда 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 - диагональный интервал между трубами, вычисленный как
Для получения дополнительной информации о вычислении числа Хагена см. [3].
Продольное и поперечное расстояния тангажа одинаковы для обоих типов расположения решетки.
Поперечное сечение двухфазной гидравлической трубы с измерениями тангажа
Когда для 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
Число Хагена используется для вычисления падения давления из-за вязкого трения. Перепад давления между портом A1 и внутренним узлом I1:
где:
μ TL является динамической вязкостью жидкости.
N R является Number of tube rows along flow direction. Это количество двухфазной гидравлической трубки, строк вдоль направления потока тепловой жидкости.
Перепад давления между портом B1 и внутренним узлом I1:
Когда для Pressure loss model задано значение Euler number per tube row
или когда Flow geometry установлено на Generic
, потеря давления из-за вязкого трения вычисляется коэффициентом падения давления, в терминах числа Эйлера, Eu:
где ξ - эмпирический коэффициент падения давления.
Перепад давления между портом A1 и внутренним узлом I1:
Перепад давления между портом B1 и внутренним узлом I1:
Общая скорость накопления массы в двухфазной жидкости определяется как:
где:
M2P - общая масса двухфазной жидкости.
A2 - массовый расход жидкости в порту A2.
B2 - массовый расход жидкости в порту B2.
Поток положителен, когда течет в блок через порт.
Уравнение сохранения энергии связывает изменение удельной внутренней энергии с теплопередачей жидкостью:
где:
u 2P является двухфазной жидкостной удельной внутренней энергией.
φ A2 является энергетической скоростью потока жидкости в порту A2.
φ B2 является энергетической скоростью потока жидкости в порту B2.
Q - скорость теплопередачи, которая положительна при выходе из двухфазного объема жидкости.
Общая скорость накопления массы в тепловой жидкости определяется как:
Уравнение энергосбережения:
где:
ϕ A1 является энергетической скоростью потока жидкости в порту A1.
ϕ B1 является энергетической скоростью потока жидкости в порту B1.
Тепло, переданное тепловой жидкости или от нее, Q, равно теплу, переданному от или в двухфазную жидкость.
[1] Справочник ASHRAE 2013 - Основные принципы. Американское общество инженеров отопления, охлаждения и кондиционирования воздуха, Inc., 2013.
[2] Ченгель, Юнус А. Теплопередача: практический подход. 3rd ed, McGraw-Hill, 2007.
[3] Шах, Р. К., и Душан П. Секулич. Основные принципы проекта теплообменника. John Wiley & Sons, 2003.
[4] Уайт, Фрэнк М. Механика Жидкости. 6-е изд, McGraw-Hill, 2009.
Condenser Evaporator (2P-MA) | E-NTU Heat Transfer | Heat Exchanger (G-TL) | Heat Exchanger (TL-MA) | Heat Exchanger (TL-TL) | Thermostatic Expansion Valve (2P)