In this paper we describe new B-spline Gaussian collocation software for solving two-dimensional parabolic partial differential equations(PDEs)defined over a rectangular region.The numerical solution is represented as...In this paper we describe new B-spline Gaussian collocation software for solving two-dimensional parabolic partial differential equations(PDEs)defined over a rectangular region.The numerical solution is represented as a bi-variate piecewise polynomial(using a tensor product B-spline basis)with time-dependent unknown coefficients.These coefficients are determined by imposing collocation conditions:the numerical solution is required to satisfy the PDE and boundary conditions at images of the Gauss points mapped onto certain subregions of the spatial domain.This leads to a large system of time-dependent differential algebraic equations(DAEs)which is solved using the DAE solver,DASPK.We provide numerical results in which we use the new software,called BACOL2D,to solve three test problems.展开更多
A novel canonical Euler splitting method is proposed for nonlinear compositestiff functional differential-algebraic equations, the stability and convergence of themethod is evidenced, theoretical results are further c...A novel canonical Euler splitting method is proposed for nonlinear compositestiff functional differential-algebraic equations, the stability and convergence of themethod is evidenced, theoretical results are further confirmed by some numerical experiments.Especially, the numerical method and its theories can be applied to specialcases, such as delay differential-algebraic equations and integral differential-algebraicequations.展开更多
文摘In this paper we describe new B-spline Gaussian collocation software for solving two-dimensional parabolic partial differential equations(PDEs)defined over a rectangular region.The numerical solution is represented as a bi-variate piecewise polynomial(using a tensor product B-spline basis)with time-dependent unknown coefficients.These coefficients are determined by imposing collocation conditions:the numerical solution is required to satisfy the PDE and boundary conditions at images of the Gauss points mapped onto certain subregions of the spatial domain.This leads to a large system of time-dependent differential algebraic equations(DAEs)which is solved using the DAE solver,DASPK.We provide numerical results in which we use the new software,called BACOL2D,to solve three test problems.
基金National Natural Science Foundation of China(Grant No.11971412)Key Project of Education Department of Hunan Province(Grant No.20A484)Project of Hunan National Center for Applied Mathematics(Grant No.2020ZYT003).
文摘A novel canonical Euler splitting method is proposed for nonlinear compositestiff functional differential-algebraic equations, the stability and convergence of themethod is evidenced, theoretical results are further confirmed by some numerical experiments.Especially, the numerical method and its theories can be applied to specialcases, such as delay differential-algebraic equations and integral differential-algebraicequations.