CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > ANSYS > CFX

Unwanted freezing with buoyancy

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

LinkBack Thread Tools Search this Thread Display Modes
Old   October 12, 2015, 00:33
Default Unwanted freezing with buoyancy
New Member
Join Date: May 2011
Posts: 21
Rep Power: 0
rae is on a distinguished road
I have a convective flow driven by the buoyancy of air, which is heated by a local energy source (approx. 1.1MW/m^3) whose size is defined by an additional variable within a much larger yet single domain.
The minimum temperature within the domain and opening boundaries is 298K. In my transient analysis the temperature plummets in a nonsensical manner to as low as 60K within a small region where inflow (about 4 m/s) is significant to the energy source volume. The attached image shows domain detail (bottom sloping line is adiabatic boundary and vertical line marks the energy source volume to its right)
I suspect the density drop (to about 0.2kg/m^3) due to temperature within the energy source is responsible for a similar density drop of adjoining elements outside the source and subsequently silly temperatures to some elements with significant inflow to the energy source???
I use Sutherlands formula for dynamic viscosity and thermal conductivity. I have not used tables and do not see where/how to clip the temperature to 298K where it should be and do not know if that is the best approach.
Alternatively perhaps the problem is a parallel partition problem ..., or then again, I have seen the tutorial for a steady state problem where gravity was introduced incrementally which seems like a great idea but I see no convenient way to do that during the iterations of each transient timestep????
Help out there would be most welcome. Thanks.
Attached Images
File Type: jpg T_low.jpg (139.4 KB, 11 views)
rae is offline   Reply With Quote

Old   October 12, 2015, 06:06
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,326
Rep Power: 138
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
It sounds like you better do some basic bug finding:
* Run it serial and parallel.
* Run it single precision, double precision.
* Run it to a tighter convergence, and a looser convergence.
* Run on a finer mesh and coarser mesh.
* Run it using constant properties versus Sutherland's
* Turn any other tricky physics on and off

Do any of these things affect it
ghorrocks is offline   Reply With Quote

Old   October 13, 2015, 20:09
New Member
Join Date: May 2011
Posts: 21
Rep Power: 0
rae is on a distinguished road
Thanks ghorrocks.
Further info: SAS model with B-CDS advection scheme being used and have diverge issue too -- last time the omega equation jumped to Rate 99.99 the step before failure on hydrodynamic equations.
Progress to date:
(1) Still running in serial mode, to one timestep beyond previous point of solver divergence so far, and temperatures still go very low.
(2) Running in double since before.
(3) Equation rms residuals bounce around 3E-4 to 2E-5 in 3 or 4 iterations, slowing the time step by factor 10 changes negligibly residuals, solver divergence or low temperature excursion.
(4) Expect smaller mesh will reduce but not remove temperature drop so still looking for other improvements.
(5) Change from Sutherland's to Material 'Air Ideal Gas' and 'Air @ 25C' (i.e. Bussinesq T dependence) hasn't changed issues.
(6) Running Expert: max continuity loops = 2, since before as found to help earlier (curiously RMS > max residual for 2nd P-mass loop occasionally); tef numerics option = 1, has allowed simulation to progress without solver divergence; and highres energy = 3, also allows simulation to progress. Still gets too cold for my liking.
rae is offline   Reply With Quote


buoyancy, convection, freezing, transient

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
How to realize Boussinesq approximation (buoyancy) in Fluent-UDF dochanxu FLUENT 1 August 12, 2022 05:26
Buoyancy flow sunilpatil CFX 1 December 17, 2014 18:30
Freezing of an Ice Cream 10mmet21 Main CFD Forum 0 February 28, 2012 07:15
Buoyancy Turbulence Chilli83 CFX 3 April 14, 2009 20:28
Buoyancy production for species concentration George Gerber FLUENT 1 May 25, 2006 16:56

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