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

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.245-256
• /
• 2014

This work is focused on the boundary layer and heat transfer characteristics of hydromagnetic flow over a stretching sheet with variable thickness. Steady, two dimensional, nonlinear, laminar flow of an incompressible, viscous and electrically conducting fluid over a stretching sheet with variable thickness and power law velocity in the presence of variable magnetic field and variable temperature is considered. Governing equations of the problem are converted into ordinary differential equations utilizing similarity transformations. The resulting non-linear differential equations are solved numerically by utilizing Nachtsheim-Swigert shooting iterative scheme for satisfaction of asymptotic boundary conditions along with fourth order Runge-Kutta integration method. Numerical computations are carried out for various values of the physical parameters and the effects over the velocity and temperature are analyzed. Numerical values of dimensionless skin friction coefficient and non-dimensional rate of heat transfer are also obtained.

• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.15-19
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• Journal of the Korean Society of Combustion
• /
• v.13 no.1
• /
• pp.17-22
• /
• 2008

A partly implicit/quasi-explicit method is introduced for the solution of detailed chemical kinetics with stiff source terms based on the standard fourth-order Runge-Kutta scheme. Present method solves implicitly only the stiff reaction rate equations, whereas the others explicitly. The stiff equations are selected based on the survey of the chemical Jaconian matrix and its Eigenvalues. As an application of the present method constant pressure combustion was analyzed by a detailed mechanism of hydrogen-air combustion with NOx chemistry. The sensitivity analysis reveals that only the 4 species in NOx chemistry has strong stiffness and should be solved implicitly among the 13 species. The implicit solution of the 4 species successfully predicts the entire process with same accuracy and efficiency at half the price.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.10a
• /
• pp.95-98
• /
• 2003

The physics of the plume-induced shock and separation particulary at a high plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with aid without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG k-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The courol methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated. The computational results show the main effect of the porous extension on plume-afterbody interactions is to in the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect on plume interference.

• The KSFM Journal of Fluid Machinery
• /
• v.3 no.3 s.8
• /
• pp.7-11
• /
• 2000

A numerical study based on the three-dimensional thin-layer Navier-Stokes solver is carried out to analyze the flowfield through a single stage transonic compressor. Explicit fout-step Runge-Kutta scheme with spatially variable time step and implicit residual smoothing is used. The governing equations we discretized with explcit finite difference method. Mired-out average method is used at the interface between rotor and stator. And, an artificial dissipation model is used to assure the stability of solution. The results with k-w turbulence model were compared to the results with Baldwin-Lomax model, and physical phenomena of transonic compressor are presented. The two turbulence models give the results that show reasonably good agreements with experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.160-166
• /
• 1998

In this paper, incompressible two-dimensional Navier-Stokes equations are numerically solved for the study of steady laminar flow around a body with the wall effect. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. To investigate the wall effect, numerical computations are carried out for the NACA 0012 section at the various blockage ratios. The pressure and skin friction on the foil surface, velocity pronto in its wake and drag coefficient are investigated as functions of the blockage ratio.

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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

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

This paper presents the result of a computational study on the wake characteristics of two tanker models. i.e HSVA and DYNE hull forms. The focus of the study is on the distributions of axial. radial and tangential velocities of the two hull forms in way of the propeller, especially over the propeller disk. The effect of bilge vortices on the velocity distribution is also concerned. For the computation of stern and wake flows of the two hull forms. the incompressible Reynolds-Averaged Navier-Stokes(RANS) equations are numerically solved by the use of a second order finite difference method, which employs a four stage Runge-Kutta scheme with a residual averaging technique and the Baldwin-Lomax model. The calculated pressure distributions on the hull surface and the axial. radial and tangential velocity distributions over the propeller disk are presented for the two hull forms. Finally, the result of wake analysis for the computed wake distribution over the propeller disk is given in comparison with those for the experimental wake distribution for the both hull forms.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2060-2065
• /
• 2003

Aeolian tone generation from a two dimensional circular cylinder is numerically investigated via direct numerical simulation and hydrodynamic-acoustic splitting method. All governing equation are spatially discretized with the sixth-order compact scheme and fourth-order Runge-Kutta method to avoid excessive numerical dissipations and dispersions of acoustic quantities. Comparisons of two results show that the previous splitting method can not accurately predict the aeroacoustic noise of wall bounded shear flow. In this study, a perturbation viscous term and a new energy equation have been developed. This modified splitting method accurately predicts aeroacoustic noise from wall-bounded shear flow. The present results agree very well with the direct numerical simulation solution.