בעיה n-גופית

בפיזיקה, הבעיה ה--גופית היא הבעיה של חיזוי תנועה של קבוצת עצמים שמימיים המקיימים אינטראקציית כבידה זה עם זה.<ref></ref> המוטיבציה לפתרון הבעיה היה הרצון להבין את תנועות השמש, הירח וכוכבי הלכת והכוכבים. במאה ה-20, הבנת הדינמיקה של מערכות כוכבי צביר כדוריות הפכה לבעיה -גופית חשובה.<ref></ref> הבעיה ה--גופית בתורת היחסות הכללית קשה הרבה יותר לפתרון, גם בשל המורכבות המתמטית של התיאוריה וגם עקב גורמים נוספים כמו עיוותי זמן ומרחב.
== היסטוריה ==
באמצעות נתונים של שלושה מיקומים של כוכב לכת במסלולו – מיקומים שקיבל סר אייזק ניוטון מהאסטרונום ג'ון פלמסטיד <ref>. </ref> – ניוטון הצליח לנסח משוואה, שחוזה את תנועתו של כוכב לכת; כלומר, נותנת את מאפייני המסלול שלו: מיקום, קוטר מסלול, זמן מחזור ומהירות מסלול.<ref name=":0">See </ref> אך במהרה התגלה כי משוואות התנועה הללו לא חזו בדיוק חלק מהמסלולים.<ref></ref> ניוטון הבין שחוסר הדיוק נובע מכוחות הכבידה שכוכבי הלכת מפעילים זה על זה.

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

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

מציאת הפתרון הכללי של הבעיה ה--גופית נחשבה חיונית ועם זאת מאתגרת מאוד. ואכן, בסוף המאה ה-19, מלך שוודיה אוסקר השני, הכריז על פרס למי שיציע פתרון לבעיה.<ref></ref>
== ניסוח כללי ==
הבעיה ה--גופית מניחה מסות נקודתיות במערכת ייחוס אינרציאלית במרחב תלת־ממדי, שנעות תחת השפעת כוחות משיכה הדדית של כבידה. לכל מסה יש וקטור מיקום . על פי החוק השני של ניוטון, מסה כפול תאוצת הגוף שווה לסך כל הכוחות שפועלים עליה. על פי חוק הכבידה העולמי של ניוטון כוח המשיכה ששתי מסות ו מפעילות זו על זו הוא

<math display=block>\mathbf{F}_{ij} = \frac{G m_i m_j}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^2} \, \frac{\left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|} = \frac{G m_i m_j \left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^3},</math>

כאשר הוא קבוע הכבידה העולמי ו <math>\left\| \mathbf{q}_j - \mathbf{q}_i\right\|</math> המרחק בין המסות.

סיכום הכוחות על כל המסות נותן את משוואות התנועה ל-n הגופים<math display="block"> m_i \frac{d^2\mathbf{q}_i}{dt^2}
 = \sum_{j=1 \atop j \ne i}^n \frac{G m_i m_j \left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^3}
 = - \frac{\partial U}{\partial \mathbf{q}_i}
 </math>

כאשר <math> U </math> היא האנרגיה הפוטנציאלית של המערכת<math display="block">U = - \sum_{1 \le i < j \le n} \frac{G m_i m_j}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|}.</math>

עם הגדרת המומנטום<math> \mathbf{p}_i = m_i \frac{d\mathbf{q}_i}{dt} </math> משוואות התנועה ההמילטוניות לבעיה ה-n-גופית הופכות ל<ref name=":1"></ref><math display="block">\frac{d\mathbf{q}_i}{dt} = \frac{\partial H}{\partial \mathbf{p}_i} \qquad \frac{d\mathbf{p}_i}{dt} = - \frac{\partial H}{\partial \mathbf{q}_i}, </math>כאשר ההמילטוניאן הוא

<math display="block">H= T + U</math>ו T האנרגיה הקינטית הכוללת של המערכת<math display="block">T = \sum_{i=1}^n \frac{\left\| \mathbf{p}_i \right\|^2}{2m_i}.</math>

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

סימטריות בבעיה ה--גופית מספקות קבועי תנועה, באופן שמפשט את הבעיה.<ref name=":1" /> היות שהבעיה סימטרית להזזה מתקיים חוק שימור התנע הכולל של המערכת, ולכן מרכז המסה <math display="block">\mathbf{C} = \frac{\displaystyle\sum_{i=1}^n m_i \mathbf{q}_i}{\displaystyle\sum_{i=1}^n m_i} </math> נע במהירות קצובה ובקו ישר. כלומר , כאשר היא המהירות הקווית ו- היא המיקום ההתחלתי. קבועי התנועה ו - מהווים סה"כ שישה קבועי תנועה. סימטריה סיבובית גורמת לכך שהתנע הזוויתי הכולל הוא קבוע <math display="block">\mathbf{A} = \sum_{i=1}^n \mathbf{q}_i \times \mathbf{p}_i,</math> כאשר <math>\times</math> הוא אופרטור המכפלה הווקטורית. שלושת המרכיבים של התנע הזוויתי הכולל מהווים שלושה קבועי תנועה נוספים. קבוע התנועה האחרון מתקבל משימור האנרגיה . מכאן שלכל בעיה -גופית יש עשרה קבועי תנועה.

מכיוון ש - ו - הם פונקציות הומוגניות מדרגה 2 ו-1, בהתאמה, משוואות התנועה אינווריאנטיות לכיול: כלומר, אם הוא פתרון, אז גם הוא פתרון לכל .<ref name="Chenciner 2007">Chenciner 2007</ref>

== מקרים מיוחדים ==
=== בעיית שני הגופים ===

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

בעיית קפלר () נפתרה על ידי יוהאן ברנולי (1667–1748) באמצעות התיאוריה הקלאסית (ולא על ידי ניוטון) על ידי הנחה שמסת הנקודה העיקרית הייתה קבועה;<ref name="Bate, Mueller, and White">See Bate, Mueller, and White, Chapter 1: "Two-Body Orbital Mechanics", pp. 1–49. These authors were from the Department of Astronautics and Computer Science, United States Air Force Academy. Their textbook is not filled with advanced mathematics.</ref> בתנועה של שני גופים, למשל השמש וכדור הארץ, כשהשמש קבועה, ואז: <math display="block">\begin{align}
m_1 \mathbf{a}_1 &= \frac{Gm_1m_2}{r_{12}^3}(\mathbf{r}_2-\mathbf{r}_1) &&\quad\text{Sun–Earth} \\
m_2 \mathbf{a}_2 &= \frac{Gm_1m_2}{r_{21}^3}(\mathbf{r}_1-\mathbf{r}_2) &&\quad\text{Earth–Sun}
\end{align}</math>המשוואה המתארת את תנועת המסה ביחס למסה מתקבלת מההפרש של שתי המשוואות הללו, ולאחר ביטול איברים משותפים מתקבל: <math display="block">\mathbf{\alpha} + \frac{\eta}{r^3} \mathbf{r} = \mathbf{0}</math> כאשר
* הוא וקטור המיקום של 
* הוא וקטור התאוצה האוילרית <math>\frac{d^2 \mathbf{r}} {dt^2}</math>
*<math>\eta = G(m_1 + m_2)
</math>
המשוואה <math> \mathbf{\alpha} + \frac{\eta}{r^3} \mathbf{r} = \mathbf{0} </math> היא המשוואה הדיפרנציאלית היסודית לבעיה הדו-גופית תחת שדה כבידה.

=== בעיית שלושת הגופים ===

סעיף זה מתייחס לפתרון הבעיה ה- גופית לאחר שנעשו הנחות מפשטות.

בעבר לא היה ידוע הרבה על הבעיה ה--גופית עבור .<ref>See Leimanis and Minorsky's historical comments.</ref> המקרה התלת-גופי () נחקר ביותר. ניסיונות רבים קודמים להבין את בעיית שלושת הגופים היו כמותיים, ומטרתם הייתה למצוא פתרונות מפורשים למצבים ספציפיים.
* ב-1687 פרסם אייזק ניוטון ב-Principia את השלבים הראשונים בחקר בעיית התנועות של שלושה גופים הכפופים למשיכה ההדדית של הכבידה שלהם.
* בשנת 1767, מצא אוילר תנועות קולינאריות, שבהן שלושה גופים נעים לאורך קו ישר קבוע. בעיית שלושת הגופים של אוילר היא מקרה מיוחד שבו שניים מהגופים מקובעים במרחב (אין לבלבל זאת עם בעיית שלושת הגופים המוגבלת מעגלית, שבה שני הגופים המסיביים מתארים מסלול מעגלי).
* בשנת 1772 גילה לגראנז' שני סוגים של פתרון מחזורי, כל אחד עבור שלושה גופים בכל מסה. בסוג אחד, הגופים נמצאים על קו ישר מסתובב. בסוג השני הגופים נמצאים על קודקודיו של של משולש שווה-צלעות מסתובב. בשני המקרים המסלולים של הגופים יהיו חתכים חרוטיים. פתרונות אלה הובילו למחקר של תצורות מרכזיות, שעבורן עבור קבוע כלשהו .
* מחקר גדול על המערכת שמש-ארץ-ירח בוצע על ידי שארל-אוג'ן דלאוני. העבודה מרמזת על הכאוטיות של הבעיה.
* בשנת 1917, פורסט ריי מולטון את ספרו מבוא למכניקה שמימית עם פתרון לבעיית שלושת הגופים המצומצמת.<ref></ref>

=== בעיית ארבעת הגופים ===
בהשראת בעיית שלושת הגופים המוגבלת המעגלית, ניתן לפשט מאוד את בעיית ארבעת הגופים על ידי התחשבות בגוף קטן יותר כבעל מסה קטנה בהשוואה לשלושת הגופים המאסיביים האחרים, שבתורם מתארים מסלולים מעגליים. זו ידועה כבעיית ארבעת הגופים המוגבלת הדו-מעגלית (הידועה גם כדגם דו-מעגלי).<ref></ref> ניסוח זה היה רלוונטי מאוד באסטרודינמיקה, בעיקר למודל של מסלולי חלליות במערכת כדור הארץ-ירח בתוספת משיכה הכבידה של השמש.

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

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

=== פתרון בטור חזקות ===
אחת הדרכים לפתור את הבעיה ה--גופית היא באמצעות "טור טיילור"<ref></ref>

נתחיל בהגדרת מערכת משוואות דיפרנציאליות של הבעיה ה--גופית:<math display="block">\frac{d^2\mathbf{x}_i(t)}{dt^2}=G \sum_{k=1 \atop k\neq i}^n \frac{m_k \left(\mathbf{x}_k(t)-\mathbf{x}_i(t)\right)}{\left|\mathbf{x}_k(t)-\mathbf{x}_i(t)\right|^{3}},</math>כיוון ש <math>\mathbf{x_i} (t0)</math> ו <math>\frac{d\mathbf{x_i}(t_0)}{dt}</math> נתונים כתנאי התחלה, גם הנגזרות השניות ידועות ב-<math>t_0</math>. הגזירה של <math>\frac{d^2\mathbf{x}_i(t)}{dt^2}</math> ביחס ל-t נותן את <math>\frac{d^3\mathbf{x}_i(t)}{dt^3}</math> שגם אותו ניתן לחשב ב-<math>t_0</math>, וכך ניתן לחשב באופן איטרטיבי את טור טיילור המתאר את ההתפתחות בזמן של המערכת.

== סימולציה נומרית ==
אמנם ישנם פתרונות אנליטיים זמינים לבעיית שני הגופים הקלאסית (כלומר הלא-יחסותית) ולקונפיגורציות ספציפיות עם , אך באופן כללי יש לפתור או לדמות בעיות -גופיות באמצעות שיטות נומריות.<ref name=":2"></ref>

=== גופים מעטים ===
עבור מספר קטן של גופים, ניתן לפתור בעיה -גופית באמצעות שיטות ישירות, הנקראות גם שיטות חלקיקים-חלקיקים. שיטות אלו משלבות באופן נומרי את משוואות התנועה הדיפרנציאליות. אינטגרציה נומרית לבעיה זו יכולה להיות מאתגרת מכמה סיבות: ראשית, פוטנציאל הכבידה הוא סינגולרי, ושואף לאינסוף כשהמרחק בין שני חלקיקים שואף לאפס. ניתן "לרכך" את פוטנציאל הכבידה על ידי הוספת איבר חיובי קטן במכנה כדי להסיר את הסינגולריות במרחקים קטנים:<ref name=":2" /><math display="block">U_\varepsilon = \sum_{1 \le i < j \le n} \frac{G m_i m_j}{ \sqrt{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^2 + \varepsilon^2} }</math> שנית, באופן כללי עבור , הבעיה ה--גופית היא כאוטית,<ref></ref> כלומר אפילו שגיאות קטנות באינטגרציה עשויות לגדול אקספוננציאלית בזמן. שלישית, סימולציה עשויה להיות על פני משכי זמן גדולים של המודל (למשל מיליוני שנים) ושגיאות מספריות מצטברות ככל שגדל מספר צעדי האינטגרציה.

ישנן מספר טכניקות לצמצום שגיאות באינטגרציה נומרית.<ref name=":0" /> מערכות קואורדינטות מקומיות משמשות להתמודדות עם סקאלות שונות בבעיות מסוימות. למשל, לסימולציה של המערכת שמש-ארץ-ירח. שיטות וריאציה או כאלה המשתמשות בתורת ההפרעות יכולות להניב מסלולים אנליטיים מקורבים ששניתן לדייק באמצעות אינטגרציה נומרית.

=== גופים מרובים ===
שיטות ישירות המשתמשות באינטגרציה נומרית דורשות סדר גודל של <math>\frac{1}{2} n^2</math>חישובים להערכת האנרגיה הפוטנציאלית. ולכן סיבוכיות הזמן שלהן היא <math>O(n^2)</math>. עבור סימולציות עם גופים מרובים, פקטור <math>O(n^2)</math> הופך את החישובים לבלתי מעשיים, ולכן פותחו שיטות מקורבות המפחיתות את סיבוכיות הזמן ביחס לשיטות הישירות.<ref name=":2" />

== ראו גם==
* בעיית כבידה של שני גופים
* בעיית קפלר
* בעיית שלושת הגופים

==קישורים חיצוניים==

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

קטגוריה:בעיות חישוביות
קטגוריה:כבידה