• Transactions of the Korean Society of Automotive Engineers
• /
• v.12 no.4
• /
• pp.42-49
• /
• 2004

When a truck cab is conveyed at a constant speed by a hanger and immersed into the painting reservoir, it may fall off from the hanger by buoyancy. In order to reduce the buoyancy, on the bottom of a cab panel are holes placed, which allow paint to flow into the inside of a cab. In this study, a differential equation is derived which can be solved numerically by using 4th-Order Runge-Kutta method to calculate transient behavior of the buoyant force with sizes and locations of the holes given. The solution is utilized to optimally determine sizes and locations of the holes.

• Journal of Astronomy and Space Sciences
• /
• v.5 no.1
• /
• pp.19-30
• /
• 1988

The existence of the j-type collision periodic orbit is examined on the condition of mass ratio 0.9878449 and Jacobian constant 2.9∼3.4. Using the Birckhoff's regularization method and the 5th order Runge-Kutta variable step-sized numerical roution introduced by Fehlberg (1968). we test their periodicities. As the results, 4 j-type collision orbits and 5 peculiar orbits are represented. There are good agreements in this collision orbits with the relationship between the period and the shape of orbit proposed by Pinotsis Zikides(1984).

• Journal of the Korean Society of Tobacco Science
• /
• v.20 no.1
• /
• pp.87-94
• /
• 1998

After we made the computer source code with mathematical model of Muramatsu et al. that was expressed by the set of simultaneous first-order ordinary differential equations in evaporation-pyrolysis zone of cigarette, we simulated the distribution profiles of temperature and density of flue-cured tobacco. Those equations were solved numerically with the Runge-Kutta-Gill algorithm assuming step size of 0.025mm by Muramatsu at at,, but in this study the advanced algorithm of Runge-Kutta 4th Order assuming step size of 0.0005mm. The initial conditions and physical parameters of Muramatsu et at. were used for solving them. The calculated values corresponded well with results of Muramatsu et al., especially the gradient of the temperature profile increased with smoldering speed and the thickness of the evaporation-pyrolysis zone decreased with increasing of smoldering speed. On the other hand, the temperature gradient decreased with increasing of the effective thermal-conductivity value and the thickness of the evaporation-pyrolysis zone increased with the effective thermal-conductivity value.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results 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 exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• The Proceedings of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.7 no.2
• /
• pp.31-35
• /
• 1993

The radial density distributions of the positive column of strongly modulated low -pressure gas discharges in mercury - rare gas mixtures at 10 tom pressure have been studied theoretically. The current was modulated inusoidally with a modulation depth of 50%. Calculations have shown that the radial profile of the excited atoms is ditferent form 0th Bessel function $J_0$(2.4r/R) and the invertion of the radial distribution of the excited atom can occur at some frequency. The hybrid method of FDM and 2nd order Runge-Kutta meth od is used for solving differenzial equations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.2
• /
• pp.101-107
• /
• 2012

In this study, Constraint Force Equation Methodology (CFEM) is used to develop a multi-body dynamic analysis program with constraints. Seven constraint models are implemented to analyze constraint motions of multiple bodies. The augmented equations with the constraints are solved with the 4th order Runge-Kutta method for higher degree of accuracy. The analysis code is verified by comparing the analysis results of the motion of bodies with various constraints to published results.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.27-37
• /
• 1991

In this paper the nonlinear free surface flows for an axisymmetric submerged body oscillating beneath the free surface are solved and the forces acting on the body are calculated. A boundary integral method is applied to solve the axisymmetric boundary value problem and the Runge-Kutta 4-th order method is used for the time stepping of the free surface location. The nonlinear forces acting on the axisymmetric body are computed and compared with published results.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.228-232
• /
• 2005

Automobiles are very important as modern society is developed. Increase of the number of the automobiles cause environmental problem, that is, air pollution. So, many countries are adopting a environmental law. Automobile manufacturing companies have developing methods to prevent air pollution with increase of the efficiency of automotive engines. PCV(Positive Crankcase Ventilation) system which is one of them is made by the closed loop that consists of combustion chamber, crankcase, manifold suction tube and manifold. PCV valve is attached on manifold tube to control the flowrate of blowby gas. PCV valve is an important part in this system but it is difficult to design PCV valve which satisfies the required flowrate of blowby gas. In this study, our purpose is to help a PCV valve designer with the development of a design program. We used 4th order Runge-Kutta method and Bernoulli's equation to analyze the spool dynamic motion. By the comparison between our program and experiment, we think that a PCV designer can use our program in their work place.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.2
• /
• pp.167-175
• /
• 2011

The purpose of this study is to examine the impact that change in speed and modeling methods has on maglevs' runnability. The study constructed equations of motion on 4-DOF, 6DOF, and 10-DOF vehicles respectively and carried out numerical analysis, applying 4th Runge Kutta method, in order to run six different model maglev as changing the vehicles speed on the same bridge that had 2000 to 1 deflection. The analysis revealed that maglev's runnability improved as speed was lower and the specific model had higher number of bogey and EMS.