מערכת לורנץ

מערכת משוואות דפרנציאליות רגילות, שנחקרו לראשונה על ידי המטאורולוג אדוארד לורנץ

בתורת הכאוס, מערכת לורנץ היא מערכת משוואות דפרנציאליות רגילות, שנחקרו לראשונה על ידי המטאורולוג אדוארד לורנץ. מערכת לורנץ מפורסמת עקב ההתנהגות הכאוטית של פתרונותיה, כתלות בתנאי ההתחלה והפרמטרים שלה המכונים מושכי לורנץ.

דוגמה למסלול במרחב הפאזה של מושך לורנץ עם הנתונים בהם השתמש לורנץ במחקרו הראשוני[1]

ההבנה של רגישות מערכת לורנץ לתנאי התחלה הובילה לשינוי בגישה של חוקרים לבעיות הנובעות ממערכות פיזיקליות שתנאי ההתחלה בהן אינם יודעים במדויק, או אינם יכולים להינתן במדויק. כך למשל, הביטוי הפופולרי של אפקט הפרפר מייצג מקרה בו משק כנפי פרפר עשוי ליצור שינויים קטנים באטמוספירה שבסופו של דבר יגרמו להופעת סופת טורנדו.

המשמעות העמוקה של הבנה זו היא שניסיון לגזור פתרון מדויק של מערכת פיזיקלית דינמית לאורך זמן נידון מראש לכישלון, גם בקני מידה גדולים בהם אפקטים קוונטיים אינם משמעותיים.

הגדרהעריכה

בשנת 1963, יצר המטאורולוג אדוארד לורנץ מודל מתמטי מופשט של הולכת החום באטמוספירה הנגזר מקירוב הציפה של בוסינסק[1].

המודל מתייחס לאטמוספירה כשכבת נוזל דליל דו ממדית, המחוממת בצורה אחידה על השפה התחתונה, ומקוררת בצורה אחידה משפתה העליונה.

המודל יוצר מערכת משוואות דיפרנציאליות רגילות הכוללת שלוש משוואות:

 

כאשר x הוא קצב הולכת החום, y הוא השינוי בטמפרטורה בציר האופקי וz הוא השינוי בטמפרטורה בציר האנכי.

הפרמטר   פרופורציוני למספר פרנטל של האטמוספירה,   לממדים הפיזיים של המערכת ו  למספר ריילי של האטמוספירה[2].

זוהי מערכת משוואות תלת-ממדית, דטרמיניסטית ולא ליניארית, אשר אין בפתרון שלה שום אלמנט רנדומלי.

אנליזהעריכה

 
המושך המוזר של לורנץ

לורנץ השתמש בגדלים המדידים   לתיאור המערכת המטאורולוגית, וקיבל כי למערכת פתרונות שמתנגים בצורה כאוטית עבור ערכים אלו.

כמעט כל הפתרונות של מערכת לורנץ עם הפרמטרים הנ"ל משתייכים לסט בעל שלוש נקודות שבת שנקרא מושך לורנץ. מושך לורנץ הוא מושך מוזר, כלומר מושך בעל מבנה פרקטלי עם מימד האוסדורף של  .

אף שהמשוואות המגדירות את מושך לורנץ הן פשוטות, הוכחת היותו מושך מוזר אינה קלה בכלל, ולמעשה היא הבעיה ה14 מבעיות שמייל, ונפתרה רק בשנת 2002 על ידי המתמטיקאי האוסטרלי וורוויק טאקר.

בפשטות, עבור   מתקבלת נקודת שבת יציבה יחידה בראשית   והיא מהווה מושך גלובלי אשר כל המסלולים מתכנסים אליו. משמעות נקודה זו היא שאין כלל הולכת חום בנוזל.

בנקודה   מתרחש פיצול קלשון, המוסיף שתי נקודות קריטיות בנקודות  , אשר משמעותן היא הולכת חום בקצב אחיד אופקית ואנכית. נקודות אלה יציבות אם ורק אם  . מכיוון שחומרים פיזיקליים הם בעלי ערכי  , נקבל דרישה מחייבת כי  [3].

סימולציית MATLABעריכה

 1 % Solve over time interval [0,100] with initial conditions [1,1,1]
 2 % ''f'' is set of differential equations
 3 % ''a'' is array containing x, y, and z variables
 4 % ''t'' is time variable
 5 
 6 sigma = 10;
 7 beta = 8/3;
 8 rho = 28;
 9 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)];
10 [t,a] = ode45(f,[0 100],[1 1 1]); % Runge-Kutta 4th/5th order ODE solver
11 plot3(a(:,1),a(:,2),a(:,3))

קישורים חיצונייםעריכה

  מדיה וקבצים בנושא מערכת לורנץ בוויקישיתוף

הערות שולייםעריכה

  1. ^ 1 2 Edward N. Lorenz, Deterministic Nonperiodic Flow., JAtS 20, 1963-03, עמ' 130–148 doi: 10.1175/1520-0469(1963)0202.0.CO;2
  2. ^ Colin Sparrow, The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors, Applied Mathematical Sciences, 1982 doi: 10.1007/978-1-4612-5767-7
  3. ^ Hirsch, Morris W., 1933-, Differential equations, dynamical systems, and an introduction to chaos, Amsterdam: Elsevier/Academic Press, 2nd ed, 2004, ISBN 978-0-08-049114-1