The object of this paper is to investigate the superconvergence properties of finite element approximations to parabolic and hyperbolic integro-differential equations. The quasi projection technique introduced earlier...The object of this paper is to investigate the superconvergence properties of finite element approximations to parabolic and hyperbolic integro-differential equations. The quasi projection technique introduced earlier by Douglas et al. is developed to derive the O(h2r) order knot superconvergence in the case of a single space variable, and to show the optimal order negative norm estimates in the case of several space variables.展开更多
The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the un...The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.展开更多
In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differen...In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differential equations with nonlinear diffusion coefficient are investigated.Sufficient conditions for oscillations of such equations are obtained.展开更多
We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the...We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the original method,i.e.,it is based on a wave decomposition at grid cell interfaces,it can be used to approximate hyperbolic problems in divergence form as well as in quasilinear form and limiting is introduced in the form of a wave limiter.展开更多
Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a...Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.展开更多
This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functio...This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
In this study, we prove the existence, uniqueness, and continuous dependence upon the data of solution to integro-differential hyperbolic equation with purely nonlocal （integral） conditions. The proofs are based on ...In this study, we prove the existence, uniqueness, and continuous dependence upon the data of solution to integro-differential hyperbolic equation with purely nonlocal （integral） conditions. The proofs are based on a priori estimates and Laplace transform method. Finally, we obtain the solution using a numerical technique （Stehfest algorithm） by inverting the Laplace transform.展开更多
We first define a kind of new function space, called the space of twice conormal distributions. With some estimates on these function spaces, we can prove that if there exist two characteristic hypersurfaces bearing s...We first define a kind of new function space, called the space of twice conormal distributions. With some estimates on these function spaces, we can prove that if there exist two characteristic hypersurfaces bearing strong and weak singularities respectively intersect transversally, then some new singularities will take place anti their strength will be no more than the original weak one.展开更多
This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a pie...This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.展开更多
This paper is concerned with the oscillations of neutral hyperbolic partial differential equations with delays. Necessary and sufficient, conditions are obtained for the oscillations of all solutions of the equations,...This paper is concerned with the oscillations of neutral hyperbolic partial differential equations with delays. Necessary and sufficient, conditions are obtained for the oscillations of all solutions of the equations, and these results are illustrated by some examples.展开更多
This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discret...This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.展开更多
During May 11–13,2022,the International Workshop on Efficient High-Order TimeDiscretization Methods for Partial Differential Equations(PDEs)took place at Villa Orlandiin Anacapri,Italy,a conference center of the Unive...During May 11–13,2022,the International Workshop on Efficient High-Order TimeDiscretization Methods for Partial Differential Equations(PDEs)took place at Villa Orlandiin Anacapri,Italy,a conference center of the University of Naples Federico II.Due to theCOVID-19 pandemic,the workshop was held in a hybrid format.Approximately 50 seniorresearchers,young scholars,and Ph.D.students attended this workshop.The purpose of theevent was to explore recent trends and directions in the area of time discretization for thenumerical solution of evolutionary PDEs with particular focus to high-order methods forhyperbolic systems with source terms and advection-diffusion-reaction equations,and withspecial emphasis on efficient time-stepping methods such as the implicit-explicit(IMEX),and the semi-implicit and strong stability preserving(SSP)time discretization.This focusedsection entitled“Efficient High-Order Time Discretization Methods for Partial DifferentialEquations”in Communications on Applied Mathematics and Computation(CAMC)consistsof six regularly peer-reviewed manuscripts,selected from submissions of works presentedduring the workshop.展开更多
Recently, Kreith, Kusano and Yoshida have studied the oscillation property of the hyperbolic equation u<sub>11</sub>-△u+c（t,x,u）=f（t,x）, （t,x）∈R<sub>+</sub>×Ωwith boundary conditi...Recently, Kreith, Kusano and Yoshida have studied the oscillation property of the hyperbolic equation u<sub>11</sub>-△u+c（t,x,u）=f（t,x）, （t,x）∈R<sub>+</sub>×Ωwith boundary condition （?）u/（?）n=g（t,x）, （t, x）∈R<sub>+</sub>×（?）Ω,and obtained some sufficient criterions for solution oscillation. In this note, we shall discuss the oscillation properties of solutions for a class of hyperbolic functional differential展开更多
For a family of linear hyperbolic damped stochastic wave equations with rapidly oscillating coefficients, we establish the homogenization result by using the sigma-convergence method. This is achieved under an abstrac...For a family of linear hyperbolic damped stochastic wave equations with rapidly oscillating coefficients, we establish the homogenization result by using the sigma-convergence method. This is achieved under an abstract assumption covering special cases like the periodicity, the almost periodicity and some others.展开更多
In this paper we develop a new closure theory for moment approximationsin kinetic gas theory and derive hyperbolic moment equations for 13 fluid variablesincluding stress and heat flux. Classical equations have either...In this paper we develop a new closure theory for moment approximationsin kinetic gas theory and derive hyperbolic moment equations for 13 fluid variablesincluding stress and heat flux. Classical equations have either restricted hyperbolicity regions like Grad’s moment equations or fail to include higher moments in apractical way like the entropy maximization approach. The new closure is based onPearson-Type-IV distributions which reduce to Maxwellians in equilibrium, but allowanisotropies and skewness in non-equilibrium. The closure relations are essentiallyexplicit and easy to evaluate. Hyperbolicity is shown numerically for a large range ofvalues. Numerical solutions of Riemann problems demonstrate the capability of thenew equations to handle strong non-equilibrium.展开更多
Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bou...Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bounded or unbounded domain in R~3 and obtain the error estimates of the approximate solution in energy norm and local maximum norm.展开更多
Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the p...Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the production line,the manufacturing layer and the workshop layer.The dynamics of re-entrant production lines are governed by hyperbolic partial differential equations(PDEs)based on the law of mass conservation.展开更多
基金The NNSF (99200204) of Liaoning Province, China.
文摘The object of this paper is to investigate the superconvergence properties of finite element approximations to parabolic and hyperbolic integro-differential equations. The quasi projection technique introduced earlier by Douglas et al. is developed to derive the O(h2r) order knot superconvergence in the case of a single space variable, and to show the optimal order negative norm estimates in the case of several space variables.
基金the National Natural Science Foundation of China(Nos.11671282 and 12171339)。
文摘The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.
基金Supported by the Natural Science Foundation of China(10471086)Supported by the Science Research Foundation of Department of Education of Hunan Province(07C164)
文摘In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differential equations with nonlinear diffusion coefficient are investigated.Sufficient conditions for oscillations of such equations are obtained.
基金This work was supported by the DFG through HE 4858/4-1
文摘We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the original method,i.e.,it is based on a wave decomposition at grid cell interfaces,it can be used to approximate hyperbolic problems in divergence form as well as in quasilinear form and limiting is introduced in the form of a wave limiter.
文摘Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.
文摘This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.
基金partially support by the University Putra Malaysia under the Grant Scheme having project number GP-IBT/2013/9420100
文摘In this study, we prove the existence, uniqueness, and continuous dependence upon the data of solution to integro-differential hyperbolic equation with purely nonlocal （integral） conditions. The proofs are based on a priori estimates and Laplace transform method. Finally, we obtain the solution using a numerical technique （Stehfest algorithm） by inverting the Laplace transform.
文摘We first define a kind of new function space, called the space of twice conormal distributions. With some estimates on these function spaces, we can prove that if there exist two characteristic hypersurfaces bearing strong and weak singularities respectively intersect transversally, then some new singularities will take place anti their strength will be no more than the original weak one.
基金partly supported by SRF for ROCS, SEMsupported by a grant from the "project 211 (phase Ⅲ)" of the Southwestern University of Finance and Economics
文摘This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.
基金Supported by Natural Science Foundation of Hebei Province（102160） and Natural Science of Education office in Hebei Province （2004123）,
文摘This paper is concerned with the oscillations of neutral hyperbolic partial differential equations with delays. Necessary and sufficient, conditions are obtained for the oscillations of all solutions of the equations, and these results are illustrated by some examples.
基金supported by the National Natural Science Foundation of China under Grant Nos.61821004 and 62250056the Natural Science Foundation of Shandong Province under Grant Nos.ZR2021ZD14 and ZR2021JQ24+1 种基金Science and Technology Project of Qingdao West Coast New Area under Grant Nos.2019-32,2020-20,2020-1-4,High-level Talent Team Project of Qingdao West Coast New Area under Grant No.RCTDJC-2019-05Key Research and Development Program of Shandong Province under Grant No.2020CXGC01208.
文摘This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.
文摘During May 11–13,2022,the International Workshop on Efficient High-Order TimeDiscretization Methods for Partial Differential Equations(PDEs)took place at Villa Orlandiin Anacapri,Italy,a conference center of the University of Naples Federico II.Due to theCOVID-19 pandemic,the workshop was held in a hybrid format.Approximately 50 seniorresearchers,young scholars,and Ph.D.students attended this workshop.The purpose of theevent was to explore recent trends and directions in the area of time discretization for thenumerical solution of evolutionary PDEs with particular focus to high-order methods forhyperbolic systems with source terms and advection-diffusion-reaction equations,and withspecial emphasis on efficient time-stepping methods such as the implicit-explicit(IMEX),and the semi-implicit and strong stability preserving(SSP)time discretization.This focusedsection entitled“Efficient High-Order Time Discretization Methods for Partial DifferentialEquations”in Communications on Applied Mathematics and Computation(CAMC)consistsof six regularly peer-reviewed manuscripts,selected from submissions of works presentedduring the workshop.
文摘Recently, Kreith, Kusano and Yoshida have studied the oscillation property of the hyperbolic equation u<sub>11</sub>-△u+c（t,x,u）=f（t,x）, （t,x）∈R<sub>+</sub>×Ωwith boundary condition （?）u/（?）n=g（t,x）, （t, x）∈R<sub>+</sub>×（?）Ω,and obtained some sufficient criterions for solution oscillation. In this note, we shall discuss the oscillation properties of solutions for a class of hyperbolic functional differential
基金the support of the CETIC(African Center of Excellence in Information and Communication Technologies)the support of the Humboldt Foundation
文摘For a family of linear hyperbolic damped stochastic wave equations with rapidly oscillating coefficients, we establish the homogenization result by using the sigma-convergence method. This is achieved under an abstract assumption covering special cases like the periodicity, the almost periodicity and some others.
文摘In this paper we develop a new closure theory for moment approximationsin kinetic gas theory and derive hyperbolic moment equations for 13 fluid variablesincluding stress and heat flux. Classical equations have either restricted hyperbolicity regions like Grad’s moment equations or fail to include higher moments in apractical way like the entropy maximization approach. The new closure is based onPearson-Type-IV distributions which reduce to Maxwellians in equilibrium, but allowanisotropies and skewness in non-equilibrium. The closure relations are essentiallyexplicit and easy to evaluate. Hyperbolicity is shown numerically for a large range ofvalues. Numerical solutions of Riemann problems demonstrate the capability of thenew equations to handle strong non-equilibrium.
基金China State Major Key Project for Basic Researches
文摘Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bounded or unbounded domain in R~3 and obtain the error estimates of the approximate solution in energy norm and local maximum norm.
文摘Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the production line,the manufacturing layer and the workshop layer.The dynamics of re-entrant production lines are governed by hyperbolic partial differential equations(PDEs)based on the law of mass conservation.
