Physics:SIMPLE algorithm

From HandWiki
Revision as of 23:58, 5 March 2023 by Carolyn (talk | contribs) (linkage)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

In computational fluid dynamics (CFD), the SIMPLE algorithm is a widely used numerical procedure to solve the Navier–Stokes equations. SIMPLE is an acronym for Semi-Implicit Method for Pressure Linked Equations. The SIMPLE algorithm was developed by Prof. Brian Spalding and his student Suhas Patankar at Imperial College, London in the early 1970s. Since then it has been extensively used by many researchers to solve different kinds of fluid flow and heat transfer problems.[1]

Many popular books on computational fluid dynamics discuss the SIMPLE algorithm in detail.[2][3] A modified variant is the SIMPLER algorithm (SIMPLE Revised), that was introduced by Patankar in 1979.[4]

Algorithm

The algorithm is iterative. The basic steps in the solution update are as follows:

  1. Set the boundary conditions.
  2. Compute the gradients of velocity and pressure.
  3. Solve the discretized momentum equation to compute the intermediate velocity field.
  4. Compute the uncorrected mass fluxes at faces.
  5. Solve the pressure correction equation to produce cell values of the pressure correction.
  6. Update the pressure field: [math]\displaystyle{ p^{k + 1} = p^k + \text{urf} \cdot p^{'} }[/math] where urf is the under-relaxation factor for pressure.
  7. Update the boundary pressure corrections [math]\displaystyle{ p_b^{'} }[/math].
  8. Correct the face mass fluxes: [math]\displaystyle{ \dot m_f^{k + 1} = \dot m_f^{*} + \dot m_f^{'} }[/math]
  9. Correct the cell velocities: [math]\displaystyle{ \vec v^{k + 1} = \vec v^{*} - \frac{{\text{Vol} \ \nabla p^{'} }}{{\vec a_P^v }} }[/math] ; where [math]\displaystyle{ {\nabla p^{'} } }[/math] is the gradient of the pressure corrections, [math]\displaystyle{ {\vec a_P^v } }[/math] is the vector of central coefficients for the discretized linear system representing the velocity equation and Vol is the cell volume.
  10. Update density due to pressure changes.

See also

References

  1. Mangani, L.; Bianchini, C. (2007). "Heat transfer applications in turbomachinery". Proceedings of the OpenFOAM International Conference 2007. https://flore.unifi.it/retrieve/handle/2158/418277/15222/OFIC-07.pdf. Retrieved 2016-03-16. 
  2. Patankar, S. V. (1980). Numerical Heat Transfer and Fluid Flow. Taylor & Francis. ISBN 978-0-89116-522-4. 
  3. Ferziger, J. H.; Peric, M. (2001). Computational Methods for Fluid Dynamics. Springer-Verlag. ISBN 978-3-540-42074-3. 
  4. Tannehill, J. C.; Anderson, D. A.; Pletcher, R. H. (1997). Computational Fluid Mechanics and Heat Transfer. Taylor & Francis. ISBN 9781560320463. https://archive.org/details/computationalflu0000tann.