例題集

2005/3/11

もくじ
１．アンモニア製造プロセスの物質収支計算 ......................................... 3
ステディステート（定常状態）シミュレーション
物質収支計算ストリーム接続（総モル流量、モル組成を伝達）

２．バッチ反応器のシミュレーション ......................................................... 9
ダイナミックシミュレーション
信号線接続
シーケンス記述（自動操作）

３．シクロヘキサン製造プロセスの物質・熱収支計算 ......................... 15
ステディステート（定常状態）シミュレーション
物質収支・熱収支計算
物性値ライブラリと物性関数の利用ストリーム接続（総モル流量、モル組成、温度、圧力を伝達）2

１．アンモニア製造プロセスの物質収支計算
プロセスの概念設計を想定した例として、アンモニア製造のプロセスの定常計算を考える。
登場する装置は反応器と分離器のみというもっとも簡単なものである。なお、総モル流量、モル組成を伝達するストリームを使った例となる。
■問題
アンモニアは以下のように窒素と水素から直接合成される。

N2 + 3H2

2NH3

原料は乾燥した水素と窒素ガスで、不純物としてメタンを含む。反応器では 20MPa、500℃の条件で触媒の存在下で反応が行われる。反応器出のガスは冷却、凝縮により分離し、未反応ガスはリサイクルされて、液体アンモニアは製品として取り出される。また、未反応
ガス中の不純物（メタン）の蓄積をパージにより防ぐ。
原料の組成（モル分率）がおのおの N2 : 0.24、H2 : 0.74、CH4 : 0.02 で、その流量が 100[kmol/h] の場合の物質収支を計算する。なお、リサイクル未反応ガスの 2%がパージされるものとする。
3

パージ
6

冷媒
反応器

5

分離器
7

N2,H2,
CH4

4

2
1

リサイクル

NH3

図 1.1 アンモニア製造プロセス■考え方
プロセスの物質収支を求めるが、ここではその処理を大まかなブロックとして考える。反
応器、分離器のほかに、合流点、分流点がある。また、バッテリーリミットとして、原料（フィード）と製品、パージがある。それぞれのブロックで入出の物質の収支式を記述し3

て、全体としての物質収支を計算する。
ブロック間には値を伝達するものとしてストリームを定義する。ここでは物質収支のみを考えるので、総モル流量（F）とモル組成（x）を考えればよい。成分は N2、H2、NH3、CH4 の４成分を登録する。
問題から図 1-2 のようなブロック・フロー・ダイアグラムを作成する。s107 ではリサイク
ル指定が必要となる点に注意されたい。
アンモニア製造プロセス－１

1

1

s103

1

Feed

1

s101

1

1

1

s104

2

s105

1

Prod

1

Purge

Separatr

Reactor

s102

Mixer
2

1

1

Splitter
2

s106
R

s107

図 1.2 ブロックフロー図
いま、総モル流量とモル組成のみを考えているので、バッテリーリミット（Feed、Prod、
Purge）や合流点（Mixer）、分流点（Splitter）の記述はほぼ決まってくる。合流点や分流点では、総モル流量と成分モル流量の収支式を記述することになる。これについては次ページの囲みを参照されたい。
ユニット
Feed

Prod / Purge

変数

入出力

F : フィード量 [kmol/h]

INPUT

x : フィード組成 [-]INPUT

F : 流量 [kmol/h]

OUTPUTx : 組成 [-]OUTPUT

Mixer
Splitterr : パージ率 [-]

INPUT

4

ソーステキスト本体
p1.F = F
p1.x = x
F = f1.F
x = f1.x
f1.F + f2.F = p1.F
f1.F*f1.x + f2.F*f2.x = p1.F*p1.x
f1.F = p1.F + p2.F
f1.F*f1.x = p1.F*p1.x + p2.F*p2.x
p1.F = f1.F * r
p1.F*p1.x = f1.F*f1.x * r

物質収支の基本式
ブロックあるいはユニット単位で物質収支をとらえると、一般に単位時間あたり次式のようになる。
（蓄積量）＝（流入量）―（流出量）＋（生成・消滅量）
定常状態を扱う場合には時間あたりの変化を考えないので、左辺は 0 となり（蓄積量）の項は消去できる。
（流入量）―（流出量）＋（生成・消滅量）= 0
つぎに流入の数、流出の数、反応のある／なしで考えると、反応器、合流点、分流点／分
離器の３タイプが基本となる。この際に総モル流量、成分流量を考慮するが、総モル流量
は Fi、モデル組成は xi,j で表すと、成分流量は Fi xi,j となる。
■反応器（流入１、流出１）
生成・消滅を反応で代表しているが、他の場合（たとえば漏れ）でも同様に考えることができる。反応生成量の総モル量を R、成分ごとの反応生成量を rj で表す。もし、生成・消滅がない場合には R = 0 と考えればよい。

F1 + R = F2

F1

F1 x1, j + r j = F2 x2, j

R

F2■合流点（流入２、流出１）
合流点は単に２つの流入が加算されて１つの流出になる。

F2

F1 + F2 = F3
F1 x1, j + F2 x2, j = F3 x3, j

F1

F3■分流点／分離器（流入１、流出２）
分流点および分離器では分配比αを導入して分配割合を表現する。分離器の場合には成分ごとの分離比を指定するので、αは１次元配列となる。

F1 = F2 + F3

F1

F1 x1, j = F2 x2, j + F3 x3, j

α

F2
F3

F2 = F1 ⋅ α
F2 x2, j = F1 x1, j ⋅ α5残る主要な装置である反応器と分離器をどのようにモデル化するかで、いろいろなアプローチが考えられよう。
■解法１
はじめにもっとも簡単なモデルで解法してみる。
・反応器での反応量はフィードされた窒素の転化率で与える。
・分離器での製品への分離は、入ってくる成分流量に対する分離比で与える。
＜反応器（Reactor）＞
転化率を conv として、種別：INPUT、型：REAL として定義する。ソーステキスト本体は反応量[kmol/h]を r とし、化学量論比を考慮して以下のようになる。
VAR r(4) "反応量 [kmol/h]"
f1.F + SUM(r) = p1.F
f1.F*f1.x + r = p1.F*p1.x-r(1) = f1.F*f1.x(1) * conv-r(2) = -r(1) * 3
r(3) = -r(1) * 2
r(4) = 0＜分離器（Separatr）＞
分離比を ratio として、種別：INPUT、型：COMP として定義する。すなわち、ratio は成分ごとに考えるので配列変数となる。
f1.F = p1.F + p2.F
f1.F*f1.x = p1.F*p1.x + p2.F*p2.x
p1.F*p1.x = f1.F*f1.x * ratio
eq: SUM(p1.x) = 1
RESET p1.F#100[0,10000] BY eq■解法１の計算結果
転化率 conv=0.16、分離比 ratio=（0.0006、0.00071、0.512、0.0010）として定常計算を行う。リサイクルをデフォルトのままにしておくが、収束は緩やかで 450 回ほど繰り返しが必要となる。

6

アンモニア製造プロセス－１

ratio(1)=0.0006
ratio(2)=0.00071
ratio(3)=0.512
ratio(4)=0.001

conv=0.16
1

1

Reactor

Feed

1

s101

1

1

Mixer
2

F=100
x(1)=0.24
x(2)=0.74
x(3)=0x(4)=0.02s103

1

F=727.92
x(1)=0.156211
x(2)=0.598859
x(3)=0.114055
x(4)=0.130875

Separatr

s104

2

s105

1

Prod

F=42.9806
x(1)=0.00158736
F=684.94
x(2)=0.00720102
x(1)=0.165914 x(3)=0.988995
x(2)=0.635986 x(4)=0.00221651
x(3)=0.0591513
x(4)=0.138949

s102

F=771.238
x(1)=0.175521
x(2)=0.649473
x(3)=0.0514818
x(4)=0.123524

1

r=0.02
1

1

Splitter
2

s106
R

s107

F=671.241x(1)=0.165914
x(2)=0.635986
x(3)=0.0591513
x(4)=0.138949

1

Purge
F=13.6988
x(1)=0.165914
x(2)=0.635986
x(3)=0.0591513
x(4)=0.138949■解法２
解法１に比べてモデル精度を上げ、つぎのように反応モデル、分離モデルを考える。
・平衡反応を考慮し、化学平衡定数を与えて計算する。
・分離器では、アンモニアは分離比で与えるが、その他の気体がアンモニアにとけ込む溶
解度[kmol/kmol]を与えて計算する。
＜反応器（Reactor）＞
化学平衡定数 K を、圧力を P として以下のように定義する。
2

K = K ' P2 =x NH 3
3
xN 2 xH 2

化学平衡定数を K として、種別：INPUT、型：REAL として定義する。また、化学量論から N
原子および H 原子について、入り出の原子数が等しいことを考慮して以下のソーステキス
トを作成する。
f1.F*f1.x(1)*2 + f1.F*f1.x(3)= p1.F*p1.x(1)*2 + p1.F*p1.x(3)f1.F*f1.x(2)*2 + f1.F*f1.x(3)*3 = p1.F*p1.x(2)*2 + p1.F*p1.x(3)*3
f1.F*f1.x(4) = p1.F*p1.x(4)eq2: p1.x(3)^2 = K*p1.x(1)*p1.x(2)^3
eq1: SUM(p1.x) = 1
RESET p1.x(1)#0.1[0,1] BY eq1
RESET p1.F#100[0,10000] BY eq2＜分離器（Separatr）＞
分離比を ratio（種別：INPUT、型：REAL）とし、溶解度を sol（種別：INPUT、型：COMP）と定義する。溶解度は、それぞれ N2 : 0.00160、H2 : 0.00740、CH4 : 0.00225 というデー7タを用いる。
f1.F = p1.F + p2.F
f1.F*f1.x = p1.F*p1.x + p2.F*p2.x
p1.F*p1.x(3) = f1.F*f1.x(3) * ratio
p1.F*p1.x(1) = p1.F*p1.x(3) * sol(1)p1.F*p1.x(2) = p1.F*p1.x(3) * sol(2)p1.F*p1.x(4) = p1.F*p1.x(3) * sol(4)eq: SUM(p1.x) = 1
RESET p1.F#1 BY eq■解法２の計算結果
化学平衡定数 K=0.35、分離比 ratio=0.512 として定常計算を行う。この場合にも収束は緩やかで、360 回あまりの繰り返しが必要となる。

アンモニア製造プロセス－２
K=0.35
1

Feed
F=100
x(1)=0.24
x(2)=0.74
x(3)=0x(4)=0.02

1

s101

1

1

Mixer

Reactorratio=0.512
1

s103

1

F=747.872
x(1)=0.158923
x(2)=0.603194
x(3)=0.110488
x(4)=0.127395

s104

Prod

1

Separatr
F=42.7828
x(1)=0.0015822
F=705.089
x(2)=0.00731768
x(1)=0.16847
x(3)=0.988875
x(2)=0.639351 x(4)=0.00222497
x(3)=0.0571895
x(4)=0.13499

2

s102

F=790.986
x(1)=0.177514
x(2)=0.652076
x(3)=0.0499595
x(4)=0.120451

2

1

s105

r=0.02
1

Splitter

1
2

s106
R

s107

F=690.988
x(1)=0.16847
x(2)=0.639351
x(3)=0.0571895
x(4)=0.13499

Purge

1

F=14.1018
x(1)=0.16847
x(2)=0.639351
x(3)=0.0571895
x(4)=0.13499■実行の操作手順
１．エクスプローラを使って、Ammonia フォルダにある Ammonia1.vmp をダブルクリッ
クして VM-EQUATRAN を起動する。
２．[ファイル]-[実行]（）で実行モードに移る。

３．[設定]-[ループオーバー回数]で回数を 500 に変更する。
４．[実行]-[起動]（）で起動し、定常計算ボタン

５ ． [フ ァ イ ル ]-[編 集 ]（をクリックして収束計算を始める。

） で 編 集 モ ー ド に 戻 り 、 [フ ァ イ ル ]-[開 く ]（Ammonia2.vmp を読み込んで同様の手順で実行する。

8）で

２．バッチ反応器のシミュレーション
ダイナミックシミュレーションの例として、ジャケット付きバッチ反応器の運転を取り上げる。なお、この例題は EQUATRAN-G ユーザーズガイド［詳細編］の例題 15 に対応している。
ユニットとしては、反応器、ジャケット、PID 制御器を考え、データの受け渡しには信号線を使う。また、バッチプロセスには運転シーケンスが付き物であるが、自動操作の機能を使って実現する。
■問題
図 2.1 のようなバッチ反応器を用いて、原料 A および B から反応
A + 3B → C
によって製品 C を製造する。

図 2.1 バッチ反応器
つぎのシーケンスによって操作を行った場合のシミュレーションを考える。
①初期状態 ······ 溶媒 V=10m3 と原料 A 10kmol が装填されており反応器およびジャケットの温度は 20℃になっている。
②昇温 ············ 熱媒(120 ℃)を用いて反応器の温度 T が 90 ℃になるまで昇温する。
③反応············ 原料 B を図 2-2 に示すようなパターンで、 積算量 30kmol になるまで供給する。反応器の温度は 90℃に制御するが、最初の 0.2h の間は手動でジ
ャケット温度を 88 ℃に設定するものとする。
④保温············ 原料の供給および、冷媒の供給を停止したまま、0.5h 保温して反応を進行させる。
⑤冷却 ············ 冷媒(20 ℃)で反応器の温度が 40℃になるまで冷却する。

9

図 2.2 運転シーケンス
計算条件は次のとおりとする。
・反応速度 r[kmol/h m3]

r = k ⋅ C A ⋅ C B3⎛−E⎞
k = A ⋅ exp⎜⎟
⎝ RT ⎠
ただし、A=2.43×1015、E=83700、R=8.3145 であり、CA 、CB は A および B の濃度[kmol/m3]である。
・反応熱∆H = －50200 kJ/kmol

・溶媒の比熱cp = 1.86 kJ/kg℃

・溶媒の密度ρ＝900 kg/m3

・ジャケットの伝熱面積

Aj = 20 m3

・ジャケットの伝熱係数

Uj = 580 W/m2 K

・PID 制御器 TC1 のパラメータ（制御動作は逆動作）
比例帯 PB =50 %積分時間 TI =1000 s
微分時間 TD =360 s
・ジャケット温度

Tj ℃

なお、ジャケット温度は TC2 の設定値 Tjset に対して、時定数 360 s の一次おくれで表わせるものとする。

10■考え方
各ユニットには以下のような外部変数を指定して接続先と信号線で結ぶ。
変数名

内容

接続先

FB

B のフィード量

Bfeed.FB

Tj

ジャケット温度

Jacket.Tj

ジャケット（Jacket） Tjset

ジャケット温度設定値

PID.OUTPID 制御器（PID）

入力値

Reactor.T

反応器（Resctor）

IN

原料 B のフィード量は 10kmol/h2 のスピードでフィード量を増やし、20kmol/h で一定量にしてから、累計 30kmol/h で止める。このためにこのフィード量を計算するフィード量の調整というユニットを用意した。
バッチ反応器
フィード量の調整

T

F

ジャケット

T

反応器

I

PID制御器図 2.3 ブロックフロー図■各ユニットモデル＜フィード量の調整＞
10kmol/h2 のスピードでフィード量を増やし、20kmol/h で一定量にしてから、累計30kmol/h で止めるという操作を表現するためのソーステキストを作成した。FBon はフィー
ド用スイッチで、１になるとフィードが開始される。
FB = min(t1*10,20) * (FBon==1)
t1' = 1;
t1#0
CHANGE t1#0 ON (FBon==1)
FBtotal' = FB

＜反応器＞
問題に与えられた反応速度式を使い、濃度変化を微分方程式で表現した。ここで、FB は B

11のフィード量である。

CA
= −r
dt
CB
F
= −3 ⋅ r + B
dt
V
Cc
=r
dt
反応器の温度はつぎの熱収支式で求められる。

Vρc p

dT
= Qr − Q jdt

ここで Qr は反応熱、Qj はジャケットによる除熱量であり、つぎのように与えられる。

Qr = ∆HVr

Q j = 3.6 ⋅ U j A j (T − T j )/* 反応 */
CA' = -rCB' = -3*r + FB/V
CC' = rr = k*CA*CB^3

k = 2.43e15*exp(-83700/8.3145/(T+273.15))/* 熱収支 */V*rho*cp*T' = Qr - Qj
Qr = -dH*V*r
Qj = 3.6*Uj*Aj*(T-Tj)
OUTPUT _time, CA, CB, CC, T, Tj STEP 0.01

＜ジャケット＞/* １次遅れ要素 */Tj' = (Tjset-Tj)/(tau/3600)＜PID 制御器＞
MAN モード、AUTO モードを切り替えて動作でき、制御動作（正動作／逆動作）の指定が可能な汎用の PID 制御器ユニットに仕上げた。PID 制御器の出力は次式で計算される。

12

⎛
1
de ⎞
OUT = K p ⎜⎜ e + ∫ edt + TD ⎟⎟
TI
dt ⎠
⎝
100
Kp =
PB
e = PV − SV
PV = IN
TIh = TI/3600
TDh = TD/3600
e = PV-SV
U0 = MV*TIh*PB/100= -MV*TIh*PB/100 WHEN INCDEC == 1
U1 # 0
CHANGE U1 # U0 ON (MODE!=1)U1' = e= 0 WHEN MODE == 1
U2 = 100/PB*(e+U1/TIh+TDh*DERIV(e,0))U3 = U2= -U2 WHEN INCDEC == 1
OUT = U3
= MV WHEN MODE == 1

＜自動操作＞
動かしながら画面から値を変更することも可能だが、問題に与えられた運転シーケンスを実現するために、以下のような操作表を作成した。PID 制御器の制御モード MODE や操作値
MV、設定値 SV とフィード用 SW である FBon に値を書き込んでいる。
[time == 00:01:00]
PID.MV = 120[time >= +00:00:00]& [Reactor.T >= 90]
Bfeed.FBon = 1
PID.MV = 88[time >= +00:12:00]
PID.MODE = 2
PID.SV = 90
[time >= +00:00:00]& [Bfeed.FBtotal >= 30]
Bfeed.FBon = 0
PID.MODE = 1
PID.MV = 90
[time >= +00:30:00]
PID.MODE = 1
PID.MV = 20
[time >= +00:00:00]

13

& [Reactor.T <= 40]
FREEZE "反応終了"■計算結果
64 倍速で約３時間半に相当するシミュレーションを行った。

図 2.4 計算結果
■実行の操作手順
１．エクスプローラを使って、BatchReact フォルダにある BatchReact.vmp をダブルクリ
ックして VM-EQUATRAN を起動する。
２．[ファイル]-[実行]（３．[ファイル]-[読み込み]（）で実行モードに移る。
）で初期状態ファイル AutoOper.vmi を開く。

４．[モニタ]-[トレンド]-[表示]（）で２つのトレンドを表示するようにしておく。

５．タイムスケールを「実時間」から「64 倍速」に変更する。
６．運転ボタンで自動実行する。

14

３．シクロヘキサン製造プロセスの物質・熱収支計算「アンモニア製造プロセス」の例題に比べて、熱収支やフラッシュ計算を取り入れた、より精密な定常状態シミュレーションである。ストリームによるデータの受け渡しには、総モル流量 F、モル組成 x のほかに、温度 T、圧力 P も考慮している（ストリーム定義でこの４つの変数にチェックを入れておく）。また、物性計算では各成分に対する物性値を標準
物性ライブラリから読み込んで利用し、エンタルピ計算やフラッシュ計算に物性計算関数を利用した。
（注）標準物性ライブラリや物性計算関数を使うためには、物性計算オプションが必要となる。
■問題
シクロヘキサン（C6H12）は以下のようにベンゼン（C6H6）の水素添加反応で合成される。
ここでは、主要装置を考慮して、メインのフローを考えることにする。
C6H6 + 3H2

C6H12

液体ベンゼンと補給水素は、リサイクル水素、リサイクル液と混合して反応器にフィードされる。補給水素には不純物として窒素が含まれる。反応器では、2.1MPa、200℃の条件下で触媒を用いて水素添加反応を行う。反応器出の飽和蒸気は冷却したのち、フラッシュ
ドラムで気液に分離する。気相は一部をパージ（8.2%）し、未反応水素はコンプレッサで加圧してリサイクルする。液相はさらに減圧してフラッシュ分離して製品シクロヘキサンを取り出す。また液相の一部はリサイクルする。
12

8.2%

11

パージ

9
7

8
フラッシュ
ドラム

4

反応器
15

ベント

10
6

H2、N2
C6H6

5

3

13

14

フラッシュ
ドラム

62%

1

2

16

図 3.1 シクロヘキサン製造プロセス

15

C6H12

原料ベンゼンは 41.83kmol/h、補給水素は 128.84kmol/h（窒素が 0.296mol%含まれる）で供給される場合の物質収支、熱収支を計算せよ。なお、反応器ではフィードされたベンゼ
ンの 99.8%が反応する（転化率が 0.998）ものとする。
出典：W.D.Seider、J.D.Seader、D.R.Lewin 著『Process Design Principles』P.232（1999）■考え方
物質収支計算の考え方は「アンモニア製造プロセス」の例題と同様である。
ここでさらに考慮したのは温度と圧力である。各成分の物性値に物性ライブラリを、温度（エンタルピ）計算やフラッシュ計算には物性計算関数を利用した。物性計算関数の引数については『VM-EQUATRAN チュートリアル』Appendix B に一覧があるので参照されたい。
まず、成分として H2、N2、ベンゼン、シクロヘキサンをライブラリから選択し、計算法は水素を含む高圧系なので SRK 式を採用した。
問題から図 3.2 のようなプロセスフロー図（PFD）を作成した。今回は各ユニットについてそれぞれユニット図形を作成して登録した。

シクロヘキサン製造プロセス

R

s104

1
2

1
1

1

s112
s111

BL03

1

SP01
1

GB01

1

s108
s109

s107

CW01

1

1

1

RA01

1

BL04

1

BL05

FA01
s115

BL02

BL01

1

1

s103

s101

1

1

GA01

s102

1

s110

1
1
2
3
4

1

VT01

2

1

s117
1

1

1
2

1

s106

MX01

FA02

s114
s113

SP02
2

1

1

R

s116
s105

GA02

図 3.2 シクロヘキサン製造プロセスのプロセスフロー図＜反応器 RA01＞
反応器では転化率と反応条件（温度、圧力）が与えられるので以下のようになる。

16

VAR xfl(_NC), xfl0(_NC) ,xfv(_NC), xfv0(_NC),xpl(_NC), xpl0(_NC) ,xpv(_NC), xpv0(_NC) ,Dh(_NC)/* 物質収支 */f1.F + SUM(r) = p1.F
f1.F*f1.x + r = p1.F*p1.x-r(1) = -r(3) * 3
r(2) = 0-r(3) = f1.F*f1.x(3) * conv
r(4) = -r(3)/* conv:INPUT *//* 出口温度・圧力 */T = p1.T/* T:INPUT */P = p1.P/* P:INPUT *//* 除熱量 */PHYHCLV( 1, f1.T, f1.P, 1, f1.x, hfl, rfl, xfl, hfv, rfv, xfv, hf )/* 入口エンタルピ */Hf = hf*1000
PHYHCLV( 1, p1.T, p1.P, 1, p1.x, hpl, rpl, xpl, hpv, rpv, xpv, hp )/* 出口エンタルピ */Hp = hp*1000
PHYGTDHV( 1, Dh )/* 標準生成熱 */Qr = (Dh(4)-Dh(3)-Dh(1)*3)*1000
Q = Hp*p1.F - Hf*f1.F - Qr*r(3)反応での生成・消滅量の計算は「アンモニア製造プロセス」の例題と同様なので省略する。
反応熱で発生する熱量は、フィード側とプロダクト側のエンタルピと、標準生成熱から計算する。そのために、物性計算関数 PHYHCLV( )を使って、フィード側とプロダクト側の気液混相のエンタルピを求め、PHYGTDHV( )で成分ごとの標準生成熱を得る。
気液混相のエンタルピ計算を行う PHYHCLV( )では、温度 f1.T、圧力 f1.P、組成 f1.x を与えると、液相エンタルピ hfl、液相率 rfl、液相組成 xfl、気相エンタルピ hfv、気相率 rfv、
気相組成 xfv、混相エンタルピ hf が得られる。ただし、使用するのは混相エンタルピのみである。
標準生成熱は PHYGTDHV( )を呼ぶことで、成分ごとの値が１次元配列で引数に渡される。これらを使って反応器で除熱する熱量を計算している。
＜クーラー CW01＞
出口の温度、圧力が与えられている。除熱量を計算するために、入口側のエンタルピを温度から換算する関数 PHYCLHVP( )で、出口側は気液混相のエンタルピ計算をする PHYHCLV( )
を使う。

17

/* 物質収支 */f1.x = p1.x
f1.F = p1.F/* 出口温度・圧力 */T = p1.T/* T:INPUT */P = p1.P/* P:INPUT *//* 除熱量 */PHYCLHVP( 1, f1.T, f1.P, f1.x, hf )/* 入口エンタルピ */PHYHCLV( 1, p1.T, p1.P, 1, f1.x, hpl, rl, xl, hpv, rv, xv, hp )/* 出口エンタルピ */Hf = hf*1000
Hp = hp*1000
Q = ( Hp - Hf )*p1.F/* 圧損 */
dP = p1.P - f1.P

＜フラッシュドラム FA01／FA02＞
フラッシュ計算は物性計算関数 PHYTFLSH( )を使用すれば容易に求められる。このときに与える温度、圧力は、フィード側の温度、圧力を使う。
/* 物質収支 */f1.F = p1.F + p2.F
p1.F/f1.F = rv
PHYTFLSH( 1, f1.T, f1.P, f1.x, rv, p2.x, p1.x )/* フラッシュ計算 *//* 入口温度・圧力 */f1.T = T
f1.P = P/* 出口温度・圧力 */T = p1.T
T = p2.T
P = p1.P
P = p2.P＜減圧弁 VT01＞
出口圧力のみが与えられており、出口温度を入口条件のエンタルピから PHYTCLV( )を使って計算する。
VAR xpl(_NC) ,xpv(_NC)/* 出口状態量 */18

f1.F = p1.F
f1.x = p1.x
T = p1.T
P = p1.P /* P は INPUT */PHYCLHLQ( 1, f1.T, f1.P, f1.x, hf )/* 入口エンタルピ */PHYTCLV( 1, P, hf, 1, f1.x, hpl, rpl, xpl, hpv, rpv, xpv, T )/* エンタルピ→温度 */f1.P = p1.P + dP

＜ポンプ

GA01／GA02＞

加える動力によるエンタルピ変化を考慮して、PHYTCL( )で出口の温度をエンタルピから換算する。
/* 物質収支 */f1.F = p1.F
f1.x = p1.x/* 熱収支 */PHYCLHLQ( 1, f1.T, f1.P, f1.x, hf )hp = hf + 3.6 * Power / f1.F
PHYTCL( 1, p1.P, hp, p1.x, T )/* 温度→エンタルピ *//* Power:INPUT *//* エンタルピ→温度 *//* 出口温度・圧力 */T = p1.T
P = p1.P/* P:INPUT */f1.P + dP = p1.P

＜コンプレッサ GB01＞
ポンプと同様なので省略する。ただし、動力からの変換では機械効率を加味してある。
＜合流点 MX01＞
入口側の条件から PHYCLHVP( )、PHYCLHLQ( )でエンタルピを計算し、PHYTCLV( )で出口側の温度に変換している。
VAR x(_NC)/* 物質収支 */f1.F + f2.F + f3.F + f4.F = p1.F
f1.F*f1.x + f2.F*f2.x + f3.F*f3.x + f4.F*f4.x = p1.F*x
x/SUM(x) = p1.x/* 熱収支 */PHYCLHVP( 1, f1.T, f1.P, f1.x, h1 )PHYCLHVP( 1, f2.T, f2.P, f2.x, h2 )/* 入口１エンタルピ *//* 入口２エンタルピ */19

PHYCLHLQ( 1, f3.T, f3.P, f3.x, h3 )/* 入口３エンタルピ */PHYCLHLQ( 1, f4.T, f4.P, f4.x, h4 )/* 入口４エンタルピ */ht = ( h1*f1.F+h2*f2.F+h3*f3.F+h4*f4.F )/ p1.F
PHYTCLV( 1, P, ht, 1, p1.x, hl, rl, xl, hv, rv, xv, T )/* エンタルピ→温度 *//* 入口圧力 */f1.P = P/* 出口温度・圧力 */P = p1.P
T = p1.T＜分流点

SP01／SP02＞

分流比の考え方は「アンモニア製造プロセス」の例題と同様。フィード側の温度、圧力でプロダクト側に伝える。
/* 物質収支 */f1.F = p1.F + p2.F
f1.F*f1.x = p1.F*p1.x + p2.F*p2.x
p1.F = f1.F * R/* R:INPUT */p1.F*p1.x = f1.F*f1.x * R/* 出口温度・圧力 */f1.T = p1.T
f1.T = p2.T
f1.P = p1.P
f1.P = p2.P

20

■計算結果F=2.94321
x(1)=0.882633
x(2)=0.0973261
x(3)=1.90505e-005
x(4)=0.020022

シクロヘキサン製造プロセス
F=33.0991 R
T=61.5264 s104
P=2234
1

Q=-3.62431e+006
rv=0.345391

1

GB01

F=104.352
T=200
P=2100

1
s107

BL02

BL01

1

1

F=41.83
T=37.8P=100
x(3)=1

1
2
3
4

s103

s101

1

1

1

GA01

s102

F=41.83
T=39.3969
P=2234

1

MX01

s106

1

BL03

F=68.3099
T=48.89
P=2000F=0.865573
x(1)=0.544475
x(2)=0.103739
x(3)=0.000337498
x(4)=0.351449
1

FA01

conv=0.9985
Q=3.56066e+00611

F=36.0423
T=48.89
P=2000s108

F=104.352
T=48.89
P=2000

CW01

1

s112
s111

SP01
s109

RA01
F=128.84
T=48.88
P=2234x(1)=0.99704
x(2)=0.00296

1

1
2

1

s115
F=42.3521
T=48.6757
1
P=100

VT01

2
1

s110

F=229.727
T=42.6792
P=2234BL04

rv=0.345391

1
2

1

1

s117

dP=1900

1

FA02

rv=0.0204375

s114
s113

SP02
2

1

1

R

s105

F=25.9578
T=51.4252
P=2234

GA02

s116

1

BL05

F=41.4866x(1)=0.000360483
x(2)=0.000122702
x(3)=0.000920939
x(4)=0.998596■実行の操作手順
１．エクスプローラを使って、Cycrohexane フォルダにある Cycrohex.vmp をダブルクリ
ックして VM-EQUATRAN を起動する。
２．[ファイル]-[実行]（）で実行モードに移る。

３．[設定]-[ループオーバー回数]で回数を 200 に変更する。
４．[実行]-[起動]（）で起動し、定常計算ボタン

21をクリックして収束計算を始める。