|
[Sponsors] |
Rhie Chow Interpolation Understanding Question |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 27, 2021, 20:52 |
Rhie Chow Interpolation Understanding Question
|
#1 |
Member
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13 |
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)? |
|
February 28, 2021, 05:00 |
|
#2 |
Senior Member
Gerry Kan
Join Date: May 2016
Posts: 361
Rep Power: 11 |
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. Last edited by Gerry Kan; February 28, 2021 at 08:40. |
|
February 28, 2021, 08:55 |
|
#3 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,840
Rep Power: 73 |
||
February 28, 2021, 13:46 |
|
#4 | |
Member
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13 |
Quote:
|
||
February 28, 2021, 13:47 |
|
#5 | |
Member
Raphael
Join Date: Nov 2012
Posts: 68
Rep Power: 13 |
Quote:
|
||
February 28, 2021, 17:09 |
|
#6 |
Senior Member
Gerry Kan
Join Date: May 2016
Posts: 361
Rep Power: 11 |
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. |
|
February 28, 2021, 17:13 |
|
#7 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,840
Rep Power: 73 |
||
March 1, 2021, 01:07 |
|
#8 | ||
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34 |
Quote:
Quote:
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) 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. |
|||
Tags |
deferred correction, rhie chow |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |