• Kyungpook Mathematical Journal
• /
• v.61 no.4
• /
• pp.831-843
• /
• 2021

In this paper we look at the three dimensional stagnation point flow problem over a rough rotating disk. We study the theoretical behaviour of the stagnation point flow, or forced flow, in the presence of a slip factor in which convective instability stationary modes appear. We make a numerical investigation of the effects of slip on the behaviour of the flow components of the stagnation point flow where the disk is rough. We provide, for the first time in the literature, a complete convective instability analysis and an energy analysis. Suitable similarity transformations are used to reduce the Navier-Stokes equations and the continuity equation into a system of highly non-linear coupled ordinary differential equations, and these are solved numerically subject to suitable boundary conditions using the bvp4c function of MATLAB. The convective instability analysis and the energy analysis are performed using the Chebyshev spectral method in order to obtain the neutral curves and the energy bars. We observe that the roughness of the disk has a destabilising effect on both Type-I and Type-II instability modes. The results obtained will be prominently treated as benchmarks for our future studies on stagnation flow.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.8
• /
• pp.764-772
• /
• 2009

In this study, fluid/structure coupled analyses have been conducted for 3-D stator and rotor configuration. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate fluid/structure responses of general stator-rotor configurations. To solve the fluid/structure coupled problems, fluid domains are modeled using the structural grid system with dynamic moving and local deforming techniques. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras(S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3-D turbine blades for fluid-structure interaction(FSI) problems. Detailed fluid/structure analysis responses for stator-rotor interaction flow conditions are presented to show the physical performance and flow characteristics.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.2
• /
• pp.150-159
• /
• 2010

Laidback fan shaped film-cooling hole is formulated numerically and optimized with the help of three-dimensional numerical analysis, surrogate methods, and the multi-objective evolutionary algorithm. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by four geometric design variables, the injection angle of the hole, the lateral expansion angle of the diffuser, the forward expansion angle of the hole, and the ratio of the length to the diameter of the hole, to maximize the film-cooling effectiveness compromising with the aerodynamic loss. The objective function values are numerically evaluated through Reynolds- averaged Navier-Stokes analysis at the designs that are selected through the Latin hypercube sampling method. Using these numerical simulation results, the Response Surface Approximation model are constructed for each objective function and a hybrid multi-objective evolutionary algorithm is applied to obtain the Pareto optimal front. The clustered points from Pareto optimal front were evaluated by flow analysis. These designs give enhanced objective function values in comparison with the experimental designs.

• International Journal of Fluid Machinery and Systems
• /
• v.8 no.3
• /
• pp.142-154
• /
• 2015

Nowadays, computational fluid dynamics is commonly used by design engineers to evaluate and compare losses in hydraulic components as it is less expensive and less time consuming than model tests. For that purpose, an automatic tool for casing and distributor analysis will be presented in this paper. An in-house mesh generator and a Reynolds Averaged Navier-Stokes equation solver using the standard $k-{\omega}$ shear stress transport (SST) turbulence model will be used to perform all computations. Two solvers based on the C++ OpenFOAM library will be used and compared to a commercial solver. The performance of the new fully coupled block solver developed by the University of Lucerne and Andritz will be compared to the standard 1.6ext segregated simpleFoam solver and to a commercial solver. In this study, relative comparisons of different geometries of casing and distributor will be performed. The present study is thus aimed at validating the block solver and the tool chain and providing design engineers with a faster and more reliable analysis tool that can be integrated into their design process.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.35-40
• /
• 2008

Generally flight vehicles have many cavities such as wheel wells, bomb bays and windows on their external surfaces and the flow around these cavities makes separation, vortex, shock and expansion waves, reattachment and other complex flow phenomenon. The flow around the cavity makes abnormal and three-dimensional noise and vibration even thought the aspect ratio (L/D) is small. The cavity giving large effects to the flow might make large noise, cause structural damage or breakage, harm the aerodynamic performance and stability, or damage the sensitive devices. In this study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's $\kappa-\omega$ turbulence model. The MPI(Message Passing Interface) parallelized code was used for calculations by PC-cluster. The cavity has the aspect ratios of 2.5, 3.5 and 4.5 with the W/D ratio of 2 for three-dimensional cavities. The Sound Pressure Level (SPL) analysis was done with FFT to check the dominant frequency of the cavity flow. The dominant frequencies were analyzed and compared with the results of Rossiter's formula and Ahuja& Mendoza's experimental datum.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.563-569
• /
• 2008

In this study, fluid/structure coupled analyses have been conducted f3r 3-D stator and rotor configuration. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate fluid/structure responses of general stator-rotor configurations. To solve the fluid/structure coupled problems, fluid domains are modeled using the structural grid system with dynamic moving and local deforming techniques. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3-D turbine blades for fluid-structure interaction (FSI) problems. Detailed fluid/structure analysis responses for stator-rotor interaction flow conditions are presented to show the physical performance and flow characteristics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.11
• /
• pp.1371-1378
• /
• 1999

Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.23 no.2
• /
• pp.51-56
• /
• 2015

This paper presents the procedure of drag prediction for EAV-1, based on a numerical analysis correlated to an in-flight test. EAV-1, developed by Korea Aerospace Research Institute, is a small-sized UAV to test a hydrogen-fuel cell power system. The long-endurance test flight of 4.5 hours provides numerous in-flight data. The thrust and drag of EAV-1 during the flight test are estimated based on the wind-tunnel test results for EAV-1's propeller performance. In addition, the CFD analysis using a commercial Navier-Stokes code is carried out for the full-scale EAV-1. The computational result suggests that the initial CFD analysis substantially under-predicts the in-flight drag in that the discrepancy is up to 27.6%. Therefore, additional investigation for more accurate drag prediction is performed; the effect of propeller slipstream is included in the CFD analysis through "fan disk" modelling. Also, the additional drag from airplane trim and load factor that actually exists during the flight test in a circular path is considered. These supplemental analyses for drag prediction turn out to be effective since the drag discrepancy reduces to 2.3%.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.5
• /
• pp.42-49
• /
• 2016

Aerodynamic analysis of a low-voltage electric motor has been performed with various inlet vent shapes. Effects of inlet vent shape on aerodynamic performance of a motor cooling fan have been investigated numerically using three-dimensional Reynolds-averaged Navier-Stokes equations. The k-${\varepsilon}$ turbulence model was used for the analysis of turbulence. The finite volume method and unstructured tetrahedral grids were used in the numerical analysis. Optimal grid system in the computational domain was selected through a grid-dependency test. From the results of the flow analysis, considerable energy loss by flow separation was observed in the flow passage. It was found that mass flow rate through the cooling fan in the low-voltage motor can be increased by modifying the inlet vent shape. And, some inlet vent shapes were suggested to improve the aerodynamic performance of the motor cooling fan.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.3
• /
• pp.229-236
• /
• 2016

Exergetic analysis was introduced in optimization of a rotating equilateral triangular internal cooling channel with staggered square ribs to maximize the net exergy gain. The objective function was defined as the net exergy gain considering the exergy gain by heat transfer and exergy losses by friction and heat transfer process. The flow field and heat transfer in the channel were analysed using three-dimensional Reynolds-averaged Navier-Stokes equations under the uniform temperature condition. Shear stress transport turbulence model has been selected as a turbulence closure through the turbulence model test. Computational results for the area-averaged Nusselt number were validated compared to the experimental data. Three design variables, i.e., the angle of rib, the rib pitch-to-hydraulic diameter ratio and the rib width-to-hydraulic diameter ratio, were selected for the optimization. The optimization was performed at Reynolds number, 20,000. Twenty-two design points were selected by Latin hypercube sampling, and the values of the objective function were evaluated by the RANS analysis at these points. Through optimization, the objective function value was improved by 22.6% compared to that of the reference geometry. Effects of the Reynolds number, rotation number, and buoyancy parameter on the heat transfer performance of the optimum design were also discussed.