
This article focuses on a class of fourth-order singularly perturbed convection diffusion equations (SPCDE) with integral boundary conditions (IBC). A numerical method based on a finite difference scheme using Shishkin mesh is presented. The proposed method is close to the first-order convergent. The discrete norm yields an error estimate and theoretical estimations are tested by numerical experiments.
Citation: V. Raja, E. Sekar, S. Shanmuga Priya, B. Unyong. Fitted mesh method for singularly perturbed fourth order differential equation of convection diffusion type with integral boundary condition[J]. AIMS Mathematics, 2023, 8(7): 16691-16707. doi: 10.3934/math.2023853
[1] | Amjid Ali, Teruya Minamoto, Rasool Shah, Kamsing Nonlaopon . A novel numerical method for solution of fractional partial differential equations involving the ψ-Caputo fractional derivative. AIMS Mathematics, 2023, 8(1): 2137-2153. doi: 10.3934/math.2023110 |
[2] | Muhammad Aslam, Muhammad Farman, Hijaz Ahmad, Tuan Nguyen Gia, Aqeel Ahmad, Sameh Askar . Fractal fractional derivative on chemistry kinetics hires problem. AIMS Mathematics, 2022, 7(1): 1155-1184. doi: 10.3934/math.2022068 |
[3] | Muhammad Bilal Riaz, Nauman Raza, Jan Martinovic, Abu Bakar, Osman Tunç . Modeling and simulations for the mitigation of atmospheric carbon dioxide through forest management programs. AIMS Mathematics, 2024, 9(8): 22712-22742. doi: 10.3934/math.20241107 |
[4] | Shabir Ahmad, Aman Ullah, Mohammad Partohaghighi, Sayed Saifullah, Ali Akgül, Fahd Jarad . Oscillatory and complex behaviour of Caputo-Fabrizio fractional order HIV-1 infection model. AIMS Mathematics, 2022, 7(3): 4778-4792. doi: 10.3934/math.2022265 |
[5] | Hashem Najafi, Abdallah Bensayah, Brahim Tellab, Sina Etemad, Sotiris K. Ntouyas, Shahram Rezapour, Jessada Tariboon . Approximate numerical algorithms and artificial neural networks for analyzing a fractal-fractional mathematical model. AIMS Mathematics, 2023, 8(12): 28280-28307. doi: 10.3934/math.20231447 |
[6] | Amnah E. Shammaky, Eslam M. Youssef . Analytical and numerical techniques for solving a fractional integro-differential equation in complex space. AIMS Mathematics, 2024, 9(11): 32138-32156. doi: 10.3934/math.20241543 |
[7] | Mohamed Akel, Muajebah Hidan, Salah Boulaaras, Mohamed Abdalla . On the solutions of certain fractional kinetic matrix equations involving Hadamard fractional integrals. AIMS Mathematics, 2022, 7(8): 15520-15531. doi: 10.3934/math.2022850 |
[8] | Kamal Shah, Hafsa Naz, Muhammad Sarwar, Thabet Abdeljawad . On spectral numerical method for variable-order partial differential equations. AIMS Mathematics, 2022, 7(6): 10422-10438. doi: 10.3934/math.2022581 |
[9] | Kottakkaran Sooppy Nisar, G Ranjith Kumar, K Ramesh . The study on the complex nature of a predator-prey model with fractional-order derivatives incorporating refuge and nonlinear prey harvesting. AIMS Mathematics, 2024, 9(5): 13492-13507. doi: 10.3934/math.2024657 |
[10] | Shao-Wen Yao, Muhammad Farman, Maryam Amin, Mustafa Inc, Ali Akgül, Aqeel Ahmad . Fractional order COVID-19 model with transmission rout infected through environment. AIMS Mathematics, 2022, 7(4): 5156-5174. doi: 10.3934/math.2022288 |
This article focuses on a class of fourth-order singularly perturbed convection diffusion equations (SPCDE) with integral boundary conditions (IBC). A numerical method based on a finite difference scheme using Shishkin mesh is presented. The proposed method is close to the first-order convergent. The discrete norm yields an error estimate and theoretical estimations are tested by numerical experiments.
Mathematical modelling is the best way to formulating problems from an application area and it is well known that several mathematical characterization of numerous growth in chemical and physical sciences is described by differential equations (DEs). In chemistry, chemical kinetics problem and CO2 with PGE problems are described by system of nonlinear (DEs) with different kind of Neumann boundary and Drichlet type conditions in different published work such as chemistry problem by Jawary and Raham [1], chemistry problem by Abbasbandy and Shirzadi [2], CO2 absorbed into PGE problem by Jawary et al. [3], Choe et al. [4], Singha et al. [5], CO2 absorbed into PGE problem by Robertson [6], chemistry problem by Matinfar et al. [7], chemistry problem by Ganji et al. [8], Dokoumetzidis et al. [9]. In the past few years, fractional calculus (FC) has found many diverse and robust applications in various research areas such as fluid dynamics, image processing, viscoelasticity and other physical phenomena. Many definitions of fractional derivatives are discovered by several mathematicians but two most famous definitions of fractional derivative are Riemann − Liouville and Caputo. Some interesting and fundamental works on various direction of the FC is given in several famous books such as by Mainardi [10], fractional differential equations by Podlubny [11], Diethelm [12], Kilbas et al. [13] and Das [14].
In the past few years, wavelets have become an increasingly newly developed famous mechanism in the several research areas of physical, chemical, computational sciences, Image manipulation, signal analysis, data compression, numerical analysis and several others research areas such as a primer on wavelets and their scinentific applications by Walker [15], wavelet: mathematics and applications by Benedetto [16], a mathematical tool for signal analysis by Chui [17], wavelet methods for dynamical problems by Gopalakrishnan and Mitra [18], Wang [19] and a wavelet operational matrix method by Wu [20]. Due to this reason, wavelets have been applied for the solution of differential equations (DEs) since the 1980s. The interesting features in this method are possibility to find-out singularities, irregular structure and transient phenomena exhibited by the analysed equations such as by Heydari et al. [21], Wang and Fan [19], Balaji [22], Rehman and Khan [23], Hosseininia [24], Pirmohabbati et al. [25], Hosseininia [26], Heydari [27] and Kumar et al. [28].
Among the several wavelet families most simple are the Haar wavelets and it has been successfully applied to several linear and nonlinear problems of physical science and other research areas such as fractional order stationary neutron transport equation, neutron point kinetics equation, fractional order nonlinear oscillatory van der pol system and fractional bagley torvik equation by Ray and Patra [29,30,31,32], a comparative study on haar wavelet and hybrid functions, nonlinear integral and integro −differential equation of first and higher order and parabolic differential equatons by Aziz et al. [33,34,35], burgers equation by Jiwari [36], fractional integral equations by Lepik [37], Poisson and biharmonic equations by Shi and Cao [38], delamination detection in composite beams by Hein and Feklistova [39], fractional order integral equations by Gao and Liao [40], lumped and distributed parameters systems by Chen and Hsiao [41], FDEs by Chen et al. [42], free vibration analysis by Xie et al. [43], fractional nonlinear differential equations by Saeed and Rehman [44], magnetohydrodynamic flow equations by Celik and Brahin [45,46], fishers equations by Hariharan et al. [47], FPDEs by Wang et al. [48], nonlinear oscillators equations Kaur et al. [49], poisson and biharmonic equations by Shi et al. [50] and free vibration analysis of functionally graded cylindrical shells by Jin et al. [51].
It is compulsory to note that the fractional chemical kinetics and condensations of CO2 and PGE problems is the first one to be solved by the Haar wavelet and generalization of Adams– Bashforth−Moulton method by us. It is also noted that there are no similar works with these methods for fractional chemical kinetics and condensations of CO2 and PGE problems available in any present published literature. It is well known by the several published research papers that the Caputo and Riemann-Liouville is most popular definition of fractional calculus.
The complete work is systematized in the following sections: Overview of basic FC are provided in section 2. Fractional Model of both Chemical Kinetics and CO2 absorbed into PHE problems are provided in section 3. In section 4, a haar wavelet and Adam Bashforth's-Moulton methods are discussed and presented for both chemistry problems. The proposed methods for solutions of both chemistry problem are provided in section 5. Numerical result and discussions are provided in section 6. Conclusion and future scope are given in sections 7.
There are numerous definition of derivative and integration are available in literature [52,53,54,55,56,57,58,59,60,61].
Definition 1. The (left sided) Riemann−Liouville fractional integral of order α>0 of a function Θ(t)∈Cα,α≥−1 is defined as,
IαtΘ(t)=1Γ(α)t∫0(t−ξ)α−1Θ(ξ)dξ,α>0,t>0; | (2.1) |
where Γ(.) is well known Gamma function.
Definition 2. The next two equations define Riemann – Liouville and Caputo fractional derivatives of order a, respectively,
RLDαtΘ(t)=dmdtm(Im−αtΘ(t))={dmΘ(t)dtm,α=m∈N,1Γ(m−α)dmdtmt∫0Θ(ξ)(t−ξ)α−m+1dξ,0≤m−1<α<m,
and,
CDαtΘ(t)=Im−αt(dmdtmΘ(t))={dmΘ(t)dtm,α=m∈N,1Γ(m−α)t∫0Θm(ξ)(t−ξ)α−m+1dξ,0≤m−1<α<m,
where t>0 and m is an integer. Two basic properties for m−1<α≤m and Θ∈L1[a,b] are given as
{(CDαtIαΘ)(t)=Θ(t),(IαCDαtΘ)(t)=Θ(t)−∑m−1k=0Θk(0+)(t−a)kk!. | (2.2) |
Let D,E and H are different location of a model of chemical process then the reactions are presented as
D⟶E, | (3.1) |
E+H⟶D+H, | (3.2) |
E+E⟶H, | (3.3) |
The concetrations of all three spaces of D,E and H are denoted by Θ1,Θ2 and Θ3 respectively. Let r1,r2 and r3 denotes the reaction rate of Eqs (3.1), (3.2) and (3.3) respectively. We consider an integer order model of chemical kinetics problem as [1,2,6,7,8]
{dΘ1(t)dt=−r1Θ1(t)+r2Θ2(t)Θ3(t),dΘ2(t)dt=r1Θ1(t)−r2Θ2(t)Θ3(t)−r3Θ22(t),dΘ3(t)dt=r3Θ22(t), | (3.4) |
with the initial conditions, Θ1(0)=1, Θ2(0)=0, Θ3(0)=0. The main target of this section is converted above inter order CK problem into fractional order CK problem. The fractional model of CK problem is presented as
{CDαtΘ1(t)=−r1Θ1(t)+r2Θ2(t)Θ3(t),0<α≤1,CDβtΘ2(t)=r1Θ1(t)−r2Θ2(t)Θ3(t)−r3Θ22(t),0<β≤1,CDγtΘ3(t)=r3Θ22(t),0<γ≤1, | (3.5) |
with the initial conditions, Θ1(0)=1, Θ2(0)=0, Θ3(0)=0 where, Dαt=dαdtα,Dβt=dβdtβ,Dγt=dγdtγ are fractional derivative with 0<α,β,γ≤1. If r1=1,r2=0, and r3=1 then
{CDαtΘ1(t)=−Θ1(t),0<α≤1,CDβtΘ2(t)=Θ1(t)−Θ22(t),0<β≤1,CDγtΘ3(t)=Θ22(t),0<γ≤1, | (3.6) |
with the initial conditions, Θ1(0)=1, Θ2(0)=0, Θ3(0)=0. The above system is representing a nonlnear reaction which was taken from litrature [2,7,8,62].
The CO2 causes in ocean acidification because it dissolves in water to form carbonic acid [63].The mathematical formulation of the concentration of CO2 and PGE is shown in Muthukaruppan et al. [64]. Now, the two nonlinear reactions equations in normalized form is presented as
{d2Υ1dt2=α1Υ1Υ21+β1Υ1+β2Υ2,d2Υ2dt2=α2Υ1Υ21+β1Υ1+β2Υ2, | (3.7) |
with boundary conditions, Υ1(0)=0, Υ1(1)=1m, Υ′2(0)=1m, Υ2(1)=1m. The whole chemistry of the above problem is given in several litratures [1,3,4]. The fractional model of the condensation of CO2 and PGE in operator form is given as,
{CDαtΥ1(t)=α1Υ1Υ21+β1Υ1+β2Υ2,1<α≤2,CDβtΥ2(t)=α2Υ1Υ21+β1Υ1+β2Υ2,1<β≤2, | (3.8) |
with the same boundary conditions Υ1(0)=0, Υ1(1)=1m and Υ′2(0)=1m,Υ2(1)=1m; where m≥3 and fractional operator is taken in Caputo sence.
The Haar functions have been discovered by Alfred Haar in 1910 and Haar wavelets are the simplest wavelet among all wavelet. The Haar sequence was also introduced by itself Alfred Haar in 1909 which is recognised as wavelet basis. The Haar wavelets are the mathematical operations which are known as Haar transform. These wavelets are build up by piecewise constant function on the real line. We used Haar wavelet operational matrix method because of its flexibility, simplicity and require very less effort of computation. Usually Haar wavelet is defined for [0, 1) but in general case we extend it up to certain interval. Haar functions are very useful in many applications as image coding, extraction of edge, binary logic design etc [20,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51]. The Haar scaling function is defined as
ϕ(x)={10≤x<1,0otherwise. | (4.1) |
The Haar wavelet mother function is defined as
ψ(x)={10≤x<12,−112≤x<1,0otherwise. | (4.2) |
The orthogonal set of Haar wavelet functions for t∈[0,1] are defined as
hi(t)=1√m{2j/2,k−12j≤t<k−0.52j,−2j/2,k−0.52j≤t<k2j,0,otherwise, | (4.3) |
where i=0,1,2,...,m−1, m=2r+1 and r is positive integer known as resolution of Harr wavelet. Also j and k represent integer decomposition of i=2j+k−1.
Any function Θ(t)∈L2([0,1)) can be expanded in terms of Haar wavelet by
Θ(t)=∞∑i=0cihi(t);Whereci=1∫0Θ(t)hi(t)dt. | (4.4) |
If we approximated as piecewise constant during each interval, Eq. 4.4 will terminated at finite terms as [65]:
Θ(t)≈m−1∑i=0cihi(t)=CTmHm(t), | (4.5) |
where Cm=[c0,c1,c2,...,cm−1]T and Hm(t)=[h0(t),h1(t),h2(t),...,hm−1(t)]T,
Using collocation points tl=(l−0.5)m, where l=0,1,...,m−1, we obtained the discrete form as
H=[h0(t0)h0(t1)⋯h0(tm−1)h1(t0)h1(t1)⋯h1(tm−1)⋮⋮⋱⋮hm−1(t0)hm−1(t1)⋯hm−1(tm−1)]. | (4.6) |
The HWOM of fractional order integration without using block pulse functions we integrate Hm(t) using Reimann-Liouville integration operator [41,66]. Then the HWOM of fractional order integration Qα is given by
QαHm(t)=IαHm(t)=[Iαh0(t),Iαh1(t),Iαh2(t),...,Iαhm−1(t)]T =[Qh0(t),Qh1(t),Qh2(t),...,Qhm−1(t)]T, | (4.7) |
where
Qh0(t)=1√mtαΓ(1+α),0≤t≤1,
Qhi(t)=1√m{0,0≤t<k−12j,2j/2ζ1(t)k−12j≤t<k−0.52j,2j/2ζ2(t)k−0.52j≤t<k2j,2j/2ζ3(t)k2j≤t<1,
where
ζ1(t)=1Γ(1+α)(t−k−12j)α,
ζ2(t)=1Γ(1+α)(t−k−12j)α−2Γ(1+α)(t−k−0.52j)α,
ζ3(t)=1Γ(1+α)(t−k−12j)α−2Γ(1+α)(t−k−0.52j)α+1Γ(1+α)(t−k2j)α.
If we take, α=1/2,m=8, then we have the operational matrix as given below:
Q1/2H8 = [0.09970.17270.22300.26390.29920.33080.35960.38630.09970.17270.22300.26390.0997−0.0147−0.0864−0.14150.24430.0333−0.1154−0.0666−0.0343−0.0223−0.0193−0.013200000.14100.24430.033−0.11540.1995−0.0534−0.0455−0.0188−0.0111−0.0075−0.0055−0.0043000.1995−0.0534−0.0455−0.0188−0.0111−0.007500000.1995−0.0534−0.0455−0.01880000000.1995−0.0534].
The above matrix is the operational matrix of Haar wavelets.
In this section we discuss about Predictor-Corrector scheme (PECE), which is the genralization of (ABM) mehod [67,68]. We obtain the numerical solution of nonlinear FDES as
DαΘ(t)=f(t,Θ(t)),0<t≤T,Θ(k)(0)=Θ(k)0, | (4.8) |
where derivative in Caputo's sense. which is equivalent to the Volterra integral equation
Θ(t)=α−1∑k=0Θk0tkk!+1Γ(α)∫t0(t−τ)α−1f(t,Θ(τ))dτ. | (4.9) |
Assume h=T/N, tn=nh, n=0,1,2,...,N ∈Z+ then the discrete form for the above equation will be
Θh(tn+1)=α−1∑k=0Θ(k)0tkn+1k!+hαΓ(α+2)f(tn+1,Θph(tn+1))+hαΓ(α+2)n∑j=0aj,n+1f(th,Θh(tj)), | (4.10) |
aj,n+1={nα+1−(n−α)(n+1)α,ifj=0,(n−j+2)α+1+(n−j)α+1−2(n−j+1)α+1,if0≤j≤n,1,ifj=1, | (4.11) |
Θph(tn+1)=α−1∑k=0Θ(k)0tkn+1k!+1Γ(α)n∑j=0bj,n+1f(tj,Θh(tj)), | (4.12) |
bj,n+1=hαα((n+1−j)α−(n−j)α). | (4.13) |
The corrector values for chemistry problem is
Θ1(n+1)=Θ1(0)+hαΓ(α+2)(−r1Θp1(n+1)+r2Θp2(n+1)Θp3(n+1))+hαΓ(α+2)n∑j=0αj,n+1(−r1Θ1(j)+r2Θ2(j)Θ3(j)),Θ2(n+1)=Θ2(0)+hβΓ(β+2)(r1Θp1(n+1)−r2Θp2(n+1)Θp3(n+1)−r3Θp2(n+1)2)+hβΓ(β+2)n∑j=0βj,n+1(r1Θ1(j)−r2Θ2(j)Θ3(j)−r3Θ22(j)),Θ3(n+1)=Θ3(0)+hγΓ(γ+2)(r3Θp2(n+1)2)+hγΓ(γ+2)n∑j=0γj,n+1r3Θ22(j). |
The corresponding predictor values are,
Θp1(n+1)=Θ1(0)+1Γ(α)n∑j=0Bj,n+1(−r1Θ1(j)+r2Θ2(j)Θ3(j)),Θp2(n+1)=Θ2(0)+1Γ(β)n∑j=0Cj,n+1(r1Θ1(j)−r2Θ2(j)Θ3(j)−r3Θ22(j)),Θp3(n+1)=Θ3(0)+1Γ(γ)n∑j=0Dj,n+1(r3Θ22(j)). |
From Eqs (4.12) and (4.14) we can calculate {αj,n+1}, {βj,n+1}, {γj,n+1}, and Bj,n+1, Cj,n+1, Dj,n+1.
Example: 1 We assume a fractional model of chemical kinetics problem is given as
{CDαtΘ1=−r1Θ1+c2Θ2Θ3,0<α≤1,CDβtΘ2=r1Θ1−r2Θ2Θ3−r3Θ22,0<β≤1,CDγtΘ3=r3Θ22,0<γ≤1, | (5.1) |
with the initial conditions, Θ1(0)=1, Θ2(0)=0 Θ3(0)=0, where r1, r2 and r3 are reaction rates. Let us assume higher derivatives in the terms of haar wavelet series.
{CDαtΘ1=CTHm(t),CDβtΘ2=GTHm(t),CDγtΘ3=KTHm(t), | (5.2) |
where C=[c0,c1,c2,...,cm−1]T, G=[g0,g1,g2,...,gm−1]T and K=[k0,k1,k2,...,km−1]T are unknown vectors. Applying Riemann-Liouville fractional integral in Eq. (5.2) and using initial conditions, we obtained
{Θ1=CTQαHm(t)+1,Θ2=GTQβHm(t),Θ3=KTQγHm(t). | (5.3) |
Now substituting the values of Θ1, Θ2 and Θ3 into the Eq. (5.1), we obtained.
{CTHm(t)=−r1(CTQαHm(t)+1)+r2(GTQβ)(KTQγ),GTHm(t)=r1(CTQαHm(t)+1)−r2(GTQβ)(KTQγ)−r3(KTQγ)(KTQγ),KTHm(t)=r3(GTQβ)2. | (5.4) |
Let r1=0.1, r2=0.02 and r3=0.009 as given in Aminikhah [69]. Now disperse the Eq. (5.4) at the collocation points tl=(l−0.5)m, where l=1,2,3,...,m. We obtained 3m nonlinear algebraic equations which can be solved by Newton iteration method, after solving we obtained the coefficients ci, gi and ki. Substitute these coefficients into the Eq. (5.3) we get desired solutions Θ1, Θ2 and Θ3.
Example 2: Consider the system of condensations of CO2 and PGE problem of arbitrary order.
{Υα1(t)=α1Υ1(t)Υ2(t)−Υα1(t)(β1Υ1(t)+β2Υ2(t)),1<α≤2,Υβ1(t)=α2Υ1(x)Υ2(x)−Υβ2(t)(β1Υ1(t)+β2Υ2(t)),1<β≤2, | (5.5) |
with boundary conditions Υ1(0)=0, Υ1(1)=1m, and Υ′2(0)=1m, Υ2(1)=1m where Υα1(t)= CDαtΥ1(t). Here for simplicity we have taken m=3 and we will take the value of α1=1, α2=2, β1=1 and β2=3 as given in Duan et al. [70], AL-jawary ad Radhi [71]. Further, we assume the higher derivative in terms of Haar wavelet series.
{Υα1(t)=CTHm(t),Υβ1(t)=KTHm(t), | (5.6) |
applying Riemann-Liouville integral operator on the above equation and using boundary conditions, we obtained
Υ1(t)−Υ′1(0)t=CTQαHm(t), | (5.7) |
substituting t=1 into Eq. (5.7) we obtained
Υ1(1)−Υ′1(0)=CTQαHm(1)
Υ′1(0)=13−CTQαHm(1), | (5.8) |
and
Υ2(0)=−KTQβHm(1). | (5.9) |
Therefore,
Υ1(t)=(13−CTQαHm(1))t+CTQαHm(t), | (5.10) |
similarly
Υ2(t)=t3−KTQβHm(1)+KTQβHm(t). | (5.11) |
Substituting the values of Υ1, Υ2 into the Eq. (5.5) and using Eq. (5.6) we obtained
CTHm(t)=(t3−CTQαHm(1)t+CTQαHm(t))(t3−KTQβHm(1)+KTQβHm(t))−CTHm(t)((t3−CTQαHm(1)t+CTQαHm(t))+3(t3−KTQβHm(1)+KTQβHm(t))).} | (5.12) |
KTHm(t)=2(t3−CTQαHm(1)t+CTQαHm(t))(t3−KTQβHm(1)+KTQβHm(t))−KTHm(t)((t3−CTQαHm(1)t+CTQαHm(t))+3(t3−KTQβHm(1)+KTQβHm(t))).} | (5.13) |
Now disperse the Eqs (5.12) and (5.13) at the collocation points tl=(l−0.5)m, where l=0,1,...,m−1. We obtained a system of nonlinear algebraic equations which can be easily solved by Newton-Iteration method using mathematical softwares, after solving we obtained the unknowns coefficients ci and ki. Substituting these coefficients into the Eqs (5.10) and (5.11) we get desired solutions Υ1 and Υ2.
All numerical simulation and graphical results of both examples are depicted through the Figures 1–14 where Figures 1–6 and Figures 7–14 are depicted for examples 1 and 2 respectively. We have depicted a comparison between numerical obtained solutions using by Haar wavelet and Adam's-Bashforth-Moulton predictor-corrector schemes through the Figures 1–3 and these figures are depicted for the values of m=64. It is clear from all figures that both obtained solutions by HWM and ABM are identical. The obtained solutions Θ1, Θ2 and Θ3 are plotted through the Figures 3–6 where the nature of solution Θ1 is of decreasing nature while other solutions Θ2 and Θ3 is of increasing nature. We plotted the resolutions Figures 7–14 for better understanding the nature of obtained solution of example 2. We plotted resolutions figures due to non-availability of its exact solution.
t | Θ1(HWM) | Θ1(ABM) | Θ2(HWM) | Θ2(ABM) |
0.1 | 0.9901 | 0.9893 | 0.0100 | 0.0107 |
0.2 | 0.9802 | 0.9794 | 0.0198 | 0.0206 |
0.3 | 0.9704 | 0.9697 | 0.0296 | 0.0303 |
0.4 | 0.9608 | 0.9600 | 0.0393 | 0.0400 |
0.5 | 0.9512 | 0.9505 | 0.0489 | 0.0495 |
0.6 | 0.9418 | 0.9410 | 0.0585 | 0.0590 |
0.7 | 0.9324 | 0.9317 | 0.0680 | 0.0683 |
0.8 | 0.9231 | 0.9224 | 0.0775 | 0.0776 |
0.9 | 0.9139 | 0.9132 | 0.0869 | 0.0868 |
1.0 | 0.9048 | 0.9041 | 0.0963 | 0.0962 |
t | Θ3(HWM) | Θ3(ABM) |
0.1 | 3.0×10−8 | 4.0×10−8 |
0.2 | 2.4×10−7 | 2.7×10−7 |
0.3 | 7.9×10−7 | 8.6×10−7 |
0.4 | 1.8×10−6 | 1.9×10−6 |
0.5 | 3.6×10−6 | 3.8×10−6 |
0.6 | 6.2×10−6 | 6.4×10−6 |
0.7 | 9.8×10−6 | 1.0×10−5 |
0.8 | 1.5×10−5 | 1.5×10−5 |
0.9 | 2.0×10−5 | 2.0×10−5 |
1.0 | 2.8×10−5 | 2.8×10−5 |
In this work, Haar wavelet operational matrix and Adam Bashforth's Moulton scheme are proposed to solve fractional chemical kinetics and another problem that relates the condensations of carbon dioxide CO2 numerically. A comparative study between fractional chemical kinetics and another problem that relates the conden sations of carbon dioxide CO2 has been done for m=64 in this work. Our tabulated and graphical results indicate that the solution will ameliorate if we will take more collocation points, i.e greater values of m. The essential advantage of HWM is that it converts problems into the system of linear or nonlinear algebraic equations so that the computation is facile and computer-oriented. Furthermore, wavelet method is much easier than other numerical methods for system of FDEs. Again, we have solved the chemistry problems at different resolutions, which produced the same results at each resolution. The precision of the solution will ameliorate if we increase the resolution. This new comparative study between the Haar wavelet operational matrix and Adam Bashforth's Moulton scheme for fractional chemical kinetics and another problem that relates the condensations of carbon dioxide CO2 indicates that both approaches can be applied successfully to the chemistry problems of chemistry science.
The first author Dr. Sunil Kumar would like to acknowledge the financial support received from the National Board for Higher Mathematics, Department of Atomic Energy, Government of India (Approval No. 2/48(20)/2016/ NBHM(R.P.)/R and D II/1014). The authors are also grateful to the editor and anonymous reviewers for their constructive comments and valuable suggestions to improve the quality of article.
The authors declare no conflict of interest in this manuscript.
[1] |
H. Chen, X. Feng, Z. Zhang, A recovery-based linear C0 finite element method for a fourth-order singularly perturbed Monge-Ampère equation, Adv. Comput. Math., 47 (2021), 21. https://doi.org/10.1007/s10444-021-09847-w doi: 10.1007/s10444-021-09847-w
![]() |
[2] |
Z. Cen, L. Liu, A. Xu, A second-order adaptive grid method for a nonlinear singularly perturbed problem with an integral boundary condition, J. Comput. Appl. Math., 385 (2021), 113205. https://doi.org/10.1016/j.cam.2020.113205 doi: 10.1016/j.cam.2020.113205
![]() |
[3] |
H. G. Debela, G. F. Duressa, Accelerated exponentially fitted operator method for singularly perturbed problems with integral boundary condition, Int. J. Differ. Equat., 2020 (2020), 9268181. https://doi.org/10.1155/2020/9268181 doi: 10.1155/2020/9268181
![]() |
[4] |
M. Denche, A. Kourta, Boundary value problem for second-order differential operators with integral conditions, Appl. Anal., 84 (2005), 1247–1266. https://doi.org/10.1080/00036810500287255 doi: 10.1080/00036810500287255
![]() |
[5] |
M. E. Durmaz, I. Amirali, G. M. Amiraliyev, An efficient numerical method for a singularly perturbed fredholm integro-differential equation with integral boundary condition, J. Appl. Math. Comput., 69 (2022), 505–528. https://doi.org/10.1007/s12190-022-01757-4 doi: 10.1007/s12190-022-01757-4
![]() |
[6] | P. A. Farrell, A. F. Hegarty, J. J. H. Miller, E. O'Riordan, G. I. Shishkin, Robust computational techniques for boundary layers, Boca Raton: CRC Press, 2000. |
[7] |
F. T. Fen, I. Y. Karaca, Positive solutions of nth-order impulsive differential equations with integral boundary conditions, An. Sti. U. Ovi. Co. Mat., 24 (2016), 243–261. https://doi.org/10.1515/auom-2016-0014 doi: 10.1515/auom-2016-0014
![]() |
[8] |
Y. Fu, H. Yao, The solution of nonlinear fourth-order differential equation with integral boundary conditions, J. Funct. Space., 2014 (2014), 890695. https://doi.org/10.1155/2014/890695 doi: 10.1155/2014/890695
![]() |
[9] |
N. Geetha, A. Tamilselvan, Parameter uniform numerical method for fourth order singularly perturbed turning point problems exhibiting boundary layers, Ain Shams Eng. J., 9 (2018), 845–853. https://doi.org/10.1016/j.asej.2016.04.018 doi: 10.1016/j.asej.2016.04.018
![]() |
[10] |
H. Guo, Z. Zhang, Q. Zou, A c0 linear finite element method for biharmonic problems, J. Sci. Comput., 74 (2018), 1397–1422. https://doi.org/10.1007/s10915-017-0501-0 doi: 10.1007/s10915-017-0501-0
![]() |
[11] |
G. J. Jayalakshmi, A. Tamilselvan, Comparative study on difference schemes for singularly perturbed boundary turning point problems with robin boundary conditions, J. Appl. Math. Comput., 62 (2020), 341–360. https://doi.org/10.1007/s12190-019-01287-6 doi: 10.1007/s12190-019-01287-6
![]() |
[12] |
T. Jankowski, Positive solutions for fourth-order differential equations with deviating arguments and integral boundary conditions, Nonlinear Anal., 73 (2010), 1289–1299. https://doi.org/10.1016/j.na.2010.04.055 doi: 10.1016/j.na.2010.04.055
![]() |
[13] |
Y. Jia, M. Xu, Y. Lin, D. Jiang, An efficient technique based on least-squares method for fractional integro-differential equations, Alex. Eng. J., 64 (2023), 97–105. https://doi.org/10.1016/j.aej.2022.08.033 doi: 10.1016/j.aej.2022.08.033
![]() |
[14] |
M. Kudu, A parameter uniform difference scheme for the parameterized singularly perturbed problem with integral boundary condition, Adv. Differ. Equ., 2018 (2018), 170. https://doi.org/10.1186/s13662-018-1620-0 doi: 10.1186/s13662-018-1620-0
![]() |
[15] |
M. Kudu, I. Amirali, G. M. Amiraliyev, A second order accurate method for a parameterized singularly perturbed problem with integral boundary condition, J. Comput. Appl. Math., 404 (2022), 113894. https://doi.org/10.1016/j.cam.2021.113894 doi: 10.1016/j.cam.2021.113894
![]() |
[16] |
Y. Li, H. Zhang, Positive solutions for a nonlinear higher order differential system with coupled integral boundary conditions, J. Appl. Math., 2014 (2014), 901094. https://doi.org/10.1155/2014/901094 doi: 10.1155/2014/901094
![]() |
[17] |
R. K. Lodhi, H. K. Mishra, Solution of a class of fourth order singular singularly perturbed boundary value problems by quintic b-spline method, J. Nigerian Math. Soc., 35 (2016), 257–265. https://doi.org/10.1016/j.jnnms.2016.03.002 doi: 10.1016/j.jnnms.2016.03.002
![]() |
[18] |
R. Mahendran, V. Subburayan, Fitted finite difference method for third order singularly perturbed delay differential equations of convection diffusion type, Int. J. Comput. Method., 16 (2019), 1840007. https://doi.org/10.1142/S0219876218400078 doi: 10.1142/S0219876218400078
![]() |
[19] | J. J. H. Miller, E. O'Riordan, G. I. Shishkin, Fitted numerical methods for singular perturbation problems, Singapore: World scientific, 2012. |
[20] |
H. K. Mishra, Fourth order singularly perturbed boundary value problems via initial value techniques, Appl. Math. Sci., 8 (2014), 619–632. https://doi.org/10.12988/AMS.2014.312676 doi: 10.12988/AMS.2014.312676
![]() |
[21] |
V. Raja, A. Tamilselvan, Fitted finite difference method for third order singularly perturbed convection diffusion equations with integral boundary condition, Arab J. Math. Sci., 25 (2019), 231–242. https://doi.org/10.1016/j.ajmsc.2018.10.002 doi: 10.1016/j.ajmsc.2018.10.002
![]() |
[22] |
V. Raja, A. Tamilselvan, Numerical method for a system of singularly perturbed convection diffusion equations with integral boundary conditions, Commun. Korean Math. Soc., 34 (2019), 1015–1027. https://doi.org/10.4134/CKMS.c180224 doi: 10.4134/CKMS.c180224
![]() |
[23] | J. Stalin Christy Roja, A. Tamilselvan, Numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type, J. Math. Model., 4 (2016), 79–102. |
[24] |
E. Sekar, Second order singularly perturbed delay differential equations with non-local boundary condition, J. Comput. Appl. Math., 417 (2023), 114498. https://doi.org/10.1016/j.cam.2022.114498 doi: 10.1016/j.cam.2022.114498
![]() |
[25] |
E. Sekar, U. Bundit, Numerical scheme for singularly perturbed mixed delay differential equation on shishkin type meshes, Fractal Fract., 7 (2023), 43. https://doi.org/10.3390/fractalfract7010043 doi: 10.3390/fractalfract7010043
![]() |
[26] |
E. Sekar, A. Tamilselvan, Finite difference scheme for third order singularly perturbed delay differential equation of convection diffusion type with integral boundary condition, J. Appl. Math. Comput., 61 (2019), 73–86. https://doi.org/10.1007/s12190-019-01239-0 doi: 10.1007/s12190-019-01239-0
![]() |
[27] |
V. Shanthi, N. Ramanujam, A numerical method for boundary value problems for singularly perturbed fourth-order ordinary differential equations, Appl. Math. Comput., 129 (2002), 269–294. https://doi.org/10.1016/S0096-3003(01)00040-6 doi: 10.1016/S0096-3003(01)00040-6
![]() |
[28] |
Q. Wang, Y. Guo, Y. Ji, Positive solutions for fourth-order nonlinear differential equation with integral boundary conditions, Discrete Dyn. Nat. Soc., 2013 (2013), 292–297. https://doi.org/10.1155/2013/684962 doi: 10.1155/2013/684962
![]() |
[29] |
M. Xu, R. Lin, Q. Zou, A c0 linear finite element method for a second-order elliptic equation in non-divergence form with cordes coefficients, Numer. Meth. Part. Differ. Equ., 39 (2023), 2244–2269. https://doi.org/10.1002/num.22965 doi: 10.1002/num.22965
![]() |
[30] |
M. Xu, C. Shi, A hessian recovery-based finite difference method for biharmonic problems, Appl. Math. Lett., 137 (2023), 108503. https://doi.org/10.1016/j.aml.2022.108503 doi: 10.1016/j.aml.2022.108503
![]() |
[31] |
M. Xu, E. Tohidi, J. Niu, Y. Fang, A new reproducing kernel-based collocation method with optimal convergence rate for some classes of bvps, Appl. Math. Comput., 432 (2022), 127343. https://doi.org/10.1016/j.amc.2022.127343 doi: 10.1016/j.amc.2022.127343
![]() |
[32] |
Y. Zhu, H. Pang, The shooting method and positive solutions of fourth-order impulsive differential equations with multi-strip integral boundary conditions, Adv. Differ. Equ., 2018 (2018), 5. https://doi.org/10.1186/s13662-017-1453-2 doi: 10.1186/s13662-017-1453-2
![]() |
1. | Ann Al Sawoor, Stability analysis of fractional-order linear neutral delay differential–algebraic system described by the Caputo–Fabrizio derivative, 2020, 2020, 1687-1847, 10.1186/s13662-020-02980-8 | |
2. | Amr MS Mahdy, Yasser Abd Elaziz Amer, Mohamed S Mohamed, Eslam Sobhy, General fractional financial models of awareness with Caputo–Fabrizio derivative, 2020, 12, 1687-8140, 168781402097552, 10.1177/1687814020975525 | |
3. | Gustavo Asumu Mboro Nchama, Angela Leon Mecias, Mariano Rodriguez Ricard, Perona-Malik Model with Diffusion Coefficient Depending on Fractional Gradient via Caputo-Fabrizio Derivative, 2020, 2020, 1085-3375, 1, 10.1155/2020/7624829 | |
4. | Juan Bory‐Reyes, Marco Antonio Pérez‐de la Rosa, Local fractional Moisil–Teodorescu operator in quaternionic setting involving Cantor‐type coordinate systems, 2021, 44, 0170-4214, 605, 10.1002/mma.6767 | |
5. | Xiuying Li, Boying Wu, Reproducing kernel functions-based meshless method for variable order fractional advection-diffusion-reaction equations, 2020, 59, 11100168, 3181, 10.1016/j.aej.2020.07.034 | |
6. | Zeyad Al-Zhour, Fundamental fractional exponential matrix: New computational formulae and electrical applications, 2021, 129, 14348411, 153557, 10.1016/j.aeue.2020.153557 | |
7. | Hu Ge-JiLe, Saima Rashid, Muhammad Aslam Noor, Arshiya Suhail, Yu-Ming Chu, Some unified bounds for exponentially tgs-convex functions governed by conformable fractional operators, 2020, 5, 2473-6988, 6108, 10.3934/math.2020392 | |
8. | Doddabhadrappla Gowda Prakasha, Naveen Sanju Malagi, Pundikala Veeresha, Ballajja Chandrappa Prasannakumara, An efficient computational technique for time‐fractional Kaup‐Kupershmidt equation , 2021, 37, 0749-159X, 1299, 10.1002/num.22580 | |
9. | Naveed Ahmad Khan, Muhammad Sulaiman, Carlos Andrés Tavera Romero, Fawaz Khaled Alarfaj, Theoretical Analysis on Absorption of Carbon Dioxide (CO2) into Solutions of Phenyl Glycidyl Ether (PGE) Using Nonlinear Autoregressive Exogenous Neural Networks, 2021, 26, 1420-3049, 6041, 10.3390/molecules26196041 | |
10. | Devendra Kumar, Hunney Nama, Dumitru Baleanu, Numerical and computational analysis of fractional order mathematical models for chemical kinetics and carbon dioxide absorbed into phenyl glycidyl ether, 2023, 53, 22113797, 107003, 10.1016/j.rinp.2023.107003 | |
11. | Khushbu Agrawal, Sunil Kumar, Ali Akgül, An algorithm for numerical study of fractional atmospheric model using Bernoulli polynomials, 2024, 70, 1598-5865, 3101, 10.1007/s12190-024-02084-6 | |
12. | Khushbu Agrawal, Sunil Kumar, Badr Saad T. Alkahtani, Sara S. Alzaid, The use of Hermite wavelet collocation method for fractional cancer dynamical system, 2024, 32, 2769-0911, 10.1080/27690911.2024.2352745 | |
13. | G. Manohara, S. Kumbinarasaiah, Fibonacci wavelets operational matrix approach for solving chemistry problems, 2023, 9, 2731-6734, 393, 10.1007/s43994-023-00046-5 | |
14. | Khushbu Agrawal, Sunil Kumar, Badr S.T. Alkahtani, Sara S. Alzaid, A numerical study on fractional order financial system with chaotic and Lyapunov stability analysis, 2024, 60, 22113797, 107685, 10.1016/j.rinp.2024.107685 |
t | Θ1(HWM) | Θ1(ABM) | Θ2(HWM) | Θ2(ABM) |
0.1 | 0.9901 | 0.9893 | 0.0100 | 0.0107 |
0.2 | 0.9802 | 0.9794 | 0.0198 | 0.0206 |
0.3 | 0.9704 | 0.9697 | 0.0296 | 0.0303 |
0.4 | 0.9608 | 0.9600 | 0.0393 | 0.0400 |
0.5 | 0.9512 | 0.9505 | 0.0489 | 0.0495 |
0.6 | 0.9418 | 0.9410 | 0.0585 | 0.0590 |
0.7 | 0.9324 | 0.9317 | 0.0680 | 0.0683 |
0.8 | 0.9231 | 0.9224 | 0.0775 | 0.0776 |
0.9 | 0.9139 | 0.9132 | 0.0869 | 0.0868 |
1.0 | 0.9048 | 0.9041 | 0.0963 | 0.0962 |
t | Θ3(HWM) | Θ3(ABM) |
0.1 | 3.0×10−8 | 4.0×10−8 |
0.2 | 2.4×10−7 | 2.7×10−7 |
0.3 | 7.9×10−7 | 8.6×10−7 |
0.4 | 1.8×10−6 | 1.9×10−6 |
0.5 | 3.6×10−6 | 3.8×10−6 |
0.6 | 6.2×10−6 | 6.4×10−6 |
0.7 | 9.8×10−6 | 1.0×10−5 |
0.8 | 1.5×10−5 | 1.5×10−5 |
0.9 | 2.0×10−5 | 2.0×10−5 |
1.0 | 2.8×10−5 | 2.8×10−5 |
t | Θ1(HWM) | Θ1(ABM) | Θ2(HWM) | Θ2(ABM) |
0.1 | 0.9901 | 0.9893 | 0.0100 | 0.0107 |
0.2 | 0.9802 | 0.9794 | 0.0198 | 0.0206 |
0.3 | 0.9704 | 0.9697 | 0.0296 | 0.0303 |
0.4 | 0.9608 | 0.9600 | 0.0393 | 0.0400 |
0.5 | 0.9512 | 0.9505 | 0.0489 | 0.0495 |
0.6 | 0.9418 | 0.9410 | 0.0585 | 0.0590 |
0.7 | 0.9324 | 0.9317 | 0.0680 | 0.0683 |
0.8 | 0.9231 | 0.9224 | 0.0775 | 0.0776 |
0.9 | 0.9139 | 0.9132 | 0.0869 | 0.0868 |
1.0 | 0.9048 | 0.9041 | 0.0963 | 0.0962 |
t | Θ3(HWM) | Θ3(ABM) |
0.1 | 3.0×10−8 | 4.0×10−8 |
0.2 | 2.4×10−7 | 2.7×10−7 |
0.3 | 7.9×10−7 | 8.6×10−7 |
0.4 | 1.8×10−6 | 1.9×10−6 |
0.5 | 3.6×10−6 | 3.8×10−6 |
0.6 | 6.2×10−6 | 6.4×10−6 |
0.7 | 9.8×10−6 | 1.0×10−5 |
0.8 | 1.5×10−5 | 1.5×10−5 |
0.9 | 2.0×10−5 | 2.0×10−5 |
1.0 | 2.8×10−5 | 2.8×10−5 |