• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.2
• /
• pp.75-83
• /
• 2016

Since a turbomachine has complex flow characteristics, which are caused by adverse pressure gradient and high speed motion, an elaborate turbulence model is needed to accurately predict the flow. Some turbulence models such as an algebraic or a two-equation eddy viscosity model have been used for in-house RANS-code, but it is difficult to obtain good result for several complex flows. In this study, Durbin's V2-F turbulence model, which has been known for better prediction for severe flow separation, is applied to T-Flow. It was validated for simple cases such as channel and compressor cascade, and its applicability to turbomachinery was shown by analyzing internal flow of a single rotor. As a result, the V2-F turbulence model shows better blade surface pressure distribution than the one-and-two equation turbulence model.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Wind and Structures
• /
• v.27 no.2
• /
• pp.101-109
• /
• 2018

Wall jet flow exists widely in engineering applications, including the simulation of thunderstorm downburst outflows, and has been investigated extensively by both experimental and numerical methods. Most previous studies focused on the scaling laws and self-similarity, while the effect of lip thickness and external stream height on mean velocity has not been examined in detail. The present work is a numerical study, using steady Reynolds-Averaged Navier Stokes (RANS) simulations at a Reynolds number of $3.5{\times}10^4$, of a turbulent plane wall jet with an external stream to investigate the influence of the wall jet domain on downstream development of the flow. The comparisons of flow characteristics simulated by the Reynolds stress turbulence model closure (Stress-omega, SWRSM) and experimental results indicate that this model may be considered reasonable for simulating the wall jet. The confined wall jet is further analyzed in a parametric study, with the results compared to the experimental data. The results indicate that the height and the width of the wind tunnel and the lip thickness of the jet nozzle have a great effect on the wall jet development. The top plate of the tunnel does not confine the development of the wall jet within 200b of the nozzle when the height of the tunnel is more than 40b (b is the height of jet nozzle). The features of the centerline flow in the mid plane of the 3D numerical model are close to those of the 2D simulated plane wall jet when the width of the tunnel is more than 20b.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.15-23
• /
• 2016

In the present study, two phase flows around a projectile vertically launched from an underwater platform have been numerically investigated by using a three dimensional multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. The relative motion between the platform and projectile was described by six degrees of freedom equations of motion with Euler angles and a chimera technique. The propulsive power of the projectile was modeled as the fluid force acting on the lower surface of the body by the compressed air emitted from the underwater platform. Various flow conditions were considered to analyze the fluid-dynamics motion parameters of the projectile. The water level of platform and the current speed around the projectile were the main parametric variables. The numerical calculations were conducted up to 0.75sec in physical time scale. The dynamics tendency of the projectile was almost identical with respect to the water level variation due to the constant buoyancy term. The moving speed of the projectile along the vertical axis inside the platform decreased when the current speed increased. This is because the inflow from outside of the platform impeded development of the compressed air emitted from the floor surface of the launch platform. As a result, the fluid force acting on the lower surface of the projectile decreased, and injection time of the projectile from the platform was delayed.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.86-92
• /
• 2014

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Advances in aircraft and spacecraft science
• /
• v.7 no.6
• /
• pp.537-551
• /
• 2020

Investigation of leading edge impingement cooling for first stage rotor blades in an aero-engine turbine, its effect on rotor temperature and trailing edge wake loss have been undertaken in this study. The rotor is modeled with the nozzle for attaining a more accurate simulation. The rotor blade is hollowed in order for the coolant to move inside. Also, plenum with the 15 jet nozzles are placed in it. The plenum is fed by compressed fresh air at the rotor hub. Engine operational and real condition is exerted as boundary condition. Rotor is inspected in two states: in existence of cooling technique and non-cooling state. Three-dimensional compressible and steady solutions of RANS equations with SST K-ω turbulent model has been performed for this numerical simulation. The results show that leading edge is one of the most critical regions because of stagnation formation in those areas. Another high temperature region is rotor blade tip for existence of tip leakage in this area and jet impingement cooling can effectively cover these regions. The rotation impact of the jet velocity from hub to tip caused a tendency in coolant streamlines to move toward the rotor blade tip. In addition, by discharging used coolant air from the trailing edge and ejecting it to the turbines main flow by means of the slot in trailing edge, which could reduce the trailing edge wake loss and a total decrease in the blade cooling loss penalty.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.6
• /
• pp.46-53
• /
• 2016

It is required to predict a limit cycle amplitude controlled by system's nonlinear behavior as well as an eigen-frequency and initial growth rate of instabilities under the linear motions, in order to fully understand combustion instabilities in a lean premixed gas turbine combustor. Special focus of the current work is placed on the limit cycle amplitude prediction using flame describing function(FDF) where the ratio of a heat release fluctuation to a given flow perturbation is expressed as a function of frequency and amplitude. In this study, the CFD modeling work based on RANS is carried out to obtain FDF, which makes that the nonlinear thermo-acoustic model is successfully developed for predicting the limit cycle amplitude of the combustion instability.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.1905-1915
• /
• 2019

This paper investigated the influence of simple support girds on flow, irrespective of having mixing vanes, in a 1 × 3 array rod bundle by using CFD methodology and the most accurate turbulence model which could reflect the actual physics of the flow was determined. In this context, a CFD model was created simulating the experimental studies on a single-phase flow [1] and the results were compared with the experimental data. In the first part of the study, influence of mesh was examined. Tetra, hybrid and poly type meshes were analyzed and convergence study was carried out on each in order to determine the most appropriate type and density. k - ε Standard and RSM LPS turbulence models were used in this section. In the second part of the study, the most appropriate turbulence model that could reflect the physics of the actual flow was investigated. RANS based turbulence models were examined using the mesh that was determined in the first part. Velocity and turbulence intensity results obtained on the upstream and downstream of the spacer grid at -3d_h, +3d_h and +40d_h locations were compared with the experimental data. In the last section of the study, the behavior of flow through the spacer grid was examined and its prominent aspects were highlighted on the most appropriate turbulence model determined in the second part. Results of the study revealed the importance of mesh type. Hybrid mesh having the largest number of structured elements performed remarkably better than the other two on results. While comparisons of numerical and experimental results showed an overall agreement within all turbulence models, RSM LPS presented better results than the others. Lastly, physical appearance of the flow through spacer grids revealed that springs has more influence on flow than dimples and induces transient flow behaviors. As a result, flow through a simple support grid was examined and the most appropriate turbulence model reflecting the actual physics of the flow was determined.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.875-882
• /
• 2019

The effects of Triangle Groove Strips (TGS) on Vortex-induced Vibration (VIV) suppression of marine riser are numerically investigated using Computational Fluid Dynamics (CFD) method. The range of Reynolds number in simulations is 4.0 × 10⁴ < Re < 1.2 × 10⁵. The two-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS) equations and Shear Stress Transport (SST) k-ω turbulence model are used to calculate the flow around marine riser. The Newmark-β method is employed for evaluating the structure dynamics of marine riser. The effect of the height ratio (ε) of TGS on VIV suppression is evaluated. The amplitude responses, frequency responses, vortex patterns and the flow around the structures are discussed in detail. With the increase of the height ratio of TGS, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When ε=0.04, the suppression effect of TGS is the best. Compared with the VIV responses of smooth marine riser, the amplitude ratio is reduced by 38.9%, the peak of the lift coefficient is reduced by 69% and the peak of the drag coefficient is reduced by 40% when Re=6.0 × 10⁴. With the increase of Reynolds number, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When the Reynolds number is 7.0 × 10⁴, the amplitude ratio can be reduced by 40.1%. As to the large-amplitude vibration cases, the TGS show nice suppression effect on VIV.