Жесткий трубопровод для потока влажного воздуха
Библиотека Simscape/Foundation/Влажный воздух/Элементы

Блок 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˙condense - скорость конденсации.
Dwfcondense - скорость потери энергии из конденсированной воды.
ФS - это скорость энергии, добавляемой источниками влаги и следового газа. и - массовые расходы воды и газа, соответственно, через порт S. Значения , и ФS определяются источниками влаги и следового газа, подключенными к порту S трубы.
Сохранение массы водяного пара связывает массовый расход водяного пара с динамикой уровня влаги во внутреннем объеме влажного воздуха:
Аналогично сохранение массы следового газа связывает массовый расход следового газа с динамикой уровня следового газа во внутреннем объеме влажного воздуха:
Сохранение массы смеси соотносит массовый расход смеси с динамикой давления, температуры и массовых долей внутреннего объема влажного воздуха:
m˙g,net−xgm˙net) =m˙net
Наконец, энергосбережение связывает расход энергии с динамикой давления, температуры и массовых долей внутреннего объема влажного воздуха:
+uIm˙net=Φnet
Уравнение состояния соотносит плотность смеси с давлением и температурой:
αIRITI
Удельная газовая константа смеси равна
xgIRg
Баланс импульса для каждой половины трубы моделирует падение давления из-за потока импульса и вязкого трения:
TIpI − TBpB) RI + ΔpBI
где:
p - давление в порте A, порте B или внутреннем узле I, как указано подстрочным индексом.
start- плотность в порте A, порте B или внутреннем узле I, как указано подстрочным индексом.
S - площадь поперечного сечения трубы.
ΔpAI и ΔpBI - потери давления из-за вязкого трения.
Потери давления из-за вязкого трения ΔpAI и ΔpBI зависят от режима потока. Числа Рейнольдса для каждой половины трубы определяются как:
где:
Dh - гидравлический диаметр трубы.
мкI - динамическая вязкость во внутреннем узле.
Если число Рейнольдса меньше значения параметра верхнего предела числа Рейнольдса потока Ламинара, то поток находится в режиме ламинарного потока. Если число Рейнольдса больше, чем значение параметра ограничения нижнего числа Рейнольдса турбулентного потока, то поток находится в режиме турбулентного потока.
В ламинарном режиме потока потери давления из-за вязкого трения составляют:
где:
fshape - значение коэффициента Shape для параметра вязкого трения ламинарного потока.
Leqv - значение параметра Агрегатная эквивалентная длина локальных сопротивлений.
В турбулентном режиме потока потери давления из-за вязкого трения составляют:
где fDarcy - коэффициент трения Дарси в портах A или B, как указано подстрочным индексом.
Коэффициенты трения Дарси вычисляются на основе корреляции Хааланда:
εrough3.7Dh) 1,11)] − 2
где αrough - значение параметра абсолютная шероховатость внутренней поверхности.
Когда число Рейнольдса находится между верхним пределом числа Рейнольдса потока Ламинара и нижними значениями предельных параметров числа Рейнольдса потока турбулентного потока, поток переходит между ламинарным потоком и турбулентным потоком. Потери давления из-за вязкого трения во время переходной области следуют за плавным соединением между теми, которые находятся в режиме ламинарного потока, и теми, которые находятся в режиме турбулентного потока.
Тепло, обмениваемое со стенкой трубы через канал Н, добавляется к энергии объема влажного воздуха, представленного внутренним узлом, посредством уравнения энергосбережения (см. массовый и энергетический баланс). Поэтому предполагается, что балансы импульса для каждой половины трубы, между портом A и внутренним узлом и между портом B и внутренним узлом, являются адиабатическими процессами. Адиабатические отношения:
[(TIpI) 2 − (TBpB) 2]
где h - специфическая энтальпия в порте A, порте B или внутреннем узле I, как указано подстрочным индексом.
Уравнение конвективной теплопередачи между стенкой трубы и внутренним объемом влажного воздуха составляет:
TH − TI)
Ssurf - площадь поверхности трубы, Ssurf = 4SL/Dh. Предполагая экспоненциальное распределение температуры вдоль трубы, конвективная теплопередача равна
1−ehcoeffSsurf'm˙avg'cpavg)
где:
Олово - температура на входе в зависимости от направления потока.
/2 - средний массовый расход от порта A к порту B.
- удельная теплота, оцениваемая при средней температуре.
Коэффициент теплопередачи, hcoeff, зависит от числа Нуссельта:
NukavgDh
где kavg - теплопроводность, оцениваемая при средней температуре. Число Нуссельта зависит от режима потока. Число Нуссельта в режиме ламинарного потока является постоянным и равно значению числа Нуссельта для параметра теплопередачи ламинарного потока. Число Нуссельта в режиме турбулентного потока вычисляется по корреляции Гниелинского:
Pravg2/3 − 1)
где Pravg - число Прандтля, оцениваемое при средней температуре. Среднее число Рейнольдса -
где мкг - динамическая вязкость, оцененная при средней температуре. Когда среднее число Рейнольдса находится между верхним пределом числа Рейнольдса потока Ламинара и нижним значением параметра предела числа Рейнольдса потока турбулента, число Нуссельта следует плавному переходу между ламинарным и турбулентным значениями числа Нуссельта.
Когда объем влажного воздуха достигает насыщения, может происходить конденсация. Удельная влажность при насыщении составляет
фwsRIRwpwsIpI
где:
αws - относительная влажность при насыщении (обычно 1).
pwsI - давление насыщения водяным паром, оцениваемое при TI.
Скорость конденсации составляет
xwsI
где startcondense - значение параметра постоянной времени конденсации.
Конденсированная вода вычитается из объема влажного воздуха, как показано в уравнениях сохранения. Энергия, связанная со сконденсированной водой, составляет
ΔhvapI)
где ΔhvapI - специфическая энтальпия испарения, оцененная при TI.
Другие количества влаги и следовых газов связаны друг с другом следующим образом:
+ xwI + xgI = 1
Отключенное давление в порте A или B - это значение соответствующей переменной Across в этом порте:
B.p
Однако переменные давления порта, используемые в уравнениях баланса импульса, pA и pB, не обязательно совпадают с давлением между переменными. A.p и B.p потому что выпускное отверстие трубы может задохнуться. Подавленный поток возникает, когда давление ниже по потоку достаточно низкое. В этот момент поток зависит только от условий на входе. Следовательно, при дросселировании давление на выходе (pA или pB, в зависимости от того, что является выходом) не может дополнительно уменьшаться, даже если давление ниже по потоку, представленное A.p или B.p, продолжает уменьшаться.
Дросселирование может происходить на выходе трубы, но не на входе трубы. Следовательно, если порт A является входным, то pA = A.p. Если порт A является выходом, то
pAchoked
Аналогично, если порт B является входом, то pB = B.p. Если порт B является выходом, то
pBchoked
Подавленные давления в портах А и В получают из баланса импульса, предполагая, что выходная скорость равна скорости звука:
cpBcvI +ΔpBI
Чтобы задать приоритет и начальные целевые значения для переменных блока перед моделированием, используйте вкладку «Переменные» в диалоговом окне блока (или раздел «Переменные» в Инспекторе свойств блока). Дополнительные сведения см. в разделе Установка приоритета и начальной цели для переменных блока и начальных условий для блоков с конечным объемом влажного воздуха.
Стенка трубы совершенно жесткая.
Поток полностью развит. Потери на трение и теплопередача не включают входные эффекты.
Влияние силы тяжести ничтожно мало.
Инерция жидкости ничтожна.
Этот блок не моделирует сверхзвуковой поток.
[1] Уайт, Ф. М., Механика жидкости. 7-е изд., раздел 6.8. Макгроу-Хилл, 2011.
[2] Cengel, Y. A., Тепло и массообмен - практический подход. 3-е изд., раздел 8.5. Макгроу-Хилл, 2007.