プログフロー反応器の温度分布 ホームページプログラミング
連立常微分方程式を解く。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