Жесткий трубопровод для потока сырого воздуха
Simscape/Библиотека фундаментов/Сырой воздух/Элементы
Блок Pipe (MA) моделирует динамику потока трубопровода во влажной воздушной сети из-за вязких потерь на трение и конвективной теплопередачи со стенкой трубопровода. Трубопровод содержит постоянный объем сырого воздуха. Давление и температура развиваются исходя из сжимаемости и теплоемкости этого объема сырого воздуха. Жидкая вода конденсируется из объема сырого воздуха, когда она достигает насыщения. Дроссельный поток происходит, когда выход достигает звукового условия.
Внимание
Поток воздуха через этот блок может дросселироваться. Если блок Mass Flow Rate Source (MA) или блок Controlled Mass Flow Rate Source (MA), соединенный с блоком Pipe (MA), задает большой массовый расход жидкости, чем возможный подавленный массовый расход жидкости, симуляция генерирует ошибку. Для получения дополнительной информации см. Раздел «Подавленный поток».
В блочных уравнениях используются эти символы. Индексы a
, w
, и g
указать свойства сухого воздуха, водяного пара и прослеживаемого газа, соответственно. Подстрочный ws
указывает водяной пар при насыщении. Индексы A
, B
, H
, и S
указать соответствующий порт. Подстрочный I
указывает свойства внутреннего объема сырого воздуха.
Массовый расход жидкости | |
Φ | Энергетическая скорость потока жидкости |
Q | Тепловая скорость потока жидкости |
p | Давление |
ρ | Плотность |
R | Удельная газовая константа |
V | Объем сырого воздуха в трубопроводе |
c v | Удельное тепло при постоянном объеме |
c p | Удельное тепло при постоянном давлении |
h | Специфическая энтальпия |
u | Удельная внутренняя энергия |
x | Массовая доля (x w - удельная влажность, что является другим термином для массовой фракции водяного пара) |
y | Мольная дробь |
φ | Относительная влажность |
r | Коэффициент влажности |
T | Температура |
t | Время |
Чистые скорости потока жидкости в объем сырого воздуха внутри трубопровода
где:
конденсация - это скорость конденсации.
Φ конденсация - это скорость потерь энергии от конденсированной воды.
Φ S - это скорость энергии, добавляемой источниками влаги и следового газа . и - массовые расходы жидкости воды и газа, соответственно, через порт S. Значения , , и Φ S определяются источниками влаги и следового газа, соединенными с портом < reservedrangesplaceholder0 > трубопровода.
Сохранение массы водяного пара относится к массовому расходу жидкости водяного пара с динамикой уровня влаги во внутреннем объеме сырого воздуха:
Точно так же сохранение массы следового газа связывает массовый расход жидкости следового газа с динамикой уровня следового газа во внутреннем объеме сырого воздуха:
Сохранение массы смеси относится к массовому расходу жидкости динамике давления, температуры и массовой доли внутреннего объема сырого воздуха:
Наконец, энергосбережение связывает энергетическую скорость потока жидкости с динамикой давления, температуры и массовой доли внутреннего объема сырого воздуха:
Уравнение состояния связывает плотность смеси с давлением и температурой:
Удельная газовая константа смеси
Баланс импульса для каждой половины трубопровода моделирует перепад давления из-за потока импульса и вязкого трения:
где:
p - давление в порте A, порте B или внутреннем узле I, как обозначено нижним индексом.
ρ - плотность порта A, порта B или внутреннего узла I, как обозначено нижним индексом.
S - площадь поперечного сечения трубопровода.
Δp AI и Δp BI - падения давления из-за вязкого трения.
Падения давления из-за вязкого трения, Δp AI и Δp BI, зависят от режима течения. Числа Рейнольдса для каждой половины трубопровода заданы как:
где:
D h - гидравлический диаметр трубопровода .
μ I является динамической вязкостью во внутреннем узле.
Если число Рейнольдса меньше значения параметра Laminar flow upper Reynolds number limit, то поток находится в ламинарном режиме течения жидкости. Если число Рейнольдса больше значения параметра Turbulent flow lower Reynolds number limit, то поток находится в турбулентном режиме течения.
В ламинарном режиме течения жидкости падения давления из-за вязкого трения:
где:
f форма является значением параметра Shape factor for laminar flow viscous friction.
L eqv является значением параметра Aggregate equivalent length of local resistances.
В турбулентном режиме течения падения давления из-за вязкого трения:
где f Darcy - коэффициент трения Дарси в порту A или B, как обозначено нижним индексом.
Коэффициенты трения Дарси вычисляются из корреляции Haaland:
где ε rough - это значение параметра Internal surface absolute roughness.
Когда число Рейнольдса находится между значениями параметров Laminar flow upper Reynolds number limit и Turbulent flow lower Reynolds number limit, поток находится на переходном этапе между ламинарным течением и турбулентным потоком. Падения давления из-за вязкого трения во время переходной области следуют за плавной связью между таковыми в ламинарном режиме течения жидкости и таковыми в турбулентном режиме течения.
Теплообмен со стенкой трубопровода через порт H добавляется к энергии объема сырого воздуха, представленного внутренним узлом, посредством уравнения энергосбережения (см. «Масса и энергетический баланс»). Поэтому балансы импульса для каждой половины трубопровода, между портом A и внутренним узлом и между портом B и внутренним узлом, приняты адиабатическими процессами. Адиабатические зависимости:
где h - специфическая энтальпия в порте A, порте B или внутреннем узле I, как обозначено нижним индексом.
Уравнение конвективной теплопередачи между стенкой трубопровода и внутренним объемом сырого воздуха является:
S surf - площадь поверхности трубопровода, S surf = 4 S L/ D h. Принимая экспоненциальное распределение температуры по трубопроводу, конвективная теплопередача
где:
T - температура на входе в зависимости от направления потока.
- средний массовый расход жидкости от порта A до порта B.
- удельное тепло, оцениваемое при средней температуре.
Коэффициент теплопередачи, h коэффициент, зависит от числа Нуссельта:
где k avg - теплопроводность, оцениваемая при средней температуре. Число Нуссельта зависит от режима течения. Число Нуссельта в ламинарном режиме течения жидкости постоянно и равно значению параметра Nusselt number for laminar flow heat transfer. Число Нуссельта в турбулентном режиме течения вычисляется из корреляции Гнилинского:
где Pr avg - число Прандтля, оцениваемое при средней температуре. Среднее число Рейнольдса
где μ avg - динамическая вязкость, оцениваемая при средней температуре. Когда среднее число Рейнольдса находится между Laminar flow upper Reynolds number limit и Turbulent flow lower Reynolds number limit значениями параметров, число Нуссельта следует плавному переходу между ламинарным и турбулентным значениями числа Нуссельта.
Когда объем сырого воздуха достигает насыщения, может происходить конденсация. Удельная влажность при насыщении
где:
φ ws - относительная влажность при насыщении (обычно 1).
p wsI является давлением насыщения водяного пара, оцениваемым в T I.
Скорость конденсации составляет
где τ condense - это значение параметра Condensation time constant.
Конденсированная вода вычитается из объема сырого воздуха, как показано в уравнениях сохранения. Энергия, связанная с конденсированной водой,
где Δh vapI - специфическая энтальпия испарения, оцениваемая в T I.
Другие количества влаги и следовых газов связаны друг с другом следующим образом:
Отмененное давление в порту A или B является значением соответствующей переменной Across в этом порту:
Однако переменные давления в порте, используемые в уравнениях баланса импульса, p A и p B, не обязательно совпадают с давлением между переменными A.p
и B.p
потому что выход трубопровода может дросселироваться. Дроссельный поток происходит, когда давление ниже по потоку достаточно низкое. На этой точке поток зависит только от условий на входном отверстии. Поэтому при дросселировании давление на выходе (p A или p B, в зависимости от того, что является выходным отверстием) не может уменьшаться дальше, даже если давление ниже по потоку, представленное A.p
или B.p
, продолжает уменьшаться.
Дросселирование может происходить на выходе трубопровода, но не на входе трубопровода. Поэтому, если порт A является входным отверстием, то p A = A.p
. Если port A является выходным отверстием, то
Точно так же, если порт B является входным отверстием, то p B = B.p
. Если port B является выходным отверстием, то
Подавленные давления в портах A и B получают из баланса импульса, принимая, что скорость выхода равна скорости звука:
Чтобы задать приоритет и начальные целевые значения для основных переменных до симуляции, используйте вкладку Variables в диалоговом окне блока (или раздел Variables в Property Inspector блоков). Для получения дополнительной информации смотрите Задать приоритет и Начальную цель для основных переменных и Начальные условия для блоков с конечным объемом влажного воздуха.
Стенки трубопровода идеально жесткие.
Поток полностью развит. Потери на трение и теплопередача не включают входные эффекты.
Эффект тяжести незначителен.
Инерция жидкости незначительна.
Этот блок не моделирует сверхзвуковой поток.
[1] Белый, F. M., Механика Жидкости. 7th Ed, раздел 6.8. Макгроу-Хилл, 2011.
[2] Cengel, Y. A., Теплопередача и массопередача - практический подход. 3-е Изд., Раздел 8.5. Макгроу-Хилл, 2007.