משוואה דיפרנציאלית חלקית

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

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

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

רקע מתמטי[עריכת קוד מקור | עריכה]

סימונים[עריכת קוד מקור | עריכה]

נגזרת חלקית של הפונקציה \ u לפי המשתנה \ x, תסומן כ-\ u_x או \ \frac{\partial u}{\partial x}.

באופן דומה נגזרת חלקית לפי \ x ואחריו לפי \ y, תסומן כ-\ u_{xy} או \ \frac{\partial^2 u}{\partial y \partial x} וכן הלאה.

תנאי שפה[עריכת קוד מקור | עריכה]

Postscript-viewer-shaded.png ערך מורחב – תנאי שפה

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

תנאי השפה מחולקים לשני סוגים עיקריים:

  • תנאי שפה מסוג דיריכלה (Dirichlet boundary condition), בהם נתון ערך הפונקציה בגבולות הבעיה (בזמן מסוים ובגבולות של המרחב).
  • תנאי שפה מסוג נוימן (Neumann boundary condition), בהם נתונה הנגזרת של הפונקציה בגבולות הבעיה (נגזרת לפי הכוון הניצב לשפה).

בנוסף ניתן לצרף אותם לסוגים חדשים של תנאי שפה:

  • תנאי שפה מסוג קושי (Cauchy boundary condition), בהם נתונים ערכי הפונקציה וערכי הנגזרת שלה, כלומר עירוב תנאי דיריכלה ונוימן.
  • תנאי שפה מסוג רובין (Robin boundary condition), בהם נתון ערך עבור סכום של הפונקציה והנגזרת שלה באותה משוואה.

פתרונות של משוואות דיפרנציאליות חלקיות[עריכת קוד מקור | עריכה]

קיום פתרון ויחידות הפתרון[עריכת קוד מקור | עריכה]

Postscript-viewer-shaded.png ערך מורחב – משפט הקיום והיחידות (משוואות דיפרנציאליות)

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

בניסוח אחר: הפתרון יציב כאשר הוא תלוי באופן רציף בתנאים הנלווים לבעיה, כלומר ההפרש בין שני פתרונות כאלה יתכנס במידה שווה ל-0, כאשר השינוי בתנאי ההתחלה שואף ל-0.

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

פתרונות אנליטיים[עריכת קוד מקור | עריכה]

עבור סוגים מסוימים של משוואות דיפרנציאליות חלקיות ניתן לפתור באופן מתימטי מלא את המשוואה ולמצוא לה פתרון אנליטי (כלומר למצוא צורה פונקציונלית של הפתרון, \ u=f(x_1, x_2,..., x_n)). שיטות עיקריות לפתרון אנליטי הן:

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

פתרונות נומריים[עריכת קוד מקור | עריכה]

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

תיאור תופעות פיזיקליות באמצעות משוואות דיפרנציאליות חלקיות[עריכת קוד מקור | עריכה]

תופעות פיזיקליות רבות מתוארות על ידי משוואות דיפרנציאליות חלקיות. הדוגמאות הנפוצות ביותר הן:

  • משוואת הגלים (בתלת-מימד): \ { \partial^2 \psi \over \partial t^2} = v^2 \nabla^2 \psi כאשר \ \psi הוא מרחק ההפרעה מנקודת שיווי המשקל במקום (x,y,z) בזמן t, ו-v היא מהירות התקדמות הגל, ו  \nabla^2 הוא אופרטור הלפלסיאן.
  • משוואת החום: \ u_t - k \nabla^2 u = 0 כאשר \ u(x,y,z,t) היא הטמפרטורה בנקודה כלשהי בזמן t, ו-k הוא מקדם הולכת החום. נקראת גם משוואת הדיפוזיה.
  • גלי הלם (משוואת אוילר): \ u_t+uu_x=0 כאשר u מתארת התנהגות ריכוז המסה בזמן, והמהירות תלויה בריכוז המסה.
  • פוטנציאל חשמלי בריק (משוואת לפלס):\ \nabla^2 \varphi = 0 כאשר \ \varphi מתאר פוטנציאל חשמלי באזור ללא מטענים. פונקציה המקיימת את משוואת לפלס נקראת פונקציה הרמונית.
  • משוואת שרדינגר: \ i \psi_t + \frac{\hbar^2}{2m}\nabla^2\psi + V \psi = 0 כאשר \ \psi (x,y,z,t) היא פונקציית הגל של חלקיק קוונטי וV היא פונקציית הפוטנציאל של החלקיק.

משוואות לינאריות מסדר שני[עריכת קוד מקור | עריכה]

תופעות רבות בטבע ניתנות לתאור על ידי משוואות לינאריות עם נגזרות עד סדר שני, בהצגת המשוואות באופן הבא: \ A(x,y)u_{xx}+B(x,y)u_{xy}+C(x,y)u_{yy}+D(x,y)u_x+E(x,y)u_y+Fu=G(x,y). ניתן למיין משוואות אלו לפי ערך הביטוי \ B^2-4AC: המשוואה היא מסוג היפרבולי כאשר הביטוי חיובי, מסוג פרבולי כאשר הביטוי מתאפס, ומסוג אליפטי כאשר הביטוי שלילי. סיבת חלוקה שכזו מובנת כאשר \ A, B, C, D, E, F הם קבועים, ו- \ G=0: במקרה זה תמיד ישנו הפתרון \ u=e^{mx+ny}, וכאשר נציבו במשוואה, נקבל את הקשר \ Am^2+Bmn+Cn^2+Dm+En+F=0, ומגאומטריה אנליטית ידוע שזהו למעשה משטח קוני במשתנים \ m,n אשר סוג חתכיו (הפרבולי, פרבולי, אליפטי) נקבעים על ידי הביטוי \ B^2-4AC.

משוואות היפרבוליות[עריכת קוד מקור | עריכה]

משוואות פרבוליות[עריכת קוד מקור | עריכה]

משוואות אליפטיות[עריכת קוד מקור | עריכה]

פתרונות אמיתיים ופתרונות מוכללים[עריכת קוד מקור | עריכה]

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

ראו גם[עריכת קוד מקור | עריכה]

לקריאה נוספת[עריכת קוד מקור | עריכה]

  • יהודה פינצ'ובר ויעקב רובינשטיין, מבוא למשוואות דיפרנציאליות חלקיות, מהדורה שלישית, הטכניון, 2006.