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

Atmospheric stability

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 10, 2005, 17:27
Default Atmospheric stability
  #1
CFD-student
Guest
 
Posts: n/a
Here is a though one (at least for me) I been trying to model hydrogen dispersion within a neutraly stratified atmosphere and I been having some trouble for more than amonth. First I tryed steady-state; this of course had convergance trouble. So I'm now trying to set-up a transient model. I'm having start-up problems. Here are the specs of my model definition. Any feedback would be inmensly apreciated.

•Domain size; L,W,H 750,350,250 meters

•Turbulent Flow: k-epsilon , 0.01-0.0002 (from ref)

•Inlet: U=7m/s and Temperature profile T/Z -2.45 K

•Density: Ideal(T)

•Hydrogen being model as a scalar from a pipe in center of domain with W=1m/s.

•Bouyancy; ON (having doubts on datum definition need help)

•Algorithm; PISO and Implicit Time-Disc

•Initialization from steady-state (no H2) run

•Time step 0.1s (having trouble defining this)

To whom ever responds... THANKS A LOT!

haoa
  Reply With Quote

Old   May 10, 2005, 17:31
Default Re: Atmospheric stability
  #2
CFD-student
Guest
 
Posts: n/a
I forgot to add this:

Boundary conditions;

•Inlet at east

•Outlet (split=1) at west

•No-Slip wall at bottom

•Symetry sides and top

USING STAR-CD

Thanks!
  Reply With Quote

Old   May 10, 2005, 20:23
Default Re: Atmospheric stability
  #3
Jim_Park
Guest
 
Posts: n/a
A couple of things you might want to look at:

I'd start with a time step that is 0.1 of the minimum cell transit time, the time it takes a packet of fluid to cross a cell.

First find min over all cells of (dx/u, dy/v, dz/w), then set your dt as 0.1 of that. The code should eventually take over and adjust (grow) the time step as needed...

You might check your references for the velocity distribution across an atmospheric boundary layer. It can be several meters thick. This would change your upstream boundary condition.

Good luck!
  Reply With Quote

Old   May 11, 2005, 02:39
Default Re: Atmospheric stability
  #4
Harry Fulmer
Guest
 
Posts: n/a
As Jim mentioned, you'll need to specify the atmospheric boundary layer velocity profile. A quick google has provided:

http://personalpages.umist.ac.uk/sta...bl/regions.pdf

Look at the section containing equation 16.

You'll need similar profile descriptions for other turbulence quantities.

More quick and lazy googling finds:

http://www.dnv.com/binaries/Applicat..._tcm4-8313.pdf

I'm sure there's lots more out there.

What I've done in the past is to start the model simple, ensure that the defined ABL velocity profile does not degrade unduly as it passes over an 'empty' flat grounded domain. To support the sheer size of the ABL (sometimes up to 300m high) you should investigate the application of the correct roughness height to the ground domain plane. Once you're happy with the ABL then introduce the contaminant source into your model.

If you've got a point source of hydrogen then you'll need bucket loads of grid and a nice high order scheme to ensure that the intial spreading is not swamped by good old numerical diffusion.

Good luck, it's tough old problem but one that you should find lots of reference to in the literature. Hey, maybe the guys at CD support will be also be able to give appropriate help and advice!
  Reply With Quote

Old   May 11, 2005, 12:49
Default Re: Atmospheric stability
  #5
CFD-student
Guest
 
Posts: n/a
Jim and Harry; Thanks for your feedback.

I used a uniform approach velocity with the idea that the BL would be fully developed by the time it got to my hydrogen source. (this is not of the same size of the actuial (ABL). But changing it to a log vel profile is a good idea.

•Do you think that the inlet vel-pprofile is truly critical?

Also my mesh is tetrahedral and is refined considerably (from 15m far from source to 4m near the source) keepping the size changes under 20%.

•The high order scheme that you suggest; do you think QUICK for all variables should work, or only for my scalar (H2?)

Thanks A LOT.
  Reply With Quote

Old   May 12, 2005, 09:39
Default Re: Atmospheric stability
  #6
WT
Guest
 
Posts: n/a
Hello,

I fully agree with Harry regarding the neccessity of having good start profiles.

Way back at my time at the university, I was dealing with atmospheric flows. I was calculating the wake of wind turbines wrt. atmospheric stratification. I remember that it was very important to have profiles of velocity, k and eps corresponding to atmosperic bl type profiles as boundary conditions, otherwise the effect of creating the bl profile was mixed with the effect of the wind turbine wake. We created the profiles by doing a 2-d simulation of a very (I remember about 2km long or so) long flat plate with appropriate roughness and stability conditions set.

Btw, the profiles are only valid in the lower part of the atmosphere (depending on stability, maybe about 200m above ground), above we assumed a constant velocity with very low turbulence level. A good starting point for information on atmospheric flow and tubulence is the book "Atmospheric turbulence" by Dutton, Panofsky and maybe "A first course in turbulence" by Tennekes, Lumley.

Hth,

-Thomas
  Reply With Quote

Old   May 12, 2005, 11:34
Default Re: Atmospheric stability
  #7
CFD-Student
Guest
 
Posts: n/a
Thanks for all your comments... I have the articles and the text books... I will revisit them.

Hope to payback the favor soon!

thnx!

Hector A. Olvera
  Reply With Quote

Old   May 13, 2005, 17:49
Default Re: Atmospheric stability (convergance)
  #8
CFD-student
Guest
 
Posts: n/a
I need some help:

I'm running the H2 dispersion in a stratified atmosphere; the scalar equation (H2) is not converging at 80 corrector steps and it has negative values (extremly small ones x*EE-175). Does any one have any suggestions why this could be hapening?

thnx!
  Reply With Quote

Old   May 14, 2005, 03:27
Default Re: Atmospheric stability (convergance)
  #9
Harry Fulmer
Guest
 
Posts: n/a
What numerical scheme are you using? Is it bounded?
  Reply With Quote

Old   May 16, 2005, 19:51
Default Re: Atmospheric stability (convergance)
  #10
CFD-student
Guest
 
Posts: n/a
I'm using PISO I have and Inlet and Outlet boundaries in east and west. no-slip on ground and symetry all the rest...

i the residuals are 0.99 approx at 30 corrector steps
  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
stability and number of processors ?? sersunzo OpenFOAM Running, Solving & CFD 0 June 9, 2009 11:45
Stability Condition for Immersed Boundary Method Dave Rudolf Main CFD Forum 0 February 6, 2007 00:36
How to generate a atmospheric boundary layer Morten Andersen CFX 3 January 16, 2007 07:48
Hydrodynamic Stability vs CFD Ben Main CFD Forum 2 October 18, 2004 02:39
Solving linear stability eq. K S Chang Main CFD Forum 3 January 19, 2004 17:01


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