• Applied Chemistry for Engineering
• /
• v.9 no.7
• /
• pp.1079-1084
• /
• 1998

Reaction kinetics between 4,4'-dihexyl methane diisocyanate($H_{12}MDI$) and n-hexanol in toluene with dibutyltin dilaurate(DBTDL) as catalyst was studied by experimental measurements and mathematical modeling. Experiments were carried out at various temperatures, catalyst concentrations and [NCO]/[OH] ratios, and the reaction kinetics were described by two second-order reactions, the one between NCO and OH leading to urethane and the other between urethane and NCO leading to allophanate. The rate constants were estimated by the Runge-Kutta 4th-order method. Experiments and mathematical simulations showed a good agreement for various experimental conditions. The [allophanate]/[urethane] ratios at 90% conversion of initial NCO were estimated to be over 20% for most conditions employed in the present study, indicating that allophanate formation might significantly affect the properties of urethane polymers.

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

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

• Journal of Biomedical Engineering Research
• /
• v.24 no.4
• /
• pp.329-337
• /
• 2003

Blood in congenital or acquired AY fistula(arteriovenous fistula) flows from arteries directly to veins. detouring peripheral micro-circulation. This makes a great effect on the hemodynamics of human cardiovascular system. In this study, a computational method using lumped parameter mode) was proposed to simulate the cardiovascular hemodynamics of patients with acute AV fistula The cardiovascular system model with a fistula compartment in left lower limb was built using 17 standard lumped compartments. Using fourth order Runge-Kutta method. we solved numerically the unsteady linear set of the ordinary differential equations resulting from application of Kirchhoff's law to the lumped parameter hemodynamic model. The baroreceptor reflex system was implemented to explain the auto-regulation effect of the cardiovascular system with acute AV fistula.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.4
• /
• pp.204-211
• /
• 2011

A dynamic analysis procedure is developed to provide a better estimation of the dynamic responses of bridge during the passage of high speed railway vehicles. Particularly, a three dimensional numerical model including the structural interaction between high speed vehicles, bridges and railway endures to analyse accurately and evaluate with in-depth parametric studies for dynamic responses of various bridge span lengths running KTX railway locomotive up to increasing maximum speed(450km/h). Three dimensional frame element is used to model the simply supported pre-stressed concrete (PSC) box bridges for four span lengths(40~25m). Track irregularity employed as a stationary random process from the given spectral density functions and irregularities of both sides of the track are assumed to have high correlation. The high-speed railway vehicle (KTX) is used as 38-degree of freedom system. Three displacements (Vertical, lateral, and longitudinal) as well as three rotational components (Pitching, rolling, and yawing) are considered in the 38-degree of freedom model. The dynamic amplification factors are evaluated by the developed procedure under various traveling conditions, such as track irregularity camber, train speed and ballast. The dynamic analysis such as Newmark-${\beta}$ and Runge-Kutta methods which are able to analyse considering the dynamic impact factors are compared and contrasted.

• Journal of the Korean Magnetics Society
• /
• v.6 no.5
• /
• pp.281-286
• /
• 1996

For the accurate analysis of motional characteristics of electrical machines, it is needed to solve the motion equations together with the electromagnetic field equations. In this paper the sequential coupling of systems, the spring mass system and the electromagnetic system, is adopted. The induced current and the magnetic fields are calculated by finite element method(FEM) with given speed. And then, with the computed elec-tromagnetic force, the mechanical equations are solved by the Runge-Kutta method. The above two processes are repeated sequentially to obtain the time domain solutions. The resultant values are applied to the energy conservation law to prove the usefulness of the proposed sequential method.

• The Korean Journal of Rheology
• /
• v.8 no.2
• /
• pp.139-146
• /
• 1996

2차원 비점성 Euler 방정식의 Stuart 해에 의하여 주어지는 주기적인 유동장, 소위 말하여 Stuart와동장내에 놓여있는 아주 작은 비대칭성 입자의 다양한 부력 변수에 따른 침 강운동을 고찰하엿다. 본 연구에서 사용한 수치방법은 매질의 준전적인 Stokes 방정식과 입 자의 평형방정식을 연계하였다. 다양한 초기 배향각 또는 형상비에 따른 궤적이나 각변화를 예측하기 위하여 힘과 토오크 관계식을 입자의 운동방정식에 적용하였으며 4차 Runge-Kutta 방정식을 이용하여 입자의 운동을 규명하였다.