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

Rhie Chow Interpolation Understanding Question

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

Like Tree4Likes
  • 1 Post By Gerry Kan
  • 1 Post By FMDenaro
  • 1 Post By FMDenaro
  • 1 Post By arjun

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 27, 2021, 20:52
Default Rhie Chow Interpolation Understanding Question
  #1
Member
 
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13
arkie87 is on a distinguished road
I've been studying Rhie Chow interpolation for fixing of checkerboard pressure profiles in collocated CFD simulations and have a question.

The way I understand it, Rhie Chow is necessary because if central differencing about the cells is used to obtain the pressure force in the momentum equations (with collocated variables), then in the continuity/pressure correction equation, a large stencil will result with 9 points (P, W, E, N, S, WW, EE, SS, NN). Not only is this hard to solve, but it can result in a checkboard pressure profile.

Conversely, if central differencing about the faces is used to make the pressure correction equation, a 5 point stencil (P, W, E, S, N) is obtained, which is easier to solve and is oscillation free, but is inconsistent with the way pressure is handled in the momentum equation and can lead to instability.

To fix this, the 9 point stencil resulting from central differencing of the cells is treated explicitly, and the 5 point stencil is added and subtracted from the continuity equation, with one term treated implicitly and one term treated explicitly. Thus, the stencil for solving the pressure correction equation has only 5 points, and an explicit term involving the difference between the 5 point and 9 point stencils is obtained.

This is essentially what Rhie Chow interpolation is, and in fact, it is very similar to a standard "deferred correction" approach, where a more favorable stencil is added and subtracted, with half being treated implicitly and half being treated explicitly, such that an convergence, the implicit and explicit terms cancel. It is typically done to arrive at the solution using a higher order discretization scheme; however, it does not remove the oscillations in the solution, but rather, just provides an convenient way to arrive at the solution.

TL;DR: thus, my question is: how does Rhie Chow interpolation remove the oscillations in the pressure field, given that it is essentially just a deferred correction of the pressure correction equation (and deferred correction approaches do not remove oscillations from high order discretization schemes)?
arkie87 is offline   Reply With Quote

Old   February 28, 2021, 05:00
Arrow
  #2
Senior Member
 
Gerry Kan's Avatar
 
Gerry Kan
Join Date: May 2016
Posts: 361
Rep Power: 11
Gerry Kan is on a distinguished road
Dear Rapheal:

The Rhie and Chow interpolation prevents pressure checkboard in collocated grid solvers by allowing face pressure to be reconstructed using face and cell values. The interpolation is, if memory serves, fourth order accurate, so in theory it will not suppress numerical oscillations beyond this.

The momentum equation is only a function of pressure gradient. So in a collocated grid elliptic solver, the gradient will be reconstructed from the pressure on the cell faces using, say, central difference. Since the grid is collocated, the face values are constructed by taking the average between neighboring cells.

Going through this derivation, you will quickly realize that the pressure gradient depends only on the pressure of the neighboring cells, and not the cell point in question. That means there will be two possible pressure fields, for the alternating 'row' and 'column' sets, that will satisfy the same pressure gradient. This is the origin of the pressure checkerboard problem.

What the Rhie and Chow interpolation does, is to use the pressure correction equation (effectively mass conservation) to reconstruct the face pressures as functions of the cell point and its immediate neighbors, so that only one pressure field will be present instead of two alternating ones, thereby eliminating the checkerboard.

Gerry.
aero_head likes this.

Last edited by Gerry Kan; February 28, 2021 at 08:40.
Gerry Kan is offline   Reply With Quote

Old   February 28, 2021, 08:55
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,840
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Question about Rhie-Chow Interpolation
aero_head likes this.
FMDenaro is online now   Reply With Quote

Old   February 28, 2021, 13:46
Default
  #4
Member
 
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13
arkie87 is on a distinguished road
Quote:
Originally Posted by Gerry Kan View Post
Dear Rapheal:

The Rhie and Chow interpolation prevents pressure checkboard in collocated grid solvers by allowing face pressure to be reconstructed using face and cell values. The interpolation is, if memory serves, fourth order accurate, so in theory it will not suppress numerical oscillations beyond this.

The momentum equation is only a function of pressure gradient. So in a collocated grid elliptic solver, the gradient will be reconstructed from the pressure on the cell faces using, say, central difference. Since the grid is collocated, the face values are constructed by taking the average between neighboring cells.

Going through this derivation, you will quickly realize that the pressure gradient depends only on the pressure of the neighboring cells, and not the cell point in question. That means there will be two possible pressure fields, for the alternating 'row' and 'column' sets, that will satisfy the same pressure gradient. This is the origin of the pressure checkerboard problem.

What the Rhie and Chow interpolation does, is to use the pressure correction equation (effectively mass conservation) to reconstruct the face pressures as functions of the cell point and its immediate neighbors, so that only one pressure field will be present instead of two alternating ones, thereby eliminating the checkerboard.

Gerry.
That is all well and good, but does not directly address my claims or question.
arkie87 is offline   Reply With Quote

Old   February 28, 2021, 13:47
Default
  #5
Member
 
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13
arkie87 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Are you saying the answer to this question is in that thread? If so, which comment?
arkie87 is offline   Reply With Quote

Old   February 28, 2021, 17:09
Default
  #6
Senior Member
 
Gerry Kan's Avatar
 
Gerry Kan
Join Date: May 2016
Posts: 361
Rep Power: 11
Gerry Kan is on a distinguished road
Quote:
Originally Posted by arkie87 View Post
That is all well and good, but does not directly address my claims or question.
You asked how the Rhie and Chow interpolation solves the pressure checkerboard. I explained to you how it works, i.e., making sure that all immediate terms are present, albeit without any math. Then you came around and told me that my explanation is no good. I don't know what other claims and questions I have missed.

Last edited by Gerry Kan; March 1, 2021 at 02:45.
Gerry Kan is offline   Reply With Quote

Old   February 28, 2021, 17:13
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,840
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by arkie87 View Post
Are you saying the answer to this question is in that thread? If so, which comment?



If you read carefully the comments and the references you will find the answers
Gerry Kan likes this.
FMDenaro is online now   Reply With Quote

Old   March 1, 2021, 01:07
Default
  #8
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by arkie87 View Post
That is all well and good, but does not directly address my claims or question.
This is your claim or question:


Quote:
Originally Posted by arkie87 View Post
thus, my question is: how does Rhie Chow interpolation remove the oscillations in the pressure field, given that it is essentially just a deferred correction of the pressure correction equation (and deferred correction approaches do not remove oscillations from high order discretization schemes)?

The basis of your question is wrong because Rhie and Chow correction is not deferred correction approach.

The deferred correction is this:

Code:
phi_face  =   phi_face_implicit +  ( phi_face_HigherOrder    - phi_face_lowerOrder)
Here phi_face_implicit is basically phi_face_lowerOrder treated implicitly.

This is not what Rhie and Chow is. Rhie and Chow is correction now and is defect correction scheme. If there is a defect then it is corrected now (not deferred correction).

The link posted above explains in detail.
Gerry Kan likes this.
arjun is offline   Reply With Quote

Reply

Tags
deferred correction, rhie chow

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
momentum interpolation for collocated grid Hadian Main CFD Forum 4 December 25, 2009 07:25
Interpolation question sdk Main CFD Forum 0 December 8, 2009 13:54
Question About solvers and interpolation schemes titio OpenFOAM Running, Solving & CFD 0 July 9, 2009 11:41
Pressure Interpolation VOF nico765 OpenFOAM Running, Solving & CFD 2 January 3, 2009 14:47
A question in understanding icoLagrangeFoam by Gschaider xiao OpenFOAM Running, Solving & CFD 7 October 27, 2008 02:17


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