我想在CFX中實(shí)現(xiàn)其默認(rèn)空化模型，當(dāng)作練習(xí)。但是不知道為什么，我用Expression寫(xiě)的Kubota模型老是算到一半的時(shí)候突然崩掉。
計(jì)算模型：ClarkY水翼；
CEL表達(dá)式：
   EXPRESSIONS:
      Anuc = 5*10^-4               #汽核體積分?jǐn)?shù)
      Cdest = 50                       #蒸發(fā)系數(shù)
      Cprod = 0.01                    #凝結(jié)系數(shù)
      Mdest = Cdest *3*Anuc *(1-av)*Rho v /Rb*dR     #蒸發(fā)質(zhì)量項(xiàng)
      Mprod = Cprod*3*av*Rho v/Rb*dR                     #凝結(jié)質(zhì)量項(xiàng)
      Pv = 2334                                                   #飽和壓力
      Rb = 0.000002                                             #氣泡半徑
      Rho l = 997                                        #液態(tài)水密度
      Rho v = 0.02558                                 #水蒸氣密度
      av = Vapor.Volume Fraction                                #水蒸氣體積分?jǐn)?shù)
      dR = (abs(Pv-pabsnc )*2.0/3.0/Rho l)^0.5
      m = Mdest*step((Pv-Pressure)/1 )-Mprod*step(( Pressure-Pv)/1)        #空化模型設(shè)置界面調(diào)用的質(zhì)量轉(zhuǎn)化項(xiàng)
    END
我在其余設(shè)置完全相同的情況下，分別使用CFX默認(rèn)的Kubota模型及我自己寫(xiě)的空化模型表達(dá)式進(jìn)行了計(jì)算。在后者沒(méi)有崩潰之前，兩者計(jì)算所得的各檢測(cè)量基本一致。計(jì)算到300-400步左右時(shí)，自帶模型正常，但是我自己寫(xiě)的總是突然崩潰，毫無(wú)征兆。
附圖說(shuō)明：圖1是計(jì)算域，圖2是自己寫(xiě)的模型調(diào)用界面，圖3是Kubota模型方程，圖4是報(bào)錯(cuò)內(nèi)容。
已經(jīng)卡在這很多天了，感激不盡。