亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Stability of Linear θ-Method for Delay Partial Functional Differential Equations with Neumann Boundary Conditions

        2022-01-19 06:23:56CHENYongtang陳永堂WANGQi王琦
        應(yīng)用數(shù)學(xué) 2022年1期
        關(guān)鍵詞:王琦

        CHEN Yongtang(陳永堂), WANG Qi(王琦)

        ( School of Mathematics and Statistics, Guangdong University of Technology,Guangzhou 510006, China)

        Abstract: This paper is mainly concerned with the numerical stability of delay partial functional differential equations with Neumann boundary conditions.Firstly, the sufficient condition of asymptotic stability of analytic solutions is obtained.Secondly, the linear θ-method is applied to discretize the above mentioned equation, and the stability of the numerical solutions is discussed for different ranges of parameter θ.Compared with the corresponding equation with Dirichlet boundary conditions, our results are more intuitive and easier to verify.Finally, some numerical examples are presented to illustrate our theoretical results.

        Key words: Delay partial functional differential equation; Neumann boundary condition;Linear θ-method; Asymptotic stability

        1.Introduction

        Delay partial functional differential equations (DPFDEs) have been widely applied because of the establishment of models to solve all kinds of problems in nature, science and life[1-5].In spite of this, there are some difficulties in obtaining the analytic solutions[6].So many researchers are constantly looking for precise numerical methods to solve DPFDEs.Among many properties for numerical solutions, stability and convergence are very important,which have aroused wide attention in recent years.Mead and Zubik-kowal[7]investigated delay partial differential equations by a Jacobi waveform relaxation method based on Chebyshev pseudo-spectral discretization, whose convergence speed gradually accelerated when parameteraincreased from 0 to 1.Meanwhile, the error of method decreased with the increase ofaand reached the minimum ata= 1.In [8], an absolutely stable difference scheme with first order accuracy was designed for the study of third-order delayed partial differential equations, while the stability estimate of the three-step difference scheme was obtained.At the same time, the authors continuously researched the existence and uniqueness of bounded solutions for semilinear DPFDEs with the same method.Similarly,ZHAO et al.[9]applied the explicit exponential Runge-Kutta methods to study semilinear DPFDEs.Stiffconvergence and conditional DN-stability of these methods were investigated and the conditions for stiffconvergence order were derived up to fourth order.In [10], Sumit and Kuldeep employed the fitted operator finite difference scheme to discretize the spatial variable and the time variable of the singularly perturbed parabolic differential equations with a time delay and the accuracy in the spatial direction was improved into second order.In addition, Bashier and Patidar[11]designed a robust fitted operator finite difference method for the numerical solution of singularly perturbed parabolic differential equations.This method was unconditionally stable and convergent with orderO(k+h2), wherekandhwere the time and space step-sizes, respectively.Compared with [10-11], Kumar and Kumar[12]used high order parameter-uniform discretization to investigate singularly perturbed parabolic partial differential equations with time delay, which was second order accuracy in time and fourth order accuracy in space.In[13-14], forward and backward Euler scheme and Crank-Nicolson scheme were used to discretize parabolic differential equations with a constant delay, respectively.It was shown that backward Euler scheme and Crank-Nicolson scheme were unconditionally delay-independently asymptotically stable while the forward Euler scheme required an additional restriction on the time and spatial stepsizes.Adam et al.[15]constructed a fitted Galerkin spectral method to solve DPFDEs.They established error estimates for the fully discrete scheme with a Galerkin spectral approximation in space and they found that the numerical solutions were very similar to those obtained by Chebyshev pseudo-spectral method.In [16], a linearized compact difference scheme was presented for a class of nonlinear delay partial differential equations with Dirichlet boundary conditions.The unique solvability, unconditional convergence and stability were discussed.

        Different from them, in this study, for the DPFDEs with Neumann boundary condition,we will investigate the stability of numerical solution, which is further extension of [17].By contrast,the stability conditions of linearθ-method in our work are simpler and easier to test,which is the main difference between the two types of boundary conditions.

        This paper is organized as follows.In Section 2, a sufficient condition of asymptotic stability of the analytic solutions is given.In Section 3, linearθ-method is used to discretize DPFDEs and the compact form is obtained.Section 4 is devoted to the stability of linearθ-method withθin different intervals.Finally, some numerical experiments are presented to validate the theoretical results in Section 5 and a brief conclusion is given in Section 6.

        2.Stability of DPFDEs

        In this section, we will give a sufficient condition of asymptotic stability for the analytic solution of the following problem

        here parametersr1>0,r2>0, diffusion coefficientsr3∈R,r4∈R, andτ >0 is the delay term.

        Definition 2.1The analytic solutionu(x,t)≡0 of Problem (2.1) is called asymptotically stable if its solutionu(x,t) close to a sufficiently differentiable functionu0(x,t) withsatisfies limt→∞u(x,t)=0.

        Theorem 2.1Assume that the solution of Problem (2.1) isu(x,t) = eλtcos(nx),whereλ ∈C,n ∈N+,x ∈[0,π] andt >0.Then the analytic solution of Problem (2.1) is asymptotically stable for

        and unstable for

        ProofLetX=B[0,π]be the Banach space equipped with the maximum norm.DefineD(A)={y ∈X:∈X,0)=(π)=0}andAy=fory ∈D(A).

        Let?r1n2(n= 1,2,···) be the eigenvalues ofA.According to Theorem 3 in [17], if all zeros of the following characteristic equations

        have negative real part,then the analytic solution is asymptotically stable.At the same time,if at least one zero has positive real part, then it is unstable.

        Denoteλ=u+vi,u,v ∈R.Letf(λ)=0, thenu+vi=r3?r1n2+(r4?r2n2)e?uτ?vτi.Separating the real and imaginary parts, we arrive at

        Now, assume the real partu <0, thenr3?r1n2< ?(r4?r2n2)e?uτcosvτ.Due tov=?(r4?r2n2)e?uτsinvτhave some roots in R, then?(r4?r2n2)>0.So we haver3?r1n2

        If(2.2)holds,then all zeros of the characteristic equations have negative real part.Therefore, the analytic solution of Problem (2.1) is asymptotically stable.

        If (2.3) holds, then there exists a zeroλ0with positive real part such thatf(λ0) = 0,which implies the analytic solution is unstable.

        Thus, the proof is completed.

        3.Linear θ-Method

        Let Δx= π/N,N ∈N+be the space step size and Δt= 1/m,m ∈N+be the time step size, respectively.Denotexj=jΔx(j= 0,1,2,··· ,N),tk=kΔt(k= 0,1,2,···), thenxjandtkconstitute a uniform space-time grid diagram.Letbe the numerical approximation ofu(xj,tk) and the linearθ-method of Problem (2.1) can be defined as:

        From the first part of (3.1), we have

        where

        Thus, (3.1) can be written as:

        where

        4.Asymptotic Stability of Linear θ-Method

        In this section, we will discuss the asymptotic stability of the numerical solution of Problem (2.1).

        Definition 4.1A numerical method applied to Section 3 of Problem (2.1) is called asymptotically stable ifclose to a sufficiently differentiable functionu0(x,t)with=0 satisfies

        According to Section 3 of [17], if we want to verify the numerical solution being asymptotically stable, we need to prove that

        is a Schur polynomial for anym ≥1, whereλjis thej-th eigenvalue of the matrixSandλj=2 cos(jΔx),j=0,1,2,··· ,N.[18]

        Therefore, it is necessary to introduce a corresponding lemma that we can prove (4.2) is a Schur polynomial .

        Lemma 4.1[17]Letγm(z) =α(z)zm ?β(z) be a polynomial, whereα(z) andβ(z) are polynomials of constant degree.Thenγm(z) is a Schur polynomial for anym ≥1 if and only if the following conditions hold

        (i)α(z)=0?|z|<1;

        (ii)|β(z)|≤|α(z)|,?z ∈C,|z|=1;

        (iii)γm(z)0,?z ∈C,|z|=1,?m ≥1.

        Denote

        then

        Theorem 4.1Under the condition(2.2),suppose that 4e1>e3,4e2>e4, e3<0,e4<0 and?e3>4e2?e4.Then (3.1) is asymptotically stable forθ ∈[0if and only if

        whereei(i=1,2,3,4) are defined in (3.2).

        Proof(Necessity) First of all, we prove the item (i) of Lemma 4.1.Fromφj(z)=0 we can derive

        In order to verify|βj(z)| ≤|αj(z)|,?z ∈C, |z|= 1, j= 0,1,2,··· ,N, we set up a function defined in complex domain

        Lettingw=x+yi, after some basic simplifications, we obtain

        In a word, whereas after (a) and (b) have been discussed, we conclude that, for allz ∈C,|z|=1,φj(z)>ψj(z),j=0,1,2,··· ,Nholds.At this point, the items (ii) and (iii) of Lemma 4.1 have been proved.

        By the above analysis, we conclude that (3.1) is asymptotically stable with the help of Lemma 4.1.

        (Sufficiency) We prove in two opposite ways:

        (a) If (4e1+4e2?e3?e4)(1?2θ)=2, we takemto be odd,j=Nandz=?1.Then for|z|=1, we observe that1)=0, which is in conflict with condition (iii) of Lemma 4.1.Thus, (3.1) is not asymptotically stable.

        (b) Assume that (4e1+4e2?e3?e4)(1?2θ)>2.Letmbe odd,j=Nandz=?1,then we get|ψj(?1)| > |φj(?1)|.This confirms that condition (ii) of Lemma 4.1 does not hold, so (3.1) is not asymptotically stable.

        Therefore, (4.5) is a necessary condition for asymptotic stability.This proof is finished.

        Next, we will prove that (3.1) is unconditionally asymptotically stable forθ ∈

        Theorem 4.2Under the condition(2.2),suppose that 4e1>e3,4e2>e4,e3<0,e4<0 and?e3>4e2?e4.Then (3.1) is unconditionally asymptotically stable for

        ProofFirst of all, due toφj(z)=0, we arrive at

        which is similar to the proof of Theorem 4.1.

        Then, in order to verify the items (ii) and (iii) of Lemma 4.1, we define the following function

        (a)θ=Setw=x+yi and|z|=1.After some reductions, we get

        By the conditions 4e1> e3,4e2> e4,e3<0,e4<0,?e3>4e2?e4and 4e1?e3≥e1(2?λj)?e3≥?e3>4e2?e4≥e2(2?λj)?e4≥?e4>0, we find that

        which confirms the items (ii) and (iii) of Lemma 4.1.

        (b)].In the same way, we arrive at

        In this case, for allz ∈C,|z|=1, we also obtain that

        This completes the proof.

        5.Numerical Experiments

        In this section, some numerical experiments are carried out to illustrate the theoretical results.Consider the following problem:

        Fig.5.1 The numerical solutions of Problem (5.1) with m=16 and θ =

        Fig.5.2 The numerical solutions of Problem(5.1) with m=25 and θ =

        Fig.5.3 The numerical solutions of Problem (5.1) with m=30 and θ =

        Fig.5.4 The numerical solutions of Problem(5.1) with m=14 and θ =

        Fig.5.5 The numerical solutions of Problem (5.1) with m=15 and θ =

        Fig.5.6 The numerical solutions of Problem(5.1) with m=5 and θ =0.5

        Fig.5.7 The numerical solutions of Problem (5.1) with m=8 and θ =0.6

        Fig.5.8 The numerical solutions of Problem(5.1) with m=15 and θ =0.6

        Fig.5.9 The numerical solutions of Problem (5.1) with m=25 and θ =0.8

        Fig.5.10 The numerical solutions of Problem(5.1) with m=50 and θ =0.8

        6.Conclusion

        The linearθ-method for solving DPFDEs with Neumann boundary conditions is proposed in this paper.The asymptotic stability condition of the analytic solutions and the numerical solutions are derived, respectively.Compared with DPFDEs with Dirichlet boundary conditions[17], it is shown that the stability condition is more intuitive and effective.In our future work, we will consider the multidimensional problem.

        猜你喜歡
        王琦
        Stability and Convergence of Non-standard Finite Difference Method for Space Fractional Partial Differential Equation
        Angle robust transmitted plasmonic colors with different surroundings utilizing localized surface plasmon resonance
        常見曲線的參數(shù)方程及其應(yīng)用
        Ultra-broadband absorber based on cascaded nanodisk arrays
        Briefly Talking About Methods Of Infiltrating Mental Health Education In Ideological And Political Teaching
        Pf- D mrt4, a potential factor in sexual development in the pearl oyster Pinctada f ucata*
        Comparison ofintestinal microbiota and activities of digestive and immune-related enzymes of sea cucumberApostichopus japonicusin two habitats*
        《皇帝的新裝》后傳
        源于現(xiàn)實之上的詩性想象
        Experimental study of flow patterns and pressure drops of heavy oil-water-gas vertical flow*
        久久精品国产91久久性色tv| 一本色道无码不卡在线观看| 亚洲色成人网站www永久| 国产福利永久在线视频无毒不卡| 国产成人久久精品一区二区三区| 中文字幕一区二区三区人妻少妇| 欧美精品在线一区| 亚洲乱码中文字幕综合69堂| 中文字幕影片免费人妻少妇| 亚洲精品无码久久久久去q| 亚洲午夜精品久久久久久人妖| 亚洲av福利天堂在线观看| 亚洲国产区中文在线观看 | 中文字幕一区二区人妻秘书| 国产精品久久久国产盗摄| 成年人黄视频大全| 国产精品国产三级国av在线观看| 国产精品午夜高潮呻吟久久av| 免费国产在线视频自拍白浆| 丰满少妇a级毛片| 成年女人毛片免费观看97| 亚欧同人精品天堂| 亚洲岛国一区二区三区| 少妇愉情理伦片高潮日本| 就去吻亚洲精品欧美日韩在线| 日本一区不卡高清在线观看| 日本淫片一区二区三区| 国产猛男猛女超爽免费视频| 在线视频观看免费视频18| 精品免费人伦一区二区三区蜜桃| 亚洲日韩AV无码美腿丝袜| 国产精品一区二区熟女不卡| 精品无码av一区二区三区| 91精品福利一区二区| 国产一区二区精品久久呦| 婷婷色婷婷开心五月四| 久久国产热这里只有精品| 人妻无码人妻有码不卡| 中文字幕乱码亚洲三区| 精品亚洲成a人片在线观看| 人妻AV无码一区二区三区奥田咲|