• Nuclear Engineering and Technology
• /
• 제46권5호
• /
• pp.655-666
• /
• 2014

The CUPID code has been developed at KAERI for a transient, three-dimensional analysis of a two-phase flow in light water nuclear reactor components. It can provide both a component-scale and a CFD-scale simulation by using a porous media or an open media model for a two-phase flow. In this paper, recent advances in the CUPID code are presented in three sections. First, the domain decomposition parallel method implemented in the CUPID code is described with the parallel efficiency test for multiple processors. Then, the coupling of CUPID-MARS via heat structure is introduced, where CUPID has been coupled with a system-scale thermal-hydraulics code, MARS, through the heat structure. The coupled code has been applied to a multi-scale thermal-hydraulic analysis of a pool mixing test. Finally, CUPID-SG is developed for analyzing two-phase flows in PWR steam generators. Physical models and validation results of CUPID-SG are discussed.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.189-192
• /
• 2002

In this paper, the submarine model, called DARPA SUBOFF model, has been numerically analyzed to investigate the aerodynamic forces variation in terms of angle of attacks and yaw angles. The SUBOFF model is consisted of the three parts : axisymmetric body, fairwater, and four symmetric stern appendages. Three dimensional unsteady incompressible Wavier-Stokes equation was used on curvilinear multi-block grid system. To validate the present code, the SUBOFF tare hull and an ellipsoid at angle of attacks of $10^{\circ}\;and\;30^{\circ}$ were simulated and a good agreement with experiments was obtained. After the code validation, the flows over SUBOFF model were simulated at three different angle of attacks and yaw angles. The variation of aerodynamic forces in terms of angle of attack and yaw angle were calculated. Also, to understand the flow features around a submarine with variation of yaw and attack angle, the pressure contours and streamlines were plotted.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.592-595
• /
• 2011

A numerical investigation on a suction vortices, free vortices and subsurface vortices behavior in the model sump system with multi-intakes is performed A test model sump and piping system were designed based on Froude similitude for the prototype of the recommended structure layout by HI-9.B Standard for Pump Intake Design of the Hydraulic Institute. A numerical analysis of three dimensional multiphase flows through the model sump is performed by using the finite volume method of the CFX code with multi-block structured grid systems. A ${\kappa}-{\omega]$ ShearStressTransportturbulencemodelandthe Rayleigh-Plesset cavitation model are used for solving turbulence cavitating flow. From the numerical analysis, several types of vortices are reproduced and their formation, growing shedding and detailed vortex structures are investigated. To reduce abnormal vortices, an anti-vortex device is considered and its effect is investigated and discussed.

• Journal of Ship and Ocean Technology
• /
• 제4권3호
• /
• pp.21-32
• /
• 2000

In the present paper, the sloshing flow in the liquid holds of a large tanker is simulated using a numerical method. In the fluid domain, the three-dimensional Navier-Stokes equation with free surface is solved using a finite difference method, and the realistic shapes of multi holds are modeled including the internal members. The time-history of the tank motion is obtained using a time-domain program for ship motion. In order to computer the impulsive pressures on internal structures, a concept of buffer zone is adopted near the tank ceiling during impact occurrence. This study demonstrates that the global fluid motion in the multi liquid holds of ships and FPSO's can be simulated using the numerical method and the corresponding local pressure can be predicted with reasonable accuracy.

• International Journal of Aeronautical and Space Sciences
• /
• 제18권4호
• /
• pp.601-613
• /
• 2017

An adjoint-based error estimation and mesh adaptation study is conducted for two-dimensional viscous flows on unstructured hybrid meshes. The error in an integral output functional of interest is estimated by a dot product of the residual vector and adjoint variable vector. Regions for the mesh to be adapted are selected based on the amount of local error at each nodal point. Triangular cells in the adaptive regions are refined by regular refinement, and quadrangular cells near viscous walls are bisected accordingly. The present procedure is applied to single-element airfoils such as the RAE2822 at a transonic regime and a diamond-shaped airfoil at a supersonic regime. Then the 30P30N multi-element airfoil at a low subsonic regime with a high incidence angle (${\alpha}=21deg.$) is analyzed. The same level of prediction accuracy for lift and drag is achieved with much less mesh points than the uniform mesh refinement approach. The detailed procedure of the adjoint-based mesh refinement for the multi-element airfoil case show that the basic flow features around the airfoil should be resolved so that the adjoint method can accurately estimate an output error.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.86-94
• /
• 2009

The non-conservative form of momentum equations is often used for some two-phase flow codes instead of a conservative form because of numerical convenience. Another non-conservative form, so called, a semi-conservative form can improve the numerical solution of these codes maintaining the numerical convenience. It is close to the conservative form but still maintains the feature of the non-conservative form. A semi-conservative form of the momentum equations and a non-conservative form of the momentum equations are implemented in CUPID[1] code. The numerical results of the semi-conservative and the non-conservative forms are compared against analytical solutions and the solutions of the FLUENT code that uses the conservative form. The results clearly showed that the semi-conservative form of the momentum equations provides better solutions than the non-conservative form, especially for heterogeneous two-phase flows.

• International Journal of Ocean System Engineering
• /
• 제1권3호
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• 한국항공우주학회지
• /
• 제30권5호
• /
• pp.32-40
• /
• 2002

희박상태나 극소장치에 관련된 기체운동을 해석하는 문제가 최근 중요한 연구주제로 부각되고 있다. 잘 알려진 DSMC와 더불어 모우멘트 기법, Chapman-Enskog 기법으로 분류되는 고차 비평형 유동 해석모델들이 이 문제에 적용되어 왔다. 본 연구에서는 Eu의 일반유체역학을 근간으로 이원자 기체에 관한 고차 해석모델을 개발하고자 한다. 회전 비평형 효과는 기체의 용적 점성계수에 관한 초과 수직응력을 고려하여 감안하였다. 개발된 계산모델을 일차원 충격파 내부구조와 단순 형상 외부의 희박 극초음속 유동장 해석에 적용하였다. 충격파 내부구조 및 전단유동 해석을 통해 회전 비평형에 의한 용적 점성계수 효과가 중요함을 확인하였다. 충격파 내부구조에 관한 이론적 예측이 실험과 잘 일치함도 확인하였다.

• 한국전산유체공학회지
• /
• 제18권1호
• /
• pp.83-90
• /
• 2013

OpenMP is becoming more and more useful as a simple parallel processing paradigm on SMP (Shared Memory Multi-Processors) computing environment with the development of multi-core processors. However, very few data is available publically regarding the OpenMP performance in CFD (Computational Fluid Dynamics). In the present study a CFD test suite is prepared for the performance evaluation of OpenMP on various multi-core systems. The test suite is composed of two-dimensional numerical simulations for inviscid/viscous and reacting/non-reacting flows using three different levels of grid systems. One to five test runs were carried out on various systems from dual-core dual threads to 16-core 32-threads systems by changing the number of threads engaged for each test up to 80. The results exhibit some interesting results and the lessons learned from the tests would be quite helpful for the further use of OpenMP for CFD studies using multi-core processor systems.