• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.1
• /
• pp.69-77
• /
• 1984

In this paper, the governing differential equations for the free vibration of uniform parabolic arches are derived on the basis of equilibrium equations of a small element of arch rib and the D'Alembert principle. A trial eigen value method is used for determining the natural frequencies and mode shapes. And the Runge-Kutta fourth order integration technique is also used in this method to perform the integration of the differential equations. A detailed study is made of the first mode for the symmetrical and anti-symmetrical vibrations of hinged arches with the Span length equal to 10 m. The effects of the rise of arch, the radius of gyration and the rotary inertia on free vibrations are presented in detail in curves and table.

• 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 Institute of Convergence Signal Processing
• /
• v.9 no.1
• /
• pp.82-88
• /
• 2008

In this paper, Runge-Kutta (RK)-Butcher algorithm is proposed to study the periodic and oscillatory problems. Simulation results obtained using RK-Butcher algorithms and the classical fourth order Runge-Kutta (RK(4)) methods are compared with the exact solutions of a few periodic and oscillatory problems to confirm the performance of the proposed algorithm. The simulation results from RK-Butcher algorithms are always found to be very close to the exact solutions of these problems. Further, it is found that the RK-Butcher algorithm is superior when compared to RK(4) methods in terms of accuracy. The RK-Butcher algorithm can be easily implemented in a programming language and a more accurate solution may be obtained for any length of time. RK-Butcher algorithm is applicable as a good numerical algorithm for studying the problems of orbit and two body as it gives the nearly identical solutions.

• Journal of the Korean Physical Society
• /
• v.73 no.9
• /
• pp.1303-1309
• /
• 2018

In the present paper, we study the entropy analysis of boundary layer flow over a slender stretching sheet under the action of a non uniform magnetic field that is acting perpendicular to the flow direction. The effects of viscous dissipation and Joule heating are included in the energy equation. Using similarity transformation technique the momentum and thermal boundary layer equations to a system of nonlinear differential equations. Numerical solutions are obtained using the shooting and fourth-order Runge-Kutta method. The expressions for the entropy generation number and Bejan number are also obtained using a suggested similarity transformation. The main objective of this article is to investigate the effects of different governing parameters such as the magnetic parameter ($M^2$), Prandtl number (Pr), Eckert number (Ec), velocity index parameter (m), wall thickness parameter (${\alpha}$), temperature difference parameter (${\Omega}$), entropy generation number (Ns) and Bejan number (Be). All these effects are portrayed graphically and discussed in detail. The analysis reveals that entropy generation reduces with decreasing wall thickness parameter and increasing temperature difference between the stretching sheet and the fluid outside the boundary layer. The viscous and magnetic irreversibilities are dominant in the vicinity of the stretching surface.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2060-2065
• /
• 2003

Aeolian tone generation from a two dimensional circular cylinder is numerically investigated via direct numerical simulation and hydrodynamic-acoustic splitting method. All governing equation are spatially discretized with the sixth-order compact scheme and fourth-order Runge-Kutta method to avoid excessive numerical dissipations and dispersions of acoustic quantities. Comparisons of two results show that the previous splitting method can not accurately predict the aeroacoustic noise of wall bounded shear flow. In this study, a perturbation viscous term and a new energy equation have been developed. This modified splitting method accurately predicts aeroacoustic noise from wall-bounded shear flow. The present results agree very well with the direct numerical simulation solution.

• Structural Engineering and Mechanics
• /
• v.65 no.5
• /
• pp.513-521
• /
• 2018

In this paper, the resonance response of spar-type floating platform in coupled heave and pitch motion is investigated using a CPU time-effective numerical method. A coupled nonlinear 2-DOF equation of motion is derived based on the potential wave theory and the rigid-body hydrodynamics. The transient responses are solved by the fourth-order Runge-Kutta (RK4) method and transformed to the frequency responses by the digital Fourier transform (DFT), and the first-order approximation of heave response is analytically derived. Through the numerical experiments, the theoretical derivation and the numerical formulation are verified from the comparison with the commercial software AQWA. And, the frequencies of resonance arising from the nonlinear coupling between heave and pitch motions are investigated and justified from the comparison with the analytically derived first-order approximation of heave response.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.197-199
• /
• 2003

In this work, the Greitzer's B-parameter model is applied for analyzing the stall and surge characteristics. The four parameters in the model are highlighted in order to establish the influence of each parameter on the system. First of all, the governing equations of stall and surge behavior are solved numerically using fourth-order Runge-Kutta method. The Taguchi method is then used to analyze the results generated to obtain the extent of effects of the parameters on the system by varying the parameters in a series of combinations. Finally, a thorough analysis is carried out on the results generated from the Taguchi method and the graphs.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.679-690
• /
• 1997

The risk assessment model for hydrlolgic safety analysis of dam and levee in developed by using Monte-Carlo and AFOSM (Advanced First-Order Second-Moment) method. The fault tree analysis and four phases approach are presented for the safety eveluation of risk of dam and levee. The risk model consists of rainfall-runoff analysis, reservoir routing and channel routing considering the variations in the model parameter. For the rainfall-runoff analysis, KRRL method is adopted with 200-year precipitation and PMP (Probable Maximum Precipitation). Reservoir routing is performed by fourth order Runge-Kutta method and channel routing by standard step method. The suggested model will contribute to safety evaluation of dam and levee and their rehabilitation decision problem.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.603-616
• /
• 2018

In this study, the acceleration vector in each time step is assumed to be a mth order time polynomial. By using the initial conditions, satisfying the equation of motion at both ends of the time step and minimizing the square of the residual vector, the m+3 unknown coefficients are determined. The order of accuracy for this approach is m+1, and it has a very low dispersion error. Moreover, the period error of the new technique is almost zero, and it is considerably smaller than the members of the Newmark method. The proposed scheme has an appropriate domain of stability, which is greater than that of the central difference and linear acceleration techniques. The numerical tests highlight the improved performance of the new algorithm over the fourth-order Runge-Kutta, central difference, linear and average acceleration methods.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.131-145
• /
• 2013

In this paper, we present an initial value technique for solving singularly perturbed differential difference equations with a boundary layer at one end point. Taylor's series is used to tackle the terms containing shift provided the shift is of small order of singular perturbation parameter and obtained a singularly perturbed boundary value problem. This singularly perturbed boundary value problem is replaced by a pair of initial value problems. Classical fourth order Runge-Kutta method is used to solve these initial value problems. The effect of small shift on the boundary layer solution in both the cases, i.e., the boundary layer on the left side as well as the right side is discussed by considering numerical experiments. Several numerical examples are solved to demonstate the applicability of the method.