משוואות אוילר (מכניקת הזורמים)

מתוך המכלול, האנציקלופדיה היהודית
קפיצה לניווט קפיצה לחיפוש
נוזל זורם דרך קובייה דמיונית. פריצת הדרך בתחום מכניקת הזורמים התרחשה בעקבות ההבנה שכדי להבין את התנהגות הזורם ניתן לחלק את הנוזל לתאים קטנים סטטיים, ולבחון את השתנות הזורם בזמן בכל תא בנפרד. כך למשל, על ידי הכלת עקרון שימור המסה על מסת הנוזל הנכנסת ויוצאת מהקוביה בתמונה, ניתן להסיק את משוואת הרצף, שהיא אחת ממשוואות אוילר.
Nuvola apps edu mathematics blue-p.svg

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


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

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

משוואות אוילר

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

המשוואות בממד אחד ומשמעותן הפיזיקלית

  1. משוואת הרצף (שימור מסה):
  2. משוואת התנועה (שימור תנע):
  3. משוואת שימור האנרגיה:

כאשר:

  • היא צפיפות הנוזל
  • היא מהירות הנוזל
  • הוא הלחץ בנוזל
  • הוא סך כל האנרגיה ליחידת נפח, כלומר , כאשר היא האנרגיה הפנימית של הנוזל.

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

משוואת הרצף

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

נראה שקילות בין המשוואה למשמעות המילולית שניתנה לעיל. נתבונן בתא סטטי בזורם בעל נפח . הביטוי משמעו קצב השינוי של הצפיפות בנפח . הביטוי הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle {\frac {\partial \left(\rho v\right)}{\partial x}}} מבטא את הקצב שבו מסה נכנסת או יוצאת מאזור זה, כאשר הביטוי מייצג את זרם המסה. המשמעות הפיזיקלית של זרם המסה דומה למשמעות של זרם חשמלי, אך עם תנועה של מסה במקום תנועת מטענים חשמליים. מכאן שהמשמעות של המשוואה היא שעל מנת למצוא את השינוי של המסה בנפח , הוא , יש להתבונן בהפרש בין זרם המסה הנכנס לתא לזרם המסה היוצא ממנו. משמעות עובדה זו היא שימור מסה, שכן היא קובעת שלא מופיע נוזל יש מאין.

משוואת התנועה

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

החוק השני של ניוטון קובע כי , והוא אכן מבטא שימור תנע. לעיתים קרובות קשה להגדיר נפח במכניקת זורמים, ועל כן לרוב עובדים עם גדלים שהם "ליחידת נפח". כך למשל נהוג לעבוד עם צפיפות במקום עם מסה. מסיבה זו, נבטא את החוק השני של ניוטון בצורה: , כאשר הוא כוח ליחידת נפח. כעת, את התאוצה בחוק ניוטון אפשר להחליף בהגדרתה - שינוי המהירות כתלות בזמן, . בנוסף, מכיוון שלחץ הוא כוח ליחידת שטח, אפשר לכתוב הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle f=-{\frac {\partial P}{\partial x}}} . מכאן מקבלים הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle -{\frac {\partial P}{\partial x}}=\rho {\frac {\partial v}{\partial t}}} , ועל ידי סידור המשוואה מתקבלת המשוואה הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle {\frac {\partial v}{\partial t}}=-{\frac {1}{\rho }}{\frac {\partial P}{\partial x}}} . משוואה זו דומה למשוואת שימור התנע שהוצגה, אך יש להסביר את האיבר האחרון שנותר - . איבר זה הוא איבר ה"הסעה", וקשור לעובדה שבחלקים שונים בנוזל ייתכנו מהירויות שונות, מה שיגרום לתאוצה. דוגמה לכך היא תנועה של נוזל דרך צינור שהולך וצר. במקרה זה, ככל שהצינור נהיה צר יותר, כך מהירות הנוזל היא גדולה יותר, ואנו מקבלים תאוצה. תופעה זו מכונה האצת הסעה (Convective acceleration).

משוואת שימור האנרגיה

משוואה זו מתארת שימור אנרגיה. בדומה למשוואת הרצף, משמעות המשוואה היא שהשינוי באנרגיה ליחידת נפח באזור בנוזל, הוא הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle {\frac {\partial E}{\partial t}}} , הוא תוצאה של זרימת אנרגיה לתוך נפח זה, המתואר על ידי זרם האנרגיה - . למעשה, ניתן להבין משוואה זו כמעין משוואת רצף עבור אנרגיה. נדגיש כאן שמכיוון שמשוואות אוילר מתארות זורמים, בחישוב האנרגיה עלינו להתחשב באנרגיה הפנימית של הנוזל, ולא באנרגיה הקינטית שלו בלבד.

המשוואות במימד אחד בזרימה רדיאלית

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

המשוואות בתלת ממד

היסטוריה

שורשן של משוואות אוילר נובע מחוקי ניוטון שנוסחו בשנת 1687[1]. בעוד חוקי ניוטון מתארים נכונה את החוקים הפיזיקליים הכלליים הפועלים על זורמים, יישום שלהם על נוזל מסוים מוביל לקשיים - שכן לא ניתן לראות נוזל רק כסך המולקולות המרכיבות אותו. זאת מכיוון שבפתרון המשוואות עבור נוזלים, רצוי להתייחס למושגים מאקרוסקופיים המתארים אותם כלחץ ונפח, שלרוב לא ניתן להגדיר היטב על מולקולות בודדות. רק כאשר מתייחסים לנוזל כאל רצף, ומשיימים את עקרונות החשבון האינפיניטסימלי עליו בכך שמחלקים אותו לגדלים קטנים ומבצעים חישוב על כל אחד בנפרד, ניתן ליישם כהלכה את חוקי ניוטון. ההבנה הזו הייתה תהליך ממושך, שהתחיל עם הניסיון הראשון של דניאל ברנולי ב־1738 ליישם עקרונות דינמיים כלליים על נוזל בתנועה, והגיעה לשיאה 17 שנה לאחר מכן עם פרסום משוואות אוילר על ידי לאונרד אוילר. המחקר שנעשה בשנים אלו לגילוי החוקים הפועלים על זורמים הוביל בין היתר לפיתוח של כלים מתמטיים חדשים, כמו כתיבה של משוואות וקטוריות, המושג של שדה מהירות וטכניקות לפתרון משוואות דיפרנציאליות חלקיות[2].

ב־1738 כתב דניאל ברנולי בספר Hydrodynamica, sive de viribus et motibus fluidorum commentarii (הידרודינמיקה, דיסרטציה על כוחות זורמים ותנועתם) את הניסיון הראשון ליישום של עקרון כללי לתנועת זורמים. ברנולי יישם את עקרון שימור האנרגיה של הויגנס על זורמים והשיג תיאור גאומטרי של חוקים הכוללים אף גרסה חד־ממדית של משוואות אוילר.

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

ב-1746 זכה המתמטיקאי ז'אן לה רון ד'אלמבר בפרס מטעם האקדמיה הפרוסית למדעים בברלין על מאמרו The Cause of winds. במקרה פרטי במאמר, ניתנות לראשונה משוואות של תנועת זורמים אי דחיסים במקרה הדו־ממדי, שהיא מקרה פרטי של משוואות אוילר. שנתיים לאחר מכן, פרסם ד'אלמבר מאמר בשם The resistance of fluid לתחרות נוספת באקדמיה הפרוסית אך לא זכה בפרס. למרות זאת, חשיבותו של המאמר אינה מוטלת בספק - הוא מכיל לראשונה את מה שכונה "משוואת הערבוליות לזורמים אי־דחיסים", תחת ההנחה שתמיד מתקיים: . הנחה זו פירושה שהערבוליות אפסית, הנחה מוטעית שאוילר עצמו ישתמש בה במאמרו שלוש שנים מאוחר יותר.

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

אך ניכר שאוילר לא היה מרוצה ממונוגרפיה זו[3]

ב־1752 כתב אוילר את המאמר Principia motus fluidorum (עברית: עקרונות התנועה של נוזלים) שהציג גרסה מוקדמת של המשוואות לזרימה אי דחיסה. ב־31 באוגוסט 1752 המאמר הוגש לאקדמיה הפרוסית למדעים בברלין[4]. המאמר כלל את הטעות שנעשתה גם על ידי המתמטיקאי ד'אלמבר בבואו למצוא חוקים לזורמים, והיא ההנחה שהערבוליות אפסית. במאמר זה פיתח אוילר לא רק גרסה דו־ממדית של המשוואות, אלא גם לראשונה גרסה תלת־ממדית. אוילר הבין אחרי זמן מה את הטעות שעשה בכך שהתעלם מפתרונות ערבוליים. בשנת 1755 כתב אוילר שלושה מונוגרפים על זורמים.

העמוד הראשון במונוגרפיה הראשונה של אוילר על מכניקת זורמים - Principes generaux de l'etat d'equilibre des fluides בהוצאת האקדמיה הפרוסית למדעים בברלין, דו"ח שנתי לשנת 1757

המונוגרפיה הראשונה, Principes generaux de l'etat d'equilibre des fluides (עברית: עקרונות כלליים בנוגע לשיווי משקל בזורמים), הוגשה לאקדמיה בברלין ב־11 באוקטובר 1753[5] והוקדשה לשיווי־משקל של זורמים, אי־דחיסים ודחיסים. במונוגרפיה זו ניסה אוילר לבסס את מדע ההידרוסטטיקה על עקרון מנחה יסודי, והוא שכל אלמנט אינפיניטסימלי בנוזל אשר נמצא בשיווי משקל יקיים את המשוואה , כאשר P הוא הלחץ בשפת האלמנט ו־f היא צפיפות הכוח הפעולת עליו. אוילר הראה כי התוצאות הידועות של הידרוסטטיקה נובעות מחוק מתמטי פשוט זה. ואכן, בתחילת המונוגרפיה נכתב:

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

המונוגרפיה השנייה, Principes generaux du mouvement des fluides (עברית: עקרונות כללים בנוגע לתנועה של זורמים), הוגשה לאקדמיה בברלין ב־4 בספטמבר 1755[7], והוקדשה ליישום העקרונות מהמונוגרפיה הקודמת לטובת הבנת תנועה של זורמים. בתחילת המונוגרפיה נכתב שהבנת תנועה של זורם תלויה בתנאי התחלה של הזורם, שמכונים "המצב הפרימיטיבי של הזורם", אחד האזכורים הראשונים לכך בספרות המתמטית. ממצב זה יש להבין את הכוחות הפועלים על הזורם וכך ניתן להסיק את מצבו הנוכחי. לאורך המונוגרפיה משתמש אוילר בעיקר בצורה התלת־ממדית של המשוואות, ומפתח אותן גם לצורות כלליות כגון עבור זורמים דחיסים:

הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle {\begin{aligned}\nabla \cdot \left(\rho {\boldsymbol {v}}\right)+{\frac {\partial p}{\partial t}}=0\end{aligned}}}

כך שההנחה ש־ אינה הכרחית.

המונוגרפיה השלישית, Continuation des recherches sur la theorie du mouvement des fluides (עברית: המשך המחקר בנוגע לתנועה של זורמים), הוגשה לאקדמיה בברלין ב־2 באוקטובר 1755[8]. המונוגרפיה היא המשך ישיר של השתיים הקודמות ופיתוח התוצאות בהן, ובתחילתה כתב אוילר על חשיבות תגליותיו:

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

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

שלושת המונוגרפיות פורסמו על ידי האקדמיה הפרוסית בברלין בשנת 1757 באוגדן Mémoires de l'académie des sciences de Berlin 11. כתבים אלו היוו פריצת דרך בהבנת מכניקת זורמים. בסוף המונוגרפיה השנייה שלו משנת 1755 תיאר אוילר את שנותר לגלות בתחום מכניקת הזורמים לאחר פיתוח משוואותיו:

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

לא כל הקהילה המדעית העריכה את אוילר כראוי על גילוי המשוואות הנושאות את שמו. בספר Méchanique analitique של ז'וזף לואי לגראנז' משנת 1788 מצוין דווקא ד'אלמבר כמייסד תחום מכניקת הזורמים, בעוד שמו של אוילר לא מצוין לכל אורך הספר. במהדורה השנייה של הספר שיצאה יותר מ־20 שנה לאחר מכן מופיע המשפט:

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

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

פיתוח של משוואות אוילר

זורם חד־ממדי, הנע ימינה במהירות v, עם הפרמטרים המשומשים בפיתוח

מספר הערות על פיתוח המשוואות:

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

משוואת הרצף

בהינתן צינור ארוך בעל שטח חתך אשר זורם בו נוזל בכיוון החיובי של ציר ה־, החידוש של אוילר היה בחלוקת הצינור לתאים בגודל אינפיניטסימלי ושימוש בכלים של החשבון אינפיניטסימלי כדי לחשב את התנועה של כל החומר הזורם באמצעות זרימה נקודתית. על כן, נחלק את הצינור לתאים אינפיניטסימליים באורך ונתבונן בתא מסוים במיקום . נסמן את הצפיפות והמהירות של הנוזל בו ב־, הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle v_{0}\left(t\right)} בהתאמה. בדומה, נסמן את הצפיפות והמהירות של הנוזל בתא משמאלו ב־ , בהתאמה (ראה איור משמאל).

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

ראשית, כתוצאה מהקידום בזמן הנוזל בתא השמאלי יתקדם מרחק של הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle v_{-}\Delta t} , ולכן סך כל המסה הנכנסת לתא ב־ היא . באותו האופן, סך כל המסה היוצאת מהתא ב־ לתא מימנו היא . מכאן שהשינוי במסה בתא ב־ הוא ההפרש בין שני הביטויים. כעת, את השינוי במסה נוכל לרשום גם במפורש בתור ההפרש , כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle m\left(t\right)=\rho_{0}\left(t\right)\cdot A\Delta x} היא המסה בתא ב־ בזמן הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle t} .

נשווה בין שני הביטויים שמצאנו ונקבל:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \rho_{0}\left(t+\Delta t\right)\cdot A\Delta x-\rho_{0}\left(t\right)\cdot A\Delta x=\rho_{-}\left(t\right)\cdot Av_{-}\Delta t-\rho_{0}\left(t\right)\cdot Av_{0}\Delta t}

נחלק ב־ וב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta t} ונקבל:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\rho_{0}\left(t+\Delta t\right)-\rho_{0}\left(t\right)}{\Delta t}=\frac{\rho_{-}\left(t\right)v_{-}-\rho_{0}\left(t\right)v_{0}}{\Delta x}}

ועל ידי לקיחת הגבולות הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \Delta x,\Delta t\rightarrow 0} נקבל מהגדרת הנגזרת את משוואת הרצף:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\partial\rho}{\partial t}=-\frac{\partial\left(\rho v\right)}{\partial x}}

ומהפיתוח מסיקים שמשוואה זו אכן מבטאת שימור מסה.

משוואת שימור תנע

כדי לפתח את משוואת שימור התנע משתמשים ברעיונות וכלים דומים. נקדם את המערכת בזמן הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta t} , ונחפש שני ביטויים לשינוי בתנע של התא ב־.

ראשית, כמות התנע הנכנס לתא האינפיניטסימלי ב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_{0}} מהתא משמאלו הוא מכפלת מסת החומר הנכנס במהירותו, כלומר: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left(Av_{-}\Delta t\cdot\rho_{-}\right)v_{-}} . בדומה, כמות התנע העוזבת את התא אל התא מימינו היא הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \left(Av_{0}\Delta t\cdot \rho _{0}\right)v_{0}} . השינוי בתנע הוא ההפרש ביניהם.

מצד שני, אפשר לכתוב במפורש את השינוי בתנע בתור ההפרש הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta p=p\left(t+\Delta t\right)-p\left(t\right)} , כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p\left(t\right)=A\Delta x\rho\left(t\right)\cdot v\left(t\right)} הוא התנע בתא בזמן t.

בנוסף על כך, בעקבות הקידום בזמן פועל מתקף על הנוזל. נסמן ב־ את הלחץ. התא השמאלי מפעיל מתקף הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle P_{-}A\Delta t} על החומר הנכנס לתא, והתא ב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_{0}} מפעיל מתקף על החומר היוצא מהתא. מכאן שהמתקף הכולל שפעל על החומר בתא ב־ הוא הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left(P_{1}-P_{0}\right)A\Delta t} .

על ידי השוואה בין שני הביטויים לשינו בתנע, מתקבלת המשוואה:

הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle A\Delta x\rho \left(t+\Delta t\right)\cdot v\left(t+\Delta t\right)-A\Delta x\rho \left(t\right)\cdot v\left(t\right)=\left(Av_{-}\Delta t\cdot \rho _{-}\right)v_{-}-\left(Av_{0}\Delta t\cdot \rho _{0}\right)v_{0}+\left(P_{1}-P_{0}\right)A\Delta t}

נחלק ב־ וב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta t} ונקבל:

על ידי לקיחת הגבולות הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta x,\Delta t\rightarrow0} נקבל מהגדרת הנגזרת את משוואת התנועה:

על מנת לקבל את המשוואה בצורתה המוכרת, נשתמש בעובדה כי ונוכל לכתוב:

כעת, נשתמש במשוואת הרצף על מנת לכתוב הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\partial\rho}{\partial t}+\frac{\partial\left(\rho v\right)}{\partial x}=0} , ונקבל:

ומהפיתוח מסיקים כי משוואה זו אכן מבטאת שימור תנע.

משוואת שימור האנרגיה

פיתוח משוואת שימור האנרגיה נעשה אף הוא בכלים וברעיונות דומים. נציג את השינוי באנרגיה, בתא ב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_{0}} בשתי דרכים שונות. הדרך הראשונה היא הפשוטה יותר, ובה נתאר את האנרגיה בתא לפני ואחרי ההתקדמות בזמן בצורה מפורשת: הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \Delta E=E\left(t+\Delta t\right)-E\left(t\right)} , כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle E\left(t\right)=A\Delta x\rho\left(t\right)\left(\frac{1}{2}v\left(t\right)^{2}+\epsilon\left(t\right)\right)} מבטאת את סך כל האנרגיה בתא בזמן t.

הדרך השנייה מתחשבת רעיונית באנרגיה שעוזבת ונכנסת לתא בעקבות הקידום בזמן. האנרגיה שנכנסת לתא היא האנרגיה של החומר שנכנס לתא מהתא ב־. סך כל החומר הנכנס לתא הוא בעל מסה של , ולכן האנרגיה שלו היא . בדומה, האנרגיה שעוזבת את התא היא הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle Au_{0}\Delta t\cdot\rho_{0}\cdot\left(\frac{1}{2}v_{0}{}^{2}+\epsilon_{0}\right)} .

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

  • אגף שמאל - הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle A\Delta x\rho\left(t+\Delta t\right)\left(\frac{1}{2}v\left(t+\Delta t\right)^{2}+\epsilon\left(t+\Delta t\right)\right)-A\Delta x\rho\left(t\right)\left(\frac{1}{2}v\left(t\right)^{2}+\epsilon\left(t\right)\right)}
  • אגף ימין -

נחלק ב־ וב־ ונקבל:

  • אגף שמאל - הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\rho\left(t+\Delta t\right)\left(\frac{1}{2}v\left(t+\Delta t\right)^{2}+\epsilon\left(t+\Delta t\right)\right)-\rho\left(t\right)\left(\frac{1}{2}v\left(t\right)^{2}+\epsilon\left(t\right)\right)}{\Delta t}}
  • אגף ימין -

ועל ידי לקיחת הגבולות נקבל מהגדרת הנגזרת:

נסמן ונקבל

מסקנות ממשוואות אוילר

משוואת ברנולי

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

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

נצא ממשוואת שימור האנרגיה:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\partial E}{\partial t}+\nabla\cdot\left(\overline{v}\left(E+p\right)\right)=0 }

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

הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \nabla \cdot \left({\overline {v}}\left(E+p\right)\right)=0\Longrightarrow \nabla \cdot \left({\overline {v}}\left({\frac {1}{2}}\rho v^{2}+\rho \epsilon +p\right)\right)=0}

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

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 0=\nabla\cdot\left(\overline{v}\rho\left(\frac{1}{2}v^{2}+\epsilon+\frac{p}{\rho}\right)\right)=\left(\frac{1}{2}v^{2}+\epsilon+\frac{p}{\rho}\right)\nabla\cdot\left(\overline{v}\rho\right)+\rho\overline{v}\cdot\nabla\left(\frac{1}{2}v^{2}+\epsilon+\frac{p}{\rho}\right) }

ממשוואת הרצף מסיקים כי הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle \nabla \cdot \left({\overline {v}}\rho \right)={\frac {\partial \rho }{\partial t}}} . בנוסף, הזרימה היא עמידה, ולכן . מכאן מתקבל:

על ידי אינטגרציה על קו מסלול (המשיק לוקטור המהירות בכל נקודה) מסיקים כי

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{1}{2}v^{2}+\epsilon+\frac{p}{\rho}=const }

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

משוואת הגלים ומהירות הקול

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

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

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

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} \rho\left(x,t\right)=\rho_{0}+\Delta\rho\left(x,t\right) ,\;\;\;\; P\left(x,t\right)=P_{0}+\Delta P\left(x,t\right) ,\;\;\;\; v\left(x,t\right)=\Delta v\left(x,t\right) \end{align} }

כאשר הם קבועים, והביטויים הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta v,\Delta P,\Delta\rho} מייצגים את ההפרעה.

נחקור את התקדמות ההפרעה עד לסדר ראשון בה, שכן לפי ההנחה מדובר בהפרעה קטנה, ולכן נזניח איברים מהצורה ואיברים מהצורה (וכאן הנחנו שגם הנגזרת של ההפרעה היא קטנה). תהליך הלינאריזציה הוא הצבת הביטויים ל־ שמופיעים למעלה במשוואות אוילר, ואז הזנחת איברים מסדר שני ב־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta} . כתוצאה מהלינאריזציה, משוואת הרצף ומשוואת שימור התנע מקבלות את הצורה הבאה:

  1. משוואת הרצף: .
  2. משוואת שימור התנע: .

כעת, נניח שהתהליך הוא אדיאבטי. הנחה זו מקובלת בפיזיקה כאשר מתכוונים לומר שהתהליך הוא "איטי". לכן ניתן לכתוב , כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle S} היא האנטרופיה. נסמן , ונראה שזוהי מהירות הקול בנוזל.

נגזור את משוואת הרצף בזמן: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{\partial^{2}\left(\Delta\rho\right)}{\partial t^{2}}+\rho_{0}\frac{\partial^{2}\left(\Delta u\right)}{\partial x\partial t}=0} .

כעת נגזור את משוואת התנועה במיקום: . נציב את הקשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta P=c_{s}^{2}\Delta\rho} , והמשוואה מקבלת את הצורה: .

על ידי חיסור שתי המשוואות מתקבלת המשוואה:

זוהי משוואת גלים עבור ההפרעה . בפרט, מסיקים כי מהירות ההפרעה היא הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle c_{s}=\left(\frac{\partial P}{\partial\rho}\right)_{S}} , ולכן זוהי מהירות הקול בנוזל.

חישוב דומה מראה שגם כאשר הנוזל אינו במנוחה אלא נע, מהירות הקול בו אינה משתנה[10].

גלי הלם

Postscript-viewer-shaded.png ערך מורחב – גל הלם

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

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

שימושים באווירודינמיקה

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

נניח שברצוננו לבנות מודל של כנף של מטוס עבור מנהרת רוח, כך שכל חלק באורך בכנף מתורגם לחלק בגודל הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x\prime=\alpha x} במודל. כך למשל אם , המודל יהיה בדיוק באורך ובגבוה שהם מאורך וגובה הכנף, בהתאמה. מגדיר את סקאלת האורך במערכת מנהרת הרוח. מכאן שגם סקאלת המהירות תשתנה, ומתקיים:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} v\prime=\frac{dx\prime}{dt}=\alpha\frac{dx}{dt}=\alpha v \end{align} }

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

ומשוואת התנועה עבור המודל (בסקאלת האורך של המודל) היא

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} \frac{\partial v\prime}{\partial t}+v\prime\frac{\partial v\prime}{\partial x\prime}=-\frac{1}{\rho}\frac{\partial P}{\partial x\prime} \end{align} }

כאשר מציבים את הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle x\prime ,v\prime } מתברר שמשוואת התנועה בסקאלה של המודל היא שונה מזו של הכנף. על מנת ששתי המשוואות יהיו זהות, יש לשנות גם את סקאלת הזמן במודל שלנו, ולהגדיר . המשמעות של השוויון היא שאם המודל מצליח לרחף במנהרת הרוח במשך שנייה, נסיק כי הכנף תצליח לרחף במשך שניות. כעת, משוואת התנועה בסקאלה של המודל היא:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} \frac{\partial v\prime}{\partial t\prime}+v\prime\frac{\partial v\prime}{\partial x\prime}=-\frac{1}{\rho}\frac{\partial P}{\partial x\prime} \end{align} }

נותר להראות שהמשוואה האחרונה זהה למשוואת התנועה עבור הכנף (כלומר, למשוואת התנועה בסקאלת הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle x,t} ). הפעם מתקיים

ולכן:

ומכאן נובע:

וכעת ניתן לראות שמשוואת התנועה בסקאלה של זהה לזו בסקאלה של הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x,t} , ולכן משוואות התנועה של הכנף והמודל יהיו זהות. מכאן שעל ידי מעברי הסקאלה הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle x\prime =\alpha x,\;\;t\prime =\alpha t} ניתן להשתמש במודל על מנת להבין את תנועת הכנף מבלי הצורך לבנות מנהרות רוח מספיק גדולות עבור הכנף בגודלה האמיתי.

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

הכללה מעבר לנוזלים קלאסיים

על נוזלים ונוזלים קוונטיים

אחת התכונות המרכזיות של על־נוזל הוא היותו חסר צמיגות. תוצאה מעניינת של תכונה זו היא היכולת של על־נוזל "לטפס" על דפנות של כלים.

בעזרת התאמות מסוימות לשפה של תורת הקוונטים, ניתן לפתח סט של משוואות אוילר עבור על נוזלים. על־נוזל הוא פאזה של חומרים מסוימים הבאה לידי ביטוי לרוב בטמפרטורות נמוכות, ומאופיינת בתכונות יוצאות דופן, כגון חוסר בצמיגות. על-נוזל מאופיין בכך שכמות מאקרוסקופית של חלקיקים מאכלסים את רמת היסוד. כתוצאה מכך, ניתן להגדיר פונקציית גל המתארת את כל הנוזל, ולא רק חלקיק בודד (כפי שנהוג לרוב בבעיות בתורת הקוונטים). פונקציית הגל היא מהצורה הבאה: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\psi\left(\overline{r},t\right)=\psi_{0}\left(\overline{r},t\right)e^{iS\left(\overline{r},t\right)}} . במקרה זה, ניתנת משמעות מעט שונה לפונקציית הגל:

  • הערך המוחלט יזוהה עם צפיפות החלקיקים הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle n_{s}=\frac{\rho_{s}}{m}} של העל־נוזל. זוהי אנלוגיה למקרה של חלקיק יחיד, שבו מזוהה עם צפיפות ההסתברות להימצאות החלקיק במקום מסוים.
  • הנגזרת של הפאזה הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle \textstyle arg\left(\psi \right)=S\left({\overline {r}},t\right)} תזוהה בתור מהירות העל־נוזל, עד כדי קבוע: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\nabla S=\frac{m}{\hbar}v_{s}} . זיהוי זה מקובל לעיתים גם במקרה של פונקציית גל של חלקיק יחיד.

מכאן, ניתן לצאת ממשוואת שרדינגר לפונקציות גל ולפתח שלוש משוואות האנלוגיות למשוואות אוילר[11]:

  1. משוואת הרצף: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\frac{\partial\left|\psi\right|^{2}}{\partial t}=-\nabla\overline{J}_{p}} , כאשר הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle {\overline {J}}_{p}} הוא זרם ההסתברות: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\overline{J}=-\frac{i\hbar}{2m}\left(\psi*\nabla\psi-\psi\nabla\psi*\right)} . משוואה זו היא משוואה מוכרת בתורת הקוונטים, ונהוג להבינה בתור משוואה המתארת שימור הסתברות[12].
  2. משוואת התנועה: , כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_{s}} הוא מהירות העל־נוזל ו־ הוא הפוטנציאל הכימי של העל־נוזל. משוואה זו מתקבלת תחת ההנחה שהצפיפות של העל-נוזל משתנה לאט במרחב. בפרט, מכיוון שמשוואת התנועה המתארת את העל־נוזל היא משוואת התנועה של אוילר ולא משוואת נאווייה-סטוקס, ניתן להסיק שהעל־נוזל הוא חסר צמיגות, שכן משוואות אוילר מתארות נוזל שהוא בפרט חסר צמיגות. זוהי תכונה חשובה של העל־נוזל.
  3. משוואת שימור האנרגיה: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\frac{1}{2}v_{s}^{2}+\frac{P}{\rho}-\sigma T=0} , זוהי משוואה הדומה למשוואת ברנולי, כאשר הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle v_{s}} היא מהירות העל־נוזל, הוא הלחץ, הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle T} הוא הטמפרטורה ו־הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \sigma} היא האנטרופיה ליחידת מסה של העל־נוזל.

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

צמיגות ומשוואות נאוויה-סטוקס

השימוש במשוואות אוילר אפשרי רק תחת ההנחה שהנוזל חסר צמיגות. הנחה זו אינה תמיד נכונה במציאות. מכיוון שצמיגות בנוזלים אנלוגית לחיכוך במצוקים, על מנת לתאר נוזלים צמיגים יש לתקן את משוואת התנועה. משוואת התנועה הכוללת את התיקון לנוזלים צמיגים נקראת משוואת נאוויה-סטוקס, והתיקון הוא הוספת האיבר לאגף ימין של משוואת התנועה של אוילר, כאשר הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \eta} הוא מקדם הצמיגות הדינמי. הסימון הוא הפעלת אופרטור הלפלסיאן על כל איבר של הווקטור .

המשמעות של האיבר הנוסף היא מעין דיפוזיה של תנע. משוואת הדיפוזיה היא מהצורה הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle\alpha\nabla^{2}n=\frac{\partial n}{\partial t}} , ותפקיד תהליך הדיפוזיה הוא למצע את פני גרדיאנטים גדולים בצפיפות הנוזל הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle n} . באנלוגיה למקרה הצפיפות בנוזל ניתן להסיק שהאיבר החדש, הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \textstyle -\frac{\eta}{\rho}\nabla^{2}\overline{v}} , ממצע על פני גרדיאנטים גדולים במהירות בנוזל, וזהו עיקר תהליך הצמיגות.

ראו גם

רקע למשוואות אוילר

מעבר למשוואות אוילר

לקריאה נוספת

  • Landau, L. D. and Lifschitz, E. M., Fluid Mechanics, Second Edition, Pergamon Press, ISBN 978-0750627672
  • Julián Simón Calero, The Genesis of Fluid Mechanics, 1640–1780, Studies in History and Philosophy of Science 22, 2008, ISBN 978-1-4020-6413-5
  • Darrigol, Olivier, Worlds of Flow. A History of Hydrodynamics from the Bernoullis to Prandtl, Oxford University Press. chapter 1, ISBN 978-0199559114

הערות שוליים

  1. ^ Julián Simón Calero, "Theoretical Constructions (II): Euler" in "The Genesis of Fluid Mechanics, 1640–1780", Studies in History and Philosophy of Science 22, 2008, ISBN 978-1-4020-6413-5, page 401
  2. ^ * Darrigol, Olivier, Worlds of Flow. A History of Hydrodynamics from the Bernoullis to Prandtl, Oxford University Press. pp. 2
  3. ^ Olivier Darrigol and Uriel Frisch, From Newton's mechanics to Euler's equations, Euler Equations: 250 years on
  4. ^ The Euler Archive, E258 -- Principia motus fluidorum
  5. ^ The Euler Archive E225 -- Principes generaux de l'etat d'equilibre des fluides
  6. ^ Leonhard Euler, Principes generaux de l'etat d'equilibre des fluides, pp. 1
  7. ^ The Euler Archive, E226 -- Principes generaux du mouvement des fluides
  8. ^ The Euler Archive, E227 -- Continuation des recherches sur la theorie du mouvement des fluides
  9. ^ Leonhard Euler, Continuation des recherches sur la theorie du mouvement des fluides, pp. 1
  10. ^ האתר של פרופסור Thayer Watkins, אוניברסיטת סן־חוזה
  11. ^ Tilley, David R. and Tilley, John., Superfluidity and Superconductivity, Third Edition. IOP Publishing. chapters 3.2-3.3
  12. ^ D. McMahon, Quantum Field Theory, Mc Graw Hill, ISBN 978-0-07-154382-8