• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.4
• /
• pp.291-306
• /
• 2011

We consider the second order Strang splitting method to approximate the solution to the KdV equation. The model equation is split into three sets of initial value problems containing convection and dispersal terms separately. TVD MUSCL or MUSCL scheme is applied to approximate the convection term and the second order centered difference method to approximate the dispersal term. In time stepping, explicit third order Runge-Kutta method is used to the equation containing convection term and implicit Crank-Nicolson method to the equation containing dispersal term to reduce the CFL restriction. Several numerical examples of weakly and strongly dispersive problems, which produce solitons or dispersive shock waves, or may show instabilities of the solution, are presented.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.343-346
• /
• 2005

Three dimensional structures of detonation wave propagating through a square-shaped duct were investigated 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 MUSCL-type TVD scheme and four stage Runge-Kutta time integration. The unsteady computational results in three dimension show the detailed mechanism of transverse mode and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record.

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

The purpose of this study is analysis of flow field where is around of injector of supersonic combustor which is bluff-body stabilized flame and hyper-mixer type of supersonic combustor injector by using hydrogen or hydrocarbon fuel. Various schemes are evaluated to supersonic backward step flow filed with massive separation region in validation step. Compounded scheme of 5th-order TVD-MUSCL, Roe FDS, S-A DES/DDES has a good performance in base and base-bleed flow.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.41-44
• /
• 2003

In spite of many researches made on the supersonic jets until now, detailed three-dimensional structures of supersonic jets are not well hewn. In the current study, the detailed structures of three-dimensional supersonic jets are numerically investigated using a CFD method. The total variation diminishing (TVD) scheme is used to solve the unsteady, three-dimensional, compressible Euler equations. Computational results are visualized to investigate the major features of supersonic jets. The three-dimensional computation results show that the structures of the supersonic jets are significantly different from those of the two-dimensional or axisymmetric supersonic jets.

• Journal of Korea Water Resources Association
• /
• v.44 no.2
• /
• pp.145-156
• /
• 2011

Two dimensional numerical model of high-order accuracy is developed to analyze complex flow including transition flow, discontinuous flow, and wave propagation to dry bed emerging at natural river flow. The bed slope term of two dimensional shallow water equation consisting of integral conservation law is treated efficiently by applying quasi-steady wave propagation scheme. In order to apply Finite Volume Method using Fractional Step Method, MUSCL scheme is applied based on HLL Riemann solver, which is second-order accurate in time and space. The TVD method is applied to prevent numerical oscillations in the second-order accurate scheme. The developed model is verified by comparing observed data of two dimenstional levee breach experiment and dam breach experiment containing structure at lower section of channel. Also effect of the source term is verified by applying to dam breach experiment considering the adverse slope channel.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1294-1299
• /
• 2004

Transient radiative heat transfer is analyzed in a one-dimensional slab using finite volume method (FVM). In this study, the step, $2^{nd}$ order upwind, and QUICK schemes are used for incident diffuse radiation and collimated beam, respectively. The results for diffuse radiation show that all schemes applied in this study give good agreements with available published results. In case of collimated beam, however, the results show deviations from the analytical solutions. To successfully describe the propagations of collimated beam, shock capturing schemes such as TVD scheme are need to be developed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.3 s.246
• /
• pp.201-207
• /
• 2006

Transient radiative heat transfer is analyzed in a one-dimensional slab using finite volume method (FVM). In this study, the step, $2^{nd}$ order upwind, and QUICK schemes are used for incident diffuse radiation and collimated beam, respectively. The results fer diffuse radiation show that all schemes applied in this study give good agreements with available published results. In case of collimated beam however, the results show deviations from the analytical solutions. To successfully describe the propagations of collimated beam shock capturing schemes such as TVD scheme are need to be developed.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.69-78
• /
• 2008

The high-order numerical method based on the adaptive mesh refinement(AMR) on the quadrilateral unstructured grids has been developed in this paper. This adaptive-grid method, originally developed with MUSCL-TVD scheme, is now extended to the WENO (weighted essentially no-oscillatory) scheme with the Runge-Kutta time integration of fifth order in spatial and temporal accuracy. The multidimensional interpolation was studied in the preliminary research, which allows us to maintain the same order of accuracy for the computation of numerical flux between two adjacent cells of different levels. Some standard benchmark tests are done to validate this method for checking the overall capacity and efficiency of the present adaptive-grid technique.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1996.11a
• /
• pp.107-116
• /
• 1996

Numerical simulations have been conducted to characterize unsteady flow structures in an axisymmetric supersonic inlet diffuser with sinusoidal pressure oscillations at the diffuser exit. The formulation is based on the unsteady Navier-Stokes equations and turbulence closure is achieved using a two-layer model with a too-Reynolds-number scheme for the near-wall treatment. The governing equations are formulated in an integral form, and are discretized by the four-stage Runge-Kutta scheme for temporal terms and the Harten-Yee upwind TVD scheme for convective terms. Results indicated that the inlet shock characteristics are significantly modified by acoustic oscillations originating from the combustor. The characteristics of shock/boundarv-layer interactions (such as the size of separation bubble, terminal shock shape, and vorticity intensity) are also greatly iufluenced by the shock oscillation due to acoustic waves.

• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.