Lorenz-System - Lorenz system

Eine Beispiellösung im Lorenz-Attraktor für ρ = 28, σ = 10 und β = 8/3

Das Lorenz-System ist ein System gewöhnlicher Differentialgleichungen, das zuerst von Edward Lorenz untersucht wurde . Es zeichnet sich durch chaotische Lösungen für bestimmte Parameterwerte und Anfangsbedingungen aus. Insbesondere ist der Lorenz-Attraktor eine Menge chaotischer Lösungen des Lorenz-Systems. In populären Medien stammt der " Schmetterlingseffekt " aus den realen Implikationen des Lorenz-Attraktors, dh in jedem physikalischen System, in Ermangelung perfekter Kenntnis der Anfangsbedingungen (sogar die winzige Störung der Luft durch einen Schmetterlingsflattern). seine Flügel) wird unsere Fähigkeit, seinen zukünftigen Kurs vorherzusagen, immer versagen. Dies unterstreicht, dass physikalische Systeme auch ohne Quanteneffekte vollständig deterministisch und dennoch von Natur aus unberechenbar sein können. Die Form des Lorenz-Attraktors selbst kann, wenn sie grafisch aufgetragen wird, auch einem Schmetterling ähneln.

Überblick

1963 entwickelte Edward Lorenz mit Hilfe von Ellen Fetter ein vereinfachtes mathematisches Modell für die atmosphärische Konvektion . Das Modell ist ein System von drei gewöhnlichen Differentialgleichungen, die heute als Lorenz-Gleichungen bekannt sind:

Die Gleichungen beziehen sich auf die Eigenschaften einer zweidimensionalen Flüssigkeitsschicht, die von unten gleichmäßig erwärmt und von oben gekühlt wird. Insbesondere beschreiben die Gleichungen die Änderungsrate von drei Größen in Bezug auf die Zeit: ist proportional zur Konvektionsrate, zur horizontalen Temperaturänderung und zur vertikalen Temperaturänderung. Die Konstanten , , und sind Systemparameter, die proportional zur Prandtl-Zahl , Rayleigh-Zahl und bestimmten physikalischen Abmessungen der Schicht selbst sind.

Die Lorenz-Gleichungen treten auch in vereinfachten Modellen für Laser , Dynamos , Thermosiphons , bürstenlose Gleichstrommotoren , elektrische Schaltungen , chemische Reaktionen und Vorwärtsosmose auf . Die Lorenz-Gleichungen sind auch die maßgebenden Gleichungen im Fourier-Raum für das Malkus-Wasserrad . Das Malkus-Wasserrad zeigt eine chaotische Bewegung, bei der sich seine Drehung, anstatt sich mit konstanter Geschwindigkeit in eine Richtung zu drehen, beschleunigt, verlangsamt, stoppt, die Richtung ändert und auf unvorhersehbare Weise zwischen Kombinationen solcher Verhaltensweisen hin und her oszilliert.

Aus technischer Sicht ist das Lorenz-System nichtlinear , nichtperiodisch, dreidimensional und deterministisch . Die Lorenz-Gleichungen waren Gegenstand von Hunderten von Forschungsartikeln und mindestens einer Studie in Buchlänge.

Analyse

Normalerweise geht man davon aus, dass die Parameter , , und positiv sind. Lorenz verwendet , um die Werte , und . Das System zeigt bei diesen (und nahegelegenen) Werten ein chaotisches Verhalten.

Wenn es dann nur einen Gleichgewichtspunkt gibt, der im Ursprung liegt. Dieser Punkt entspricht keiner Konvektion. Alle Bahnen konvergieren zum Ursprung, der ein globaler Attraktor ist , wenn .

Eine Heugabelgabelung tritt bei auf , und zwei weitere kritische Punkte erscheinen bei: und Diese entsprechen einer stetigen Konvektion. Dieses Gleichgewichtspunktpaar ist nur dann stabil, wenn

die nur positiv gelten kann, wenn . Beim kritischen Wert verlieren beide Gleichgewichtspunkte durch eine unterkritische Hopf-Bifurkation an Stabilität .

Wenn , , und , hat das Lorenz-System chaotische Lösungen (aber nicht alle Lösungen sind chaotisch). Fast alle Anfangspunkte tendieren zu einer invarianten Menge – dem Lorenz-Attraktor – einem seltsamen Attraktor , einem Fraktal und einem selbsterregten Attraktor in Bezug auf alle drei Gleichgewichte. Seine Hausdorff-Dimension wird von oben durch die Lyapunov-Dimension (Kaplan-Yorke-Dimension) auf 2,06 ± 0,01 geschätzt , und die Korrelationsdimension wird auf 2,05 ± 0,01 geschätzt. Die exakte Lyapunov-Dimensionsformel des globalen Attraktors lässt sich analytisch unter klassischen Einschränkungen der Parameter ermitteln:

Der Lorenz-Attraktor ist schwer zu analysieren, aber die Wirkung der Differentialgleichung auf den Attraktor wird durch ein ziemlich einfaches geometrisches Modell beschrieben. Der Beweis, dass dies tatsächlich der Fall ist, ist das vierzehnte Problem auf der Liste der Probleme von Smale . Dieses Problem wurde 2002 von Warwick Tucker als erstes gelöst .

Für andere Werte von zeigt das System verknotete periodische Bahnen an. Damit wird zum Beispiel ein T (3,2) Torusknoten .

Beispiellösungen des Lorenz-Systems für verschiedene Werte von ρ
Lorenz Ro14 20 41 20-200px.png Lorenz Ro13-200px.png
ρ = 14, σ = 10, β = 8/3 (Vergrößern) ρ = 13, σ = 10, β = 8/3 (Vergrößern)
Lorenz Ro15-200px.png Lorenz Ro28-200px.png
ρ = 15, σ = 10, β = 8/3 (Vergrößern) ρ = 28, σ = 10, β = 8/3 (Vergrößern)
Für kleine Werte von ρ ist das System stabil und entwickelt sich zu einem von zwei Fixpunktattraktoren. Wenn ρ größer als 24,74 ist, werden die Fixpunkte zu Repulsoren und die Flugbahn wird von ihnen auf sehr komplexe Weise abgestoßen.
Sensible Abhängigkeit von der Anfangsbedingung
Zeit t = 1 (Vergrößern) Zeit t = 2 (Vergrößern) Zeit t = 3 (Vergrößern)
Lorenz caos1-175.png Lorenz caos2-175.png Lorenz caos3-175.png
Diese Abbildungen – erstellt mit ρ = 28, σ = 10 und β = 8/3 – zeigen drei Zeitabschnitte der 3-D-Entwicklung zweier Trajektorien (eine in Blau, die andere in Gelb) im Lorenz-Attraktor beginnend bei zwei anfänglichen Punkte, die sich nur um 10 −5 in der x- Koordinate unterscheiden. Anfangs scheinen die beiden Trajektorien zusammenzufallen (nur die gelbe ist zu sehen, da sie über die blaue gezeichnet wird), aber nach einiger Zeit ist die Divergenz offensichtlich.

Verbindung zum Zeltplan

Eine Nachbildung der Ergebnisse von Lorenz, die auf Mathematica erstellt wurde . Punkte über der roten Linie entsprechen den Schaltkeulen des Systems.

In Abbildung 4 seiner Arbeit hat Lorenz den vom System erreichten relativen Maximalwert in z-Richtung gegen den vorherigen relativen Maximalwert in z-Richtung aufgetragen. Dieses Verfahren wurde später als Lorenz-Karte bekannt (nicht zu verwechseln mit einem Poincaré-Plot , der die Schnittpunkte einer Trajektorie mit einer vorgegebenen Oberfläche darstellt). Das resultierende Grundstück hat eine Form, die der Zeltkarte sehr ähnlich ist . Lorenz fand auch heraus, dass das System auf die nächste Keule umschaltet, wenn der maximale z-Wert über einem bestimmten Grenzwert liegt. In Kombination mit dem bekannten Chaos der Zeltkarte zeigte er, dass das System chaotisch zwischen den beiden Lappen wechselt.

Simulationen

MATLAB-Simulation

% Solve over time interval [0,100] with initial conditions [1,1,1]
% ''f'' is set of differential equations
% ''a'' is array containing x, y, and z variables
% ''t'' is time variable

sigma = 10;
beta = 8/3;
rho = 28;
f = @(t,a) [-sigma*a(1) + sigma*a(2); rho*a(1) - a(2) - a(1)*a(3); -beta*a(3) + a(1)*a(2)];
[t,a] = ode45(f,[0 100],[1 1 1]);     % Runge-Kutta 4th/5th order ODE solver
plot3(a(:,1),a(:,2),a(:,3))

Mathematica-Simulation

Standardweg:

tend = 50;
eq = {x'[t] == σ (y[t] - x[t]), 
      y'[t] == x[t] (ρ - z[t]) - y[t], 
      z'[t] == x[t] y[t] - β z[t]};
init = {x[0] == 10, y[0] == 10, z[0] == 10};
pars = {σ->10, ρ->28, β->8/3};
{xs, ys, zs} = 
  NDSolveValue[{eq /. pars, init}, {x, y, z}, {t, 0, tend}];
ParametricPlot3D[{xs[t], ys[t], zs[t]}, {t, 0, tend}]

Weniger ausführlich:

lorenz = NonlinearStateSpaceModel[{{σ (y - x), x (ρ - z) - y, x y - β z}, {}}, {x, y, z}, {σ, ρ, β}];
soln[t_] = StateResponse[{lorenz, {10, 10, 10}}, {10, 28, 8/3}, {t, 0, 50}];
ParametricPlot3D[soln[t], {t, 0, 50}]

Python-Simulation

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
from mpl_toolkits.mplot3d import Axes3D

rho = 28.0
sigma = 10.0
beta = 8.0 / 3.0

def f(state, t):
    x, y, z = state  # Unpack the state vector
    return sigma * (y - x), x * (rho - z) - y, x * y - beta * z  # Derivatives

state0 = [1.0, 1.0, 1.0]
t = np.arange(0.0, 40.0, 0.01)

states = odeint(f, state0, t)

fig = plt.figure()
ax = fig.gca(projection="3d")
ax.plot(states[:, 0], states[:, 1], states[:, 2])
plt.draw()
plt.show()

Herleitung der Lorenz-Gleichungen als Modell für atmosphärische Konvektion

Die Lorenz-Gleichungen sind aus der Oberbeck-Boussinesq-Näherung der Gleichungen abgeleitet , die die Flüssigkeitszirkulation in einer flachen Flüssigkeitsschicht beschreiben, die von unten gleichmäßig erhitzt und von oben gleichmäßig gekühlt wird. Diese Flüssigkeitszirkulation wird als Rayleigh-Bénard-Konvektion bezeichnet . Es wird angenommen, dass die Flüssigkeit in zwei Dimensionen (vertikal und horizontal) mit periodischen rechteckigen Randbedingungen zirkuliert.

Die partiellen Differentialgleichungen, die die Strömungsfunktion und die Temperatur des Systems modellieren , werden einer spektralen Galerkin-Näherung unterzogen : Die hydrodynamischen Felder werden in Fourier-Reihen erweitert, die dann stark auf einen einzigen Term für die Strömungsfunktion und zwei Terme für die Temperatur verkürzt werden. Dies reduziert die Modellgleichungen auf einen Satz von drei gekoppelten, nichtlinearen gewöhnlichen Differentialgleichungen. Eine detaillierte Herleitung findet sich beispielsweise in nichtlinearen Dynamiktexten. Das Lorenz-System ist eine reduzierte Version eines größeren Systems, das zuvor von Barry Saltzman untersucht wurde.

Lösung von Smales 14. Problem

Smales 14. Problem lautet: „Haben die Eigenschaften des Lorenz-Attraktors die eines seltsamen Attraktors ?“, wurde 2002 von Warwick Tucker bejaht . Um dieses Ergebnis zu beweisen, verwendete Tucker rigorose numerische Methoden wie Intervallarithmetik und Normalformen . Zunächst definiert Tucker einen Querschnitt , der von den Strömungstrajektorien quer geschnitten wird. Daraus kann man die First-Return-Map definieren , die jedem den Punkt zuordnet, an dem sich die Trajektorie von First schneidet .

Dann wird der Beweis in drei Hauptpunkte aufgeteilt, die bewiesen sind und die Existenz eines seltsamen Attraktors implizieren. Die drei Punkte sind:

  • Unter der First-Return-Karte existiert eine Regionsinvariante, d. h
  • Die Rückkehrabbildung lässt ein vorwärtsinvariantes Kegelfeld zu
  • Vektoren innerhalb dieses invarianten Kegelfeldes werden durch die Ableitung der Rückkehrabbildung gleichförmig erweitert .

Um den ersten Punkt zu beweisen, bemerken wir, dass der Querschnitt von zwei Bögen geschnitten wird, die von (siehe ) gebildet werden. Tucker bedeckt die Lage dieser beiden Bögen durch kleine Rechtecke , die Vereinigung dieser Rechtecke ergibt . Das Ziel ist nun zu beweisen, dass für alle Punkte in , der Fluss die Punkte in , in zurückbringt . Um dies zu tun, nehmen wir einen Plan unten in einem Abstand klein ist , dann durch das Zentrum unter dem und Euler Integrationsverfahren verwendet wird , kann man sich ausrechnen , wo die Strömung bringt in denen gibt uns einen neuen Punkt . Dann kann man mit Hilfe der Taylor-Entwicklung abschätzen, wo die Punkte in abgebildet werden. Dies gibt uns ein neues Rechteck mit dem Mittelpunkt von . Somit wissen wir, dass alle Punkte in in abgebildet werden . Das Ziel ist , diese Methode rekursiv zu tun , bis der Ablauf zu zurückkommt und wir erhalten ein Rechteck in so dass wir wissen , dass . Das Problem ist, dass unsere Schätzung nach mehreren Iterationen ungenau werden kann, daher teilt Tucker in kleinere Rechtecke auf und wendet den Prozess dann rekursiv an. Ein weiteres Problem besteht darin, dass der Fluss bei der Anwendung dieses Algorithmus „horizontaler“ wird (siehe ), was zu einer dramatischen Zunahme der Ungenauigkeit führt. Um dies zu verhindern, ändert der Algorithmus die Ausrichtung der Querschnitte, sodass sie entweder horizontal oder vertikal werden.

Beiträge

Lorenz würdigt die Beiträge von Ellen Fetter in seinem Beitrag, die für die numerischen Simulationen und Zahlen verantwortlich ist. Außerdem half Margaret Hamilton bei den ersten numerischen Berechnungen, die zu den Ergebnissen des Lorenz-Modells führten.

Galerie

Siehe auch

Anmerkungen

Verweise

Weiterlesen

Externe Links