Твердый кабелепровод для потока жидкости в двухфазных жидких системах
Двухфазная Жидкость/Элементы
Блок Pipe (2P) моделирует динамику потока двухфазной жидкости в твердом канале. Динамическая сжимаемость и тепловая мощность производства жидкости приняты ненезначительные. Двухфазные жидкие порты A и B сохранения представляют входные отверстия канала. Тепловой порт H сохранения представляет стену канала, через которую происходит теплопередача со средой канала.
Блок предоставляет возможность моделировать жидкую инерцию, сопротивление внезапным изменениям в массовой скорости потока жидкости. По умолчанию жидкое моделирование инерции выключено. Эта установка является соответствующей, когда силы давления, управляющие потоком далеко, превышают инерционные силы, действующие на поток.
Настройка по умолчанию уменьшает вычислительные затраты и рекомендуется для большинства моделей. Однако жидкая инерция может стать важной, если массовая скорость потока жидкости изменяется быстро. В таких случаях, включая жидкое моделирование инерции может помочь улучшить точность симуляции.
Энергосбережение в канале наблюдается посредством уравнения:
где:
M является жидкой массой в канале.
u я - определенная внутренняя энергия жидкости в канале.
ϕ A является энергетической скоростью потока жидкости в канал через порт A.
ϕ B является энергетической скоростью потока жидкости в канал через порт B.
Q H является уровнем теплового потока в канал через стену канала, представленную портом H.
Теплопередача между стеной канала и внутренним жидким объемом моделируется как конвективный процесс с уровнем теплового потока, вычисленным как:
где:
Коэффициент h является средним коэффициентом теплопередачи в канале.
Surf S является площадью поверхности канала.
T H является температурой канала стенки.
T я - температура жидкости в канале.
Вычисление коэффициента теплопередачи зависит от жидкой фазы. В подохлажденных жидких и перегретых фазах пара коэффициент:
где звездочка обозначает значение, характерное для рассмотренной фазы (жидкость или пар) и:
Nu
является средним номером Nusselt в канале.
k
я - средняя теплопроводность в канале.
D h является гидравлическим диаметром канала (то, что имело бы сечение общей формы, если бы это было сделано круговым).
В двухфазной смеси тот же коэффициент:
где нижний M
обозначает значение, характерное для двухфазной смеси, и индекс SL
указывает на значение, полученное для влажной жидкости.
В ламинарных течениях номер Nusselt принят постоянный и равный значению, заданному в диалоговом окне блока. Ламинарное течение, которое применяет номер Nusselt, когда число Рейнольдса меньше, чем значение, ввело для параметра Laminar flow upper Reynolds number limit.
Турбулентное течение, которое применяет номер Nusselt, когда число Рейнольдса больше, чем значение, ввело для параметра Laminar flow upper Reynolds number limit. В транзитной области между ламинарным и турбулентным течением функция кубического полинома смешивает два числа Nusselt. Это смешивание гарантирует плавный переход между режимами потока.
В жидкости и фазах пара, номер Nusselt для турбулентного течения следует из корреляции Гниелинского:
где, как прежде, звездочка обозначает рассмотренную фазу и:
f является фактором трения канала.
Re
является числом Рейнольдса.
Pr я - число Прандтля.
Фактор трения вычисляется как:
где ε r является шероховатостью канала. Число Рейнольдса вычисляется как:
где нижний Avg
обозначает среднее значение между портами и:
S является площадью поперечного сечения канала.
v я - определенный объем.
ν я - кинематическая вязкость.
В двухфазной смеси номер Nusselt для турбулентного течения следует из корреляции Каваллини и Зеччина:
где нижний SL
обозначает значение для влажной жидкости, SV
преобразовывают значение в нижний индекс для влажного пара, и:
x я - качество пара.
v является определенным объемом.
Число Рейнольдса влажной жидкости вычисляется как:
Массовое сохранение в канале наблюдается посредством уравнения:
где:
ρ является жидкой плотностью.
p я - давление в канале.
V является объемом жидкости в канале.
массовая скорость потока жидкости в канал через порт A.
массовая скорость потока жидкости в канал через порт B.
∊ M является сроком исправления, который составляет сглаживание частных производных плотности через контуры перехода фазы.
Блок смешивает частные производные плотности различных областей с помощью функции кубического полинома. В качестве пара 0–0.1, эта функция смешивает производные подохлажденных жидких и двухфазных областей смеси. В качестве пара 0.9–1, это смешивает те из двухфазной смеси и перегретых областей пара. Срок исправления в массовом уравнении сохранения,
добавляется, чтобы исправить для числовых ошибок, введенных функцией кубического полинома, с:
M как жидкая масса в канале, вычисленном из уравнения:
v I как определенный объем жидкости в канале.
τ как временная константа фазового перехода — характеристическая длительность события фазового перехода. Эта константа гарантирует, что фазовые переходы не происходят мгновенно, эффективно вводя задержку каждый раз, когда они происходят.
Уравнения баланса импульса определены отдельно для каждого раздела хаф-пайпа. В хаф-пайпе, смежном с портом A:
где:
p A является давлением в порте A.
S является площадью поперечного сечения канала.
ν A является определенным объемом жидкости в порте A.
F visc, A является вязкой силой трения в хаф-пайпе, смежном с портом A.
I A является жидкой инерцией в порте A:
Параметр L является длиной канала.
В хаф-пайпе, смежном с портом B:
где:
p B является давлением в порте B.
ν B является определенным объемом жидкости в порте B.
F visc, B является вязкой силой трения в хаф-пайпе, смежном с портом B.
I B является жидкой инерцией в порте B:
Жидкие условия инерции, I A и I B, являются нулем, когда параметр Fluid inertia устанавливается на Off
. Вычисление вязких сил трения, F visc, A и F visc, B зависит от режима потока, пластинчатого или бурного.
В пластинчатом режиме — то есть, когда число Рейнольдса меньше, чем значение Laminar flow upper Reynolds number limit, заданное в диалоговом окне блока — вязкая сила трения в хаф-пайпе, смежном с портом A,
в то время как в хаф-пайпе, смежном с портом B, это
где:
Форма f является форм-фактором канала.
Эффективность L является эффективной длиной канала — сумма длины канала и совокупная эквивалентная продолжительность локальных сопротивлений.
D h является гидравлическим диаметром канала.
В бурном режиме — то есть, когда число Рейнольдса больше, чем значение Turbulent flow lower Reynolds number limit, заданное в диалоговом окне блока — вязкая сила трения в хаф-пайпе, смежном с портом A,
в то время как в хаф-пайпе, смежном с портом B, это
где:
f A является фактором трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом A.
f B является фактором трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом B.
Фактор трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом A, следует от уравнения Haaland как
и в хаф-пайпе, смежном с портом B как
где:
∊ r является относительной шероховатостью канала.
Re A является числом Рейнольдса в хаф-пайпе, смежном с портом A,
Re B является числом Рейнольдса в хаф-пайпе, смежном с портом B,
Функция кубического полинома используется, чтобы смешать потери трения в области перехода между ламинарным течением и турбулентным течением.
Диалоговое окно блока не имеет ссылки Source code. Чтобы просмотреть исходный код для различных настроек блока, откройте следующие файлы в редакторе MATLAB®:
Передайте по каналу с потерями и динамической сжимаемостью (жидкая инерция Off
) — pipe_resistive_dynamic.ssc
Передайте по каналу с потерями, динамической сжимаемостью и жидкой инерцией (жидкая инерция On
) — pipe_resistive_dynamic_inertia.ssc
Стена канала тверда.
Поток полностью разрабатывается.
Эффект силы тяжести незначителен.
Теплопередача вычисляется относительно температуры жидкого объема в канале. Чтобы смоделировать градиент температуры из-за теплопередачи вдоль длинного канала, соедините несколько Канал (2P) блоки последовательно.
Расстояние между входным отверстием канала и выходом. Значением по умолчанию является 5
m.
Внутренняя область канала, нормальная к направлению потока. Эта область является постоянной вдоль канала. Значением по умолчанию является 0.01
m^2.
Диаметр эквивалентного канала с круговым сечением. В цилиндрическом канале гидравлический диаметр совпадает со своим фактическим диаметром. Значением по умолчанию является 0.1
m.
Падение давления из-за локальных сопротивлений, таких как изгибы, входные отверстия и подборы кривой, выраженные как эквивалентная продолжительность этих сопротивлений. Значением по умолчанию является 0.1
m.
Средняя высота всей поверхности дезертирует на внутренней поверхности канала. Этот параметр включает вычисление трения, включают режим турбулентного течения. Значением по умолчанию является 1.5e-5
m.
Самое большое значение числа Рейнольдса, соответствующего полностью разработанному ламинарному течению. Когда число Рейнольдса повышается выше этого предела, поток постепенно переходы от пластинчатого до бурного. Значением по умолчанию является 2000
.
Наименьшее значение числа Рейнольдса, соответствующего полностью разработанному турбулентному течению. Когда число Рейнольдса падает ниже этого предела, поток постепенно переходы от бурного до пластинчатого. Значением по умолчанию является 4000
.
Полуэмпирический параметр, кодирующий эффект геометрии канала на вязких потерях трения, понесен в пластинчатом режиме. Соответствующее значение, чтобы использовать зависит от перекрестной частной формы канала.
Типичные значения включают 56
для квадратного сечения, 62
для прямоугольного сечения и 96
для концентрического сечения кольца [1]. Значением по умолчанию, соответствуя круговому сечению, является 64
.
Коэффициент пропорциональности между конвективной и проводящей теплопередачей в пластинчатом режиме. Этот параметр включает вычисление конвективной теплопередачи в ламинарных течениях. Его значение изменяется с площадью поперечного сечения канала и тепловыми граничными условиями, например, постоянный температурный или постоянный поток тепла в стене канала. Значением по умолчанию, соответствуя круговому сечению канала, является 3.66
.
Опция, чтобы смоделировать жидкую инерцию, сопротивление жидкости к быстрому ускорению. Значением по умолчанию является Off
.
Термодинамическая переменная, с точки зрения которой можно задать начальные условия компонента. Настройкой по умолчанию является Temperature
.
Давление в камере в начале симуляции, заданной против абсолютного нуля. Значением по умолчанию является 0.101325
MPa.
Температура в камере в начале симуляции, заданной против абсолютного нуля. Этот параметр активен, когда опция Initial fluid energy specification установлена в Temperature
. Значением по умолчанию является 293.15
K.
Массовая часть пара в камере в начале симуляции. Этот параметр активен, когда опция Initial fluid energy specification установлена в Vapor quality
. Значением по умолчанию является 0.5
.
Часть объема пара в камере в начале симуляции. Этот параметр активен, когда опция Initial fluid energy specification установлена в Vapor void fraction
. Значением по умолчанию является 0.5
.
Определенная энтальпия жидкости в камере в начале симуляции. Этот параметр активен, когда опция Initial fluid energy specification установлена в Specific enthalpy
. Значение по умолчанию является 1500
kJ/kg.
Определенная внутренняя энергия жидкости в камере в начале симуляции. Этот параметр активен, когда опция Initial fluid energy specification установлена в Specific internal energy
. Значение по умолчанию является 1500
kJ/kg.
Характеристическая длительность события фазового перехода. Эта константа вводит задержку в переход между фазами. Значением по умолчанию является 0.1
s.
Блок имеет два двухфазных жидких порта сохранения, A и B. Порт H является тепловым портом сохранения, представляющим стену канала, через которую происходит теплообмен.
[1] Белый, F.M., вязкий поток жидкости, McGraw-Hill, 1991
Локальное ограничение (2P) | Переменное локальное ограничение (2P)