advanced30 min studyLesson 331 in Mathematics

Flows, orbits, and phase portraits

Treat an autonomous differential equation as motion through state space, distinguishing a geometric orbit from its clock and its existence interval.

On this page

What you’ll understand

  • Use uniqueness to derive the local flow composition law.
  • Read invariant regions and equilibrium barriers in a phase portrait.
  • Distinguish smooth local solvability from existence for all time.

A differential equation specifies a velocity at every allowed state. Instead of asking only for a formula, follow how those velocities organize the whole state space: which regions are preserved, where motion slows, and which states can be approached. The resulting phase portrait records possible motion without assigning every point a time stamp. To reason correctly about it, keep three objects separate: the vector field, the time-parametrized solution, and the geometric orbit traced by that solution.

A local flow is a family of compatible solutions

Let f be continuously differentiable on an open subset U of Euclidean space. The autonomous equation x dot equals f(x) has a unique local solution through every initial state. Write its value after elapsed time t as phi subscript t of x. Uniqueness says that evolving for s and then for t agrees with evolving for s plus t, whenever all the required solution segments exist. The identity at time zero and this composition law make solution curves into a local flow.

Autonomous means that the velocity depends only on state, not on the external clock. A periodically forced oscillator needs a forcing phase added to its state before it becomes autonomous; dropping that phase can make projected curves cross. Distinct geometric orbits of a unique autonomous flow cannot cross at a shared state. If two parametrized solutions reach that same state at different times, uniqueness identifies their future and past portions after shifting the clock, rather than creating two competing velocity choices.

ϕ0(x)=x,ϕt(ϕs(x))=ϕt+s(x)\phi_0(x)=x,\qquad \phi_t(\phi_s(x))=\phi_{t+s}(x)
The composition law holds on its natural domain: each indicated solution segment must exist and remain in the vector field's domain.↔ Scroll wide equations horizontally.

An orbit forgets speed, but preserves direction

The orbit through a point is the set of states visited over its maximal existence interval. Equilibria are singleton orbits where f vanishes. In one dimension, the sign of f supplies arrows between equilibria, and uniqueness prevents a solution from passing through an equilibrium in finite time and then leaving it. A positive scalar multiple of f has the same oriented geometric trajectories, while generally changing their speed. A time trace and a phase portrait therefore answer different questions about the same model.

A set is positively invariant if every solution starting there stays there for all future times for which the model is considered. Full invariance additionally includes backward evolution. A bounded interval whose boundary arrows point inward is a useful positive trapping region, but checking a picture at a handful of boundary points is not a proof for a higher-dimensional boundary. Long-time accumulation is described by an omega-limit set; compact containment supplies convergent subsequences, while an escaping orbit can have no finite accumulation point at all.

ω(x)={y:∃tj→∞ with ϕtj(x)→y}\omega(x)=\{y:\exists t_j\to\infty\ \text{with}\ \phi_{t_j}(x)\to y\}
This future limit-set definition presumes the solution exists for arbitrarily large positive times. A precompact forward orbit inside U has a nonempty compact omega-limit set.↔ Scroll wide equations horizontally.

Smoothness is local, completeness is additional

Even a polynomial vector field need not generate a flow for every real time. The scalar equation x dot equals x squared has solutions that blow up in finite forward time from positive initial states. Local existence and uniqueness are fully compatible with this failure of completeness. A standard continuation result says a maximal solution cannot stop at a finite time while remaining in a compact subset of U. It must instead escape every such compact subset, through unbounded growth or approach to the domain boundary.

This distinction matters in numerical work. A plot that stops at a finite axis limit may show a plotting choice, a solver failure, or actual mathematical blow-up. A trajectory that seems to settle near a stationary point may simply move very slowly. Check invariance, sign information, and existence estimates alongside the computation. Phase-space coordinates also matter: a second-order scalar equation requires both position and velocity as state, because position alone does not determine which direction the system will move next.

Read a logistic orbit and its clock

For x dot = 3x(1−x/4) with x(0)=1, find the future orbit, its limit, and the state at time (log 3)/3.

  1. The equilibria are zero and four. The vector field is positive between them, so uniqueness and the equilibrium barriers keep this solution in (0,4) and moving upward.

  2. Separating variables gives x(t)=4/(1+3 exp(−3t)). It has the prescribed initial value and differentiating verifies the differential equation. The formula exists for every real t and ranges through the entire open interval (0,4).

  3. At t=(log 3)/3, the exponential equals one third, giving x=2. As t tends to positive infinity the solution approaches four; as t tends to negative infinity it approaches zero without reaching either endpoint.

The result

The full geometric orbit is (0,4), while its forward portion from the given initial state is [1,4). Four is the future limit; the clock determines when the orbit passes through two.

Explore the model

Interactive lab

Interactive lab · fields

Separate a velocity field from a time trace

Pause and reset, then set both controls to source strength q = 0 and rotation Ω = 1. Predict arrows tangent to circles about the origin. For a fresh comparison, pause and reset, then set q = 0 and Ω = 2: the directions stay tangential while the field magnitudes double at a fixed point. If interpreted as an ODE, these two fields have the same circular orbits with different speeds. The instrument draws the softened vector field; its moving marks are illustrative and are not a numerical integration of arbitrary initial trajectories.

Separate a velocity field from a time traceThe softened illustration is q·r/(r² + 0.2) + Ω(−y, x). Moving markers indicate local direction, not integrated particle paths or a Maxwell solution. Current readings: Fx at (1, 0): 0.833; Fy at (1, 0): 0.5; Field magnitude there: 0.972. All parameters are adjustable with the labeled controls below.−3.2−1.601.63.2−1.5−0.7500.751.5xyq = 1
Local field directionSource strength
Fx at (1, 0)
0.833
Fy at (1, 0)
0.5
Field magnitude there
0.972
t = 0 s

The softened illustration is q·r/(r² + 0.2) + Ω(−y, x). Moving markers indicate local direction, not integrated particle paths or a Maxwell solution. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA phase portrait specifies geometric motion and direction; the vector-field magnitude supplies the clock.

Check your understanding

What does continuous differentiability of f on an open domain guarantee for x dot = f(x)?

Try it yourself

Solve x dot = x² with x(0)=2. Find the maximal real-time interval containing zero and explain why the local flow law remains valid despite finite-time blow-up.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Differentiate the reciprocal of a nonzero solution, then locate the denominator's zero.

Compare with the solution

The reciprocal satisfies (1/x) dot = −1, so x(t)=2/(1−2t). The denominator vanishes at t=1/2, and the maximal interval containing zero is (−infinity,1/2). The vector field is smooth everywhere, but the solution escapes every bounded set as that endpoint approaches. Composition still holds whenever both time segments and their sum lie within the corresponding solution domains. The failure is global completeness, not uniqueness or the local composition identity.

Take it with you

  • An autonomous flow composes elapsed-time evolutions wherever they exist.
  • An orbit is a geometric set, while a solution also records speed and time.
  • Compact containment supports continuation and limit-set conclusions that local smoothness alone cannot supply.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.