• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.61-66
• /
• 2005

Natural convection flows in a cubical, air-filled cavity that has one pair of opposing faces isothermal at different temperatures, Th and Tc, the remaining faces having a linear variation from Tc to Th are numerically simulated by a new solution code(PowerCFD) using unstructured cell-centered method. Solutions are obtained for configurations with a Rayleigh number as high as 105 and three inclination angles ${\theta}$ of the isothermal faces from horizontal: namely ${\theta}=0$, 45 and $90^{\circ}$. Interesting features are presented in detail and comparisons are made with benchmark solutions and experimental results found in the literature. It is found that the code is capable of producing accurately the nature of the laminar convection in a cubical, air-filled cavity with differentially heated walls.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.71-78
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.71-78
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.1-11
• /
• 2004

A comprehensive study has been made for the investigation of the convergence and stability characteristics of the LU scheme for solving the Navier-Stokes equations on unstructured meshes. For this purpose the characteristics of the LU scheme was initially studied for a scalar model equation. Then the analysis was extended to the Navier-Stokes equations. It was shown that the LU scheme has an inherent stiffness in the streamwise direction. This stiffness increases when the grid aspect ratio becomes high and the cell Reynolds number becomes small. It was also shown that the stiffness related to the grid aspect ratio can be effectively eliminated by performing proper subiteration. The results were validated for a flat-plate turbulent flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.65-70
• /
• 2009

In this paper a hybrid grid generation program for general 2-D region is introduced. The program is developed by using JAVA programming language, and it can be used either as an application program on a local computer or as an applet in the network environment. The hybrid grid system for a 2-D problem means a combination of triangular cells and quadrilateral cells, and it can offer both of the high flexibility of triangular cells and the high accuracy and efficiency of structured-type quadrilateral cells. To accommodate a quadrilateral-cell region and a triangular-cell region into one computational domain, it is importance to take good care of the interface between two different regions so that overall good grid quality can be maintained. In this research advancing layer method(ALM) augmented by elliptic smoothing method is used for the quadrilateral-cell region and advancing front method(AFM) is used for the triangular-cell region. A special treatment technique for the interface between those two regions is also developed. The interface treatment technique is basically to prevent the propagation of small cell size due to ALM method into the triangular region and maintain the smooth transition of cell-size scale between two different regions. By applying current technique high-quality hybrid grids for general 2-D regions can be easily generated, and typical grid generation results and flow solutions are demonstrated.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.211-217
• /
• 2009

The numerical program has been developed for the purpose of the complicate geometries application using CIP method. The non-staggered, non-orthogonal, and unstructured grid system can be also used for the various geometries in the program. For validating CIP solver, the lid-driven cavity flow and solitary wave propagation flow are carried out. Test results show a good agreement with the verified results. The dynamic solver was used for the behavior of moving body. Interface process between the two solvers is introduced. The research was performed on the flow problem around torpedo and log and the flow problem in a tank in order to analyze the three phase flow problem Although the comparison to the verified results was not quantitatively performed, the trend of the results was reasonable.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.60-65
• /
• 2000

An aerodynamic design method has been developed by using a three-dimensional unstructured Euler code and an adjoint code with a discrete approach. The resulting adjoint code is applied to a wing design problem of super-sonic transport with a wing-body-nacelle configuration. Hicks-Henne shape functions are adopted far the surface geometry perturbation, and the elliptic equation method is employed fer the interior grid modification during the design process. Interior grid sensitivities are neglected except those for design parameters associated with nacelle translation. The Sequential Quadratic Programming method is used to minimize the drag with constraints on the lift and airfoil thickness. Successful design results confirm validity and efficiency of the present design method.

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

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.175-177
• /
• 2003

The convergence characteristics of the LV scheme for the Euler equations have been investigated by using the Von Neumann stability analysis. The results indicated that the convergence rate is governed by a specific combination of CFD parameters. Based on this insight, it is shown that the convergence characteristics of the LV scheme is not deteriorated at any grid aspect-ratio as long as the local time step is defined based on the parameter combination. The numerical results demonstrated that this time step definition provide a uniform convergence for grid aspect-ratios between one to$1{\times}10^{4}$.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.156-157
• /
• 2003

A simple method is proposed to split the flux vector of the Euler equations by introducing two artificial wave speeds. The direction of wave propagation can be adjusted by these two wave speeds. This idea greatly simplifies the upwinding, and leads to a new family of upwind schemes. Numerical flux function for multi-dimensional Euler equations is formulated for any grid system, structured or unstructured. A remarkable simplicity of the scheme is that it successfully achieves one-sided approximation for all waves without recourse to any matrix operation. Moreover, its accuracy is comparable with the exact Riemann solver. For 1-D Euler equations, the scheme actually surpasses the exact solver in avoiding expansion shocks without any additional entropy fix. The scheme can exactly resolve stationary contact discontinuities, and it is also freed of the carbuncle problem in multi­dimensional computations.