CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Rhie & Chow interpolation

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

Like Tree2Likes
  • 2 Post By frederic felten

LinkBack Thread Tools Search this Thread Display Modes
Old   January 17, 2003, 06:41
Default Rhie & Chow interpolation
Posts: n/a
My CFD friends...

Can someone explain to me in plain and simple English what the idea behind Rhie and Chow interpolation is? I know that its purpose is to couple pressure and velocity on co-located meshes, but how does it work?

How do I calculate the pressure gradient (face pressure) before I solve the momentum equation? Do we use linear interpolation for the pressure or some other scheme? And what about the pressure correction equation? Where does it fit in?

Thank you. Looking forward to your replies! Barry
  Reply With Quote

Old   January 17, 2003, 06:57
Default Re: Rhie & Chow interpolation
Posts: n/a
i am also a fresh in CFD.And does your Rhie and Chow mean some similarity as ROE?
  Reply With Quote

Old   January 17, 2003, 12:18
Default Re: Rhie & Chow interpolation
frederic felten
Posts: n/a
Hi there,

Here is my understanding of the Rhie & Chow interpolation procedure and some of the issues associated with it.

The Cartesian velocity components u_i (or u, v, w) are stored with the pressure, p, at the cell center. In addition, contravariant volume fluxes, F_i ( i= 1, 2, 3), are defined at the cell face in a manner analogous to the staggered-mesh system. The volume fluxes are not solution variables, but rather are determined through interpolation of the cell-centered u_i values plus a projection operation that guarantees exact conservation of mass. Use of the mass-conserving volume fluxes results in a pressure equation identical to that in the staggered-mesh system and thus also leads to fully-coupled velocity and pressure fields. While the pressure field determined in this manner leads to mass conserving volume fluxes, it leaves the primary solution variables, u_i, only approximately divergence free. As pointed out by Morinishi et al.(J Comput. Phys., 143, 1998, pp 90-124), this defect leads to one source of kinetic energy conservation error.

In my case, the pressure gradient on the face, is obtain just by differentiating the pressure values on both side of the face. In 1D, you then only need 3 points (i-1,i,i+1) to write the laplace operator.

Now as far as the correction step for the cartesian cell center velocities is concern, this step will allow you to reduce the divergence level of the primary variables u_i. I did some testing involving this correction step and whether I did it or not, the code (turbulent channel flow) was still running fine, providing good results in both cases. Still, when the correction step was performed, the level of divergence came back half the non-corrected one. My advice to you is: do perform this correction step so that you'll avoid any kind of funky behavior that could eventually lead to instability.

Hope this helps.


Frederic Felten.

utkunun and galleon like this.
  Reply With Quote

Old   January 20, 2003, 09:11
Default Re: Rhie & Chow interpolation
Posts: n/a
Dear Frderic

Thank you for the explanation. It goes a long way in improving my understanding of the matter!

Regards Barry
  Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Collocated grid (Rhie and Chow) C.Robin Main CFD Forum 5 February 29, 2008 03:19
Rhie and Chow James Houghton Main CFD Forum 0 October 9, 2005 07:02
Rhie & Chow method/flame propagation? yang Main CFD Forum 0 May 19, 2005 08:00
Rhie - Chow Interpolation ? Javier Goicochea Main CFD Forum 0 November 27, 2003 18:01
Rhie - Chow speed up factors Goicox Main CFD Forum 0 November 22, 2003 00:00

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