CFD Online Logo CFD Online URL
Home > Forums > CFX

CFX CEL Expression for Cavitation Modeling

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   February 14, 2013, 17:37
Default CFX CEL Expression for Cavitation Modeling
New Member
Join Date: Jun 2012
Posts: 12
Rep Power: 7
kimotbwb is on a distinguished road
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!

Kind regards
Attached Images
File Type: jpg KUBOTA_CEL_Expressions.jpg (66.3 KB, 143 views)
File Type: jpg KUBOTA_origin_Formulation.jpg (19.2 KB, 93 views)
kimotbwb is offline   Reply With Quote

Old   February 15, 2013, 06:51
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 14,068
Rep Power: 109
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
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.
ghorrocks is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
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

All times are GMT -4. The time now is 03:03.