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

How to enforce the continuity equation?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 13, 2005, 12:29
Default How to enforce the continuity equation?
  #1
dragon
Guest
 
Posts: n/a
1)Q1:

For pure implicit method,the simple method is always employed to decouple pressure and velocity, mass conservation is enforced only at the end of the outer iteration.

My question is :

When the convective term is treated explicitly (so there is no nonlinearity here),and viscous term implicitly,then how to enforce the continuity equation? Is the following right?

a)calculate u* form the previous pressure p0.

b)solve pressure correction p' from u*

c)correct p=p0+p'

d)use p0=p go to 1) again until the convergent solution is reached. Is that right?

2)Q2: Runge-kutta methods are always explicit,how to enforce continuity using RK method?

Any comment or referrence is greatly appreciated.

  Reply With Quote

Old   March 13, 2005, 13:25
Default Re: How to enforce the continuity equation?
  #2
Hrvoje Jasak
Guest
 
Posts: n/a
You forgot two things:

1) After solving the pressure correction equation, correct the face fluxes with the calculated pressure gradient - it is the pressure equation that guaranteen that your fluxes will be conservative

2) You will need some under-relaxation - explicit for p and built into the equation for U (if steady, for p alays the same).

Enjoy,

Hrv
  Reply With Quote

Old   March 13, 2005, 20:59
Default Re: How to enforce the continuity equation?
  #3
dragon
Guest
 
Posts: n/a
thanks dear Hrv But since the convective term is treated explicitly,so the coeff of this term is determined only by the face fluxes of previous time step,while the viscous term doesn't use fluxes at all. So why do we correct the face fluxes ? It seems we don't use it at all.
  Reply With Quote

Old   March 14, 2005, 02:55
Default Re: How to enforce the continuity equation?
  #4
Hrvoje Jasak
Guest
 
Posts: n/a
OK, you are doing Finite Volume, right? Therefore, the convection term is evaluated by multiplying the face flux by the face value. So:

- for the passive scalar theta it would flux F times theta on the face:

F*\theta_f

- for velocity it would be flux times the velocity on the face:

F*\U_f

In both cases the "on the face" bit is calculated using your favourite differencing scheme but the face flux is always used!

Implicit and explicit:

- for implicit, you will express theta on the face as alpha times theta left plus beta times theta right (and maybe fiddle with it a bit more) and then get the matrix coefficients as alpha times the flux and beta times the flux:

\theta_f = \alpha*\theta_P + \beta*\theta_N

a_P = \alpha*F

a_N = \beta*F

(typically \alpha + \beta = 1)

- for explicit, you will evaluate the flux and put the whole thing into the source/residual instead of the matrix, but the idea is the same, i.e. you always need the flux.

Better?

Hrv
  Reply With Quote

Old   March 14, 2005, 03:47
Default Re: How to enforce the continuity equation?
  #5
dragon
Guest
 
Posts: n/a
yes ,I use finite volume. I know we always need to use the flux,what I mean is the following: 1)For explicit convection term, the flux is evaluated from the previous time step(and put into the source), and so there is no need to correct the flux after the P' equation is solved and the velocity need not be corrected either. The P' is only used for pressure correction. 2)For implicit convection term, the flux and velocity are corrected after the P' is solved.

Is that right?

thanks dear hrv

dragon
  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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
mixture mass continuity v.s. volumetric mixture continuity in pressure equation kaifu OpenFOAM 0 June 9, 2011 08:14
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
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 14:53.