Моделирует теплообмен между сырой воздушной сетью и тепловой гидравлической сетью
Simscape/Жидкости/Интерфейсы Гидравлической сети/Теплообменники
Блок Heat Exchanger (TL-MA) моделирует теплообменник с одной сырой воздушной сетью, которая течет между портами A2 и B2, и одной тепловой гидравлической сетью, которая течет между портами A1 и B1. Потоки жидкости могут быть выровнены в параллельных встречных или поперечных строениях потока.
Тепловой жидкостный воздушный теплообменник не подходит для систем охлаждения. Смотрите Испаритель Конденсатора (2P-MA) или Испаритель Конденсатора (TL-2P) для теплообменников, которые могут использоваться в холодильных системах.
Можно смоделировать сторону сырого воздуха как поток в трубах, поток вокруг трубы с тепловой жидкостью или с помощью эмпирической, типовой параметризации. Сырой воздух содержит воздух, следящий газ и водяной пар, которые могут конденсироваться в течение всего цикла теплообмена. Модель блока учитывает скрытое тепло, которое выделяется, когда вода конденсируется на поверхности теплопередачи. Этот жидкий слой не собирается на поверхности и, как принято, полностью удаляется из последующего потока сырого воздуха. Скорость конденсации влаги возвращается как физический сигнал в порт W.
Блок использует метод Effectiveness-NTU (E-NTU), чтобы смоделировать теплопередачу через общую стенку. Загрязнение на стенках теплообменника, которое увеличивает тепловое сопротивление и уменьшает теплообмен между двумя жидкостями, также моделируется. Можно также опционально смоделировать ребра как на стороне сырого воздуха, так и на стороне тепловой жидкости. Падение давления из-за вязкого трения с обеих сторон теплообменника может быть смоделировано аналитически или путем типовой параметризации, которую можно использовать, чтобы настроить на собственные данные.
Можно смоделировать сторону тепловой жидкости как поток в трубах, поток вокруг трубки с сырым воздухом или с помощью эмпирической, типовой параметризации.
Эффективность теплообменника основана на выбранном строении теплообменника, свойствах жидкости, геометрии трубы и строения потока на каждой стороне теплообменника, а также использовании и размере ребер.
Параметр Flow arrangement присваивает относительные пути потока между двумя сторонами:
Parallel flow
указывает, что жидкости движутся в том же направлении.
Counter flow
указывает, что жидкости движутся параллельно, но в противоположных направлениях.
Cross flow
указывает, что жидкости движутся перпендикулярно друг другу.
Когда Flow arrangement установлено на Cross flow
используйте параметр Cross flow arrangement, чтобы указать, разделены ли потоки тепловой жидкости или сырого воздуха на несколько путей перегородками или стенками. Без этих разделений поток может свободно перемешиваться и рассматривается как смешанный. Обе жидкости, одна жидкость или ни одна жидкость не могут быть смешаны в поперечном потоке. Смешивание гомогенизирует температуру жидкости вдоль направления потока второй жидкости и изменяется перпендикулярно второму потоку жидкости.
Несмешанные потоки изменяются по температуре как вдоль, так и перпендикулярно пути потока второй жидкости.
Примеры строений кросс-потока
Обратите внимание, что направление потока во время симуляции не влияет на выбранную настройку расположения потока. Порты на блоке не отражают физические положения портов в системе физического теплообмена.
Все устройства для потока являются однопроходными, что означает, что жидкости не делают нескольких поворотов в теплообменнике для дополнительных точек теплопередачи. Чтобы смоделировать многоходовой теплообменник, можно расположить несколько блоков теплообменника (TL-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 показан на Шаггерном рисунке.
Только одна жидкость может Flow geometry установлена на Flow perpendicular to bank of circular tubes
за раз. Другая жидкость должна быть сконфигурирована таким образом Flow inside one or more tubes
или Generic
. Если Flow geometry для обеих жидкостей установлено в Flow perpendicular to bank of circular tubes
, вы получите ошибку.
Строение теплообменника без ребер, когда параметр Total fin surface area установлен в 0 m^2
. Ребра вводят дополнительную площадь поверхности для дополнительного теплопередачи. Каждая сторона жидкости имеет отдельную площадь плавника.
Скорость теплопередачи вычисляется по усредненным свойствам обеих жидкостей.
Теплопередача рассчитывается как:
где:
C Min является меньшим из скоростей теплоемкости двух жидкостей. Скорость теплоемкости является продуктом удельного тепла жидкости, c p и массового расхода жидкости. C Min всегда положительный.
T In, TL является входной температурой тепловой жидкости.
T In, MA является температурой на входе сырого воздуха.
ε - эффективность теплообменника.
Эффективность является функцией скорости теплоемкости и количества передаточных модулей, NTU, а также изменяется на основе расположения потока теплообменника, которое более подробно обсуждается в Effectiveness by Flow Organization. Значение NTU определяется как:
где 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.
Эффективность теплообменника изменяется в зависимости от его строения потока и смешивания в каждой жидкости.
Когда 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 & Moist Air 2 unmixed
или Thermal Liquid 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 - жидкое число Прандтля.
Для ламинарных течений число Нуссельта задается параметром Laminar flow Nusselt number.
Для переходных потоков число Нуссельта является смесью между ламинарным и турбулентным числами Нуссельта.
Когда 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].
Измерения l L и l T показаны в поперечном сечении трубчатого блока ниже. Эти расстояния одинаковы для обоих типов расположения банка сетки.
Поперечное сечение труб с измерениями тангажа
Когда параметр Heat transfer coefficient model установлен в Colburn equation
или когда Flow geometry установлено на Generic
число Нуссельта вычисляется эмпирическим уравнением Колберна:
где a, b и c заданы в параметре Coefficients [a, b, c] for a*Re^b*Pr^c.
Приведенные ниже уравнения применяются как к стороне тепловой жидкости, так и к стороне сырого воздуха и используют соответствующие свойства жидкости.
Падение давления из-за вязкого трения изменяется в зависимости от режима течения и строения.
Для турбулентных потоков, когда число Рейнольдса выше Turbulent flow lower Reynolds number limit, и когда Pressure loss model установлено на Correlations for flow inside tubes
, падение давления на трение вычисляется в терминах коэффициента трения Дарси.
Для стороны тепловой жидкости перепадом давления между портом 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, и когда Pressure loss model установлено на Correlations for flow inside tubes
, падение давления на трение вычисляется в терминах Laminar friction constant for Darcy friction factor, λ. λ является пользовательским параметром, когда Tube cross-section задано значение Generic
в противном случае значение вычисляется внутренне.
Перепад давления между портом A1 и внутренним узлом I1:
where - динамическая вязкость жидкости. Перепад давления между портом B1 и внутренним узлом I1:
Для переходных потоков, когда Pressure loss model установлено в Correlations for flow inside tubes
перепадом давления из-за вязкого трения является сглаженная смесь между значениями для ламинарного и турбулентного падения давления.
Когда Pressure loss model установлено на Pressure loss coefficient
или когда Flow geometry установлено на Generic
, падения давления из-за вязкого трения вычисляются эмпирическим коэффициентом падения давления, ξ. Те же уравнения применяются как к стороне сырого воздуха, так и к стороне тепловой жидкости и используют соответствующие свойства жидкости.
Для стороны тепловой жидкости перепадом давления между портом A1 и внутренним узлом I1 является:
Перепад давления между портом B1 и внутренним узлом I1:
Когда Flow geometry установлено на Flow perpendicular to bank of circular tubes
Число Хагена используется для вычисления падения давления из-за вязкого трения. Те же уравнения применяются как к стороне сырого воздуха, так и к стороне тепловой жидкости и используют соответствующие свойства жидкости.
Для стороны сырого воздуха перепадом давления между портом A2 и внутренним узлом I2 является:
где:
μ - динамическая вязкость сырой воздушной жидкости.
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:
Общая скорость накопления массы в тепловой жидкости определяется как:
где:
MTL - общая масса тепловой жидкости.
A1 - массовый расход жидкости в порту A1.
B1 - массовый расход жидкости в порту B1.
Поток положителен, когда течет в блок через порт.
Уравнение сохранения энергии связывает изменение удельной внутренней энергии с теплопередачей жидкостью:
где:
u TL является удельной внутренней энергией тепловой жидкости.
φ A1 является энергетической скоростью потока жидкости в порту A1.
φ B1 является энергетической скоростью потока жидкости в порту B1.
Q - скорость теплопередачи, которая положительна при выходе из объема тепловой жидкости.
Существует три уравнения для сохранения массы на стороне сырого воздуха: одно для смеси сырого воздуха, одно для конденсированного водяного пара и одно для следового газа.
Примечание
Если для Trace gas model задано значение None
в Moist Air Properties (MA) блоке следовой газ не моделируется блоками во влажной воздушной сети. В блоке Heat Exchanger (TL-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 (2P-MA) | Condenser Evaporator (TL-2P) | E-NTU Heat Transfer | Heat Exchanger (G-TL) | Heat Exchanger (TL-TL)