|
: <math>{\partial u\over \partial t} =
ck \left({\partial^2 u\over \partial x^2 } +
{\partial^2 u\over \partial y^2 } +
{\partial^2 u\over \partial z^2 }\right)
<br /><br />
==פתרון משוואתכללי החום בממדלממד אחד ==
פתרון המשוואה, בממד אחד, על ידי [[הפרדת משתנים]] הוא:
:<math> \ u(t,x) = X(x) T(t). \quad</math>
<br />
:<math>\frac{T'(t)}{kT(t)} = \frac{X''(x)}{X(x)}. \quad </math>
<br />שני הצדדים של המשוואה הם משוואות התלויות במשתנים שונים, לכן הם חייבים להיות שווים לקבוע מספרי. הקבוע חייב להיות שלילי מכיוון שאחרת הטמפרטורה תגיע לאינסוף, ונסמנו λ²-. הפתרון הסופי המתקבל הוא:
:<math>T(t) = A e^{-\lambda^2 k t} \quad </math>
ו
:<math>X(x) = B \sin(\lambda \, x) + C \cos(\lambda \, x).</math>
<br />
אחרי שגילנו ש
<math>\ \frac{dT}{dt}= c\,{\delta}T</math>
שווה ל
<math>\ T(t) = e^{-ct}</math>
עכשיו אנחנו יכולים להתמודד עם הבעיה של גילוי הטמפ' בכל מקום ובכל זמן בתוך גוף.
לצורך העניין נסתכל על גוף שיש בו רק ציר אחד כלומר: צינור ארוך שמבודד לחלוטין מצדדיו.
<!-- עזרה בהוספת תמונה -->
קצוות הצינור מוחזקים בטמפ' קבועה של 0 מעלות צלזיוס. בתחילת הניסוי הטמפ' בכל נקודה בצינור מיוצגת ע"י הגרף שמצורף מעליו. כלומר פונקציה כלשהי של <math>\ T_0</math>
אנחנו מחפשים את הטמפרטורה בכל מקום בכל זמן. עלינו להגיע למשוואה דיפרנציאלית שתתאר את השינוי בטמפ'.כדי לעשות זאת נחלק את הצינור בדמיוננו לאינספור פרוסות דקות מאוד (ברוחב dX ) כך שבכל פרוסה הטמפ' יכולה להיחשב כאחידה.
<!-- עזרה בהוספת תמונה -->
את הפרמטר λ נקבל מתנאי השפה של הבעיה, והמשך הפתרון על ידי [[טור פורייה]].
נסתכל על מיקום מסוים X. הטמפ' בפרוסה זו היא <math>\ T(X)</math> והטמפ' הזו תשתנה בגלל שבסביבה ישנן שתי פרוסות שהטמפ' בהן שונה. אחת היא הפרוסה שמיקומה
<math>\ X+dX</math>
והטמפ' שלה:
<math>\ T(t+dx)</math>
והשניה שמיקומה
<math>\ x-dx</math>
והטמפ' בה <math>\ T(x-dx)</math>. קצב השינוי של הטמפ' תלוי בהפרשים של <math>\ T(x+dx)-T(x)</math> (הנגזרת) ו<math>\ T(x-dx)-T(x)</math>
הפרשים אלה קשורים לנגזרת של T לפי X <math>\ \frac{dT}{dx} = \ \frac{T(x+dx)-T(x)}{dx}</math>
ו<math>\ \frac{dT}{dx} = \ \frac{T(x-dx)-T(x)}{dx}</math>
למי שלא יודע זוהי הגדרת הנגזרת של T לפי X כאשר השתמשתי בסימן + כדי לייצג את הנגזרת בקצה הימני של הפרוסה ובסימן – כדי לייצג את הנגזרת בקצה השמאלי של הפרוסה.
==דוגמה לאילוץ תנאי התחלה ושפה==
כדי להבין זאת יותר טוב יש להסתכל על גרף הטמפ' לפי המקום: נבחר נקודה על הגרף. נסמן אותה X. הנגזרת שסימנו ב + מייצגת את שיפוע הגרף בכיוון ימינה מהנקודה וזו שסימנו ב – את שיפוע הגרף בכיוון שמאלה מהנקודה. שימו לב שהם הפוכים: אם כשהולכים ימינה הטמפ' יורדת אז כשהולכים שמאלה הטמפ' עולה. כלומר: הטמפ' בפרוסה x תגדל בגלל שבצד אחד שלה הטמפ' חמה יותר אבל תקטן בגלל שבצד השני הטמפ' קרה יותר. לכן השינוי בטמפ' בסופו של דבר יהיה תלוי בהפרש שבין הנגזרת מצד אחד לנגזרת מהצד השני:<math>\ \frac{dT}{dx} = \ \frac{d^2T}{dx^2}</math>
כלומר לקצב השינוי של הנגזרת כשזזים קצת ב X או במילים אחרות לנגזרת של הנגזרת.
קיבלנו אם כן שקצב השינוי של הטמפ' בכל נקודה פרופורציוני לנגזרת השניה של הטמפ' לפי המקום בנקודה זו.
מתמטית נייצג זאת כך:
ניקח מוט באורך L, המבודד כולו פרט לקצה אחד שלו, שם הוא מוחזק בטמפרטורה קבועה. התנאים שנקבל:
<math>\ \frac{dT}{dt} = \ c \frac{d^2T}{dx^2}</math>
זוהי משוואה דיפרנציאלית שכוללת נגזרות לפי שני משתנים שונים. כשגוזרים פונקציה של כמה משתנים כל פעם לפי חלק ממשתניה הנגזרות נקראות 'נגזרות חלקיות' ולכן משוואה כזו מכונה 'משוואה דיפרנציאלית חלקית'.
* '''תנאי התחלה:''' בזמן t=0 כל המוט בטמפ' החדר: <math> u(0,x) = T_{0} </math>
== פתרון לצינור באורך פאי ==
* '''תנאי שפה א':''' הטמפרטורה במקום x=0 קבועה תמיד: <math> u(t,0) = T_{max} </math>
* '''תנאי שפה ב':''' כל המוט מבודד, כך שנוכל לכתוב לגבי קצה המוט: <math> {\partial u\over \partial x}(t,L) = 0 </math>
מכיוון שבמשוואה הגדלים דיפרנציאליים, נוכל לבחור את נקודת האפס כרצוננו. נבחר את נקודת האפס של הטמפרטורה ב- <math> T_{max} </math>. למרות שהטמפרטורה במוט תמיד שלילית בסקלה זו, היא הטובה ביותר להתייחס בה לבעיה.
קודם עלינו לפתור משוואה של גוף חד ממדי לפני שנעבור לדו ממדי ותלת ממדי.
נאלץ את תנאי שפה א':
נתון צינור באורך פאי. הטמפ' בקצוות שלו היא 0 מעלות. ו <math>\ T_0 = 30 </math> מעלות.
<!-- עזרה בהוספת תמונה -->
עכשיו אנחנו צריכים לחפש איזושהיא פונקצייה שבנקודה-0 ובנקודה <math>\ \pi</math> שווה ל-0.(כמו נתוני הקצה בצינור...)
הפונקצייה אשר תתאר מצב זה היא פונקציית הסינוס. כיוון שבפאי וב0 היא מתאפסת.
חילקנו קודם את הצינור לdxים קטנים, שבכל אחד מהם הטמפ' כרגע היא 30 מעלות. אחנו צריכים לסכום את הטמפ' בכל הצינור. הפעולה שעושה זאת היא אינטגראל. את האינטגראל אנו עושים לפי אורך הצינור (מאפס עד פאי).
לפי טור פורייה, כל פונקצייה מחזורית ניתן להציג כסכום של פונקציות סינוס וקוסינוס שארכי הגל שלהם הם שברים שלמים של אורך הגל של הפונקצייה המקורית.
לכן משוואת החום החדשה היא:
<math>\ A_n\,e^{-cn^2t}\,Sin(nx)</math>
קודם כל, ה <math>\ A_n = \ A_1 + \ A_2 + \ A_3... \ lim_{A \to \infty}</math>
*כשA מייצג את מכפלת הסיונסים\קוסינוסים.
לכן תיאורטית לא ניתן להגיע לפונ' המדויקת אלא לטמפ' אשר שואפת לפונקצייה זו כמו סדרה מתכנסת.
<math>u(t,0) = A e^{-\lambda^2 k t} * C = 0 \quad </math>
כעת נעשה אינטרגראל מ0 עד פאי. הסיבה לכך היא שכאשר נציב בסינוס(ב-X) פאי או 0 יצא לנו כפולה של פאי או 0 שהסינוס שלהם הוא 0 ומקיים את תנאי הקצה.
ומכאן
נכפיל את הצינור בשתיים (אפשר לעשות זאת לפי משפט פורייה)
במקום זאת, פשוט נחלק את אורך הצינור ב2.
זאתי הנוסחא הכללית (כשאורך הצינור הוא פאי)
<math>\C = \frac{2}{\pi}\int_{0}^{\pi}\, \quad T_0\,Sin(nx)\,dx</math>
נאלץ את תנאי ב':
וכעת נחזור לתרגיל:
<math>\ \frac{2}{\pi}\int_{0}^{\pi}\, \ 30\,Sin(nx)\,dx</math>
בגלל ש30 קבוע, ניתן להוציאו מחוץ לאינטגראל ולכפול אותו בשתיים חלקי פאי.
<math>\ \frac{60}{\pi}\int_{0}^{\pi}\, \ Sin(nx)\,dx</math>
האינטגראל של Sin הוא Cos-.
<math>\ \frac{60}{\pi}\, \ \frac{-Cos(nx)}{n}</math>
כעת נציב פאי ואז ם (גבולות האינטגראל)
<math>\sum_{0}^{\pi} \frac{60}{\pi}\, \ \frac{-Cos(nx)}{n}</math>
<math> {\partial u\over \partial x}(t,L) = D e^{-\lambda^2 k t} \cos(\lambda \, L) = 0 </math>
מה שיוצא לאחר ההצבה הוא:
כאשר D כולל בתוכו מספר פרמטרים שהיו קודם.
<math>\ \frac{60}{\pi}\,[ \ \frac{-Cos({\pi}x)}{n}] - \ \frac{60}{\pi}\,[ \ \frac{-Cos(0)}{n}] </math>
למשוואה זו יש שתי פתרונות-כשN זוגי וכN אי זוגי.
כאשר N זוגי <math>\ N_{even} = 0</math> אם תציבו מספר זוגי בN תגלו זאת.
ונקבל כי λ יכול להיכתב בסדרה של ערכים אפשריים:
כאשר N אי זוגי <math>\ N_{odd} = \ \frac{60}{\pi}\, [\frac{-Cos{\pi}}{1} + \frac{cos0}{1}] = \frac{120}{\pi}</math>
<math>\lambda T(x\,t)_{n} = \frac{120}{2{\pi}} \,e^{-c4t \over L}\,Sin (2x)n + \frac{1201}{4{\pi2}}\,e^{-16t}\,Sin(4x)...\ \frac{120}{n{\pi}}\,e^{-cn^2t}\,Sin(nx)</math>
והרי לכם הפתרון של בעייה בממד אחד.
==כיצד יש לגשת לבעיה?==
א.נאפס תנאי שפה ואורך בצורה הבאה:
<math>\ Sin[\frac{n{\pi}x}{Length}]</math>
*הערה:רק כאשר אורך הצינור הוא אינו פאי או כפולה של פאי.
משתמשים בנוסחא זו בכדי לסדר את תנאי השפה-אם נציב X=0 ו X=length תצא כפולה של פאי אשר מתאפסת.(לפי מחזוריות Sin)
*הערה: כאשר משנים תנאי שפה במשוואה שמעלים את N בריבוע, מעלים את מה שהצבנו בSin לא כולל הX והN!
*הערה: כאשר באחד מן הקצוות הטמפ' אינה 0, באינטגראל יש להוריד את פונ' הצינור המתאימה, לפרק לשתי טורי פורייה, ולהוסיף את הפונ' לתוצאה.
כשיש מצב בו אחד הצדדים של הצינור מבודד-נכפיל את הצינור בשתיים, ואז תיווצר בעיה חדשה עם אורך <u>כפול</u> . וזו הופכת לבעיה רגילה.
כאשר שני הצדדים מבודדים נעשה אותו הליך של סידור תנאי שפה רק עם קוסינוס, כיוון שקוסינוס של פאי הוא 1- ושל שני פאי הוא 1. לכן בפתרון יצא או 1- או 1 כפול משוואת החום.
ב. עושים פוריה(עם ה An וכו'), עושים אינטגראל מ0 עד אורך הצינור, הצבה באינטגראל, הצבה של זוגי ואי זוגי ואורך הצינור אשר מתאפס בכל כפולה שלו. ומגלים את משוואת החום בכל מקום ובכל זמן.
מכאן שגם הפרמטר החופשי A יכול להיות מספר ערכים אפשריים, נקבל:
<math> u(t,x) = \sum_{n=0}^\infty A_{n} e^{-\lambda_{n}^2 k t} \ \sin(\lambda_{n} \, x) </math>
כשנשאר לנו למצוא את A<sub>n</sub>. נאלץ את תנאי השפה באמצעות [[טור פורייה]], וכפל בפונקציה <math>\sin(\lambda_{m} \, x)</math> עבור מספר שלם כלשהו m. אחרי ביצוע אינטגרל על כל המוט, נקבל:
<math> A_{n} = \frac{4 \ (T_{0}-T_{max})}{\pi \, (2n + 1)} </math>
כך שהפתרון למקרה אחרי אילוץ כל התנאים הוא:
<math> u(t,x) = T_{max} - \frac{4 \ (T_{max} - T_{0})}{\pi} \ \sum_{n=0}^\infty \frac{1}{2n+1} \ \sin \left( \pi \ (n+\frac{1}{2}) \ \frac{x}{L} \right) \ \exp \left( -\frac{k \ \pi^{2}}{L^{2}} \ (n+\frac{1}{2})^{2} \ t \right) </math>
[[קטגוריה:תרמודינמיקה]]
[[sv:Värmeledningsekvationen]]
[[zh:熱傳導方程式]]
כל הזכויות שמורות לעודדי, יותם ודוראל.
| 