• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.51-56
• /
• 2006

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are applied to two typical test cases of a translation test and a collapsing water column problem for the accurate capturing of fluid interfaces. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method, which is based on the finite-volume technique and fully conservative. The calculated results are found to show that SMART scheme gives the best performance with respect to accuracy and robustness.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.193-200
• /
• 1999

An unstructured implicit Euler solver is parallelized on a Cray T3E. Spatial discretization is accomplished by a cell-centered finite volume formulation using an unpwind flux differencing. Time is advanced by the Gauss-Seidel implicit scheme. Domain decomposition is accomplished by using the k-way N-partitioning method developed by Karypis. In order to analyze the parallel performance of the solver, flows over a 2-D NACA 0012 airfoil and a 3-D F-5 wing were investigated.

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

The Newton-Krylov method on the unstructured grid flow solver using the cell-centered spatial discretization oi compressible Euler equations is presented. This flow solver uses the reconstructed primitive variables to get the higher order solutions. To get the quadratic convergence of Newton method with this solver, the careful linearization of face flux is performed with the reconstructed flow variables. The GMRES method is used to solve large sparse matrix and to improve the performance ILU preconditioner is adopted and vectorized with level scheduling algorithm. To get the quadratic convergence with the higher order schemes and to reduce the memory storage. the matrix-free implementation and Barth's matrix-vector method are implemented and compared with the traditional matrix-vector method. The convergence and computing times are compared with each other.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.98-105
• /
• 1999

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

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

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.13.2-13.2
• /
• 2005

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.61.1-61.1
• /
• 2005

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.