# momentum interpolation for collocated grid

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

 September 14, 2004, 09:11 momentum interpolation for collocated grid #1 Hadian Guest   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. Regards, Hadian

 September 19, 2004, 20:05 Re: momentum interpolation for collocated grid #2 Halim Choi Guest   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

 September 20, 2004, 01:08 Re: momentum interpolation for collocated grid #3 Hadian Guest   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

 September 20, 2004, 04:14 Re: momentum interpolation for collocated grid #4 Halim Choi Guest   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

 December 25, 2009, 08:25 bc for colocated grid #5 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 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. thanks jyothish

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ztdep Main CFD Forum 0 August 9, 2011 21:44 ztdep Main CFD Forum 1 June 6, 2007 15:00 David FLUENT 4 June 8, 2005 16:20 Suresh FLUENT 0 October 15, 2003 13:18 Hansong Hang Main CFD Forum 12 September 16, 1998 22:26

All times are GMT -4. The time now is 09:39.