Jump to content

SIMPLEC algorithm

From Wikipedia, the free encyclopedia
This is an old revision of this page, as edited by Kshtzsrswt (talk | contribs) at 19:21, 11 November 2014. The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

The SIMPLEC (Semi-Implicit Method for Pressure Linked Equations-Consistent) algorithm; a modified form of SIMPLE algorithm; is a commonly used numerical procedure in Computational Fluid Dynamics to solve the very famous Navier-Stokes equations.
This algorithm was developed by Van Doormal and Raithby in 1984. The algorithm follows the same steps like the SIMPLE algorithm with a little variation that the momentum equations are manipulated which allows SIMPLEC velocity correction equations to omit terms that are less significant than those omitted in SIMPLE. Basically, it tries to prevent the effects of dropping velocity neighbor correction terms.
[1]

Algorithm

The steps involved are same as the SIMPLE algorithm and the algorithm is iterative in nature. Steps are as follows:

  1. Specify the boundary conditions and guess the initial values.
  2. Determine the velocity and pressure gradients.
  3. Calculate the pseudo velocities.
  4. Solve for the pressure equation and get the p.
  5. Set p*=p.
  6. Using p* solve the discretized momentum equation and get u* and v*.
  7. Solve the pressure correction equation.
  8. Get the pressure correction term and evaluate the corrected velocities and get p, u, v, Φ*.
  9. Solve all other discretized transport equations.
  10. If Φ shows convergence, then STOP and if not, then set p*=p, u*=u, v*=v, Φ*=Φ and start the iteration again.
    [2][3]
Pseudo velocity equations in x and y dir.
Pressure equation
Discretized momentum equations
Pressure correction equation
Velocity Correction equations in X and Y dir.
Transport equations

Peculiar features

  • The discretized pressure correction equation is same as in the SIMPLE algorithm, except for the d terms which are used in momentum equations.
  • p=p*+p' which tells that the under relaxing factor is not there in SIMPLEC as it was in SIMPLE.
  • SIMPLEC algorithm is seen to converge 1.2-1.3 times faster than the SIMPLE algorithm
  • It doesn't solve extra equations like SIMPLER algorithm.
  • The cost per iteration is same as in the case of SIMPLE.
  • It also destroys a good velocity field; alike SIMPLE; unless there is also a good pressure guess.[4]

See also

References

  1. ^ "Variants of SIMPLE algorithm" (PDF). engineering.purdue.edu. Retrieved 11 November 2014.
  2. ^ Versteeg, H.K.; Malalasekera, W. An introduction to Computational Fluid Dynamics- The finite volume method (1st edition, 1995 ed.). Longman Group Ltd. p. 149-151.
  3. ^ Patankar, S. V. (1980). Numerical Heat Transfer and Fluid Flow. Taylor & Francis. ISBN 978-0-89116-522-4.
  4. ^ "SIMPLE solver for driven cavity problem" (PDF). engineering.purdue.edu.

Kshtzsrswt (talk)Kshitiz Saraswat, Kushal Mundle