🧪 演習6: メタノール-水の精留塔演習

1段の気液平衡(フラッシュ)から、12段の精留塔へ ― 純度を「作り込む」感覚をCOCOで確かめる

この演習について

演習1では「1回だけ加熱・蒸発させる」フラッシュ蒸留を扱い、フィード20 mol%が蒸気側で約51 mol%まで濃縮されることを確認しました。今回はその発展として、12段の平衡段を積み重ねた精留塔(COCOのColumn(ChemSep)ユニット)を使い、メタノール50 wt%のフィードから塔頂98.4 wt%・塔底はほぼ純水(0.14 wt%)という分離を作り込みます。

✅ この演習で確認すること

フィード(メタノール約36 mol%)を12段の塔に通すと、還流比3・フィード全体のうち63%を塔底側に回収する条件で、塔頂は約97.2 mol%(98.4 wt%)まで濃縮され、塔底はメタノールが約0.08 mol%(0.14 wt%)まで下がる。段数を積み重ね、塔底への回収割合を調整することで、フラッシュ1回では届かなかった高純度に到達できることを確認する。

結果の数値(先に答えを見ておく)

項目 フィード F 留出 D(塔頂) 缶出 W(塔底)
流量 [mol/s] 約 0.181 約 0.0669 約 0.1138
流量 [kg/h] 15.0 7.61 7.39
メタノール モル分率 0.360 0.972 約 0.0008
水 モル分率 0.640 0.028 約 0.9992
重量分率 50.0 wt% 98.4 wt% 0.14 wt%
温度 86.77°C(供給時) 66.90°C 100.37°C
圧力 [Pa] 101325 101325 101325

フィードの流量・組成・圧力はChemSepのFeedストリーム編集画面(Stream > Connections > Overall)に実際に入力した値(流量15 kg/h、圧力101325 Pa、モル分率Methanol 0.36/Water 0.64)です。温度はバブルポイント条件相当の86.77°Cでした。留出(D)・缶出(W)はSolve後のストリーム一覧表に表示された実測値(D:7.61444 kg/h、Mass frac Methanol 0.983864、66.9046°C/W:7.38556 kg/h、Mass frac Methanol 0.00137408、100.365°C)を、モル分率・mol/sに換算したものです。

フィード
36.0%
留出(塔頂)
97.2%
缶出(塔底)
≈0%

メタノールのモル分率をバーの長さで比較。フラッシュ1回(演習1)では蒸気側が約51%止まりだったのに対し、12段の精留塔では97.2%まで濃縮が進み、塔底側はほぼ0%(実質メタノールをほとんど含まない)まで下がることに注目。

なぜこの系・この条件を選ぶか

演習1で扱った水-エタノール系は約89 mol%エタノールで共沸点を持ち、普通の蒸留ではそれ以上濃縮できませんでした。一方、メタノール-水系は常圧で共沸点を作らない非共沸系です。そのため理論上は段数と還流比を増やせばいくらでも純度を上げられます。今回の12段・還流比3という条件は、「共沸の壁がない系では、段数を積み重ねることで狙った純度まで作り込める」という精留塔の基本を体感させるのに適した設定です。

フィード濃度50 wt%は、化学工業での粗メタノール精製や実験室スケールの回収工程に近い値で、塔頂98.4 wt%・塔底はほぼ純水という結果も、教材として「共沸の壁がなければどこまで分離が進むか」を実感させる現実的な数値です。

COCO(COFE)操作手順

Step 1: 新規フローシートを開く

  1. COFEを起動する。
  2. File > New で新規フローシートを作成する。

Step 2: Column(ChemSep)ユニットを配置する

  1. Insertメニュー(またはユニット操作挿入ボタン)からUnit Operationを選ぶ。
  2. Select Unit Operationウィンドウで、蒸留塔(Column/ChemSep Column)に該当する項目を選択し、フローシート上に配置する。

⚠️ カテゴリ名・項目名はバージョンにより異なる場合があります

COCOのバージョンによってUnit Operation一覧でのカテゴリ表記が異なることがあります。お手元の環境で実際に選んだカテゴリ名・項目名に読み替えてご利用ください。

Step 3: ChemSep側で熱力学モデルを設定する

  1. 配置したColumnユニットをダブルクリックし、ChemSep(ChemSep(TM))の設定画面を開く。
  2. 左側のツリーでComponentsを選び、MethanolWaterを検索して追加する。
  3. 左側のツリーでPropertiesを選び、上部のThermodynamicsタブを開く。
  4. Select Thermodynamic Modelsの各項目を次のように設定する:K-value = DECHEMA、Activity coefficient = Wilson、Vapor pressure = Extended Antoine、Enthalpy = Ideal(Equation of stateは選択不可のまま)。
  5. 右上のBIP estimationにチェックを入れる(BIP T est. = 350 K、T dependence = Noneのまま)。下部のWilsonタブ・Extended Antoineタブに、Methanol-Waterの二成分パラメータ(Aij = 858.975、Aji = 2017.36)が自動的に読み込まれる。

💡 演習1(NRTL)との違い

演習1のFlashユニットではCOCOのFlowsheet Configuration(Settings)でフローシート全体の活量係数モデルとしてNRTLを設定しましたが、今回のChemSepカラムユニットでは自分自身の設定画面(Properties > Thermodynamics)の中で個別に指定し、活量係数モデルにはWilsonを使います。フローシート全体のProperty packとは別枠になっている点、モデルの種類が異なる点の両方に注意してください。

⚠️ 画面構成・項目名はバージョンにより異なる場合があります

ChemSepの画面構成はバージョンによって多少異なることがあります。お手元の環境で実際の画面に読み替えてご利用ください。

Step 4: ストリームを接続する

  1. 演習1と同じ要領でInsert > Stream(Ctrl+I)を使い、フィードストリーム(F)をカラムのフィード段(6段目)に接続する。
  2. 塔頂(留出)ストリーム(D)をTotal Condenser側のDistillate出口に、塔底(缶出)ストリーム(W)をReboiler側のBottoms出口に接続する。

💡 ストリーム名について

塔底側のストリームはデフォルトでWという名前になります(Bottomsの意味ですが、表示上はWです)。以降の説明でも塔底ストリームはWと表記します。

Step 5: フィードストリームの条件を入力する

  1. フィードストリーム(F)を開くと、左側にStream > Connections > Overall というツリーが表示される。Overallを開き、pressure・temperature・mole fraction・flowの各欄に入力する。
  2. pressure = 101325 Pa、mole fraction [Methanol] = 0.36、mole fraction [Water] = 0.64flow = 15 kg/h(単位はmol/sではなくkg/h)を入力する。温度はバブルポイント(飽和液)条件とし、直接入力はしていない(この組成では約86.77°Cになる)。
  3. 入力後、MW欄に平均分子量(約0.0231 kg/mol)が自動計算されて表示される。この値が0.023付近になっていればモル分率の入力が正しい。

⚠️ 単位がFlashユニットと異なる

演習1のFlashユニットでは温度をK、流量をmol/sで入力しましたが、ChemSepのストリーム編集画面ではデフォルトで温度は°C、流量はkg/hで表示・入力する点が異なります。同じCOCO上でもユニットによって既定の単位系が違うので、入力のたびに単位欄を確認する習慣をつけてください。

Step 6: カラムの仕様を入力する

  1. 左ツリーのOperationで、理論段数 = 12フィード段 = 6段目、塔内圧力 = 101325 Pa(全段一定)を入力する。
  2. 左ツリーのSpecifications > Column specsを開く。Top specification(塔頂側の指定)はReflux ratioを選び、値に3を入力する。Bottom specification(塔底側の指定)はFraction of combined feeds recoveredを選び、値に0.630000を入力する。
  3. Top product name/Bottom product nameや、Condenser duty name(Qcondenser)/Reboiler duty name(Qreboiler)はデフォルトのままでよい。

✅ 入力するのはこの5つだけ

この演習で生徒が実際に触る数値は、フィード条件(流量・組成・圧力)・理論段数・フィード段・還流比・塔底への回収割合(Fraction of combined feeds recovered)の5つです。この5つを入れて計算を回すと、留出(D)は98.4 wt%、缶出(W)は0.14 wt%という高い純度の分離が得られます。

⚠️ 「Fraction of combined feeds recovered」とは

これは組成の目標値ではなく、フィード全体(モルベース)のうち何割を塔底側に回収するかという分割比の指定です。今回の0.63という値は、フィードの63%を塔底(W)へ、残り37%を塔頂(D)へ向かわせる設定です。塔頂から取り出す量を絞り込むことで、その中身(純度)を濃くしています。塔頂・塔底の実際の純度(wt%)は、この分割比・段数・還流比・フィード組成から結果として決まる値であり、直接指定しているわけではありません。

Step 7: フローシート上にストリーム一覧表を追加する

  1. 結果を見やすくするため、フローシート上にF・D・Wの3ストリームをまとめた一覧表を追加しておく。
  2. 表の行には Pressure・Temperature・Flow rate・Mass frac Methanol・Mass frac Water を並べる。実行前はD・Wの列がN/A(未計算)と表示される。

💡 なぜ実行前に表を作るのか

実行前にF・D・Wを横に並べた表を用意しておくと、Solveした瞬間にD・WのN/Aが実際の数値に置き換わり、フィードに対して留出・缶出がどう変化したかを一目で比較できます。フラッシュ蒸留(3-1)のフィード・蒸気・液の対比と同じ考え方です。

Step 8: 計算を実行する

  1. ツールバーの▶(Solve/Run)をクリックする。
  2. 精留塔は段数が多い分、フラッシュ単独より収束に時間がかかることがある。収束しない場合は還流比や初期値を見直す。

Step 9: 結果を確認する

  1. Step 7で作ったストリーム一覧表を見て、D・WのN/Aが実際の数値に置き換わっていることを確認する。
  2. 留出(D)でMass frac Methanolが約0.984(98.4 wt%相当)になっていることを確認する。
  3. 缶出(W)でMass frac Methanolが約0.0014(0.14 wt%相当、ほぼ純水)になっていることを確認する。
  4. フィード50 wt% → 塔頂98.4 wt%(純度が上がる)、塔底はほぼ純水(メタノールがほとんど残らない)という対比を確認する。
メタノール-水精留塔演習のシミュレーション結果画面

図:COCOでの実行結果

授業でのトーク例

Q1. 演習1のフラッシュ蒸留では蒸気側は約51%までしか濃縮できなかった。今回は98.4 wt%まで濃縮でき、塔底はほぼ純水になったのはなぜ?

→ 平衡段を1段でなく12段積み重ね、さらに還流で液を塔内に戻しているため。1段ごとに少しずつ濃縮が進み、それが積み重なる。

Q2. 水-エタノール系(演習1)ではどれだけ段数を増やしても89 mol%より濃くならなかった。メタノール-水系ではどうなる?

→ メタノール-水系は常圧で共沸点を作らないため、理論上は段数と還流比を増やせばさらに高純度(例えば99%以上)にも近づけられる。

Q3. 還流比を3より大きくすると、塔頂の純度と塔底の純度はどう変わる?

→ 一般に還流比を上げるほど分離性能は上がる(同じ段数でより高純度、または同じ純度をより少ない段数で達成できる)が、還流に使う熱量(凝縮器・再沸器の負荷)は増える。

Q4. フィード段を6段目でなく、もっと上(塔頂寄り)や下(塔底寄り)にずらすとどうなる?

→ フィード組成に対して最適な段からずれるため、同じ段数・還流比でも到達できる純度が下がることがある。実際にCOCOでフィード段を変えて比較させると効果的。

Q5. 塔底への回収割合(Fraction of combined feeds recovered = 0.63)を、もっと大きく(例えば0.8)するとどうなる?

→ 塔頂から取り出す留出の量がさらに絞られ、その中身は塔の一番上の「最も濃い部分」だけになるため純度はさらに上がる方向に働く。ただしその代わりに、缶出側にわずかに残るメタノールの量は増える傾向になる。純度を上げようとすると、取り出せる量や回収率とのトレードオフが生じることを実感させるポイント。

補足

フィードの流量・組成・圧力はChemSepの画面に実際に入力した値(15 kg/h、101325 Pa、モル分率Methanol 0.36/Water 0.64)そのものです。留出・缶出の数値もSolve後に画面に表示された実測値です。生徒には「50 wt% → 塔頂98.4 wt% / 塔底ほぼ0 wt%」という定性的な結論を中心に説明するとよいでしょう。