• 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.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.96-101
• /
• 1998

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• Journal of the Korean earth science society
• /
• v.39 no.4
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.212-218
• /
• 1995

Turbulent flow and heat transfer characteristics around staggered tube bundles were studied using a non-orthogonal boundary fitted coordinate system and the low Reynolds .kappa. - .epsilon. turbulence model suggested by Lam and Bremhorst. The predicted flow characteristics for two tube pitches and tube arrangement showed good agreement with the experimental data except the strongly curved region. The predicted Nusselt number was compared with measurements obtained in the staggered rough bundles and it revealed the similar trend to measurements, but the location of the maximum and minimum heat transfer differed somewhat from the measurements.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.59-70
• /
• 1992

A Finite volume method for the computation of the two-dimensional, incompressible, steady, laminar Navier-Stokes equation is developed using a non-staggered grid system in a general curvilinear coordinate. The numerical pressure fluctuations, usually encountered when the non-staggered grid system is used, is suppressed by the momentum interpolation method. Flows around a NACA0012 foil section have been computed by the present method and the results show good agreements with other experimental and numerical ones.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.