Spectral submanifold

Last updated
Schematic illustration of a spectral submanifold
W
(
E
)
{\displaystyle {\mathcal {W}}(E)}
emanating from a spectral subspace
E
{\displaystyle E}
. A trajectory
p
(
t
)
{\displaystyle p(t)}
in the reduced coordinates is mapped to the phase space via the manifold parametrization
W
(
p
)
{\displaystyle W(p)}
. Spectral submanifold.JPG
Schematic illustration of a spectral submanifold emanating from a spectral subspace . A trajectory in the reduced coordinates is mapped to the phase space via the manifold parametrization .

In dynamical systems, a spectral submanifold (SSM) is the unique smoothest invariant manifold serving as the nonlinear extension of a spectral subspace of a linear dynamical system under the addition of nonlinearities. [2] SSM theory provides conditions for when invariant properties of eigenspaces of a linear dynamical system can be extended to a nonlinear system, and therefore motivates the use of SSMs in nonlinear dimensionality reduction.

Contents

SSMs are chiefly employed for the exact model reduction of dynamical systems. For the automated computation of SSMs and the analysis of the reduced dynamics, open source online software packages such as SSMTool and SSMLearn have been published. These tools allow to study system dynamics either from the underlying equations of motion or from trajectory data, supporting both analytical and data-driven approaches. [3] [4] Detailed documentation for SSMTool is provided online. [5]

Definition

Consider a nonlinear ordinary differential equation of the form

with constant matrix and the nonlinearities contained in the smooth function .

Assume that for all eigenvalues of , that is, the origin is an asymptotically stable fixed point. Now select a span of eigenvectors of . Then, the eigenspace is an invariant subspace of the linearized system

Under addition of the nonlinearity to the linear system, generally perturbs into infinitely many invariant manifolds. Among these invariant manifolds, the unique smoothest one is referred to as the spectral submanifold.

An equivalent result for unstable SSMs holds for .

Existence

The spectral submanifold tangent to at the origin is guaranteed to exist provided that certain non-resonance conditions are satisfied by the eigenvalues in the spectrum of . [6] In particular, there can be no linear combination of equal to one of the eigenvalues of outside of the spectral subspace. If there is such an outer resonance, one can include the resonant mode into and extend the analysis to a higher-dimensional SSM pertaining to the extended spectral subspace.

Non-autonomous extension

The theory on spectral submanifolds extends to nonlinear non-autonomous systems of the form

with a quasiperiodic forcing term. [7]

Significance

Spectral submanifolds are useful for rigorous nonlinear dimensionality reduction in dynamical systems. The reduction of a high-dimensional phase space to a lower-dimensional manifold can lead to major simplifications by allowing for an accurate description of the system's main asymptotic behaviour. [8] For a known dynamical system, SSMs can be computed analytically by solving the invariance equations, and reduced models on SSMs may be employed for prediction of the response to forcing. [9]

Furthermore these manifolds may also be extracted directly from trajectory data of a dynamical system with the use of machine learning algorithms. [10]

See also

Related Research Articles

In mathematics, particularly linear algebra and functional analysis, a spectral theorem is a result about when a linear operator or matrix can be diagonalized. This is extremely useful because computations involving a diagonalizable matrix can often be reduced to much simpler computations involving the corresponding diagonal matrix. The concept of diagonalization is relatively straightforward for operators on finite-dimensional vector spaces but requires some modification for operators on infinite-dimensional spaces. In general, the spectral theorem identifies a class of linear operators that can be modeled by multiplication operators, which are as simple as one can hope to find. In more abstract language, the spectral theorem is a statement about commutative C*-algebras. See also spectral theory for a historical perspective.

<span class="mw-page-title-main">Affine space</span> Euclidean space without distance and angles

In mathematics, an affine space is a geometric structure that generalizes some of the properties of Euclidean spaces in such a way that these are independent of the concepts of distance and measure of angles, keeping only the properties related to parallelism and ratio of lengths for parallel line segments. Affine space is the setting for affine geometry.

<span class="mw-page-title-main">Nonlinear dimensionality reduction</span> Projection of data onto lower-dimensional manifolds

Nonlinear dimensionality reduction, also known as manifold learning, is any of various related techniques that aim to project high-dimensional data, potentially existing across non-linear manifolds which cannot be adequately captured by linear decomposition methods, onto lower-dimensional latent manifolds, with the goal of either visualizing the data in the low-dimensional space, or learning the mapping itself. The techniques described below can be understood as generalizations of linear decomposition methods used for dimensionality reduction, such as singular value decomposition and principal component analysis.

<span class="mw-page-title-main">Contact geometry</span> Branch of geometry

In mathematics, contact geometry is the study of a geometric structure on smooth manifolds given by a hyperplane distribution in the tangent bundle satisfying a condition called 'complete non-integrability'. Equivalently, such a distribution may be given as the kernel of a differential one-form, and the non-integrability condition translates into a maximal non-degeneracy condition on the form. These conditions are opposite to two equivalent conditions for 'complete integrability' of a hyperplane distribution, i.e. that it be tangent to a codimension one foliation on the manifold, whose equivalence is the content of the Frobenius theorem.

In mathematics, in particular in algebraic topology and differential geometry, the Stiefel–Whitney classes are a set of topological invariants of a real vector bundle that describe the obstructions to constructing everywhere independent sets of sections of the vector bundle. Stiefel–Whitney classes are indexed from 0 to n, where n is the rank of the vector bundle. If the Stiefel–Whitney class of index i is nonzero, then there cannot exist everywhere linearly independent sections of the vector bundle. A nonzero nth Stiefel–Whitney class indicates that every section of the bundle must vanish at some point. A nonzero first Stiefel–Whitney class indicates that the vector bundle is not orientable. For example, the first Stiefel–Whitney class of the Möbius strip, as a line bundle over the circle, is not zero, whereas the first Stiefel–Whitney class of the trivial line bundle over the circle, , is zero.

In differential geometry, a field in mathematics, a Poisson manifold is a smooth manifold endowed with a Poisson structure. The notion of Poisson manifold generalises that of symplectic manifold, which in turn generalises the phase space from Hamiltonian mechanics.

<span class="mw-page-title-main">Semisimple Lie algebra</span> Direct sum of simple Lie algebras

In mathematics, a Lie algebra is semisimple if it is a direct sum of simple Lie algebras.

<span class="mw-page-title-main">Bifurcation theory</span> Study of sudden qualitative behavior changes caused by small parameter changes

Bifurcation theory is the mathematical study of changes in the qualitative or topological structure of a given family of curves, such as the integral curves of a family of vector fields, and the solutions of a family of differential equations. Most commonly applied to the mathematical study of dynamical systems, a bifurcation occurs when a small smooth change made to the parameter values of a system causes a sudden 'qualitative' or topological change in its behavior. Bifurcations occur in both continuous systems and discrete systems.

In mathematics, the Lagrangian Grassmannian is the smooth manifold of Lagrangian subspaces of a real symplectic vector space V. Its dimension is 1/2n(n + 1) (where the dimension of V is 2n). It may be identified with the homogeneous space

<span class="mw-page-title-main">Hyperbolic equilibrium point</span>

In the study of dynamical systems, a hyperbolic equilibrium point or hyperbolic fixed point is a fixed point that does not have any center manifolds. Near a hyperbolic point the orbits of a two-dimensional, non-dissipative system resemble hyperbolas. This fails to hold in general. Strogatz notes that "hyperbolic is an unfortunate name—it sounds like it should mean 'saddle point'—but it has become standard." Several properties hold about a neighborhood of a hyperbolic point, notably

In the mathematics of evolving systems, the concept of a center manifold was originally developed to determine stability of degenerate equilibria. Subsequently, the concept of center manifolds was realised to be fundamental to mathematical modelling.

In the mathematical discipline of matrix theory, a Jordan matrix, named after Camille Jordan, is a block diagonal matrix over a ring R, where each block along the diagonal, called a Jordan block, has the following form:

Numerical continuation is a method of computing approximate solutions of a system of parameterized nonlinear equations,

A normally hyperbolic invariant manifold (NHIM) is a natural generalization of a hyperbolic fixed point and a hyperbolic set. The difference can be described heuristically as follows: For a manifold to be normally hyperbolic we are allowed to assume that the dynamics of itself is neutral compared with the dynamics nearby, which is not allowed for a hyperbolic set. NHIMs were introduced by Neil Fenichel in 1972. In this and subsequent papers, Fenichel proves that NHIMs possess stable and unstable manifolds and more importantly, NHIMs and their stable and unstable manifolds persist under small perturbations. Thus, in problems involving perturbation theory, invariant manifolds exist with certain hyperbolicity properties, which can in turn be used to obtain qualitative information about a dynamical system.

In mathematics, the slow manifold of an equilibrium point of a dynamical system occurs as the most common example of a center manifold. One of the main methods of simplifying dynamical systems, is to reduce the dimension of the system to that of the slow manifold—center manifold theory rigorously justifies the modelling. For example, some global and regional models of the atmosphere or oceans resolve the so-called quasi-geostrophic flow dynamics on the slow manifold of the atmosphere/oceanic dynamics, and is thus crucial to forecasting with a climate model.

In dynamical systems, a branch of mathematics, an invariant manifold is a topological manifold that is invariant under the action of the dynamical system. Examples include the slow manifold, center manifold, stable manifold, unstable manifold, subcenter manifold and inertial manifold.

In mathematics, inertial manifolds are concerned with the long term behavior of the solutions of dissipative dynamical systems. Inertial manifolds are finite-dimensional, smooth, invariant manifolds that contain the global attractor and attract all solutions exponentially quickly. Since an inertial manifold is finite-dimensional even if the original system is infinite-dimensional, and because most of the dynamics for the system takes place on the inertial manifold, studying the dynamics on an inertial manifold produces a considerable simplification in the study of the dynamics of the original system.

In geometric measure theory an approximate tangent space is a measure theoretic generalization of the concept of a tangent space for a differentiable manifold.

In the mathematics of dynamical systems, the concept of Lyapunov dimension was suggested by Kaplan and Yorke for estimating the Hausdorff dimension of attractors. Further the concept has been developed and rigorously justified in a number of papers, and nowadays various different approaches to the definition of Lyapunov dimension are used. Remark that the attractors with noninteger Hausdorff dimension are called strange attractors. Since the direct numerical computation of the Hausdorff dimension of attractors is often a problem of high numerical complexity, estimations via the Lyapunov dimension became widely spread. The Lyapunov dimension was named after the Russian mathematician Aleksandr Lyapunov because of the close connection with the Lyapunov exponents.

In mathematics, specifically in spectral theory, an eigenvalue of a closed linear operator is called normal if the space admits a decomposition into a direct sum of a finite-dimensional generalized eigenspace and an invariant subspace where has a bounded inverse. The set of normal eigenvalues coincides with the discrete spectrum.

References

  1. Jain, Shobhit; Haller, George (2022). "How to compute invariant manifolds and their reduced dynamics in high-dimensional finite element models". Nonlinear Dynamics. 107 (2): 1417–1450. doi: 10.1007/s11071-021-06957-4 . hdl: 20.500.11850/519249 . S2CID   232269982.
  2. Haller, George; Ponsioen, Sten (2016). "Nonlinear normal modes and spectral submanifolds: Existence, uniqueness and use in model reduction". Nonlinear Dynamics. 86 (3): 1493–1534. arXiv: 1602.00560 . doi:10.1007/s11071-016-2974-z. S2CID   44074026.
  3. Jain, Shobhit; Li, Mingwu; Thurnher, Thomas; Haller, George (2024-10-09), SSMTool 2.6: Computation of invariant manifolds in high-dimensional mechanics problems, [object Object], doi:10.5281/zenodo.13909795 , retrieved 2024-11-12
  4. haller-group/SSMLearn, Haller Group, Nonlinear Dynamics, 2024-11-12, retrieved 2024-11-12
  5. "SSMTool". ssmtool.mae.sustech.edu.cn. Retrieved 2024-11-12.
  6. Cabré, P.; Fontich, E.; de la Llave, R. (2003). "The parametrization method for invariant manifolds I: manifolds associated to non-resonant spectral subspaces". Indiana Univ. Math. J. 52: 283–328. doi:10.1512/iumj.2003.52.2245. hdl: 2117/876 .
  7. Haro, A.; de la Llave, R. (2006). "A parameterisation method for the computation of invariant tori and their whiskers in quasiperiodic maps: Rigorous results". Differ. Equ. 228 (2): 530–579. Bibcode:2006JDE...228..530H. doi:10.1016/j.jde.2005.10.005.
  8. Rega, Giuseppe; Troger, Hans (2005). "Dimension Reduction of Dynamical Systems: Methods, Models, Applications". Nonlinear Dynamics. 41 (1–3): 1–15. doi:10.1007/s11071-005-2790-3. S2CID   14728580.
  9. Ponsioen, Sten; Pedergnana, Tiemo; Haller, George (2018). "Automated computation of autonomous spectral submanifolds for nonlinear modal analysis". Journal of Sound and Vibration. 420: 269–295. arXiv: 1709.00886 . Bibcode:2018JSV...420..269P. doi:10.1016/j.jsv.2018.01.048. S2CID   44186335.
  10. Cenedese, Mattia; Axås, Joar; Bäuerlein, Bastian; Avila, Kerstin; Haller, George (2022). "Data-driven modeling and prediction of non-linearizable dynamics via spectral submanifolds". Nature Communications. 13 (1): 872. arXiv: 2201.04976 . Bibcode:2022NatCo..13..872C. doi:10.1038/s41467-022-28518-y. PMC   8847615 . PMID   35169152.