Rhie & Chow interpolation
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 colocated 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 
Re: Rhie & Chow interpolation
i am also a fresh in CFD.And does your Rhie and Chow mean some similarity as ROE?

Re: Rhie & Chow interpolation
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 staggeredmesh system. The volume fluxes are not solution variables, but rather are determined through interpolation of the cellcentered u_i values plus a projection operation that guarantees exact conservation of mass. Use of the massconserving volume fluxes results in a pressure equation identical to that in the staggeredmesh system and thus also leads to fullycoupled 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 90124), 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 (i1,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 noncorrected 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. Sincerely, Frederic Felten. 
Re: Rhie & Chow interpolation
Dear Frderic
Thank you for the explanation. It goes a long way in improving my understanding of the matter! Regards Barry 
All times are GMT 4. The time now is 13:31. 