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

variable density in single fluid simulation not running with total energy

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 18, 2012, 06:34
Default variable density in single fluid simulation not running with total energy
  #1
New Member
 
viking's Avatar
 
Sebastian
Join Date: Jul 2010
Location: Cologne
Posts: 10
Rep Power: 7
viking is on a distinguished road
Dear CFX-Users,

i have a problem simulating a single fluid case with variable density. i used an expression in CEL for modelling the density. the density decreases for low pressure values and should represent the increasing bubble friction (desorption of non condensable gases) in areas of pressure below the atmospheric pressure.

in my simulation model there is couette and poseuille flow through a gap (rotor-stator). the boundary conditions are total pressure at the inlet and opening entrainment at the outlet.

i allready have done some other simulations with these boundaries and the total energy model but different fluid propperties as variable densities on multiphase flow. they did all work. for the single fluid simulation the solver writes an error "divide by zero" and even an initialisation from a constant density simulation does not solve the problem for total energy model. i also tried some oter timesteps.

furthermore the obove described CEL density expression workes fine with an thermal energy model or isotherm as well.

can anyone tell me, why i get the divide by zero error?
the error does not occure whithin my CEL expressions, in this case the solver would tell me the problematic CEL.

kind regards
yours viking
viking is offline   Reply With Quote

Old   September 18, 2012, 20:17
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,923
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
When you put weird kinks and bends in the density curve that always makes convergence difficult. Your problem is simply difficult numerics (assuming there is no gross error in your equations).

One hint is to make sure that the density is continuous against pressure (I assume you have density as a function of pressure), and preferably keep the first derivative continuous as well.
ghorrocks is offline   Reply With Quote

Old   September 18, 2012, 20:19
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,923
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
And also consider using the cavitation model for this as well. The physics are quite similar so you may be able to adapt the cavitation model to this. That model has stuff built in to handle the convergence difficulties.
ghorrocks 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
Please Help.....Insufficient Catalogue Size Paresh Jain CFX 17 March 11, 2014 12:40
error in COMSOL:'ERROR:6164 Duplicate Variable' bhushas Main CFD Forum 1 May 30, 2008 04:35
What is the total energy for incompressible fluid? Harry Dong Main CFD Forum 12 February 4, 2006 01:55
Energy dissipation and Drag coefficient Freeman Main CFD Forum 10 January 27, 2006 08:42
Why FVM for high-Re flows? Zhong Lei Main CFD Forum 23 May 14, 1999 13:22


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