Газовая сторона теплообменника
Simscape/Жидкости/Интерфейсы Гидравлической сети/Теплообменники/Основные компоненты
Блок Heat Exchanger Interface (G) моделирует теплопередачу газообразным потоком в теплообменнике. Используйте второй блок теплообменника, чтобы смоделировать пару жидкостей. Интерфейсы могут быть в различных жидких областях, таких как один в жидкости и один в газе. Используйте блок E-NTU Heat Transfer, чтобы связать интерфейсы и захватить теплообмен между жидкостями.
Конструкция блока с фиксированным объемом позволяет вам захватывать изменения массовых расходов жидкости жидкости из-за сжимаемости. Общая скорость накопления массы равна сумме массовых расходов жидкости через порты:
где - скорость накопления массы и - массовый расход жидкости. Нижние индексы обозначают порты A и B. Массовый расход жидкости положителен, когда он направлен в газовый канал. Изменения плотности отражаются в скорости накопления массы:
где:
ρ - плотность.
p - давление.
u является специфической внутренней энергией.
V - объем.
Импульс балансировки между входным и выходным портами теплообменника определяет направление потока и скорость внутри теплообменника. Изменения импульса вызваны в основном потерями на трение от поворотов трубопровода, которые переходят к изменениям давления. Локальные сопротивления, такие как повороты, колена и тройники, могут привести к разделению потока, что приводит к незначительным дополнительным падениям давления. Для устойчивых потоков массовый расход жидкости остается постоянным.
Баланс импульса применяется к каждому сегменту объема газа (трубопровода). Этот рисунок показывает блок труб, разделенный на два объема и три узла. Узлы соответствуют портам A, B и объему жидкости I. Состояния жидкости, такие как давление и температура, и свойства жидкости, такие как плотность и вязкость, заданы в этих узлах.
Обратите внимание, что инерция потока незначительна, и поток рассматривается как квазистационарное состояние. Перемещение переходных процессов к массовым расходам жидкости может быть смещено: из-за связи между плотностью, давлением и температурой распространение изменений по всей системе не мгновенно. Другие источники и стоки импульса, такие как различия в головке между портами или радиальные деформации стенки канала, не рассматриваются. Баланс импульса для половины объема в порту A:
где p - давление в узле, указанном в индексе. Δp f, A является общим падением давления между узлом порта и внутренним узлом из-за трения. Общее падение давления включает как основные, так и незначительные потери. Для половины объема в порту B, баланс импульса:
Скачки давления от трения изменяются с квадратом массового расхода жидкости для турбулентных потоков и с величиной массового расхода жидкости для ламинарных течений. Этот скачок давления характеризуется трехмерными параметрами: коэффициент трения Дарси, коэффициент падения давления и число Эйлера. Эти числа вычисляются из эмпирических корреляций или оцениваются из интерполяционных таблиц, в зависимости от параметра Pressure loss parameterization.
Классификация «ламинарного» или «турбулентного» потока основана на числе Рейнольдса. Когда число Рейнольдса выше параметра Turbulent flow lower Reynolds number limit, поток полностью турбулентен. Ниже параметра Laminar flow upper Reynolds number limit поток полностью ламинарен. Числа Рейнольдса, промежуточные этим значениям, указывают на переходный поток. Переходные потоки показывают характеристики как ламинарного, так и турбулентного потоков. На Simscape™ Fluids™ языке между этими ограничивающими значениями применяется численное смешение.
Correlations for tubes
Для труб используется коэффициент трения Дарси, f D. В половине объема в порту A баланс импульса :
где L - длина трубки, а L Add - добавленная длина трубки, которая воспроизводит незначительные вязкие потери, если используется вместо колен, тройников, объединений или других локальных сопротивлений. A - площадь поперечного сечения трубы; в случае неоднородной площади поперечного сечения следует использовать A мин. D H является гидравлическим диаметром трубы или диаметром круга, равным по площади сечению трубы :
Если трубка имеет круглое сечение, гидравлический диаметр и диаметр трубки совпадают.
Для половины объема в порту B, баланс импульса:
Для турбулентных потоков коэффициент трения Дарси вычисляется корреляцией Haaland. Число Рейнольдса устанавливается в ограничивающем порте:
где ε R - шероховатость стенки, принятая за характерную высоту. Этот параметр задан в параметре Internal surface absolute roughness.
Для ламинарных течений коэффициент трения зависит от формы трубы и вычисляется с масштабный фактор трубы:
где λ - масштабный фактор. Число Рейнольдса вычисляется в ограничивающем порте как:
Подставляя Re в уравнение падения давления в порту A, баланс импульса переформулируется как:
Точно так же для половины объема в порту B, баланс импульса:
Для каналов, отличных от труб, используйте коэффициент падения давления, ξ. Для турбулентных потоков в половине объема в порту A, баланс импульса:
Для турбулентных потоков в половине объема в порту B, баланс импульса:
Для ламинарных течений в половине объема в порту A, баланс импульса:
где ReL является Laminar flow upper Reynolds number limit параметров блоков. Для ламинарных течений в половине объема в порту B, баланс импульса:
Tabulated data - Darcy friction factor vs. Reynolds number
Можно использовать табличные данные для определения коэффициента трения Дарси на основе числа Рейнольдса для потоков труб. Для половины объема в порту A, баланс импульса:
Для половины объема в порту B, баланс импульса:
Для турбулентного режима коэффициент трения определяется из табличной функции числа Рейнольдса:
Точки прерывания табличной функции получают из вектора параметров блоков. Параметр Reynolds number vector for Darcy friction factor задает независимую переменную, а параметр Darcy friction factor vector задает зависимую переменную. Между точками останова применяется линейная интерполяция. За пределами табличных данных области значений ближайшая точка останова определяет коэффициент трения.
В ламинарном режиме коэффициент трения вычисляется из масштабного фактора,
Tabulated data - Euler number vs. Reynolds number
Можно использовать табличные данные для определения числа Эйлера на основе числа Рейнольдса. Это вычисление зависит от режима течения, и число Эйлера сформулировано как табличная функция числа Рейнольдса:
Точки останова в Tabulated data - Euler number vs. Reynolds number
заданы векторами числа Рейнольдса и числа Эйлера. Параметр Reynolds number vector for Euler number задает независимые переменные, числа Рейнольдса, а параметр Euler number vector задает зависимую переменную, число Эйлера, для каждого числа Рейнольдса. Линейная интерполяция используется для определения значений между точками останова. Вне табличных данных области значений используется значение в ближайшей точке останова.
Для турбулентных потоков баланс импульса для половины объема в порте A равен:
где Eu - число Эйлера в порту A. Для турбулентных потоков баланс импульса для половины объема в порте B равен:
Для ламинарного течения в половине объема в порту A, баланс импульса:
где ReL является параметром Laminar flow upper Reynolds number limit, а EuL - числом Эйлера, оцененным из табличных данных при этом числе Рейнольдса. Для ламинарного течения в половине объема в порту B, баланс импульса:
Энергетический баланс внутри объема газа является суммой его скоростей потока жидкости через контуры канала и связанной с этим теплопередачей. Энергия может быть передана путем консультирования в портах и путем конвекции в стенке. Хотя проводимость способствует энергетическому балансу в портах, она часто незначительна по сравнению с advection. Однако проводимость не незначительна в почти стационарных жидкостях, таких как, когда жидкости застойны или изменяют направление. Уравнение энергетического баланса:
где:
- частная производная внутренней энергии объема газа от давления при постоянной температуре и объеме.
p I - давление объема газа .
- частная производная внутренней энергии объема газа от температуры при постоянном давлении и объеме.
T I - температура объема газа .
Φ A и Φ B являются скоростями потока энергии в портах A и B, соответственно.
Q - скорость теплопередачи.
Приключение и проводимость учитываются в Φ, а конвекция - в Q. Скорость теплопередачи положительная при направлении в объем газа.
Теплопередача между двумя жидкостями теплообменника происходит несколькими способами: посредством конвекции на границах раздела жидкости, проводимости через слои заросшего загрязнения и проводимости через толщину стенки.
Теплопередача проходит за пределы газового канала и, следовательно, требует других блоков для моделирования всей системы теплообменника. Второй блок интерфейса теплообменника моделирует второй канал потока, в то время как блок E-NTU Heat Transfer моделирует тепловой поток через стенку. Параметры теплопередачи, которые специфичны для газового канала, но требуются блоком E-NTU Heat Transfer, доступны через порты физического сигнала:
Порт C выводит скорость теплоемкости, которая является мерой способности газа поглощать тепло и необходима для вычисления количества модулей (NTU). Скорость теплоемкости вычисляется как:
где C R - скорость теплоемкости, а c p - удельная теплота.
Порт HC выводит коэффициент теплопередачи, U.
Если коэффициент теплопередачи рассматривается как константа, его значение равномерно по каналу потока. Если коэффициент теплопередачи переменен, он вычисляется для каждого порта из выражения:
где Nu - число Нуссельта, k - теплопроводность, а D H, Q - гидравлический диаметр для теплопередачи. Гидравлический диаметр D H, Q вычисляется как:
где S Q является параметром Heat transfer surface area, а L Q является параметром Length of flow path for heat transfer.
Нижняя граница среднего коэффициента теплопередачи является параметром Minimum gas-wall heat transfer coefficient.
Число Нуссельта определяется эмпирическими корреляциями с числами Рейнольдса и Прандтля. Используйте параметр Heat transfer parametrization, чтобы выбрать наиболее подходящую формулировку для вашей симуляции.
Простейшая параметризация, Constant heat transfer coefficient
, получает коэффициент теплопередачи непосредственно от значения параметра Gas-wall heat transfer coefficient . Correlations for tubes
использует аналитические выражения с постоянными или вычисленными параметрами, чтобы захватить зависимость числа Нуссельта от режима течения для потоков труб.
Остальные параметризации являются табличными функциями числа Рейнольдса. Они полезны для изменения чисел Нуссельта или коэффициентов теплопередачи между режимами течения. Функции сгенерированы из экспериментальных данных, связывающих число Рейнольдса с коэффициентом Колборна или числа Рейнольдса и Прандтля с числом Нуссельта.
Constant heat transfer coefficient
Использование Constant heat transfer coefficient
, заданный в параметре Gas-wall heat transfer coefficient, устанавливает коэффициент теплопередачи как константу, и не использует число Нуссельта в вычислениях. Используйте эту параметризацию как простое приближение для газовых потоков, ограниченных ламинарным режимом.
Correlation for tubes
Число Нуссельта зависит от режима течения при использовании Correlation for tubes
. Для турбулентных потоков его значение изменяется пропорционально числу Рейнольдса и вычисляется из корреляции Гнилинского:
где Re - число Рейнольдса, Nu - число Нуссельта, а Pr - число Прандтля. Коэффициент трения, f, тот же, что и коэффициент, используемый в вычислениях падения давления в трубе. Для ламинарных течений число Нуссельта является константой. Его значение получается из параметра Nusselt number for laminar flow heat transfer, NuL:
Tabulated data - Colburn factor vs. Reynolds number
Можно использовать табличные данные для определения коэффициента Колборна на основе числа Рейнольдса. Уравнение Колберна используется для определения числа Нуссельта, которое изменяется пропорционально числу Рейнольдса. Colburn j -factor является мерой пропорциональности между числами Рейнольдса, Прандтля и Нуссельта:
ReQ - число Рейнольдса, основанное на гидравлическом диаметре теплопередачи, D H, Q и от минимальной площади свободного потока канала, A Min:
Tabulated data - Nusselt number vs. Reynolds number & Prandtl number
Можно использовать табличную функцию, чтобы определить число Нуссельта из чисел Прандтля и Рейнольдса. Линейная интерполяция используется для определения значений между точками останова. Число Нуссельта является функцией как Re, так и Pr, и поэтому параметры Reynolds number vector for Nusselt number, Prandtl number vector for Nusselt number и Nusselt number table, Nu(Re,Pr) определяют точки останова таблицы:
Табличное число Рейнольдса должно быть вычислено с помощью гидравлического диаметра для теплопередачи, D H, Q.