scholarly journals Qualitative behavior of stiff ODEs through a stochastic approach

2020 ◽  
Vol 10 (2) ◽  
pp. 181-187
Author(s):  
Hande Uslu ◽  
Murat Sari ◽  
Tahir Cosgun
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Real Life ◽  
Stochastic Approach ◽  
Computational Time ◽  
Qualitative Behavior ◽  
Computational Error ◽  
Linear Ordinary Differential Equations ◽  
Stiff Ordinary Differential Equations ◽  
Simulation Techniques

In the last few decades, stiff differential equations have attracted a great deal of interest from academic society, because much of the real life is covered by stiff behavior. In addition to importance of producing model equations, capturing an exact behavior of the problem by dealing with a solution method is also handling issue. Although there are many explicit and implicit numerical methods for solving them, those methods cannot be properly applied due to their computational time, computational error or effort spent for construction of a structure. Therefore, simulation techniques can be taken into account in capturing the stiff behavior. In this respect, this study aims at analyzing stiff processes through stochastic approaches. Thus, a Monte Carlo based algorithm has been presented for solving some stiff ordinary differential equations and system of stiff linear ordinary differential equations. The produced results have been qualitatively and quantitatively discussed.

scholarly journals BBDF-Alpha for solving stiff ordinary differential equations with oscillating solutions

2020 ◽  
Vol 51 (2) ◽  
pp. 123-136
Author(s):  
Iskandar Shah Mohd Zawawi
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Numerical Experiments ◽  
Stability Region ◽  
Newton Iteration ◽  
Computational Time ◽  
First Order ◽  
Stiff Ordinary Differential Equations ◽  
Oscillating Solutions ◽  
The Stability

In this paper, the block backward differentiation α formulas (BBDF-α) is derived for solving first order stiff ordinary differential equations with oscillating solutions. The consistency and zero stability conditions are investigated to prove the convergence of the method. The stability region in the entire negative half plane shows that the derived method is A-stable for certain values of α. The implementation of the method using Newton iteration is also discussed. Several numerical experiments are conducted to demonstrate the performance of the method in terms of accuracy and computational time.

scholarly journals Formulative Visualization of Numerical Methods for Solving Non-Linear Ordinary Differential Equations

2021 ◽  
Vol 2 (2) ◽  
pp. 79-88
Author(s):  
Jeevan Kafle ◽  
Bhogendra Kumar Thakur ◽  
Grishma Acharya
Keyword(s):  
Numerical Methods ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Real Life ◽  
Euler Method ◽  
Linear Ordinary Differential Equations ◽  
Life Problems ◽  
Backward Euler ◽  
Non Linear ◽  
Runge Kutta

Many physical problems in the real world are frequently modeled by ordinary differential equations (ODEs). Real-life problems are usually non-linear, numerical methods are therefore needed to approximate their solution. We consider different numerical methods viz., Explicit (Forward) and Implicit (Backward) Euler method, Classical second-order Runge-Kutta (RK2) method (Heun’s method or Improved Euler method), Third-order Runge-Kutta (RK3) method, Fourth-order Runge-Kutta (RK4) method, and Butcher fifth-order Runge-Kutta (BRK5) method which are popular classical iteration methods of approximating solutions of ODEs. Moreover, an intuitive explanation of those methods is also be presented, comparing among them and also with exact solutions with necessary visualizations. Finally, we analyze the error and accuracy of these methods with the help of suitable mathematical programming software.

scholarly journals Fixed Coefficient A ( α ) Stable Block Backward Differentiation Formulas for Stiff Ordinary Differential Equations

Symmetry ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 846 ◽  
Author(s):  
Zarina Ibrahim ◽  
Nursyazwani Mohd Noor ◽  
Khairil Othman
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Execution Time ◽  
Stability Region ◽  
Scientific Literature ◽  
Necessary Conditions ◽  
Computational Time ◽  
Stiff Ordinary Differential Equations ◽  
Backward Differentiation Formulas ◽  
Differentiation Formulas

The main contribution in this paper is to construct an implicit fixed coefficient Block Backward Differentiation Formulas denoted as A ( α ) -BBDF with equal intervals for solving stiff ordinary differential equations (ODEs). To avoid calculating the differentiation coefficients at each step of the integration, the coefficients of the formulas will be stored, with the intention of optimizing the performance in terms of precision and computational time. The plots of their A ( α ) stability region are provided, and the order of the method is also verified. The necessary conditions for convergence, such as the consistency and zero stability of the method, are also discussed. The numerical results clearly showed the efficiency of the method in terms of accuracy and execution time as compared to other existing methods in the scientific literature.

scholarly journals Stability Analysis of Singly Diagonally Implicit Block Backward Differentiation Formulas for Stiff Ordinary Differential Equations

Mathematics ◽  
2019 ◽  
Vol 7 (2) ◽  
pp. 211 ◽  
Author(s):  
Saufianim Jana Aksah ◽  
Zarina Ibrahim ◽  
Iskandar Mohd Zawawi
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Triangular Matrix ◽  
Practical Significance ◽  
Computational Time ◽  
Step Size ◽  
Time Step ◽  
Stiff Ordinary Differential Equations ◽  
Backward Differentiation Formulas ◽  
Differentiation Formulas

In this research, a singly diagonally implicit block backward differentiation formulas (SDIBBDF) for solving stiff ordinary differential equations (ODEs) is proposed. The formula reduced a fully implicit method to lower triangular matrix with equal diagonal elements which will results in only one evaluation of the Jacobian and one LU decomposition for each time step. For the SDIBBDF method to have practical significance in solving stiff problems, its stability region must at least cover almost the whole of the negative half plane. Step size restriction of the proposed method have to be considered in order to ensure stability of the method in computing numerical results. Efficiency of the SDIBBDF method in solving stiff ODEs is justified when it managed to outperform the existing methods for both accuracy and computational time.

Normal form method in search for periodic solutions of ordinary differential equations. Case of the fourth order

2018 ◽  
pp. 24-34
Author(s):  
V. F. Edneral ◽  
O. D. Timofeevskaya
Keyword(s):  
Normal Form ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Periodic Solutions ◽  
Fourth Order ◽  
Numerical Solutions ◽  
Numerical Calculations ◽  
Computational Time ◽  
Resonant Normal Form ◽  
Normal Form Method

Introduction:The method of resonant normal form is based on reducing a system of nonlinear ordinary differential equations to a simpler form, easier to explore. Moreover, for a number of autonomous nonlinear problems, it is possible to obtain explicit formulas which approximate numerical calculations of families of their periodic solutions. Replacing numerical calculations with their precalculated formulas leads to significant savings in computational time. Similar calculations were made earlier, but their accuracy was insufficient, and their complexity was very high.Purpose:Application of the resonant normal form method and a software package developed for these purposes to fourth-order systems in order to increase the calculation speed.Results:It has been shown that with the help of a single algorithm it is possible to study equations of high orders (4th and higher). Comparing the tabulation of the obtained formulas with the numerical solutions of the corresponding equations shows good quantitative agreement. Moreover, the speed of calculation by prepared approximating formulas is orders of magnitude greater than the numerical calculation speed. The obtained approximations can also be successfully applied to unstable solutions. For example, in the Henon — Heyles system, periodic solutions are surrounded by chaotic solutions and, when numerically integrated, the algorithms are often unstable on them.Practical relevance:The developed approach can be used in the simulation of physical and biological systems.

Sign in / Sign up

Export Citation Format

Share Document