Μετάβαση στο περιεχόμενο

Σύστημα Λόρεντζ

Από τη Βικιπαίδεια, την ελεύθερη εγκυκλοπαίδεια
Ο ελκυστής του Λόρεντζ

Το Σύστημα Λόρεντζ[1] είναι μια μορφοκλασματική δομή που αντιστοιχεί στη μακροχρόνια συμπεριφορά του ταλαντωτή Λόρεντζ. Ο ελκυστής δείχνει πώς οι διάφορες μεταβλητές του δυναμικού συστήματος εξελίσσονται με την πάροδο του χρόνου σε μια μη περιοδική τροχιά.

Το 1963, ο μετεωρολόγος Έντουαρντ Λόρεντζ ήταν ο πρώτος που τόνισε την πιθανή χαοτική φύση της μετεωρολογίας. Το μοντέλο Λόρεντζ, επίσης γνωστό ως δυναμικό σύστημα Λόρεντζ ή ταλαντωτής Λόρεντζ, είναι ένα απλουστευμένο πρότυπο των μετεωρολογικών φαινομένων που βασίζεται στη ρευστομηχανική. Το μοντέλο αυτό είναι ένα τρισδιάστατο δυναμικό σύστημα που παράγει χαοτική συμπεριφορά υπό ορισμένες συνθήκες.

Το πρότυπο του Λόρεντζ είχε σημαντικό αντίκτυπο στην ανάδειξη των πιθανών ορίων της ικανότητας πρόβλεψης της μακροπρόθεσμης κλιματικής και μετεωρολογικής εξέλιξης. Αποτελεί σημαντικό στοιχείο της θεωρίας ότι η ατμόσφαιρα των πλανητών και των αστέρων μπορούν να περιλαμβάνουν μια μεγάλη ποικιλία οιονεί περιοδικών καθεστώτων και υπόκεινται σε απότομες και φαινομενικά τυχαίες μεταβολές. Αποτελεί επίσης ένα χρήσιμο παράδειγμα στη θεωρία των δυναμικών συστημάτων, χρησιμεύοντας ως πηγή νέων μαθηματικών εννοιών[2].

Ο ελκυστής και οι σχετικές εξισώσεις δημοσιοποιήθηκαν το 1963 από τον Έντουαρντ Λόρεντζ.[3]

Από μαθηματική άποψη, η σύζευξη της γήινης ατμόσφαιρας με τον ωκεανό περιγράφεται από το σύστημα συζευγμένων μερικών διαφορικών εξισώσεων Ναβιέρ-Στόκες της ρευστομηχανικής. Αυτό το σύστημα εξισώσεων ήταν πολύ περίπλοκο για να επιλυθεί αριθμητικά με τους υπολογιστές που υπήρχαν την εποχή του Λόρεντζ. Έτσι, ο Λόρεντζ είχε την ιδέα να αναζητήσει ένα εξαιρετικά απλουστευμένο πρότυπο αυτών των εξισώσεων για τη μελέτη μιας συγκεκριμένης φυσικής κατάστασης: του φαινομένου της συναγωγής Ρέιλεϊ-Μπενάρ. Το αποτέλεσμα ήταν ένα δυναμικό διαφορικό σύστημα με τρεις μόνο βαθμούς ελευθερίας, το οποίο ήταν πολύ πιο απλό να ολοκληρωθεί αριθμητικά από τις αρχικές εξισώσεις.

Δυναμικό διαφορικό σύστημα Λόρεντζ

[Επεξεργασία | επεξεργασία κώδικα]

Αυτό το διαφορικό σύστημα γράφεται ως εξής[4]:

όπου σ, ρ και β είναι τρεις αυστηρά θετικές, σταθερές παράμετροι. Δεδομένου ότι το σύστημα προκύπτει από την απλοποίηση των εξισώσεων που διέπουν τη συναγωγή Rayleigh-Mπενάρ, οι παράμετροι ονομάζονται από τη φυσική τους προέλευση: σ είναι ο αριθμός Prandtl[5], και ρ, που καταχρηστικά αναφέρεται ως "αριθμός Rayleigh"[6], είναι στην πραγματικότητα ο λόγος του αριθμός Rayleigh Ra προς τον κρίσιμο αριθμό Rayleigh[6] Rac (Rac είναι η τιμή του Ra πάνω από την οποία το φυσικό σύστημα δεν μπορεί να παραμείνει αναλλοίωτο με την πάροδο του χρόνου, όπου επομένως υπόκειται σε κινήσεις συναγωγής).

Οι δυναμικές μεταβλητές x, y και z αντιπροσωπεύουν την κατάσταση του συστήματος σε κάθε χρονική στιγμή. Η φυσική ερµηνεία είναι η εξής: x(t) είναι ανάλογη της έντασης της κίνησης συναγωγής, y(t) είναι ανάλογη της διαφοράς θερµοκρασίας µεταξύ των ανοδικών και καθοδικών ρευµάτων και z(t) είναι ανάλογη της απόκλισης του κατακόρυφου προφίλ θερµοκρασίας από ένα γραµµικό προφίλ[7]..

Συχνά υποθέτουμε σ = 10 και β = 8/3, με το ρ να παραμένει μεταβλητό. Το σύστημα παρουσιάζει χαοτική συμπεριφορά για ρ = 28, μετά από κάποιο χρονικό διάστημα στο ψευδο-περιοδικό καθεστώς (της τάξης των 30 δευτερολέπτων). Ο χρόνος αυτός εξαρτάται από το ρ : μπορεί να φτάσει τα 180 δευτερόλεπτα για ρ = 24,3. Το κατώφλι για την εγκαθίδρυση του χαοτικού καθεστώτος είναι μεταξύ 24,29 και 24,30 για αυτές τις τιμές σ και β[8].

Σημεία ισορροπίας

[Επεξεργασία | επεξεργασία κώδικα]

Τα σημεία ισορροπίας ή σταθερά σημεία του συστήματος είναι οι σταθερές (x,y,z) λύσεις του διαφορικού συστήματος. Υπάρχουν ένα ή τρία τέτοια σημεία:

  • το σταθερό σημείο (0, 0, 0) υπάρχει ανεξάρτητα από τις τιμές των πραγματικών παραμέτρων σ, ρ και β. Είναι σταθερό όταν (σύγκλιση στην κατάσταση ηρεμίας, χωρίς συναγωγή), ασταθές διαφορετικά (συναγωγή) ,
  • τα δύο συμμετρικά σταθερά σημεία και υπάρχουν μόνο για την . Ανάλογα με τις τιμές των σ και ρ, είναι σταθερές (σταθερή συναγωγή) ή ασταθείς (περιοδικά κυμαινόμενη συναγωγή ή χαοτική).
Απορρόφηση 25.000 σημείων από έναν ελκυστή Λόρεντζ.

Ο ελκυστής Λόρεντζ ορίζεται ως το σύνολο των μακροχρόνιων τροχιών του παραπάνω δυναμικού συστήματος Λόρεντζ.

Πράγματι, όταν οι παράμετροι σ, ρ και β παίρνουν τις ακόλουθες τιμές: σ = 10, ρ = 28 και β = 8/3, το διαφορικό δυναμικό σύστημα Λόρεντζ παρουσιάζει έναν παράξενο ελκυστή με τη μορφή φτερών πεταλούδας, που απεικονίζεται στο διπλανό σχήμα.

Για σχεδόν όλες τις αρχικές συνθήκες (εκτός από εκείνες για τα σταθερά σημεία), η τροχιά του συστήματος πλησιάζει γρήγορα τον ελκυστή, με την περιφορά της τροχιάς αρχικά να τυλίγεται γύρω από τη μία πτέρυγα, στη συνέχεια να πηδά από τη μία πτέρυγα στην άλλη για να αρχίσει να τυλίγεται γύρω από την άλλη πτέρυγα, και ούτω καθεξής, με έναν προφανώς ακανόνιστο τρόπο.

Η ύπαρξη ενός παράξενου ελκυστή για ορισμένες τιμές παραμέτρων είχε πιθανολογηθεί από τον Έντουαρντ Λόρεντζ ήδη από το 1963 βάσει αριθμητικών προσομοιώσεων. Ωστόσο, έπρεπε να περιμένουμε μέχρι το 2001 για μια αυστηρή απόδειξη του γεγονότος αυτού από τον Ουόρικ Τάκερ[9].

Σε αυτές τις περιπτώσεις, ο ελκυστής είναι ένα φράκταλ με διάσταση Χάουστορφ μεταξύ 2 και 3.[10]

Προκειμένου να μελετηθεί μαθηματικά ο ελκυστής Λόρεντζ, εμφανίστηκαν στην επιστημονική βιβλιογραφία αυστηρά αλλά διακριτά πρότυπα [11]. Το κύριο πλεονέκτημά τους ήταν ότι ήταν ευκολότερο να μελετηθούν και, κυρίως, ότι μπορούσαμε να είμαστε σίγουροι για την ύπαρξή τους (αν και δεν γνωρίζαμε αν ο ελκυστής Λόρεντζ περιγράφεται σωστά από αυτά τα πρότυπα, βλ. παραπάνω).

Προσομοίωση στο Julia.

Προσομοίωση στο Julia

[Επεξεργασία | επεξεργασία κώδικα]
using Plots
# define the Lorenz attractor
@kwdef mutable struct Lorenz
    dt::Float64 = 0.02
    σ::Float64 = 10
    ρ::Float64 = 28
    β::Float64 = 8/3
    x::Float64 = 2
    y::Float64 = 1
    z::Float64 = 1
end

function step!(l::Lorenz)
    dx = l.σ * (l.y - l.x);         l.x += l.dt * dx
    dy = l.x * (l.ρ - l.z) - l.y;   l.y += l.dt * dy
    dz = l.x * l.y - l.β * l.z;     l.z += l.dt * dz
end

attractor = Lorenz()

# initialize a 3D plot with 1 empty series
plt = plot3d(
    1,
    xlim = (-30, 30),
    ylim = (-30, 30),
    zlim = (0, 60),
    title = "Lorenz Attractor",
    marker = 2,
)

# build an animated gif by pushing new points to the plot, saving every 10th frame
@gif for i=1:1500
    step!(attractor)
    push!(plt, attractor.x, attractor.y, attractor.z)
end every 10

Προσομοίωση με το MATLAB

[Επεξεργασία | επεξεργασία κώδικα]
% Resolver para el intervalo de tiempo [0,100] con condiciones iniciales [1,1,1]
% ''f'' es un conjunto de ecuaciones diferenciales
% ''a'' es un arreglo que contiene variables x, y, z
% ''t'' es la variable de tiempo

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]);     % Solución de EDO de Runge-Kutta de 4.º/5.º orden
plot3(a(:,1),a(:,2),a(:,3))

Προσομοίωση με το Python

[Επεξεργασία | επεξεργασία κώδικα]
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  # Desempaqueta el vector de estado
    return sigma * (y - x), x * (rho - z) - y, x * y - beta * z  # Derivadas

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.add_subplot(projection='3d')
ax.plot(states[:, 0], states[:, 1], states[:, 2])
plt.show()
  1. «Σύστημα Lorenz». www.hellenicaworld.com. Ανακτήθηκε στις 13 Δεκεμβρίου 2023. 
  2. «L'attracteur de Lorenz, paradigme du chaos - Séminaire Poincaré - CNRS» (PDF). 
  3. Gulick, Denny (26 Απριλίου 2012). Encounters with Chaos and Fractals. CRC Press. ISBN 978-1-4665-5875-5. 
  4. Magnit?ski?, Nikola? Aleksandrovich· Sidorov, Sergey Vasilevich (2006). New Methods for Chaotic Dynamics. World Scientific. ISBN 978-981-256-817-5. 
  5. «Prandtl Number - an overview | ScienceDirect Topics». www.sciencedirect.com. Ανακτήθηκε στις 13 Δεκεμβρίου 2023. 
  6. 6,0 6,1 «Rayleigh Number - an overview | ScienceDirect Topics». www.sciencedirect.com. Ανακτήθηκε στις 13 Δεκεμβρίου 2023. 
  7. Tels que définis dans Lorenz 1963, σελ. 135.
  8. La valeur du seuil peut dépendre des paramètres de résolution numérique du système.
  9. «A Rigorous ODE solver and Smale's 14th Problem». archive.wikiwix.com. Ανακτήθηκε στις 14 Δεκεμβρίου 2023. 
  10. Grassberger, Peter; Procaccia, Itamar (1983-10-01). «Measuring the strangeness of strange attractors». Physica D: Nonlinear Phenomena 9 (1): 189–208. doi:10.1016/0167-2789(83)90298-1. ISSN 0167-2789. https://www.sciencedirect.com/science/article/pii/0167278983902981. 
  11. J. Guckenheimer et R. F. Williams, « Structural stability of Lorenz attractors », Inst. Hautes études Sci. Publ. Math., vol. 50,‎ 1979, p. 59-72.