Creating IB-Solver for turbulent channel flow
 

I expanded the icoIbFoam solver in foam extend 3.1 by the pressure-correction used in the solver channelFoam for a constant mass flow rate so i can use the new solver as DNS-solver for turbulent channel flow with IB. The problem is that the pressure- correction seems to be performed in the whole channel including the IB. This leads to a wrong mass flow. Therfore I was wondering if there is an easy way to change the code so that the pressure-correction is only performed in the cells outside the IB?

Link to used solver:
https://www.dropbox.com/sh/wfgf194ul...u3_FkDEAa?dl=0

In IB method pressure correction is usually solved in all domain including the marked regions.

Note that this is not hard and fast rule, I in my own code did solve only in flow regions due to problems in accuracy (Not openfoam code). It needs lots of change in code to do because matrix will now be smaller and different.

If you are facing accuracy issues mesh refinement is only friend you have.