ConceptComplete

Maximal Solutions and Continuation

The local existence theorem guarantees solutions on short intervals. The theory of maximal solutions addresses how far solutions can be extended and what happens at the boundary of their existence interval.


Maximal Interval of Existence

Definition8.5Maximal solution

A solution x:(tβˆ’,t+)β†’Rn\mathbf{x}: (t^-, t^+) \to \mathbb{R}^n of xβ€²=F(t,x)\mathbf{x}' = \mathbf{F}(t, \mathbf{x}) is maximal (or non-continuable) if it cannot be extended to a larger interval. The interval (tβˆ’,t+)(t^-, t^+) is called the maximal interval of existence. When t+<∞t^+ < \infty, we say the solution has a finite-time blowup or reaches the boundary of the domain Ξ©\Omega.

Theorem8.3Continuation theorem

Let F\mathbf{F} be continuous and locally Lipschitz on an open set Ξ©βŠ‚RΓ—Rn\Omega \subset \mathbb{R} \times \mathbb{R}^n. Then every solution of the IVP can be extended to a maximal solution. Moreover, if (tβˆ’,t+)(t^-, t^+) is the maximal interval and t+<∞t^+ < \infty, then for every compact set KβŠ‚Ξ©K \subset \Omega, there exists tK<t+t_K < t^+ such that (t,x(t))βˆ‰K(t, \mathbf{x}(t)) \notin K for all t∈(tK,t+)t \in (t_K, t^+).

In other words, as t→t+t \to t^+, the solution either leaves every compact subset of Ω\Omega (blows up or reaches the boundary).


Finite-Time Blowup

ExampleFinite-time blowup
  1. xβ€²=x2x' = x^2, x(0)=1x(0) = 1: Solution x(t)=1/(1βˆ’t)x(t) = 1/(1-t) blows up at t+=1t^+ = 1. Here x(t)β†’+∞x(t) \to +\infty.

  2. xβ€²=1+x2x' = 1 + x^2, x(0)=0x(0) = 0: Solution x(t)=tan⁑tx(t) = \tan t blows up at t+=Ο€/2t^+ = \pi/2.

  3. The system xβ€²=yx' = y, yβ€²=x+x2y' = x + x^2 with appropriate initial conditions can exhibit finite-time blowup in the norm βˆ₯x(t)βˆ₯β†’βˆž\|\mathbf{x}(t)\| \to \infty.

RemarkCriteria for global existence

If F\mathbf{F} satisfies a linear growth condition βˆ₯F(t,x)βˆ₯≀A+Bβˆ₯xβˆ₯\|\mathbf{F}(t, \mathbf{x})\| \leq A + B\|\mathbf{x}\| on RΓ—Rn\mathbb{R} \times \mathbb{R}^n, then by Gronwall's inequality, solutions satisfy:

βˆ₯x(t)βˆ₯≀(βˆ₯x0βˆ₯+A∣tβˆ’t0∣)eB∣tβˆ’t0∣\|\mathbf{x}(t)\| \leq (\|\mathbf{x}_0\| + A|t-t_0|)e^{B|t-t_0|}

which is finite for all finite tt. Therefore all solutions exist for all time. This is a sufficient condition for global existence.


Gronwall's Inequality

Theorem8.4Gronwall's inequality

Let u,Ξ±,Ξ²:[t0,T]β†’Ru, \alpha, \beta: [t_0, T] \to \mathbb{R} be continuous with Ξ²β‰₯0\beta \geq 0. If

u(t)≀α(t)+∫t0tΞ²(s)u(s) dsforΒ allΒ t∈[t0,T],u(t) \leq \alpha(t) + \int_{t_0}^{t} \beta(s)u(s)\,ds \quad \text{for all } t \in [t_0, T],

then

u(t)≀α(t)+∫t0tΞ±(s)Ξ²(s)exp⁑ ⁣(∫stΞ²(r) dr)ds.u(t) \leq \alpha(t) + \int_{t_0}^{t} \alpha(s)\beta(s)\exp\!\left(\int_s^t \beta(r)\,dr\right)ds.

In particular, if Ξ±\alpha is non-decreasing, then u(t)≀α(t)exp⁑ ⁣(∫t0tΞ²(s) ds)u(t) \leq \alpha(t)\exp\!\left(\int_{t_0}^t \beta(s)\,ds\right).

ExampleContinuous dependence via Gronwall

Let x(t)\mathbf{x}(t) and y(t)\mathbf{y}(t) solve xβ€²=F(t,x)\mathbf{x}' = \mathbf{F}(t,\mathbf{x}) with initial conditions x0\mathbf{x}_0 and y0\mathbf{y}_0 respectively. If F\mathbf{F} is Lipschitz with constant LL, then:

βˆ₯x(t)βˆ’y(t)βˆ₯≀βˆ₯x0βˆ’y0βˆ₯+∫t0tLβˆ₯x(s)βˆ’y(s)βˆ₯ ds.\|\mathbf{x}(t) - \mathbf{y}(t)\| \leq \|\mathbf{x}_0 - \mathbf{y}_0\| + \int_{t_0}^{t} L\|\mathbf{x}(s) - \mathbf{y}(s)\|\,ds.

By Gronwall: βˆ₯x(t)βˆ’y(t)βˆ₯≀βˆ₯x0βˆ’y0βˆ₯eL(tβˆ’t0)\|\mathbf{x}(t) - \mathbf{y}(t)\| \leq \|\mathbf{x}_0 - \mathbf{y}_0\| e^{L(t-t_0)}. This proves continuous dependence on initial conditions: nearby initial data produce nearby solutions on bounded time intervals.

RemarkSensitivity to initial conditions

The exponential factor eL(tβˆ’t0)e^{L(t-t_0)} in the continuous dependence estimate shows that errors can grow exponentially. This is the mathematical basis of sensitive dependence on initial conditions in chaotic systems. The bound is sharp: for the linear system xβ€²=Lxx' = Lx, the ratio ∣x(t;x0)βˆ’x(t;y0)∣/∣x0βˆ’y0∣=eLt|x(t; x_0) - x(t; y_0)| / |x_0 - y_0| = e^{Lt} exactly.