Fitted Parameter Exponential Spline Method for Singularly Perturbed Delay Differential Equations with a Large Delay

Prasad, E. Siva; Omkar, R.; Phaneendra, Kolloju

doi:https://doi.org/10.1155/2022/9291834

Computational and Mathematical Methods

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 9291834 | https://doi.org/10.1155/2022/9291834

Fitted Parameter Exponential Spline Method for Singularly Perturbed Delay Differential Equations with a Large Delay

E. Siva Prasad,¹R. Omkar,¹and Kolloju Phaneendra¹

Academic Editor: Higinio Ramos

Received23 Mar 2022

Revised09 Jun 2022

Accepted26 Jun 2022

Published16 Aug 2022

Abstract

In this paper, we present a new computational method based on an exponential spline for solving a class of delay differential equations with a negative shift in the differentiated term. When the shift parameter is , the proposed method works well and also controls the oscillations in the solution’s layer region. To accomplish this, we included a parameter in the proposed numerical scheme that is based on a special type of mesh, and the parameter is evaluated using the theory of singular perturbation. Maximum absolute errors and convergences of numerical solutions are tabulated to demonstrate the efficiency of the proposed computational method and to support the convergence analysis of the presented scheme.

1. Introduction

In this article, we examine a section of differential equations with a delay in the first-order derivative term and a small parameter multiplied by the second-order derivative. These equations are widely used in the simulation of a wide range of physical applications, such as construction of logical circuits by means of a network built up by 1D excitable media [1], hybrid optical systems with delayed feedback [2], population dynamics [3], nonlinear delay differential equations relating to physiological control systems [4], red blood cell system [5], predator-prey models [6], a Hamiltonian boundary value problem associated with optimal control theory [7], and neuronal variability problems connected to patterns of nerve action potentials formed by unit quantal inputs occurring at random which are studied [8].

Bellman et al. [9], Doolan et al. [10], Driver [11], Norkin [12], Kokotovic et al. [13], Miller et al. [14], Smith [15], and O’Malley [16] are just a few of the books in the collection that can be used for further investigations of precise characteristics of the aforementioned class of models and perturbation problems. Lange and Miura [17, 18] devised an asymptotic method for solving equations with layer behaviour. Researchers presented a mathematical model in [17] for predicting the time it takes nerve cells with arbitrary synaptic inputs in dendrites to produce potential action. The researchers in [18] demonstrated the issues by using solutions that exhibit rapid oscillations. Kadalbajoo and his colleagues began widespread numerical work using finite differentiation techniques as well as the B-spline technique with fitted mesh [19–22]. The researchers in [23] suggested a numerical approach via nonpolynomial spline for the solution of differential-difference equations with small and large delay.

In [24], researchers proposed an unconditionally stable collocation method for solving singularly perturbed one-dimensional time dependent problems. The method employs the implicit finite difference method on a uniform mesh in the temporal direction and the B-spline collocation method on a nonuniform Shishkin mesh in the spatial direction. Kadalbajoo and Gupta [25] discussed a B-spline collocation approach on a piecewise uniform Shishkin mesh for the numerical solution of a singularly perturbed problem with a boundary layer on the right side of the domain. Furthermore, it is demonstrated that the approach is uniformly convergent in the maximum norm of the second order. A brief survey of various computational techniques on different type of singularly perturbed problems is given in [26]. Gupta and Kadalbajoo [27] proposed a parameter uniform method for solving a singularly perturbed one-dimensional parabolic equation with multiple boundary turning points on a rectangular domain that combines the Euler method for time variables on a uniform mesh and the B-spline collocation method for space variables on a Shishkin mesh. The authors of [28] developed a numerical scheme for time-dependent singularly perturbed problems based on Mickens nonstandard finite difference method. The method is shown to be unconditionally stable and parameter uniformly convergent for both large and small shift arguments. Das [29] presented an a posteriori-based convergence analysis on an adaptive mesh for a nonlinear singularly perturbed system of delay differential equations. The layer-adaptive solution on the a posteriori-generated mesh has been shown to converge uniformly to the exact solution with optimal order accuracy.

The authors of [30] considered a moving mesh-adaptive algorithm that adapts meshes to boundary layers for the numerical solution of singularly perturbed convection-diffusion-reaction problems with two small parameters. A special positive monitor function is used to generate the mesh and demonstrate that the proposed method is parameter independent and convergent for a class of model problems. Das and Natesan [31] used cubic spline approximation on a piecewise uniform Shishkin mesh to solve a system of reaction diffusion differential equations for Robin type boundary-value problems. This method is a second-order uniform convergence method that uses a central difference scheme for the outer region of the boundary layers and a cubic spline approximation for the inner region. The authors of [32] presented two adaptive methods for a system of reaction diffusion problems using the -refinement strategy, which moves the mesh points toward the boundary layers. On the adaptively generated equidistributed mesh generated by a curvature-based monitor function, first-order uniform solutions are obtained using forward and backward schemes. To improve the discrete solution to second-order uniform accuracy, a cubic spline-based scheme is used. The authors of [33] investigated a numerical scheme for solving singularly perturbed parabolic partial differential equations using the Euler method and the upwind scheme on an adaptive mesh. They also examined the convergence analysis of the proposed scheme. The authors of [34] investigated the approximation of solutions to a class of fractional order Volterra-Fredholm integro-differential equations. The detailed analysis of the existence and uniqueness of the solution as well as the solution’s bounds was addressed. A perturbation approach based on homotopy analysis is used to solve the same problem. The authors of [35] considered the initial and boundary value problems of a class of fractional order Volterra integro-differential equations of the first kind. Using Leibniz’s rule, the same problem is reduced to a Volterra integro-equation of the second kind, and an operative-based method is used to approximate their solutions. Das et al. [36] investigate the homotopy perturbation method for approximating solutions of Volterra-Fredholm integro-differential equations using semianalytical approaches.

The following is the paper’s outline: the problem is described in Section 2 along with the conditions for the layer’s behaviour, and the exponential spline function is defined in Section 3. Section 4 discusses the numerical scheme with a large delay. Section 5 consists of truncation error and convergence analysis. Numerical experiments are addressed in Section 6. Finally, in Section 7, the article’s discussions and conclusions are addressed.

2. Statement of the Problem

Consider a singularly perturbed delay differential equation with a delay term, that is, a negative shift in the first-order derivative in accordance with the conditions. where are expected to be smooth, bounded functions in [0, 1], is the delay parameter, and is a finite constant. When , Equation (1) is reduced to a singularly perturbed equation with boundary layer behaviour and turning points that depend on the sign of . Throughout the interval [0, 1], the solution has a boundary layer on the left end side when is positive and on the right end side if is negative. When the delay parameter is of , the behaviour of the layer can change and even be ruined, or the solution exhibits oscillatory behaviour. To monitor these oscillations in the solution profile, we devise a numerical method that combined a special mesh introduced in [21] with the exponential spline method’s fitting parameter.

2.1. Properties of Continuous Problem

Let be the differential operator for the problem Equation (1) which is defined for any smooth function as .

Lemma 1 (Continuous maximum principle). Let be the smooth function satisfying and Then, implies that .

Proof. Let be such that and assume that . Clearly, ; therefore, and. Now, in order to prove , we consider different cases.
Case 1:
Case 2:
Combining the above two cases, we get which contradicts the hypothesis that Hence, our assumption that is wrong, and thus,

Lemma 2. Under the assumption that , where , are positive constants, the solution of Equation (1) with boundary conditions Equation (2) exists and satisfies

Proof. Let us construct the two-barrier function defined by Then, we have

Case 1. For , since , i.e., . Since , and we choose the constant so that the sum of the first two terms dominates the third term in the above inequality (7), we then obtain

Case 2. For , since , i.e., . From the inequality (8) and (10), we get
Therefore, using Lemma 1, we obtain , for all , which gives required estimate.
Lemma 1 implies that the solution is unique, and since the problem under consideration is linear, the existence of the solution is implied by its uniqueness. Additionally, Lemma 2 gives the boundedness of the solution [22].

Lemma 3. The bounds on the derivatives of the solution of Equation (1) with respect to is provided by

Proof. One can refer to [28].

3. Exponential Spline

We divide [0, 1] into equal subdomains of grid size , so that with . Let be the precise solution and be an estimate to by the exponential spline passing through the points and . The exponential spline function has the following form for each ^th segment, satisfying the condition of first derivative continuity at the joint nodes . where , and are constants and is a parameter. Here, is of the class interpolates at mesh points depends on a parameter and reduces to the cubic spline in [0,1] as

Define , to calculate the coefficient values in Equation (12) in terms of , and Following the exercise in algebra, we get

where , for .

The continuity condition of the first derivative, i.e., at , allows us to obtain the following relations for where .

The local truncation error of the scheme in Equation (14) is for .

As a result, we obtain different orders for for different values of and in Equation (14): (i) fourth order for any arbitrary and with ; (ii) six orders for .

4. Numerical Approach with Fitting Parameter for Large Delay

The solution’s layer behaviour is preserved, and precise results are obtained when compared to the perturbation parameter; the shift parameter is smaller [22, 23]. However, if is of order , the solution’s layer behaviour is no longer well-preserved, and oscillations becomes visible. As a result, we are developing a numerical scheme with fitting parameters that is based on an exponential spline method and a special type of mesh developed in [21]. In addition, we examine the solution layer behaviour graphically for large delays in order to demonstrate the significance of the fitting parameter in our proposed scheme.

4.1. Numerical Scheme for Left-End Boundary Layer

Assume that and in [0, 1], where is a positive constant. As a result of this theory, Equations (1) and (2) show layer behaviour at for small values of . Now, we consider a mesh developed in [21], i.e., chooses the mesh as , where , is the mantissa of , and is a positive integer. The discretization of the boundary value problem Equations (1) and (2) with a special mesh result in

with

Using the difference scheme in Equation (14) and the difference estimates of the first-order derivatives listed below, we get

Now, to enhance the scheme’s effectiveness, we incorporate a fitting parameter, which is determined using the theory of singular perturbations.

By multiplying Equation (19) by and assuming limit as [4], we get

Based on singular perturbation theory pertaining to the layer at the interval’s left end [16], we have

Now, by plugging Equation (21) into Equation (20) and exercising, we have which is the parameter that has been fitted for the layer on the left-end of the domain.

4.2. Numerical Scheme for Right-End Boundary Layer

Assume that and in [0, 1], where is a positive constant. As a result of this theory, Equations (1) and (2) shows layer behaviour at for small values of . Based on singular perturbation theory for the layer at the right end of the interval [16], we have

In the exponential spline difference scheme Equation (18), we insert a fitting parameter and obtain the fitting parameter employing the similar procedure as in the case of the left boundary layer. Now, by simplifying the difference scheme in Equation (19), we get where

Following Equation (16), the scheme of Equation (27) can be expressed as

The Gauss elimination scheme with partial pivoting is used to solve the aforementioned system of equations.

5. Truncation Error

The local truncation error for the scheme Equation (27) is obtained using Taylor’s series expansion, as shown below. where . Clearly, for any random choice of and whose sum is equal to and .

5.1. Convergence Analysis

The matrix equation for the left-end boundary layer problem, including the specified boundary conditions in Equation (24), is as follows: and , where .

and , are related vectors of Equation (19).

Let satisfies the equation

Let be the error so that .

Subtracting Equation (30) from Equation (33), we obtain the error equation

Let be the sum of the row elements of the matrix . Then, we have

Let , , and .

Since and with sufficiently small it is verified that is irreducible and monotone [37, 38]. Therefore,

So, we can deduce from Equation (34) that

For small , we have .

Let , and we define , since.

Hence,

Furthermore,

Using Equations (36), (38), and (39), we get

As a result, it demonstrates that the proposed exponential spline scheme is second-order convergent. The same procedure can also be used to examine the convergence analysis for the right end boundary layer.

6. Numerical Experiments

Three problems are selected to validate the competence of the proposed scheme. The maximum absolute errors are obtained using the principle of double mesh and the order of convergence for the examples, which are tabulated in Tables 1–3. Furthermore, the calculated arithmetic solutions to the considered experiments are depicted by graphs, with and without fitting parameters for various values of and

Example 1. .

Example 2.

Example 3. .

7. Discussions and Conclusion

The delay differential equation of order two is considered, which has a large delay in the convention term and layers at the left and right ends of the interval. An exponential spline finite difference scheme is constructed using the continuity of its first-order derivative condition at the joint nodes. We considered the case when is of order , layer’s behaviour can change and even be destroyed, or the solution can exhibit oscillatory behaviour. To address these drawbacks in solutions, we experimented with a new scheme based on the special mesh introduced in [21], which involves adding a parameter to the exponential spline method. Tables 1–3 illustrate the maximum errors as well as the order of convergence for Examples 1–3. Table 4 compares the maximum absolute errors in Example 1 to other methods for different shift parameter values. Figures 1–6 depict diagrams of the solutions to the experiments for various values of , and we compared the graphs with the fitting parameter for various values of to the graphs without the fitting parameter. When the value is larger than the , the solution includes oscillations, as shown in Figures 2 and 4 (without fitting parameter), whereas in Figures 1 and 3 (with parameter), the oscillations are controlled in solutions, and behaviour of the layer is well-preserved. However, when dealing with a right-end layer, as shown in Figures 5 and 6, layer behaviour of the solution is preserved even for Finally, we found that exponential spline fitting parameter method with special mesh had a significant improvement in terms of regulating the oscillations in the solutions of delay differential equations.

Data Availability

Data is available in the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

With the motivation of the numerical scheme derived in the manuscript titled “Numerical Approach for Differential-Difference Equations having Layer Behaviour with Small or Large Delay using Non-Polynomial Spline,” https://assets.researchsquare.com/files/rs-261904/v1_covered.pdf (Ref. [23]), we have derived a scheme for the problem having large delay in this manuscript. We thank the reviewers for their constructive comments/suggestions to strengthen the manuscript.

References

M. Bestehornand York and E. V. Grigorieva, “Formation and propagation of localized states in extended systems,” Annalen der Physik, vol. 13, no. 78, pp. 423–431, 2004.
View at: Publisher Site | Google Scholar
M. W. Derstine, H. M. Gibbs, F. A. Hopf, and D. L. Kaplan, “Bifurcation gap in a hybrid optically bistable system,” Physical Review A, vol. 26, no. 6, pp. 3720–3722, 1982.
View at: Publisher Site | Google Scholar
Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Academic Press, New York, 1993.
M. C. Mackey and L. Glass, “Oscillation and chaos in physiological control systems,” Science, vol. 197, no. 4300, pp. 287–289, 1977.
View at: Publisher Site | Google Scholar
A. Lasota and M. Walewska, “Mathematical models of the red blood cell system,” Roczniki Polskiego Towarzystwa Matematycznego, vol. 6, pp. 25–40, 1976.
View at: Google Scholar
A. Martin and S. Raun, “Predator-prey models with delay and prey harvesting,” Journal of Mathematical Biology, vol. 43, no. 3, pp. 247–267, 2001.
View at: Publisher Site | Google Scholar
V. Y. Glizer, “Asymptotic solution of a boundary-value problem for linear singularly-perturbed functional differential equations arising in optimal control theory,” Journal of Optimization Theory and Applications, vol. 106, no. 2, pp. 309–335, 2000.
View at: Publisher Site | Google Scholar
R. B. Stein, “Some models of neuronal variability,” Biophysical Journal, vol. 7, no. 1, pp. 37–68, 1967.
View at: Publisher Site | Google Scholar
R. Bellman, K. L. Cooke, R. Bellman, and J. Gillis, Differential-Difference Equations, Academic Press, New York, vol. 16, no. 12, 1963.
View at: Publisher Site
E. P. Doolan, J. J. H. Miller, and W. H. A. Schilders, Uniform Numerical Methods for Problems with Initial and Boundary Layers, Boole Press, Dublin, 1980.
R. D. Driver, Ordinary and Delay Differential Equations, Springer, Belin-Heidelberg, New York, 1977.
View at: Publisher Site
S. B. Norkin, Introduction to the Theory and Application of Differential Equations with Deviating Arguments, Academic Press, 1973.
P. V. Kokotovic, H. K. Khalil, and J. O’Reilly, Singular Perturbation Methods in Control Analysis and Design, Academic Press, New York, 1986.
J. J. Miller, E. O'Riordan, G. I. Shishkin, and R. B. Kellogg, Fitted Numerical Methods for Singular Perturbation Problems, World Scientific, Singapore, 1996.
View at: Publisher Site
H. Smith, An Introduction to Delay Differential Equations with Applications to the Life Sciences, Springer, Berlin, 2011.
R. E. O’Malley, Introduction to Singular Perturbations, Academic Press, New York, 1974.
C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary-value problems for differential-difference equations. v. Small shifts with layer behavior,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 249–272, 1994.
View at: Publisher Site | Google Scholar
C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary-value problems for differential-difference equations. vi. Small shifts with rapid oscillations,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 273–283, 1994.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of singularly perturbed delay differential equations with layer behavior,” Applied Mathematics and Computation, vol. 157, no. 1, pp. 11–28, 2004.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and K. K. Sharma, “<![CDATA[<B>Numerical treatment of boundary value problems for second order singularly perturbed delay differential equations</B>]]>,” Journal of Computational and Applied Mathematics, vol. 24, no. 2, pp. 151–172, 2005.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and K. K. Sharma, “A numerical method based on finite difference for boundary value problems for singularly perturbed delay differential equations,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 692–707, 2008.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and D. Kumar, “Fitted mesh _B_ -spline collocation method for singularly perturbed differential-difference equations with small delay,” Applied Mathematics and Computation, vol. 204, no. 1, pp. 90–98, 2008.
View at: Publisher Site | Google Scholar
M. Lalu, K. Phaneendra, and S. P. Emineni, “Numerical approach for differential-difference equations having layer behaviour with small or large delay using non-polynomial spline,” Soft Computing, vol. 25, no. 21, pp. 13709–13722, 2021.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo, V. Gupta, and A. Awasthi, “A uniformly convergent B-spline collocation method on a nonuniform mesh for singularly perturbed one-dimensional time-dependent linear convection-diffusion problem,” Journal of Applied Mathematics and Computing, vol. 220, no. 1-2, pp. 271–289, 2008.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and V. Gupta, “Numerical solution of singularly perturbed convection-diffusion problem using parameter uniform B-spline collocation method,” Journal of Mathematical Analysis and Applications, vol. 355, no. 1, pp. 439–452, 2009.
View at: Publisher Site | Google Scholar
M. K. Kadalbajoo and V. Gupta, “A brief survey on numerical methods for solving singularly perturbed problems,” Applied Mathematics and Computation, vol. 217, no. 8, pp. 3641–3716, 2010.
View at: Publisher Site | Google Scholar
V. Gupta and M. K. Kadalbajoo, “A layer adaptive B-spline collocation method for singularly perturbed one-dimensional parabolic problem with a boundary turning point,” Numerical Methods for Partial Differential Equations, vol. 27, no. 5, pp. 1143–1164, 2011.
View at: Publisher Site | Google Scholar
K. Bansal and K. K. Sharma, “Parameter uniform numerical scheme for time dependent singularly perturbed convection-diffusion-reaction problems with general shift arguments,” Numerical Algorithms, vol. 75, no. 1, pp. 113–145, 2017.
View at: Publisher Site | Google Scholar
P. Das, “An a posteriori-based convergence analysis for a nonlinear singularly perturbed system of delay differential equations on an adaptive mesh,” Numerical Algorithms, vol. 81, no. 2, pp. 465–487, 2019.
View at: Publisher Site | Google Scholar
P. Das and V. Mehrmann, “Numerical solution of singularly perturbed convection- diffusion-reaction problems with two small parameters,” BIT Numerical Mathematics, vol. 56, no. 1, pp. 51–76, 2015.
View at: Publisher Site | Google Scholar
P. Das and S. Natesan, “A uniformly convergent hybrid scheme for singularly perturbed system of reaction-diffusion Robin type boundary-value problems,” Journal of Applied Mathematics and Computing, vol. 41, no. 1-2, pp. 447–471, 2013.
View at: Publisher Site | Google Scholar
P. Das, S. Rana, and J. Vigo-Aguiar, “Higher order accurate approximations on equidistributed meshes for boundary layer originated mixed type reaction diffusion systems with multiple scale nature,” Applied Numerical Mathematics, vol. 148, pp. 79–97, 2020.
View at: Publisher Site | Google Scholar
P. Das, “A higher order difference method for singularly perturbed parabolic partial differential equations,” Journal of Difference Equations and Applications, vol. 24, no. 3, pp. 452–477, 2018.
View at: Publisher Site | Google Scholar
P. Das, S. Rana, and H. Ramos, “A perturbation-based approach for solving fractional-order Volterra-Fredholm integro differential equations and its convergence analysis,” International Journal of Computer Mathematics, vol. 97, no. 10, pp. 1994–2014, 2020.
View at: Publisher Site | Google Scholar
P. Das, S. Rana, and H. Ramos, “On the approximate solutions of a class of fractional order nonlinear Volterra integro-differential initial value problems and boundary value problems of first kind and their convergence analysis,” Journal of Computational and Applied Mathematics, vol. 404, article 113116, 2022.
View at: Publisher Site | Google Scholar
P. Das, S. Rana, and H. Ramos, “Homotopy perturbation method for solving Caputo-type fractional-order Volterra-Fredholm integro-differential equations,” Computational and Mathematical Methods, vol. 1, no. 5, pp. 1–9, 2019.
View at: Publisher Site | Google Scholar
R. S. Varga, Matrix Iterative Analysis, Prentice-Hall, Englewood Cliffs, New Jersey, 1962.
D. M. Young, Iterative Solutions of Large Linear Systems, Academic press, New York, 1971.

Copyright

Copyright © 2022 E. Siva Prasad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

428

Downloads

419

Citations