CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Velocity Under-relaxation in SIMPLE type methods

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

LinkBack Thread Tools Display Modes
Old   June 27, 2005, 15:21
Default Velocity Under-relaxation in SIMPLE type methods
Matt U.
Posts: n/a
A (somewhat general) description of the SIMPLE and the SIMPLEC algorithm for pressure correction for incompressible flows for implicit time integration is given by:

0. Use values from previous time step to start outer iteration.

1. Form coefficients and solve momentum equations for intermediate velocity or u*.

2. Solve continuity equation (as the pressure correction equation) for p'.

3. Correct the pressure using p'.

4. Correct the velocity using u'.

5. Use corrected velocities and pressures as improved estimates and repeat from Step 1 until values of p' and u' corrections are "small enough".

(note this description has been described by numerous authors; the current description was loosely taken from page 168 of Ferziger's Computational Methods for Fluid Dynamics).

NOW MY QUESTION: Where do you apply velocity under-relaxation? Do you apply it in Step 1, where the solution of the momentum equations is done for u*, OR in Step 4, where the velocities are updated before the next outer iteration OR does it even matter?

Note that if you under-relax velocity in Step 1 it can be done by incorporating the under-relaxation factor in the matrix coefficients in the discrete momentum equation; otherwise the velocity under- relaxation could take place after the velocity correction in Step 4, where the new velocity would be a combination of the corrected velocity and the velocity from the previous outer iteration.


Matt U.

  Reply With Quote

Old   June 27, 2005, 23:08
Default Re: Velocity Under-relaxation in SIMPLE type metho
Posts: n/a
As far as I know, the under-relaxation in SIMPLE is integrated into the matrix calculation, in other words, it is embodied in the coefficiants AIP,AIM,CON, and so on.
  Reply With Quote

Old   June 28, 2005, 03:31
Default Re: Velocity Under-relaxation in SIMPLE type metho
Posts: n/a
If you apply it at 4 the corrected velocity field will not satisfy continuity (because you have reduced the correction). The pressure correction at 3 is usually reduced relative to the velocity correction.

The relaxation is added at 1 and, usually, included in the coefficients of the pressure correction equation.

The word usually is occurring because there are lots of variations in the details of pressure correction schemes.

  Reply With Quote

Old   June 28, 2005, 04:44
Default Re: Velocity Under-relaxation in SIMPLE type metho
Halim Choi
Posts: n/a

(1) The momentum equation can be written as follows with the underrelaxation factor written explicitly at the east face, for example;

(Ap/alpa)e Ue=(sigma(AnbUnb))e+be+(1-alpa)(Ap/alpa)e Ue,n-1---(1)

where n-1 is the previous iteration level and alpa is the underrelaxation factor.

Redefine (Ap)e=(Ap/alpa)e -----(2)

then the above equation can be written as follows;

(Ap)e Ue=(sigma(AnbUnb))e+be+(1-alpa)(Ap)e Ue,n-1 ---(3)

Solve above equation to get Ue*.

(2) In stage 4, we need velocity correction equations. The above equation can be written as follows,

Ue=(Hu)e+(Du)e(Pp-PE)+other terms,

where (Du)e=(dy/Ap)e. (Hu)e=(sigma(AnbUnb)+ScdV)/(Ap)e,

Note that the underrelaxationfactor is embodied in (Du)e by equation (2) and ScdV is the source term (be) without the pressure gradient term.

The velocity correction equation in SIMPLE can be derived as follows,

Ue=Ue* + (Du)e(P'p-P'E).

Note that Ue* is the velocity obtained by Eq.(3), not the velocity at the previous interation level.

I hope this helps,

Halim Choi
  Reply With Quote

Old   June 29, 2005, 13:33
Default Re: Velocity Under-relaxation in SIMPLE type metho
Posts: n/a
yes, the under/over-relaxation is applied during coefficient generation, before solver
  Reply With Quote

Old   June 30, 2005, 08:55
Default Re: Velocity Under-relaxation in SIMPLE type metho
Matt U.
Posts: n/a
Good point about the velocity field and the continuity equation. Thanks.
  Reply With Quote

Old   July 4, 2005, 05:29
Default Any difference for relaxation in couple solver?
Posts: n/a
Do you know if there's any difference to apply under-relaxation in couple solver (velocity & pressure are solved at the same time)? Someone told me before that relaxation has little effect on coupler solver. I don't understand why. Is this really true?

Thanks for your advice.
  Reply With Quote


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
compilation problem with "fvPatch::lookupPatchField" Ya_Squall2010 OpenFOAM Programming & Development 7 October 27, 2014 10:24
Instable natural convection case Peter88 OpenFOAM 5 August 18, 2011 01:23
turbulent jet simulation antonio_ing OpenFOAM Running, Solving & CFD 5 September 16, 2010 02:31
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 04:05
RasInterFoam STRANGE RESULTS AT BOUNDARY kumar2 OpenFOAM Running, Solving & CFD 8 March 24, 2008 18:38

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