Adaptive Time Steps Runge-Kutta Methods: Comparative Analysis of Simulation Time in Nonlinear and Harmonically Excited Pendulum and Duffing Oscillators
Journal Title: Current Journal of Applied Science and Technology - Year 2017, Vol 19, Issue 5
Abstract
Time management without integrity compromise is an integral part of good engineering practice. The present study investigated for the required computation time in nonlinear and harmonically excited oscillators (Pendulum and Duffing). Adaptive time steps simulation of their governing equations with personal computer were performed by Runge-Kutta methods (RK2, RK3, RK4, RK5, RK5M) and one blend (RKB) comprising unsteady and steady solutions. The respective Pendulum and Duffing Poincare result at damp quality (4, 0.168), excitation amplitude (1.5, 0.21) and excitation frequency (2/3, 1.0) were used to validate the codes developed in FORTRAN environment. Actual simulation time was monitored at three different lengths of excitation periods (40000, 80000 and 120000) using the current time subroutine call command. Except for RK2, the validation Poincaré results compare well with the counterpart available in the literature for the oscillators. The actual computation time decrease rapidly with increasing order of Runge-Kutta method, but suffered relative increase for the blended method. The difference in computation time required between RK5 and RK5M is negligible for all studied cases. The Pendulum required longer actual computation time (4-115 seconds) than Duffing (3-52 seconds). The respective normalized range of time step for Pendulum and Duffing formed a simple average ratios {(1.0): (7.5): (13.3): (26.0): (24.0): (29.7)} and {(1.0): (3.7): (5.1): (9.1): (8.0): (11.2)} for RK2, RK3, RK4, RK5, RK5M and RKB. It is concluded that substantial time management can be achieved by the Runge-Kutta methods except RK2 that permitted strictly shorter simulation time steps leading to longer actual simulation time and consequently simulated largely an unacceptable Poincare result.
Authors and Affiliations
T. A. O. Salau, O. O. Ajide
Keywords
Related Articles
- EP ID EP318318
- DOI 10.9734/BJAST/2017/32388
- Views 129
- Downloads 0
Numerical Evaluation of the Axial Resistance Force in Steel Column Design
The finite element method (FEM) is a form of computational analysis that provides approximated results with acceptable accuracy. By using the FEM, developing theoretical models capable of properly analyzing the effects o...
Electric Field Induced Resistance Change in PEDOT-PSS Films
Resistance of PEDOT-PSS films changes considerably with time after application of DC voltage. Since, variation of resistance over time can lead to serious errors when used as a strain gauge or gas sensors or any device,...
Adoption of Cloud Computing Services in Healthcare Sectors: Special Attention to Private Hospitals in Colombo District, Sri Lanka
Cloud computing is the most important part of information technology that provides solutions to improve healthcare sectors. This innovation helps organizations to enhance or redesign the services at less expensive charge...
Biostratigraphy of Ochigbo–1 Well, Offshore Niger Delta; Evidence from Foraminifera, Spores and Pollen
This study involved the palynological and micropalaeontological studies of Ochigbo – 1 well, offshore Niger Delta. A total of Sixty-two (62) ditch cutting samples were composited at intervals of 100 ft covering a total d...
Land-Use/ Land-Cover Change Detection and Analysis in Aglar Watershed, Uttarakhand
Land-Use and Land-Cover (LULC) change has a significant impact on watershed processes such as hydrology, soil erosion, carbon sequestration, etc. While analyzing the interactions of various natural and human activities a...