我想在CFX中实现其默认空化模型，当作练习。但是不知道为什么，我用Expression写的Kubota模型老是算到一半的时候突然崩掉。
计算模型：ClarkY水翼；
CEL表达式：
   EXPRESSIONS:
      Anuc = 5*10^-4               #汽核体积分数
      Cdest = 50                       #蒸发系数
      Cprod = 0.01                    #凝结系数
      Mdest = Cdest *3*Anuc *(1-av)*Rho v /Rb*dR     #蒸发质量项
      Mprod = Cprod*3*av*Rho v/Rb*dR                     #凝结质量项
      Pv = 2334                                                   #饱和压力
      Rb = 0.000002                                             #气泡半径
      Rho l = 997                                        #液态水密度
      Rho v = 0.02558                                 #水蒸气密度
      av = Vapor.Volume Fraction                                #水蒸气体积分数
      dR = (abs(Pv-pabsnc )*2.0/3.0/Rho l)^0.5
      m = Mdest*step((Pv-Pressure)/1 )-Mprod*step(( Pressure-Pv)/1)        #空化模型设置界面调用的质量转化项
    END
我在其余设置完全相同的情况下，分别使用CFX默认的Kubota模型及我自己写的空化模型表达式进行了计算。在后者没有崩溃之前，两者计算所得的各检测量基本一致。计算到300-400步左右时，自带模型正常，但是我自己写的总是突然崩溃，毫无征兆。
附图说明：图1是计算域，图2是自己写的模型调用界面，图3是Kubota模型方程，图4是报错内容。
已经卡在这很多天了，感激不尽。