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

momentum interpolation for collocated grid

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

LinkBack Thread Tools Search this Thread Display Modes
Old   September 14, 2004, 09:11
Default momentum interpolation for collocated grid
Posts: n/a
Dear friends

I am developing a 2D implicit unsteady code using collocated grid. As you know, to avoid the Checkerboard pressure prediction, the momentum interpolation is needed. since I use the under-relaxation factor in solving momentum equation the original Rhie&Chow method can not be used directly. To have the results independent of under-relaxation and time increment a modified Rhie&Chow interpolation should be used. One of them has been explained in the following paper: Lien ,F.S.; Leschziner, M.A. "A General Non-Orthogonal Collocated Finite Volume Algorithm for Turbulent flow at all Speeds…..", Comput. Methods Appl. Mech. Engrg, 114(1994) 123-148 Now my question is about implementation of such a method. 1- In this type of momentum interpolation the Ap at cell face is needed. should I interpolate the Ap by linear interpolation or 1/Ap should be interpolated? 2- In the momentum interpolation equation there are the cell face velocity at previous iteration and previous time step. So we should store 2 velocity component at each face. So in this case the required memory will be more than what is needed for staggered grid. Is this true? Or there is a trick not to store the velocities at cell faces? 3- Should I use the cell face velocity resulted from momentum interpolation only in pressure correction or it should be used in momentum equation too?

Thanks in advance for your replies.



  Reply With Quote

Old   September 19, 2004, 20:05
Default Re: momentum interpolation for collocated grid
Halim Choi
Posts: n/a
(1) Yes, usually 1/Ap is linearly interpolated. (2) Yes, the cell face velocities at the previous iteration and time step should be stored in order to get solutions which are independent of under-relaxation factor and time. (3) Usually the mass fluxes (density*velocity*area) are stored at the cell-faces instead of the velocities. The cell face mass fluxes are used not only for the source term of pressure correction equation, but also for treating the convection terms of transport equations such as the momentum, energy and various turbulent quantities.

I hope this helps,

Halim Choi
  Reply With Quote

Old   September 20, 2004, 01:08
Default Re: momentum interpolation for collocated grid
Posts: n/a
Thank you for your reply.

As you stated, we need the velocity of last iteration and last time step to calculate the new cell face velocity. this need to store 2 velocity component (u and v) at each face in non-orthogonal grid. if i store only the mass flux at each face, i can not calculate the velocity components which are needed for the momentum interpolation. how should i treat it?

if we assume the same underrelaxation factor for u and v, it is possible to combine the equations and interpolate the contravariant velocity components directly on the cell face. so i will need only one velocity for each face. what is your suggestion for this trick?

Thanks again for your attention. Hadian
  Reply With Quote

Old   September 20, 2004, 04:14
Default Re: momentum interpolation for collocated grid
Halim Choi
Posts: n/a
Dear Hadian;

(1) Store the Cartesian velocity components (u,v) at the cell-centers and store mass flux only at the cell-faces. The velocity components which are needed for the momentum interpolation to get the cell-face mass flux (in nonorthogonal grid the cell-face mass flux is density*contravariant velocity component) are the cell-centered Cartesian velocity components and the mass flux at the previous time step and iteration. Combining the expression for the cell-face velocities from the momentum interpolation, you can calculate the contravariant velocities at the cell-faces. Instead of storing the mass flux at the cell-faces, you can store both Cartesian velocity components at the cell-faces. In this case, you need one more storage and a routine to calculate the cell-face mass flux which are needed to treat the convection terms in the transport equation.

(2) Yes. You got it. It is generally assumed that the under-relaxation factors for both cell-centered velocity components are same in the momentum interpolation method. If you do not like this, store both Cartesian velocity components at the cell-faces. Then, you do not have any problems, but it is messy and need more computer storage and programing.

Could you please read my paper "Use of the momentum interpolation method for numerical solution of incompressible flows in complex geometries : choosing cell face velocities" Numerical Heat Transfer, Part B, vol.23, pp.21-41, 1993.

Good luck

Halim Choi
  Reply With Quote

Old   December 25, 2009, 07:25
Default bc for colocated grid
Join Date: Jun 2009
Posts: 64
Rep Power: 16
jyothishkumar is on a distinguished road
Dear Halim Choi,

In the case of colocated grid, can you please tell me the boundary condition (right side, consider a one dimensional case) at the right side.
Basically i am using cell centered CV. Assume I have 7 points (cell centered). I have 1 and 7 as my ghost point. i have on the west face of point 2 my velocity inlet boundary condition. At 7 i have my pressure boundary condition. I dont know how to give the bc at the exit. I tried with some extrapolations at the exit but nothing worked out well.

In the momentum equation, for finding the convective quantity (F) whether i have to use average nodal velocity or directly the face velocity.
I am interpolating the ap(i) value at the faces. Please help me in this regard.


jyothishkumar is offline   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
On grid independent solution for pulsatile flow David FLUENT 5 March 25, 2022 03:33
momentum wall boundary condition in unstructured grid ztdep Main CFD Forum 0 August 9, 2011 21:44
"grid points" or "grid interface" ztdep Main CFD Forum 1 June 6, 2007 15:00
Grid Adaptation Suresh FLUENT 0 October 15, 2003 13:18
Numerical methods for discontinuous grid interfaces? Hansong Hang Main CFD Forum 12 September 16, 1998 22:26

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