• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.57-64
• /
• 2015

This paper describes a volume of fluid(VOF) method, mRHRIC for the simulation of free surface flows around ship hulls and provides its validation against benchmark test cases. The VOF method is developed on the basis of RHRIC method developed by Park et al. that uses high resolution differencing schemes to algebraically preserve both the sharpness of interface and the boundedness of volume fraction. A finite volume method is used to solve the governing equations, while the realizable ${\kappa}-{\varepsilon}$ model is used for turbulence closure. The present numerical results of the resistance performance tests for DTMB5415 and KCS hull forms show a good agreement with available experimental data and those of other free surface methods.

• Journal of Energy Engineering
• /
• v.24 no.1
• /
• pp.149-163
• /
• 2015

In this study, we conducted study on the confirmation of thermal-hydraulic safety for Mast assembly with Computational Fluid Dynamics(CFD) analysis. Before performing the natural convection analysis for the Mast assembly by using CFD code, we validated the CFD code against two benchmark natural convection data for the evaluation of turbulence models and confirmation of its applicability to the natural convection flow. From the first benchmark test which was performed by Betts et al. in the simple rectangular channel, we selected standard k-omega turbulence model for natural convection. And then, calculation performance of CFD code was also investigated in the sub-channel of rod bundle by comparing with PNL(Pacific Northwest Laboratory) experimental data and prediction results by MATRA and Fluent 12.0 which were performed by Kwon et al.. Finally, we performed main natural convection analysis for fuel assembly inside the Mast assembly by using validated turbulence model. From the calculation, we observed stable natural circulation flow between the mast assembly and pool side and evaluated the thermal-hydraulic safety by calculating the departure from nucleate boiling ratio.

• International Journal of High-Rise Buildings
• /
• v.1 no.2
• /
• pp.87-97
• /
• 2012

The indoor environments in high-rise buildings are generally well enclosed by defined boundary conditions. Here, a numerical simulation method based on the Lattice Boltzmann method (LBM), which aims to model and simulate the turbulent flow accurately in an enclosed environment, and its comparison with traditional computational fluid dynamics (CFD) results, are presented in this paper. CFD has become a powerful tool for predicting and evaluating enclosed airflows with the rapid advance in computer capacity and speed, and various types of CFD turbulence modeling and its application and validation have been reported. The LBM is a relatively new method; it involves solving of the discrete Boltzmann equation to simulate the fluid flow with a collision model instead of solving Navier-Stokes equations. In this study, the LBM-based scheme of flow pattern and particle dispersion analyses are validated using the benchmark test case of two- and three-dimensional and isothermal conditions (IEA/Annex 20 case); the prediction accuracy and advantages are also discussed by comparison with the results of CFD.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1318-1325
• /
• 2017

Direct Numerical Simulation (DNS) serves as an irreplaceable tool to probe the complexities of multiphase flow and identify turbulent mechanisms that elude conventional experimental measurement techniques. The insights unlocked via its careful analysis can be used to guide the formulation and development of turbulence models used in multiphase computational fluid dynamics simulations of nuclear reactor applications. Here, we perform statistical analyses of DNS bubbly flow data generated by Bolotnov ($Re_{\tau}=400$) and LueTryggvason ($Re_{\tau}=150$), examining single-point statistics of mean and turbulent liquid properties, turbulent kinetic energy budgets, and two-point correlations in space and time. Deformability of the bubble interface is shown to have a dramatic impact on the liquid turbulent stresses and energy budgets. A reduction in temporal and spatial correlations for the streamwise turbulent stress (uu) is also observed at wall-normal distances of $y^+=15$, $y/{\delta}=0.5$, and $y/{\delta}=1.0$. These observations motivate the need for adaptation of length and time scales for bubble-induced turbulence models and serve as guidelines for future analyses of DNS bubbly flow data.

• Wind and Structures
• /
• v.13 no.2
• /
• pp.169-189
• /
• 2010

This paper describes the use of coupled Computational Fluid Dynamics (CFD) and Rigid Body Dynamics (RBD) in modelling the aerodynamic behaviour of wind-borne plate type objects. Unsteady 2D and 3D Reynolds Averaged Navier-Stokes (RANS) CFD models are used to simulate the unsteady and non-uniform flow field surrounding static, forced rotating, auto-rotating and free-flying plates. The auto-rotation phenomenon itself is strongly influenced by vortex shedding, and the realisable k-epsilon turbulence modelling approach is used, with a second order implicit time advancement scheme and equal or higher order advection schemes for the flow variables. Sequentially coupling the CFD code with a RBD solver allows a more detailed modelling of the Fluid-Structure Interaction (FSI) behaviour of the plate and how this influences plate motion. The results are compared against wind tunnel experiments on auto-rotating plates and an existing 3D analytical model.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• Journal of the Korean Society of Safety
• /
• v.33 no.3
• /
• pp.92-98
• /
• 2018

The flow of cooling water in a passive containment cooling system (PCCS), used to remove heat released in design basis accidents from a concrete containment of light water nuclear power plant, was conducted in order to investigate the thermo-fluid equilibrium among many parallel tubes of PCCS. Numerical simulations of the subcooled boiling flow within a coolant loop of a PCCS, which will be installed in innovative pressurized-water reactor (PWR), were conducted using the commercially available computational fluid dynamics (CFD) software ANSYS-CFX. Shear stress transport (SST) and the RPI model were used for turbulence closure and subcooled flow boiling, respectively. As the first step, the simplified geometry of PCCS with 36 tubes was modeled in order to reduce computational resource. Even and uneven thermal loading conditions were applied at the outer walls of parallel tubes for the simulation of the coolant flow in the PCCS at the initial phase of accident. It was observed that the natural circulation maintained in single-phase for all even and uneven thermal loading cases. For uneven thermal loading cases, coolant velocity in each tube were increased according to the applied heat flux. However, the flows were mixed well in the header and natural circulation of the whole cooling loop was not affected by uneven thermal loading significantly.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.04a
• /
• pp.652-653
• /
• 2009

In this study, flow-induced vibration (FIV) analyses have been conducted for a 3D compressor blade model. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed dynamic responses of designed compressor blades. Fluid domains are modeled using the computational grid system with local grid deforming and remeshing techniques. Reynolds-averaged Navier-Stokes equations with $\kappa-\varepsilon$ turbulence model are solved for unsteady flow problems of the rotating compressor model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D compressor blade for fluid-structure interaction (FSI) problems. Detailed dynamic responses and instantaneous pressure contours on the blade surfaces considering flow-separation effects are presented to show the multi-physical phenomenon of the rotating compressor blade.

• 유체기계공업학회:학술대회논문집
• /
• 2005.12a
• /
• pp.244-249
• /
• 2005

Flow distribution and pressure drop analysis for an inlet plenum of a Pebble Bed Modular Reactor (PBMR) have been performed using Computational Fluid Dynamics. Three-dimensional Navier-Stokes equations have been solved in conjunction with $k-{\epsilon}$ model as a turbulence closure. Non-uniformity in flow distribution is assessed for the reference case and parametric studies have been performed for rising channels diameter, Reynolds number and angle between the inlet ports. Also, two different shapes of the inlet plenum namely, rectangular shape and oval shape, have been analysed. The relative flow mal-distribution parameter shows that the flow distribution in the rising channels for the reference case is strongly non-uniform. As the rising channels diameter decreases, the uniformity in the flow distribution as well as the pressure drop inside the inlet plenum increases. Reynolds number is found to have no effect on the flow distribution in the rising channels for both the shapes of the inlet plenum. The increase in angle between the inlet ports makes the flow distribution in the rising channels more uniform.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.6
• /
• pp.551-559
• /
• 2009

In this study, flow-induced vibration(FIV) analyses have been conducted for a 3D compressor blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responses of designed compressor blades. Fluid domains are modeled using the computational grid system with local grid deforming and remeshing techniques. Reynolds-averaged Navier-Stokes equations with $\kappa-\epsilon$ turbulence model are solved for unsteady flow problems of the rotating compressor model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D compressor blade for fluid-structure interaction(FSI) problems. Detailed dynamic responses and instantaneous pressure contours on the blade surfaces considering flow-separation effects are presented to show the multi-physical phenomenon of the rotating compressor blade.