SIMPLE algorithm: Difference between revisions

Revision as of 18:45, 30 August 2013

Latest revision as of 21:20, 7 June 2024

(26 intermediate revisions by 20 users not shown)

Line 1:

{{Short description|Computational fluid dynamics algorithm}}

In [[computational fluid dynamics]] (CFD), the '''SIMPLE algorithm''' is a widely used [[numerical algorithm|numerical procedure]] to solve the [[Navier-StokesNavier–Stokes equation]]s. ''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|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.<ref>{{cite webconference |urllast1=https://engineeringMangani |first1=L.purdue.edu/ME608/webpage/project-reports/SIMPLE-DrivenCavity.pdf |titlelast2=SIMPLEBianchini solver|first2=C. for|conference=[[Proceedings drivenof cavitythe flowOpenFOAM problemInternational |format=PDFConference |date=2007]] |accessdateyear=2011-08-21}}</ref><ref>{{cite2007 web|url=httphttps://wwwflore.opensourcecfdunifi.comit/conference2008retrieve/2007handle/media2158/proceedings418277/15222/OFIC-07_ManganiLuca07.pdf?phpMyAdmin=f7f6d5b5be922e9db00b8b978b532798 |title=Heat transfer applications in turbomachinery |access-date= |accessdate=20112016-0803-2116}}</ref>

Many popular books on computational fluid dynamics discuss the SIMPLE algorithm in detail.<ref>{{cite book |last=Patankar |first=S. V. | authorlinkauthor-link = Suhas Patankar |title=Numerical Heat Transfer and Fluid Flow |publisher=[[Taylor & Francis]] |year=1980 |isbn=978-0-89116-522-4}}</ref><ref>{{cite book |last=Ferziger |first=J. H. | authorlinkauthor-link = J. H. Ferziger |coauthors author2= Peric, M. |title=Computational Methods for Fluid Dynamics|publisher=[[ Springer-Verlag]] |year=2001 |isbn= 978-3-540-42074-3}}</ref> This algorithm forms the basis of Commercial CFD packages.<ref>{{cite web|url=http://www.salihnet.freeservers.com/engineering/cfd/cfd_history.html |title=History of CFD |publisher=Salihnet.freeservers.com |date= |accessdate=2011-08-21}}</ref>

A modified variant is the ''SIMPLER'' algorithm (SIMPLE Revised), that was introduced by Patankar in 1979.<ref>{{cite book |last=Tannehill|first=J. C.| authorlinkauthor2 = |coauthorsAnderson, = [[DaleD. A. Anderson|Anderson,author2-link D.= Dale A.]]; Anderson |author3=Pletcher, R. H. |title=Computational Fluid Mechanics and Heat Transfer |url=https://archive.org/details/computationalflu0000tann|url-access=registration|publisher=[[Taylor & Francis]] |year=1997 |isbn=9781560320463 }}</ref>

A modified variant is the ''SIMPLER'' algorithm (SIMPLE Revised), that was introduced by Patankar in 1979.<ref>{{cite book |last=Tannehill|first=J. C.| authorlink = |coauthors = [[Dale A. Anderson|Anderson, D. A.]]; Pletcher, R. H. |title=Computational Fluid Mechanics and Heat Transfer |publisher=[[Taylor & Francis]] |year=1997 |isbn=}}</ref>

== Algorithm ==

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

# Set the boundary conditions.

# Compute the gradients of velocity and pressure.

# Solve the discretized momentum equation to compute the intermediate velocity field.

# Compute the uncorrected mass fluxes at faces.

# Solve the pressure correction equation to produce cell values of the pressure correction.

# Update the pressure field: <math> p^{k + 1} = p^k + \text{urf} \bulletcdot p^{'} </math> where urf is the under-relaxation factor for pressure.

# Update the boundary pressure corrections <math> p_b^{'} </math>.

# Correct the face mass fluxes: <math>\dot m_f^{k + 1} = \dot m_f^{*} + \dot m_f^{'} </math>

# Correct the cell velocities: <math> \vec v^{k + 1} = \vec v^{*} - \frac{{\text{Vol} \ \nabla p^{'} }}{{\vec a_P^v }} </math> ; where <math> {\nabla p^{'} } </math> is the gradient of the pressure corrections, <math> {\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.

# Update density due to pressure changes.

== See also ==

* [[PISO algorithm]]

* [[SIMPLEC algorithm]]

==References==

{{reflistReflist}}

[[Category:Computational fluid dynamics]]