MATLAB code to solve the Channel Flow Problem

Matlab code to solve the 2D steady state incompressible NavierStokes equations using SIMPLE algorithm for a Channel

This is a MATLAB code that uses Finite Volume Method to discretize the channel flow domain to solve the continuity and the X,Y momentum equations using the Semi-Implicit Method for Pressure Linked Equation (SIMPLE). U,V velocities are declared and solved along the staggered mesh while the pressure uses the normal mesh generated after discretization. The upwind differencing scheme of interpolation is used for obtaining the face velocities at the staggered mesh faces. In addition the same code has certain lines commented out that can be un-commented to introduce a block or an obstacle in the centre of the pipe and obtain the corresponding contours. Pressure correction equations are overrelaxed. X,Y momentum under relaxed. After velocity update mass imbalance monitor is displayed to see it die out and satisfy mass conservation. Post-processing is done using contourf and quiver functions