# variable density in single fluid simulation not running with total energy

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

 September 18, 2012, 06:34 variable density in single fluid simulation not running with total energy #1 New Member     Sebastian Join Date: Jul 2010 Location: Cologne Posts: 10 Rep Power: 7 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

 September 18, 2012, 20:17 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 10,923 Rep Power: 85 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.

 September 18, 2012, 20:19 #3 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 10,923 Rep Power: 85 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Paresh Jain CFX 17 March 11, 2014 12:40 bhushas Main CFD Forum 1 May 30, 2008 04:35 Harry Dong Main CFD Forum 12 February 4, 2006 01:55 Freeman Main CFD Forum 10 January 27, 2006 08:42 Zhong Lei Main CFD Forum 23 May 14, 1999 13:22

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