3-Zone Pipe (2P)

Трубопровод для транспортировки изменяющей фазу жидкости с теплопередачей

  • Библиотека:
  • Simscape/Жидкости/Двухфазная жидкость/Трубы и подборы кривой

  • 3-Zone Pipe (2P) block

Описание

Блок Pipe (2P) моделирует трубопровод с изменяющей фазу жидкостью. Каждая жидкая фаза называется зоной, которая является дробным значением между 0 и 1. Зоны не смешиваются. 3-Zone Pipe (2P) использует следующую за ним модель для отслеживания переохлаждаемой жидкости (L), смеси пар-жидкость (M) и перегретого пара (V) в трех зонах. Относительное количество пространства, которое зона занимает в системе, называется дробью длины зоны в системе.

Порт H является тепловым портом, который представляет температуру окружающей среды. Скорость теплопередачи между жидкостью и окружением зависит от жидкой фазы каждой зоны. Стенки трубопровода моделируются внутри блока, и температура стенки трубопровода в каждой зоне может быть различной. Давление и температура зависят от динамической сжимаемости жидкости и теплоемкости зоны жидкости.

Теплопередача между жидкостью и стенкой

Коэффициент конвективной теплопередачи между жидкостью и стенкой, αF, изменяется на зону в соответствии с числом Нуссельта:

αF=NukDH,

где:

  • Nu - номер зоны Нуссельта.

  • k среднюю теплопроводность жидкости.

  • D H - Hydraulic diameter трубопровода, эквивалентный диаметр некруглого сечения.

Число Нуссельта, используемое в коэффициенте теплопередачи, является большим из числа Нуссельта турбулентного и ламинарного потока.

Для турбулентных потоков в переохлажденных жидких или перегретых паровых зонах число Нуссельта вычисляется корреляцией Гнилинского:

Nu=f8(Re1000)PR1+12.7( f8)1/2(Pr2/31),

где:

  • Re - среднее Число Рейнольдса зоны.

  • Pr - среднее значение Prandtl в зоне.

  • f - коэффициент трения Дарси, вычисленный из корреляции Haaland:

    1f=1.8журнал[(εDH3.7)1.11+6.9Re],

    где ε - стенка Internal surface absolute roughness.

Для турбулентных потоков в зоне смеси пар-жидкости число Нуссельта вычисляется корреляцией Каваллини-Цекхина:

Nu=являютсяSLbPRSLc{[(ρSLρSV1)xOut+1]1+b[(ρSLρSV1)xIn+1]1+b}(1+b)(ρSLρSV1)(xOutxIn).

Где:

  • Re SL является числом Рейнольдса насыщенной жидкости.

  • Pr SL является числом Прандтля насыщенной жидкости.

  • ρ SL является плотностью насыщенной жидкости.

  • ρ SV является плотностью насыщенного пара.

  • a = 0,05, b = 0,8 и c = 0,33.

Когда ребра моделируются на внутренней поверхности трубопровода, коэффициент теплопередачи равен:

αF=NukDH(1+ηIntsInt),

где:

  • η Int - это Internal fin efficiency.

  • s Int - это Ratio of internal fins surface area to no-fin surface area.

Для ламинарных течений число Нуссельта задается параметром Laminar flow Nusselt number.

Специфическая энтальпия

Скорость теплопередачи от жидкости основана на изменении конкретной энтальпии в каждой зоне:

Q=m˙Q(ΔhL+ΔhM+ΔhV),

где m˙Q - массовый расход жидкости теплопередачи. Это входное отверстие трубопровода массового расхода жидкости, или m˙A или m˙B, в зависимости от направления потока жидкости.

В жидкой и паровой зонах изменение специфической энтальпии определяется как:

Δh=cp(THTI)[1exp(zSWm˙Qcp(αF1+αE1))],

где:

  • cp - удельная теплота жидкости или пара.

  • TH - температура окружающей среды.

  • TI - температура входного отверстия жидкости.

  • z - фракция длины зоны жидкости.

  • αE - коэффициент теплопередачи между стенкой и окружением.

  • SW - площадь поверхности стенки:

    SW=4ADHL,

    где:

    • 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 означает отсутствие ребер на стенке трубопровода.

В зоне смеси пар-жидкость изменение специфической энтальпии рассчитывается как:

Δh=(THTS)zSWm˙Q(αF1+αE1),

где TS - температура насыщения жидкости. Принято, что смесь пар-жидкость всегда находится при этой температуре.

Скорость теплопередачи

Общая теплопередача между жидкостью и стенкой трубопровода является суммой теплопередачи в каждой жидкой фазе:

QF=QF,L+QF,V+QF,M.

Скорость теплопередачи между жидкостью и трубопроводом в жидкой зоне:

QF,L=m˙Qcp,L[TW,Lминута(TI,TS)][1exp(zLSWαLm˙Qcp,L)].

где TW,L - температура стенки, окружающей жидкую зону.

Скорость теплопередачи между жидкостью и трубопроводом в зоне смеси:

QF,M=(THTSat)zMSWαM.

Скорость теплопередачи между жидкостью и трубопроводом в паровой зоне:

QF,V=m˙Qcp,V[TW,Vминута(TI,TSat)][1exp(zVSWαVm˙Qcp,V)],

где TW,V - температура стенки, окружающей паровую зону.

Теплопередача между Стенкой и Окружением

Если стенка трубопровода имеет конечную толщину, коэффициент теплопередачи между стеной и окружением, αE, определяется:

1αE=1αW+1αExt(1+ηExtsExt),

где α W - коэффициент теплопередачи, обусловленный проводимостью через стенку:

αW=kWDHln(1+tWDH),

и где:

  • 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.

Скорость теплопередачи

Скорость теплопередачи между каждой зоной стенки и окружения:

QH,zone=(THTW)zSWαE.

Общая теплопередача между стенкой и окружением составляет:

QH=QH,L+QH,V+QH,M.

Управляющие дифференциальные уравнения

Скорость теплопередачи зависит от тепловой массы стенки, CW:

CW=cp,WρWSW(tW+tW2DH),

где:

  • c p, W является Wall specific heat.

  • ρW является Wall density.

Руководящими уравнениями для теплопередачи между жидкостью и внешним окружением являются, для жидкой зоны:

QH,LQF,L=CW[zLdTW,Ldt+макс.(dzLdt,0)(TW,LTW,M)],

для зоны смеси:

QH,MQF,M=CW[zMdTW,Mdt+минута(dzLdt,0)(TW,LTW,M)+минута(dzVdt,0)(TW,VTW,M)],

и для паровой зоны:

QH,VQF,V=CW[zVdTW,Vdt+макс.(dzVdt,0)(TW,VTW,M)].

Баланс импульса

Перепад давления по трубе состоит из двух факторов: изменения давления из-за изменений плотности и изменения давления из-за трения на стенках трубопровода.

Для турбулентных потоков, когда число Рейнольдса выше Turbulent flow lower Reynolds number limit, падение давления вычисляется с точки зрения коэффициента трения Дарси. Перепад давления между портом A и внутренним узлом I:

pApЯ=(1ρI1ρA*)(m˙AS)2+fAm˙A|m˙A|2ρIDHS2(L+LAdd2),

где:

  • ρI - плотность жидкости во внутреннем узле I.

  • ρA* - плотность жидкости в порте A. Это то же самое, что и ρA, когда поток является установившимся; когда поток является переходным, он вычисляется из внутреннего состояния жидкости с помощью адиабатического выражения:

    uA*+pAρA*+12(m˙AρA*S)2=h+12(m˙AρS)2,

    где:

    • h - средняя специфическая энтальпия, h=hLzL+hVzV+hMzM.

    • ρ - средняя плотность, ρ=ρLzL+ρMzM+ρVzV.

    Это связано с тем, что расчет теплопередачи происходит на внутреннем узле I.

  • m˙A - массовый расход жидкости через порт A.

  • L является Pipe length.

  • L Add является Aggregate equivalent length of local resistances, которая является эквивалентной длиной трубки, которая вводит такое же количество потерь, как и сумма потерь из-за других локальных сопротивлений.

Обратите внимание, что коэффициент трения Дарси зависит от числа Рейнольдса и вычисляется в обоих портах.

Перепад давления между портом B и внутренним узлом I:

pBpЯ=(1ρI1ρB*)(m˙BS)2+fBm˙B|m˙B|2ρIDHS2(L+LAdd2).

где:

  • ρB* - плотность жидкости в порте B. Это то же самое, что и ρB, когда поток является установившимся; когда поток является переходным, он вычисляется из внутреннего состояния жидкости с помощью адиабатического выражения:

    uB*+pBρB*+12(m˙BρB*S)2=h+12(m˙BρS)2.

  • m˙B - массовый расход жидкости через порт B.

Для ламинарных течений, когда число Рейнольдса ниже Laminar flow upper Reynolds number limit, падение давления на трение вычисляется в терминах Laminar friction constant for Darcy friction factor, λ. Перепад давления между портом A и внутренним узлом I:

pApЯ=(1ρI1ρA*)(m˙AS)2+λμm˙A2ρIDH2S(L+LAdd2),

where - средняя динамическая вязкость жидкости:

μ=μLzL+μMzM+μVzV.

Перепад давления между портом B и внутренним узлом I:

pBpЯ=(1ρI1ρB*)(m˙BS)2+λμm˙B2ρIDH2S(L+LAdd2).

Для переходных потоков перепадом давления из-за вязкого трения является сглаженная смесь между значениями для ламинарного и турбулентного падения давления.

Баланс массы

Общая скорость накопления массы определяется как:

dMdt=m˙A+m˙B,

где M - общая масса жидкости в трубопроводе. Что касается зон жидкости, скорость накопления массы является функцией изменения плотности, ρ, относительно давления, p и удельной внутренней энергии жидкости, u:

dMdt=[(dρdp)udpdt+(dρdu)pduoutdt+ρLdzLdt+ρMdzMdt+ρVdzVdt]V,

где:

  • uout является специфической внутренней энергией после того, как произошла вся теплопередача.

  • V - общий объем жидкости или объем трубопровода.

Энергетический баланс

Уравнение энергетического разговора:

Mduoutdt+(m˙A+m˙B)uout=ϕA+ϕB+QF,

где:

  • ϕ A является расходом энергии в порту A.

  • ϕ B - расход энергии в порту B.

  • Q F - скорость теплопередачи между жидкостью и стенкой.

Допущения и ограничения

  • Стенки трубопровода идеально жесткие.

  • Поток полностью развит. Потери на трение и теплопередача не включают входные эффекты.

  • Инерция жидкости незначительна.

  • Эффект тяжести незначителен.

  • Когда давление выше критического давления жидкости, большие значения свойств тепловой жидкости (таких как число Прандтля, теплопроводность и удельное тепло) могут не точно отражать теплообмен в трубопроводе.

Порты

Выход

расширить все

Вектор с долями длин жидкой, смешанной фазы и паровых зон в трубопроводе.

Сохранение

расширить все

Открытие, посредством которого двухфазная жидкость течет внутрь или из трубопровода. Порты A и B могут каждый функционировать как вход или выход. Теплопроводность допускается между двухфазными портами жидкости и жидкостью, внутренней к трубопроводу (хотя ее влияние обычно актуально только при почти нулевых скоростях потока жидкости).

Открытие, посредством которого двухфазная жидкость течет внутрь или из трубопровода. Порты A и B могут каждый функционировать как вход или выход. Теплопроводность допускается между двухфазными портами жидкости и жидкостью, внутренней к трубопроводу (хотя ее влияние обычно актуально только при почти нулевых скоростях потока жидкости).

Тепловое граничное условие на внешней поверхности стенки трубопровода. Используйте этот порт для захвата теплообмена различных видов - например, проводящего, конвективного или радиационного - между стенкой трубопровода и окружением. Теплообмен между внутренней поверхностью стенки и жидкостью захватывается непосредственно в блоке.

Параметры

расширить все

Геометрия

Расстояние между портами трубопровода. Каждая из жидких, двухфазных и паровых зон составляет часть этого расстояния. Зональные дроби могут изменяться, но их совокупная длина, совпадающая с расстоянием трубопровода, является фиксированной.

Внутренняя область трубопровода, перпендикулярная направлению потока. Поперечное сечение трубопровода принято постоянным на протяжении всей ее длины.

Отношение площади открытия поперечного сечения трубопровода к периметру этой площади. Трубопровод не должен быть цилиндрическим, а его поперечное сечение имеет любую форму. Этот параметр задает диаметр, который общее сечение имело бы, если бы оно было округлым.

Вязкое трение

Комбинированная длина всех локальных сопротивлений, присутствующих в трубопроводе. Локальные сопротивления включают повороты, подборы кривой, якоря, а также входные и выходные отверстия трубопровода. Эффект локальных сопротивлений состоит в увеличении эффективной длины трубопровода. Эта длина добавляется к геометрической длине трубопровода для вычисления трения.

Средняя глубина всех поверхностных дефектов на внутренней поверхности трубопровода. Поверхностные дефекты влияют на падение давления через трубопровод в турбулентном режиме течения.

Число Рейнольдса, выше которого поток трубопровода начинает переход от ламинарного к турбулентному. Это значение является максимальным значением числа Рейнольдса, соответствующим полностью развитому ламинарному течению.

Число Рейнольдса, ниже которого поток трубопровода начинает переходить от турбулентного к ламинарному. Это значение является минимальным числом Рейнольдса, соответствующим полностью развивающемуся турбулентному потоку.

Безразмерный коэффициент, используемый для захвата эффектов геометрии поперечного сечения на вязкие потери на трение, происходящие в ламинарном режиме течения жидкости. Типичные значения 64 для круглого сечения, 57 для квадратного сечения и 62 для прямоугольного сечения с соотношением сторон 2.

Теплопередача

Коэффициент теплопередачи для теплообмена между окружением (в порту H) и внешней поверхностью стенки трубопровода. Если для этого параметра задано значение inf, тогда тепловое сопротивление между окружением и стенкой принято равным нулю. Затем стенка трубопровода имеет равномерную температуру, равную температуре, связанной с портом H.

Средняя толщина материала стенки трубопровода. Если для этого значения задано значение 0, тогда и тепловое сопротивление, обусловленное проводимостью через стенку трубопровода, и тепловое хранение, обусловленное тепловой массой стенки трубопровода, приняты незначительными.

Теплопроводность материала стенки трубопровода. Если для этого параметра задано значение inf, тогда тепловое сопротивление от проводимости через стенку трубопровода принято незначительным.

Теплоемкость на единицу массы материала стенки трубопровода. Если для этого параметра задано значение 0, тогда тепловое хранение из-за тепловой массы трубопровода принято незначительным.

Массовая плотность материала стенки трубопровода. Если для этого параметра задано значение 0, тогда тепловое хранение из-за тепловой массы стенки трубопровода принято незначительным.

Отношение общей площади теплопередающей поверхности ребер на внешней стороне стенки трубопровода к площади поверхности стенки трубопровода без каких-либо ребер. Наличие ребер служит для улучшения конвективной теплопередачи между стенкой трубопровода и окружением.

Отношение фактической скорости теплообмена между внешними ребрами и окружением к ее идеальному значению (если ребра полностью удерживались при температуре стенки трубопровода). Этот параметр является функцией геометрии плавника.

Отношение общей площади теплопередающей поверхности ребер на внутренней стороне стенки трубопровода к площади поверхности стенки трубопровода без каких-либо ребер. Наличие ребер служит для улучшения конвективной теплопередачи между стенкой трубопровода и жидкостью.

Отношение фактической скорости теплообмена между внутренними ребрами и окружением к ее идеальному значению (если ребра полностью удерживались при температуре стенки трубопровода). Этот параметр является функцией геометрии плавника.

Отношение конвективного и проводящего теплообмена в ламинарном режиме течения жидкости. Этот параметр является функцией сечения канала геометрии. Типичные значения 3.66 для круглого сечения, 2.98 для квадратного сечения и 3.39 для прямоугольного сечения с соотношением сторон 2.

Вкладка Эффекты и начальные условия

Термодинамическая переменная, в терминах которой можно задать начальные условия компонента.

Давление в трубопроводе в начале симуляции, заданное относительно нуля.

Температура жидкости в трубопроводе в начале симуляции. Этот параметр может быть скаляром или вектором с 2 элементами. Если это скаляр, начальная температура принимается равномерной по всему трубопроводу. Если это вектор, начальная температура принимается линейно изменяющейся между портами. Первый векторный элемент задает начальную температуру во входном отверстии и втором векторном элементе в выходном отверстии.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Temperature.

Качество пара, или массовая доля пара, в трубопроводе в начале симуляции. Этот параметр может быть скаляром или вектором с 2 элементами. Если это скаляр, начальное качество пара принято равномерным по всему трубопроводу. Если это вектор, начальное качество пара принимается линейно изменяющимся между портами. Первый векторный элемент дает начальное качество пара во входном отверстии и втором векторном элементе в выходном отверстии.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Vapor quality.

Паровая пустая фракция, или паровая объемная фракция, в трубопроводе в начале симуляции. Этот параметр может быть скаляром или вектором с 2 элементами. Если это скаляр, начальная доля паров в пустоте принимается равномерной по всему трубопроводу. Если это вектор, начальная паровая пустая фракция изменяется линейно между портами. Первый векторный элемент дает начальную долю паров в входном отверстии и второй векторный элемент в выходном отверстии..

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Vapor void fraction.

Специфическая энтальпия жидкости в трубопроводе в начале симуляции. Этот параметр может быть скаляром или вектором с 2 элементами. Если это скаляр, начальная специфическая энтальпия принята равномерной по всему трубопроводу. Если это вектор, начальная специфическая энтальпия принимается линейно изменяющейся между портами. Первый векторный элемент дает начальную специфическую энтальпию во входном отверстии и второй векторный элемент в выходном отверстии.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Specific enthalpy.

Удельная внутренняя энергия жидкости в трубопроводе в начале симуляции. Этот параметр может быть скаляром или вектором с 2 элементами. Если это скаляр, начальная удельная внутренняя энергия принята равномерной по всему трубопроводу. Если это вектор, начальная удельная внутренняя энергия изменяется линейно между портами. Первый векторный элемент дает начальную удельную внутреннюю энергию во входном отверстии и втором векторном элементе в выходном отверстии.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Specific internal energy.

Ссылки

[1] Белый, F.M., Механика жидкости, 7th Ред, раздел 6.8. Макгроу-Хилл, 2011.

[2] Qengel, Y.A., Теплопередача и массопередача - практический подход, 3rd Ред, раздел 8.5. Макгроу-Хилл, 2007.

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ Simulink ®

.
Введенный в R2018b
Для просмотра документации необходимо авторизоваться на сайте