The replication of Zwart cavitation model using CEL languages
 

I want to reproduce the mass transfer rate of Zwart cavitation model using the CEL language. the following languages is CEL which were used to replicate the cavitation model of Zwart.

If = if(pabsnc>=3574[kg m^-1 s^-2], pabsnc-3574[kg m^-1 s^-2], 3574[kg m^-1 s^-2]-pabsnc)
Rob = 0.000001[m]
F1 = sqrt((2/3)*(If /Water.density))
Cond = 0.01*3*Vapor.vf*Vapor.densitync*F1 /Rob
Vap = 50*3*0.0005*(1-Vapor.vf)*Vapor.densitync*F1 /Rob
Mass transfer = Vap -Cond

but when i do the solve-process, the computation will diverage quickly. i do not know why. please give me some help. thanks a lot

First of all, does the simulation run correctly when you use the built in cavitation model?

If that is OK then do some basic debugging:

Send each of these variables to a user variable, defined to be CEL and set to the value of each of your expressions. Change your output file to write the results file every iteration. I would also add the residuals to the results file as that might have useful information as well.

Then rerun the simulation until it crashes. Look at the results file and check all your variables are the expected values.