Твердый кабелепровод для потока жидкости в двухфазных гидросистемах
Двухфазная Жидкость/Элементы
Блок Pipe (2P) моделирует динамику потока двухфазной жидкости в твердом трубопроводе. Динамическая сжимаемость и тепловая мощность производства жидкости приняты ненезначительные. Двухфазные жидкие порты А сохранения и B представляют входы трубопровода. Тепловой порт H сохранения представляет стену трубопровода, через которую происходит теплопередача со средой трубопровода.
Блок предоставляет возможность к инерции жидкости модели, сопротивление внезапным изменениям в массовом расходе жидкости. По умолчанию моделирование инерции жидкости выключено. Эта установка является соответствующей, когда силы давления, управляющие потоком далеко, превышают инерционные силы, действующие на поток.
Настройка по умолчанию уменьшает вычислительные затраты и рекомендуется для большинства моделей. Однако инерция жидкости может стать важной, если массовый расход жидкости изменяется быстро. В таких случаях, включая моделирование инерции жидкости может помочь улучшить точность симуляции.
Энергосбережение в трубопроводе наблюдается посредством уравнения:
где:
M является жидкой массой в трубопроводе.
u я - определенная внутренняя энергия жидкости в трубопроводе.
ϕ 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 является объемом жидкости в трубопроводе.
массовый расход жидкости в трубопровод через порт А.
массовый расход жидкости в трубопровод через порт B.
∊ M является термином коррекции, который составляет сглаживание частных производных плотности через контуры перехода фазы.
Блок смешивает частные производные плотности различных областей с помощью функции кубического полинома. В качестве пара 0–0.1, эта функция смешивает производные подохлажденных жидких и двухфазных областей смеси. В качестве пара 0.9–1, это смешивает те из двухфазной смеси и перегретых областей пара. Термин коррекции в массовом уравнении сохранения,
добавляется, чтобы откорректировать для числовых ошибок, введенных функцией кубического полинома, с:
M как жидкая масса в трубопроводе, вычисленном из уравнения:
v I как определенный объем жидкости в трубопроводе.
τ как постоянная времени фазового перехода — характеристическая длительность события фазового перехода. Эта константа гарантирует, что фазовые переходы не происходят мгновенно, эффективно вводя задержку каждый раз, когда они происходят.
Уравнения баланса импульса определены отдельно для каждого раздела хаф-пайпа. В хаф-пайпе, смежном с портом А:
где:
p A является давлением в порте А.
S является площадью поперечного сечения трубопровода.
ν A является определенным объемом жидкости в порте А.
F visc, A является вязкой силой трения в хаф-пайпе, смежном с портом А.
I 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, заданное в диалоговом окне блока — вязкая сила трения в хаф-пайпе, смежном с портом А,
в то время как в хаф-пайпе, смежном с портом B, это
где:
Форма f является масштабным фактором трубопровода.
Эффективность L является эффективной длиной трубопровода — сумма длины трубопровода и совокупная эквивалентная продолжительность локальных сопротивлений.
D h является гидравлическим диаметром трубопровода.
В турбулентном режиме — то есть, когда число Рейнольдса больше значения Turbulent flow lower Reynolds number limit, заданного в диалоговом окне блока — вязкая сила трения в хаф-пайпе, смежном с портом А,
в то время как в хаф-пайпе, смежном с портом B, это
где:
f A является коэффициентом трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом А.
f B является коэффициентом трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом B.
Коэффициент трения Дарси для турбулентного течения в хаф-пайпе, смежном с портом А, следует из уравнения Haaland как
и в хаф-пайпе, смежном с портом B как
где:
∊ r является относительной шероховатостью трубопровода.
Re A является числом Рейнольдса в хаф-пайпе, смежном с портом А,
Re B является числом Рейнольдса в хаф-пайпе, смежном с портом B,
Функция кубического полинома используется, чтобы смешать потери на трение в области перехода между ламинарным течением и турбулентным течением.
Стена трубопровода тверда.
Поток полностью разрабатывается.
Эффект силы тяжести незначителен.
Теплопередача вычисляется относительно температуры объема жидкости в трубопроводе. Чтобы смоделировать градиент температуры из-за теплопередачи вдоль длинного трубопровода, соедините несколько блоков Pipe (2P) последовательно.
Расстояние между входом трубопровода и выходом. Значением по умолчанию является 5 m.
Внутренняя область трубопровода, нормальная к направлению потока. Эта область является постоянной вдоль трубопровода. Значением по умолчанию является 0.01 м^2.
Диаметр эквивалентной трубы круглого сечения. В цилиндрическом трубопроводе гидравлический диаметр совпадает со своим фактическим диаметром. Значением по умолчанию является 0.1 m.
Падение давления из-за локальных сопротивлений, таких как повороты, входы и подборы кривой, описанные как эквивалентная продолжительность этих сопротивлений. Значением по умолчанию является 0.1 m.
Средняя высота всей поверхности дезертирует на внутренней поверхности трубопровода. Этот параметр включает вычисление коэффициента трения в режиме турбулентного течения. Значением по умолчанию является 1.5e-5m.
Самое большое значение числа Рейнольдса, соответствующего полностью разработанному ламинарному течению. Когда число Рейнольдса повышается выше этого предела, поток постепенно переходы от ламинарного к турбулентному. Значением по умолчанию является 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
