Сопротивление потока (2P)

Общее сопротивление в двухфазном жидком ответвлении

  • Библиотека:
  • Simscape / Библиотека Основы / Двухфазная Жидкость / Элементы

Описание

Блок Flow Resistance (2P) моделирует общий перепад давления в двухфазном жидком сетевом ответвлении. Перепад давления пропорционален квадрату массовой скорости потока жидкости и к плотности двухфазной жидкости. Коэффициент пропорциональности определяется от номинальных условий работы, заданных в диалоговом окне блока.

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

Массовый баланс

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

m˙A+m˙B=0,

где m˙A и m˙B заданы как массовые скорости потока жидкости на компонент через порты A и B, соответственно.

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

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

ϕA+ϕB=0,

где ϕ A и ϕ B является энергетическими скоростями потока жидкости в сопротивление потока через порты A и B.

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

Соответствующие внешние силы на жидкости включают тех из-за давления в портах и тех из-за вязкого трения в стенах компонента. Сила тяжести проигнорирована, как другие массовые силы. При выражении фрикционных сил с точки зрения коэффициента потерь ξ приводит к полуэмпирическому выражению:

Δp=ξνm˙22S2,

где:

  • Δp является перепадом давления от порта к порту B — то есть, p + p B.

  • ξ является коэффициентом потерь.

  • ν является определенным объемом, инверсией массовой плотности ρ — то есть, 1/ ρ.

  • S является областью потока.

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

Δp=ξνm˙|m˙|2S2,

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

Δp=ξνm˙m˙2+m˙Th 22S2,

где m˙Th пороговая скорость потока жидкости массы, ниже которой линеаризуется перепад давления. Данные показывают измененный перепад давления против локальной массовой скорости потока жидкости (изогнитесь I):

  • Выше m˙Th, перепад давления аппроксимирует, который выразил в исходном уравнении (изогнитесь II), и это меняется m˙2. Эта зависимость соразмерна с наблюдаемым в турбулентных течениях.

  • Ниже m˙Th, перепад давления аппроксимирует прямую линию с наклоном, частично зависящим от m˙Th (изогнитесь III), и это меняется m˙. Эта зависимость соразмерна с наблюдаемым в ламинарных течениях.

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

ξ2S2=Δp*ν*m˙*2,

где звездочка (*) обозначает значение в номинальных условиях работы. Лежание в основе всех этих вычислений является предположением что пороговая скорость потока жидкости массы m˙Th намного меньше, чем номинальная стоимость m˙*. Заменяя дробный ξ / (2S2) в выражении для урожаев перепада давления:

Δp=νΔp*ν*m˙*2(m˙m˙2+m˙Th 2).

или, эквивалентно:

Δp=Cνm˙m˙2+m˙Th 2,

где C является коэффициентом пропорциональности между перепадом давления через сопротивление потока и локальной массовой скоростью потока жидкости. Это задано как:

C=Δp*ν*m˙*2.

Если определенный объем — и поэтому массовая плотность — принята, чтобы быть инвариантной, то ее номинальная стоимость и фактические значения должны всегда быть равными. Дело обстоит так каждый раз, когда номинальная стоимость задана в диалоговом окне блока как 0 — специальное значение раньше сигнализировало к блоку, что определенный объем является константой. Отношением этих двух является затем 1 и продукт C, до которого уменьшает ν:

Cν=Δp*m˙*2.

Порты

Сохранение

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

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

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

Параметры

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

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

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

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

Область потока в портах сопротивления потока. Порты приняты, чтобы быть идентичными в размере.

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

Переменные

Массовая скорость потока жидкости в сопротивление через порт в начале симуляции.

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

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

Введенный в R2017b

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