Next Article in Journal
Applying Artificial Intelligence in Cryptocurrency Markets: A Survey
Next Article in Special Issue
Some Inexact Version of the Iterative Exponential Method for Solving Nonsmooth Equations
Previous Article in Journal
Special Issue “Selected Algorithmic Papers From CSR 2020”
Previous Article in Special Issue
Dynamics and Stability on a Family of Optimal Fourth-Order Iterative Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems

1
School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X01, Scottsville, Pietermaritzburg 3201, South Africa
2
Mathematics Department, University of Eswatini, Private Bag 4, Kwaluseni M201, Eswatini
Algorithms 2022, 15(11), 427; https://doi.org/10.3390/a15110427
Submission received: 19 October 2022 / Revised: 5 November 2022 / Accepted: 9 November 2022 / Published: 14 November 2022

Abstract

:
This study presents a new variant of the hybrid block methods (HBMs) for solving initial value problems (IVPs). The overlapping hybrid block technique is developed by changing each integrating block of the HBM to incorporate the penultimate intra-step point of the previous block. In this paper, we present preliminary results obtained by applying the overlapping HBM to IVPs of the first order, utilizing equally spaced grid points and optimal points that maximize the local truncation errors of the main formulas at the intersection of each integration block. It is proven that the novel method reduces the local truncation error by at least one order of the integration step size, O(h). In order to demonstrate the superiority of the suggested method, numerical experimentation results were compared to the corresponding HBM based on the standard non-overlapping grid. It is established that the proposed method is more accurate than HBM versions of the same order that have been published in the literature.

1. Introduction

This study presents a variation of the hybrid block methods (HBMs) for integrating systems of initial value problems of the first order (IVPs). It is now well established that HBMs are effective at solving IVPs that may or may not be stiff. Researchers are now concentrating on enhancing these approaches to increase their efficacy and extend them to other forms of differential equations, such as boundary value problems. Several authors have developed optimised hybrid block methods in recent years by selecting intra-step points that minimise the local truncation error of one or more of the formulas in the set of hybrid block methods (see, for instance, [1,2,3,4,5]). The hybrid block methods were applied with an adaptive strategy in [6,7] in an effort to increase their efficiency and circumvent the difficulties associated with short step sizes. Other researchers (see [8,9]) have implemented the optimised hybrid block methods in nonlinear partial differential equations using cubic spline basis functions. In this study, the hybrid block method is modified by overlapping adjacent integration blocks at the last intra-step point. This adjustment produces a A -stable version with a truncation error at least O(h) greater than the comparable HBM without overlap.
The performance of the proposed overlapping HBM is tested utilising an error analysis of numerically challenging solutions to different nonlinear IVPs. In order to determine the approach’s robustness, it is also evaluated on IVPs with extreme stiffness.
The remainder of this paper is organised into four sections. The derivation of the overlapping hybrid block methods is covered in Section 2. Section 3 examines the error and stability properties of the proposed method. Section 4 presents the experimental and numerical results, and Section 5 provides a summary and conclusion.

2. Derivation of the Method

Consider the initial value problem defined on t [ 0 , T ] by
y ( t ) = f ( t , y ( t ) ) , y ( t 0 ) = y 0
where f ( t , y ) is a general, continuous nonlinear function.
Let h be a prescribed step size on a mesh I defined by
I = { t n : t 0 < t 1 < < t N = T } .
We choose m + 1 collocation points, defined by
t n + p i = t 0 + h p i , i = 0 , 1 , 2 , , m
in the first interval [ t 0 , t 1 ] . For the intervals, [ t n , t n + 1 ] , where 1 n N 1 , we choose m + 2 collocation points:
{ t n p 1 = t n h p 1 } t n + p i = t n + h p i : i = 0 , 1 , 2 , , m .
The parameters p i are called collocation parameters, and they satisfy the condition
0 = p 0 < p 1 < p 2 < < p m 2 < p m 1 < p m = 1 .
An illustrative representation of the grid is shown in Figure 1, for the special case when m = 3 and N = 4 . It is worth noting that, for n > 0 , the grid I n overlaps with one grid point ( t n 1 + p m 1 = t n p 1 ) from the previous grid I n 1 . For this reason, the method will be referred to as the overlapping hybrid block method.
In the non-overlapping block I 1 , we suppose that the solution can be approximated by the polynomial:
y ( t 0 + s h ) , 0 s 1 ,
that satisfies the collocation conditions:
y ( t p i ) = f p i .
By the Lagrange interpolation formula, the function y ( t 0 + s h ) can be approximated as
y ( t 0 + s h ) = j = 0 m y ( t p j ) j ( s )
where
j ( s ) = i = 0 , i j m s p i p j p i .
Using the collocation Equation (3), we obtain
y ( t 0 + s h ) = j = 0 m f p j j ( s )
which is integrated from 0 to p i to give
y p i = y 0 + h j = 0 m β i , j f p j
where
β i , j = 0 p i j ( τ ) d τ , j = 0 , 1 , , m .
The matrix form of (6) is
A 1 Y 1 = A 0 Y 0 + h ( B 0 F 0 + B 1 F 1 )
where
A 1 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 , A 0 = 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
B 0 = β 1 , 0 0 0 0 β 2 , 0 0 0 0 β 3 , 0 0 0 0 β m , 0 0 0 0 , B 1 = β 1 , 1 β 1 , 2 β 1 , 3 β 1 , m β 2 , 1 β 2 , 2 β 2 , 3 β 2 , m β 3 , 1 β 3 , 2 β 3 , 3 β 3 , m β m , 1 β m , 2 β m , 3 β m , m
and
F 0 = f 0 , f p 1 , , f p m 1 , F 1 = f p 1 , f p 2 , , f p m , Y 0 = y 0 , y p 1 , , y p m 1 , Y 1 = y p 1 , y p 2 , , y p m .
In the block I n = [ t n , t n + 1 ] ( n > 0 ), we consider the collocation parameters
p 1 < 0 < p 1 < p 2 < < p m 1 < p m
and assume the following collocation conditions:
y ( t n + p i ) = f n + p i , y ( t n p 1 ) = f n p 1
for i = 0 , 1 , , m . By the Lagrange interpolation formula, the function y ( t n + s h ) can be approximated as
y ( t n + s h ) = j = 1 m y ( t n + p j ) j ( s ) .
Integrating from 0 to p i gives
y n + p i = y n + h j = 1 m β ^ i , j f n + p j
where
β ^ i , j = 0 p i j ( τ ) d τ , j = 1 , 0 , 1 , , m .
Equation (10) can be written in matrix form as
A 1 Y n + 1 = A 0 Y n + h ( B ^ 0 F n + B ^ 1 F n + 1 )
where
Y n + 1 = y n + p 1 y n + p 2 y n + p m , Y n = y n y n p 1 y n p m 1 , F n + 1 = f n + p 1 f n + p 2 f n + p m , F n = f n f n p 1 f n p m 1
B ^ 0 = β ^ 1 , 0 β ^ 1 , 1 0 β ^ 2 , 0 β ^ 2 , 1 0 β ^ m , 0 β ^ m , 1 0 B ^ 1 = β ^ 1 , 1 β ^ 1 , 2 β ^ 1 , m β ^ 2 , 1 β ^ 2 , 2 β ^ 2 , m β ^ m , 1 β ^ m , 2 β ^ m , m

3. Analysis of the Method

In this section, we present the truncation error and stability properties of the proposed overlapping hybrid block method.

3.1. Local Truncation Error

Theorem 1. 
The local truncation error, in the integration block [ t n , t n + 1 ] , for the hybrid block method defined as
y n + p i = y n + h j = 0 m β i , j f n + p j , β i , j = 0 p i j ( τ ) d τ
is given by
L i [ z ( t n ) ; h ] = h m + 2 ( m + 2 ) ! p i m + 2 ( m + 2 ) j = 0 m β i , j p j m + 1 z ( m + 2 ) ( t n ) + O ( h m + 3 )
for i = 1 , 2 , , m , when considering the grid points:
t n < t n + p 1 < t n + p 1 < < t n + p m .
Proof. 
Given a sufficiently differentiable function z ( t ) , the truncation error of the hybrid block method can be written in terms of a linear operator L i as
L i [ z ( t n ) ; h ] = z ( t n + h p i ) z ( t n ) h j = 0 m β i , j z ( t n + h p j ) , i = 1 , 2 , , m .
Using Taylor series to expand (14) gives
L i [ z ( t n ) ; h ] = k = 1 K p i k k ! h k z ( k ) ( t n ) k = 1 K h k k k ! j = 0 m β i , j p j k 1 z ( k ) ( t n ) + O ( h K + 1 ) = k = 2 m + 1 h k k ! p i k k j = 0 m β i , j p j k 1 z ( k ) ( t n ) + k = m + 2 K h k k ! p i k k j = 0 m β i , j p j k 1 z ( k ) ( t n ) + O ( h K + 1 )
where K m + 2 is a positive integer.
Using the identity:
j = 0 m j ( x ) = 1
we note that
j = 0 m β i , j = j = 0 m 0 p i j ( x ) d x = p i
and
j = 0 m β i , j p j k 1 = j = 0 m 0 p i j ( x ) d x p j k 1 = p i k k , f o r k = 2 , , m + 1 .
In a related study on implicit Runge–Kutta processes, the identity (16) was also reported by [10]. Substituting the identities in (15) gives the result in (13). □
As an immediate consequence, we obtain the following result, for the truncation error in the integration blocks [ t n , t n + 1 ] ( n > 0 ):
Corollary 1. 
The local truncation error, in the integration block [ t n , t n + 1 ] , for the one-point overlapping hybrid block method:
y n + p i = y n + h j = 1 m β ^ i , j f n + p j , β ^ i , j = 0 p i j ( τ ) d τ
is
L i [ z ( t n ) ; h ] = h m + 3 ( m + 3 ) ! p i m + 3 ( m + 3 ) j = 0 m + 1 β ^ i , j 1 p j 1 m + 2 z ( m + 3 ) ( t n ) + O ( h m + 4 ) .
Proof. 
The proof is similar to the one for Theorem 1, given above, and uses the identity:
j = 0 m + 1 β ^ i , j 1 p j 1 k 1 = j = 0 m + 1 0 p i j 1 ( x ) d x p j 1 k 1 = p i k k , f o r   k = 2 , , m + 2
after considering the m + 2 grid points:
t n p 1 < t n < t n + p 1 < t n + p 1 < < t n + p m .
By minimizing the local truncation error of the main method at t n + 1 (i.e., p m = 1 ), it was shown in [11] that the optimal intra-step points are the roots of the polynomial:
Q i ( p ) = 1 i ( i + 1 ) k = 1 i k ( 1 ) i + k p k 1 ( i + k ) ! ( k ! ) 2 ( i k ) !
which is proportional to the first derivative of the shifted Legendre polynomials of degree m given by
P i * ( t ) = k = 0 i ( 1 ) i + k t k ( i + k ) ! ( k ! ) 2 ( i k ) ! .
The shifted Legendre polynomials are defined in [ 0 , 1 ] and have the property:
P i * ( 1 t ) = ( 1 ) i P i * ( t ) , i = 0 , 1 , 2 , .
Furthermore, P i * ( t ) has i distinct roots in ( 0 , 1 ) . By expanding (19), we obtain the polynomials that give optimal intra-step points with m = 2 , 3 , 4 , 5 as follows:
Q 2 ( p ) = 2 p 1 , Q 3 ( p ) = 5 p 2 5 p + 1 , Q 4 ( p ) = 14 p 3 21 p 2 + 9 p 1 , Q 5 ( p ) = 42 p 4 84 p 3 + 56 p 2 14 p + 1 , Q 6 ( p ) = 132 p 5 330 p 4 + 300 p 3 120 p 2 + 20 p 1 .
As a consequence of Equation (21), the optimal roots (intra-step points) have the form
p 1 , p 2 , , , p m 2 1 2 , 1 p m 2 1 2 , , 1 p 2 , 1 p 1
when m is odd and
p 1 , p 2 , , , p m 2 1 , 1 2 , 1 p m 2 1 , , 1 p 2 , 1 p 1
when m is even.
The truncation errors for overlapping and non-overlapping HBM algorithms with equally spaced and optimally spaced intra-step points are compared in Table 1 for m = 3 , 4 , 5 . The table reveals that the truncation errors of the overlapping HBM are at least h orders of magnitude smaller than those of the non-overlapping HBM. This is consistent with the observations made in the preceding theorems regarding the explicit formulas for local truncation errors.

3.2. Zero Stability

In the limit when h 0 , the zero stability of the hybrid methods for all blocks ( n = 0 , 1 , , N 1 ) can be analysed from the zeros of the polynomial that is computed from
ρ ( λ ) = d e t A 1 λ A 0 .
Zero stability is confirmed if the roots of ρ ( λ ) satisfy | λ j | 1 and all the roots with
| λ j | = 1 have a multiplicity of 1. The characteristic polynomial for the hybrid block methods considered in this study is
ρ ( λ ) = λ m ( λ 1 ) .
Thus, the proposed overlapping hybrid block methods are zero stable for any selection of intra-step points. In a related investigation, Reference [1] also validated the zero stability of the HBMs.

3.3. Absolute Stability

If errors introduced in one time step do not accumulate in successive time steps, the general stability of a numerical approach is confirmed. A region of absolute stability is defined as the set of points z C such that the roots of the characteristic equation associated with the test equation y = λ y ( λ < 0 ) lie within a unit circle. For hybrid block methods, the absolute stability in the first non-overlapping interval is well established (see, for example, [5]). When this test problem is applied to the overlapping grid methods (11), the following matrix equation is obtained:
Y n + 1 = R ( z ) Y n
where z = λ h and R ( z ) is a matrix equation given by
R ( z ) = ( A 1 z B ^ 1 ) 1 ( A 0 + z B ^ 0 ) .
By finding the eigenvalues of the matrix R ( z ) , one can derive the stability function H ( z ) . Figure 2 depicts, for illustrative purposes, the absolute stability region of the proposed overlapping grid hybrid block method created from optimally spaced and equally spaced intra-step points for m = 3 . Given that the shaded area encompasses the entirety of the left half-plane, the absolute stability of the overlapping HBM is confirmed for both equally spaced and optimal intra-step points.

4. Numerical Experimentation and Results

In this section, we implement the proposed overlapping hybrid block method to generate numerical solutions for problems that are typically difficult to solve numerically using standard methods. Before employing the hybrid block approach to integrate the nonlinear IVPs, they are first linearized. For a general nonlinear system of first-order IVPs, the linearisation method is detailed below.
Consider the system of s nonlinear first-order IVPs of the form
y k = f k ( t , y 1 , , y s ) = L k ( t , y 1 , , y k 1 , y k + 1 , , y s ) y k + N k ( t , y 1 , , y k 1 , y k , y k + 1 , y s )
where L k ( ) is the coefficient of the state variable y k , and N k ( ) is the component of the k-th equation that is not linear in y k , for k = 1 , 2 , , s . The hybrid block method is sequentially applied iteratively to each equation in the system (25), with the solution of the k-th equation used immediately when solving the ( k + 1 ) th equation. Consequently, the iterative scheme is developed as
y k , r + 1 = L k ( t , y 1 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , , y s , r ) y k , r + 1 + N k ( t , y 1 , r + 1 , , y k 1 , r + 1 , y k , r , , y s , r ) .
This linearization method employs the Gauss–Seidel method for decoupling and solving large systems of nonlinear equations. The method has been adapted from the waveform relaxation method that has been reported in the literature (see, for example, [12,13]).
Noting that (26) is linear, it can be applied to the HBM method after setting
f ( t , y ) = ϕ ( t ) y + ψ ( t )
where ϕ = L ( ) and ψ = N ( ) are known functions of t. As a result, the matrix form of the HBM, (11), becomes
A 1 Y n + 1 = A 0 Y n + h B ^ 0 ( ϕ n ( d ) Y n + ψ n ) + h B ^ 1 ( ϕ n + 1 ( d ) Y n + 1 + ψ n + 1 )
or
P n Y n + 1 = Q n
where
P n = A 1 h B ^ 1 ϕ n + 1 ( d ) , Q n = A 0 + h B ^ 0 ϕ n ( d ) Y n + h B ^ 0 ψ n + B ^ 1 ψ n + 1
ψ n = ψ n , ψ n p 1 , , ψ n p m 1 T , ψ n + 1 = ψ n + p 1 , ψ n + p 2 , , ψ n + p m T ,
ϕ n ( d ) = ϕ n ϕ n p 1 ϕ n p m 1 , ϕ n + 1 ( d ) = ϕ n + p 1 ϕ n + p 2 ϕ n + p m
All numerical examples in this study were solved using (28) applied iteratively over 20 iterations for each integration block [ t n , t n + 1 ] . Through numerical experimentation, it was determined that sufficient convergence would have occurred prior to the 20th iteration level. All numerical results presented in this study were generated by iterating the scheme (28) up to r = 20 for all n values. In light of the linearisation of nonlinear IVPs via partitioning, the hybrid block method defined by (28) is referred to as the linear partitioning overlapping hybrid block method. The method is computationally efficient and easy to implement. Due to the lack of derivatives and Jacobian matrices, it circumvents the challenges of Newton-based iterative techniques that rely on derivatives.
In the examples, given below, the linear partitioning overlapping hybrid block method was applied using equally spaced and the shifted Legendre-based optimal intra-step points. In both cases, m = 3 and m = 5 were used for illustrative purposes. The intra-step points corresponding to these values of m are given in Table 2 below.
Using m = 3 as an example, the matrix parameters for the optimal overlapping hybrid block approach are as follows:
B 0 = 5 + 11 120 0 0 11 5 120 0 0 1 12 0 0 , B 1 = 25 5 120 25 13 5 120 5 1 120 13 5 + 25 120 5 + 25 120 5 1 120 5 12 5 12 1 12
which is applied on the scheme:
A 1 Y 1 = A 0 Y 0 + h ( B 0 F 0 + B 1 F 1 )
in the first interval [ t 0 , t 1 ] and
B ^ 0 = 11 5 + 61 600 5 5 9 2200 0 49 11 5 600 5 5 + 9 2200 0 1 12 0 0 , B ^ 1 = 55 4 5 300 137 65 5 600 37 5 61 6600 17 5 + 35 150 5 5 + 113 600 37 5 49 6600 5 12 5 12 1 12
which is applied in [ t n , t n + 1 ] , for the overlapping hybrid block scheme:
A 1 Y 1 = A 0 Y 0 + h ( B ^ 0 F 0 + B ^ 1 F 1 )
for n = 1 , 2 , , N 1 . The proposed overlapping grid approach is referred to as the overlapping optimal hybrid block method (OOHBM) for the purposes of this paper. The original optimal hybrid block method (OHBM), without overlap, is constructed by applying Equation (29) to all the integration blocks [ t n , t n + 1 ] for n = 0 , 1 , 2 , , N 1 . It is important to note that the size of the matrices in both the original HBM scheme and the overlapping HBM scheme is the same. This indicates that the computational effort necessary to solve the IVPs will be roughly equivalent for both methods.

4.1. Example 1

Consider the following example that belongs to the class of problems with the so-called Runge phenomenon:
y = 144 ( 1 2 t ) y 2 , y ( 0 ) = 1 37 .
The exact solution of Example (31) is y ( t ) = 1 36 ( 2 t 1 ) 2 + 1 . In this example, the linear partitioning iterative method with
ϕ ( t , y r ) = 144 ( 1 2 t ) y r 2 , ψ ( t , y r ) = 0
was used with the HBM schemes to solve (31) using h = 0.01 and up to 20 iterations, within the integration domain [0,1]. The convergence plot, computed as the maximum error per iteration in the last block [ t N 1 , t N ] , is shown in Figure 3. The convergence of the optimal overlapping (OOHBM) and the original optimal hybrid block technique (OHBM) with no overlapping grid for m = 3 and m = 5 are compared. The graph demonstrates that the overlapping HBM plateaus at a higher degree of accuracy than the non-overlapping HBM. This is one indicator of the overlapping HBM’s higher precision. It is also apparent from Figure 3 that the linear partitioning HBM converges swiftly, with full convergence reached within 10 iterations for all HBM variants evaluated in this scenario.
A comparison of the maximum absolute errors, computed from the exact solution, is displayed in Table 3. It can be seen from the table that the maximum error improves when the overlapping grid is applied. In this example, the improvement in the maximum error when comparing equally spaced and optimal hybrid block methods is not as significant as the difference between the original and proposed overlapping HBM. It can also be observed that the maximum error obtained using the overlapping HBM with equally spaced points is greater than the error of the original optimal HBM, for both m = 3 and m = 5 . This appears to indicate that introducing grid overlap into HBM systems reduces the maximum error more than adjusting the type of intra-step points.
The maximum error in each block [ t n , t n + 1 ] is defined as
E n = max 1 i m | y ( x n + p i ) y n + p i | .
Plotting the profile of the maximum error in each block for n = 0 , 1 , , N 1 reveals the distribution of the error as it traverses the complete domain of integration. Figure 4 depicts the maximum error profile E n for 100 intervals (blocks) used during the deployment of the overlapping optimal HBM and original optimal HBM to Example 31. The graph illustrates the significant reduction in the maximum error when the overlapping grid method is implemented. The maximum error begins in the same location when n = 0 . This is expected given that the initial interval is identical for both optimal overlapping (OOHBMs) and standard optimal non-overlapping HBMs (OHBMs). After n = 0 , the maximum errors of the OOHBM are less than those of the OHBM when m = 3 and when m = 5 . The reduction in the computed error is consistent with the theorems for local truncation errors derived in the preceding section, which demonstrates that the truncation error of the overlapping HBM is O ( h ) more than that of the corresponding non-overlapping HBM.

4.2. Example 2

Consider the logistic equation:
y = κ y ( 1 y ) , y ( 0 ) = y 0
with the exact solution y ( t ) = y 0 ( 1 y 0 ) e κ t + y 0 . The problem was integrated on [ 0 , 10 ] using a step size of h = 0.1 for the case when κ = 2 .
Table 4 compares the maximum errors produced by the OOHBM versus the OHBM. Included also are the maximum errors of the corresponding hybrid block methods constructed using intra-step points equally spaced apart. The table clearly illustrates the increase in accuracy brought about by the introduction of overlapping in grids that use both equally spaced and optimal intra-step points.
Figure 5 provides a comparison of the convergence plots for the standard OHBM and the proposed OOHBM. The illustration demonstrates the speedy convergence of the linear partitioning iterative approach and the superior accuracy of the proposed OOHBM schemes in comparison to the OHBM schemes.
Figure 6 illustrates the maximum error profile E n across all blocks as a function of the block index n. The decrease in error shown in Figure 6 when switching from the OHBM to the OOHBM is consistent with the pattern observed in Figure 4 of the preceding case. This is additional proof that the concept of overlapping grids improves the accuracy of hybrid methods.
Next, we discuss the implementation of the overlapping HBM on systems of nonlinear equations that are notoriously difficult to numerically solve due to their stiff characteristics or their propensity to change rapidly over short time intervals.

4.3. Example 3

Consider the the van der Pol oscillator with nonlinear damping:
y 1 = y 2 , y 1 ( 0 ) = 0.1 , y 2 = μ ( 1 y 1 2 ) y 2 y 1 , y 2 ( 0 ) = 0
where μ is a constant parameter. In this example, we set
ψ 1 ( t , y r , y ˙ r ) = y 2 , r , ϕ 1 ( t , y r , y ˙ r ) = 0 , ψ 2 ( t , y r , y ˙ r ) = y 1 , r + 1 , ϕ 2 ( t , y r , y ˙ r ) = μ ( 1 y 1 , r + 1 2 ) .
This problem was solved using the OOHBM-3 with h = 0.01 in the domain of integration [0,50]. Since there is no analytic solution to the problem, we compared the results of the solution profiles to those that have been previously reported in the literature. In Figure 7, the solution profiles and phase portrait for the van der Pol oscillator when μ = 5 are depicted. Similar results have been reported in the literature using different numerical solution techniques (see, for example, [14]).

4.4. Example 4

Consider the highly stiff system of IVPs below, which has applications in chemical kinetics [15].
y 1 = k 1 y 1 + k 2 y 2 y 3 , y 1 ( 0 ) = 1 , y 2 = k 1 y 1 k 2 y 2 y 3 k 3 y 2 2 , y 2 ( 0 ) = 0 , y 3 = k 3 y 2 2 , y 3 ( 0 ) = 0 ,
For specific values of the reaction constants, this problem becomes highly stiff and unsolvable using ordinary numerical approaches. In this study, we employed the values that were also used to assess the robustness of the extended backward difference methods described in [16] and the implicit Runge–Kutta method as summarized in [17]. Consequently, we set
k 1 = 0.04 , k 2 = 10 4 , k 3 = 3 × 10 7 .
Figure 8 depicts the solution profiles generated by OOHBM-3 with h = 0.001 . The nature of the solution profile for y 2 provides insight into the numerical difficulties of solving this problem. Over a brief period of time, the concentration of y 2 rises abruptly to a maximum and then falls sharply. Compared to the concentrations of y 1 and y 3 represented in Figure 8a, the size of the concentration for y 2 is less. Due to the widely different concentration amounts, the l o g plot on the x-axis is used to compare all the concentration profiles over a longer period of time, as depicted in Figure 8c. These profiles are comparable to those published in [17].

4.5. Example 5

Consider the Lorenz system [18] given by
y 1 = a ( y 2 y 1 ) , y 1 ( 0 ) = 1 y 2 = y 1 y 3 + b y 1 y 2 , y 2 ( 0 ) = 5 y 3 = y 1 y 2 c y 3 , y 3 ( 0 ) = 10
The constants are a = 10 , b = 28 , c = 8 / 3 . With these selected constants, the Lorenz system exhibits chaotic behaviour, with the solution profiles fluctuating rapidly over extremely brief time intervals. To precisely and efficiently capture the solution, robust numerical approaches are necessary. The proposed overlapping HBM with the fewest intra-step points was employed to tackle this problem.
The parameters of the iteration scheme are given below.
ϕ 1 ( t , y 1 , y 2 , y 3 ) = a , ψ 1 ( t , y 1 , y 2 , y 3 ) = a y 2 , ϕ 2 ( t , y 1 , y 2 , y 3 ) = 1 , ψ 2 ( t , y 1 , y 2 , y 3 ) = b y 1 y 1 y 3 , ϕ 3 ( t , y 1 , y 2 , y 3 ) = c , ψ 3 ( t , y 1 , y 2 , y 3 ) = y 1 y 2
In Figure 9 and Figure 10, respectively, the time series solutions and phase portraits for the Lorenz equation are displayed. These results were produced using OOHBM-3 with h = 0.01 across the integrating domain [0,50]. The results are identical to those reported in numerous publications on the Lorenz system’s solutions (see, for example, [19]).

5. Conclusions

In this study, an enhancement to the conventional hybrid block methods for solving first-order equations was developed. This was accomplished by overlapping successive integration blocks after the first block. When executed on a variety of initial value problems, the resulting method, called the overlapping hybrid block method, demonstrated a significant gain in accuracy. In addition, theorems were established that demonstrate an improvement in local truncation error by an order of the step size, O ( h ) . It was also demonstrated that the proposed method is A-stable. These findings show that the concept of inserting overlapping blocks in hybrid block methods may be a profitable field of study. Further research is required to extend the applicability of the overlapping technique to other types of differential equations, such as boundary value problems, fractional-order differential equations, and differential-algebraic systems of equations (DAEs).

Funding

This research was funded by the University of KwaZulu-Natal.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ramos, H. An optimized two-step hybrid block method for solving first-order initial-value problems in ODEs. Differ. Geom.-Dyn. Syst. 2017, 19, 107–118. [Google Scholar]
  2. Kashkari, B.S.; Syam, M.I. Optimization of one step block method with three hybrid points for solving first-order ordinary differential equations. Results Phys. 2019, 12, 592–596. [Google Scholar] [CrossRef]
  3. Muideen, O.; Yahaya, H.; Raphael, A.; Morufu, O. Optimized Three-step Hybrid Block Method for Stiff Problems in Ordinary Differential Equations. J. Sci. Eng. 2020, 17, 80–95. [Google Scholar]
  4. Singh, G.; Garg, A.; Singla, R.; Kanwar, V. A novel two-parameter class of optimized hybrid block methods for integrating differential systems numerically. Comput. Math. Methods 2021, 3, e1214. [Google Scholar] [CrossRef]
  5. Singla, R.; Singh, G.; Ramos, H.; Kanwar, V. A Family of A-Stable Optimized Hybrid Block Methods for Integrating Stiff Differential Systems. Math. Probl. Eng. 2022, 2022, 5576891. [Google Scholar] [CrossRef]
  6. Singh, G.; Garg, A.; Kanwar, V.; Ramos, H. An efficient optimized adaptive step-size hybrid block method for integrating differential systems. Appl. Math. Comput. 2019, 362, 124–567. [Google Scholar] [CrossRef]
  7. Ramos, H.; Rufai, M.A. An adaptive pair of one-step hybrid block Nyström methods for singular initial-value problems of Lane–Emden–Fowler type. Math. Comput. Simul. 2022, 193, 497–508. [Google Scholar] [CrossRef]
  8. Ramos, H.; Kaur, A.; Kanwar, V. Using a cubic B-spline method in conjunction with a one-step optimized hybrid block approach to solve nonlinear partial differential equations. Comput. Appl. Math. 2022, 41, 34. [Google Scholar] [CrossRef]
  9. Kaur, A.; Kanwar, V. Numerical Solution of Generalized Kuramoto–Sivashinsky Equation Using Cubic Trigonometric B-Spline Based Differential Quadrature Method and One-Step Optimized Hybrid Block Method. Int. J. Appl. Comput. Math. Vol. 2022, 8, 20. [Google Scholar] [CrossRef]
  10. Butcher, J.C. Implicit runge-kutta processes. Math. Comput. 1964, 18, 50–64. [Google Scholar] [CrossRef]
  11. Orakwelu, M.G. Generalised Implicit Block Hybrid Algorithms for Initial Value Problems. Ph.D. Thesis, University of KwaZulu-Natal, Durban, South Africa, 2019. [Google Scholar]
  12. Pasetto, M.; Waisman, H.; Chen, J.S. A waveform relaxation Newmark method for structural dynamics problems. Comput. Mech. 2019, 63, 1223–1242. [Google Scholar] [CrossRef]
  13. Barbara, Z.K.; Vandewalle, S. Waveform Relaxation for Functional-Differential Equations. SIAM J. Sci. Comput. 1999, 21, 207–226. [Google Scholar]
  14. Trefethen, L.N.; Birkisson, Á.; Driscoll, T.A. Exploring ODEs; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2017. [Google Scholar]
  15. Robertson, H. The solution of a set of reaction rate equations. In Numerical Analysis: An Introduction; Academic Press: Cambridge, MA, USA, 1966; pp. 178–182. [Google Scholar]
  16. Cash, J. Second derivative extended backward differentiation formulas for the numerical integration of stiff systems. SIAM J. Numer. Anal. 1981, 18, 21–36. [Google Scholar] [CrossRef]
  17. Butcher, J.C. Numerical Methods for Ordinary Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  18. Lorenz, E.N. Deterministic nonperiodic flow. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  19. Hirsch, M.W.; Smale, S.; Devaney, R.L. Differential Equations, Dynamical Systems, and an Introduction to Chaos; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar]
Figure 1. Illustration of the overlapping grid when N = 4 and m = 3 .
Figure 1. Illustration of the overlapping grid when N = 4 and m = 3 .
Algorithms 15 00427 g001
Figure 2. Stability regions for the overlapping hybrid block method. (a) Equally spaced intra-step points. (b) Optimal intra-step points.
Figure 2. Stability regions for the overlapping hybrid block method. (a) Equally spaced intra-step points. (b) Optimal intra-step points.
Algorithms 15 00427 g002
Figure 3. Comparison of the convergence plots for Example 1 when h = 0.01 .
Figure 3. Comparison of the convergence plots for Example 1 when h = 0.01 .
Algorithms 15 00427 g003
Figure 4. Maximum error profile for Example 1.
Figure 4. Maximum error profile for Example 1.
Algorithms 15 00427 g004
Figure 5. Comparison of the convergence plots for Example 2 when h = 0.1 .
Figure 5. Comparison of the convergence plots for Example 2 when h = 0.1 .
Algorithms 15 00427 g005
Figure 6. Maximum error profile for Example 2.
Figure 6. Maximum error profile for Example 2.
Algorithms 15 00427 g006
Figure 7. Solution profile of the van der Pol oscillator example when h = 0.01 . (a) y 1 ( t ) vs. t. (b) y 2 ( t ) vs. t. (c) Phase portrait.
Figure 7. Solution profile of the van der Pol oscillator example when h = 0.01 . (a) y 1 ( t ) vs. t. (b) y 2 ( t ) vs. t. (c) Phase portrait.
Algorithms 15 00427 g007
Figure 8. Solution profiles of Example 4 with m = 3 and h = 0.001 . (a) concentrations of y 1 and y 3 (b) concentrations of y 2 (c) concentration y 1 , y 2 and y 3 .
Figure 8. Solution profiles of Example 4 with m = 3 and h = 0.001 . (a) concentrations of y 1 and y 3 (b) concentrations of y 2 (c) concentration y 1 , y 2 and y 3 .
Algorithms 15 00427 g008
Figure 9. Time series profiles of the Lorenz system. (a) y 1 ; (b) y 2 ; (c) y 3 .
Figure 9. Time series profiles of the Lorenz system. (a) y 1 ; (b) y 2 ; (c) y 3 .
Algorithms 15 00427 g009
Figure 10. Phase portraits of the Lorenz system. (a) y 1 and y 2 ; (b) y 1 and y 3 ; (c) y 2 and y 3 .
Figure 10. Phase portraits of the Lorenz system. (a) y 1 and y 2 ; (b) y 1 and y 3 ; (c) y 2 and y 3 .
Algorithms 15 00427 g010
Table 1. Truncation errors for the overlapping and non-overlapping HBM with equally spaced and optimal intra-step points.
Table 1. Truncation errors for the overlapping and non-overlapping HBM with equally spaced and optimal intra-step points.
Equally SpacedOptimal Points
L i Non-OverlappingOverlappingNon-OverlappingOverlapping
p 1 19 y ( 5 ) ( t n ) h 5 174960 11 y ( 6 ) t n h 6 1049760 y ( 5 ) ( t n ) h 5 6000 5 4 3 5 y ( 6 ) t n h 6 450000
m = 3 p 2 y ( 5 ) ( t n ) h 5 21870 y ( 7 ) t n h 7 1653372 y ( 5 ) ( t n ) h 5 6000 5 1 + 3 5 y ( 6 ) t n h 6 450000
p 3 y ( 5 ) ( t n ) h 5 6480 y ( 6 ) t n h 6 38880 y ( 7 ) ( t n ) h 7 1512000 y ( 7 ) t n h 7 1512000
p 1 3 y ( 6 ) ( t n ) h 6 655360 271 y ( 7 ) t n h 7 990904320 y ( 6 ) ( t n ) h 6 493920 49 8 21 y ( 7 ) t n h 7 145212480
m = 4 p 2 y ( 6 ) ( t n ) h 6 368640 y ( 7 ) t n h 7 12386304 y ( 6 ) ( t n ) h 6 322560 14 + 21 y ( 7 ) t n h 7 27095040
p 3 3 y ( 6 ) ( t n ) h 6 655360 13 y ( 7 ) t n h 7 36700160 y ( 6 ) ( t n ) h 6 493920 49 + 21 y ( 7 ) t n h 7 145212480
p 4 y ( 7 ) ( t n ) h 7 1935360 y ( 7 ) t n h 7 1935360 y ( 9 ) ( t n ) h 9 1422489600 y ( 9 ) t n h 9 1422489600
p 1 863 y ( 7 ) ( t n ) h 7 4725000000 13 y ( 8 ) t n h 8 1750000000 147 30 7 y ( 7 ) ( t n ) h 7 160030080 12 147 30 7 + 26 7 + 13 y ( 8 ) t n h 8 13442526720
p 2 37 y ( 7 ) ( t n ) h 7 295312500 y ( 8 ) t n h 8 295312500 147 + 30 7 y ( 7 ) ( t n ) h 7 160030080 6.05861 × 10 9 y ( 8 ) ( t n ) h 8
m = 5 p 3 29 y ( 7 ) ( t n ) h 7 175000000 13 y ( 8 ) t n h 8 1750000000 147 + 30 7 y ( 7 ) ( t n ) h 7 160030080 1.05281 × 10 8 y ( 8 ) ) h 8
p 4 8 y ( 7 ) ( t n ) h 7 73828125 13 y ( 9 ) t n h 9 27685546875 147 30 7 y ( 7 ) ( t n ) h 7 160030080 12 147 30 7 + 2 7 + 1 y ( 8 ) t n h 8 13442526720
p 5 11 y ( 7 ) ( t n ) h 7 37800000 11 y ( 8 ) t n h 8 378000000 y ( 11 ) h 11 2112397056000 y ( 11 ) t n h 11 2112397056000
Table 2. Intra-step points used in numerical experimentation.
Table 2. Intra-step points used in numerical experimentation.
mEquidistantOptimal Points
3 1 3 , 2 3 1 10 ( 5 ± 5 )
5 1 5 , 2 5 , 3 5 , 4 5 1 2 1 ± 1 21 7 ± 2 7
Table 3. Maximum errors for Example 1 using different choices of intra-step points.
Table 3. Maximum errors for Example 1 using different choices of intra-step points.
m = 3 m = 5
Original HBMProposed HBMOriginal HBMProposed HBM
Equally spaced points 9.9970 × 10 7 6.7542 × 10 8 6.0825 × 10 10 4.2735 × 10 11
Optimal points 1.8106 × 10 7 3.2112 × 10 8 1.4701 × 10 10 1.5855 × 10 11
Table 4. Maximum errors for Example 2 using different choices of intra-step points.
Table 4. Maximum errors for Example 2 using different choices of intra-step points.
m = 3 m = 5
Original HBMProposed HBMOriginal HBMProposed HBM
Equally spaced points 1.3918 × 10 7 1.5292 × 10 8 7.2577 × 10 12 8.5256 × 10 13
Optimal points 6.3022 × 10 9 4.4041 × 10 10 1.2751 × 10 12 6.2343 × 10 14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Motsa, S. Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems. Algorithms 2022, 15, 427. https://doi.org/10.3390/a15110427

AMA Style

Motsa S. Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems. Algorithms. 2022; 15(11):427. https://doi.org/10.3390/a15110427

Chicago/Turabian Style

Motsa, Sandile. 2022. "Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems" Algorithms 15, no. 11: 427. https://doi.org/10.3390/a15110427

APA Style

Motsa, S. (2022). Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems. Algorithms, 15(11), 427. https://doi.org/10.3390/a15110427

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop