Pipe (MA)

Жесткий трубопровод для потока сырого воздуха

  • Библиотека:
  • Simscape/Библиотека фундаментов/Сырой воздух/Элементы

  • Pipe (MA) block

Описание

Блок Pipe (MA) моделирует динамику потока трубопровода во влажной воздушной сети из-за вязких потерь на трение и конвективной теплопередачи со стенкой трубопровода. Трубопровод содержит постоянный объем сырого воздуха. Давление и температура развиваются исходя из сжимаемости и теплоемкости этого объема сырого воздуха. Жидкая вода конденсируется из объема сырого воздуха, когда она достигает насыщения. Дроссельный поток происходит, когда выход достигает звукового условия.

Внимание

Поток воздуха через этот блок может дросселироваться. Если блок Mass Flow Rate Source (MA) или блок Controlled Mass Flow Rate Source (MA), соединенный с блоком Pipe (MA), задает большой массовый расход жидкости, чем возможный подавленный массовый расход жидкости, симуляция генерирует ошибку. Для получения дополнительной информации см. Раздел «Подавленный поток».

В блочных уравнениях используются эти символы. Индексы a, w, и g указать свойства сухого воздуха, водяного пара и прослеживаемого газа, соответственно. Подстрочный ws указывает водяной пар при насыщении. Индексы A, B, H, и S указать соответствующий порт. Подстрочный I указывает свойства внутреннего объема сырого воздуха.

m˙Массовый расход жидкости
ΦЭнергетическая скорость потока жидкости
QТепловая скорость потока жидкости
pДавление
ρПлотность
RУдельная газовая константа
VОбъем сырого воздуха в трубопроводе
c vУдельное тепло при постоянном объеме
c pУдельное тепло при постоянном давлении
hСпецифическая энтальпия
uУдельная внутренняя энергия
xМассовая доля (x w - удельная влажность, что является другим термином для массовой фракции водяного пара)
yМольная дробь
φОтносительная влажность
rКоэффициент влажности
TТемпература
tВремя

Массовый и энергетический баланс

Чистые скорости потока жидкости в объем сырого воздуха внутри трубопровода

m˙net=m˙A+m˙Bm˙condense+m˙wS+m˙gSΦnet=ΦA+ΦB+QHΦcondense+ΦSm˙w,net=m˙wA+m˙wBm˙condense+m˙wSm˙g,net=m˙gA+m˙gB+m˙gS

где:

  • m˙конденсация - это скорость конденсации.

  • Φ конденсация - это скорость потерь энергии от конденсированной воды.

  • Φ S - это скорость энергии, добавляемой источниками влаги и следового газа .m˙wS и m˙gS - массовые расходы жидкости воды и газа, соответственно, через порт S. Значения m˙wS, m˙gS, и Φ S определяются источниками влаги и следового газа, соединенными с портом < reservedrangesplaceholder0 > трубопровода.

Сохранение массы водяного пара относится к массовому расходу жидкости водяного пара с динамикой уровня влаги во внутреннем объеме сырого воздуха:

dxwIdtρIV+xwIm˙net=m˙w,net

Точно так же сохранение массы следового газа связывает массовый расход жидкости следового газа с динамикой уровня следового газа во внутреннем объеме сырого воздуха:

dxgIdtρIV+xgIm˙net=m˙g,net

Сохранение массы смеси относится к массовому расходу жидкости динамике давления, температуры и массовой доли внутреннего объема сырого воздуха:

(1pIdpIdt1TIdTIdt)ρIV+RaRwRI(m˙w,netxwm˙net)+RaRgRI(m˙g,netxgm˙net)=m˙net

Наконец, энергосбережение связывает энергетическую скорость потока жидкости с динамикой давления, температуры и массовой доли внутреннего объема сырого воздуха:

ρIcvIVdTIdt+(uwIuaI)(m˙w,netxwm˙net)+(ugIuaI)(m˙g,netxgm˙net)+uIm˙net=Φnet

Уравнение состояния связывает плотность смеси с давлением и температурой:

pI=ρIRITI

Удельная газовая константа смеси

RI=xaIRa+xwIRw+xgIRg

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

Баланс импульса для каждой половины трубопровода моделирует перепад давления из-за потока импульса и вязкого трения:

pApI=(m˙AS)2(TIpITApA)RI+ΔpAIpBpI=(m˙BS)2(TIpITBpB)RI+ΔpBI

где:

  • p - давление в порте A, порте B или внутреннем узле I, как обозначено нижним индексом.

  • ρ - плотность порта A, порта B или внутреннего узла I, как обозначено нижним индексом.

  • S - площадь поперечного сечения трубопровода.

  • Δp AI и Δp BI - падения давления из-за вязкого трения.

Падения давления из-за вязкого трения, Δp AI и Δp BI, зависят от режима течения. Числа Рейнольдса для каждой половины трубопровода заданы как:

ReA=|m˙A|DhSμIReB=|m˙B|DhSμI

где:

  • D h - гидравлический диаметр трубопровода .

  • μ I является динамической вязкостью во внутреннем узле.

Если число Рейнольдса меньше значения параметра Laminar flow upper Reynolds number limit, то поток находится в ламинарном режиме течения жидкости. Если число Рейнольдса больше значения параметра Turbulent flow lower Reynolds number limit, то поток находится в турбулентном режиме течения.

В ламинарном режиме течения жидкости падения давления из-за вязкого трения:

ΔpAIlam=fshapem˙AμI2ρIDh2SL+Leqv2ΔpBIlam=fshapem˙BμI2ρIDh2SL+Leqv2

где:

  • f форма является значением параметра Shape factor for laminar flow viscous friction.

  • L eqv является значением параметра Aggregate equivalent length of local resistances.

В турбулентном режиме течения падения давления из-за вязкого трения:

ΔpAItur=fDarcyAm˙A|m˙A|2ρIDhS2L+Leqv2ΔpBItur=fDarcyBm˙B|m˙B|2ρIDhS2L+Leqv2

где f Darcy - коэффициент трения Дарси в порту A или B, как обозначено нижним индексом.

Коэффициенты трения Дарси вычисляются из корреляции Haaland:

fDarcyA=[1.8log(6.9ReA+(εrough3.7Dh)1.11)]2fDarcyB=[1.8log(6.9ReB+(εrough3.7Dh)1.11)]2

где ε rough - это значение параметра Internal surface absolute roughness.

Когда число Рейнольдса находится между значениями параметров Laminar flow upper Reynolds number limit и Turbulent flow lower Reynolds number limit, поток находится на переходном этапе между ламинарным течением и турбулентным потоком. Падения давления из-за вязкого трения во время переходной области следуют за плавной связью между таковыми в ламинарном режиме течения жидкости и таковыми в турбулентном режиме течения.

Теплообмен со стенкой трубопровода через порт H добавляется к энергии объема сырого воздуха, представленного внутренним узлом, посредством уравнения энергосбережения (см. «Масса и энергетический баланс»). Поэтому балансы импульса для каждой половины трубопровода, между портом A и внутренним узлом и между портом B и внутренним узлом, приняты адиабатическими процессами. Адиабатические зависимости:

hAhI=(RIm˙AS)2[(TIpI)2(TApA)2]hBhI=(RIm˙BS)2[(TIpI)2(TBpB)2]

где h - специфическая энтальпия в порте A, порте B или внутреннем узле I, как обозначено нижним индексом.

Конвективная теплопередача

Уравнение конвективной теплопередачи между стенкой трубопровода и внутренним объемом сырого воздуха является:

QH=Qconv+kISsurfDh(THTI)

S surf - площадь поверхности трубопровода, S surf = 4 S L/ D h. Принимая экспоненциальное распределение температуры по трубопроводу, конвективная теплопередача

Qconv=|m˙avg|cpavg(THTin)(1ehcoeffSsurf|m˙avg|cpavg)

где:

  • T - температура на входе в зависимости от направления потока.

  • m˙avg=(m˙Am˙B)/2 - средний массовый расход жидкости от порта A до порта B.

  • cpavg - удельное тепло, оцениваемое при средней температуре.

Коэффициент теплопередачи, h коэффициент, зависит от числа Нуссельта:

hcoeff=NukavgDh

где k avg - теплопроводность, оцениваемая при средней температуре. Число Нуссельта зависит от режима течения. Число Нуссельта в ламинарном режиме течения жидкости постоянно и равно значению параметра Nusselt number for laminar flow heat transfer. Число Нуссельта в турбулентном режиме течения вычисляется из корреляции Гнилинского:

Nutur=fDarcy8(Reavg1000)Pravg1+12.7fDarcy8(Pravg2/31)

где Pr avg - число Прандтля, оцениваемое при средней температуре. Среднее число Рейнольдса

Reavg=|m˙avg|DhSμavg

где μ avg - динамическая вязкость, оцениваемая при средней температуре. Когда среднее число Рейнольдса находится между Laminar flow upper Reynolds number limit и Turbulent flow lower Reynolds number limit значениями параметров, число Нуссельта следует плавному переходу между ламинарным и турбулентным значениями числа Нуссельта.

Насыщение и конденсация

Когда объем сырого воздуха достигает насыщения, может происходить конденсация. Удельная влажность при насыщении

xwsI=φwsRIRwpwsIpI

где:

  • φ ws - относительная влажность при насыщении (обычно 1).

  • p wsI является давлением насыщения водяного пара, оцениваемым в T I.

Скорость конденсации составляет

m˙condense={0,если xwIxwsIxwIxwsIτcondenseρIV,если xwI>xwsI

где τ condense - это значение параметра Condensation time constant.

Конденсированная вода вычитается из объема сырого воздуха, как показано в уравнениях сохранения. Энергия, связанная с конденсированной водой,

Φcondense=m˙condense(hwIΔhvapI)

где Δh vapI - специфическая энтальпия испарения, оцениваемая в T I.

Другие количества влаги и следовых газов связаны друг с другом следующим образом:

φwI=ywIpIpwsIywI=xwIRwRIrwI=xwI1xwIygI=xgIRgRIxaI+xwI+xgI=1

Дроссельный поток

Отмененное давление в порту A или B является значением соответствующей переменной Across в этом порту:

pAunchoked=A.ppBunchoked=B.p

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

Дросселирование может происходить на выходе трубопровода, но не на входе трубопровода. Поэтому, если порт A является входным отверстием, то p A = A.p. Если port A является выходным отверстием, то

pA={A.p,если A.ppAchokedpAchoked,если A.p<pAchoked 

Точно так же, если порт B является входным отверстием, то p B = B.p. Если port B является выходным отверстием, то

pB={B.p,если B.ppBchokedpBchoked,если B.p<pBchoked 

Подавленные давления в портах A и B получают из баланса импульса, принимая, что скорость выхода равна скорости звука:

pAchokedpI=pAchoked(pAchokedTIpITA1)cpAcvI+ΔpAIpBchokedpI=pBchoked(pBchokedTIpITB1)cpBcvI+ΔpBI

Переменные

Чтобы задать приоритет и начальные целевые значения для основных переменных до симуляции, используйте вкладку Variables в диалоговом окне блока (или раздел Variables в Property Inspector блоков). Для получения дополнительной информации смотрите Задать приоритет и Начальную цель для основных переменных и Начальные условия для блоков с конечным объемом влажного воздуха.

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

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

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

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

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

  • Этот блок не моделирует сверхзвуковой поток.

Порты

Выход

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

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

Выходной порт физического сигнала, который выводит сигнал вектора. Вектор содержит давление (в Pa), температуру (в K), уровень влаги и проследить измерения уровня газа внутри компонента. Используйте блок Measurement Selector (MA), чтобы распаковать этот векторный сигнал.

Сохранение

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

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

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

Тепловой порт сопоставлен с температурой стенки трубопровода. Блок включает конвективную теплопередачу между влажной воздушной смесью внутри трубопровода и стенкой трубопровода.

Подключите этот порт к S портов блока из библиотеки «Источники влаги и прослеживания газа», чтобы добавить или удалить влагу и проследить газ. Для получения дополнительной информации см. Использование источников влаги и прослеживания газа.

Зависимости

Этот порт видим, только если вы задаете значение параметра Moisture and trace gas source Controlled.

Параметры

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

Главный

Длина трубопровода по направлению потока.

Внутренняя область трубопровода, перпендикулярная направлению потока.

Диаметр эквивалентного цилиндрического трубопровода с той же площадью поперечного сечения.

Трение и теплопередача

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

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

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

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

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

Отношение конвективного к проводящему теплопередаче в ламинарном режиме течения жидкости. Его значение зависит от геометрии поперечного сечения трубопровода и тепловых граничных условий стенки трубопровода, таких как постоянная температура или постоянный тепловой поток. Типичное значение 3,66 для круглого сечения с постоянной температурой стенки [2].

Влага и прослеживание газа

Относительная влажность, выше которой происходит конденсация.

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

Этот параметр контролирует видимость портовых S и предоставляет следующие опции для моделирования уровня влаги и прослеживания газа внутри компонента:

  • None - Никакая влага или следовой газ не впрыскивается или не извлекается из блока. Порт S скрыт. Это значение по умолчанию.

  • Constant - Влага и следовой газ впрыскиваются в блок или извлекаются из него с постоянной скоростью. Те же параметры, что и в блоках Moisture Source (MA) и Trace Gas Source (MA), становятся доступными в Moisture and Trace Gas разделе интерфейса блока. Порт S скрыт.

  • Controlled - Влага и следовой газ впрыскиваются в блок или извлекаются из него с изменяющейся во времени скоростью. Порт S доступен. Соедините Controlled Moisture Source (MA) и Controlled Trace Gas Source (MA) блоки с этим портом.

Выберите, будет ли блок добавлять или удалять влагу в виде водяного пара или жидкой воды:

  • Vapor Энтальпия добавленной или удаленной влаги соответствует энтальпии водяного пара, которая больше, чем у жидкой воды.

  • Liquid Энтальпия добавленной или удаленной влаги соответствует энтальпии жидкой воды, которая меньше, чем у водяного пара.

Зависимости

Активируется, когда для параметра Moisture and trace gas source задано значение Constant.

Водяной пар массового расхода жидкости через блок. Положительное значение добавляет влагу в связанный объем сырого воздуха. Отрицательное значение извлекает влагу из этого объема.

Зависимости

Активируется, когда для параметра Moisture and trace gas source задано значение Constant.

Выберите метод спецификации температуры влаги:

  • Atmospheric temperature - Используйте температуру атмосферы, заданную блоком Moist Air Properties (MA), соединенным с цепью.

  • Specified temperature - Задайте значение при помощи параметра Temperature of added moisture.

Зависимости

Активируется, когда для параметра Moisture and trace gas source задано значение Constant.

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

Зависимости

Активируется, когда для параметра Added moisture temperature specification задано значение Specified temperature.

Проследите массовый расход жидкости газа через блок. Положительное значение добавляет следовой газ в связанный объем сырого воздуха. Отрицательное значение извлекает следовой газ из этого объема.

Зависимости

Активируется, когда для параметра Moisture and trace gas source задано значение Constant.

Выберите метод спецификации для температуры прослеживаемого газа:

  • Atmospheric temperature - Используйте температуру атмосферы, заданную блоком Moist Air Properties (MA), соединенным с цепью.

  • Specified temperature - Задайте значение при помощи параметра Temperature of added trace gas.

Зависимости

Активируется, когда для параметра Moisture and trace gas source задано значение Constant.

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

Зависимости

Активируется, когда для параметра Added trace gas temperature specification задано значение Specified temperature.

Примеры моделей

Ссылки

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

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

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

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

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