• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1752-1756
• /
• 2008

A simple, accurate and efficient mesh generation technique, the cut-cell method, is able to represent an arbitrarily complex geometry. Both structured and unstructured grid meshes are used in this method. First, the numerical domain is constructed with regular Cartesian grids as a background grid and then the solid boundaries or bodies are cut out of the background Cartesian grids. As a result, some boundary cells can be contained two numerical conditions such as the flow and solid conditions, where the special treatment is needed to simulate such physical characteristics. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. Also, the TVD-WAF method is applied on the Cartesian cut-cell grids to stabilize numerical results. Present method is validated for the rectangular dam break problems. Initially, a conventional grid is constructed with the Cartesian regular mesh only and then applied to the dam-break flow simulation. As a comparative simulation, a cut-cell grids are applied to represent the flow domain rotated with arbitrary angles. Numerical results from this study are compared with the results from the case of the Cartesian regular mesh only. A good agreement is achieved with other numerical results presented in the literature.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.05a
• /
• pp.275-280
• /
• 2008

Shock tube flow measurement has been often troubled with 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 specify 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.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.161-164
• /
• 2002

In the present study, a passive control method, using the 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. Baldwind-Lomax turbulence model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure losses 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/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 influences the control effectiveness.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.151-154
• /
• 2002

The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee's total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number $M_{s}$, of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin-impulse waves and its location strongly depend upon the distance between two parallel tubes, L/d and the incident shock Mach number, $M_{s}$. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.5
• /
• pp.666-674
• /
• 2002

A rapid expansion of the moist air or stream through transonic nozzle often leads to not-equilibrium condensation shock, causing a considerable amount of energy loss to the entire flow field. Depending on amount of heat released, condensation shock wave occurs in the nozzle and interacts with the boundary layer flow. In the current study, a passive control technique using a porous wall with a plenum cavity underneath is applied for purpose of alleviation the condensation shock wave in a transonic nozzle. A droplet growth equation is incorporated into two-dimensional wavier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. An experiment using an indraft transonic wind tunnel is made to validate the present computational results. The results obtained show that the magnitude of condensation shock wave is reduced by the current passive control method.

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

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.432-437
• /
• 2001

A computational work of the impulse wave which is discharged from the open end of a pipe is compared to the Lighthill's aeroacoustics theory. The second-order total variation diminishing(TVD) scheme is employed to solve the axisymmetric, compressible, unsteady Euler equations. The relationship between the initial compressure wave form and the resulting impulse wave is characterized in terms of the peak pressure. The overpressure, pressure gradient and wavelength of the initial compression wave are changed to investigate the influence of the initial compressure wave form on the peak pressure of impulse wave. The results obtained show that for the initial compression wave of a large wavelength and small pressure gradient the peak pressure of the impulse wave depends upon the wavelength and pressure gradient of compression wave, but for the initial compression wave of a short wavelength and large pressure gradient the peak pressure of the impulse wave is almost constant regardless of the wavelength and pressure gradient of compression wave. The peak pressure of the impulse wave is increased with an increase in the overpressure of the initial compression wave. The results from the numerical analysis are well compared to the results from the aeroacoutics theory with a good agreement.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.640-645
• /
• 2001

The present study is to investigate the characteristics of the impulse waves discharged from the exit of the convergent and divergent pipes. An experiment is carried out using a shock tube with an open end and is compared to the computation of the axisymmetric, compressible, unsteady Euler equations, which are solved by the second-order total variation diminishing(TVD) scheme. For the computational work, some initial compression waves are assumed inside the pipe so that those are identical to the experimental ones of the shock tube. The results show that the peak pressures of the impulse waves discharged from the exit of convergent and divergent pipes decrease with an increase in the wavelength of the initial compression wave. All of the impulse waves have a strong directivity toward the pipe axis, regardless of the exit type of the pipe employed. The impulse waves discharged from the divergent pipe are stronger than those from the straight pipe, while the impulse waves of the convergent pipe are weaker than those from the straight pipe. It is believed that the convergent pipe can playa role of a passive control to reduce the peak pressure of the impulse wave. The present computations represent the experimented impulse waves with a good accuracy.