• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.387-393
• /
• 2019

In this study, the free motion of a riser due to vortex shedding was numerically simulated with Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) turbulence models. A numerical simulation program was developed by applying the Rhie-Chow interpolation method to the pressure correction of the OpenFOAM standard solver pimpleDyMFoam. To verify the developed program, the vortex shedding around the fixed riser at Re = 3900 was calculated, and the results were compared with the existing experimental and numerical data. Moreover, the vortex-induced vibration of a riser supported by a linear spring was numerically simulated while varying the spring constant. The results are compared with published direct numerical simulation (DNS) results. The present calculation results show that the numerical method is appropriate for simulating the vortex-induced motion of a riser, including lock-in phenomena.

• Journal of computational fluids engineering
• /
• v.6 no.1
• /
• pp.1-13
• /
• 2001

The natural convection heat transfer and solidification in a stratified pool are studied. The flow and heat transfer characteristics in a heat generating pool are compared between single-layered and double-layered pools. And local Nusselt number distributions on outer walls are obtained to consider thermal loads on a vessel wall. The cooling and solidification of Al₂O₃/Fe melt in a hemispherical vessel are simulated to study the mechanism of heat transfer and temperature distribution. A unstructured mesh is chosen for this study because of the non-orthogonality originated from the boundaries of double-layered pool. Interface between the layers is modeled to be fixed. With this assumption mass flux across the interface is neglected, but shear force and heat flux are considered by boundary conditions. The colocated cell-centered finite volume method is used with the Rhie-Chow interpolation to compute cell face velocity. To prevent non-physical solutions near walls in case body force is large the wall pressure is extrapolated by the way to include body force. The numerical solutions calculated by current method show that averaged downward heat flux of the double-layered pool increases compared to single-layered pool and maximum temperature occurs right below the interface of the layers.