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

The 'other' half of the continuity equation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 2, 2007, 10:32
Default The 'other' half of the continuity equation
  #1
desA
Guest
 
Posts: n/a
I've been having a few thoughts about how we formulate the 'incompressible Navier-Stokes' system & would value a few insights.

Background:

1. We have 3 momentum equations - x,y,z cartesian directions

2. We begin with the transient form of the continuity eqn

(all derivs are partials, r=density, N = nabla, v = velocity, D = substantial deriv):

dr/dt + v.Nr + r*(N.v) = 0

or

dr/dt + v.Nr + r*div(v) = 0

massage into form:

div(v) = -(1/r)*(dr/dt + v.Nr) = -Dr/r

For incompressible flow, we claim that -Dr/r, & discard this to end up with:

div(v) = 0

our well-known 'friend'...

3. We use div(v)=0 to drop the velocity cross-term form the compressible momentum equations, with some manipulation.

4. We now form a system of 4 equations, re-using div(v)=0:

3 x fn(u,v,w,p) -> 3 incompressible momentum eqns

1 x fn(u,v,w) -> div(v)=0

------------------

A few provocative thoughts:

1. Is it a 'wise' idea to try use a closing equation deficient in p? Would it not be preferable to use an equation which contains p? We do not have 4 equations in 4 unknowns, but rather, 3 equations in 4 unknowns & one semi-closing equation (deficient in p).

2. In our restriction to obtain div(v)=0, without making specific conditions for all terms of Dr/r to be zero, have we not created a 'density wave' that is going to enter our solution space?

3. Further to point 2, if ' r=c*P ', as is the case for almost incompressible fluids, with 'c' being very, very small, we have then created a 'pressure wave' in our solution space.

4. With so-called 'steady solvers', are we not merely substituting 'time' with 'iteration number'?

5. It is interesting that the pressure 'blowup' is an often seen malady in cfd solutions.

-----------------

A little food for thought.

desA
  Reply With Quote

Old   February 2, 2007, 12:05
Default Re: The 'other' half of the continuity equation
  #2
andy
Guest
 
Posts: n/a
> I've been having a few thoughts about how we formulate the
: 'incompressible Navier-Stokes' system & would value a few insights.

The assumption of incompressibility does not necessarily mean density=constant. There are many incompressible flow combustion codes with strong density gradients.

Incompressibility is normally invoked by splitting the pressure the into two components: a constant in space but possibly time varying background thermodynamic pressure and a small spatially varying component that is not involved in determining the fluid state but does enter in the momentum balance. Because the thermodynamic pressure is constant the fast moving acoustic waves have been eliminated which is usually the reason for invoking the assumption of incompressibiliy.

If the flow is also assumed to be isentropic then the energy equation will have nothing useful to contribute and can be dropped. With a constant thermodynamic pressure and constant entropy then the density will also be constant.

  Reply With Quote

Old   February 2, 2007, 14:17
Default Re: The 'other' half of the continuity equation
  #3
desA
Guest
 
Posts: n/a
Thanks Andy for your input - it is very much appreciated.

desA wrote:
: I've been having a few thoughts about how we formulate the
: 'incompressible Navier-Stokes' system & would value a few insights.


Andy replied:
:The assumption of incompressibility does not necessarily mean density=constant. There are many incompressible flow combustion codes with strong density gradients.

Incompressibility is normally invoked by splitting the pressure the into two components: a constant in space but possibly time varying background thermodynamic pressure and a small spatially varying component that is not involved in determining the fluid state but does enter in the momentum balance. Because the thermodynamic pressure is constant the fast moving acoustic waves have been eliminated which is usually the reason for invoking the assumption of incompressibiliy.

If the flow is also assumed to be isentropic then the energy equation will have nothing useful to contribute and can be dropped. With a constant thermodynamic pressure and constant entropy then the density will also be constant.


desA's reply:
:Now this is very interesting. Can I ask in which field you work - it would appear to be combustion judging form your early statement? I find the approach you're using to make complete sense.

Now, I've seen a number of ideas around this incompressibility issue, including:

density = const wrt x,y,z,t => div(v)=0 is direct

constant fluid properties - no mention of spatial, or temporal.

Also mention of a compressibility factor relating to density change. Less that 30% being practically considered 'incompressible'.

Many of the engineering fluids & convection books are written this way. The mathemagicians seem to often simply use div(v)=0 as the guarantee of incompressibility. At very slow flows, would we still expect to see a pressure-wave? I'm talking about as low at the onset of instability, say Re 40 odd over a cylinder?

Thanks again Andy, for your input into the debate.

desA

  Reply With Quote

Old   February 3, 2007, 22:25
Default Re: The 'other' half of the continuity equation
  #4
Ahmed
Guest
 
Posts: n/a
I just would like to add a simple clarification: Incompressibility means the effect of pressure on the density is small enough to be considered negligible (in Engineering Applications where the Mach number is less than 0.3 this has be proven experimentally), This does not exclude the effect of Temperature on density which Andy is referring to from his observations on combustion processes. Diaw, The incompressible Navier Stokes equations are an Engineering Assumption that is nice and good for design purposes only, but not for a research of the wave phenomnun you are talking about. Computer codes blow up because of rounding errors not because of the lack of pressure effects on the density.
  Reply With Quote

Old   February 4, 2007, 06:06
Default Re: The 'other' half of the continuity equation
  #5
desA
Guest
 
Posts: n/a
Welcome, Ahmed - great to have you aboard...

Ahmed wrote:

I just would like to add a simple clarification: Incompressibility means the effect of pressure on the density is small enough to be considered negligible (in Engineering Applications where the Mach number is less than 0.3 this has be proven experimentally), This does not exclude the effect of Temperature on density which Andy is referring to from his observations on combustion processes.


desA replies: I agree with this observation. So, how would engineers implement the 'continuity' closure equation, in a way that best expresses their observations?

Ahmed wrote:

Diaw, The incompressible Navier Stokes equations are an Engineering Assumption that is nice and good for design purposes only, but not for a research of the wave phenomnun you are talking about. Computer codes blow up because of rounding errors not because of the lack of pressure effects on the density.


There is a lot of current research going into these 'blow-up' solutions at present. Something is definitely there, in the physics, beyond the early numerical roundoff argument. Ironically with my current work, it is a very bad day should I come cross a numerical 'blowup'& I really torture the N-S.

I have managed to isolate what I call 'momentum waveforms' in the N-S equations. I can show them via tensor & vector calculus - as well as in various simulations. In current work, I've been able to predict, via tensor calculus, where the solution bifurcations should occur. These waveforms are there alright - & in reasonably large dimension, especially at Reynolds numbers climb, with almost negligible compressibility & bulk modulus in the order of that for water.

My current debate is what is the 'correct' continuity closure equation to use which can serve physical, mathematical & computing requirements.

desA
  Reply With Quote

Old   February 6, 2007, 08:07
Default Re: The 'other' half of the continuity equation
  #6
bystander
Guest
 
Posts: n/a
Why don't you publish your findings?
  Reply With Quote

Old   February 6, 2007, 08:36
Default Re: The 'other' half of the continuity equation
  #7
desA
Guest
 
Posts: n/a
bystander wrote:

Why don't you publish your findings?

desA replies:

Valid question. All in good time. At the moment, it is somewhat of a rapid emergence of discovery & I'm fast-tracking, with mostly hand-written notes as I go (over 30 kg of paper I have published some papers in my university seminars. You can read a little between the lines on this forum as to what I'm exploring.

I'll follow up whether it would be acceptable to put these local papers up for review on a website. This stuff has been fun & I honestly wish I had a 48-hour day.

As things stand, I'm in the process of writing up another local paper & am in the process of possibly changing universities. I plan a series of papers during the next few months & am busy collecting suitable material.

You would be most welcome, as would anyone else, to e-mail me directly in the meantime & we can compare notes.

desA

  Reply With Quote

Old   February 6, 2007, 16:36
Default Re: The 'other' half of the continuity equation
  #8
Ahmed
Guest
 
Posts: n/a
Diaw I am not going to write about the merits of Engineers in our modern life, you can see the end result around you by yourself. But the point is when Engineers make statements, these are generally correct statements because they are based on real life experiences. Now when you try to apply their statements it is your responsibility to make sure that you are having the right conditions to apply these statements, this is a very tricky point. In my previous posting I said when the flow Mach number is below 0.3, then the incompressibility assumption is valid, if you check all books written about Aerodynamics you will find this same assumption and a lot of experimental work has validated it. Now if you take this assumption outside the field of aerodynamics, it is your responsibility to make sure you meet the right conditions. Industrial accidents do happen everywhere no matter what the level of technology is (Just remember the Chernobeel reactor and who was conducting the tests at the time) I was a wittness of one of those accidents where life and property loss was enormous (to say the least), that accident was motivated by that same incompressibilty assumption, water, as you can check all technical references written by engineers and scientists, is treated as incompressible fluid, but the accident I am talking about happened because that assumption was taken blindly to its full extent, but neglecting the effect of temperature drop and how the whole physical process will take place no matter how you treat it mathematically. Excluding the compressibilty of matter because the Mach number is below 0.3 is valid only in Aerodynamic applications, surje waves do happen in Hydraulic applications at much lower speeds, and they develop because Water has compressibility Enjoy the reading of what is not written. My computer suffered from one of these accidents, and I had to install a new hard drive, I lost all my email system, I wish you could send me your email address How the preparations for the world cup are going on
  Reply With Quote

Old   February 6, 2007, 18:13
Default Re: The 'other' half of the continuity equation
  #9
Ahmed
Guest
 
Posts: n/a
I would like you to log on the following web site: http://web.mit.edu/fluids/www/Shapiro/ncfmf.html Download the file entitled Pressure fields and acceleration If you are a carefull listener, you will notice how professor Shapiro referes to the water as "nearly incompressible" fluid. Enjoy your time
  Reply With Quote

Old   February 6, 2007, 22:51
Default Re: The 'other' half of the continuity equation
  #10
diaw (Des_Aubery)
Guest
 
Posts: n/a
Hi Ahmed,

Thanks for your comments. I am an Engineer by the way, on a Physics platform, & so I understand precisely where you are coming from. I've also been around the block a few times too many.

I have also had the same approach to understanding incompressibility & that is why I'm working on gaining an understanding of the physical mechanisms at work, in order to bring the mathematical-modeling techniques we rely on, closer to reality. This incompressible / nearly-incompressible issue is something not currently well-enough understood - it has been easier for folks to work at high-speed than at low speed.

Thanks for the Shapiro link - I do have the set on file. I'll work through the link you've suggested - thanks so much for that.

If you draw an analogy to how electromagnetic waves are generated, the physics books mention that these occur when a particle is accelerated - we should, thus, by inference 'expect' to see wave phenomena in fluid flows - we haven't really spent time looking.

The point I'm currently trying to settle is the appropriate closure form of the 'Continuity Equation' to best represent the physics of the situation. I have shown concern for the traditional mathematical approach of a deficient in 'p' closure of form div(v)=0.

Thanks for your input. My e-mail is <momentumwaves@gmail.com> I look forward to your direct e-mail.

desA
  Reply With Quote

Old   February 7, 2007, 08:35
Default Re: The 'other' half of the continuity equation
  #11
desA
Guest
 
Posts: n/a
Hi Ahmed,

I had another look at Shapiro this afternoon. A few things were pretty interesting, & these were not directly covered in his talk. They were the 'unsteady effects' in the following situatons:

1. Hydrogen bubbles coming onto the nose of an obstruction. Observe closely the bubble patterns, from some distance before impact.

2. Hydrogen bubbles moving into constricted flow region (1/2 converging nozzle). Observe the bubble packets as they accelerate towards the rights of the screen.

3. Hydrogen streamline continuous trace through same 1/2 nozzle. Observe bubble pattern to right of screen & tell me if it looks as smooth as the entry pattern. If not, what do you see?

4. The water test nozzle at startup. Two shots of this. Very interesting unsteady flow-patterns.

Very interesting stuff.

desA
  Reply With Quote

Reply

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 write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 16:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58
continuity equation Rafal Main CFD Forum 4 November 29, 2006 10:27
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 08:52.