2-Port Constant Volume Chamber (2P)

Ёмкость с двумя портами и фиксированным объемом двухфазной жидкости

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

  • 2-Port Constant Volume Chamber (2P) block

Описание

Блок 2-Port Constant Volume Chamber (2P) моделирует накопление массы и энергии в ёмкости, содержащей фиксированный объем двухфазной жидкости. Ёмкость имеет два входных отверстия, обозначенных A и B, через которые может течь жидкость. Объем жидкости может обмениваться теплом с тепловой сетью, например с той, которая представляет окружающую ёмкость, через тепловой порт, маркированный H.

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

Сопротивление потоку между каждым входным отверстием и внутренней частью ёмкости принято незначительным. Поэтому давление во внутреннем пространстве равно давлению во входных отверстиях. Аналогично, тепловое сопротивление между тепловым портом и внутренней частью ёмкости принято незначительным. Температура в интерьере равна температуре в тепловом порту.

Баланс массы

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

[(ρp)udpdt+(ρu)pdudt]V=m˙A+m˙B+ϵM,

где левая сторона является скоростью накопления массы и:

  • ρ - плотность.

  • p - давление.

  • u - специфическая внутренняя энергия.

  • V - объем.

  • m˙ - массовый расход жидкости.

  • ϵ M является членом коррекции, введенным для расчета числовой ошибки, вызванной сглаживанием частных производных.

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

Частные производные в уравнении баланса массы вычисляются путем применения метода конечной разности к сведенным в табличным данным в блоке Two-Phase Fluid Properties (2P) и интерполяции результатов. Частные производные затем сглаживаются на контурах фазового перехода с помощью кубических полиномиальных функций. Эти функции применяются между:

  • Переохлажденная жидкость и двухфазная смесь фазы домены, когда качество пара находится в 0-0,1 областях значений.

  • Двухфазная смесь и перегретые пары фазы домены, когда качество пара находится в 0-0,9 области значений.

Сглаживание вводит небольшую числовую ошибку, которую блок корректирует путем добавления к балансу массы термина коррекции ϵ M, заданного как :

ϵM=MV/ντ.

где:

  • M - масса жидкости в ёмкости.

  • ν - конкретный объем.

  • τ - характеристическая длительность события изменения фазы.

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

dMdt=m˙A+m˙B.

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

Энергия может войти и выйти из ёмкости двумя способами: с потоком жидкости через порты А и B и с потоком тепла через порт H. Никакая работа не выполняется ни с жидкостью внутри ёмкости, ни с ней. Скорость накопления энергии во внутреннем объеме жидкости должна равняться сумме энергии, скоростей потока жидкости через порты А, B и H:

E˙=ϕA+ϕB+QH,

где:

  • ϕ - скорость потока жидкости энергии.

  • Q - тепловая скорость потока жидкости.

  • E - полная энергия.

Пренебрегая кинетической энергией жидкости, общая энергия в ёмкости:

E=Mu.

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

Перепад давления из-за вязкого трения между отдельными портами и внутренней частью ёмкости принимается незначительным. Гравитация игнорируется, как и другие силы тела. Давление во внутреннем объеме жидкости должно тогда равняться давлению в порте А и в порту B:

p=pA=pB.

Предположения

  • Ёмкость имеет фиксированный объем жидкости.

  • Сопротивление потоку между входным и внутренним отверстиями ёмкости незначительно.

  • Тепловое сопротивление между тепловым портом и внутренней частью ёмкости незначительно.

  • Кинетическая энергия жидкости в ёмкости незначительна.

Порты

Сохранение

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

Открытие, посредством которого жидкость может войти и выйти из ёмкости.

Открытие, посредством которого жидкость может войти и выйти из ёмкости.

Интерфейс, посредством которого жидкость в ёмкости обменивается теплом с тепловой сетью.

Параметры

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

Вкладка « параметры»

Объем жидкости в ёмкости. Этот объем является постоянным во время симуляции.

Площадь входного отверстия, перпендикулярная направлению потока.

Площадь входного отверстия, перпендикулярная направлению потока.

Вкладка Эффекты и начальные условия

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

Давление в ёмкости в начале симуляции, заданное относительно нуля.

Температура в ёмкости в начале симуляции, заданная при абсолютном нуле.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Temperature.

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

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Vapor quality.

Объемная доля пара в ёмкости в начале симуляции.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Vapor void fraction.

Специфическая энтальпия жидкости в ёмкости в начале симуляции.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Specific enthalpy.

Удельная внутренняя энергия жидкости в ёмкости в начале симуляции.

Зависимости

Этот параметр активен, когда опция Initial fluid energy specification установлена на Specific internal energy.

Характерное время до равновесия события изменения фазы, происходящего в ёмкости. Увеличьте этот параметр, чтобы замедлить скорость изменения фазы или уменьшите его, чтобы ускорить скорость.

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

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

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