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.展开更多
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.展开更多
This paper presents a system representation that can be applied to the description of the interaction between systems connected through common boundaries. The systems consist of partial differential equations that are...This paper presents a system representation that can be applied to the description of the interaction between systems connected through common boundaries. The systems consist of partial differential equations that are first order with respect to time, but spatially higher order. The representation is derived from the instantaneous multisymplectic Hamiltonian formalism;therefore, it possesses the physical consistency with respect to energy. In the interconnection, particular pairs of control inputs and observing outputs, called port variables, defined on the boundaries are used. The port variables are systematically introduced from the representation.展开更多
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.展开更多
A symbolic computation method to decide whether the solutions to the system Of linear partial differential equation is complete via using differential algebra and characteristic set is presented. This is a mechanizati...A symbolic computation method to decide whether the solutions to the system Of linear partial differential equation is complete via using differential algebra and characteristic set is presented. This is a mechanization method, and it can be carried out on the computer in the Maple environment.展开更多
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.展开更多
In this paper, Laplace decomposition method (LDM) and Pade approximant are employed to find approximate solutions for the Whitham-Broer-Kaup shallow water model, the coupled nonlinear reaction diffusion equations and ...In this paper, Laplace decomposition method (LDM) and Pade approximant are employed to find approximate solutions for the Whitham-Broer-Kaup shallow water model, the coupled nonlinear reaction diffusion equations and the system of Hirota-Satsuma coupled KdV. In addition, the results obtained from Laplace decomposition method (LDM) and Pade approximant are compared with corresponding exact analytical solutions.展开更多
Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.
In this paper, we introduce the notion of a (2+1)-dimenslonal differential equation describing three-dimensional hyperbolic spaces (3-h.s.). The (2+1)-dimensional coupled nonlinear Schrodinger equation and its...In this paper, we introduce the notion of a (2+1)-dimenslonal differential equation describing three-dimensional hyperbolic spaces (3-h.s.). The (2+1)-dimensional coupled nonlinear Schrodinger equation and its sister equation, the (2+1)-dimensional coupled derivative nonlinear Schrodinger equation, are shown to describe 3-h.s, The (2 + 1 )-dimensional generalized HF model:St=(1/2i[S,Sy]+2iσS)x,σx=-1/4i tr(SSxSy), in which S ∈ GLc(2)/GLc(1)×GLc(1),provides another example of (2+1)-dimensional differential equations describing 3-h.s. As a direct con-sequence, the geometric construction of an infinire number of conservation lairs of such equations is illustrated. Furthermore we display a new infinite number of conservation lairs of the (2+1)-dimensional nonlinear Schrodinger equation and the (2+1)-dimensional derivative nonlinear Schrodinger equation by a geometric way.展开更多
We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide...We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide a novel general structure of infinite-dimensional port-Hamiltonian systems (pHs) for which the Partitioned Finite Element Method (PFEM) straightforwardly applies. The proposed strategy is applied to abstract multidimensional linear hyperbolic and parabolic systems of PDEs. Then we show that instructional model problems based on the wave equation, Mindlin equation and heat equation fit within this unified framework. Secondly, we introduce the ongoing project SCRIMP (Simulation and Control of Interactions in Multi-Physics) developed for the numerical simulation of infinite-dimensional pHs. SCRIMP notably relies on the FEniCS open-source computing platform for the finite element spatial discretization. Finally, we illustrate how to solve the considered model problems within this framework by carefully explaining the methodology. As additional support, companion interactive Jupyter notebooks are available.展开更多
Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of ...Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.展开更多
In this paper we investigate the properties of the solutions of a class of second order neutral differential inequalities with distributed type deviating arguments . We apply the properties on the inequalities to obta...In this paper we investigate the properties of the solutions of a class of second order neutral differential inequalities with distributed type deviating arguments . We apply the properties on the inequalities to obtain some new criteria for all solutions of certain neutral hyperbolic partial functions differential equations to be oscillatory.展开更多
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 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.展开更多
In this paper, we study the oscillation of solutions to the systems of impulsive neutral delay parabolic partial differential equations. Under two different boundary conditions, we obtain some sufficient conditions fo...In this paper, we study the oscillation of solutions to the systems of impulsive neutral delay parabolic partial differential equations. Under two different boundary conditions, we obtain some sufficient conditions for oscillation of all solutions to the systems.展开更多
In this paper, some sufficient conditions for the oscillation for solutions to systems of n-th order partial functional differential equations are obtained.
区域综合能源系统(regional integrated energy system,RIES)的最优能流计算是求解RIES的设备配置、优化调度、故障分析等问题的基础。考虑供冷/热和供气管道传输能量的动态特性,建立RIES动态最优能流计算模型,其中基于特征线法获得了供...区域综合能源系统(regional integrated energy system,RIES)的最优能流计算是求解RIES的设备配置、优化调度、故障分析等问题的基础。考虑供冷/热和供气管道传输能量的动态特性,建立RIES动态最优能流计算模型,其中基于特征线法获得了供冷/热管道和供气管道动态偏微分方程的代数解析解。针对基于供冷/热系统质–量调节模式下管道能量传输时滞变量造成RIES的动态能流计算模型难以求解的问题,提出采用分段插值法获得供冷/热管道两端节点温度之间关系的近似表达式并加入动态最优能流计算模型中。此外,针对优化模型中供冷/热系统的流量与温度相乘的双线性项,提出一种能够缩紧松弛间隙的分段凸包络松弛方法将原混合整数非线性优化模型转化为混合整数二次约束规划模型,能够在保证计算精度的同时实现高效求解。最后以某个RIES算例进行分析,验证了所提方法的计算准确性和高效性,并与常用的质调节模式相比,表明在供冷/热系统质–量调节模式下能找到经济性更优的RIES运行点。展开更多
基金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.
文摘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.
文摘This paper presents a system representation that can be applied to the description of the interaction between systems connected through common boundaries. The systems consist of partial differential equations that are first order with respect to time, but spatially higher order. The representation is derived from the instantaneous multisymplectic Hamiltonian formalism;therefore, it possesses the physical consistency with respect to energy. In the interconnection, particular pairs of control inputs and observing outputs, called port variables, defined on the boundaries are used. The port variables are systematically introduced from the representation.
基金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.
文摘A symbolic computation method to decide whether the solutions to the system Of linear partial differential equation is complete via using differential algebra and characteristic set is presented. This is a mechanization method, and it can be carried out on the computer in the Maple environment.
基金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.
文摘In this paper, Laplace decomposition method (LDM) and Pade approximant are employed to find approximate solutions for the Whitham-Broer-Kaup shallow water model, the coupled nonlinear reaction diffusion equations and the system of Hirota-Satsuma coupled KdV. In addition, the results obtained from Laplace decomposition method (LDM) and Pade approximant are compared with corresponding exact analytical solutions.
基金This work is supported in part by NNSF of China(10571126)and in part by Program for New Century Excellent Talents in University.
文摘Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.
基金The project partially supported by National Natural Science Foundation of China
文摘In this paper, we introduce the notion of a (2+1)-dimenslonal differential equation describing three-dimensional hyperbolic spaces (3-h.s.). The (2+1)-dimensional coupled nonlinear Schrodinger equation and its sister equation, the (2+1)-dimensional coupled derivative nonlinear Schrodinger equation, are shown to describe 3-h.s, The (2 + 1 )-dimensional generalized HF model:St=(1/2i[S,Sy]+2iσS)x,σx=-1/4i tr(SSxSy), in which S ∈ GLc(2)/GLc(1)×GLc(1),provides another example of (2+1)-dimensional differential equations describing 3-h.s. As a direct con-sequence, the geometric construction of an infinire number of conservation lairs of such equations is illustrated. Furthermore we display a new infinite number of conservation lairs of the (2+1)-dimensional nonlinear Schrodinger equation and the (2+1)-dimensional derivative nonlinear Schrodinger equation by a geometric way.
文摘We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide a novel general structure of infinite-dimensional port-Hamiltonian systems (pHs) for which the Partitioned Finite Element Method (PFEM) straightforwardly applies. The proposed strategy is applied to abstract multidimensional linear hyperbolic and parabolic systems of PDEs. Then we show that instructional model problems based on the wave equation, Mindlin equation and heat equation fit within this unified framework. Secondly, we introduce the ongoing project SCRIMP (Simulation and Control of Interactions in Multi-Physics) developed for the numerical simulation of infinite-dimensional pHs. SCRIMP notably relies on the FEniCS open-source computing platform for the finite element spatial discretization. Finally, we illustrate how to solve the considered model problems within this framework by carefully explaining the methodology. As additional support, companion interactive Jupyter notebooks are available.
文摘Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.
文摘In this paper we investigate the properties of the solutions of a class of second order neutral differential inequalities with distributed type deviating arguments . We apply the properties on the inequalities to obtain some new criteria for all solutions of certain neutral hyperbolic partial functions differential equations to be oscillatory.
文摘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.
基金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(10471086).
文摘In this paper, we study the oscillation of solutions to the systems of impulsive neutral delay parabolic partial differential equations. Under two different boundary conditions, we obtain some sufficient conditions for oscillation of all solutions to the systems.
文摘In this paper, some sufficient conditions for the oscillation for solutions to systems of n-th order partial functional differential equations are obtained.
文摘区域综合能源系统(regional integrated energy system,RIES)的最优能流计算是求解RIES的设备配置、优化调度、故障分析等问题的基础。考虑供冷/热和供气管道传输能量的动态特性,建立RIES动态最优能流计算模型,其中基于特征线法获得了供冷/热管道和供气管道动态偏微分方程的代数解析解。针对基于供冷/热系统质–量调节模式下管道能量传输时滞变量造成RIES的动态能流计算模型难以求解的问题,提出采用分段插值法获得供冷/热管道两端节点温度之间关系的近似表达式并加入动态最优能流计算模型中。此外,针对优化模型中供冷/热系统的流量与温度相乘的双线性项,提出一种能够缩紧松弛间隙的分段凸包络松弛方法将原混合整数非线性优化模型转化为混合整数二次约束规划模型,能够在保证计算精度的同时实现高效求解。最后以某个RIES算例进行分析,验证了所提方法的计算准确性和高效性,并与常用的质调节模式相比,表明在供冷/热系统质–量调节模式下能找到经济性更优的RIES运行点。