• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.1
• /
• pp.40-48
• /
• 2003

Large eddy simulation is used to investigate the compressible flow over a open cavity, The sub-grid scale stresses are modeled using the dynamic model. The compressible Navier-Stokes equations are solved with the sixth order accurate compact finite difference scheme in the space and the 4th order Runge-Kutta scheme in the time. The results show a typical flow pattern of the shear layer mode of oscillation over the cavity. The votical disturbances, the roll-up of vorticity, and impingement and scattering of vorticity at the downstream cavity edge can be seen in the shear layer. Predicted acoustic resonant frequency is in good agreement with that of the empirical formula. The mean flow streamlines are nearly horizontal along the mouth of the cavity. The pressure has its minimum value in the vortex core 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.

• Journal of Astronomy and Space Sciences
• /
• v.26 no.2
• /
• pp.199-210
• /
• 2009

Even though there are several Global Navigation Satellite Systems under development, only GPS and GLONASS are currently available for satellite positioning. In this study, GLONASS orbits were predicted from broadcast ephemeris using the 4th-order Runge-Kutta numerical integration. For accuracy validation, predicted orbits were compared with precise ephemeris. The RMS(Root Mean Square) and maximum 3-D errors were 14.3 km and 17.4 km for one-day predictions. In case of 7-day predictions, the RMS and maximum 3-D errors were 15.7 and 40.1 km, respectively. Also, the GLONASS satellite visibility predictions were compared with real observations, and they agree perfectly except for several epochs when the satellite signal was blocked by nearby buildings.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.2A
• /
• pp.43-49
• /
• 2011

In this study, to analyze the collision behaviors of the bridge super-structure and the vessel which the collision point is located far from its rotation center such as bridge of a vessel and equipments on a barge, the simplified collision model was proposed. The model was configured to denote the mass, stiffness and the nonlinear behaviors of the bridge and the vessel. The nonlinear equation of motions of the proposed model were numerically solved by 4th order Runge-Kutta method. The parametric studies were performed for various collision conditions by the standardized Korean barge vessel in term of barge width, and its effects to the maximum collision load of bridge were analyzed.

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

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

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