• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.12
• /
• pp.785-792
• /
• 2016

In the present study, the effects of non-equilibrium condensation on the drag divergence Mach number with the angle of attack in a transonic 2D moist air flow of NACA0012 are investigated using the TVD finite difference scheme. For the same ${\alpha}$, the maximum upstream Mach number of the shock wave, Mmax, and the size of supersonic bubble decrease with the increase in ${\Phi}_0$. For the same $M_{\infty}$, ${\Phi}_0$, and $T_0$, the length of the non-equilibrium condensation zone ${\Delta}_z$ decreases with increasing ${\Phi}_0$. On the other hand, because of the attenuating effect of non-equilibrium condensation on wave drag, which is related to the interaction between the shock wave and the boundary layer, the drag coefficient $C_D$ decreases with an increase in ${\Phi}_0$ for the same $M_{\infty}$ and ${\alpha}$. For the same ${\alpha}$, $M_D$ increases with increasing ${\Phi}_0$, while $M_D$ decreases with an increase in ${\alpha}$.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.12
• /
• pp.1009-1015
• /
• 2014

The present study investigated the effects of non-equilibrium condensation with the angle of attack on the coefficients of pressure, lift, and drag in the transonic 2-D flow of NACA0012 by numerical analysis of the total variation diminishing (TVD) scheme. At $T_0=298k$ and ${\alpha}=3^{\circ}$, the lift coefficients for $M_{\infty}=0.78$ and 0.81 decreased monotonically with increasing ${\Phi}_0$. In contrast, for $M_{\infty}$ corresponding to the Mach number of the force break, $C_L$ increased with ${\Phi}_0$. For ${\alpha}=3^{\circ}$ and ${\Phi}_0=0%$, $C_D$ increased markedly as $M_{\infty}$ increased. However, at ${\Phi}_0=60%$ and ${\alpha}=3^{\circ}$, which corresponded to the case of the condensation having a large influence, $C_D$ increased slightly as $M_{\infty}$ increased. The decrease in profile drag by non-equilibrium condensation grew as the angle of attack and stagnation relative humidity increased for the same free stream transonic Mach number. At ${\Phi}_0=0%$, the coefficient of the wave drag increased with the attack angle and free stream Mach number. When ${\Phi}_0$ > 50%, the coefficient of the wave drag decreased as ${\alpha}$ and $M_{\infty}$ increased. Lowering ${\Phi}_0$ and increasing $M_{\infty}$ increased the maximum Mach number.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.3
• /
• pp.248-253
• /
• 2009

This paper presents a numerical study for unsteady flows in a high-pressure steam turbine with a partial admission stage. Compressible Navier-Stokes equations are solved by the high-order high-resolution finite-difference method based on the fourth-order compact MUSCL TVD scheme, Roe's approximate Riemann solver, and the LU-SGS scheme. The SST-model is also solved for evaluating the eddy-viscosity. The unsteady two-dimensional flows through whole nozzle-rotor cascade channels considering a partial admission are numerically investigated. 108 nozzle passages with two blockages and 60 rotor passages are simultaneously calculated. The influence of the flange in the nozzle box to the lift of rotors is predicted. Also the efficiency of the partial admission stage changing the number of blockages and the number of nozzles is parametrically predicted.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1079-1089
• /
• 2009

The objective of this study is to develop the accurate, robust and high resolution two-dimensional numerical model that solves the computationally difficult hydraulic problems, including the wave front propagation over dry bed and abrupt change in bathymetry. The developed model in this study solves the conservative form of the two-dimensional shallow water equations using an unsplit finite volume scheme and HLLC approximate Riemann solvers to compute the interface fluxes. Bed-slope term is discretized by the divergence theorem in the framework of FVM for application of unsplit scheme. Accurate and stable SGM, in conjunction with the MUSCL which is second-order-accurate both in space and time, is adopted to balance with fluxes and source terms. The exact C-property is shown to be satisfied for balancing the fluxes and the source terms. Since the spurious oscillations in second-order schemes are inherent, an efficient slope limiting technique is used to supply TVD property. The accuracy, conservation property and application of developed model are verified by comparing numerical solution with analytical solution and experimental data through the simulations of one-dimensional dam break flow without bed slope, steady transcritical flow over a hump and two-dimensional dam break flow with a constriction.

• 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 of Propulsion Engineers
• /
• v.13 no.3
• /
• pp.27-33
• /
• 2009

A computational study has been performed out to evaluate the effect of a vertical porous barrier on the pressure oscillations in a supersonic cavity. The porous barriers with different perforations were vertically installed into a rectangular cavity at Mach numbers 1.50, 1.83 and 2.50. TVD finite difference MUSCL scheme was employed to solve the two-dimensional, unsteady, compressible Navier-Stokes equations. The present vertical porous barrier considerably altered the characteristics of the time-dependent shear layers that occur at the upstream edge of cavity and remarkably reduced the pressure oscillations inside the supersonic cavity. The present results showed that the effectiveness of passive control using the present porous vertical barrier is dependent on Mach number and the perforation of the porous barrier.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.1
• /
• pp.27-33
• /
• 2009

Shock tube flow measurement has been often hampered a finite opening time of diaphragm, but there is no systematic work to investigate its effect on the shock tube flow. In the present study, both the experimental and computational works have been performed on the shock tube flows at low pressure ratios. The computational analysis has been performed using the two-dimensional, unsteady, compressible Navier-Stokes equations, based upon a TVD MUSCL finite difference scheme. It is known that the present computational results reproduce the experimental data with good accuracy and simulate successfully the process of diaphragm opening as a function of time. The concept of an imaginary center is introduced to quantify the non-centered expansion wave due to a finite opening time of diaphragm. The results obtained show that the diaphragm opening time is reduced as the initial pressure ratio of shock tube increases, leading to the effect of a finite opening time of diaphragm to be more remarkable at low pressure ratios.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.21-29
• /
• 2002

In the present study, a passive control method, using a porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible, Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwin-Lomax model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure loss of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock wave/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil as well. It is also found that the location of the porous ventilation significantly affects the control effectiveness.