exponenta event banner

Труба (MA)

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

  • Библиотека:
  • Библиотека Simscape/Foundation/Влажный воздух/Элементы

  • Pipe (MA) block

Описание

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

Внимание

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

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

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

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

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

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

где:

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

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

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

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

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

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

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

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

(1pIdpIdt 1TIdTIdt) αIV + Ra RwRI (m˙w,net−xwm˙net) + Ra RgRI (m˙g,net−xgm˙net) =m˙net

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

(uwI uaI) (m˙w,net−xwm˙net) + (ugI uaI) (m˙g,net−xgm˙net) +uIm˙net=Φnet

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

pI = αIRITI

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

RI = xaIRa + xwIRw + xgIRg

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

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

pA pI = (m˙AS) 2⋅ (TIpI TApA) RI + ΔpAIpB pI = (m˙BS) 2⋅ (TIpI − TBpB) RI + ΔpBI

где:

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

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

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

  • ΔpAI и ΔpBI - потери давления из-за вязкого трения.

Потери давления из-за вязкого трения ΔpAI и ΔpBI зависят от режима потока. Числа Рейнольдса для каждой половины трубы определяются как:

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

где:

  • Dh - гидравлический диаметр трубы.

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

Если число Рейнольдса меньше значения параметра верхнего предела числа Рейнольдса потока Ламинара, то поток находится в режиме ламинарного потока. Если число Рейнольдса больше, чем значение параметра ограничения нижнего числа Рейнольдса турбулентного потока, то поток находится в режиме турбулентного потока.

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

ΔpAIlam=fshapem˙A⋅μI2ρI⋅Dh2⋅S⋅L+Leqv2ΔpBIlam=fshapem˙B⋅μI2ρI⋅Dh2⋅S⋅L+Leqv2

где:

  • fshape - значение коэффициента Shape для параметра вязкого трения ламинарного потока.

  • Leqv - значение параметра Агрегатная эквивалентная длина локальных сопротивлений.

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

ΔpAItur=fDarcyAm˙A⋅|m˙A|2ρI⋅Dh⋅S2⋅L+Leqv2ΔpBItur=fDarcyBm˙B⋅|m˙B|2ρI⋅Dh⋅S2⋅L+Leqv2

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

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

fDarcyA = [− 1 .8log (6 .9ReA + (εrough3.7Dh) 1,11)] 2fDarcyB = [1 .8log (6 .9ReB + (εrough3.7Dh) 1,11)] − 2

где αrough - значение параметра абсолютная шероховатость внутренней поверхности.

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

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

hA hI = (RIm˙AS) 2 [(TIpI) 2 (TApA) 2] hB hI = (RIm˙BS) 2 [(TIpI) 2 − (TBpB) 2]

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

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

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

QH = Qconv + kISsurfDh (TH − TI)

Ssurf - площадь поверхности трубы, Ssurf = 4SL/Dh. Предполагая экспоненциальное распределение температуры вдоль трубы, конвективная теплопередача равна

Qconv=|m˙avg'cpavg (TH олово) (1−ehcoeffSsurf'm˙avg'cpavg)

где:

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

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

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

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

hcoeff = NukavgDh

где kavg - теплопроводность, оцениваемая при средней температуре. Число Нуссельта зависит от режима потока. Число Нуссельта в режиме ламинарного потока является постоянным и равно значению числа Нуссельта для параметра теплопередачи ламинарного потока. Число Нуссельта в режиме турбулентного потока вычисляется по корреляции Гниелинского:

Нутур = fDarcy8 (Реавг − 1000) Pravg1 + 12 .7fDarcy8 (Pravg2/3 − 1)

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

Reavg=|m˙avg'DhSμavg

где мкг - динамическая вязкость, оцененная при средней температуре. Когда среднее число Рейнольдса находится между верхним пределом числа Рейнольдса потока Ламинара и нижним значением параметра предела числа Рейнольдса потока турбулента, число Нуссельта следует плавному переходу между ламинарным и турбулентным значениями числа Нуссельта.

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

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

xwsI = фwsRIRwpwsIpI

где:

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

  • pwsI - давление насыщения водяным паром, оцениваемое при TI.

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

m˙condense={0,if xwI≤xwsIxwI−xwsIτcondenseρIV,if xwI > xwsI

где startcondense - значение параметра постоянной времени конденсации.

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

Φcondense=m˙condense (hwI ΔhvapI)

где ΔhvapI - специфическая энтальпия испарения, оцененная при TI.

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

фwI = ywIpIpwsIywI = xwIRwRIrwI = xwI1 xwIygI = xgIRgRIxaI + xwI + xgI = 1

Подавленный поток

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

pAunchoked = A.ppBunchoked = B.p

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

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

pA = {A.p, если  A.p≥pAchokedpAchoked,if A.p <  pAchoked

Аналогично, если порт B является входом, то pB = B.p. Если порт B является выходом, то

pB = {B.p, если  B.p≥pBchokedpBchoked,if B.p <  pBchoked

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

pAchoked−pI=pAchoked (pAchokedTIpITA−1) cpAcvI +ΔpAIpBchoked−pI=pBchoked (pBchokedTIpITB−1) cpBcvI +ΔpBI

Переменные

Чтобы задать приоритет и начальные целевые значения для переменных блока перед моделированием, используйте вкладку «Переменные» в диалоговом окне блока (или раздел «Переменные» в Инспекторе свойств блока). Дополнительные сведения см. в разделе Установка приоритета и начальной цели для переменных блока и начальных условий для блоков с конечным объемом влажного воздуха.

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

  • Стенка трубы совершенно жесткая.

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

  • Влияние силы тяжести ничтожно мало.

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

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

Порты

Продукция

развернуть все

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

Порт вывода физического сигнала, который выводит векторный сигнал. Вектор содержит измерения давления (в Па), температуры (в К), уровня влаги и уровня следового газа внутри компонента. Для распаковки этого векторного сигнала используется блок выбора измерений (MA).

Сохранение

развернуть все

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

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

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

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

Зависимости

Этот порт отображается только в том случае, если для параметра Источник влаги и газа трассировки задано значение Controlled.

Параметры

развернуть все

Главный

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

Внутренняя площадь трубы по нормали к направлению потока.

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

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

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

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

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

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

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

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

Влага и следовой газ

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

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

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

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

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

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

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

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

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

Зависимости

Включено, если для параметра Источник влаги и газа трассировки установлено значение Constant.

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

Зависимости

Включено, если для параметра Источник влаги и газа трассировки установлено значение Constant.

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

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

  • Specified temperature - укажите значение с помощью параметра Температура добавленной влаги.

Зависимости

Включено, если для параметра Источник влаги и газа трассировки установлено значение Constant.

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

Зависимости

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

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

Зависимости

Включено, если для параметра Источник влаги и газа трассировки установлено значение Constant.

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

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

  • Specified temperature - укажите значение с помощью параметра Температура добавленного газа.

Зависимости

Включено, если для параметра Источник влаги и газа трассировки установлено значение Constant.

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

Зависимости

Активируется, если для параметра Спецификация добавленной температуры газа трассировки установлено значение Specified temperature.

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

Ссылки

[1] Уайт, Ф. М., Механика жидкости. 7-е изд., раздел 6.8. Макгроу-Хилл, 2011.

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

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

Создание кода C/C + +
Создайте код C и C++ с помощью Simulink ® Coder™

.
Представлен в R2018a