פרדוקס ד'אלמבר

ז'אן לה רון ד'אלמבר, הוגה הפרדוקס
ערכים אמפיריים עבור מקדם גרר כפונקציה של מספר ריינולדס. הקו המלא מייצג ספירה חלקה והקו השבור מייצג ספירה מחוספסת. המספרים לאורך הקווים מייצגים שינוי במשטר הזרימה והשינויים במקדם הגרר בהתאם: 2: זרימה צמודה (זרימת סטוקס) וזרימה מנותקת יציבה 3: זרימה מנותקת עם שכבת גבול למינרית במעלה ההיתנקות היוצרת שובל ערבולים סדורים 4: זרימה מנותקת עם שיכבת גבול במעלה ההיתנתקות והזרימה בשובל היא [[זרימה טורבולנטית 5: זרימה מנותקת על קריטית עם שכבת גבול טורבולנטית]]

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

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

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

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

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

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

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

== היפרדות של זרימה לא צמיגה: קירכהוף וריילי ==
Batchelor (2000), p. 499, eq. (6.13.12).}} עם לחץ קבוע בשני קווי זרם המופרדים על גבי הגוף

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

עם זאת, התנגדויות בסיסיות התעוררו כנגד גישה זו: קלווין ציין שאם משטח נע במהירות קבועה דרך הנוזל, והמהירות בשובל שווה לזו של המשטח. ההיקף האינסופי של השובל - המתרחב ככל שהמרחק מהמשטח מתרחק כפי שמתקבל מהתאוריה - משמעותו שיש אנרגיה קינטית אינסופית בשובל. משמעות שאנחנו חייבים לדחות על בסיס פיזיקלי. יתר על כן, הפרש הלחצים הנצפה בין החלק הקדמי והחלק האחורי של המשטח וכוחות הגרר כתוצאה מכך, הם הרבה יותר גדולים מהצפוי: למשטח שטוח בניצב לזרימה מקדם הגרר המחושב הוא <math>C_D=0.88</math>, ואילו בניסויים <math>C_D=2.0 </math>. זאת, בעיקר בשל שאיבה בצד השובל של המשטח, הנגרם על ידי הזרימה הבלתי יציבה בשובל האמיתי (בניגוד לתאוריה שמניחה מהירות זרימה קבועה ששווה למהירות של המשטח). לכן, תאוריה זו אינה מספקת כדי להסביר את הגרר שפועל על גוף בזרימה.
== שכבות גבול דקות: פרנטל ==
חלוקת הלחץ בזרימה סביב צילינדר מעגלי. הקווים הכחולים המקווקווים מראים את חלוקת הלחץ לפי האוריית הזרימה הפוטנציאלית, הגורמת לפרדוקס ד'אלמבר. הקווים הכחולים הרציפים מגיעים בניסויים שנעשו במספרי ריינולדס גבוהים.
הפיזיקאי הגרמני לודוויג פרנטל הציע בשנת 1904 כי ההשפעות של שכבת גבול צמיגה דקה אולי יכולה להיות המקור של גרר משמעותי. פרנטל שיער, שבמהירויות גבוהות ומספרי ריינולדס גבוהים, תנאי שפה של אי-החלקה גורם לווריאציה משמעותית של מהירויות הזרימה בשכבה דקה בקרבת הגוף. זה מוביל ליצירה של ערבוליות ופיזור צמיג של אנרגיה קינטית בשכבת הגבול. פיזור האנרגיה, אשר חסרה בתאוריות של זרימה לא צמיגה, מסבירה את ההפרדה של הזרימה על גבי גופים שאינם פשוטים. הלחץ הנמוך באזור השובל יוצרות גרר צורה, אשר יכול להיות גדול יותר מאשר גרר החיכוך בשל לחץ הגזירה הצמיגי על הקיר.

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

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

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

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

ב-1951 המתמטיקאי ג'יימס ג'יי סטוקר מבקר בחריפות את הפרק הראשון בספרו של בירקהוף:

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

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

== הוכחה של אפס גרר בזרימה פוטנציאלית תמידית ==
=== זרימה פוטנציאלית ===
קווי זרם של זרימה פוטנציאלית מסביב לצילינדר עגול בזרימה אחידה.
שלוש ההנחות העיקריות בגזירה של פרדוקס ד'אלמבר הן שהזרימה הקבועה היא בלתי דחיסה, לא צמיגה, ואי רוטציונית. נוזל לא צמיג מתואר על ידי משוואות אוילר, שבזרימה בלתי דחיסה נראות כדלהלן:
:<math>\begin{align} & \boldsymbol{\nabla} \cdot \boldsymbol{u} = 0 && \text{(Conservation of mass)} \\ & \frac{\partial}{\partial t} \boldsymbol{u} + \left(\boldsymbol{u} \cdot \boldsymbol{\nabla}\right) \boldsymbol{u} = - \frac{1}{\rho} \boldsymbol{\nabla} p && \text{(Conservation of momentum)} \end{align}</math>

כאשר <math>u</math> הוא מהירות הזורם, <math>p</math> הלחץ, <math>\rho</math> הצפיפות, ו<math>\nabla</math> הוא אופרטור הגרדיאנט. ההנחה שהזרימה היא רוטציונית גוררת שהמהירות מקיימת <math>\nabla\times u=0</math> ולכן אנחנו מקבלים:
:<math>\left(\boldsymbol{u} \cdot \boldsymbol{\nabla}\right) \boldsymbol{u} = \tfrac12 \boldsymbol{\nabla} \left(\boldsymbol{u} \cdot \boldsymbol{u}\right) - \boldsymbol{u} \times \boldsymbol{\nabla} \times \boldsymbol{u} = \tfrac12 \boldsymbol{\nabla} \left(\boldsymbol{u} \cdot \boldsymbol{u}\right) \qquad (1)</math>
כאשר השוויון הראשון הוא זהות והשוויון השני נובע מהעובדה שהזרימה אי-רוטציונית. בנוסף, לכל זרימה אי-רוטציונית, קיימת פונקציית מהירות <math>\varphi</math> כך ש<math>u=\nabla \varphi</math>. הכנסת כל זה למשוואת שימור התנע נותן:
:<math>\boldsymbol{\nabla} \left( \frac{\partial\varphi}{\partial t} + \tfrac12 \boldsymbol{u} \cdot \boldsymbol{u} + \frac p\rho \right) = \boldsymbol{0}</math>
לכן, החלק שבסוגריים חייב להיות קבוע (כל תלות שיש בזמן יכולה להתבטל על ידי הגדרה מגדש של <math>\varphi</math>). בהנחה שהזורם נמצא במנוחה באינסוף ושהלחץ מוגדר כאפס שם, הקבוע שווה לאפס, ולכן:
<math> \frac{\partial\varphi}{\partial t} + \tfrac12 \boldsymbol{u} \cdot \boldsymbol{u} + \frac p\rho = 0, \qquad (2) </math>:

שזה משוואת ברנולי לזרימה פוטנציאלית לא יציבה.

=== אפס גרר ===
אם גוף נע במהירות קבועה <math>v</math> דרך הנוזל שנמצא במנוחה רחוק מאוד - באינסוף. לכן שדה המהירות של הנוזל חייב להקבע על ידי הגוף, ולכן הוא מהצורה <math>u(x,t)=u(x-v\cdot t,0)</math> כאשר <math>x</math> הוא הקוארדינטה של הווקטור המרחבי, ולכן:
:<math>\frac{\partial \boldsymbol{u}}{\partial t} + \left( \boldsymbol{v} \cdot \boldsymbol{\nabla} \right) \boldsymbol{u} = \boldsymbol{0}. </math>

מכיוון ש<math>u=\nabla \varphi </math>, אפשר לעשות אינטגרציה לפי <math>x</math> ולקבל:
:<math>\frac{\partial\varphi}{\partial t} = -\boldsymbol{v} \cdot \boldsymbol{\nabla} \varphi + R(t) = -\boldsymbol{v} \cdot \boldsymbol{u} + R(t). </math>

הכוח <math>F</math> שהנוזל מפעיל על הגוף מחושב באמצעות אינטגרל משטחי:
:<math>\boldsymbol{F} = - \int_A p\, \boldsymbol{n}\; \mathrm{d} S </math>

כאשר <math>A</math> הוא שטח הגוף, ו<math>n</math> הוא וקטור הנורמל למישור הגוף. אבל ממשואה (2) אנחנו מקבלים:
:<math> p = - \rho \Bigl( \frac{\partial\varphi}{\partial t} + \tfrac12 \boldsymbol{u} \cdot \boldsymbol{u} \Bigr) = \rho \Bigl( \boldsymbol{v} \cdot \boldsymbol{u} - \tfrac12 \boldsymbol{u} \cdot \boldsymbol{u} - R(t) \Bigr), </math>

לכן:
:<math> \boldsymbol{F} = - \int_A p\, \boldsymbol{n}\; \mathrm{d} S = \rho \int_A \left(\tfrac12 \boldsymbol{u} \cdot \boldsymbol{u} - \boldsymbol{v} \cdot \boldsymbol{u}\right) \boldsymbol{n}\; \mathrm{d} S, </math>

כאשר התרומה של <math>R(t)</math> לאינטגרל שווה לאפס.
בשלב זה נהיה הרבה יותר פשוט לעבוד עם רכיבי הווקטור. הרכיב ה<math>k</math> של המשוואה הזאת הוא:
:<math> F_k = \rho \int_A \sum_i (\tfrac12 u_i^2 - u_i v_i) n_k \, \mathrm{d} S. \qquad (3) </math>

אם <math>V</math> הוא הנפח של הנוזל, אזי על פי משפט גאוס:
:<math>\frac12 \int_A \sum_i u_i^2 n_k \, \mathrm{d} S = - \frac12 \int_V \frac{\partial}{\partial x_k} \left( \sum_i u_i^2 \right) \,\mathrm{d} V. </math>

צד ימין הוא אינטגרל על נפל אינסופי ולכן צריך הצדקה, אשר אפשר להשיגה מתאוריית הזרימה הפוטנציאלית אשר אומרת כי המהירות דועכת כמו <math>r^{-3}</math> - כמו פוטנציאל חשמלי של דיפול בגוף סופי תלת־ממדי. האינטגרנד בתוך האינטגרל הנפחי יכול להכתב כדלהלן:
:<math>\frac12 \frac{\partial}{\partial x_k} \left( \sum_i u_i^2 \right) = \sum_i u_i \frac{\partial u_k}{\partial x_i} = \sum_i \frac{\partial(u_iu_k)}{\partial x_i}</math>

כאשר שוויון (1) והאי-דחיסות של הנוזל ממומשים. כאשר מחזירים את הנוסחא הזו בחזרה לתוך האינטגר הנפחי ומשתמשים בחוק גאוס שוס, מקבלים:
:<math> - \frac12 \int_V \frac{\partial}{\partial x_k} \left( \sum_i u_i^2 \right) \,\mathrm{d} V = -\int_V \sum_i \frac{\partial(u_iu_k)}{\partial x_i} \,\mathrm{d} V = \int_A u_k \sum_i u_i n_i \,\mathrm{d} S.</math>

כאשר מכניסים את זה למשוואה (3) אנו מוצאים ש:
:<math> F_k = \rho \int_A \sum_i (u_k u_i n_i - v_i u_i n_k) \, \mathrm{d} S. </math>

הנוזל אינו יכול לחדו אל תוך הגוף, ולכן <math>n\cdot u=n\cdot v</math> על משטח הגוף, ולכן:
:<math> F_k = \rho \int_A \sum_i (u_k v_i n_i - v_i u_i n_k) \, \mathrm{d} S. </math>
ולבסוף, כוח הגרר הפועל על הגוף חייב להיות בכיוון התנועה, ולכן:
:<math> F_{drag}= \boldsymbol{v} \cdot \boldsymbol{F} = \sum_i v_i F_i = 0.</math>

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

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

קטגוריה:ערכים שבהם תבנית בריטניקה אינה מתאימה
ד'אלמבר
קטגוריה:מכניקת הזורמים