SIMPLE algorithm for variable density flow
 

Hello everybody,

I would like to apply the SIMPLE algorithm to a laminar flow in a straight channel. Along the channel there is gain/loss of mass set by boundary conditions. Usually the SIMPLE algorithm is presented for constant density (Patankar, Ferziger Peric books and different articles).

To take into account the variable density, shall I just include the density value in the appropriate node, both in the momentum discreitzation and then in the continuity equation?

thanks in advance

Mauri

Look for papers by the group of Wei Shyy. They have been working on extending SIMPLE algorithm over the entire mach number regime 0-1.