• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

• Journal of the Korean Society of Combustion
• /
• v.1 no.2
• /
• pp.41-49
• /
• 1996

A numerical study is conducted to investigate the periodically unstable shock induced combustion around blunt bodies in stoichiometric hydrogen-air mixtures. Euler equations are spatially discretized by upwind-biased third order scheme and temporally integrated by Runge-Kutta method. Chemistry model used in this study involves 8 elementary kinetics steps and 7 species. At a constant Mach number, the effects of projectile size, inflow pressure and inflow temperature are examined with Lehr#s experimental condition as a reference. In addition to oscillation frequency, characteristic distances and time averaged values are found from the result to find an relation with dimensionless parameters. As a result, it is found that the effects of inflow pressure and body size are very similar and $Damk{\ddot{o}}hler$ number plays an important role in determining the instability characteristics.

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

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

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

• Water for future
• /
• v.17 no.1
• /
• pp.45-56
• /
• 1984

The whole process from the model design to the results of the test, of hydraulic model test of restricted entry surge tank of Hapcheon dam, is reviewed with the respect to the flowchart of the experiment. And the experimental results are compared with the numerical values which are calculated by Runge-Kutta-Gill scheme. The comparision show a reasonable agreement. In final design, it doesn't matter that only numerical values are considered in case of the short design period, or difficulties of budget, and or the comparably simple type surge tank as Hapcheon dam.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.4
• /
• pp.350-356
• /
• 1993

We carry out a numerical calculation to understand the nonlinear wave deformation around breakwaters using the Boussinesq equation, which is weakly nonlinear and weakly dispersive shallow water equation. A numerical method based on a finite element scheme and fourth order Runge-Kutta algorithm is employed to investigate the diffraction of incident waves by the breakwater. As a computational model, two-dimensional wave flume is treated. The breakwaters is perpendicular to the side wall of a channel. From the numerical results, the wave deformations according to the change of the length and the thickness of breakwaters are investigated. We also investigate the effect of the nonlinearity by comparing the results with the linear solutions.

• Bulletin of the Korean Mathematical Society
• /
• v.48 no.4
• /
• pp.873-896
• /
• 2011

The semi-discrete central scheme and central upwind scheme use Runge-Kutta (RK) time discretization. We do the Lax-Wendroff (LW) type time discretization for both schemes. We perform numerical experiments for various problems including two dimensional Riemann problems for Burgers' equation and Euler equations. The results show that the LW time discretization is more efficient in CPU time than the RK time discretization while maintaining the same order of accuracy.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.310-313
• /
• 2006

The present study aims at investigating the effectiveness of a new passive cavity flow control technique, sub-cavity. The characteristics of cavity flow oscillation with the device are compared with those with other control techniques tested previously, including a triangular bump and blowing jet. In the computation, the three-dimensional, unsteady Navier-Stokes equations governing the supersonic cavity flow are solved based on an implicit finite volume scheme spatially and multi-stage Runge-Kutta scheme temporally. Large eddy simulation (LES) is carried out to properly predict the turbulent features of cavity flow. The present results show that the pressure oscillation near the downstream edge dominates overall time-dependent cavity pressure variations, and the amplitude of the pressure oscillation can be reduced in the presence of a sub-cavity.

• Water for future
• /
• v.29 no.5
• /
• pp.173-184
• /
• 1996

In this study, a two-dimensional shallow-water equation was used to develop the mathematical model for computing water levels and flow distribution. In the discretization equations, based on the finite volume method (FVM), the third order Runge-Kutta method and the third order upwind scheme were introduced to handle the unsteady and vconvective terms in the governing equations. To determine the accuracy of the developed model, it was applied to the rectangular horizontal channel in a frictionless flow. The water depth and velocity obtained by the numerical model were found to agree closely with the exact solution. The model was also applied to the rectangular channel with both the symmetric and the non symmetric constriction. The velocity distribution of the flow and the propagation of the flood wave were simulated and the results well described the flow characteristics.