• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Journal of the Korean Society of Safety
• /
• v.17 no.3
• /
• pp.50-55
• /
• 2002

In this numerical study, characteristics of swirl generation by the fan and selection of the location of the fan was studied theoretically by the modified TEACH code. The governing equations for the system are solved by means of the k dimensional version of the SIMPLE method and STAGGERED grid. From the present results, the optimal position of the fm is 0.625(h/H).

• Proceedings of the Korean Institute of Industrial Safety Conference
• /
• 2000.06a
• /
• pp.42-46
• /
• 2000

In this numerical study, characteristics of swirl generation by the fan and selection of the location of the fan was studied theoretically by the PHOENICS soft ware. The governing equations for the system are solved by means of the three dimensional version of the SIMPLE algorithm and STAGGERED grid. From the present results, the optimal position of the fan is 0.625($\ell$/L). Here we can survey the big swirl near the fan.

• Journal of Ship and Ocean Technology
• /
• v.7 no.4
• /
• pp.32-40
• /
• 2003

Numerical simulation code is developed to study the formation of liquid drops from a nozzle into another quiescent liquid. The Navier-Stokes equations for two immiscible, incompressible, Newtonian fluids are solved on a fixed, staggered grid of cylindrical axisymmetric coordinates. Interfacial motion is captured using a Front-Tracking Method. The time variation of interfacial shape simulated by the code is in excellent agreement with experiments. Simulation results show that the viscosity ratio affects the size of the satellite drops.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.03a
• /
• pp.50-58
• /
• 2003

We have developed a locally variable time-step scheme matching with discontinuous grids in the flute-difference method for the efficient simulation of seismic wave propagation. The first-order velocity-stress formulations are used to obtain the spatial derivatives using finite-difference operators on a staggered grid. A three-times coarser grid in the high-velocity region compared with the grid in the low-velocity region is used to avoid spatial oversampling. Temporal steps corresponding to the spatial sampling ratio between both regions are determined based on proper stability criteria. The wavefield in the margin of the region with smaller time-step are linearly interpolated in time using the values calculated in the region with larger one. The accuracy of the proposed scheme is tested through comparisons with analytic solutions and conventional finite-difference scheme with constant grid spacing and time step. The use of the locally variable time-step scheme with discontinuous grids results in remarkable saving of the computation time and memory requirement with dependency of the efficiency on the simulation model. This implies that ground motion for a realistic velocity structures including near-surface sediments can be modeled to high frequency (several Hz) without requiring severe computer memory

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.1
• /
• pp.122-130
• /
• 1998

The present work analyzes the pressure drop and heat transfer characteristics of the plate heat exchanger with staggered dimples. The flow is assumed to be constant property, three dimensional and laminar. A thermal boundary condition is uniform wall temperature and it is assumed that the flow is periodically fully developed. Elliptic grid generation is used for proper modelling of the internal tube geometry with dimples. Computations have been carried out for a variety of geometric parameters, Reynolds number, and Prandtl number. The friction factor ratio and the ratio of a module average Nusselt number are presented for the cases considered in this study. It is found that the distance between dimples has a substantial effect on the pressure drop and heat transfer.

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

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.51-61
• /
• 2008

The SMAC (Simplified Marker And Cell) algorithm is extended for an application to thermal non-equilibrium two-phase flows in light water nuclear reactors (LWRs). A two-fluid three-field model is adopted and a multi-dimensional unstructured grid is used for complicated geometries. The phase change and the time derivative terms appearing in the continuity equations are implemented implicitly in a pressure correction equation. The energy equations are decoupled from the momentum equations for faster convergence. The verification of the present numerical method was carried out against a set of test problems which includes the single and the two-phase flows. The results are also compared to those of the semi-implicit ICE method, where the energy equations are coupled with the momentum equation for pressure correction.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.37-50
• /
• 2016

A large eddy simulation (LES) of a turbulent channel flow is performed by using the third order low-storage Runge-Kutta method in time and second order finite difference formulation in space with staggered grid at a Reynolds number, $Re_{\tau}=590$ based on the channel half width, ${\delta}$ and wall shear velocity, $u_{\tau}$. To reduce the calculation cost of LES, algebraic wall model (AWM) is applied to approximate the near-wall region. The computation is performed in a domain of $2{\pi}{\delta}{\times}2{\delta}{\times}{\pi}{\delta}$ with $32{\times}20{\times}32$ grid points. Standard Smagorinsky model is used for subgrid-scale (SGS) modeling. Essential turbulence statistics of the flow field are computed and compared with Direct Numerical Simulation (DNS) data and LES data using no wall model. Agreements as well as discrepancies are discussed. The flow structures in the computed flow field have also been discussed and compared with LES data using no wall model.