File Exchange

image thumbnail


version (2.18 KB) by Hassan Khawaja
CFD solution using SIMPLE


Updated 09 May 2019

View License

This work presents a method of solution of fundamental governing equations of computational fluid dynamics (CFD) using Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) in MATLAB®. The fundamental governing equations of fluid mechanics are based on three laws of conservation, referred to the law of conservation of mass, the law of conservation of momentum and law of conservation of energy. Continuity equation represents the law of conservation of mass, Navier-Stokes equations represents the law of conservation of momentum, and energy equation represents the law of conservation of energy. In SIMPLE, the continuity and Navier-Stokes equations are required to be discretized and solved in a semi-implicit way. This article presents discretization and method of solution applied to the flow around a 2-D square body. Code is written in MATLAB®. The results show the pressure and velocity fields of the converged solution.

Cite As

Hassan Khawaja (2019). SIMPLE (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (2)

Line 155 throws an error if Nx and Ny are not equal. Because dp1 is [Nx Nx] dimension, but p has dimension [Ny Nx]


Thanks for the comments. The code has been corrected. It now computes flow around a square. Nx and Ny can be varied.

Revised abstract

MATLAB Release Compatibility
Created with R2016b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.

Learn About Live Editor