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

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

Почему выбор неизвестных определяет линейность

Нелинейное представление потока

Расход компонента часто записывается как:

[ \dot{m}_i = \dot{m}w_i ]

где (\dot{m}) — общий массовый расход, а (w_i) — массовая доля компонента.

Если и (\dot{m}), и (w_i) неизвестны, их произведение содержит две неизвестные переменные. Полученный баланс становится нелинейным, несмотря на то, что лежащий в его основе закон сохранения является линейным.

Линейное представление через компонентные расходы

Определите сам расход компонента как неизвестную величину:

[ \dot{m}_i = x_i ]

Тогда баланс по компоненту можно записать непосредственно как:

[ \sum_j a_{ij}x_j = b_i ]

где коэффициенты (a_{ij}) представляют собой известные соотношения потоков, коэффициенты распределения или стехиометрические коэффициенты. Вектор неизвестных содержит расходы компонентов, а не произведения общих расходов на составы.

Последующее определение характеристик потока

После того как компонентные расходы рассчитаны, определите общий расход потока:

[ \dot{m}_{\text{total}} = \sum_i \dot{m}_i ]

Затем рассчитайте массовую долю компонента, используя:

[ w_i = \frac{\dot{m}i}{\dot{m}_{\text{total}}} ]

Тот же подход применяется на молярной основе:

[ \dot{n}_{\text{total}} = \sum_i \dot{n}_i \qquad y_i = \frac{\dot{n}i}{\dot{n}_{\text{total}}} ]

Это отделяет расчет линейного баланса от расчета состава.

Как построить модель баланса

Установите базис расчета

Выберите базис перед определением неизвестных. Для непрерывных процессов в установившемся режиме это обычно расход в единицу времени; для периодических процессов это может быть одна партия или заданная загрузка.

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

Начертите схему процесса и пронумеруйте потоки

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

Например, двухкомпонентный выходной поток может быть представлен как:

[ x_1 = \dot{m}{A,\text{out}} \qquad x_2 = \dot{m}{B,\text{out}} ]

а не как пара, состоящая из общего выходного расхода и неизвестной доли состава.

Запишите независимые балансы компонентов

Для нереакторного узла в установившемся состоянии:

[ \text{поступление компонента} = \text{выход компонента} ]

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

Сформируйте матричный вид

После приведения каждого уравнения к виду, где неизвестные члены находятся слева, а известные — справа, запишите:

[ A\mathbf{x}=\mathbf{b} ]

Например:

[ \begin{bmatrix} a_{11} & a_{12}\ a_{21} & a_{22} \end{bmatrix} \begin{bmatrix} x_1\ x_2 \end{bmatrix}

\begin{bmatrix} b_1\ b_2 \end{bmatrix} ]

Здесь (A) содержит известные коэффициенты, (\mathbf{x}) содержит неизвестные компонентные расходы, а (\mathbf{b}) содержит известные расходы питания, спецификации или источники.

Проверка разрешимости системы

Выполните анализ степеней свободы

Рассчитайте:

[ \text{Степени свободы}

\text{количество неизвестных}

\text{количество независимых уравнений} ]

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

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

Проверьте независимость уравнений

Простого подсчета уравнений недостаточно, так как некоторые балансы могут быть линейно зависимыми. Проверьте ранг матрицы (A) перед решением.

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

import numpy as np

x = np.linalg.solve(A, b)

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

Проверьте сохранение после решения

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

[ \mathbf{r}=A\mathbf{x}-\mathbf{b} ]

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

Понимание компромиссов

Компонентные неизвестные увеличивают количество переменных

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

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

Измерения состава все еще могут быть необходимы

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

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

Некоторые модели процессов принципиально нелинейны

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

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

Отрицательные рассчитанные расходы требуют расследования

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

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

Правильный выбор для вашей цели

Используйте формулировку, соответствующую цели модели:

Если ваша основная задача — надежный линейный решатель: Определите индивидуальные массовые или молярные расходы компонентов как вектор неизвестных и решите полученную систему (A\mathbf{x}=\mathbf{b}).
Если ваша основная задача — определение состава потоков: Сначала найдите компонентные расходы, затем получите общие расходы и массовые или мольные доли путем суммирования и нормализации.
Если ваша основная задача — согласование данных пилотной установки: Включите измеренные данные о потоках и составах как независимые уравнения, затем проверьте ранг матрицы и невязки баланса.
Если ваша основная задача — детальное моделирование процесса: Сохраните переменные компонентных расходов для балансов сохранения, но используйте нелинейный решатель, когда этого требуют равновесие, кинетика или соотношения свойств.

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

Сводная таблица:

Элемент моделирования	Рекомендуемый подход	Преимущество
Основные неизвестные	Индивидуальные массовые или молярные расходы компонентов	Сохраняет линейность уравнений
Матричная модель	Выражение балансов как `A x = b`	Позволяет использовать стандартные линейные решатели
Общий расход	Суммирование компонентных расходов после решения	Отделяет расчет от состава
Состав	Деление расхода каждого компонента на общий расход	Обеспечивает массовые или мольные доли
Разрешимость	Проверка степеней свободы и ранга матрицы	Подтверждает возможность получения единственного решения
Валидация	Оценка невязок баланса и исследование отрицательных расходов	Выявляет проблемы моделирования или данных
Нелинейное поведение	Использование нелинейных решателей для равновесия, кинетики и свойств	Соответствует полной модели процесса

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