CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Stationary Stokes flow with FVM

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 4, 2009, 11:00
Default Stationary Stokes flow with FVM
  #1
New Member
 
Christoph Heinrich
Join Date: Mar 2009
Posts: 12
Rep Power: 7
christoph_heinrich is on a distinguished road
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?

Thanks in advance!

Kind regards
Christoph
christoph_heinrich is offline   Reply With Quote

Old   May 5, 2009, 12:30
Default
  #2
Senior Member
 
Henrik Rusche
Join Date: Mar 2009
Location: Braunschweig, Niedersachsen, Germany
Posts: 272
Rep Power: 8
henrik is on a distinguished road
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
henrik is offline   Reply With Quote

Old   August 27, 2012, 08:50
Default
  #3
New Member
 
Lun Yu
Join Date: Aug 2012
Posts: 8
Rep Power: 3
victorconan is on a distinguished road
Quote:
Originally Posted by henrik View Post
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()
+ 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?
victorconan is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Stokes flow Rajneesh FLUENT 2 June 18, 2013 07:40
Stokes flow Shuo Main CFD Forum 0 May 4, 2008 21:51
Flow laminar and stationary of water in a pipe manuel OpenFOAM Running, Solving & CFD 6 March 24, 2007 19:23
Stokes flow on cylinder Khurram Main CFD Forum 0 August 22, 2004 19:10
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 12:19


All times are GMT -4. The time now is 23:23.