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

divergerence free initial conditions for 2d flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 1, 2005, 23:05
Default divergerence free initial conditions for 2d flow
  #1
quarkz
Guest
 
Posts: n/a
Dear all,

I need to generate initial conditions for a 2d flow past an airfoil for my FVM NS solver.

What I've done is simply specify u(infinity)=1, v=0 at the left most grid for the initial conditions. The rest of the grids have zero velocity.

I've read some paper which talks about getting a divergence free velocity field. Can someone explain that?

Do I need to specify some pressure values as well?

Thank you.

  Reply With Quote

Old   August 1, 2005, 23:27
Default Re: divergerence free initial conditions for 2d fl
  #2
Ahmed
Guest
 
Posts: n/a
Thermodynamic state is defined by two independent properties, what you have is just velocity (2 components of velocity) you need the total pressure too
  Reply With Quote

Old   August 2, 2005, 04:47
Default Re: divergerence free initial conditions for 2d fl
  #3
quarkz
Guest
 
Posts: n/a
Hi,

can you explain a bit more? I assume gauge pressure to be zero. I'm not interested in the temperature and I am simulating low Re laminar incompressible flow.

In that case, how do I specify the total pressure?

Thanks

  Reply With Quote

Old   August 2, 2005, 06:51
Default Re: divergerence free initial conditions for 2d fl
  #4
Vinod Dhiman
Guest
 
Posts: n/a
Velocity fields are always required to be divergence free, since that is continuity equation for incompressible flow. So, if you've given u=constant, v=0 is right. And usually that's the initial condition for free stream flow. Is your grid starting from stagnation point, or before that. Beginning of grid should be enough upstream of the airfoil to give the condition you are giving.
  Reply With Quote

Old   August 2, 2005, 07:40
Default Re: divergerence free initial conditions for 2d fl
  #5
quarkz
Guest
 
Posts: n/a
Hi,

thanks alot. My starting pt is quite a few chords away so I believe it should be ok.
  Reply With Quote

Old   August 4, 2005, 02:53
Default Re: divergerence free initial conditions for 2d fl
  #6
Ahmed
Guest
 
Posts: n/a
Allow me to be honest with you, Could you please explain what do you mean by "gauge pressure to be zero." and how this is related to the "p" that appears in the Navier Stokes equation and called "pressure"( equilibrium pressure or thermodynamic pressure) and how are you calculating that "p" in your code. I would recommend that you review the derivation of the momentum equation in the following books 1- Continuum Mechanics by mace and mace 2- Incompressible flow by Panton 3- Boundary Layer Theory by Schlichting. Your mistake is very common among a lot of researchers and practicing engineers, I mean you are not alone in this just review other postings on this forum. The pressure that is measured by a gauge attached to the walls of your device or domain is either the total pressure or the stagnation, not the p that appears in the Navier Stokes equation (which can only be measured if your measuring device is travelling with a fluid particle in the Lagrangian Sense of description of motion).
  Reply With Quote

Old   August 17, 2005, 11:52
Default Re: divergerence free initial conditions for 2d fl
  #7
quarkz
Guest
 
Posts: n/a
hi Ahmed. in that case, can u explain a bit on what u said?

i'm solving NS in non-dimensional form. hence, the pressure is divided by p(infinity). in other words, if i multiply the p i get for solving the NS eqns by p(infinity), what is this new "pressure" then?

Thanks alot!
  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
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM 14 October 27, 2010 12:13
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21
Negative value of k causing simulation to stop velan OpenFOAM Running, Solving & CFD 1 October 17, 2008 06:36
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 04:17.