• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.86-87
• /
• 2003

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be changed and deformed with time if we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.24-30
• /
• 2005

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be time-dependently changing and deforming according to the movement of the boundaries when we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.337-345
• /
• 2002

The present study addresses a computational result of unsteady gas flow through a critical nozzle. The axisymmetric, unsteady, compressible, Wavier-Stokes equations are solved using a finite volume method that makes use of the second order upwind scheme for spatial derivatives and the multi-stage Runge-Kutta integral scheme for time derivatives. The steady solutions of the governing equation system are validated with the previous experimental data to ensure that the present computational method is valid to predict the critical nozzle flows. In order to simulate the effects of back pressure fluctuations on the critical nozzle flows, an excited pressure oscillation with an amplitude and frequency is assumed downstream of the exit of the critical nozzle. The results obtained show that for low Reynolds numbers, the unsteady effects of the pressure fluctuations can propagate upstream of the throat of critical nozzle, and thus giving rise to the applicable fluctuations in mass flow rate through the critical nozzle, while for high Reynolds numbers, the pressure signals occurring at the exit of the critical nozzle do not propagate upstream beyond the nozzle throat. For very low Reynolds number, it is found that the sonic line near the throat of the critical nozzle remarkably fluctuateswith time, providing an important mechanism for pressure signals to propagate upstream of the nozzle throat, even in choked flow conditions. The present study is the first investigation to clarify the unsteady effects on the critical nozzle flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.145-150
• /
• 2005

Generally, Compressible Navier-Stokes codes are used to solve high mach number flows. But, Most of high mach number flows embrace low mach number flows. This phenomenon results in low convergence rate and non-physical solution in CFD analysis. So Many researchers developed preconditioning technique to solve these problems. This Study presents how to modify previous compressible N-S computer code with little changes of structure into preconditioned compressible N-S code applying Roe's Approximate Riemann Solver. And this study show developed preconditioning code is very well operated at all mach number flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.6
• /
• pp.788-802
• /
• 1998

A numerical method using dual time stepping and preconditioning procedure for efficient computations of unsteady low speed compressible flow problems is developed. The time-derivative preconditioning method which is valid at low speed flow conditions cannot maintain temporal accuracy because of the modification of the time-derivative term in Navier-Stokes equations. The dual time stepping procedure is incorporated to enable the time accurate computations and this procedure introduces a pseudo-time derivative in addition to the physical time derivative. At a given physical time, an inner iteration can be carried out until a steady state in pseudo-time is achieved. This will effectively yield a time accurate solution. Computational capabilities of the above algorithm are demonstrated through computation of a variety of practical fluid flows and it is shown that the algorithms is efficient in the essentially incompressible flows and low Mach number compressible flows with heat source.

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

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.556-561
• /
• 2003

Flight bodies are subject to highly unstable and severe flow conditions during taking-off and landing periods. In this situation, the flight bodies essentially experience accelerating or decelerating flows, and the aerodynamic characteristics can be completely different from those of steady flows. In the present study, unsteady aerodynamic characteristics of an aerofoil accelerating at subsonic speeds are investigated using a computational method. Two-dimensional, unsteady, compressible Navier-Stokes simulations are conducted with a one-equation turbulence model, Spalart-Allmaras, and a fully implicit finite volume scheme. An acceleration factor is defined to specify the unsteady aerodynamics of the aerofoil. The results show that the acceleration of the subsonic aerofoil generally leads to a variation in aerodynamic characteristics and it is more significant at angles of attack.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.60-68
• /
• 2005

For study on the unsteady wall interference effect, flows around a forced oscillating airfoil in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The Spalart-Allmaras one-equation model is employed for the turbulence effect. The computed results of the oscillating airfoil having a thin wake showed that the lift curve slope is increased and the magnitude of hysteresis loop is reduced by the interference effects. Since the vortex around the airfoil is generated and convected downstream faster than the free-air condition, the phase of lift, drag and pitching moment coefficients was shifted. The pressure on the test section wall shows harmonic terms having the oscillating frequency contained in the wail effect.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1511-1521
• /
• 2002

High-speed moist air or steam flow has long been of important subject in engineering and industrial applications. Of many complicated gas dynamics problems involved in moist air flows, the most challenging task is to understand the nonequilibrium condensation phenomenon when the moist air rapidly expands through a flow device. Many theoretical and experimental studies using supersonic wind tunnels have devoted to the understanding of the nonequilibrium condensation flow physics so far. However, the nonequilibrium condensation can be also generated in the subsonic flows induced by the unsteady expansion waves in shock tube. The major flow physics of the nonequilibrium condensation in this application may be different from those obtained in the supersonic wind tunnels. In the current study, the nonequilibrium condensation phenomenon caused by the unsteady expansion waves in a shock tube is analyzed by using the two-dimensional, unsteady, Navier-Stokes equations, which are fully coupled with a droplet growth equation. The third-order TVD MUSCL scheme is applied to solve the governing equation systems. The computational results are compared with the previous experimental data. The time-dependent behavior of nonequilibrium condensation of moist air in shock tube is investigated in details. The results show that the major characteristics of the nonequilibrium condensation phenomenon in shock tube are very different from those in the supersonic wind tunnels.