Трубопровод для транспортировки изменяющей фазу жидкости с теплопередачей
Simscape/Жидкости/Двухфазная жидкость/Трубы и подборы кривой
Блок Pipe (2P) моделирует трубопровод с изменяющей фазу жидкостью. Каждая жидкая фаза называется зоной, которая является дробным значением между 0 и 1. Зоны не смешиваются. 3-Zone Pipe (2P) использует следующую за ним модель для отслеживания переохлаждаемой жидкости (L), смеси пар-жидкость (M) и перегретого пара (V) в трех зонах. Относительное количество пространства, которое зона занимает в системе, называется дробью длины зоны в системе.
Порт H является тепловым портом, который представляет температуру окружающей среды. Скорость теплопередачи между жидкостью и окружением зависит от жидкой фазы каждой зоны. Стенки трубопровода моделируются внутри блока, и температура стенки трубопровода в каждой зоне может быть различной. Давление и температура зависят от динамической сжимаемости жидкости и теплоемкости зоны жидкости.
Коэффициент конвективной теплопередачи между жидкостью и стенкой, αF, изменяется на зону в соответствии с числом Нуссельта:
где:
Nu
- номер зоны Нуссельта.
k среднюю теплопроводность жидкости.
D H - Hydraulic diameter трубопровода, эквивалентный диаметр некруглого сечения.
Число Нуссельта, используемое в коэффициенте теплопередачи, является большим из числа Нуссельта турбулентного и ламинарного потока.
Для турбулентных потоков в переохлажденных жидких или перегретых паровых зонах число Нуссельта вычисляется корреляцией Гнилинского:
где:
Re
- среднее Число Рейнольдса зоны.
Pr
- среднее значение Prandtl в зоне.
f - коэффициент трения Дарси, вычисленный из корреляции Haaland:
где ε - стенка Internal surface absolute roughness.
Для турбулентных потоков в зоне смеси пар-жидкости число Нуссельта вычисляется корреляцией Каваллини-Цекхина:
Где:
Re SL является числом Рейнольдса насыщенной жидкости.
Pr SL является числом Прандтля насыщенной жидкости.
ρ SL является плотностью насыщенной жидкости.
ρ SV является плотностью насыщенного пара.
a = 0,05, b = 0,8 и c = 0,33.
Когда ребра моделируются на внутренней поверхности трубопровода, коэффициент теплопередачи равен:
где:
η Int - это Internal fin efficiency.
s Int - это Ratio of internal fins surface area to no-fin surface area.
Для ламинарных течений число Нуссельта задается параметром Laminar flow Nusselt number.
Скорость теплопередачи от жидкости основана на изменении конкретной энтальпии в каждой зоне:
где - массовый расход жидкости теплопередачи. Это входное отверстие трубопровода массового расхода жидкости, или A или B, в зависимости от направления потока жидкости.
В жидкой и паровой зонах изменение специфической энтальпии определяется как:
где:
cp - удельная теплота жидкости или пара.
TH - температура окружающей среды.
TI - температура входного отверстия жидкости.
z - фракция длины зоны жидкости.
αE - коэффициент теплопередачи между стенкой и окружением.
SW - площадь поверхности стенки:
где:
A - Cross-sectional area трубопровода.
L является Pipe length.
Обратите внимание, что эта площадь поверхности стенки не включает площадь ребра, которая определяется параметрами Ratio of external fins surface area to no-fin surface area и Ratio of internal fins surface area to no-fin surface area. Ребра расположены пропорционально площади поверхности стенки. Значение 0
означает отсутствие ребер на стенке трубопровода.
В зоне смеси пар-жидкость изменение специфической энтальпии рассчитывается как:
где TS - температура насыщения жидкости. Принято, что смесь пар-жидкость всегда находится при этой температуре.
Общая теплопередача между жидкостью и стенкой трубопровода является суммой теплопередачи в каждой жидкой фазе:
Скорость теплопередачи между жидкостью и трубопроводом в жидкой зоне:
где TW,L - температура стенки, окружающей жидкую зону.
Скорость теплопередачи между жидкостью и трубопроводом в зоне смеси:
Скорость теплопередачи между жидкостью и трубопроводом в паровой зоне:
где TW,V - температура стенки, окружающей паровую зону.
Если стенка трубопровода имеет конечную толщину, коэффициент теплопередачи между стеной и окружением, αE, определяется:
где α W - коэффициент теплопередачи, обусловленный проводимостью через стенку:
и где:
kW является Wall thermal conductivity.
tW является Wall thickness.
α Ext - это External environment heat transfer coefficient.
η Ext - это External fin efficiency.
s Ext - это Ratio of external fins surface area to no-fin surface area.
Если стенка не имеет тепловой массы, коэффициент теплопередачи между стеной и окружением равен коэффициенту теплопередачи окружения, α Ext.
Скорость теплопередачи между каждой зоной стенки и окружения:
Общая теплопередача между стенкой и окружением составляет:
Скорость теплопередачи зависит от тепловой массы стенки, CW:
где:
c p, W является Wall specific heat.
ρW является Wall density.
Руководящими уравнениями для теплопередачи между жидкостью и внешним окружением являются, для жидкой зоны:
для зоны смеси:
и для паровой зоны:
Перепад давления по трубе состоит из двух факторов: изменения давления из-за изменений плотности и изменения давления из-за трения на стенках трубопровода.
Для турбулентных потоков, когда число Рейнольдса выше Turbulent flow lower Reynolds number limit, падение давления вычисляется с точки зрения коэффициента трения Дарси. Перепад давления между портом A и внутренним узлом I:
где:
ρI - плотность жидкости во внутреннем узле I.
ρA* - плотность жидкости в порте A. Это то же самое, что и ρA, когда поток является установившимся; когда поток является переходным, он вычисляется из внутреннего состояния жидкости с помощью адиабатического выражения:
где:
h - средняя специфическая энтальпия,
ρ - средняя плотность,
Это связано с тем, что расчет теплопередачи происходит на внутреннем узле I.
A - массовый расход жидкости через порт A.
L является Pipe length.
L Add является Aggregate equivalent length of local resistances, которая является эквивалентной длиной трубки, которая вводит такое же количество потерь, как и сумма потерь из-за других локальных сопротивлений.
Обратите внимание, что коэффициент трения Дарси зависит от числа Рейнольдса и вычисляется в обоих портах.
Перепад давления между портом B и внутренним узлом I:
где:
ρB* - плотность жидкости в порте B. Это то же самое, что и ρB, когда поток является установившимся; когда поток является переходным, он вычисляется из внутреннего состояния жидкости с помощью адиабатического выражения:
B - массовый расход жидкости через порт B.
Для ламинарных течений, когда число Рейнольдса ниже Laminar flow upper Reynolds number limit, падение давления на трение вычисляется в терминах Laminar friction constant for Darcy friction factor, λ. Перепад давления между портом A и внутренним узлом I:
where - средняя динамическая вязкость жидкости:
Перепад давления между портом B и внутренним узлом I:
Для переходных потоков перепадом давления из-за вязкого трения является сглаженная смесь между значениями для ламинарного и турбулентного падения давления.
Общая скорость накопления массы определяется как:
где M - общая масса жидкости в трубопроводе. Что касается зон жидкости, скорость накопления массы является функцией изменения плотности, ρ, относительно давления, p и удельной внутренней энергии жидкости, u:
где:
uout является специфической внутренней энергией после того, как произошла вся теплопередача.
V - общий объем жидкости или объем трубопровода.
Уравнение энергетического разговора:
где:
ϕ A является расходом энергии в порту A.
ϕ B - расход энергии в порту B.
Q F - скорость теплопередачи между жидкостью и стенкой.
Стенки трубопровода идеально жесткие.
Поток полностью развит. Потери на трение и теплопередача не включают входные эффекты.
Инерция жидкости незначительна.
Эффект тяжести незначителен.
Когда давление выше критического давления жидкости, большие значения свойств тепловой жидкости (таких как число Прандтля, теплопроводность и удельное тепло) могут не точно отражать теплообмен в трубопроводе.
[1] Белый, F.M., Механика жидкости, 7th Ред, раздел 6.8. Макгроу-Хилл, 2011.
[2] Qengel, Y.A., Теплопередача и массопередача - практический подход, 3rd Ред, раздел 8.5. Макгроу-Хилл, 2007.