• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.18 no.7
• /
• pp.578-586
• /
• 2006

The present study has been conducted to examine the effect of obstacles around a cooling tower and an air-guide to prevent recirculation. In order to analyze the interaction between external flow and cooling tower exit flow, the external region as well as the cooling, tower are included in computational domain. Two dimensional analysis is performed using the finite volume method with non-orthogonal and unstructured grid system. The standard ${\kappa}-{\varepsilon}$ turbulence model is used. To investigate the recirculation phenomena, flow and temperature fields are calculated with three approaches such as, the distance between cooling tower and obstacle, the allocated geometrical type, and the effect of height of obstacle. In addition, the air-guide is considered in the current computation. The mean recirculation rate increases with the height of obstacle. The effect of air-guide to reduce the mean recirculation rate is obviously observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.3
• /
• pp.254-263
• /
• 2017

The cavitation erosion remains an industrial issue for many applications. This paper deals with the cavitation intensity, which can be described as the fluid mechanical loading leading to cavitation damage. The estimation of this quantity is a challenging problem both in terms of modeling the cavitating flow and predicting the erosion due to cavitation. For this purpose, a numerical methodology was proposed to estimate cavitation intensity from 3D unsteady cavitating flow simulations. CFD calculations were carried out using Code_Saturne, which enables U-RANS equations resolution for a homogeneous fluid mixture using the Merkle's model, coupled to a $k-{\varepsilon}$ turbulence model with the Reboud's correction. A post-process cavitation intensity prediction model was developed based on pressure and void fraction derivatives. This model is applied on a flow around a hydrofoil using different physical (inlet velocities) and numerical (meshes and time steps) parameters. The article presents the cavitation intensity model as well as the comparison of this model with experimental results. The numerical predictions of cavitation damage are in good agreement with experimental results obtained by pitting test.

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

The finite volume based multi-block RANS code, WAVIS developed at KRISO, is used to simulate the turbulent flows around a submarine with the realizable $\textsc{k}-\varepsilon$ turbulence model. RANS methods are verified and validated at the level of validation uncertainty 1.54% of the stagnation pressure coefficient for the solution of the turbulent flows around SUBOFF submarine model without appendages. Another SUBOFF configuration, axisymmetric body with four identical stem appendages, is also computed and validated with the experimental data of the nominal wake and hydrodynamic coefficients. The hydrodynamic forces and moments for SUBOFF model and a practical submarine are predicted at several drift and pitch angles. The computed results are in extremely good agreement with experimental data. Furthermore, it is noteworthy that all the computations at the present study were carried out in a PC and the CPU time required for 2.8 million grids was about 20 hours to get fully converged solution. The current study shows that CFD can be a very useful and cost effective tool for the prediction of the hydrodynamic performance of a submarine in the basic design stage.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.2
• /
• pp.48-59
• /
• 2017

In this study, an analysis of the thermo-fluid dynamic and missile-motion performance was carried out through a numerical simulation inside the missile canister. Calculation was made in an analytical volume using dynamic grid and evaporated water was used as a coolant. To analyze the interaction among the hot gas, coolant, and mixture flow, Realizable $k-{\varepsilon}$ turbulence and VOF (Volume Of Fluid) model were chosen and a parametric study was performed with the change of coolant flow rate. As a result of the analysis, pressure of the canister showed a large difference depending on the presence or absence of the coolant, and also showed a dependancy on the amount of coolant. Velocity and acceleration were dependent on the canister pressure.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped SMEP (Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• Journal of Biomedical Engineering Research
• /
• v.13 no.2
• /
• pp.87-96
• /
• 1992

In thIns paper, a numerical simulation of steady laminar and turbulent flow in a two dimensional model for the total artificial heart is'presented. A trlleaflet polyurethane valve was simulated at the outflow orifice while the Inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-$\varepsilon$-E model. The SIMPLER algo rithm was used to solve the problem in primitive variables. The numerical solutions of the slulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the Inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller wlth the flap valve than with the trlleaflet valve. These resu1ts also suggest a correlation be- tween high turbulent stresses and the presence of thrombus In the vicinity of the valves in the total artificial hearts. The computed velocity vectors and trubulent stresses were comparable with previ ously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solo talons suggests that geometries similar to the flap valve(or a tilting disc valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.6
• /
• pp.750-757
• /
• 2016

This study presents the wave run-up height around single and multiple surface-piercing cylinders according to wave period and steepness. In order to simulate 3D incompressible viscous two-phase turbulent flow, the present study employed a volume of fluid (VOF) method with realizable $k-{\varepsilon}$ turbulence model based on commercial Computational Fluid Dynamics (CFD) software, "STAR-CCM". The wave periods at model scale were 1.269s and 1.692s for a single cylinder and 1.716s for multiple cylinders. In each case, wave steepness of has 1/30 and 1/16 were used, respectively. Consequently, the results for wave run-up height with regard to wave steepness and period were compared with those of relevant previous experimental studies. The numerical simulation results showed a good qualitative agreement with experiments.

• Journal of Korean Society for Atmospheric Environment
• /
• v.21 no.5
• /
• pp.525-535
• /
• 2005

Using a three-dimensional computational fluid dynamics (CFD) model with the $k-{\varepsilon}$ turbulence closure scheme based on the renormalization group theory, flow regimes in urban street canyons are classified according to the building and street aspect ratios. The transition between skimming flow (SF) and wake interference flow (WIF) is determined with the size of double-eddy circulation generated behind the upwind building. The transition between WIF and isolated roughness flow (IRF) is determined with the flow reattachment distance from the upwind building. The critical aspect ratios at which the flow transition occurs are found and compared with those in previous studies. The results show that the flow-regime classification method used in this study is quite reasonable and that the values of the critical aspect ratios are generally consistent with those in fluid experiments or large-eddy simulation. The regression equation describing a relation between the building and street aspect ratios at the flow-regime transition is presented.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.1
• /
• pp.9-16
• /
• 2012

This paper presents the numerical results of a simulation of the free surface flow around a blunt bow ship model and focuses on the validation of the proposed method with a brief investigation of the relation between the resistance and free surface behavior. A finite volume method based on the Reynolds Averaged Navier-Stokes (RANS) approach is used to solve the governing flow equations, where the free surface, including wave breaking,is captured by using a two-phase Level-Set (LS) method. For turbulence closure, a two equation k-${\varepsilon}$ model with the standard wall function technique is used. Finally, the numerical results are compared with the available experimental data, showing good agreement.