Communication Special Issues

Cold helium plasma as a modifier of free radical processes in the blood: in vitro study

  • We studied the influence of various exposures (1, 2, 3, 5, 10 and 15 min.) effects of cold helium plasma on the state of free radical processes in the biological fluid (human blood). It was found that the nature of the response of antioxidant systems of blood to cold plasma treatment is directly determined by exposure, as evidenced by the results of chemiluminescence assessment of the intensity of free radical processes and total antioxidant activity, as well as the dynamics of the concentration of malonic dialdehyde. At the same time, the short-term processing mode of biological fluid allows us to establish the presence of a two-phase influence of the factor on oxidative processes (“the phenomenon of the antioxidant window”).

    Citation: Andrew K. Martusevich, Alexander G. Galka, Konstantin A. Karuzin, Alexander N. Tuzhilkin, Svetlana L. Malinovskaya. Cold helium plasma as a modifier of free radical processes in the blood: in vitro study[J]. AIMS Biophysics, 2021, 8(1): 34-40. doi: 10.3934/biophy.2021002

    Related Papers:

    [1] Yinlin Ye, Hongtao Fan, Yajing Li, Ao Huang, Weiheng He . An artificial neural network approach for a class of time-fractional diffusion and diffusion-wave equations. Networks and Heterogeneous Media, 2023, 18(3): 1083-1104. doi: 10.3934/nhm.2023047
    [2] Kexin Li, Hu Chen, Shusen Xie . Error estimate of L1-ADI scheme for two-dimensional multi-term time fractional diffusion equation. Networks and Heterogeneous Media, 2023, 18(4): 1454-1470. doi: 10.3934/nhm.2023064
    [3] Narcisa Apreutesei, Vitaly Volpert . Reaction-diffusion waves with nonlinear boundary conditions. Networks and Heterogeneous Media, 2013, 8(1): 23-35. doi: 10.3934/nhm.2013.8.23
    [4] Yaxin Hou, Cao Wen, Yang Liu, Hong Li . A two-grid ADI finite element approximation for a nonlinear distributed-order fractional sub-diffusion equation. Networks and Heterogeneous Media, 2023, 18(2): 855-876. doi: 10.3934/nhm.2023037
    [5] Zhe Pu, Maohua Ran, Hong Luo . Effective difference methods for solving the variable coefficient fourth-order fractional sub-diffusion equations. Networks and Heterogeneous Media, 2023, 18(1): 291-309. doi: 10.3934/nhm.2023011
    [6] Feiyang Peng, Yanbin Tang . Inverse problem of determining diffusion matrix between different structures for time fractional diffusion equation. Networks and Heterogeneous Media, 2024, 19(1): 291-304. doi: 10.3934/nhm.2024013
    [7] Panpan Xu, Yongbin Ge, Lin Zhang . High-order finite difference approximation of the Keller-Segel model with additional self- and cross-diffusion terms and a logistic source. Networks and Heterogeneous Media, 2023, 18(4): 1471-1492. doi: 10.3934/nhm.2023065
    [8] Aslam Khan, Abdul Ghafoor, Emel Khan, Kamal Shah, Thabet Abdeljawad . Solving scalar reaction diffusion equations with cubic non-linearity having time-dependent coefficients by the wavelet method of lines. Networks and Heterogeneous Media, 2024, 19(2): 634-654. doi: 10.3934/nhm.2024028
    [9] Caihong Gu, Yanbin Tang . Global solution to the Cauchy problem of fractional drift diffusion system with power-law nonlinearity. Networks and Heterogeneous Media, 2023, 18(1): 109-139. doi: 10.3934/nhm.2023005
    [10] Verónica Anaya, Mostafa Bendahmane, David Mora, Ricardo Ruiz Baier . On a vorticity-based formulation for reaction-diffusion-Brinkman systems. Networks and Heterogeneous Media, 2018, 13(1): 69-94. doi: 10.3934/nhm.2018004
  • We studied the influence of various exposures (1, 2, 3, 5, 10 and 15 min.) effects of cold helium plasma on the state of free radical processes in the biological fluid (human blood). It was found that the nature of the response of antioxidant systems of blood to cold plasma treatment is directly determined by exposure, as evidenced by the results of chemiluminescence assessment of the intensity of free radical processes and total antioxidant activity, as well as the dynamics of the concentration of malonic dialdehyde. At the same time, the short-term processing mode of biological fluid allows us to establish the presence of a two-phase influence of the factor on oxidative processes (“the phenomenon of the antioxidant window”).


    Time-fractional partial differential equations (TFPDEs) have attracted considerable attention due to their ability to model memory and nonlocal properties. It has been established that TFPDEs are crucial mathematical and physical models for describing a wide range of anomalous phenomena and complex systems in the fields of natural science and engineering [1,2,3,4,5]. Successful applications of TFPDEs include signal processing [6], Powell-Eyring fluid [7], fluid mechanics [8,9], and robotics [10]. The multi-term TFPDEs have been found useful in describing many complex natural phenomena, such as magnetic resonance imaging (MRI) [11], viscoelastic mechanical models [12], and oxygen delivery through capillaries to a tissue (see [13] and the references therein). The multi-term time-fractional mixed sub-diffusion and diffusion wave equations belong to an important class of multi-term TFPDEs that have three major subclasses: the time-fractional sub-diffusion equations (TFSDEs), whose fractional order belongs to (0,1); the time fractional diffusion-wave equations (TFDWEs), whose fractional order belongs to (1,2); and TFDEs whose fractional order belongs to (0,2). Some researchers studied the analytical solutions of these equations, such as the authors of [14] who studied the solution of multi-term TFDWEs using the method of separation of variables. In [15], the authors used the Laguerre polynomials for obtaining the approximate solution of 1D TFDWE. The authors of [16] obtained the analytic solution of the multi-term space-time fractional advection-diffusion equation using a method based on spectral representation of fractional Laplacian operator. In most cases, the analytical solution of the these equations is provided in terms of special functions such as the multivariate Mittag-Leffler function and the Fox H-function [17,18], which are very complex and challenging to evaluate. Therefore, to solve these equations, numerical methods would be preferable, particularly in situations where analytical solutions are not present. Various numerical schemes have been developed to solve TFPDEs. Liu et al.[19] developed a first-order finite difference approach with stability constraints for both time and spatial derivatives. The multi-term time-fractional mixed sub-diffusion with a variable coefficient was introduced by Zhao [20] using the finite element method. The solution of the wave-diffusion equation based on the Legendre spectral method and finite difference method for discretization of space and time derivatives was established in [21]. Shen and Gu [22] developed two novel finite difference schemes for 2D time fractional mixed diffusion and wave-diffusion equations. Agrawal [23] utilized the Laplace transform and finite sine transform for obtaining the solution of the wave-diffusion equation. In [24], the authors coupled spatial extrapolation method with the Crank-Nicholson method to obtain the solution to the fractional diffusion problem. In [25], a fully discrete spectral method was employed for the solution of a novel multi-term time-fractional mixed diffusion and diffusion-wave equation. The mesh-less method based on radial basis function was employed by [26] for the approximate solution of multi-term time-fractional mixed sub-diffusion. In [27], the authors coupled the dual reciprocity method and the improved singular boundary method for multi-term fractional wave-diffusion equations. Ye et al. [28] solved the 2D and 3D multi-term time-space fractional diffusion equations via a series expansion. Li et al.[29] analyzed the temporal asymptotic behavior with well-posedness for the solution of the multi-term time-fractional diffusion equation. The authors of [30] employed the spectral method for the numerical solution of the two-dimensional multi-term time-fractional mixed sub-diffusion. Sun [31] examined the simultaneous inversion of the potential term and the fractional orders in a multi-term time-fractional diffusion equation. Spectral methods are among the most popular techniques for discretizing spatial variables in partial differential equations, which have the reputation of producing extremely accurate approximations of sufficiently smooth solutions [32,33]. Over time, spectral methods have been widely used in many fields, such as quantum mechanics [34], fluid dynamics [35], weather forecasting [36], and heat conduction [37], due to their high-order accuracy. These techniques haven't, however, been applied to many problems where the finite-difference and finite-element approaches are still frequently used because of some disadvantages. One disadvantage is that using spectral methods to discretize partial differential equations results in the solution of large systems of linear or nonlinear equations that require full matrices. On the other hand, finite-difference and finite-element approaches produce sparse matrices, which can be handled by suitable techniques to significantly reduce the computational complexity. Another drawback of spectral methods is that they face difficulties for problems defined on complex geometries. Although, some authors have attempted to use spectral methods in complex geometries[38,39]. In this work, our aim is to obtain the approximate solution of the multi-term time-fractional mixed sub-diffusion using the pseudospectral method in space coupled with the Laplace transform in time. The pseudospectral method is an easy method to implement and offers low computational cost with high accuracy. Also, it has been proven that for time discretization, the Laplace transform is suitable for solving various initial and boundary value problems and can be utilized as an alternative to the classical time stepping techniques [40]. In finite difference time stepping techniques, the accuracy can be obtained for extremely short time steps, which results in high computational cost and crucial stability constraints. A large number of valuable works are available in the literature, which couples Laplace transform in time with another method in space (see [41,42] and references therein). In this article, we consider a two-dimensional multi-term time-fractional mixed sub-diffusion and diffusion-wave equation of the following form [43]:

    pk=1ϱ1,kαku(ˉξ,t)tαk+ϱ2Dtu(ˉξ,t)+qr=1ϱ3,rβru(ˉξ,t)tβr+ϱ4u(ˉξ,t)=ϱ5Δu(ˉξ,t)+ϱ6γΔu(ˉξ,t)tγ+G(ˉξ,t),ˉξΓ,t[0,1] (1.1)

    with boundary conditions

    LBu(ˉξ,t)=h(ˉξ,t),ˉξΓ,t[0,1] (1.2)

    and initial condition

    u(ˉξ,0)=g1(ˉξ),ut(ˉξ,0)=g2(ˉξ),ˉξΓ (1.3)

    where G(ˉξ,t),g1(ˉξ),g2(ˉξ),h(ˉξ,t) are given continuous functions, ϱ1,k>0,ϱ3,r>0,ϱl>0, l=2,4,5,6, the coefficients are not simultaneously equal to zero, 1<α1<α2<...<αp<2,0<β1<β2<...<βq<1,0<γ<1, ˉξ=(ξ,ζ),Dt=t,Δ=2ξ2+2ζ2,LB is the boundary operator, ΓR2 is the domain, Γ is its boundary, and αu(ˉξ,t)tα, βu(ˉξ,t)tβ, and γu(ˉξ,t)tγ are fractional derivatives in Caputo's sense, defined as in [44].

    αu(ˉξ,t)tα=1Γ(nα)t0nu(ˉξ,ν)νn(tν)αn+1dν,n1<αn

    In literature, there are many physical processes whose characteristics cannot be described by utilizing the time fractional sub-diffusion equations or the time fractional wave diffusion equation independently. In order to avoid such limitations, the multi-term time-fractional mixed sub-diffusion can be utilized as an effective mathematical model [43]. These models are particularly beneficial at representing anomalous diffusion processes and capturing the properties of media, including power-law frequency dependence [30]. Additionally, they are proficient in the modeling of various types of viscoelastic damping, modeling the unsteady flow of a fractional Maxwell fluid [45], and describing the behavior of an Oldroyd-B fluid, which has been used to simulate the response of many dilute polymeric liquids [22].

    In this section, the Laplace transform is employed on the multi-term time fractional mixed sub-diffusion and diffusion-wave equation in order to reduce the problem into an equivalent time-independent problem. The Laplace transform of u(ˉξ,t) is denoted and defined as in [44]

    ˆu(ˉξ,s)=0estu(ˉξ,t)dt

    The Laplace transform of Caputo's derivative αu(ˉξ,t)tα of fractional order α(m1,m] is defined as in [44]

    L{αu(ˉξ,t)tα}=sαˆu(ˉξ,s)m1j=0sαj1u(j)(ˉξ,0)

    Applying Laplace transform on Eqs (1.1)–(1.3), we have

    pk=1ϱ1,k(sαkˆu(ˉξ,s)sαk1u(ˉξ,0)sαk2ut(ˉξ,0))+ϱ2(sˆu(ˉξ,s)u(ˉξ,0))+qr=1ϱ3,r(sβrˆu(ˉξ,s)sβr1u(ˉξ,0))+ϱ4ˆu(ˉξ,s)=ϱ5Lˆu(ˉξ,s)+ϱ6(sγLˆu(ˉξ,s)sγ1Lu(ˉξ,0))+ˆG(ˉξ,s)

    and

    LBˆu(ˉξ,s)=ˆh(ˉξ,s)

    which implies

    (pk=1ϱ1,ksαkI+ϱ2sI+qr=1ϱ3,rsβrI+ϱ4Iϱ5Lϱ6sγL)ˆu(ˉξ,s)=ˆH(ˉξ,s) (1.4)
    LBˆu(ˉξ,s)=ˆh(ˉξ,s) (1.5)

    where ˆH(ˉξ,s)=pk=1ϱ1,ksαk1g1(ˉξ)+pk=1ϱ1,ksαk2g2(ˉξ)+ϱ2g1(ˉξ)+qr=1ϱ3,rsβr1g1(ˉξ)ϱ6sγ1Lg1(ˉξ)+ˆG(ˉξ,s), L=Δ, and I is the identity operator. In the proposed approach, first we discretize the linear differential operator L via the pseudospectral method, then the full-discrete system (Eqs (1.4) and (1.5)) is solved in Laplace transform space. Finally the solution of Eqs (1.1)–(1.3) is obtained via inversion of the Laplace transform. Our solution method is highly suitable for parallel computations; it belongs to the class of parallel in time methods for TFPDEs (see [46]). The pseudospectral method is discussed in next section.

    The pseudospectral method is an accurate and extremely precise approach for the numerical solution of TFPDEs. The pseudospectral method uses the Chebyshev points for collocation that are nonzero over the entire domain, whereas finite element methods use basis functions that are nonzero only on small subintervals [35]. In the pseudospectral method, the expansion of solutions is characterized by some global basis functions, such as Lagrange's interpolation polynomials. The concept of interpolation and differentiation matrices [47] is very important in the pseudospectral method. To provide a comprehensive description of the pseudospectral method based on Lagrange's interpolation polynomial basis, the main aspects of the differentiation matrices are reviewed in the next section.

    In the pseudospectral method, the solution is considered over [1,1] and the data {(ξj,ˆu(ξj))}mj=0 is interpolated by the Lagrange polynomial (LP) ηj(ξ) of degree m [33,48]

    Im(ξ)=mj=0ηj(ξ)ˆuj

    where ηj(ξ) is the LP at the point ξj(j=0,1,...,m), which is

    ηj(ξ)=(ξξ0)...(ξξj1)(ξξj+1)...(ξξm)(ξjξ0)...(ξjξj1)(ξjξj+1)...(ξjξm) (1.6)

    where ˆuj=ˆu(ξj). The Chebyshev nodes are defined as

    ξj=cos(jπm),j=0,1,...,m (1.7)

    and are used to discretize the domain [1,1]. The 1st derivative ˆu(ξ)ξ is approximated as

    ˆu(ξ)ξΦmˆu

    and the matrix Φm has elements in the following form

    [Φm]i,j=ηj(ξi),i,j=1,2,...,m

    The non-diagonal elements of [Φm]i,j have the following form

    [Φm]i,j=αjαi(ξiξj),ij

    where α1j=mij(ξiξj), and the elements on diagonal entries are expressed as

    [Φm]i,j=mj=0,ji[Φm]i,j,i=0,1,2,..,m

    The elements of the μth-order derivative of Φμm are analytically obtained as

    [Φμm]i,j=ημj(ξi),i,j=1,2,...,m

    In [47,49], more accurate evaluation of the differentiation matrices is elaborated. For the differentiation matrices, Welfert [47] deduced a handy recursion relation as follows

    [Φμm]kj=μξkξj(αjαk[Φ(μ1)m]kk[Φ(μ1)m]kj),kj

    Chebyshev points are remarkable for their ability to enable spectral convergence when appropriate analytic assumptions are made.

    Let us consider the square domain Γ=[1,1]2. Generally, the point in Γ is denoted by ˉξ and is expressed as

    ˉξij=(cos(πim),cos(πjm)),i,j=0,1,2,...,m

    The LPs associated to ˉξij can be written as

    ηij(ˉξ)=ηi(ξ)ηj(ζ)

    where ηij(ˉξij)=δij. The 2nd-order derivatives of the LPs Eq (1.6) are given as

    2ηij(ˉξrs)ξ2=ηi(ξr)ηj(ζs)=[Φ2m]riδjs,
    2ηij(ˉξrs)ζ2=ηi(ξr)ηj(ζs)=δri[Φ2m]sj

    where Φ2m represents the second order differentiation matrix based on Chebyshev points. Applying the operator L on the LPs Eq (1.6) based on the points {ˉξrs} is

    L(ηij(ˉξrs))=([Φ2m]riδjs+δri[Φ2m]sj)

    Consequently, the approximation L via the pseudospectral method is obtained as

    LM=(ImΦ2m+Φ2mIm)

    The conditions in Eq (1.5) are incorporating by considering the interpolation matrix LM and considering all points ˉξ. Furthermore, the rows of LM in correspondence with boundary nodes are replaced with unit vectors that have a one in accordance with the diagonal elements of LM. Hence, the boundary conditions LBˆu(ˉξ,s)=ˆh(ˉξ,s) in Eq (1.5) will be implemented directly [33]. By rearranging the columns and rows of the matrix LM, the following block matrix is obtained

    LΓ=[WR0I]

    where W and I are the nonzero square blocks of order (mmB)×(mmB) and (mB×mB), respectively. However, mB represents the boundary points. The solution of the system Eqs (1.4) and (1.5) can be attained by solving the linear block system

    LΓˆu(ˉξ,s)=[ˆH(ˉξ,s)ˆh(ˉξ,s)]

    where ˆH and ˆh collect the values at interior and boundary collocation points, respectively. In the final step, the inverse Laplace transform is used to obtain the solution of Eqs (1.1)–(1.3).

    The interpolation operator is based on Chebyshev collocation points Eq (1.7) and Lagrange polynomials Eq (1.6) as follows:

    Im:C(Γ)Pm,Im(u)=mj=0u(ξj)ηj(ξ).

    The steps proposed in [50] will be followed for the construction of the interpolation polynomial error bound, for a constant Mm, which satisfies the following estimate

    Im(u)Mmu,uC[1,1] (1.8)

    Additionally,

    Im(u)=u,uPm

    It is possible to show that for Chebyshev interpolation,

    Mm=ln(1+m)π2+1(m+1)

    Based on m, the constant of stability increases very slowly. For uCm+1[1,1], we have the following bound [50]

    uIm(u)12mΓ(m+2)u(m+1) (1.9)

    Theorem 1. [50] If Eqs (1.8) and (1.9) hold, and uCm+1[1,1], we have

    u()Im(u)()2(m)(M()m+1)Γ(m+2)u(m+1),=0,1,2,...,m (1.10)

    with the stability constant

    M()m=1Γ(+1)[MmΓ(m+1)Γ(m+1)]

    Proof: By using Eqs (1.9) and (1.10) for Eq (1.1) in the one-dimension case, the governing operator is written as L=2ξ2; thus, the error estimate is given as

    E=(pk=1ϱ1,kDαktu+ϱ2Dtu+qr=1ϱ3,rDβrtu+ϱ4uϱ5Luϱ6Dγt(Lu))(pk=1ϱ1,kDαktIm(u)+ϱ2Im(Dtu)+qr=1ϱ3,rDβrtIm(u)+ϱ4Im(u)ϱ5LIm(u)ϱ6Dγt(LIm(u)))=pk=1ϱ1,kDαkt(uIm(u))+ϱ2(Dt(uIm(u)))+qr=1ϱ3,rDβrt(uIm(u))+ϱ4(uIm(u))ϱ5L(uIm(u))ϱ6Dγt(L(uIm(u)))pk=1ϱ1,kDαkt(uIm(u))+ϱ2Dt(uIm(u))+qr=1ϱ3,rDβrt(uIm(u))+ϱ4(uIm(u))+ϱ5L(uIm(u))+ϱ6Dγt(L(uIm(u)))pk=1|ϱ1,k|Dαkt(uIm(u))+|ϱ2|Dt(uIm(u))+qr=1|ϱ3,r|Dβrt(uIm(u))+|ϱ4|uIm(u)+|ϱ5|uξξIm(u)ξξ+|ϱ6|Dγt(uξξIm(u)ξξ)Epk=1|ϱ1,k|Dαkt(uIm(u))+|ϱ2|Dt(uIm(u))+qr=1|ϱ3,r|Dβrt(uIm(u))+|ϱ4|12mΓ(m+2)u(m+1)+|ϱ5|2(M(2)m+1)Γ(m)(12)m1u(m+1)+|ϱ6|Dγt(uξξIm(u)ξξ)

    since the time derivatives are accurately computed. Therefore, the error bound of Dαt(uIm(u)),Dt(uIm(u)),Dβt(uIm(u)) and uIm(u) are of the same order and the error bound of Dγt(uξξIm(u)ξξ) and uξξIm(u)ξξ are of the same order. Thus, the following stability bound is obtained:

    ENu(m+1)

    where the constant N is determined by computing the coefficients of u(m+1). A similar stability estimate can be obtained for two-dimensional problems by using the tensor product interpolation operators.

    In this section, we utilize numerical inverse Laplace methods to transform the pseudospectral method solution ˆu(ˉξ,s) from the Laplace domain to the time domain u(ˉξ,t) as follows:

    u(ˉξ,t)=12πiρ+iρiesτˆu(ˉξ,s)ds,ρ>ρ0 (1.11)

    Here, ˆu(ˉξ,s) needs to be inverted and the converging abscissa is ρ0 and ρ>ρ0. Thus, the open half plane Re(s)<ρ contains all the singularities of ˆu(ˉξ,s). Our goal is to approximate the integral defined in Eq (1.11). In most cases, the integral defined in Eq (1.11) is quite difficult to be evaluated analytically; therefore, a numerical method needs to be employed. The integral defined in Eq (1.11) may be evaluated by several numerical algorithms available in literature. Every approach has a distinct application and is appropriate for a certain type of function. In this article, we use the two popular inversion algorithms, the improved Talbot's method [51] and the Stehfest's method [52], which are presented as follows.

    Here, we utilize the Talbot's method to approximate u(ˉξ,t)

    u(ˉξ,t)=12πiρ+iρiestˆu(ˉξ,s)ds=12πiΛestˆu(ˉξ,s)ds,Re(s)>0 (1.12)

    where Λ is a suitably chosen contour. In the Talbot's method, the numerical quadrature is applied to the integral in Eq (1.12). The trapezoidal and midpoint rule are the two effective rules used in conjunction with the deformation of contour [53]. The purpose of contour deformation is to handle the exponential factor. Particularly, the path of integration can be deformed to a Hankel contour, i.e., a contour whose real part starts at in the 3rd quadrant and winds around all the singularities of the transform function going again to in the 2nd quadrant. The exponential component decays rapidly on such contour, making the Bromwich integral appropriate for approximation and using the trapezoidal and midpoint rule. Cauchy's theorem justifies such deformation, provided that the contour remains in the region where the transformed function ˆu(ˉξ,s) is analytic. Furthermore, some mild restrictions are required in the left complex plane on the decay of the transformed function ˆu(ˉξ,s)[54,55]. We consider the contour in the following parametric form as[51]

    Λ:s=s(ω),πωπ,Res(±π)=

    We have

    s(ω)=mTtϖ(ω),ϖ(ω)=η1+η2ωcot(η3ω)+η4iω (1.13)

    where η1,η2,η3, and η4 are to be selected by the user. Plugging Eq (1.13) in Eq (1.12), we have

    u(ˉξ,t)=12πiππes(ω)tˆu(ˉξ,s)(s(ω))s(ω)dω (1.14)

    Midpoint rule with uniform step h=2πmT is utilized to approximate Eq (1.14) as

    uApp(ˉξ,t)1mTimTk=1es(ωk)tˆu(s(ωk))s(ωk),ωk=π+(k12)h (1.15)

    The error analysis of the improved Talbot's approach is based on the following theorem.

    Theorem 1.1. [51] Let ωk be defined as in Eq (1.15). Let g:ΞC be an analytic function in the set

    Ξ={ωC:π<Reω<π,ands<Imω<r}

    when r,s>0, then

    ππg(ω)dω2πmTmTj=1g(ωk)=Ψ(ϕ)+Ψ+(ψ)

    where

    Ψ+(ϕ)=12(π+iϕπ+π+iϕπ+iϕ+ππ+iϕ)(1+itan(mTω2))g(ω)dω

    and

    Ψ(ψ)=12(πiψπ+πiψπiψ+ππiψ)(1itan(mTω2))g(ω)dω

    0<ϕ<r,0<ψ<s,andmTeven if mT is an odd number; we can replace tan(mTω2) with cot(mTω2) if g(ω) is real valued, that is, g(ˉω)=¯g(ω) and if r and s can be taken to be equal, then

    Ψ(ψ)=Ψ+(ψ)+Ψ(ψ)=Reπ+iψπ+iψ(1+itan(mTω2))g(ω)dω

    By examining the complex tangent function's behavior, this may be bounded as

    |Ψ(ψ)|4πCexp(rmT)1

    The above process has been done for even mT and C and r are some positive constants. A similar approach can be used for an odd mT.

    The optimal values of the parameters included in Eq (1.13) can be utilized to find the optimal contour of integration, which is necessary to obtain the best results. The authors of [51] have obtained the optimal values of the parameters as follows

    η1=0.61220,η2=0.50170,η3=0.64070,andη4=0.26450

    The corresponding error estimate is given as

    Errest=| uApp(ˉξ,t)u(ˉξ,t)|=O(e(1.358)mT)

    One of the most effective and straightforward techniques for Laplace transform inversion is the Gaver-Stehfest approach. The latter part of the 1960s saw its design. Because of its simplicity and efficacy, it has become more and more popular in a variety of fields, including computational physics, finance, economics, and chemistry. The basis of the Gaver-Stehfest approach is the series of Gaver approximants, as found by Gaver [56]. Since the convergence of the Gaver approximants was essentially logarithmic, acceleration was necessary. A linear acceleration approach was provided by Stehfest [52] using the Salzer acceleration method. Using a series of functions, the Gaver-Stehfest method approximates u(ˉξ,t) as

    uApp(ˉξ,t)=ln2tmSk=1θkˆu(ˉξ,ln2tk) (1.16)

    where θk are given as

    θk=(1)mS2+kmin(k,mS2)n=k+12nmS2(2n)!(mS2n)!n!(n1)!(kn)!(2nk)! (1.17)

    Solving the system Eqs (1.4) and (1.5) for the corresponding Laplace parameters s=ln2tk,k=1,2,3,...,mS, the approximate solution uApp(ˉξ,t) of the problem in Eq (1.1) can be obtained via Eq (1.16). The Gaver-Stehfest algorithm has a few noteworthy qualities, such as: (ⅰ) u(ˉξ,t) are linear in the context of values of ˆu(ˉξ,s); (ⅱ) the values of ˆu(ˉξ,s) are required merely for real value of s; (ⅲ) the procedure of determining the coefficients is fairly effortless; (ⅳ) for constant functions, this approach results in significantly precise approximations, i.e., if uc,thenuAppc for all mS1. In literature, this methodology has been employed by many researchers in [57,58], in which it is revealed that this strategy converges promptly to uApp(ˉξ,t) (given u(ˉξ,t) is non-oscillatory).

    Convergence

    The convergence of uApp(ˉξ,t) has been derived by the author in [57]. The results are based on the following theorem.

    Theorem 1.2. Let u:(0,)R be a locally integrable function. Let s>0, define the Laplace transform ˆu(ˉξ,s), and let uApp(ˉξ,t) be the numerical solution as given by Eq (1.16).

    1. uApp(ˉξ,t) converges given u(ˉξ,t) near t.

    2. Let for some real number ϕ and 0<ψ<0.25,

    ψ0|u(ˉξ,tlog2(1/2+η))+u(ˉξ,tlog2(1/2η))2ϕ|η1dη<

    then uApp(ˉξ,t)ϕ as mS+.

    3. Let u(ˉξ,t) be of bounded variation near t, then

    uApp(ˉξ,t)u(ˉξ,t+0)+u(ˉξ,t0)2,asmS+

    Corollary 1.3. Using the assumptions of the above theorem, if

    u(ˉξ,t+η)u(ˉξ,t)=O(|η|ϑ)

    η, and some ϑ, then uApp(ˉξ,t)u(ˉξ,t) as mS+.

    Moreover, the authors in [59] conducted a number of experiments to determine how parameters affected the numerical scheme's correctness. Their conclusions are as follows: "For ν1 significant digits, select a mS positive integer 2.2ν1: After setting the system precision to ν2=1.1mS, compute θi,1imS for a given mS using Eq (1.17). Next, compute the uApp(ˉξ,t) in Eq (1.16) for the provided transformed function ˆu(ˉξ,s) and the argument t." These conclusions indicate that the error is 10(ν1+1)uApp(ˉξ,t)u(ˉξ,t)u(ˉξ,t)10ν1, where mS=2.2ν1 [40].

    In this section, three test problems are considered to assess the effectiveness and accuracy of the Laplace transformed pseudospectral method. We compute the maximum absolute error (Err), the relative error (Err2), and the root mean square error (Errrms) between the numerical and the exact solutions, which are defined as follows:

    Err=max1km|u(ˉξk,t)uApp(ˉξk,t)|,
    Err2=mk=1(u(ˉξk,t)uApp(ˉξk,t))2mk=1(u(ˉξk,t))2,
    Errrms=mk=1(u(ˉξk,t)uApp(ˉξk,t))2m

    where uApp(ˉξk,t) is the approximate solution and u(ˉξk,t) is the exact solution.

    Problem 1

    We consider Eq (1.1) with exact solution

    u(ξ,ζ,t)=(t3+1)sin(πξ)sin(πζ)

    and the source term is

    G(ξ,ζ,t)=sin(πξ)sin(πζ)[ϱ1,16Γ(4α)t3α+3ϱ2t2+ϱ3,16Γ(4β)t3β+(ϱ4+2ϱ5π2)(t3+1)+2ϱ6π26Γ(4γ)t3γ]

    where ϱ1,1=ϱ2=ϱ3,1=ϱ4=ϱ5=ϱ6=1 and (ξ,ζ)[1,1]2. The initial-boundary conditions are obtained from the analytical solution. The Err2,Err, and Errrms errors of the proposed method for problem 1 obtained via the Stehfest's and the Talbot's methods by using various values of m,mT,mS,α,β,γ, and t=1 are presented in Tables 1 and 2, respectively. The numerical solution to problem 1 is shown in Figure 1. Figure 2 presents a comparison of error norms Err2,Err, and Errrms computed using the Stehfest's and the Talbot's methods for various values of mS and mT with α=1.5,β=0.7,γ=0.3,m=30, and t=1. The comparison of Err2,Err, and Errrms computed by the Stehfest's and the Talbot's methods for different values of t with α=1.5,β=0.7,γ=0.3,mS=14,mT=30, and m=30 is shown in Figure 3. The comparison of Err2,Err, and Errrms computed using the Stehfest's and the Talbot's methods versus α with β=0.7,γ=0.3,m=30,mS=14,mT=30, and t=1 is presented in Figure 4. The comparison of Err2,Err, and Errrms computed via the Stehfest's and the Talbot's methods versus β with α=1.5,γ=0.7,m=30,mS=14,mT=30, and t=1 is presented in Figure 5. The comparison of Err2,Err, and Errrms computing by the Stehfest's and the Talbot's methods versus γ with α=1.5,β=0.7,m=30,mS=14,mT=30, and t=1 is presented in Figure 6. The plot of Err2 obtained via the Stehfest's and the Talbot's methods for α[1,2] and β[0,1] with γ=0.3,t=1,mS=16,mT=30, and m=30 is shown in Figure 7. In Figure 8, the plots of Err obtained using the Stehfest's and the Talbot's methods for α[1,2] and γ[0,1] with β=0.7,t=1,mS=16,mT=30, and m=30 are shown. Similarly, Figure 9 presents the plots of Errrms obtained using the Stehfest's and the Talbot's methods for α[1,2] and t[0,1] with β=0.7,γ=0.3,mS=16,mT=30, and m=30. It is observed that the accuracy of the Stehfest's method steadily decreases for mS14. From the results presented in tables and figures, we conclude that the proposed method is accurate, stable, and efficient. The numerical results demonstrate that the Talbot's method is more accurate than the Stehfest's method.

    Table 1.  The Err2,Err,Errrms computed by the Stehfest's method for problem 1.
    (α,β,γ) m mS Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 14 4.7216×1004 9.6894×1005 2.1462×1005 0.199250
    25 5.3662×1004 1.0751×1004 2.1465×1005 0.326975
    28 7.5644×1004 1.0128×1004 2.7016×1005 0.556057
    30 10 6.9931×1002 5.9791×1003 2.3310×1003 0.575029
    12 1.5653×1003 1.3408×1004 5.2178×1005 0.663883
    14 7.8051×1004 1.3790×1004 2.6017×1005 0.755536
    (1.6,0.5,0.7) 22 14 4.0428×1004 4.4012×1005 1.8376×1005 0.180055
    25 4.8390×1004 4.4246×1005 1.9356×1005 0.324530
    28 6.4225×1004 8.5592×1005 2.2938×1005 0.576450
    30 10 6.9931×1002 5.9791×1003 2.3310×1003 0.581491
    12 1.5680×1003 1.3409×1004 5.2266×1005 0.688463
    14 7.4051×1004 1.3304×1004 2.4684×1005 0.719391
    (1.9,0.8,0.5) 22 14 3.8530×1004 4.2171×1005 1.7514×1005 0.184399
    25 4.5504×1004 4.3708×1005 1.8202×1005 0.329158
    28 5.2746×1004 6.0001×1005 1.8838×1005 0.556485
    30 10 6.9931×1002 5.9791×1003 2.3310×1003 0.514627
    12 1.5660×1003 1.3402×1004 5.2201×1005 0.630183
    14 6.5840×1004 1.3505×1004 2.1947×1005 0.760430

     | Show Table
    DownLoad: CSV
    Table 2.  The Err2,Err,Errrms computed by the Talbot's method for problem 1.
    (α,β,γ) m mT Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 30 1.4423×1010 3.8758×1011 6.5561×1012 1.431363
    25 2.2738×1010 6.1642×1011 9.0953×1012 2.823492
    28 3.0634×1010 5.6744×1011 1.0941×1011 4.853607
    30 20 5.1049×1006 4.3647×1007 1.7016×1007 4.409938
    25 1.1297×1008 9.6706×1010 3.7658×1010 5.576429
    30 3.7908×1010 8.4039×1011 1.2636×1011 6.435613
    (1.6,0.5,0.7) 22 30 8.8800×1011 1.6472×1011 4.0364×1012 1.432723
    25 1.3312×1010 3.5030×1011 5.3248×1012 2.644643
    28 1.9449×1010 3.9329×1011 6.9462×1012 4.498548
    30 20 5.1049×1006 4.3647×1007 1.7016×1007 4.512210
    25 1.1277×1008 9.6519×1010 3.7589×1010 5.250997
    30 2.5877×1010 8.3971×1011 8.6256×1012 6.315079
    (1.9,0.8,0.5) 22 30 6.3703×1011 1.7695×1011 2.8956×1012 1.465253
    25 9.4756×1011 2.4410×1011 3.7902×1012 2.624655
    28 1.5070×1010 5.0718×1011 5.3821×1012 4.675835
    30 20 5.1049×1006 4.3647×1007 1.7016×1007 4.469602
    25 1.1281×1008 9.6487×1010 3.7602×1010 5.300393
    30 2.0811×1010 4.6520×1011 6.9370×1012 6.311030

     | Show Table
    DownLoad: CSV
    Figure 1.  Numerical solution of problem 1 obtained by proposed scheme.
    Figure 2.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus mS and for α=1.5,β=0.7,γ=0.3,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus mT for α=1.5,β=0.7,γ=0.3,t=1, and m=30.
    Figure 3.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus t for α=1.5,β=0.7,γ=0.3,mS=14, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus t for α=1.5,β=0.7,γ=0.3,mT=30, and m=30.
    Figure 4.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus α for β=0.7,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus α for β=0.7,γ=0.3,mT=30,t=1, and m=30.
    Figure 5.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus β for α=1.5,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus β for α=1.5,γ=0.3,mT=30,t=1, and m=30.
    Figure 6.  (a) Plots of Err2,Err, and Errrms versus γ obtained by Stehfest's method for α=1.5,β=0.7,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms versus γ obtained by Talbot's method for α=1.5,β=0.7,mT=30,t=1, and m=30.
    Figure 7.  (a) Plot of Err2 versus α and β obtained by Stehfest's method for γ=0.3,t=1,mS=16, and m=30. (b) Plot of Err2 versus α and β obtained by Talbot's method for γ=0.3,t=1,mT=24, and m=30.
    Figure 8.  (a) Plot of Err versus α and γ obtained by Stehfest's method for β=0.7,t=1,mS=16, and m=30. (b) Plot of Err versus α and γ obtained by Talbot's method for β=0.7,t=1,mT=24, and m=30.
    Figure 9.  (a) Plot of Errrms versus α and t obtained by Stehfest's method for β=0.7,γ=0.3,mS=16, and m=30. (b) Plot of Errrms versus α and t obtained by Talbot's method for β=0.7,γ=0.3,mT=24, and m=30}.

    Problem 2

    We consider Eq (1.1) with exact solution

    u(ξ,ζ,t)=(t3+1)(1ξ2ζ2)

    and the source term is

    G(ξ,ζ,t)=(1ξ2ζ2)[ϱ1,16Γ(4α)t3α+3ϱ2t2+ϱ3,16Γ(4β)t3β+ϱ4(t3+1)]+4ϱ5(t3+1)+ϱ66Γ(4γ)t3γ

    where ϱ1,1=ϱ2=ϱ3,1=ϱ4=ϱ5=ϱ6=1, (ξ,ζ)[1,1]2 and the initial-boundary conditions are obtained from the analytical solution.

    The Err2,Err, and Errrms errors of the proposed method for problem 2 obtained via the Stehfest's and the Talbot's methods by using various values of m,mT,mS,α,β,γ, and t=1 are presented in Tables 3 and 4, respectively.

    Table 3.  The Err2,Err,Errrms computed by Stehfest's method for problem 2.
    (α,β,γ) m mS Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 14 4.4852×1003 1.6695×1003 2.0387×1004 0.188688
    25 6.8895×1003 1.5726×1003 2.7558×1004 0.394322
    28 1.7303×1002 6.2685×1003 6.1797×1004 0.524445
    30 10 9.4263×1002 5.9793×1003 3.1421×1003 0.568557
    12 2.1084×1003 3.2664×1004 7.0279×1005 0.658355
    14 2.4313×1002 9.5089×1003 8.1042×1004 0.818657
    (1.6,0.5,0.7) 22 14 2.2303×1003 1.0201×1003 1.0138×1004 0.200221
    25 6.1739×1003 2.3026×1003 2.4695×1004 0.316503
    28 9.0251×1003 2.3816×1003 3.2232×1004 0.510553
    30 10 9.4257×1002 5.9791×1003 3.1419×1003 0.550572
    12 2.0884×1003 2.1157×1004 6.9612×1005 0.663605
    14 1.3078×1002 3.9249×1003 4.3594×1004 0.735535
    (1.9,0.8,0.5) 22 14 1.7605×1003 4.8804×1004 8.0021×1005 0.181867
    25 3.2429×1003 1.5817×1003 1.2971×1004 0.318878
    28 7.1727×1003 2.6411×1003 2.5617×1004 0.512775
    30 10 9.4250×1002 5.9791×1003 3.1417×1003 0.585836
    12 2.1968×1003 2.7998×1004 7.3227×1005 0.644230
    14 1.1617×1002 4.9286×1003 3.8722×1004 0.708296

     | Show Table
    DownLoad: CSV
    Table 4.  The Err2,Err,Errrms computed by Talbot's method for problem 2.
    (α,β,γ) m mT Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 30 2.2778×1009 8.4692×1010 1.0354×1010 1.387058
    25 5.2077×1009 2.4683×1009 2.0831×1010 2.893077
    28 5.9240×1009 1.8102×1009 2.1157×1010 4.540104
    30 20 6.8804×1006 4.3647×1007 2.2935×1007 4.290310
    25 1.5504×1008 1.7167×1009 5.1681×1010 5.406541
    30 9.2636×1009 2.3673×1009 3.0879×1010 6.253001
    (1.6,0.5,0.7) 22 30 1.1655×1009 4.5895×1010 5.2979×1011 1.413894
    25 2.4327×1009 8.7358×1010 9.7309×1011 2.551745
    28 5.9798×1009 1.5282×1009 2.1357×1010 4.802169
    30 20 6.8807×1006 4.3647×1007 2.2936×1007 4.319773
    25 1.5690×1008 2.3208×1009 5.2301×1010 5.494691
    30 9.9438×1009 3.3832×1009 3.3146×1010 6.327990
    (1.9,0.8,0.5) 22 30 9.5101×1010 3.0484×1010 4.3228×1011 1.500904
    25 1.6298×1009 5.9318×1010 6.5192×1011 2.605735
    28 3.3339×1009 1.4952×1010 1.1907×1010 4.379358
    30 20 6.8805×1006 4.3647×1007 2.2935×1007 4.246011
    25 1.5398×1008 1.2926×1009 5.1325×1010 5.247654
    30 6.2433×1009 3.0606×1009 2.0811×1010 6.329043

     | Show Table
    DownLoad: CSV

    The numerical solution of problem 2 is shown in Figure 10. Figure 11 presents a comparison of error norms Err2,Err, and Errrms computed via the Stehfest's and the Talbot's methods versus mS and mT with α=1.5,β=0.7,γ=0.3,m=30, and t=1. The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus t with α=1.5,β=0.7,γ=0.3,mS=14,mT=30, and m=30 is shown in Figure 12.

    Figure 10.  Numerical solution of problem 2 obtained by proposed scheme.
    Figure 11.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus mS and for α=1.5,β=0.7,γ=0.3,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus mT for α=1.5,β=0.7,γ=0.3,t=1, and m=30.
    Figure 12.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus t for α=1.5,β=0.7,γ=0.3,mS=14, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus t for α=1.5,β=0.7,γ=0.3,mT=30, and m=30.

    The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus α with β=0.7,γ=0.3,m=30,mS=14,mT=30, and t=1 is presented in Figure 13. The comparison of Err2,Err, and Errrms computed using the Stehfest's and the Talbot's methods versus β with α=1.5,γ=0.7,m=30,mS=14,mT=30, and t=1 is shown in Figure 14.

    Figure 13.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus α for β=0.7,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus α for β=0.7,γ=0.3,mT=30,t=1, and m=30.
    Figure 14.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus β for α=1.5,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus β for α=1.5,γ=0.3,mT=30,t=1, and m=30.

    The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus γ with α=1.5,β=0.7,m=30,mS=14,mT=30, and t=1 is presented in Figure 15. The plots of Err2 obtained using the Stehfest's and the Talbot's methods for α[1,2] and β[0,1] with γ=0.3,t=1,mS=16,mT=30, and m=30 are shown in Figure 16.

    Figure 15.  (a) Plots of Err2,Err, and Errrms versus γ obtained by Stehfest's method for α=1.5,β=0.7,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms versus γ obtained by Talbot's method for α=1.5,β=0.7,mT=30,t=1, and m=30.
    Figure 16.  (a) Plot of Err2 versus α and β obtained by Stehfest's method for γ=0.5,t=1,mS=16, and m=30. (b) Plot of Err2 versus α and β obtained by Talbot's method for γ=0.5,t=1,mT=24, and m=30.

    In Figure 17, the plots of Err obtained using the Stehfest's and the Talbot's methods for α[1,2] and γ[0,1] with β=0.7,t=1,mS=16,mT=30, and m=30 are shown. Similarly, Figure 18 shows the plot of Errrms obtained via the Stehfest's and the Talbot's methods for α[1,2] and t[0,1] with β=0.7,γ=0.3,mS=16,mT=30, and m=30. Comparable performances like the one we witnessed in Problem 2 are observed.

    Figure 17.  (a) Plot of Err versus α and γ obtained by Stehfest's method for β=0.8,t=1,mS=16, and m=30. (b) Plot of Err versus α and γ obtained by Talbot's method for β=0.8,t=1,mT=24, and m=30.
    Figure 18.  (a) Plot of Errrms versus α and t obtained by Stehfest's method for β=0.8,γ=0.5,mS=16, and m=30. (b) Plot of Errrms versus α and t obtained by Talbot's method for β=0.8,γ=0.5,mT=24, and m=30.

    Problem 3

    We consider Eq (1.1) with analytical solution as

    u(ξ,ζ,t)=t2sin(1ξ)(eξ1)sin(1ζ)(eζ1)

    and the source term is

    G(ξ,ζ,t)=sin(1ξ)(eξ1)sin(1ζ)(eζ1)[ϱ1,12Γ(3α)t2α+2ϱ2t+ϱ3,12Γ(3β)t2β+ϱ4t2](ϱ5t2+ϱ62Γ(3γ)t2γ)[sin(ζ1)(eζ1)(2cos(ξ1)eξ+sin(ξ1))+sin(ξ1)(eξ1)(2cos(ζ1)eζ+sin(ζ1))]

    where ϱ1,1=ϱ2=ϱ3,1=ϱ4=ϱ5=ϱ6=1 and ˉξ=(ξ,ζ)[1,1]2. The initial-boundary conditions are obtained from the analytical solution. The Err2,Err, and Errrms errors obtained using the proposed method for problem 3 obtained via the Stehfest's and the Talbot's methods by using various values of m,mT,mS,α,β,γ, and t=1 are presented in Tables 5 and 6, respectively.

    Table 5.  The Err2,Err,Errrms computed by Stehfest's method for problem 3.
    (α,β,γ) m mS Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 14 3.3480×1005 8.6337×1006 1.5218×1006 0.254196
    25 9.7449×1005 3.7684×1005 3.8980×1006 0.363822
    28 2.3165×1004 1.0850×1004 8.2733×1006 0.585754
    30 10 2.2116×1004 1.8488×1005 7.3719×1006 0.542104
    12 1.7804×1004 1.4958×1005 5.9348×1006 0.687321
    14 3.1175×1004 7.7178×1005 1.0392×1006 0.802970
    (1.6,0.5,0.7) 22 14 3.0799×1005 1.0421×1005 1.4000×1006 0.191420
    25 8.6136×1005 3.5275×1005 3.4454×1006 0.339131
    28 1.4053×1004 4.7890×1005 5.0190×1006 0.559366
    30 10 2.2101×1004 1.8488×1005 7.3670×1006 0.583434
    12 1.8068×1004 1.4953×1005 6.0227×1006 0.673867
    14 1.2604×1004 4.4833×1005 4.2014×1006 0.733930
    (1.9,0.8,0.5) 22 14 2.0609×1005 5.6164×1006 9.3675×1007 0.185813
    25 5.6168×1005 1.6845×1005 2.2467×1006 0.337030
    28 4.7054×1005 2.1159×1005 1.6805×1006 0.513103
    30 10 2.2103×1004 1.8488×1005 7.3676×1006 0.648050
    12 1.7979×1004 1.4953×1005 5.9930×1006 0.627418
    14 1.3026×1004 4.7709×1005 4.3420×1006 0.718273

     | Show Table
    DownLoad: CSV
    Table 6.  The Err2,Err,Errrms computed by Talbot's method for problem 3.
    (α,β,γ) m mT Err2 Err Errrms CPU (sec.)
    (1.3,0.2,0.9) 22 30 2.7733×1011 1.2506×1011 1.2606×1012 1.406847
    25 6.7413×1011 2.5993×1011 2.6965×1012 2.649332
    28 1.5151×1011 4.4017×1011 5.4112×1012 4.415065
    30 20 6.8451×1008 5.7242×1009 2.2817×1009 4.279803
    25 1.6659×1010 3.2206×1011 5.5531×1012 5.292333
    30 1.5071×1010 4.2789×1011 5.0236×1012 6.272503
    (1.6,0.5,0.7) 22 30 2.4943×1011 1.2897×1011 1.1338×1012 1.503324
    25 5.4494×1011 1.8470×1011 2.1798×1012 2.687497
    28 1.1024×1011 3.6371×1011 3.9370×1012 4.498761
    30 20 6.8453×1008 5.7242×1009 2.2818×1009 4.177472
    25 1.4326×1010 3.6404×1011 4.7752×1012 5.148521
    30 1.3951×1010 6.3866×1011 4.6503×1012 6.443956
    (1.9,0.8,0.5) 22 30 7.6747×1012 2.4080×1012 3.4885×1013 1.476408
    25 3.0024×1011 1.6923×1011 1.2010×1012 2.663378
    28 4.6899×1011 1.8591×1011 1.6750×1012 4.540820
    30 20 6.8448×1008 5.7242×1009 2.2816×1009 4.224953
    25 1.3981×1010 3.3329×1011 4.6602×1012 5.341039
    30 5.1002×1011 2.4527×1011 1.7001×1012 6.258475

     | Show Table
    DownLoad: CSV

    The numerical method solution to problem 3 is shown in Figure 19. Figure 20 represents a comparison of error norms Err2,Err, and Errrms computed by the Stehfest's and the Talbot's methods versus mS and mT with α=1.5,β=0.7,γ=0.3,m=30, and t=1. The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus t with α=1.5,β=0.7,γ=0.3,mS=14,mT=30, and m=30 is shown in Figure 21.

    Figure 19.  Numerical solution of problem 3 obtained by proposed scheme.
    Figure 20.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus mS and for α=1.5,β=0.7,γ=0.3,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus mT for α=1.5,β=0.7,γ=0.3,t=1, and m=30.
    Figure 21.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus t for α=1.5,β=0.7,γ=0.3,mS=14, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus t for α=1.5,β=0.7,γ=0.3,mT=30, and m=30.

    The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods for different values of α with β=0.7,γ=0.3,m=30,mS=14,mT=30, and t=1 is shown in Figure 22. The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods for different values of β with α=1.5,γ=0.7,m=30,mS=14,mT=30, and t=1 is presented in Figure 23.

    Figure 22.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus α for β=0.7,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus α for β=0.7,γ=0.3,mT=30,t=1, and m=30.
    Figure 23.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus β for α=1.5,γ=0.3,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus β for α=1.5,γ=0.3,mT=30,t=1, and m=30.

    The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus γ with α=1.5,β=0.7,m=30,mS=14,mT=30, and t=1 is presented in Figure 24. The plot of Err2 using the Stehfest's and the Talbot's methods for α[1,2] and β[0,1] with γ=0.3,t=1,mS=16,mT=30, and m=30 is shown in Figure 25.

    Figure 24.  (a) Plots of Err2,Err, and Errrms versus γ obtained by Stehfest's method for α=1.5,β=0.7,mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms versus γ obtained by Talbot's method for α=1.5,β=0.7,mT=30,t=1, and m=30.
    Figure 25.  (a) Plot of Err2 versus α and β obtained by Stehfest's method for γ=0.5,t=1,mS=16, and m=30. (b) Plot of Err2 versus α and β obtained by Talbot's method for γ=0.5,t=1,mT=24, and m=30.

    In Figure 26, the plots of Err obtained using the Stehfest's and the Talbot's methods for α[1,2] and γ[0,1] with β=0.7,t=1,mS=16,mT=30, and m=30 are shown. Similarly, Figure 27 shows the plots of Errrms using the Stehfest's and the Talbot's methods for α[1,2] and t[0,1] with β=0.7,γ=0.3,mS=16,mT=30, and m=30.

    Figure 26.  (a) Plot of Err versus α and γ obtained by Stehfest's method for β=0.8,t=1,mS=16, and m=30. (b) Plot of Err versus α and γ obtained by Talbot's method for β=0.8,t=1,mT=24, and m=30.
    Figure 27.  (a) Plot of Errrms versus α and t obtained by Stehfest's method for β=0.8,γ=0.5,mS=16, and m=30. (b) Plot of Errrms versus α and t obtained by Talbot's method for β=0.8,γ=0.5,mT=24, and m=30.

    It is observed that the accuracy of the Stehfest's method steadily decreases for mS14. From the results presented in tables and figures, we conclude that the proposed method is accurate, stable, and efficient. The numerical results demonstrates that the Talbot's method is more accurate than the Stehfest's method.

    Problem 4

    We consider a limiting case of Eq (1.1) with analytical solution as

    u(ξ,ζ,t)=Eγ(π22tγ)cos(π2ξ)cos(π2ζ)

    and the source term is

    G(ξ,ζ,t)=Eγ(π22tγ)cos(π2ξ)cos(π2ζ)

    where ϱ1,1=ϱ3,1=ϱ2=0, ϱ6=ϱ4=ϱ5=1, (ξ,ζ)[1,1]2, and the initial-boundary conditions are obtained from the analytical solution. Eγ(t)=n=0tnΓ(nγ+1),(γC,Re(γ)>0) is the Mittag Leffler function. The Err2,Err, and Errrms errors of the proposed method for problem 4 by using various values of m,mT,mS,γ, and t=1 are presented in Tables 7 and 8, respectively. The numerical solution of problem 4 is shown in Figure 28. Figure 29 presents a comparison of error norms Err2,Err, and Errrms obtained via the Stehfest's and the Talbot's methods versus mS and mT with γ=0.7,m=30, and t=1. The comparison of Err2,Err, and Errrms obtained using the Stehfest's and the Talbot's methods versus t with γ=0.7,mS=14,mT=25, and m=30 is shown in Figure 30. The comparison of Err2,Err, and Errrms using the Stehfest's and the Talbot's methods versus γ with m=30,mS=14,mT=25, and t=1 is presented in Figure 31. The plot of Err2 using the Stehfest's and the Talbot's methods for γ[0.5,1] and t[0,1] with mS=14,mT=25, and m=30 is shown in Figure 32. In Figure 33, the plots of Err computed using the Stehfest's and the Talbot's methods with γ[0.5,1],t[0,1],mS=14,mT=25, and m=30 are shown. Similarly, Figure 34 shows the plot of Errrms obtained using the Stehfest's and the Talbot's methods with γ[0.5,1],t[0,1],mS=14,mT=25, and m=30. The computational results undeniably demonstrate that Talbot's method is greater in precision than Stehfest's method.

    Table 7.  The Err2,Err,Errrms computed by Stehfest's method for problem 4.
    m mS Err2 Err Errrms CPU (sec.)
    γ=0.5 21 14 7.5454×1005 1.0208×1005 3.5930×1006 0.155965
    24 8.7095×1005 1.0340×1005 3.6290×1006 0.263790
    27 9.6493×1005 1.0212×1005 3.5738×1006 0.475612
    30 08 3.3892×1005 3.2475×1006 1.1297×1006 0.433870
    10 1.1589×1004 1.1104×1005 3.8630×1006 0.562005
    12 1.1026×1004 1.0576×1005 3.6754×1006 0.658822
    γ=0.7 21 14 8.2377×1006 1.7617×1006 3.9227×1007 0.156614
    24 7.5928×1006 1.5342×1006 3.1637×1007 0.263522
    27 1.4674×1005 3.3973×1006 5.4349×1007 0.459460
    30 08 2.2064×1004 2.1141×1005 7.3547×1006 0.428611
    10 2.3103×1004 2.2138×1005 7.7011×1006 0.567407
    12 3.1044×1005 2.9706×1006 1.0348×1006 0.623242
    γ=0.9 21 14 9.4376×1005 1.2762×1005 4.4941×1006 0.160795
    24 1.0828×1004 1.2942×1005 4.5116×1006 0.288674
    27 1.2360×1004 1.2863×1005 4.5779×1006 0.424751
    30 08 1.2649×1003 1.2120×1004 4.2162×1005 0.462303
    10 1.8471×1003 1.7699×1004 6.1570×1005 0.521154
    12 6.5184×1004 6.2457×1005 2.1728×1005 0.645028

     | Show Table
    DownLoad: CSV
    Table 8.  The Err2,Err,Errrms computed by Talbot's method for problem 4.
    m mT Err2 Err Errrms CPU (sec.)
    γ=0.5 21 30 7.5633×1005 1.0211×1005 3.6016×1006 0.565917
    25 9.0040×1005 1.0253×1005 3.6016×1006 1.337206
    29 1.0445×1004 1.0278×1005 3.6016×1006 2.600535
    30 10 1.0893×1004 1.0438×1005 3.6311×1006 1.103470
    12 1.0800×1004 1.0348×1005 3.6000×1006 1.386152
    14 1.0804×1004 1.0353×1005 3.6015×1006 1.538342
    γ=0.7 21 30 3.2260×1011 4.3311×1012 1.5362×1012 0.598598
    25 3.8574×1011 4.3535×1012 1.5430×1012 1.348554
    29 4.3999×1011 4.2745×1012 1.5172×1012 2.666525
    30 12 7.2355×1008 6.9329×1009 2.4118×1009 1.266167
    16 2.8863×1010 2.7649×1011 9.6211×1012 1.678309
    20 4.5249×1011 4.3492×1012 1.5083×1012 2.078595
    γ=0.9 21 30 5.4576×1012 1.7432×1012 2.5989×1013 0.579733
    25 6.9966×1012 1.5988×1012 2.7987×1013 1.395903
    29 1.3880×1011 4.5894×1012 4.7862×1013 2.534978
    30 12 8.2775×1008 7.9314×1009 2.7592×1009 1.266089
    16 3.3689×1010 3.2289×1011 1.1230×1011 1.622688
    20 3.5520×1012 8.5037×1013 1.1840×1013 1.994124

     | Show Table
    DownLoad: CSV
    Figure 28.  Numerical solution of problem 4 obtained by proposed scheme.
    Figure 29.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus mS and for γ=0.7,t=1, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus mT for γ=0.7,t=1, and m=30. It can be seen that the observed errors are in good agreement with theoretical error.
    Figure 30.  (a) Plots of Err2,Err, and Errrms obtained by Stehfest's method versus t for γ=0.7,mS=14, and m=30. (b) Plots of Err2,Err, and Errrms obtained by Talbot's method versus t for γ=0.7,mT=25, and m=30.
    Figure 31.  (a) Plots of Err2,Err, and Errrms versus γ obtained by Stehfest's method for mS=14,t=1, and m=30. (b) Plots of Err2,Err, and Errrms versus γ obtained by Talbot's method for mT=25,t=1, and m=30.
    Figure 32.  (a) Plot of Err2 versus γ and t obtained by Stehfest's method for mS=14, and m=30. (b) Plot of Err2 versus γ and t obtained by Talbot's method for mT=25, and m=30.
    Figure 33.  (a) Plot of Err versus γ and t obtained by Stehfest's method for mS=14, and m=30. (b) Plot of Err versus γ and t obtained by Talbot's method for mT=25, and m=30.
    Figure 34.  (a) Plot of Errrms versus γ and t obtained by Stehfest's method for mS=14, and m=30. (b) Plot of Errrms versus γ and t obtained by Talbot's method for mT=25, and m=30.

    In the current work, we have developed an efficient and stable numerical method, which combines the numerical Laplace transform method in time with the pseudospectral method in space for the numerical solution of the two-dimensional time fractional multi-term mixed sub-diffusion and diffusion-wave equation. The proposed method offers an excellent approach to the solution process of the considered equation. In our technique, first, the Laplace transform was employed, which reduced the problem into an elliptic problem in the Laplace transform domain, then the pseudospectral method was utilized to obtain the approximate solution to the transformed problem. Finally, the improved Talbot's method and the Stehfest's method were used to convert the obtained solution in Laplace transform domain back into time domain. The improved Talbot's method and the Stehfest's method provide a numerical inversion process that is accurate, stable, easy to implement, and does not suffer from stability issues which occur with finite difference time stepping methods. From the results presented in tables and figures, it can be observed that the Stehfest's method provides optimal results for mS<14, and for mS14, its accuracy decreases. Furthermore, the obtained numerical results clearly demonstrate that the Talbot's method is considerably more precise than the Stehfest's method.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors D. Santina and N. Mlaiki would like to thank the Prince Sultan University for paying the publication fees for this work through TAS LAB. This research is supported by Prince Sultan University, Saudi Arabia.

    The authors declare that they have no conflicts of interest.



    Conflict of interest



    The authors declare no conflict of interest.

    [1] Alkawareek MY, Gorman SP, Graham WG, et al. (2014) Potential cellular targets and antibacterial efficacy of atmospheric pressure non-thermal plasma. Int J Antimicrob Ag 43: 154-160.
    [2] Dobrynin D, Fridman G, Friedman G, et al. (2009) Physical and biological mechanisms of direct plasma interaction with living tissue. New J Phys 11: 115020.
    [3] Flynn PB, Busetti A, Wielogorska E, et al. (2016) Non-thermal plasma exposure rapidly attenuates bacterial AHL-dependent quorum sensing and virulence. Sci Rep 6: 26320.
    [4] Martusevich AK, Solov'eva AG, Galka AG, et al. (2019) Effects of helium cold plasma on erythrocyte metabolism. B Exp Biol Med 167: 198-200.
    [5] Wiegand C, Fink S, Beier O, et al. (2016) Dose-and time-dependent cellular effects of cold atmospheric pressure plasma evaluated in 3D skin models. Skin Pharmacol Phys 29: 257-265.
    [6] Ermolaeva SA, Varfolomeev AF, Chernukha MY, et al. (2011) Bactericidal effects of non-thermal argon plasma in vitro, in biofilms and in the animal model of infected wounds. J Med Microbiol 60: 75-83.
    [7] Lee C, Subhadra B, Choi HG, et al. (2019) Inactivation of mycobacteria by radicals from non-thermal plasma jet. J Microbiol Biotechn 29: 1401-1411.
    [8] Ahn HJ, Kim KI, Kim G, et al. (2011) Atmospheric-pressure plasma jet induces apoptosis involving mitochondria via generation of free radicals. PloS one 6: e28154.
    [9] Brun P, Pathak S, Castagliuolo I, et al. (2014) Helium generated cold plasma finely regulates activation of human fibroblast-like primary cells. PloS one 9: e104397.
    [10] Jawaid P, Rehman MU, Zhao QL, et al. (2016) Helium-based cold atmospheric plasma-induced reactive oxygen species-mediated apoptotic pathway attenuated by platinum nanoparticles. J Cell Mol Med 20: 1737-1748.
    [11] Szili EJ, Hong SH, Oh JS, et al. (2018) Tracking the penetration of plasma reactive species in tissue models. Trends Biotechnol 36: 594-602.
    [12] Yan D, Cui H, Zhu W, et al. (2017) The strong cell-based hydrogen peroxide generation triggered by cold atmospheric plasma. Sci Rep 7: 10831.
    [13] Lackmann JW, Schneider S, Edengeiser E, et al. (2013) Photons and particles emitted from cold atmospheric-pressure plasma inactivate bacteria and biomolecules independently and synergistically. J R Soc Interface 10: 20130591.
    [14] Martusevich AA, Solov'Ieva AG, Martusevich AK (2013) Influence of singlet oxygen inhalation on the state of blood pro-and antioxidant systems and energy metabolism. B Exp Biol Med 156: 41-43.
    [15] Soloveva AG, Peretyagin SP (2016) The effect of subchronic inhalations of nitric oxide on metabolic processes in blood of experimental animals. Biomeditsinskaia Khimiia 62: 212-214.
    [16] Kim SY, Lee SY, Min SC (2019) Improvement of the antioxidant activity, water solubility, and dispersion stability of prickly pear cactus fruit extracts using argon cold plasma treatment. J Food Sci 84: 2876-2882.
    [17] Sheteiwy MS, An J, Yin M, et al. (2019) Cold plasma treatment and exogenous salicylic acid priming enhances salinity tolerance of Oryza sativa seedlings. Protoplasma 256: 79-99.
    [18] Peretyagin SP, Martusevich AK, Solovyeva AG, et al. (2013) Enzymological evaluation of hepatotropic effect of ozone in a subchronic experiment. B Exp Biol Med 154: 789-791.
    [19] Smith NL, Wilson AL, Gandhi J, et al. (2017) Ozone therapy: an overview of pharmacodynamics, current research, and clinical utility. Med Gas Res 7: 212-219.
    [20] Zanardi I, Borrelli E, Valacchi G, et al. (2016) Ozone: a multifaceted molecule with unexpected therapeutic activity. Curr Med Chem 23: 304-314.
    [21] Zeng J, Lu J (2018) Mechanisms of action involved in ozone-therapy in skin diseases. Int Immunopharmacol 56: 235-241.
    [22] Vanin AF (2019) What is the mechanism of nitric oxide conversion into nitrosonium ions ensuring S-nitrosating processes in living organisms. Cell Biochem Biophys 77: 279-292.
  • This article has been cited by:

    1. Ujala Gul, Zareen A. Khan, Salma Haque, Nabil Mlaiki, A Local Radial Basis Function Method for Numerical Approximation of Multidimensional Multi-Term Time-Fractional Mixed Wave-Diffusion and Subdiffusion Equation Arising in Fluid Mechanics, 2024, 8, 2504-3110, 639, 10.3390/fractalfract8110639
    2. Yubing Jiang, Hu Chen, Error estimate of GL‐ADI scheme for 2D multiterm nonlinear time‐fractional subdiffusion equation, 2024, 47, 0170-4214, 14588, 10.1002/mma.10290
    3. Salma Aljawi, Ahmad Aloqaily, Nabil Mlaiki, A local meshless method for the numerical solution of multi-term time-fractional generalized Oldroyd-B fluid model, 2024, 10, 24058440, e34061, 10.1016/j.heliyon.2024.e34061
    4. Jia Tan, Tian-jun Wang, Spectral collocation method coupled with domain decomposition for exterior problems of the Fisher equation, 2024, 175, 08981221, 289, 10.1016/j.camwa.2024.10.003
    5. Zeeshan Ali, Kamal Shah, Thabet Abdeljawad, Amjad ali, Analyzing wave dynamics of Burger-Poisson fractional partial differential equation, 2025, 26668181, 101153, 10.1016/j.padiff.2025.101153
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3914) PDF downloads(155) Cited by(4)

Figures and Tables

Figures(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog