• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.47-51
• /
• 2004

Accurate Prediction of a supersonic missile base drag continues to defy even well-rounded CFD codes. In an effort to address the accuracy and predictability of the base drags, the influence of grid system and competitive turbulence models on the base drag is analyzed. Characteristics of some turbulence models is reviewed through incompressible turbulent flow over a flat plate, and performance for the base drag prediction of several turbulence models such as Baldwin-Lomax(B-L), Spalart-Allmaras(S-A), $\kappa-\epsilon$, $\kappa-\omega$ model is assessed. When compressibility correction is injected into the S-A model, prediction accuracy of the base drag is enhanced. The NSWC wind tunnel test data are utilized for comparison of CFD and semi-empirical codes on the accuracy of base drag predictability: they are about equal, but CFD tends to perform better. It is also found that, as angle of attack of a missile with control (ins increases, even the best CFD analysis tool we have lacks the accuracy needed for the base drag prediction.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.299-304
• /
• 2001

The current study describes experimental and computational work on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order to control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is incorporated into the two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order tractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indraft transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is significantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.55-63
• /
• 1995

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.450-455
• /
• 2001

The time-dependent behavior of nonequilibrium condensation of moist air through the Ludwieg tube is investigated with a computational fluid dynamics(CFD) method. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme, with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The computational results are compared with the previous experiments using the Ludwieg tube with a downstream diaphragm. The results clearly show that for an initial relative humidity below 30% there is no periodic oscillation of the condensation shock wave, but for an initial relative humidity over 40% the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to nonequilibrium condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity, and the periodic excursions of the condensation shock wave are responsible for the total pressure loss.

• Journal of Mechanical Science and Technology
• /
• v.17 no.12
• /
• pp.2066-2077
• /
• 2003

The time-dependent behavior of unsteady condensation of moist air through the Ludwieg tube is investigated by using a computational fluid dynamics (CFD) work. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme, with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The predicted results are compared with the previous experiments using the Ludwieg tube with a diaphragm downstream. The present computations represent the experimental flows well. The time-dependent unsteady condensation characteristics are discussed based upon the present predicted results. The results obtained clearly show that for an initial relative humidity below 30% there is no periodic oscillation of the condensation shock wave, but for an initial relative humidity over 40% the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to unsteady condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity and it results from the periodic excursions of the condensation shock wave.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.997-1004
• /
• 2002

The current study describes experimental and computational works on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order 10 control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is copuled with two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indrafi transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is signi ficantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

• Advances in aircraft and spacecraft science
• /
• v.7 no.4
• /
• pp.371-385
• /
• 2020

A forward-facing aerospike attached to a payload fairing of a satellite launch vehicle significantly alters its flowfield and decreases the aerodynamic drag in transonic and low supersonic speeds. The present payload fairing is an axisymmetric configuration and consists of a blunt-nosed body along with a conical section, payload shroud, boat tail and followed by a booster. The main purpose of the present numerical simulations is to evaluate flowfield and assess the performance of aerodynamic drag coefficient with and without aerospike attached to a payload fairing of a typical satellite launch vehicle in freestream Mach number range 0.8 ≤ M_∞ ≤ 3.0 and freestream Reynolds number range 33.35 × 10⁶/m ≤ Re_∞ ≤ 46.75 × 10⁶/m whichincludes the maximum aerodynamic drag and maximum dynamic conditions during ascent flight trajectory of the satellite launch vehicle. A numerical simulation has been carried out to solve time-dependent compressible turbulent axisymmetric Reynolds-averaged Navier-Stokes equations. The closure of the system of equations is achieved using the Baldwin-Lomax turbulence model. The aerodynamic drag reduction mechanism is analysed employing numerical results such as velocity vector plots, density and Mach contours in conjunction with the experimental flow visualization pictures. The variations of wall pressure coefficient over the payload fairing with and without aerospike are exhibiting different kind of flowfield characteristics in the transonic and low supersonic speeds. The numerically computed results are compared with schlieren pictures, oil flow patterns and measured wall pressure distributions and exhibit good agreement between them.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.123-132
• /
• 2002

A numerical study was conducted to investigate the combustion phenomena of normal start and unstart processes based on ISL's RAMAC 30 experiments with different diluent amounts and fill pressures in a ram accelerator. The initial projectile launching speed was 1.8 km/s which corresponded to the superdetonative speed of the stoichiometric $H_2/O_2$ mixture diluted with 5 $CO_2$ or 4 $CO_2$. Experiments with same condition except for projectile surface material demonstrated that ignition was successful with an aluminum projectile, but no combustion was observed in case of a steel projectile. In this study, it was found that neither shock nor viscous heating was sufficient to ignite the mixture at a low speed of 1.8 km/s, as was found in the experiments using a steel projectile. However, we could succeed in igniting the mixtures by imposing a minimal amount of additional heat to the combustor section and simulate the normal start and unstart processes found in the experiments with an aluminum projectile. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations coupled with a Baldwin-Lomax turbulence model and detailed chemistry reaction equations of $H_2/O_2/CO_2$ suitable for high-pressure gaseous combustion were considered. The governing equations were discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit, time accurate integration method. The numerical results matched almost exactly to the experimental results. As a result, it was found that the normal start and unstart processes depended on the strength of gas mixture, development of shock-induced combustion wave stabilized by the first separation bubble, and its size and location.