Astrophysical Stellar Convection (ODE) First Described: 1966

Moore-Spiegel Attractor

Discovered by D.W. Moore & Edward A. Spiegel

Executive Summary

Pulsating stellar envelopes and chaotic convection waves occurring in variable giant stars.

System of Equations

Phase Space Formulation
dxdt=y\frac{dx}{dt} = y
dydt=z\frac{dy}{dt} = z
dzdt=z(TR+Rx2)yTx\frac{dz}{dt} = -z - (T - R + R x^2) y - T x

Governing Parameters

TT
Thermal Relaxation Parameter = 26.0

Ratio of cooling time to dynamical expansion time in the pulsating stellar gas shell.

RR
Reynolds / Buoyancy Ratio = 100.0

Measures radiative buoyancy forces driving convective instabilities against gravity.

Overview

Proposed in 1966 by astrophysicists Derek W. Moore and Edward A. Spiegel, this third-order differential equation system models the chaotic vertical oscillations of a compressible gas parcel in a variable star’s atmosphere.

In stars undergoing thermodynamic pulsations (such as Cepheid variables or red supergiants), buoyancy forces can destabilize the stellar envelope. As gas rises, it cools and radiates heat away; as it plunges back inward, it heats up and re-expands. Moore and Spiegel proved that non-linear coupling between thermal dissipation and acoustic resonance drives the stellar surface into turbulent, chaotic pulsations.

The Mathematical Model

Represented as a third-order autonomous ODE x...+x¨+(TR+Rx2)x˙+Tx=0\dddot{x} + \ddot{x} + (T - R + Rx^2)\dot{x} + Tx = 0, rewritten in state-space coordinates (x,y,z)(x, y, z):

dxdt=y\frac{dx}{dt} = y

dydt=z\frac{dy}{dt} = z

dzdt=z(TR+Rx2)yTx\frac{dz}{dt} = -z - (T - R + Rx^2)y - Tx

Coordinates

  • xx: Vertical displacement of the stellar fluid element relative to equilibrium.
  • yy: Radial velocity of the oscillating stellar gas.
  • zz: Acceleration and rate of thermal buoyant energy transfer.

Phase Space Geometry

  • Funnel-and-Wave Morphology: The non-linear damping term (TR+Rx2)(T - R + Rx^2) acts like a Van der Pol oscillator where dissipation depends directly on the square of displacement x2x^2.
  • Explosive Outbursts: The trajectory winds through a tight spiral at small amplitudes before violently expanding outward into wide, sweeping parabolic arcs, mimicking flares and irregular stellar eruptions.
Previous Type
Rikitake Dynamo Attractor
All Attractor Types
Next Type
Polynomial Attractor (Sprott 3D System)