Module

for

The Lorenz Attractor

Lorenz Attractor

The Lorenz attractor is a set of differential equations which are popular in the field of Chaos. The equations describe the flow of fluid in a box which is heated along the bottom. This model was intended to simulate medium-scale atmospheric convection.  Lorenz simplified some of the Navier-Stokes equations in the area of fluid dynamics and obtained three ordinary differential equations

,

.

The parameter p is the Prandtl number, is the quotient of the Rayleigh number  and critical Rayleigh number  and b is a geometric factor.  Lorenz is attributed to using the values .

There are three critical points (0,0,0)  corresponds to no convection,  and the two points

The latter two points are to be stable, only if the following equation holds

.

Computer Programs  Lorenz Attractor  Lorenz Attractor

Program (Runge-Kutta Method in 3D space)  To compute a numerical approximation for the solution of the initial value problem

with  ,

with  ,

with  ,

over the interval
at a discrete set of points.

Mathematica Subroutine (Runge-Kutta Method in n-dimensions).

``````

``````

Example 1.  Solve the Lorenz I. V. P.
and
and

and
Solution 1.

Example 2.  Solve the Lorenz I. V. P.
and
and

and
Use Mathematica's NDSolve procedure.
Solution 2.

Example 3.  Solve the Lorenz I. V. P.

Solution 3.

Example 4.  Solve the Lorenz I. V. P.

Use Mathematica's NDSolve procedure.
Solution 4.

Rössler Attractor

In 1976 the Swiss mathematician Otto Rössler was studying oscillations in chemical reactions and discovered another set of equations with attractor.  They also are involved with the study of  Navier-Stokes equations.

,

,

.

Rössler is acclaimed to use the parameters where a = 0.2, b = 0.2, and c = 5.7.  Screw chaos occurs when the parameter values are a = 0.343,  b = 1.82  and  c = 9.75.  Rössler wanted to find a minimum system which would exhibit chaos. This system of equations looks easier than the Lorenz system, but it is harder to analyze.

Example 5.  Solve the Rössler attractor I. V. P.
with  ,
with  ,

with  .
Solution 5.

Example 6.  Solve the Rössler attractor I. V. P.
with  ,
with  ,
with  .
Solution 6.

Lorenz Attractor  Lorenz Attractor  Internet hyperlinks to web sites and a bibliography of articles.

(c) John H. Mathews 2004