• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• 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년도 학술대회
• /
• 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년 추계학술대회논문집
• /
• 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.

• 한국항공우주학회지
• /
• 제32권8호
• /
• pp.1-11
• /
• 2004

본 연구에서는 비정렬 격자계에서 가장 많이 쓰이는 근사 해법 중에 하나인 LU 기법의 Navier-Stokse 방정식에 대한 수렴성 및 안정성에 관한 연구를 수행하였다. 적절한 스칼라 모델 방정식을 사용하여 LU 기법이 갖는 고유한 특성에 관한 해석적 논의를 수행하였으며, 이를 Navier-Stokes 방정식으로 확장하여 해석하였다. 그 결과 LU 기법의 강성도는 격자 종횡비가 높아짐에 띠라, 그리고 격자 레이놀즈 수 감소함에 따라 증가하게 된다. 또한 내부반복계산을 통해서 이러한 강성도가 부분적으로 극복될 수 있음을 보였으며, 평판 난류 유동 해석을 통해서 해석 결과를 검증하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• 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년 추계학술대회논문집
• /
• 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년도 춘계 학술대회논문집
• /
• 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년도 추계 학술대회논문집
• /
• 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년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• 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년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• 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.