CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Stationary Stokes flow with FVM (https://www.cfd-online.com/Forums/openfoam/64231-stationary-stokes-flow-fvm.html)

 christoph_heinrich May 4, 2009 11:00

Stationary Stokes flow with FVM

Hi everybody,

at the moment I am trying to implement a stationary Stokes solver in Matlab for a simple geometry. I am using the FVM with midpoint rule and CDS. In the pressure field I get oscillations. How can I circumvent this problem on a collocated grid? Any idea?

Kind regards
Christoph

 henrik May 5, 2009 12:30

Dear Christoph,

you need to do Rhie & Chow interpolation to avoid pressure-velocity decoupling on collocated grids.

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/rhiechow.pdf

Reference to original paper:
http://powerlab.fsb.hr/ped/kturbo/Op...jeJasakPhD.pdf

Henrik

 victorconan August 27, 2012 08:50

Quote:
 Originally Posted by henrik (Post 215216) Dear Christoph, you need to do Rhie & Chow interpolation to avoid pressure-velocity decoupling on collocated grids. http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/rhiechow.pdf Reference to original paper: http://powerlab.fsb.hr/ped/kturbo/Op...jeJasakPhD.pdf Henrik
Hi, Henrik
I suppose the Rhie-Chow interpolation mentioned by Fabian Peng Karrholm does not exist in Openfoam at all. According to Rhie-Chow interpolation, we should have interpolation for U and p on surface respectively. I have felt puzzled for a long time about how to make it in openfoam. I tried this:
volScalarField Vmesh=mesh.V()
phi=fvc::interpolate(U) & mesh.Sf()
- pEqn.flux()