# Laminar flow between parallel plates with SIMPLE algorithm

 Register Blogs Members List Search Today's Posts Mark Forums Read

 November 18, 2016, 15:21 Laminar flow between parallel plates with SIMPLE algorithm #1 New Member   Alper Join Date: Jul 2013 Posts: 8 Rep Power: 11 Hello Folks, I'm trying to code in FORTRAN incompressible laminar flow between two parallel plates. I put uniform inlet velocity. I use finite volume method. I'm using 1-D convection-diffusion equation and representing it by the power law. I'm linking velocity and pressure fields using the SIMPLE algorithm. Dimensions of the channel are Length x Height is 2x0.02. Reynolds number is 150. So my question is, I need to divide this domain to 160x40 control volumes but when I do it, solution blows up. When I divide domain to 40x160 it runs without any problems. So I wonder if anyone has any comments why Length (2m) cannot be divided into 160 control volumes.

 July 19, 2017, 02:02 #2 New Member   Alper Join Date: Jul 2013 Posts: 8 Rep Power: 11 It's been a long time I forgot about this thread. Just in case if anyone will have the similar problem I solved the problem. I realized that in the first few iterations solution was not propagating to the end of the channel when I was dividing length to the 160 meshes. Solution is to add global mass conservation to the code which basically forces mass to be conserved for the first few iterations until solution reaches to the end of the channel. Global mass conservation added to the code by checking if the ratio of total outlet mass to the total inlet mass is equal to 1 after each iteration. If not, it is forced to be. Sent from my iPhone using CFD Online Forum mobile app

 Tags convection-diffusion, laminar flow, parallel plates, simple