• The KSFM Journal of Fluid Machinery
• /
• v.17 no.6
• /
• pp.52-58
• /
• 2014

In the present work, design optimization of film-cooling hole array on the pressure side of high pressure turbine nozzle was conducted. There are four rows of fan-shaped film cooling holes on the nozzle pressure side surface and each row has a straight array of holes in the spanwise direction for baseline model. For design optimization, hole distributions in streamwise and spanwise directions for three rows of holes except first row are parameterized as a 2nd-order shape function. Three-dimensional compressible RANS equations are used for flow and thermal analysis around the nozzle surface and optimization technique using Design of Experiment, Kriging surrogate model and Genetic Algorithm is used. The results shows that averaged adiabatic wall temperature at the whole nozzle surface decreases about 2.7% and averaged film cooling effectiveness at the pressure side of nozzle increased about 8.2%.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.259-267
• /
• 2017

Flow fields around a KARI-11-180 airfoil, SDM and transonic body are numerically simulated by using an unstructured meshes based compressible flow solver developed at KAIST. RANS equations are solved to analyse the flow fields and Roe's FDS method is adopted to evaluate convective fluxes. Turbulence effect of the flow fields is modeled by a SA model, SST model and ${\gamma}-{\widetilde{Re}}_{{\theta}t}$ model. It is found that smaller drag coefficients are predicted for the KARI-11-180 airfoil when a transition phenomenon is considered and small deviations exist between CFD and EFD results. For the SDM, flow separation is observed at a leading edge and calculated aerodynamic properties show similar tendencies to experimental results. A shock wave on main wings of the transonic body is successfully captured by the present flow solver at a Mach number 0.9. Estimated pressure profiles by means of the present CFD method also agree well with those of wind tunnel results.

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

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

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

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.435-449
• /
• 2014

The effects of the gravity field and the free surface on the cavity shape and the drag are investigated through a numerical analysis for the steady supercavitating flow past a simple two-dimensional body underneath the free surface. The continuity and the RANS equations are numerically solved for an incompressible fluid using a $k-{\epsilon}$ turbulence model and a mixture fluid model has been applied for calculating the multiphase flow of air, water and vapor using the method of volume of fluid and the Schnerr-Sauer cavitation model. Numerical solutions have been obtained for the supercavitating flow about a two-dimensional $30^{\circ}$ wedge in wide range of depths of submergence and inflow velocities. The results are presented for the cavity shape, especially the length and the width, and the drag of the wedge in comparison with those of the case for the infinite fluid flow neglecting the gravity and the free surface. The influences of the gravity field and the free surface on the aforementioned quantities are discussed. The length and the width of the supercavity are reduced and the centerline of the cavity rises toward the free surface due to the effects of the gravity field and the free surface. The drag coefficient of the wedge, however, is about the same except for shallow depths of submergence. As the supercavitating wedge is approaching very close to the free surface, it is found the length and the width of a cavity are shorten even though the cavitation number is reduced. Also the present result suggests that, under the influence of the gravity field and the free surface, the length of the supercavity for a certain cavitation number varies and moreover is proportional to the inverse of the submergence depth Froude number.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.5
• /
• pp.355-363
• /
• 2013

Panel methods are essential tools for analyzing a fluid-flow problem around complex three dimensional bodies, but they lack ability to solve viscous effects. On the other hand, CFD methods are considered as powerful tools for analyzing fluid-flow characteristics including viscosity. However, they also have short falls, requiring more computing time and showing different results depending on the selection of turbulence models and grid systems. In this paper a hybrid scheme combining a panel method and a CFD method is suggested. The scheme adopts a panel method for far-field solution where viscous effects are negligible and a CFD method for the solution of RANS equations in near-field where viscous effects are relatively strong. The intermediate region between the far-field and near-field is introduced where the calculated field point velocities by the panel method are given as boundary velocities for the CFD method. To verify the scheme, calculated results, by a panel method, a CFD method and the hybrid scheme, for a two dimensional foil section are compared. The suggested hybrid scheme gives reasonable results, while computation time and memory can be dramatically reduced. By using the hybrid scheme efforts can be concentrated for the local flow near the leading and trailing edges, by providing more dense grid system, where detailed flow characteristics are required.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.5
• /
• pp.60-66
• /
• 2012

One of commonly physical phenomena encountered in pump sump systems in which its significant influence to the hydraulic performance of pump system plays an important role in the field of fluid engineering, is the appearance of free surface and submerged vortices. In this paper, a study of the vortices behavior and their formative mechanism of asymmetry is considered in this paper by using numerical approach. The Reynolds-Averaged Navier-Stokes (RANS) equations and k-omega Shear Stress Transport turbulence model used to describe the properties of turbulent flows, in company with VOF multiphase model, are implemented by Fluent code with multi-block structured grid system. In the numerical simulation, the calculated elevation of air-water interface and vortex core contours are used to classify visually surface vortices as well as submerged vortices. It is shown that the free surface vortex is identified by the concavity of liquid region from the free surface and swirling flow at that own plane. To investigate the distinctive behavior of these vortices corresponding to each given flow rate at the same water level, some numerical testing of them are considered here in such a manner that the flow pattern of surface vortex are obtained similarly to the obtained results from experiment. Furthermore, the influence due to the change of grid refinement and the variation of depth of the concavity are also considered in this paper. From that, these influential factors will be implemented to design a good pump sump with higher performance in the future.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.1-7
• /
• 2014

The present study numerically investigated the effect of the advance ratio on the wake characteristics of the marine propeller in the propeller open water test. Therefore, a wide range of the advance ratio(0.2${\kappa}-{\omega}$SST Model are considered. The three-dimensional vortical structures of tip vortices are visualized by the swirl strength, resulting in fast decay of the tip vortices with increasing the advance ratio. Furthermore, to better understanding of the wake evolution, the contraction ratio of the slip stream for different advance ratios is extracted from the velocity fields. Consequently, the slip stream contraction ratio decreases with increasing the advance ratio and successively the difference of the slip stream contraction ratio between J=0.2 and J=0.8 is about 0.1R.