Моделирование газовых систем

Применения по назначению

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

  • Пневматическое приведение в действие механических систем

  • Транспорт природного газа через сети канала

  • Газовые турбины для производства электроэнергии

  • Воздушное охлаждение тепловых компонентов

Вы задаете газовые свойства в связанном цикле при помощи блока Gas Properties (G). Этот блок позволяет вам выбрать между тремя уровнями идеализации: совершенный газ, полусовершенный газ или действительный газ (см. Газовые Модели Свойства).

Сетевые переменные

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

Газовые модели свойства

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

Вы выбираете газовую модель свойства при помощи блока Gas Properties (G), который задает газовые свойства в подключенной схеме.

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

  • Тепловое уравнение состояния указывает на отношение плотности с температурой и давлением.

  • Тепловое уравнение состояния указывает на отношение удельной теплоемкости с температурой и давлением.

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

Газовая модель свойстваТепловое уравнение состоянияТепловое уравнение состоянияТранспортные свойства
СовершенныйИдеальный газовый законПостоянныйПостоянный
ПолусовершенныйИдеальный газовый закон1D поиск по таблице температурой1D поиск по таблице температурой
Действительный2D поиск по таблице температурой и давление2D поиск по таблице температурой и давление2D поиск по таблице температурой и давление

Идеальный газовый закон реализован в библиотеке Simscape™ Foundation Gas как

p = ZρRT

где:

  • p является давлением.

  • Z является фактором сжимаемости.

  • R является определенной газовой константой.

  • T является температурой.

Фактором сжимаемости, Z, обычно является функция давления и температуры. Это составляет отклонение от идеального газового поведения. Газ идеален когда Z = 1. В совершенных и полусовершенных газовых моделях свойства Z должен быть постоянным, но это не должно быть равно 1. Например, если вы моделируете неидеальный газ (Z ≠ 1), но температура и давление системы значительно не отличаются, можно использовать совершенную газовую модель и задать соответствующее значение Z. В следующей таблице перечислены фактор сжимаемости Z для различных газов в 293.15 K и 0,101325 МПа:

ГазФактор сжимаемости
Сухой воздух0.99962
Углекислый газ0.99467
Кислород0.99930
Водород1.00060
Гелий1.00049
Метан0.99814
Природный газ0.99797
Аммиак0.98871
R-134a0.97814

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

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

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

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

Блоки с газовым объемом

Компоненты в газовой области моделируются с помощью объемов управления. Объем управления охватывает газ в компоненте и разделяет его от окружающей среды и других компонентов. Потоки газа и тепловые потоки через поверхность управления представлены портами. Газовый объем в компоненте представлен с помощью внутреннего узла, который обеспечивает давление газа и температуру в компоненте. Этот внутренний узел не видим, но можно получить доступ к его параметрам и переменным с помощью регистрации данных Simscape. Для получения дополнительной информации займитесь Журналированием Данных моделирования.

Следующие блоки в библиотеке Gas моделируются как компоненты с газовым объемом. В случае Управляемого Водохранилища (G) и Водохранилища (G), объем принят, чтобы быть бесконечно большим.

БлокГазовый объем
Постоянная камера объема с 2 портами (G)Конечный
Постоянная камера объема с 3 портами (G)Конечный
Постоянная камера объема (G)Конечный
Передайте (G) по каналуКонечный
Вращательный механический конвертер (G)Конечный
Переводный механический конвертер (G)Конечный
Водохранилище (G)Бог
Управляемое водохранилище (G)Бог

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

Ссылочный узел и основывающиеся правила

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

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

Библиотека Foundation Gas содержит блок Absolute Reference (G), но, в отличие от других областей, вы не используете его для основания газовых схем. Цель блока Absolute Reference (G) состоит в том, чтобы предоставить ссылку Pressure & Temperature Sensor (G). Если вы будете использовать блок Absolute Reference (G) в другом месте в газовой сети, он инициирует утверждение симуляции, потому что давление газа и температура не могут быть в абсолютном нуле.

Начальные условия для блоков с конечным газовым объемом

В этом разделе рассматриваются определенные требования инициализации для блоков, смоделированных с конечным газовым объемом. Эти блоки перечислены в Блоках с Газовым Объемом.

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

Состояние газового объема представлено дифференциальными переменными во внутреннем узле блока. Как дифференциальные переменные, они требуют, чтобы начальные условия были заданы до запуска симуляции. Диалоговое окно каждого блока, смоделированного с конечным газовым объемом, имеет вкладку Variables, которая перечисляет три переменные:

  • Pressure of gas volume

  • Temperature of gas volume

  • Density of gas volume

По умолчанию Pressure of gas volume и Temperature of gas volume имеют высокий приоритет с целевыми значениями, равными стандартному условию (0.101325 MPa и 293.15 K). Можно настроить целевые значения, чтобы представлять соответствующее начальное состояние газового объема для блока. Density of gas volume имеет приоритет по умолчанию None, потому что только начальные условия двух из этих трех переменных необходимы, чтобы полностью определить начальное состояние газового объема. При желании альтернативный способ задать начальные условия состоит в том, чтобы изменить Density of gas volume на высокий приоритет с соответствующим целевым значением, и затем изменить или Pressure of gas volume или Temperature of gas volume к приоритету ни одного.

Важно, чтобы только двум из этих трех переменных установили их приоритеты на High для каждого блока с конечным газовым объемом. Размещение высокоприоритетных ограничений на все три переменные приводит к сверхспецификации с решателем, не могущим найти решение для инициализации, которое удовлетворяет желаемые начальные значения. С другой стороны размещение высокоприоритетного ограничения только на одну переменную делает систему под-заданным, и решатель может разрешить переменные с произвольными и неожиданными начальными значениями. Для получения дополнительной информации о переменной инициализации и контакте со сверхспецификацией, смотрите, Инициализируют Переменные для Системы Массового Демпфера Spring.

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

Дросселируемый поток

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

Следующая модель иллюстрирует дросселируемый поток. В этой модели блок Ramp имеет наклон 0,005 и время начала 10. Блоку Simulink-PS Converter установили Input signal unit на Mpa. Все другие блоки имеют значения параметров по умолчанию. Время симуляции составляет 50 с. Когда вы моделируете модель, давление в порте блока Local Restriction (G) увеличивается линейно с атмосферного давления, запускающегося в 10 с. Давление в порте B фиксируется в атмосферном давлении.

Следующий рисунок показывает регистрируемые данные моделирования для блока Local Restriction (G). Число Маха в ограничении (Mach_R) достигает 1 приблизительно в 20 с, указывая, что поток дросселируется. Массовая скорость потока жидкости (mdot_A) перед потоком дросселируется, следует за типичным квадратичным поведением относительно увеличивающегося перепада давлений. Однако массовая скорость потока жидкости после потока дросселируется, становится линейным, потому что дросселируемая массовая скорость потока жидкости зависит только от восходящего давления и температуры, и восходящее давление увеличивается линейно.

То, что дросселируемая массовая скорость потока жидкости зависит только от восходящих условий, может вызвать несовместимость с Массовым Источником Скорости потока жидкости (G) или Управляемым Массовым Источником Скорости потока жидкости (G) соединенный в нисходящем направлении дросселируемого блока. Считайте модель показанной на следующем рисунке, который содержит блок Controlled Mass Flow Rate Source (G) вместо Управляемого Источника Давления (G).

Если бы источник управлял увеличивающейся массовой скоростью потока жидкости слева направо через Локальное Ограничение (G), симуляция успешно выполнилась бы, даже если бы поток стал дросселируемым, потому что Управляемый Массовый Источник Скорости потока жидкости (G) будет восходящим из дросселируемого блока. Однако в этой модели блок Gain инвертирует направление потока, так, чтобы Управляемый Массовый Источник Скорости потока жидкости (G) был нисходящим из дросселируемого блока. Давление в восходящем направлении Локального Ограничения (G) фиксируется в атмосферном давлении. Поэтому дросселируемая массовая скорость потока жидкости в этой ситуации является постоянной. Когда массовая скорость потока жидкости, которой управляют, увеличивается, в конечном счете это станет больше, чем это постоянное значение дросселируемой массовой скорости потока жидкости. На данном этапе массовая скорость потока жидкости, которой управляют, и дросселируемая массовая скорость потока жидкости не могут быть согласованы и сбои симуляции. Просмотр регистрируемых данных моделирования в Проводнике Результатов Simscape показывает, что симуляция перестала работать только в точке, когда Число Маха достигает 1, и поток становится дросселируемым.

В целом, если модель, вероятно, будет дросселировать, используйте источники давления, а не массовые источники скорости потока жидкости. Если модель содержит массовые исходные блоки скорости потока жидкости и сбои симуляции, используйте Проводник Результатов Simscape, чтобы осмотреть переменные Числа Маха во всем Локальном Ограничении (G), Переменном Локальном Ограничении (G) и Канале (G) блоки, соединенные вдоль того же ответвления как массовый источник скорости потока жидкости. Если отказ симуляции происходит, когда Число Маха достигает 1, вероятно, что существует нисходящий массовый источник скорости потока жидкости, пытающийся управлять массовой скоростью потока жидкости, больше, чем возможная дросселируемая массовая скорость потока жидкости.

Переменная Числа Маха для блоков ограничения называется Mach_R. Блок Pipe (G) имеет две переменные Числа Маха, Mach_A и Mach_B, представляя Число Маха в порте A и порте B, соответственно.

Реверсирование потока

Поток газа через схему несет энергию от одного газового объема до другого газового объема. Поэтому энергетическая скорость потока жидкости между двумя связанными блоками зависит от направления потока. Если потоки газа от блока A до блока B, то энергетическая скорость потока жидкости между двумя блоками основана на определенной общей энтальпии блока A. С другой стороны, если потоки газа от блока B до блока A, то энергетическая скорость потока жидкости между двумя блоками основана на определенной общей энтальпии блока B. Чтобы сглаживать переход для робастности симуляции, энергетическая скорость потока жидкости также включает вклад на основе различия в определенных общих энтальпиях двух блоков в скоростях потока жидкости малой массы. Областью сглаживания управляет Gas Properties (G) параметры блоков Mach number threshold for flow reversal.

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

Робастность симуляции может быть сложной для моделей, которые показывают быстрые реверсирования потока и большой перепад температур между блоками. Быстрые реверсирования потока могут быть результатом наличия низких сопротивлений потока (например, коротких каналов) между большими газовыми объемами. Большим перепадом температур может быть результат энергии, добавленной источниками, чтобы поддержать большой перепад давлений в модели с небольшой теплоотдачей. В этих моделях может быть необходимо увеличить значение параметров Mach number threshold for flow reversal, чтобы избежать отказа симуляции.

Похожие темы