In the work of Colliander et al. (2020) a minimal lattice model was constructed describing the transfer of energy to high frequencies in the defocusing nonlinear Schrödinger equation. In the present ...work, we present a systematic study of the coherent structures, both standing and traveling, that arise in the context of this model. We find that the nonlinearly dispersive nature of the model is responsible for standing waves in the form of discrete compactons. On the other hand, analysis of the dynamical features of the simplest nontrivial variant of the model, namely the dimer case, yields both solutions where the intensity is trapped in a single site and solutions where the intensity moves between the two sites, which suggests the possibility of moving excitations in larger lattices. Such excitations are also suggested by the dynamical evolution associated with modulational instability. Our numerical computations confirm this expectation, and we systematically construct such traveling states as exact solutions in lattices of varying size, as well as explore their stability. A remarkable feature of these traveling lattice waves is that they are of “antidark” type, i.e., they are mounted on top of a non-vanishing background. These studies shed light on the existence, stability and dynamics of such standing and traveling states in 1+1 dimensions, and pave the way for exploration of corresponding configurations in higher dimensions.
•A minimal lattice model that describes the transfer of energy to high frequencies in the defocusing nonlinear Schrödinger equation supports both standing and traveling coherent structures.•Standing coherent structures take the form of compactons, which are solitary waves with compact support.•The dimer system can be solved exactly, and yields both solutions where the intensity is trapped in a single site and solutions where the intensity moves between the two lattice sites.•Lattice traveling waves, which take the form of a localized excitation superimposed on a nonzero background, can be constructed numerically for larger lattices.
We consider the Hamilton-Jacobi equationH(x,Du)+λ(x)u=c,x∈M, where M is a connected, closed and smooth Riemannian manifold. The functions H(x,p) and λ(x) are continuous. H(x,p) is convex, coercive ...with respect to p, and λ(x) changes the signs. The first breakthrough to this model was achieved by Jin-Yan-Zhao 11 under the Tonelli conditions. In this paper, we consider more detailed structure of the viscosity solution set and large time behavior of the viscosity solution on the Cauchy problem. To the best of our knowledge, it is the first detailed description of the large time behavior of the HJ equations with non-monotone dependence on the unknown function.
Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE ...discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.
This article introduces a novel systematic methodology for modeling a class of multidimensional linear mechanical systems that directly allows to obtain their infinite-dimensional port-Hamiltonian ...representation. While the approach is tailored to systems governed by specific kinematic assumptions, it encompasses a wide range of models found in current literature, including ℓ-dimensional elasticity models (where ℓ = 1, 2, 3), vibrating strings, torsion in circular bars, classical beam and plate models, among others. The methodology involves formulating the displacement field using primary generalized coordinates via a linear algebraic relation. The non-zero components of the strain tensor are then calculated and expressed using secondary generalized coordinates, enabling the characterization of the skew-adjoint differential operator associated with the port-Hamiltonian representation. By applying Hamilton's principle and employing a specially developed integration by parts formula for the considered class of differential operators, the port-Hamiltonian model is directly obtained, along with the definition of boundary inputs and outputs. To illustrate the methodology, the plate modeling process based on Reddy's third-order shear deformation theory is presented as an example. To the best of our knowledge, this is the first time that a port-Hamiltonian representation of this system is presented in the literature.
•Systematic modeling of multidimensional flexible linear mechanical systems using the port-Hamiltonian framework.•Explicit definition of energy and co-energy variables for this class of port Hamiltonian systems.•Explicit definition of distributed and boundary power conjugated ports suitable for control or interconnection purposes.
We propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of ...SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.
We show the Godbillon-Vey invariant arises as a 'restricted Casimir' invariant for three-dimensional ideal fluids associated to a foliation. We compare to a finite-dimensional system, the rattleback, ...where analogous phenomena occur.
Exponential/extended Runge–Kutta–Nyström (ERKN) methods have been validated by numerous theoretical and numerical results to be more efficient and suitable than classical Runge–Kutta–Nyström (RKN) ...methods in dealing with highly oscillatory Hamiltonian systems. Once numerical solutions are required to achieve a prescribed local error tolerance for practical problems, the embedded ERKN pairs with adaptive stepsize are needed. Given the higher efficiency of high-order methods than low-order methods, this paper focuses on high-order embedded ERKN pairs. Using the particular mapping from RKN methods into ERKN methods, we establish an approach to constructing high-order embedded ERKN pairs. By diagonalizing the frequency matrix, an implementation algorithm with less calculation amount than the direct calculation procedure is proposed for high-dimensional problems. In addition, the dispersion and dissipation of the ERKN pairs ERKN4(3) and ERKN8(6) are analyzed. Furthermore, the application of ERKN pairs to an important highly oscillatory problem, i.e., the wave equation prescribed with different boundary conditions is discussed. For the periodic boundary condition, a special implementation algorithm incorporated with FFT techniques is presented, which decreases the calculation amount in one time step from O(N2) to O(NlogN). Finally, numerical results with the FPU problem, the Klein–Gordon equation in the nonrelativistic limit regime, and a two-dimensional wave equation demonstrate the superiority of ERKN pairs over classical RKN pairs.
•An effective approach is presented to establish high-order embedded ERKN pairs.•A special implementation algorithm by diagonalizing the frequency matrix is proposed to reduce the calculation amount.•The application of the embedded ERKN pairs to wave equations prescribed with different boundary conditions is discussed.
The design of an observer-based state feedback (OBSF) controller with guaranteed passivity properties for port-Hamiltonian systems (PHS) is addressed using linear matrix inequalities (LMIs). The ...observer gain is freely chosen and the LMIs conditions such that the state feedback is equivalent to control by interconnection with an input strictly passive (ISP) and/or an output strictly passive (OSP) and zero state detectable (ZSD) port-Hamiltonian controller are established. It is shown that the proposed controller exponentially stabilizes a class of infinite-dimensional PHS and asymptotically stabilizes a class of finite-dimensional non-linear PHS. A Timoshenko beam model and a microelectromechanical system are used to illustrate the proposed approach.