CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

CFX Cavitation Model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 26, 2013, 05:38
Default CFX Cavitation Model
  #1
New Member
 
Join Date: Jun 2012
Posts: 12
Rep Power: 5
kimotbwb is on a distinguished road
Hi there,

i have a couple of questions regarding the CFX Cavitation Model.
CFX uses the simplified Rayleigh-Plesset model per default. The formulation of the source term (transport based equation model) is done by Zwart (this model is similar to the Kubota formulation)
I tried to implement the MERKLE model by using CEL expressions, but stability problems occur. I even tried to implement the original formulation used by CFX, but even there i got stability problems. The solver fails and produces a fatal overflow in the linear equations after a couple of iterations.
Does anybody has some experiences with the implementation of user defined cavitation models ?
Perhaps somebody could sent me his CEL-formulation as a template?

king regards!
kimotbwb is offline   Reply With Quote

Old   February 26, 2013, 17:58
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,646
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Cavitation models are always difficult to get stable as they involve such large changes in density.

Implementing a model in CEL can add additional stability problems as the model is not implicitly included in the numerics but is updated separately. This means the implementation of existing models with CEL will be less stable (as you found).

So I would try user fortran to implement your model. I am not sure but hopefully this is better linked into the solver and will result in better numerical stability.
ghorrocks is offline   Reply With Quote

Old   February 26, 2013, 18:24
Default
  #3
New Member
 
Join Date: Jun 2012
Posts: 12
Rep Power: 5
kimotbwb is on a distinguished road
Thanks for your reply. We have already thought of this. But the first impression is that there isn't any relative difference regarding the stability problems elsewhere as for CEL and for Fortran routine. We tryied to reproduce the original formulation for the source and the sink term that CFX model uses both with CEL and Fortran, but we obtained stability problems in both cases.
There are some literature references available from authors who implemented the KUBOTA and the MERKLE model for their purposes and got stable and accurate results.
We tried to stabilize the CEL and the Fortran routine by some tricks like if/else statements for condensation and for vaporization, but actualy the linear solver exceeds the bounds and fails out with an overflow.

The funny think is, that we even tryied to reproduce the original CFX model (see KUBOTA) on a working testcase. The original CFX model gets very good results in a close range depending on the cavitation number, but we failed trying to reproduce the same testcase by using CEL and fortran routine.

It would be very interesting to know, if somebody has got a CEL-template or maybe a routine.f with a working and stable user defined cavitation model.
kimotbwb is offline   Reply With Quote

Old   February 27, 2013, 06:51
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,646
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Oh yes, and did you use source term coefficients with your source terms?
ghorrocks is offline   Reply With Quote

Old   February 27, 2013, 07:51
Default
  #5
New Member
 
Join Date: Jun 2012
Posts: 12
Rep Power: 5
kimotbwb is on a distinguished road
yes of cource. we used the default coefficients of the CFX source term model. means: Fvap=50.0 and Fcond=0.01.
did you ment that?
kimotbwb is offline   Reply With Quote

Old   February 27, 2013, 17:50
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,646
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
No. When you use a source term you are recommended to define a linearisation coefficient to assist with convergence. The source term coefficient is general taken as the derivative of the source term WRT to equation variable. In several cases I have done I could not get a source term to converge, but the addition of a source term coefficient made it converge quickly and easily. It makes a massive difference.

Have a look in the documentation about source terms, especially the linearisation coefficients.
ghorrocks is offline   Reply With Quote

Old   June 18, 2013, 02:08
Default
  #7
New Member
 
aerodung
Join Date: Apr 2013
Location: Canada
Posts: 15
Rep Power: 4
hydroaero is on a distinguished road
Hi kimotbwb,
Did you finish the implementation of Merkle model in CFX? Your results are good?
I plan to implement Merkle model in CFX by User Fortran soon.
Any advice from your part would be appreciated!
See you!
hydroaero is offline   Reply With Quote

Old   July 10, 2013, 13:32
Default
  #8
New Member
 
aerodung
Join Date: Apr 2013
Location: Canada
Posts: 15
Rep Power: 4
hydroaero is on a distinguished road
Dear Glenn Horrocks,
I’m implemeting new cavitation models using CFX user routine with Fortran code.
The model is running but doesn’t converge.
I want to use linearization coefficients to improve the convergence.
I read the CFX documentation about linearization coefficients for source terms.
As you said, linearization coefficient = derivative of source term to equation variable.
Do you know where to specify these linearization coefficients for cavitation modelling in CFX?
Thanks for your help!
See you soon!

Quote:
Originally Posted by ghorrocks View Post
No. When you use a source term you are recommended to define a linearisation coefficient to assist with convergence. The source term coefficient is general taken as the derivative of the source term WRT to equation variable. In several cases I have done I could not get a source term to converge, but the addition of a source term coefficient made it converge quickly and easily. It makes a massive difference.

Have a look in the documentation about source terms, especially the linearisation coefficients.

Last edited by hydroaero; July 11, 2013 at 13:37.
hydroaero is offline   Reply With Quote

Old   July 13, 2013, 07:20
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,646
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
I think you need to talk to CFX support for this type of detail.
ghorrocks is offline   Reply With Quote

Old   November 12, 2013, 22:20
Default
  #10
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 263
Rep Power: 4
huangxianbei is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
Cavitation models are always difficult to get stable as they involve such large changes in density.

Implementing a model in CEL can add additional stability problems as the model is not implicitly included in the numerics but is updated separately. This means the implementation of existing models with CEL will be less stable (as you found).

So I would try user fortran to implement your model. I am not sure but hopefully this is better linked into the solver and will result in better numerical stability.
Hi
I'm trying to implement the ZGB cavitation model just to ensure the way through Fortran is right. But I'm confused about the volume fraction definition of the vapor, for the vapor volume fraction is the main variable in the cavitation model. The user_get_var or use_get_gvar seems powerless in getting the vapor volume fraction data. Could you make some suggestion about how to get this variable?
huangxianbei is offline   Reply With Quote

Old   November 12, 2013, 22:24
Default
  #11
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 263
Rep Power: 4
huangxianbei is on a distinguished road
Quote:
Originally Posted by hydroaero View Post
Dear Glenn Horrocks,
Iím implemeting new cavitation models using CFX user routine with Fortran code.
The model is running but doesnít converge.
I want to use linearization coefficients to improve the convergence.
I read the CFX documentation about linearization coefficients for source terms.
As you said, linearization coefficient = derivative of source term to equation variable.
Do you know where to specify these linearization coefficients for cavitation modelling in CFX?
Thanks for your help!
See you soon!
Hi
I'm trying to implement a cavitation model into CFX too. But the problem I'm facing is the variable-vapor volume fraction. Can you give me the guidance about how to get the vapor volume fraction data during the calculation?
huangxianbei is offline   Reply With Quote

Reply

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
"Full" buoyancy model in CFX Jade M Main CFD Forum 2 March 18, 2010 13:08
CFX Liquid Evaporation Model Jinx CFX 3 January 28, 2010 17:31
v2-f turbulence model in CFX? flga CFX 14 November 23, 2006 07:12
software new version CAVITATION model! ROOZBEH CFX 1 October 29, 2005 12:32
Use of 1 equation turbulence model in CFX 4.3 Niels Deen CFX 0 July 19, 2000 08:50


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