Data-driven partial differential equation identification is a potential breakthrough to solve the lack of physical equations in complex dynamic systems.However,existing equation identification methods still cannot eff...Data-driven partial differential equation identification is a potential breakthrough to solve the lack of physical equations in complex dynamic systems.However,existing equation identification methods still cannot effectively identify equations from multivariable complex systems.In this work,we combine physical constraints such as dimension and direction of equation with data-driven method,and successfully identify the Navier-Stocks equations from the flow field data of Karman vortex street.This method provides an effective approach to identify partial differential equations of multivariable complex systems.展开更多
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.展开更多
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.展开更多
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 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.
A new theory on the construction of optimal truncated Low-Dimensional Dynamical Systems (LDDSs) with different physical meanings has been developed, The physical properties of the optimal bases are reflected in the us...A new theory on the construction of optimal truncated Low-Dimensional Dynamical Systems (LDDSs) with different physical meanings has been developed, The physical properties of the optimal bases are reflected in the user-defined optimal conditions, Through the analysis of linear and nonlinear examples, it is shown that the LDDSs constructed by using the Proper Orthogonal Decomposition (POD) method are not the optimum. After comparing the errors of LDDSs based on the new theory POD and Fourier methods, it is concluded that the LDDSs based on the new theory are optimally truncated and catch the desired physical properties of the systems.展开更多
This paper presents an algorithm to reduce a nonlinear algebraic partial differential equation system into the involutive characteristic set with respect to an abstract involutive prolongation direction, which covers ...This paper presents an algorithm to reduce a nonlinear algebraic partial differential equation system into the involutive characteristic set with respect to an abstract involutive prolongation direction, which covers the existing algorithms based on Riquier method, Thomas method, and Pommaret method. It also provides new algorithms for computing involutive characteristic sets due to the existence of new involutive directions. Experiments show that these new algorithms may be used to significantly reduce the computational steps in Wu-Ritt's characteristic set method for algebraic partial differential equations.展开更多
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.展开更多
The present work is based on the third-order partial differential equation (PDE) of acoustics of viscoelastic solids for the quasi-equilibrium (QE) component of the average normal stress. This PDE includes the stress-...The present work is based on the third-order partial differential equation (PDE) of acoustics of viscoelastic solids for the quasi-equilibrium (QE) component of the average normal stress. This PDE includes the stress-relaxation time (SRT) for the material and is applicable at any value of the SRT. The notion of a smart deicing system (SDS) for blade shells (BSs) of a wind turbine is specified. The work considers the stress in a BS as the one caused by the operational load on the BS. The work develops key design issues of a prospective ice-detection system (IDS) able to supply an array of the heating elements of an SDS with the element-individual spatiotemporal data and procedures for identification of the material parameters of atmospheric-ice (AI) layer accreted on the outer surfaces of the BSs. Both the SDS and IDS flexibly allow for complex, curvilinear and space-time-varying shapes of BSs. The proposed IDS presumes monitoring of the QE components of the normal stresses in BSs. The IDS is supposed to include an array of pressure-sensing resistors, also known as force-sensing resistors (FSRs), and communication hardware, as well as the parameter-identification software package (PISP), which provides the identification on the basis of the aforementioned PDE and the data measured by the FSRs. The IDS does not have hardware components located outside the outer surfaces of, or implanted in, BSs. The FSR array and communication hardware are reliable, and both cost- and energy-efficient. The present work extends methods of structural-health/operational-load monitoring (SH/OL-M) with measurements of the operational-load-caused stress in closed solid shells and, if the prospective PISP is used, endows the methods with identification of material parameters of the shells. The identification algorithms that can underlie the PISP are computationally efficient and suitable for implementation in the real-time mode. The identification model and algorithms can deal with not only the single-layer systems such as the BS layer without the AI layer or two-layer systems but also multi-layer systems. The outcomes can be applied to not only BSs of wind turbines but also non-QE closed single- or multi-layer deformable solid shells of various engineering systems (e.g., the shells of driver or passenger compartments of ships, cars, busses, airplanes, and other vehicles). The proposed monitoring of the normal-stress QE component in the mentioned shells extends the methods of SH/OL-M. The topic for the nearest research is a better adjustment of the settings for the FSR-based measurement of the mentioned components and a calibration of the parameter-identification model and algorithms, as well as the resulting improvement of the PISP.展开更多
In this paper, oscillation of solutions to a class of impulsive delay parabolic partial differential equations system with higher order Laplace operator is studied. Under two different boundary value conditions, we es...In this paper, oscillation of solutions to a class of impulsive delay parabolic partial differential equations system with higher order Laplace operator is studied. Under two different boundary value conditions, we establish some sufficient criteria with respect to the oscillations of such systems, employing first-order impulsive delay differential inequalities. The results fully reflect the influence action of impulsive and delay in oscillation.展开更多
A conceptual model for microscopic-macroscopic slow-fast stochastic systems is considered. A dynamical reduction procedure is presented in order to extract effective dynamics for this kind of systems. Under appropriat...A conceptual model for microscopic-macroscopic slow-fast stochastic systems is considered. A dynamical reduction procedure is presented in order to extract effective dynamics for this kind of systems. Under appropriate assumptions, the effective system is shown to approximate the original system, in the sense of a probabilistic convergence.展开更多
The optimal control problem for a nonlinear elliptic population system is considered.First, under certain hypotheses, the existence and uniqueness of coexistence state solutions are shown. Then the existence of the op...The optimal control problem for a nonlinear elliptic population system is considered.First, under certain hypotheses, the existence and uniqueness of coexistence state solutions are shown. Then the existence of the optimal control is given and the optimality system is established.展开更多
Trial equation method is a powerful tool for obtaining exact solutions of nonlinear differential equations. In this paper, the improved Boussinesq is reduced to an ordinary differential equation under the travelling w...Trial equation method is a powerful tool for obtaining exact solutions of nonlinear differential equations. In this paper, the improved Boussinesq is reduced to an ordinary differential equation under the travelling wave transformation. Trial equation method and the theory of complete discrimination system for polynomial are used to establish exact solutions of the improved Boussinesq equation.展开更多
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.展开更多
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.
The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference ...The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference scheme is constructed. The stability and convergence of that scheme are studied. Numerical experiments are carried out. The appropriate graphical illustrations and tables are given.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China(No.92152301).
文摘Data-driven partial differential equation identification is a potential breakthrough to solve the lack of physical equations in complex dynamic systems.However,existing equation identification methods still cannot effectively identify equations from multivariable complex systems.In this work,we combine physical constraints such as dimension and direction of equation with data-driven method,and successfully identify the Navier-Stocks equations from the flow field data of Karman vortex street.This method provides an effective approach to identify partial differential equations of multivariable complex systems.
文摘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.
文摘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.
文摘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.
基金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.
基金The project supported by the National Natural Science Foundation of ChinaLNM,Institute of Mechanics,CAS
文摘A new theory on the construction of optimal truncated Low-Dimensional Dynamical Systems (LDDSs) with different physical meanings has been developed, The physical properties of the optimal bases are reflected in the user-defined optimal conditions, Through the analysis of linear and nonlinear examples, it is shown that the LDDSs constructed by using the Proper Orthogonal Decomposition (POD) method are not the optimum. After comparing the errors of LDDSs based on the new theory POD and Fourier methods, it is concluded that the LDDSs based on the new theory are optimally truncated and catch the desired physical properties of the systems.
基金This work was supported in part by a National Key Basic Research Project of China(No.G19980306)the National Natural Science Foundation of China(Grant No.69725002).
文摘This paper presents an algorithm to reduce a nonlinear algebraic partial differential equation system into the involutive characteristic set with respect to an abstract involutive prolongation direction, which covers the existing algorithms based on Riquier method, Thomas method, and Pommaret method. It also provides new algorithms for computing involutive characteristic sets due to the existence of new involutive directions. Experiments show that these new algorithms may be used to significantly reduce the computational steps in Wu-Ritt's characteristic set method for algebraic partial differential equations.
文摘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.
文摘The present work is based on the third-order partial differential equation (PDE) of acoustics of viscoelastic solids for the quasi-equilibrium (QE) component of the average normal stress. This PDE includes the stress-relaxation time (SRT) for the material and is applicable at any value of the SRT. The notion of a smart deicing system (SDS) for blade shells (BSs) of a wind turbine is specified. The work considers the stress in a BS as the one caused by the operational load on the BS. The work develops key design issues of a prospective ice-detection system (IDS) able to supply an array of the heating elements of an SDS with the element-individual spatiotemporal data and procedures for identification of the material parameters of atmospheric-ice (AI) layer accreted on the outer surfaces of the BSs. Both the SDS and IDS flexibly allow for complex, curvilinear and space-time-varying shapes of BSs. The proposed IDS presumes monitoring of the QE components of the normal stresses in BSs. The IDS is supposed to include an array of pressure-sensing resistors, also known as force-sensing resistors (FSRs), and communication hardware, as well as the parameter-identification software package (PISP), which provides the identification on the basis of the aforementioned PDE and the data measured by the FSRs. The IDS does not have hardware components located outside the outer surfaces of, or implanted in, BSs. The FSR array and communication hardware are reliable, and both cost- and energy-efficient. The present work extends methods of structural-health/operational-load monitoring (SH/OL-M) with measurements of the operational-load-caused stress in closed solid shells and, if the prospective PISP is used, endows the methods with identification of material parameters of the shells. The identification algorithms that can underlie the PISP are computationally efficient and suitable for implementation in the real-time mode. The identification model and algorithms can deal with not only the single-layer systems such as the BS layer without the AI layer or two-layer systems but also multi-layer systems. The outcomes can be applied to not only BSs of wind turbines but also non-QE closed single- or multi-layer deformable solid shells of various engineering systems (e.g., the shells of driver or passenger compartments of ships, cars, busses, airplanes, and other vehicles). The proposed monitoring of the normal-stress QE component in the mentioned shells extends the methods of SH/OL-M. The topic for the nearest research is a better adjustment of the settings for the FSR-based measurement of the mentioned components and a calibration of the parameter-identification model and algorithms, as well as the resulting improvement of the PISP.
基金the Natural Science Foundation of Hunan Province under Grant 05JJ40008.
文摘In this paper, oscillation of solutions to a class of impulsive delay parabolic partial differential equations system with higher order Laplace operator is studied. Under two different boundary value conditions, we establish some sufficient criteria with respect to the oscillations of such systems, employing first-order impulsive delay differential inequalities. The results fully reflect the influence action of impulsive and delay in oscillation.
基金supported by NSF of China (10901065, 10971225, and11028102)the NSF Grants 1025422 and 0731201the Cheung Kong Scholars Program, and an open research grant from the State Key Laboratory for Nonlinear Mechanics at the Chinese Academy of Sciences
文摘A conceptual model for microscopic-macroscopic slow-fast stochastic systems is considered. A dynamical reduction procedure is presented in order to extract effective dynamics for this kind of systems. Under appropriate assumptions, the effective system is shown to approximate the original system, in the sense of a probabilistic convergence.
基金Supported by National Natural Science Foundation of P. R. China (60334040, 60225003)
文摘The optimal control problem for a nonlinear elliptic population system is considered.First, under certain hypotheses, the existence and uniqueness of coexistence state solutions are shown. Then the existence of the optimal control is given and the optimality system is established.
文摘Trial equation method is a powerful tool for obtaining exact solutions of nonlinear differential equations. In this paper, the improved Boussinesq is reduced to an ordinary differential equation under the travelling wave transformation. Trial equation method and the theory of complete discrimination system for polynomial are used to establish exact solutions of the improved Boussinesq equation.
文摘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.
文摘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.
文摘The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference scheme is constructed. The stability and convergence of that scheme are studied. Numerical experiments are carried out. The appropriate graphical illustrations and tables are given.
文摘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.