• 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.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.1
• /
• pp.72-83
• /
• 2017

A thermo-fluid dynamic analysis was performed through the numerical simulation of a missile canister. Calculation was made in a fixed analytical volume and fully 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 parametric study was performed with the change of coolant flow rate. It could be found that the pressure on the canister top nonlinearly increased with the increase of coolant flow rate. Temperature and coolant distribution were closely related to the flow behavior in canister. Temperature on the canister bottom indicated a decrease being proportional to coolant flow rate in early times but after a specific time, the temperature increased with the tendency being reversed. In addition, the early part of temperature showed a fluctuating phenomenon because of the overall circulatory flow of mixture gas.

• 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.

• Wind and Structures
• /
• v.17 no.4
• /
• pp.361-377
• /
• 2013

Wall functions have been widely used in computational fluid dynamics (CFD) simulations and can save significant computational costs compared to other near-wall flow treatment strategies. However, most of the existing wall functions were based on the asymptotic characteristics of near-wall flow quantities, which are inapplicable in complex and non-equilibrium flows. A modified wall function is thus derived in this study based on flow over a plate at zero-pressure gradient, instead of on the basis of asymptotic formulations. Turbulent kinetic energy generation ($G_P$), dissipation rate (${\varepsilon}$) and shear stress (${\tau}_{\omega}$) are composed together as the near-wall expressions. Performances of the modified wall function combined with the nonlinear realizable k-${\varepsilon}$ turbulence model are investigated in homogeneous equilibrium atmosphere boundary layer (ABL) and flow around a 6 m cube. The computational results and associated comparisons to available full-scale measurements show a clear improvement over the standard wall function, especially in reproducing the boundary layer flow. It is demonstrated through the two case studies that the modified wall function is indeed adaptive and can yield accurate prediction results, in spite of its simplicity.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.14 no.4
• /
• pp.273-279
• /
• 2011

Turbulent flow analysis around a wind turbine blade was performed to evaluate the power performance of offshore wind turbine. Fluent package was utilized to solve the Reynolds-averaged Navier-Stokes equations in non-inertial rotating coordinates. The realizable k-$\varepsilon$ model was used for turbulence closure and the grid system combining structured and unstructured grids was generated. In the first, lift and drag forces of 2-D foil section were calculated and compared with existing experimental data for the validation. Then torque and thrust of the wind turbine blade having NACA 4-series sections were calculated with fixed pitch angle and rpm. Tip speed ratio was varied by changing wind speed. In the next, three kinds of end plate were attached at the tip of blade in order to increase the power of the wind turbine. Among them the end plate attached at the suction side of the blade was found to be most effective. Furthermore, performance analysis with tilt angle and rake was also performed.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.12-24
• /
• 2016

The standard ${\kappa}-{\varepsilon}$ and realizable ${\kappa}-{\varepsilon}$ models are adopted to improve the prediction performance on the recirculating flow. In this paper, the backward facing step flows are used to assess the prediction performance of the recirculation zone. The model constants of turbulence model are obtained by the experimental results and they have a different value according to the flow. In the case of an isotropic flow situation, decaying of turbulent kinetic energy should follow a power law behavior. In accordance with the power law, the coefficients for the dissipation rate of turbulent kinetic energy are not universal. Also, the other coefficients as well as the dissipation coefficient are not constant. As a result, a suitable coefficients can be varied according to each of the flow. The changes of flow over the backward facing step in accordance with model constants of the ${\kappa}-{\varepsilon}$ models show that the reattachment length is dependent on the growth rate(${\lambda}$) and the ${\kappa}-{\varepsilon}$ models can be improved the prediction performance by changing the model constants about the recirculating flow. In addition, it was investigated for the curvature correction effect of the ${\kappa}-{\varepsilon}$ models in the recirculating flow. Overall, the curvature corrected ${\kappa}-{\varepsilon}$ models showed an excellent prediction performance.

• Fire Science and Engineering
• /
• v.28 no.4
• /
• pp.35-43
• /
• 2014

The prediction performance of 9 model sets, which combine 3 turbulent models and 3 combustion models, was investigated numerically for turbulent partially-premixed jet flame. The standard ${\kappa}-{\varepsilon}$ (SKE), Realizable ${\kappa}-{\varepsilon}$ (RKE) and Reynolds stress model (RSM) were used as a turbulence model, and the eddy dissipation concept (EDC), steady laminar flamelet (SLF) and unsteady laminar flamelet model (ULF) were also adopted as a combustion model. The prediction performance of those 9 model sets was evaluated quantitatively and qualitatively for Sandia D flame of which flame structure was measured precisely. The flame length was predicted as, from longest to shortest, RSM > SKE > RKE, and the RKE predicted the flame length of the jet flame much shorter than experiment. The flame temperature was over predicted by the combination of RSM + SLF or RSM + ULF while the flame length obtained by RSM + SLF and RSM + ULF was well agreed with the experiment. The combination of SKE + SLF and SKE + ULF predicts well the flame length as well as the temperature distribution. The SKE turbulence model was most superior to the other turbulent models, and SKE + ULF showed the best prediction performance for the structure of turbulent partially-premixed jet flame.

• Journal of Power System Engineering
• /
• v.21 no.6
• /
• pp.5-12
• /
• 2017

The steady-state, incompressible and three-dimensional numerical analysis was carried out to investigate the internal flow fields characteristics according to wind tunnel contraction type. The turbulence model used in this study is a realizable $k-{\varepsilon}$ modified from the standard $k-{\varepsilon}$ model. As a results, the distribution of the axial mean velocity components along the central axis of the flow model is very similar to the ASME and BE types, and the cubic and cosine types. When the flow passes through the interior space of the analytical models, the flow resistance at the inlet of the plenum chamber is the largest at BS type contraction, but the smallest at cubic type contraction. The boundary layer thickness is the smallest in the cosine type contraction as the axial distance increases. The maximum turbulent kinetic energy in the test section is the smallest in the order of the contraction of cubic type and cosine type. Comprehensively, cubic type contraction is the best choice for wind tunnel performance, and cosine type contraction can be the next best solution.

• Journal of Aerospace System Engineering
• /
• v.13 no.5
• /
• pp.39-48
• /
• 2019

The turbulent flow in the rotor-stator mixer is based on shear characteristics generated by the interaction of the stator with the rotor rotating at high speed. In this study, the flow characteristics analysis of the unsteady state generated by the interaction of the rotor and the stator in the prototype model of the emulsion-fuel related mixer development was performed with the MRF and SMM by applying the ANSYS FLUENT $k-{\varepsilon}$ (RKE) turbulence model. The behavior and shear characteristics of the flow particles generated at the interface between the designed rotor and stator, and trends such as velocity distribution and turbulence eddy dissipation, were predicted and verified using the CFD analysis.