Numerical computation plays an important role in the study of differential equations with time-delay,because a simple and explicit analytic solution is usually un available.Time-stepping methods based on discretizing ...Numerical computation plays an important role in the study of differential equations with time-delay,because a simple and explicit analytic solution is usually un available.Time-stepping methods based on discretizing the temporal derivative with some step-size∆t are the main tools for this task.To get accurate numerical solutions,in many cases it is necessary to require∆t<τand this will be a rather unwelcome restriction whenτ,the quantity of time-delay,is small.In this paper,we propose a method for a class of time-delay problems,which is completely meshless.The idea lies in representing the solution by its Laplace inverse transform along a carefully de-signed contour in the complex plane and then approximating the contour integral by the Filon-Clenshaw-Curtis(FCC)quadrature in a few fast growing subintervals.The computations of the solution for all time points of interest are naturally parallelizable and for each time point the implementations of the FCC quadrature in all subintervals are also parallelizable.For each time point and each subinterval,the FCC quadrature can be implemented by fast Fourier transform.Numerical results are given to check the efficiency of the proposed method.展开更多
基金The first author was supported by NSFC(No.11771313,No.61573010)the Project of China Postdoctoral Science Foundation(No.2015M580777,No.2016T90841)The second author was supported by NSFC(No.11771163).
文摘Numerical computation plays an important role in the study of differential equations with time-delay,because a simple and explicit analytic solution is usually un available.Time-stepping methods based on discretizing the temporal derivative with some step-size∆t are the main tools for this task.To get accurate numerical solutions,in many cases it is necessary to require∆t<τand this will be a rather unwelcome restriction whenτ,the quantity of time-delay,is small.In this paper,we propose a method for a class of time-delay problems,which is completely meshless.The idea lies in representing the solution by its Laplace inverse transform along a carefully de-signed contour in the complex plane and then approximating the contour integral by the Filon-Clenshaw-Curtis(FCC)quadrature in a few fast growing subintervals.The computations of the solution for all time points of interest are naturally parallelizable and for each time point the implementations of the FCC quadrature in all subintervals are also parallelizable.For each time point and each subinterval,the FCC quadrature can be implemented by fast Fourier transform.Numerical results are given to check the efficiency of the proposed method.
