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

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.

Mathematics ◽  
2019 ◽  
Vol 7 (2) ◽  
pp. 211 ◽  
Author(s):  
Saufianim Jana Aksah ◽  
Zarina Ibrahim ◽  
Iskandar Mohd Zawawi

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.


2020 ◽  
Vol 51 (2) ◽  
pp. 123-136
Author(s):  
Iskandar Shah Mohd Zawawi

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.


Mathematics ◽  
2020 ◽  
Vol 8 (6) ◽  
pp. 914
Author(s):  
Zarina Bibi Ibrahim ◽  
Amiratul Ashikin Nasarudin

Recently, block backward differentiation formulas (BBDFs) are used successfully for solving stiff differential equations. In this article, a class of hybrid block backward differentiation formulas (HBBDFs) methods that possessed A –stability are constructed by reformulating the BBDFs for the numerical solution of stiff ordinary differential equations (ODEs). The stability and convergence of the new method are investigated. The methods are found to be zero-stable and consistent, hence the method is convergent. Comparisons between the proposed method with exact solutions and existing methods of similar type show that the new extension of the BBDFs improved the stability with acceptable degree of accuracy.


Author(s):  
Hande Uslu ◽  
Murat Sari ◽  
Tahir Cosgun

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.


Symmetry ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 952
Author(s):  
Amiratul Ashikin Nasarudin ◽  
Zarina Bibi Ibrahim ◽  
Haliza Rosali

In this research, a six-order, fully implicit Block Backward Differentiation Formula with two off-step points (BBDFO(6)), for the integration of first-order ordinary differential equations (ODEs) that exhibit stiffness, is proposed. The order, consistency and stability properties of the method are discussed, and the method is found to be zero stable and consistent. Hence, the method is convergent. The numerical comparisons with the existing methods of a similar type are given to demonstrate the accuracy of the derived method.


Author(s):  
I. S. M. Zawawi ◽  
Z. B. Ibrahim ◽  
F. Ismail ◽  
Z. A. Majid

This paper focuses on the derivation of diagonally implicit two-point block backward differentiation formulas (DI2BBDF) for solving first-order initial value problem (IVP) with two fixed points. The method approximates the solution at two points simultaneously. The implementation and the stability of the proposed method are also discussed. A performance of the DI2BBDF is compared with the existing methods.


Author(s):  
V. F. Edneral ◽  
O. D. Timofeevskaya

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