Stationary Stokes flow with FVM
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?
Thanks in advance!
you need to do Rhie & Chow interpolation to avoid pressure-velocity decoupling on collocated grids.
Reference to original paper:
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:
phi=fvc::interpolate(U) & mesh.Sf()
+ fvc::snGrad(p) * fvc::interpolate(Vmesh);
I suppose these three above items are what Rhie-Chow interpolation is, but the compiling of it failed. I don't know why and how to correct it. Could you do me a favor?
|All times are GMT -4. The time now is 17:30.|