Жесткий трубопровод для потока жидкости в двухфазных гидросистемах
Двухфазная жидкость/элементы
Блок Pipe (2P) моделирует динамику потока двухфазной жидкости внутри жёсткого трубопровода. Динамическая сжимаемость и теплоемкость жидкости приняты незначительными. Двухфазные порты А и B для сохранения жидкости представляют входные отверстия трубопровода. Тепловой порт H представляет стенку трубопровода, через которую происходит теплопередача с окружающей трубой.
Блок предоставляет опцию для моделирования инерции жидкости, сопротивления внезапным изменениям в массовом расходе жидкости. По умолчанию моделирование инерции жидкости отключено. Эта настройка подходит, когда силы давления, движущие поток далеко, превышают инерционные силы, действующие на поток.
Настройка по умолчанию снижает вычислительные затраты и рекомендуется для большинства моделей. Однако инерция жидкости может стать важной, если массовый расход жидкости изменяется быстро. В таких случаях включение моделирования инерции жидкости может помочь улучшить точность симуляции.
Энергосбережение в трубопроводе наблюдается уравнением:
где:
M - масса жидкости в трубопроводе.
u I является специфической внутренней энергией жидкости в трубопроводе.
ϕ A - это расход энергии в трубопровод через порт A.
ϕ B является расходом энергии в трубопровод через порт B.
Q H - скорость теплового потока в трубу через стенку трубопровода, представленная портом H.
Теплопередача между стенкой трубопровода и внутренним объемом жидкости моделируется как конвективный процесс с скоростью потока жидкости тепла, рассчитанным как:
где:
h коэффициент является средним коэффициентом теплопередачи в трубопроводе.
S Surf является площадью поверхности трубопровода.
T H - температура стенки трубопровода.
T I - температура жидкости в трубопроводе.
Вычисление коэффициента теплопередачи зависит от фазы жидкости. В переохлажденной жидкой и перегретой паровой фазах коэффициент равен:
где звездочка обозначает значение, характерное для рассматриваемой фазы (жидкость или пар) и:
Nu - среднее число Нуссельта в трубопроводе.
kЯ - средняя теплопроводность в трубопроводе.
D h - гидравлический диаметр трубопровода (то, что сечение общей формы имело бы, если бы оно было выполнено круглым).
В двухфазной смеси тот же коэффициент является:
где индекса M обозначает значение, характерное для двухфазной смеси и SL индекс указывает значение, полученное для насыщенной жидкости.
В ламинарных течениях число Нуссельта принимается постоянным и равным значению, заданному в диалоговом окне блока. Ламинарное течение применяется, когда число Рейнольдса меньше значения, введенного для параметра Laminar flow upper Reynolds number limit.
Число Нуссельта турбулентного потока применяется, когда число Рейнольдса больше значения, введенного для параметра Laminar flow upper Reynolds number limit. В переходной области между ламинарным и турбулентным потоками кубическая полиномиальная функция смешивает два числа Нуссельта. Это смешение обеспечивает плавный переход между режимами течения.
В жидкой и паровой фазах число Нуссельта для турбулентного потока следует из корреляции Гнилинского:
где, как и прежде, звездочка обозначает рассматриваемую фазу и:
f - коэффициент трения трубопровода.
Re - число Рейнольдса.
Pr я - это число Прандтля.
Коэффициент трения вычисляется как:
где ε r - шероховатость трубопровода. Число Рейнольдса вычисляется как:
где индекса Avg обозначает среднее значение между портами и:
S - площадь поперечного сечения трубопровода.
v I - конкретный объем.
ν I - кинематическая вязкость.
В двухфазной смеси число Нуссельта для турбулентного потока следует из корреляции Каваллини и Цекхина:
где индекса SL обозначает значение для насыщенной жидкости, SV опустить значение насыщенного пара и:
x I - качество пара.
v - конкретный объем.
Число Рейнольдса насыщенной жидкости вычисляется как:
Сохранение массы в трубопроводе наблюдается уравнением:
где:
ρ - плотность жидкости.
p я давление в трубопроводе.
V - объем жидкости в трубопроводе.
массовый расход жидкости в трубопровод через порт А.
- массовый расход жидкости в трубопровод через порт B.
∊ M является поправочным термином, который учитывает сглаживание частных производных плотности через контуры фазового перехода.
Блок смешивает частные производные плотности различных областей с помощью кубической полиномиальной функции. При качестве пара 0-0,1 эта функция смешивает производные переохлажденной жидкой и двухфазной смесевых областей. При качестве пара 0,9-1 он смешивает пары двухфазной смеси и перегретых паровых областей. Член коррекции в уравнении сохранения массы,
добавляется для исправления числовых ошибок, введенных кубической полиномиальной функцией, с:
M как масса жидкости в трубопроводе, вычисленная из уравнения:
v I как конкретный объем жидкости в трубопроводе.
τ как постоянная времени изменения фазы - характеристическая длительность события изменения фазы. Эта константа гарантирует, что изменения фазы не происходят мгновенно, эффективно вводя временную задержку каждый раз, когда они происходят.
Уравнения баланса импульса заданы отдельно для каждой половины участка трубопровода. В половине трубопровода, примыкающего к порту А:
где:
p A - давление в порту A.
S - площадь поперечного сечения трубопровода.
ν A - удельный объем жидкости в порту A.
F visc, A является вязкой силой трения в половине трубопровода, примыкающей к порту 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 eff - эффективная длина трубопровода - сумма длины трубопровода и совокупная эквивалентная длина локальных сопротивлений.
D h - гидравлический диаметр трубопровода.
В турбулентном режиме - то есть, когда число Рейнольдса больше Turbulent flow lower Reynolds number limit значения, заданного в диалоговом окне блока - вязкая сила трения в половине трубопровода, соседствующей с портом А,
в то время как в половине трубопровода, соседствующего с портом B, это
где:
f A является коэффициентом трения Дарси для турбулентного потока в половине трубопровода, примыкающего к порту A.
f B является коэффициентом трения Дарси для турбулентного потока в половине трубопровода, примыкающего к порту B.
Коэффициент трения Дарси для турбулентного потока в половине трубопровода, соседствующего с портом А, следует из уравнения Haaland как
и в половине трубопровода, примыкающего к порту B, как
где:
∊ r - относительная шероховатость трубопровода.
Re A - число Рейнольдса в полупроводнике, соседствующем с портом А,
Re B - число Рейнольдса в полупроводнике, соседствующем с портом B,
Кубическая полиномиальная функция используется, чтобы смешать потери на трение в переходной области между ламинарным течением и турбулентным потоком.
Стенки трубопровода жесткие.
Поток полностью развит.
Эффект тяжести незначителен.
Теплопередача вычисляется относительно температуры объема жидкости в трубопроводе. Чтобы смоделировать градиент температуры из-за теплопередачи по длинному трубопроводу, соедините несколько блоков Pipe (2P) последовательно.
Расстояние между входным и выходным отверстиями трубопровода. Значение по умолчанию 5 м.
Площадь внутреннего трубопровода, перпендикулярная направлению потока. Эта площадь постоянна по длине трубопровода. Значение по умолчанию 0.01 м ^ 2.
Диаметр эквивалентной трубы круглого сечения. В цилиндрическом трубопроводе гидравлический диаметр совпадает с его фактическим диаметром. Значение по умолчанию 0.1 м.
Падение давления из-за локальных сопротивлений, таких как повороты, входные отверстия и подборы кривой, выраженное в виде эквивалентной длины этих сопротивлений. Значение по умолчанию 0.1 м.
Средняя высота всех поверхностных дефектов на внутренней поверхности трубопровода. Этот параметр включает вычисление коэффициента трения в турбулентном режиме течения. Значение по умолчанию 1.5e-5м.
Наибольшее значение числа Рейнольдса, соответствующего полностью развитому ламинарному течению. Когда число Рейнольдса увеличивается выше этого предела, поток постепенно переходит от ламинарного к турбулентному. Значение по умолчанию 2000.
Наименьшее значение числа Рейнольдса, соответствующее полностью развитому турбулентному потоку. Когда число Рейнольдса падает ниже этого предела, поток постепенно переходит от турбулентного к ламинарному. Значение по умолчанию 4000.
Полу-эмпирический параметр, кодирующий эффект геометрии трубопровода на вязкие потери на трение, происходящие в ламинарном режиме. Соответствующее значение для использования зависит от формы поперечного сечения трубопровода.
Типичные значения включают 56 для квадратного сечения, 62 для прямоугольного сечения и 96 для концентрического поперечного сечения кольцевого пространства [1]. Значение по умолчанию, соответствующее круглое сечение, является 64.
Коэффициент пропорциональности конвективной и проводящей теплопередачи в ламинарном режиме. Этот параметр позволяет вычислять конвективную теплопередачу в ламинарных течениях. Его значение изменяется с площадью поперечного сечения по каналу и тепловыми граничными условиями, например, постоянной температурой или постоянным тепловым потоком на стенке трубопровода. Значение по умолчанию, соответствующее круглое сечение сечению, является 3.66.
Опция для моделирования инерции жидкости, сопротивления жидкости быстрому ускорению. Значение по умолчанию является Off.
Термодинамическая переменная, в терминах которой можно задать начальные условия компонента. Настройкой по умолчанию является Temperature.
Давление в ёмкости в начале симуляции, заданное относительно нуля. Значение по умолчанию 0.101325 МПа.
Температура в ёмкости в начале симуляции, заданная при абсолютном нуле. Этот параметр активен, когда опция 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 кДж/кг.
Удельная внутренняя энергия жидкости в ёмкости в начале симуляции. Этот параметр активен, когда опция Initial fluid energy specification установлена на Specific internal energy. Значение по умолчанию 1500 кДж/кг.
Характеристическая длительность события изменения фазы. Эта константа вводит временную задержку в переход между фазами. Значение по умолчанию 0.1 с.
Блок имеет два двухфазных порта для жидкости, A и B. Порт H является тепловым портом, представляющим стенку трубопровода, через которую происходит теплообмен.
[1] Белый, F.M., Вязкий Поток Жидкости, Макгроу-Хилл, 1991
