Misplaced Pages

Inertial manifold

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.

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 many physical applications, inertial manifolds express an interaction law between the small and large wavelength structures. Some say that the small wavelengths are enslaved by the large (e.g. synergetics). Inertial manifolds may also appear as slow manifolds common in meteorology, or as the center manifold in any bifurcation. Computationally, numerical schemes for partial differential equations seek to capture the long term dynamics and so such numerical schemes form an approximate inertial manifold.

Introductory Example

Consider the dynamical system in just two variables  p ( t ) {\displaystyle p(t)} and  q ( t ) {\displaystyle q(t)} and with parameter  a {\displaystyle a} :

d p d t = a p p q , d q d t = q + p 2 2 q 2 . {\displaystyle {\frac {dp}{dt}}=ap-pq\,,\qquad {\frac {dq}{dt}}=-q+p^{2}-2q^{2}.}
  • It possesses the one dimensional inertial manifold  M {\displaystyle {\mathcal {M}}} of q = p 2 / ( 1 + 2 a ) {\displaystyle q=p^{2}/(1+2a)} (a parabola).
  • This manifold is invariant under the dynamics because on the manifold  M {\displaystyle {\mathcal {M}}}
  d q d t = d d t p 2 1 + 2 a = 2 p d p d t 1 + 2 a = 2 a p 2 1 + 2 a 2 p 4 ( 1 + 2 a ) 2 {\displaystyle {\frac {dq}{dt}}={\frac {d}{dt}}{\frac {p^{2}}{1+2a}}={\frac {2p{\frac {dp}{dt}}}{1+2a}}={\frac {2ap^{2}}{1+2a}}-{\frac {2p^{4}}{(1+2a)^{2}}}} which is the same as
  q + p 2 2 q 2 = p 2 1 + 2 a + p 2 2 ( p 2 1 + 2 a ) 2 = 2 a p 2 1 + 2 a 2 p 4 ( 1 + 2 a ) 2 . {\displaystyle -q+p^{2}-2q^{2}=-{\frac {p^{2}}{1+2a}}+p^{2}-2\left({\frac {p^{2}}{1+2a}}\right)^{2}={\frac {2ap^{2}}{1+2a}}-{\frac {2p^{4}}{(1+2a)^{2}}}.}
  • The manifold  M {\displaystyle {\mathcal {M}}} attracts all trajectories in some finite domain around the origin because near the origin  d q d t q {\displaystyle {\frac {dq}{dt}}\approx -q} (although the strict definition below requires attraction from all initial conditions).

Hence the long term behavior of the original two dimensional dynamical system is given by the 'simpler' one dimensional dynamics on the inertial manifold  M {\displaystyle {\mathcal {M}}} , namely  d p d t = a p 1 1 + 2 a p 3 {\displaystyle {\frac {dp}{dt}}=ap-{\frac {1}{1+2a}}p^{3}} .

Definition

Let u ( t ) {\displaystyle u(t)} denote a solution of a dynamical system. The solution  u ( t ) {\displaystyle u(t)} may be an evolving vector in H = R n {\displaystyle H=\mathbb {R} ^{n}} or may be an evolving function in an infinite-dimensional Banach space  H {\displaystyle H} .

In many cases of interest the evolution of  u ( t ) {\displaystyle u(t)} is determined as the solution of a differential equation in  H {\displaystyle H} , say d u / d t = F ( u ( t ) ) {\displaystyle {du}/{dt}=F(u(t))} with initial value u ( 0 ) = u 0 {\displaystyle u(0)=u_{0}} . In any case, we assume the solution of the dynamical system can be written in terms of a semigroup operator, or state transition matrix, S : H H {\displaystyle S:H\to H} such that u ( t ) = S ( t ) u 0 {\displaystyle u(t)=S(t)u_{0}} for all times t 0 {\displaystyle t\geq 0} and all initial values  u 0 {\displaystyle u_{0}} . In some situations we might consider only discrete values of time as in the dynamics of a map.

An inertial manifold for a dynamical semigroup  S ( t ) {\displaystyle S(t)} is a smooth manifold  M {\displaystyle {\mathcal {M}}} such that

  1. M {\displaystyle {\mathcal {M}}} is of finite dimension,
  2. S ( t ) M M {\displaystyle S(t){\mathcal {M}}\subset {\mathcal {M}}} for all times  t 0 {\displaystyle t\geq 0} ,
  3. M {\displaystyle {\mathcal {M}}} attracts all solutions exponentially quickly, that is, for every initial value  u 0 H {\displaystyle u_{0}\in H} there exist constants  c j > 0 {\displaystyle c_{j}>0} such that dist ( S ( t ) u 0 , M ) c 1 e c 2 t {\displaystyle {\text{dist}}(S(t)u_{0},{\mathcal {M}})\leq c_{1}e^{-c_{2}t}} .

The restriction of the differential equation  d u / d t = F ( u ) {\displaystyle du/dt=F(u)} to the inertial manifold  M {\displaystyle {\mathcal {M}}} is therefore a well defined finite-dimensional system called the inertial system. Subtly, there is a difference between a manifold being attractive, and solutions on the manifold being attractive. Nonetheless, under appropriate conditions the inertial system possesses so-called asymptotic completeness: that is, every solution of the differential equation has a companion solution lying in  M {\displaystyle {\mathcal {M}}} and producing the same behavior for large time; in mathematics, for all  u 0 {\displaystyle u_{0}} there exists  v 0 M {\displaystyle v_{0}\in {\mathcal {M}}} and possibly a time shift  τ 0 {\displaystyle \tau \geq 0} such that dist ( S ( t ) u 0 , S ( t + τ ) v 0 ) 0 {\displaystyle {\text{dist}}(S(t)u_{0},S(t+\tau )v_{0})\to 0} as  t {\displaystyle t\to \infty } .

Researchers in the 2000s generalized such inertial manifolds to time dependent (nonautonomous) and/or stochastic dynamical systems (e.g.)

Existence

Existence results that have been proved address inertial manifolds that are expressible as a graph. The governing differential equation is rewritten more specifically in the form d u / d t + A u + f ( u ) = 0 {\displaystyle du/dt+Au+f(u)=0} for unbounded self-adjoint closed operator  A {\displaystyle A} with domain  D ( A ) H {\displaystyle D(A)\subset H} , and nonlinear operator  f : D ( A ) H {\displaystyle f:D(A)\to H} . Typically, elementary spectral theory gives an orthonormal basis of  H {\displaystyle H} consisting of eigenvectors  v j {\displaystyle v_{j}} : A v j = λ j v j {\displaystyle Av_{j}=\lambda _{j}v_{j}} , j = 1 , 2 , {\displaystyle j=1,2,\ldots } , for ordered eigenvalues 0 < λ 1 λ 2 {\displaystyle 0<\lambda _{1}\leq \lambda _{2}\leq \cdots } .

For some given number  m {\displaystyle m} of modes, P {\displaystyle P}  denotes the projection of  H {\displaystyle H} onto the space spanned by  v 1 , , v m {\displaystyle v_{1},\ldots ,v_{m}} , and Q = I P {\displaystyle Q=I-P}  denotes the orthogonal projection onto the space spanned by  v m + 1 , v m + 2 , {\displaystyle v_{m+1},v_{m+2},\ldots } . We look for an inertial manifold expressed as the graph  Φ : P H Q H {\displaystyle \Phi :PH\to QH} . For this graph to exist the most restrictive requirement is the spectral gap condition  λ m + 1 λ m c ( λ m + 1 + λ m ) {\displaystyle \lambda _{m+1}-\lambda _{m}\geq c({\sqrt {\lambda _{m+1}}}+{\sqrt {\lambda _{m}}})} where the constant  c {\displaystyle c} depends upon the system. This spectral gap condition requires that the spectrum of  A {\displaystyle A} must contain large gaps to be guaranteed of existence.

Approximate inertial manifolds

Several methods are proposed to construct approximations to inertial manifolds, including the so-called intrinsic low-dimensional manifolds.

The most popular way to approximate follows from the existence of a graph. Define the  m {\displaystyle m} slow variables  p ( t ) = P u ( t ) {\displaystyle p(t)=Pu(t)} , and the 'infinite' fast variables  q ( t ) = Q u ( t ) {\displaystyle q(t)=Qu(t)} . Then project the differential equation d u / d t + A u + f ( u ) = 0 {\displaystyle du/dt+Au+f(u)=0} onto both  P H {\displaystyle PH} and  Q H {\displaystyle QH} to obtain the coupled system d p / d t + A p + P f ( p + q ) = 0 {\displaystyle dp/dt+Ap+Pf(p+q)=0} and d q / d t + A q + Q f ( p + q ) = 0 {\displaystyle dq/dt+Aq+Qf(p+q)=0} .

For trajectories on the graph of an inertial manifold  M {\displaystyle M} , the fast variable  q ( t ) = Φ ( p ( t ) ) {\displaystyle q(t)=\Phi (p(t))} . Differentiating and using the coupled system form gives the differential equation for the graph:

d Φ d p [ A p + P f ( p + Φ ( p ) ) ] + A Φ ( p ) + Q f ( p + Φ ( p ) ) = 0. {\displaystyle -{\frac {d\Phi }{dp}}\left+A\Phi (p)+Qf(p+\Phi (p))=0.}

This differential equation is typically solved approximately in an asymptotic expansion in 'small'  p {\displaystyle p} to give an invariant manifold model, or a nonlinear Galerkin method, both of which use a global basis whereas the so-called holistic discretisation uses a local basis. Such approaches to approximation of inertial manifolds are very closely related to approximating center manifolds for which a web service exists to construct approximations for systems input by a user.

See also

References

  1. ^ R. Temam. Inertial manifolds. Mathematical Intelligencer, 12:68–74, 1990
  2. Roberts, A. J. (1985). "Simple examples of the derivation of amplitude equations for systems of equations possessing bifurcations". Journal of the Australian Mathematical Society, Series B. 27 (1). Cambridge University Press (CUP): 48–65. doi:10.1017/s0334270000004756. ISSN 0334-2700.
  3. Robinson, James C (1996-09-01). "The asymptotic completeness of inertial manifolds". Nonlinearity. 9 (5). IOP Publishing: 1325–1340. Bibcode:1996Nonli...9.1325R. doi:10.1088/0951-7715/9/5/013. ISSN 0951-7715. S2CID 250890338.
  4. Schmalfuss, Björn; Schneider, Klaus R. (2007-09-18). "Invariant Manifolds for Random Dynamical Systems with Slow and Fast Variables". Journal of Dynamics and Differential Equations. 20 (1). Springer Science and Business Media LLC: 133–164. Bibcode:2008JDDE...20..133S. doi:10.1007/s10884-007-9089-7. ISSN 1040-7294. S2CID 123477654.
  5. Pötzsche, Christian; Rasmussen, Martin (2009-02-18). "Computation of nonautonomous invariant and inertial manifolds" (PDF). Numerische Mathematik. 112 (3). Springer Science and Business Media LLC: 449–483. doi:10.1007/s00211-009-0215-9. ISSN 0029-599X. S2CID 6111461.
  6. Maas, U.; Pope, S.B. (1992). "Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space". Combustion and Flame. 88 (3–4). Elsevier BV: 239–264. doi:10.1016/0010-2180(92)90034-m. ISSN 0010-2180.
  7. Bykov, Viatcheslav; Goldfarb, Igor; Gol'dshtein, Vladimir; Maas, Ulrich (2006-06-01). "On a modified version of ILDM approach: asymptotic analysis based on integral manifolds". IMA Journal of Applied Mathematics. 71 (3). Oxford University Press (OUP): 359–382. doi:10.1093/imamat/hxh100. ISSN 1464-3634.
  8. Roberts, A. J. (1989). "The Utility of an Invariant Manifold Description of the Evolution of a Dynamical System". SIAM Journal on Mathematical Analysis. 20 (6). Society for Industrial & Applied Mathematics (SIAM): 1447–1458. doi:10.1137/0520094. ISSN 0036-1410.
  9. Foias, C.; Jolly, M.S.; Kevrekidis, I.G.; Sell, G.R.; Titi, E.S. (1988). "On the computation of inertial manifolds". Physics Letters A. 131 (7–8). Elsevier BV: 433–436. Bibcode:1988PhLA..131..433F. doi:10.1016/0375-9601(88)90295-2. ISSN 0375-9601.
  10. Roberts, A. J. (2002-06-04). "A holistic finite difference approach models linear dynamics consistently". Mathematics of Computation. 72 (241): 247–262. CiteSeerX 10.1.1.207.4820. doi:10.1090/S0025-5718-02-01448-5. S2CID 11525980.
  11. "Construct centre manifolds of ordinary or delay differential equations (autonomous)".
Category: