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.展开更多
Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.
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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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 order to get the exact traveling wave solutions to nonlinear partial differential equation, the complete discrimination system for polynomial and direct integral method are applied to the considered equation. All s...In order to get the exact traveling wave solutions to nonlinear partial differential equation, the complete discrimination system for polynomial and direct integral method are applied to the considered equation. All single traveling wave solutions to the equation can be obtained. As an example, we give the solutions to (3 + 1)-dimensional breaking soliton 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.展开更多
The numerical analysis of the approximate inertial manifold in,weakly damped forced KdV equation is given. The results of numerical analysis under five models is the same as that of nonlinear spectral analysis.
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.展开更多
Under the travelling wave transformation, some nonlinear partial differential equations such as the Getmanou equation are transformed to ordinary differential equation. Then using trial equation method and combing com...Under the travelling wave transformation, some nonlinear partial differential equations such as the Getmanou equation are transformed to ordinary differential equation. Then using trial equation method and combing complete discrimination system for polynomial, the classifications of all single traveling wave solution to this equation are obtained.展开更多
In the article, the nonlinear equation is reduced to an ordinary differential equation under the travelling wave transformation. Using trial equation method, the ODE is reduced to the elementary integral form. In the ...In the article, the nonlinear equation is reduced to an ordinary differential equation under the travelling wave transformation. Using trial equation method, the ODE is reduced to the elementary integral form. In the end, complete discrimination system for polynomial is used to solve the corresponding integrals and obtain the classification of all single travelling wave solutions to the equation.展开更多
In order to get the traveling wave solutions of the Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZK-BBM) equation, it is reduced to an ordinary differential equation (ODE) under the travelling wave transformation first. T...In order to get the traveling wave solutions of the Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZK-BBM) equation, it is reduced to an ordinary differential equation (ODE) under the travelling wave transformation first. Then complete discrimination system for polynomial is applied to the ZK-BBM equation. The traveling wave solutions of the equation can be obtained.展开更多
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.展开更多
基金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 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.
文摘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.
文摘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.
文摘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.
基金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.
文摘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.
文摘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.
基金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 order to get the exact traveling wave solutions to nonlinear partial differential equation, the complete discrimination system for polynomial and direct integral method are applied to the considered equation. All single traveling wave solutions to the equation can be obtained. As an example, we give the solutions to (3 + 1)-dimensional breaking soliton 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.
文摘The numerical analysis of the approximate inertial manifold in,weakly damped forced KdV equation is given. The results of numerical analysis under five models is the same as that of nonlinear spectral analysis.
文摘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.
文摘Under the travelling wave transformation, some nonlinear partial differential equations such as the Getmanou equation are transformed to ordinary differential equation. Then using trial equation method and combing complete discrimination system for polynomial, the classifications of all single traveling wave solution to this equation are obtained.
文摘In the article, the nonlinear equation is reduced to an ordinary differential equation under the travelling wave transformation. Using trial equation method, the ODE is reduced to the elementary integral form. In the end, complete discrimination system for polynomial is used to solve the corresponding integrals and obtain the classification of all single travelling wave solutions to the equation.
文摘In order to get the traveling wave solutions of the Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZK-BBM) equation, it is reduced to an ordinary differential equation (ODE) under the travelling wave transformation first. Then complete discrimination system for polynomial is applied to the ZK-BBM equation. The traveling wave solutions of the equation can be obtained.
基金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.
