プログフロー反応器の温度分布 ホームページプログラミング 


連立常微分方程式を解く。diffeq3cat.html

スクリプト部分の記述

(VBScript の微分方程式の部分)
Sub SUB2000
'--微分方程式 変数の割り当て:転化率 x → Y(0) , 温度 T →Y(1)
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