|February 14, 2013, 17:37||
CFX CEL Expression for Cavitation Modeling
Join Date: Jun 2012
Posts: 12Rep Power: 6
I want to modify the cavitation model used in CFX by adding some specific user defined mass transfer soure terms for vaporization and condensation. This can be done by adding some expressions using CEL or by user CEL function which uses a Fortran subroutine. In my case i am using some simple CEL expressions.
What I want do to, is to adopte the so caled KUBOTA-model (sometimes also caled ZWART-model).
The formulation is based on a modification of the Rayleigh-Plesset equation for individual bubbles.
My problem is, to correctly adopte the KUBOTA-model in terms of CEL expressions, as some equations reproduce a negative sqrt! so the solver crashes.
You can see the entire fomulation and the used CEL expressions in the attached file. The basic equation which is used as the user defined cavitation model, is called KUBOTA.
As you can see, the two equations "PTermVap" and "PTermCond" can obtain negative values depending on if vaporization (p < pvap) or if condensation (p > pavp) occurs.
Did i adopte the fomulations correctly ? How can i avoid that the sqrt gets negative values ?
The original formulation is also attached as a *.JPG. So you can see the entire equations and compare it with the KUBOTA equation that I am trying to implement.
I also used the variable "pasbnc" instead of "Pressure" as the modeling guide suggests that the pressure in (pvap-p) should be the nonclipped absolute pressure (variable "pabsnc").
Hopefully somebody can help. Thanks a lot!
|February 15, 2013, 06:51||
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,559Rep Power: 97
You are evaluating both PTermVap and PTermCond, and at all times at least one of them will be undefined with a sqrt of a negative number. So put them in your if statement for the KOBUTA function and then only the valid one gets evaluated.
|Thread||Thread Starter||Forum||Replies||Last Post|
|problem in modeling thin wedge in cfx||m.salemi||CFX||1||July 25, 2011 18:49|
|Help needed for writing CEL expression||kmgraju||CFX||5||September 20, 2010 18:30|
|CEL in CFX??||Vishal||CFX||2||April 5, 2008 11:48|
|Help to create a CEL Expression||Juan Maria Campos||CFX||2||December 3, 2007 08:52|
|Lift, Drag Vs time chart,calculations||Jamesd69climber||CFX||8||February 17, 2005 18:23|