• 대한기계학회논문집B
• /
• 제22권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.557-560
• /
• 2002

In this study, simulation of weak shock waves are peformed by a two-dimensional thermal fluid or compressible fluid model of the lattice Boltzmann method. The shock wave represents an abrupt change in fluids properties, in which finite variations in pressure, internal energies, and density occur over the shock thickness. The characteristics of the proposed model with a simple distribution function is verified by calculation of the sound speeds, and the shock tube problem. The reflection of a weak shock wave by wedge propagating in a channel is performed. The results agree well with those by finite difference method or by experiment. In the simulation of unsteady shock wave diffraction around a sharp corner, we show a flow field of vortical structure near the comer.

• Journal of Mechanical Science and Technology
• /
• 제20권2호
• /
• pp.286-294
• /
• 2006

Flows in a ramjet inlet is simulated for the study of the rocket-ramjet transition. The flow is unsteady, two-dimensional axisymmetric, compressible and turbulent. Double time marching method is used for the unsteady calculation and HLLC method is used as a higher order MUSCL method. As for turbulent calculation, $\kappa-\omega$ SST model is used for more accurate viscous calculations. Sinusoidal pressure perturbation is given at the exit and the flow fields at the inlet is studied. The cruise condition as well as the ground test condition are considered. The pressure level for the ground test condition is relatively low and the effect of the pressure perturbation at the combustion chamber is small. The normal shock at the cruise condition is very sensitive to the pressure perturbation and can be easily detached from the cowl when the exit pressure is relatively high. The sudden decrease in the mass flux is observed when the inlet flow becomes subcritical, which can make the inlet incapable. The amplitude of travelling pressure waves becomes larger as the downstream pressure increases, and the wavelength becomes shorter as Mach number increases. The phase difference of the travelling perturbed pressure wave in space is 180 degree.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.414-419
• /
• 2011

Missile am fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 60 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure distribution, etc. are analyzed according to the angle of attack. The results at a low angle of attack are compared with other results before a stall condition. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. Unsteady velocity field, periodic vortex shedding, the unsteady pressure distribution on the airfoil surface, and the acoustic fields are analyzed. The effects of these unsteady characteristics in the aerodynamic coefficients are analyzed.

• 한국전산유체공학회지
• /
• 제18권1호
• /
• pp.1-6
• /
• 2013

Missile and fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 50 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure, entropy distribution, etc. are analyzed according to the angle of attack. The results of average lift coefficients are compared with other results according to the angle of attack. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. The primary and secondary oscillating frequencies are analyzed by the effects of these unsteady aerodynamic characteristics.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1742-1747
• /
• 2003

When a shock wave is discharged from the exit of a duct, complicated flow is formed near the duct exit. The flow field is much more complicated under the ground effects or any other objects near the exit of a duct, such as the circumstance near the exit of the high-speed railway tunnel. The resulting flow is essentially three-dimensional unsteady with the effects of strong compressibility. In the current study, three-dimensional flow fields of the weak shock wave which is discharged from the exit of a duct are numerically investigated using a CFD method. Computations are performed for the weak shock wave in the range below 1.5. The results obtained show that the directivity and magnitude of the weak shock discharged strongly depend upon the Mach number of initial shock wave and are significantly influenced by the ground effects.

• 대한기계학회논문집B
• /
• 제21권12호
• /
• pp.1624-1634
• /
• 1997

A numerical study on two-dimensional unsteady transonic cascade flow has been performed by adopting dual time stepping and the k-.omega. turbulence model. An explicit 4 stage Runge-Kutta scheme for the compressible Navier-Stokes equations and an implicit Gauss-Seidel iteration scheme for the k-.omega. turbulence model are proposed for fictitious time stepping. This mixed time stepping scheme ensures the stability of numerical computation and exhibits a good convergence property with less computation time. Typical steady-state convergence accelerating schemes such as local time stepping, residual smoothing and multigrid combined with dual time stepping shows good convergence properties. Numerical results are presented for unsteady laminar flow past a cylinder and turbulent shock buffeting problem for bicircular arc cascade flow is discussed.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
• /
• pp.387-390
• /
• 2006

Numerical study has been conducted in two dimensions about a NACA0012 airfoil with an oscillating microflap on the surface. We show that this microflap is effective in controlling the unsteady stall at high angles of attack. We solve the compressible Navier-Stokes equations for the Reynolds numbers with an extensible chimera grid fitted to the oscillatory microflap. For turbulent calculation, we adopt the SST $k-{\omega}$ model. We investigate the parametric effect of angle of attacks, Reynolds number, and the location where the microflap is installed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.297-301
• /
• 2008

The flight vehicles have cavities such as wheel wells and bomb bays. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves. Resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. The flow field is observed to oscillate in the "shear layer mode" with low aspect ratio. In the present 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 flow field is observed to oscillate in the shear layer mode" with large aspect ratio. Based on the SPL(Sound Pressure Level) analysis of the pressure variation at the cavity trailing edge, the dominant frequency was analyzed and compared with the results of Rossiter's formul. The aero-acoustic wave analyzed with CPD(Correlation of Pressure Distribution).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.297-301
• /
• 2008

The flight vehicles have cavities such as wheel wells and bomb bays. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves. Resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. The flow field is observed to oscillate in the "shear layer mode" with low aspect ratio. In the present study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's ${\kappa}-{\varepsilon}$ turbulence model. The flow field is observed to oscillate in the "shear layer mode" with large aspect ratio. Based on the SPL(Sound Pressure Level) analysis of the pressure variation at the cavity trailing edge, the dominant frequency was analyzed and compared with the results of Rossiter's formul. The aero-acoustic wave analyzed with CPD(Correlation of Pressure Distribution).