• Journal of the Korean Society of Propulsion Engineers
• /
• v.10 no.2
• /
• pp.1-14
• /
• 2006

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma} $ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Journal of Korea Water Resources Association
• /
• v.38 no.6 s.155
• /
• pp.429-438
• /
• 2005

In this study, a couple of ordinary differential equations which can describe random waves are derived from the Boussinesq equations. Incident random waves are generated by using the TMA(TEXEL storm, MARSEN, ARSLOE) shallow-water spectrum. The governing equations are integrated with the 4-th order Runge-Kutta method. By using newly derived wave equations, nonlinear energy interaction of propagating waves in constant depth is studied. The characteristics of random waves propagate over a sinusoidally varying topography lying on a sloping beach are also investigated numerically. Transmission and reflection of random waves are considerably affected by nonlinearity.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Tribology and Lubricants
• /
• v.33 no.6
• /
• pp.309-314
• /
• 2017

This paper presents a theoretical investigation of the dynamic characteristics of externally pressurized air pad bearings with closed loop grooves. These grooves are made on the surface of bearings to reduce the number of supply holes so that manufacturing costs can be reduced. The semi-implicit method is applied to calculate the time varying pressure profile on the air bearing surface owing to the advantages of numerical stability and fast time tracing characteristics. The static pressure of the groove bearings is much higher than that without grooves, so the groove bearings can provide high load carrying capacity. The equation of motion considering vertical motion and tilting motion are also solved using the Runge-Kutta 4th order method. By combining the semi-implicit method and the Runge-Kutta method, fast calculations of the dynamic behavior of the air bearing can be achieved. The variations of bearing reaction force, air film reaction moment, height, and tilting angle are investigated for the step force input, which is 20% higher than the bearing reaction, when the nominal clearance is 6 mm. The effect of the groove width and the groove depth are investigated by calculating the dynamic behavior. The possibility of the air hammering with the depth of the groove is found and discussed.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.1 no.1
• /
• pp.27-38
• /
• 1995

In this paper, a program for the calculation of GZ curve for a ship in waves is developed and GZ curves for a ferry in the still water and in waves are calculated. And the added mass, damping, restoring forces and wave exciting forces are calculated by using the strip theory given by Salvesen, Tuck, Faltinsen. Capsizing simulations are perfoned in consideration if the nonlinear restoring forces of the ship in waves by using the Runge-Kutta 4-th method.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1997.04a
• /
• pp.204-209
• /
• 1997

This paper reports on the theoretical study of the vane motions in a positive displacement vane pump. Vane detachment cause the pressure fluctuation, noise, wear in cam ring, and decrease the volumetric efficiency. Dynamic equation of vane motion and flow continuity equation have been modeled and solved simultaneously using 4th order Runge-Kutta method. As results of analysis, vane detachment occurs due to pressure overshoot by excess compression in the pumping chamber. Amount of vane detachment has been reduced by decreasing the pressure overshoot.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.39 no.5
• /
• pp.516-518
• /
• 1990

We have derived the equation which involves the variation of electron energy with time in a lowly ionized plasma when a low alternating electric field is applied. We consider only elastic collisions between electrons and neutral atoms. This equation is solved using the 4th-order Runge-Kutta method, and applied to argon gas discharge which is driven by source frequency of 100, 1K, 10K, 100K, and 1M (Hz). The results show that the variation of electron energy becomes flat with higher frequencies.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. 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, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. 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, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.