気相触媒反応層の温度分布 表計算
Excelで連立常微分方程式に挑戦！ルンゲクッタ法で積分をおこなう。プログラムはホームページプラグミングと同じ。微分方程式の記述は以下のよう。
Sub SUB2000()
'---------ﾊﾟﾗﾒｰﾀ
F0 = 40000
ZA0 = 0.05
ZB0 = 0.2
P = 101.3
DELTAA = (1 - 1 - 1) / 1
CP = 31.4
Q = -142400
TW = 673
ROB = 1100
UAS = Cells(1, 2)
PA = ZA0 * (1 - Y(0)) * P / (1 + DELTAA * ZA0 * Y(0))
PB = (ZB0 - ZA0 * Y(0)) * P / (1 + DELTAA * ZA0 * Y(0))
K = 24400 * Exp(-75370 / 8.314 / Y(1))
R = -K * PA * PB
'---------ﾋﾞﾌﾞﾝ ﾎｳﾃｲｼｷ
Select Case NK
Case 0
' [ DYDX=Y(0)' ]
DYDX = -R * ROB / (ZA0 * F0)
Case 1
' [ DYDX=Y(1)' ]
DYDX = (Q * R * ROB - UAS * (Y(1) - TW)) / (CP * F0)
End Select
End Sub
