未知数として成分流量を使用してください。 化学工学の単位操作やパイロットプラントの物質収支では、各未知数を個別の成分質量流量、モル流量、または成分量として定式化します。未知のストリーム全流量に未知の質量分率やモル分率を掛け合わせるような表現は避けてください。そのような積は非線形方程式を生み出し、標準的な線形行列法による直接的な解法を妨げるためです。

最も信頼性の高い線形定式化は、まず成分流量を解き、その結果から全流量と組成を計算することです。これにより Ax = b 形式の方程式系が生成され、方程式が独立しており自由度がゼロであれば、標準的な数値線形代数を用いて効率的に解くことができます。

未知数の選択が線形性を決定する理由

非線形なストリーム表現

成分流量はしばしば次のように書かれます：

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

ここで、(\dot{m}) は全質量流量、(w_i) は成分質量分率です。

(\dot{m}) と (w_i) の両方が未知である場合、その積には2つの未知変数が含まれます。根底にある保存則は線形であっても、結果として得られる収支式は非線形になります。

線形な成分流量表現

成分流量そのものを未知数として定義します：

[ \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}}} ]

これにより、線形収支計算と組成計算を分離できます。

収支モデルの構築方法

計算基準（ベース）の確立

未知数を定義する前に基準を選択してください。連続定常運転の場合、通常は単位時間あたりの流量となります。バッチ運転の場合は、1バッチまたは指定された仕込み量となります。

基準間の変換が明示的に含まれていない限り、モデル全体を通して質量流量またはモル流量を一貫して使用してください。

プロセスを描き、ストリームに番号を振る

プロセスフロー図には、すべての流入、流出、リサイクル、分岐、および中間ストリームを特定する必要があります。未知のストリームに関連する各成分流量が候補変数となります。

例えば、2成分の出口ストリームは次のように表せます：

[ 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は、大学、研究機関、企業が、実用的なパイロットプラントシステム、成分ベースのモデリングサポート、および実験教育、プロセス開発、スケールアップ研究に適した機器を通じて、より明確なトレーニングおよびテストワークフローを構築できるよう支援します。