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

Buoyancy driven flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 14, 2012, 09:07
Arrow Buoyancy driven flow
  #1
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
Hi everybody,

I'm going to solve a simple flow over flat plate and I'm going to have heat equation and buoyancy effects in my flow. My domain size is 1500(height)x6000x6000 [m] and my first grid center is located 0.03 away from the ground. I have periodic lateral BC, Symmetry top BC and no-slip wall BC.

I can simulate the flow when I have a neutral case but with temperature it diverges very fast. Do you have any idea what could be wrong?

Note that with this size Grashof number becomes enormous, but I don't know if that's the problem.

I appreciate all your suggestions
chivaee is offline   Reply With Quote

Old   March 14, 2012, 14:07
Default
  #2
New Member
 
Join Date: Jun 2009
Posts: 9
Rep Power: 8
sgrshukla is on a distinguished road
How about first trying with reducing your URF and making your solution more stable.
sgrshukla is offline   Reply With Quote

Old   March 14, 2012, 14:26
Default
  #3
Senior Member
 
Andrew
Join Date: Mar 2009
Location: Washington, DC
Posts: 194
Rep Power: 8
mettler is on a distinguished road
it sounds like a grid problem. Especially with buoyancy driven flow. Most of the interaction will be close to the wall, so you need a good grid in that region. I am assuming your plate is vertical and not horizontal.
mettler is offline   Reply With Quote

Old   March 15, 2012, 08:02
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,581
Rep Power: 20
FMDenaro will become famous soon enough
maybe, it can be a case of numerical instability due to the high Gr number...
FMDenaro is offline   Reply With Quote

Old   March 15, 2012, 10:03
Default
  #5
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
Quote:
Originally Posted by sgrshukla View Post
How about first trying with reducing your URF and making your solution more stable.
I see your point, but I've tried with URFs of down to 0.15 with no success ...
chivaee is offline   Reply With Quote

Old   March 15, 2012, 10:11
Default
  #6
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
Quote:
Originally Posted by mettler View Post
it sounds like a grid problem. Especially with buoyancy driven flow. Most of the interaction will be close to the wall, so you need a good grid in that region. I am assuming your plate is vertical and not horizontal.
Hmmmm it could be. I've used a simple tanh function to stretch in the wall normal direction and uniform in the other two directions ... the plate is horizontal, in fact its an attempt to model atmospheric boundary layer simulation with 1.5 km height ... my issue is that the stretching rate is more than 30% but it looks quite smooth.

Now, even the temperature equation works fine but only if I disregard the gravitation.
chivaee is offline   Reply With Quote

Old   March 15, 2012, 10:13
Default
  #7
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
I might also take a look at the block communications (it's a multi-block structured code I'm using) to see if something goes strange ...
chivaee is offline   Reply With Quote

Old   March 15, 2012, 10:15
Default
  #8
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
maybe, it can be a case of numerical instability due to the high Gr number...
mmm, yeah, meybe, since it scales up as L^3 (L being characteristic length scale) ... but isn't it what normally happens? I mean I'm just simulating a small part of the real atmosphere with temperature stratification ...
chivaee is offline   Reply With Quote

Old   March 19, 2012, 11:39
Default problem partly solved, not sure if it's true ...
  #9
New Member
 
Sarlak
Join Date: Mar 2012
Posts: 17
Rep Power: 5
chivaee is on a distinguished road
Hi again,
regarding the issue with buoyancy driven flow, I tried to decouple temperature and momentum equations completely i.e. after the solution for pressure and velocities converge, I solve the heat equation and most importantly "add the buoyancy contributions".

Seems it solved the issue ...

I know it's a general practice to treat the heat equation like this (quasi-steady) but not sure if it's true when you have buoyancy? since buoyant force affects the flow field directly and might make it inaccurate ...
chivaee is offline   Reply With Quote

Old   March 19, 2012, 11:56
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,581
Rep Power: 20
FMDenaro will become famous soon enough
Quote:
Originally Posted by chivaee View Post
Hi again,
regarding the issue with buoyancy driven flow, I tried to decouple temperature and momentum equations completely i.e. after the solution for pressure and velocities converge, I solve the heat equation and most importantly "add the buoyancy contributions".

Seems it solved the issue ...

I know it's a general practice to treat the heat equation like this (quasi-steady) but not sure if it's true when you have buoyancy? since buoyant force affects the flow field directly and might make it inaccurate ...
are you searching for the steady solution? then your method can be quite good, otherwise you need more accuracy in time, for example doing a predictor-corrector integration
FMDenaro is offline   Reply With Quote

Old   April 19, 2012, 04:46
Default show me by CFX
  #11
Senior Member
 
Safia
Join Date: Oct 2010
Location: Australia
Posts: 161
Rep Power: 6
happy is on a distinguished road
Quote:
Originally Posted by chivaee View Post
Hi again,
regarding the issue with buoyancy driven flow, I tried to decouple temperature and momentum equations completely i.e. after the solution for pressure and velocities converge, I solve the heat equation and most importantly "add the buoyancy contributions".

Seems it solved the issue ...

I know it's a general practice to treat the heat equation like this (quasi-steady) but not sure if it's true when you have buoyancy? since buoyant force affects the flow field directly and might make it inaccurate ...
could you please show me how you could do it by CFX?
My problem is the same as yours. I have fluid flow with heat transfer and the buoyancy takes place through my model.
Regards
happy is offline   Reply With Quote

Old   May 14, 2012, 23:56
Default
  #12
Senior Member
 
Safia
Join Date: Oct 2010
Location: Australia
Posts: 161
Rep Power: 6
happy is on a distinguished road
Quote:
Originally Posted by happy View Post
could you please show me how you could do it by CFX?
My problem is the same as yours. I have fluid flow with heat transfer and the buoyancy takes place through my model.
Regards
finally, I got it!!
thanks for your good thread.However, I'm still sad because I asked to show me how until I do it after 2 months!!..some guys do not care about the others. you will feel better whan youbecome helpful member
happy is offline   Reply With Quote

Reply

Tags
buoyancy flow, divergence, grashof number

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
Flow meter Design CD adapco Group Marketing CD-adapco 3 June 21, 2011 08:33
is there any parallel code for the famous Lid Driven Cavity flow? gholamghar Main CFD Forum 0 August 1, 2010 01:55
2D Lid Driven Cavity Flow simulation using MATLAB josephlm Main CFD Forum 3 June 24, 2010 01:58
steam flow in a pipe driven by a pressure gradient between inlet and outlet SalvoCalvo Main CFD Forum 0 March 11, 2010 07:52
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


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