Creating IB-Solver for turbulent channel flow

urion · May 3, 2016, 10:42

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

arjun · May 12, 2016, 10:21

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.

May 3, 2016, 10:42	Creating IB-Solver for turbulent channel flow	#1
urion New Member K.H Join Date: Feb 2016 Posts: 15 Rep Power: 10	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

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
compressible flow calculation error using rhoSimpleFoam solver	student4326	OpenFOAM Running, Solving & CFD	7	November 2, 2015 12:34
Using pimpleFoam to model pulsatile flow in a micro channel	isaacp	OpenFOAM Running, Solving & CFD	0	January 30, 2012 01:03
ATTENTION! Reliability problems in CFX 5.7	Joseph	CFX	14	April 20, 2010 16:45
Wrong Pressure countors for channel flow driven by mass flow rate or press gradient	ash-khan	FLUENT	0	June 21, 2009 11:02
Open Channel Flow using InterFoam type solver	sxhdhi	OpenFOAM Running, Solving & CFD	3	May 5, 2009 22:58

May 12, 2016, 10:21		#2
arjun Senior Member Arjun Join Date: Mar 2009 Location: Nurenberg, Germany Posts: 1,290 Rep Power: 34	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.