מבוא
עד כה עסקנו בבעיות הולכה מהסוג:
ופתרנו את משוואות אלו כדי למצוא את פילוג הטמפרטורה, ומשם את קצב מעבר החום. הייחוס היחידי שלנו להסעה היה כאשר התחשבנו בו כתנאי שפה למוצק עליו אנו כותבים את המשוואה. נזכיר כי הסעה כוללת את המעבר אנרגיה גם ע”י תנועה מקרוסקופית של הזורם, וגם ע”י תנועה רנדומלית של החלקיקים.
נקודה מעניינת בהסעה היא שהאופן הסעה של מסה אנלוגית לאופן הסעה של מעבר חום. כלומר, כמו בזרימה, נעסוק כאן בשכבות גבול, זרימה טורבולנטית שכבתית (למינרית), זרימה מול פלטה שטוחה וכו’.
שכבות גבול
עיקרון השכבות גבול הוא חיוני להבנת ההסעה של מסה וחום בין משטח וזורם. לכן, נעבור בקצרה על שכבת גבול מהירות, לפני שנעבור לשכבת גבול תרמית.
שכבת גבול מהירות
כאשר זרימה מציפה פוגשת בפלטה שטוחה, על הפלטה עצמה, הזרם נעצר - לפי תנאי האי-החלקה. חלקיקים אלו לאחר מכן פועלים כנגד המהירות של החלקיקים בשכבה מעליהם, וחלקיקים אלו מאטים גם כן את החלקיקים בשכבה מעליהם - מה שיוצר בסוף שכבת גבול.
התפתחות שכבת גבול מהירות על פלטה שטוחה. (Bergman & Lavine, 2017).
מקדם החיכוך המקומי על הפלטה מוגדר בעזרת המאמץ גזירה
כאשר מאמץ הגזירה מוגדר בעזרת הצמיגות הדינמית של הזורם:
כדאי לעבור על שאר הנושאים בשכבת גבול בזרימה (הנחות לשכבת גבול, המשוואות, הפתרונות למשוואות וכו’) לפני שעוברים לשכבת גבול תרמית.
שכבת גבול תרמית
כמו ששכבת גבול מהירות מתפתחת כאשר יש זרימת זורם על משטח, שכבת גבול תרמית מתפתחת אם יש הפרש טמפרטורות בין הזרימה המציפה של הזורם והמשטח. נביט באיור הבא:
שכבת גבול תרמית המתפתחת על משטח איזותרמי. (Bergman & Lavine, 2017).
בדופן הלוח, פרופיל הטמפרטורה הוא אחיד, עם
נפתח כעת את הקשר בין התנאים בשכבת גבול זו ומקדם ההסעה
כאשר
בשפה, אין תנועת זורם, והמעבר חום מתרחש אך ורק ע”י הולכה. מחוק הקירור של ניוטון:
ולכן:
לפיכך, התנאים בשכבת גבול, המשפיעים על גרדיאנט הטמפרטורה
זרימה טורבולנטית ולמינרית
בזרימה התמקדנו באזור הלמינרי של שכבות הגבול, אבל חשוב לנו כאן להתייחס גם לשכבת הגבול הטורבולנטית.
התפתחות שכבת גבול מהירות על פלטה שטוחה. (Bergman & Lavine, 2017).
בשכבת גבול הלמינרית, הזרם מסודר באופן שכבתי וניתן לזהות קווי זרם לאורכם החלקיקי זורם נעים. הסדר הזה ממשיך עד שמגיעים ל-אזור מעבר, שלאורכו מתרחש מעבר מזרימה למינרית לטורבולנטית.
הזרם בשכבת גבול הטורבולנטית, בכללי, הוא מאוד מפוזר, והוא מאופיין בתנועה תלת ממדית רנדומלית. הערבוב בתוך השכבת גבול לוקח זורם במהירות גבוהה לשפת המוצק ומעביר זורם במהירות נמוכה הרחק יותר אל הזרימה המציפה.
השוואה בין פרופיל זרימה של שכבת גבול למינרית וטורבולנטית עבור אותה הזרימה המציפה. מאחר והמהירות משתנה בזמן בזרימה טורבולנטית, המהירות הממוצעת,
, משורטטת בגרף. (Bergman & Lavine, 2017).
כדי לחשב האם השכבת גבול טורבולנטית או למינרית, נהוג להניח שהמעבר מתרחש במיקום קריטי
אלא אם כן מצוין אחרת.
משוואות שכבת הגבול
בזרימה חישבנו את עובי שכבת הגבול, המקדם חיכוך ועוד פרמטרים נוספים המאפיינים את שכבת הגבול בעזרת משוואות שכבת הגבול. הגענו במהלך הדרך למאזן תנע קווי בכיוון
באותו אופן, אנו יכולים גם לפתח את המשוואה הבאה משימור אנרגיה:
הביטויים בצד שמאל מתארים את סך הקצב בו אנרגיה תרמית עוזבת את נפח הבקרה כתוצאה מתנועה של זורם. הביטוי הראשון בצד ימין מתאר את סך האנרגיה התרמית הנכנסת כתוצאה מהולכה בכיוון
בדרך כלל מזניחים איבוד לחיכוך בתוך הזרם, כי
^eq-*6-29
משוואות שכבת הגבול המנורמלות
כמו בזרימה, נוח לנו לעבוד עם גדלים חסרי ממד, ולכן נסמן:
תנאי שפה במונחים של
נציב לתוך משוואת שימור האנרגיה (*6.29):
כאשר:
הגדרה: מספר פרנדטל
הביטוי
נקרא מספר פרנדטל (Prandtl), והוא מוגדר כ: כאשר
הוא ההולכה התרמית של ה-זורם ( עבור fluid). מספר זה בעצם מתאר את היחס בין כמה מהר התנע מתפזר, וכמה מהר הטמפרטורה מתפזרת.
משימור האנרגיה נסיק שהפתרון למשוואה (IH6.36) הוא מהצורה:
מהגדרת מקדם ההסעה, משוואה (IH6.5):
במונחים חסרי ממד:
ביטוי זה מרמז על תלות במספר חסר ממד הנקרא:
הגדרה: מספר נוסלט
מספר נוסלט מוגדר באופן הבא:
מספר זה שווה לגרדיאנט הטמפרטורה החסר ממד בשפה, והוא מספק מידה לכמות המעבר חום בהסעה המתרחש בשפה.
ממשוואה (IH6.47), נסיק ש-עבור גאומטרייה נתונה:
כמו ש-
אנלוגיות שכבות גבול
כמהנדסים, עיקר העניין שלנו בשכבות גבול הוא במציאת המספרים החסרי ממד המאפיינים אותם כמו
אנלוגיית ריינולדס
מטבלה 6.1 ניתן לראות שעבור
נגדיר עוד מספר חסר ממד:
הגדרה:
מספר סטנטון (Stanton) מוגדר כ:
וכעת נוכל לרשום את (IH6.66) בצורה:
משוואה זו נקראת אנלוגיית ריינולדס, והיא מקשרת בין השכבת גבול מהירות לשכבת הגבול התרמית. משוואה זו נכונה עבור המקרה של
קיים קשר אמפירי (מבוסס על ניסויים) שאומר שעבור
פלטה שטוחה תחת זרימה מציפה
למרות פשטותה, זרימה מציפה על פלטה שטוחה מתרחשת במספר רב של בעיות הנדסיות. כפי שכבר הזכרנו, זרימת למינרית מתחרשת בדופן הפלטה, ומתפתחת לזרימה טורבולנטית במיקום
זרימה מציפה על פלטה שטוחה. (Bergman & Lavine, 2017).
ננתח את הזרימות השונות מבחינת מהירות וטמפרטורה בעזרת פתרון בלסיוס:
זרימה למינרית בפלטה שטוחה
מפתרון בלסיוס לזרימה, עובי שכבת הגבול הוא
ומקדם החיכוך:
הממוצע:
באותו אופן, ניתן למצוא בעזרת פתרון בלסיוס את
והממוצע:
נזכור כי מספר פרנדטל מתאר את היחס בין יעילות המעבר אנרגיה ומומנטום ע”י דיפוזיה בשכבת המהירות לבין יעילות המעבר אנרגיה ומומנטום ע”י דיפוזיה בשכבה התרמית. לכן, נצפה שתהיה לו השפעה לקשר בין עובי שכבת גבול המהירות לעובי שכבת הגבול התרמית.
ואכן, מהשוואה בין הפתרונות, עבור פלטה שטוחה:
זרימה טורבולנטית בפלטה שטוחה
עבור זרימה טורבולנטית אין תקווה למצוא פתרון באופן אנליטי, ולכן פונים לניסויים.
עבור זרימה טורבולנטית ישרה, כאשר
עובי שכבת הגבול ניתן לקירוב כ:
מהשוואת ביטויים אלו למשוואות המתאימות בזרימה למינרית, (IH7.19) ו-(IH7.20), נסיק כי גדילת השכבת גבול באזור הטורבולנטי הרבה יותר מהירה, והירידה במקדם החיכוך יותר מתונה. עבור זרימה טורבולנטית, התפתחות שכבת הגבול מושפעת משינויים רנדומליים בזורם, ולא מדיפוזיה מולקולרית. לכן, קצב גדילת השכבת גבול התרמית לא תלוי ב-
מאנלוגיית ריינולדס ומשוואה (IH7.34), מספר נוסלט המקומי הופך להיות:
והממוצע:
זרימה מעורבת
עבור זרימה למינרית המתרחשת לכל אורך הפלטה, ניתן להשתמש במשוואות (IH7.29) ו-(IH7.30) לחישוב המקדמים הממוצעים. למעשה, גם אם מתרחש מעבר לזרימה טורבולנטית בסוף הפלטה, למשל בטווח
במקרה המעורב, נחשב את המקדם מעבר חום בהסעה הממוצע באופן הבא:
כאשר, לפי משוואות (IH7.23) ו-(IH7.36):
מאינטגרציה נקבל:
כאשר משוואה זו תקפה כאשר
באותו אופן ניתן לחשב כוח חיכוך ממוצע. אם
לשכבת גבול שהיא טורבולנטית לחלוטין, מתקיים
דוגמה: 7.2 בספר
נתון לוח שטוח ברוחב מטר המוחזק בטמפרטורה אחידה וקבועה
ע”י סדרת מחממים באורך . הלוח חשוף לאוויר במהירות ובטמפרטורה .
איזה מהמחממים דורש הספק מקסימאלי וכמה?פתרון:
עבור אוויר, מטבלה A.4 בספר:מספר ריינולדס לחלק המחמם הראשון, באורך
: במידה וריינולדס הקריטי הוא
, נסיק שהמעבר לטורבולנציה יתרחש ב: נסיק שהמעבר מתרחש במחמם החמישי. נסיק שישנם שלוש אפשרויות למחמם קריטי הדורש את ההספק המקסימלי:
- המחמם הראשון
- המחמם החמישי
- המחמם השישי
בכללי, עבור המחמם ה-
-י, ההספק חום צריך להיות: לפי שימור אנרגיה, ההספק הדרוש למחמם
ניתן לחישוב ע”י חיסור קצב איבוד החום של ה- מחממים הראשונים מקצב איבוד החום של כלל ה- מחממים. אם מגדירים את כמקדם החום להסעה הממוצע למחממים עד , ההספק חום הדרוש למחמם שווה לקצב מעבר החום בהסעה מהמחמם, שניתן לביטוי כ: עבור המחמם הראשון, מספר נוסלט:
ולכן מקדם מעבר החום:
עבור המחמם החמישי, לפי משוואה (IH7.38):
לכן המקדם מעבר חום במחממים
עד : את הערך של
נוכל למצוא באותו אופן שמצאנו את , ואז נמצא כי: ממשוואה (3):
באותו אופן מחשבים עבור
ונקבל בסוף: מאחר ו-
, נסיק שהמחמם השישי ידרוש הכי הרבה הספק חום, שהוא:
אורך התחלה ללא חימום
מספרי נוסלט לעיל תקפים למקרים בהם טמפרטורת השפה
יש מספר מקרים יוצאי דופן שנפוצים גם הם, כמו למשל הקיום של אורך התחלה ללא חימום (
פלטה שטוחה במקביל לזרם עם אורך התחלה ללא חימום. (Bergman & Lavine, 2017).
לכן, אין מעבר חום ב-
כאשר
מבחינת מספר נוסלט ממוצע, עבור זרימה למינרית או טורבולנטית:
כאשר
פלטה שטוחה בשטף חום קבוע
קיימת גם האפשרות של שטף חום קבוע במקום טמפרטורה קבועה. במקרים אלו, עבור זרימה למינרית:
ולזרימה למינרית:
מבחינת מספר נוסלט ממוצע, ניתן להראות עבור זרימה למינרית ש:
שזה גדול רק ב-
הסעה חיצונית על משטחים עקומים
בזרימה דיברנו על ניתוק זרימה - תופעה המתרחשת בגאומטריות לא ישרות, שזה תאכלס כמעט כל הגאומטריות.
התפתחות שכבות גבול והתנתקות על גליל עגול. (Bergman & Lavine, 2017).
לאזור מאחורי הגליל באיור לעיל אנו קוראים שובל (wake).
נקודת ההתנתקות תלויה ב-
ואנו נניח שעבור גליל, כאשר
הסעה על גליל וכדור
בגרף הבא מוצגים תוצאות ניסויים למספר נוסלט לפי
מספר נוסלט מקומי לזרימה הנורמלית לגליל עגול. (Bergman & Lavine, 2017).
לא במפתיע, תוצאות אלו משופעות מטבע התפתחות השכבות גבול על שטח הפנים של הגליל. נביט למשל בגרף של של
לא ניתן לקבל פתרון אנליטי עבור גלילים במספר
כאשר
ישנה עוד קורלציה שמכילה את כלל טווחי ה-
כאשר כל הערכים מחושבים בטמפרטורה הממוצעת
עבור כדור:
כאשר
זרימה לאורך מערך גלילים
מעבר חום ממערך גלילים בזרימה מציפה רלוונטי למספר יישומים תעשייתיים, כמו יצירת קיטור בדוד או קירור בצינורות של מזגן. גאומטריית מערך גלילים טיפוסית היא מהצורה:
סכמה של צינורות בזרימה מציפה. (Bergman & Lavine, 2017).
בדרך כלל, זורם אחד עובר על הצינורות, בעוד זורם שוני בטמפרטורה שונה זורם דרך הצינורות. בפרק זה אנו דנים ספציפית בהסעה של הזרימה המציפה על הצינורות.
המערך גלילים יכול להיות מסודר או בצורה ישרה (aligned) או מדורגת (staggered) לכיוון הזרימה שבמהירות
אופן סידור מערך הגלילים. (a) מיושר. (b) מדורג. (Bergman & Lavine, 2017).
התצורה מאופיינת ע”י קוטר גליל
זרימה סביב הגלילים בשורה הראשונה של המערך דומה לזו סביב גליל אחד. לכן, מקדם מעבר החום עבור המערך בשורה הראשונה הוא בערך שווה למקדם עבור גליל אחד. עבור השורות הבאות, זה כבר תלוי באופן הסידור של מערך הגלילים.
תנאי זרימה עבור (a) סידור ישר ו-(b) סידור מדורג. (Bergman & Lavine, 2017).
גלילים מיושרים מעבר לשורה הראשונה נמצאים בשובלים של הגלילים במעלה הזרימה, ועבור ערכי
עבור גלילים מדורגים, מסלול הזרימה יותר כאוטי, וערבוב הזרימה גדל ביחס למערך הגלילים המיושר. בכללי, מעבר החום יותר טוב עבור הזרימה היותר כאוטית של מערך מדורג, בייחוד עבור ערכי ריינולדס קטנים (
לרוב, נרצה לדעת את מקדם מעבר החום ה-ממוצע של כלל מערך הגלילים. מניסויים, נמצא הקשר הבא:
כאשר
אם יש
כאשר
מספר ריינולדס המקסימלי
למערך ישר,
למערך מדורג, המהירות המקסימלית מתרחשת או ב-
לכן,
ואז:
אם
מאחר והזורם יכול לחוות שינוי גדול בטמפרטורה כאשר הוא נע דרך מערך הצינורות, קצב מעבר החום כבר לא פרופורציונלי ישירות ל-
כאשר
כאשר
תרגילים
תרגיל 1
נתונים ביטויים למקדמי מעבר חום בהסעה מקומיים עבור זרימה מקבילה ללוח ישר, גם למינרית וגם טורבולנטית:
סעיף א’
מהו מקדם מעבר החום הממוצע בהסעה על גבי פלטה באורך
פתרון:
לפי זרימה מעורבת:
נחשב כל אינטגרל בנפרד:
נקבל:
סעיף ב’
האם קיים אורך
פתרון:
נמצא מתי
נציב נתונים ונקבל:
תרגיל 2
כדי למנוע היווצרות קרח על כנפי מטוס, מוצע להתקין גופי חימום על הכנפיים. האורך האופייני של חתך הכנף הוא
תכונות האוויר:
פתרון:
שטף החום הממוצע הדרוש הוא:
בהזנחת גרדיאנט לחץ בסביבת המטוס, ניתן להעריך את מקדם מעבר החום בהסעה באמצעות אנלוגיית ריינולדס:
נחלץ את מקדם מעבר החום הממוצע בהסעה:
נשים לב כי:
ולכן:
מה שאומר שמקדם מעבר החום הממוצע:
ומכאן:
שאלה 3
אוויר בלחץ אטמוספרי וטמפרטורה
סעיף א’
מצאו את עובי שכבת הגבול התרמית בקצה הלוח.
פתרון:
ראשית, נמצא את תכונות האוויר:. עבור הטמפרטורה הממוצעת:
התכונות הן:
נמצא אם אנחנו בתחום טורבולנטי או למינרי:
לפי פתרון בלסיוס, עבור האזור הלמינרי, משוואה (IH7.19):
ולכן בקצה הלוח:
את עובי שכבת הגבול התרמית נמצא בעזרת משוואה (IH7.24):
נציב ערכים ונקבל:
סעיף ב’
מהו שטף החום הממוצע על פני הלוח?
פתרון:
לפי הגדרת שטף החום:
את מספר נוסלט נמצא ממשוואה (IH7.30):
ולכן שטף החום:
סעיף ג’
בניסוי נמדד על הלוח כוח הגרר
פתרון:
מתוך הגדרת כוח הגרר ומקדם החיכוך:
מנגד, מאנלוגיית ריינולדס, משוואה (IH6.70):
נשווה בין שני הביטויים ל-
לכן שטף החום:
נציב ערכים ונקבל:
שאלה 4
משטח של גוף חימום מוחזק בטמפרטורה של
סכמת המשטחים
סעיף א’
חשבו את קצב מעבר החום הכולל ליחידת רוחב מן המשטח כאשר הפרוסה המחוספסת ממוקמת ב-מעלה הזרם.
פתרון:
תכונות האוויר עבור
לפי זרימה מעורבת, משוואה (IH7.38):
מאחר והזרימה תהיה טורבולנטית מההתחלה, מהגדרת
ולכן קצב המעבר חום ליחידת רוחב:
נקבל:
סעיף ב’
חשבו את קצב מעבר החום הכולל ליחידת רוחב מן המשטח כאשר הפרוסה המחוספסת ממוקמת ב-מורד הזרם.
פתרון:
באותו אופן כמו סעיף קודם, רק הפעם נשים לב שאנחנו מתחילים מזרימה למינרית, ועלינו לבדוק אם אנו בכלל מגיעים לזרימה טורבולנטית:
קיבלנו ש-
כלומר, קיבלנו שהזרימה משתנה לטורבולנטית איפשהו באמצע, לפני המעבר למשטח המחוספס. משוואה (IH7.38) עדיין תקפה, אבל הפעם:
נציב ונקבל:
ולכן קצב המעבר חום:
סעיף ג’
חזרו על הסעיף הקודם אבל עבור זרימת אוויר במהירות
פתרון:
נבדוק שוב אם עד חצי מהאורך הכולל
נסיק שעד מחצית מהאורך הכולל אין מעבר לזרימה טורבולנטית. אבל, מאחר ואנו עוברים לפלטה המחוספסת, נעבור בכל זאת לאזור טורבולנטי ב-
לפי משוואה (IH7.39):
שוב, לפי משוואה (IH7.38):
לכן מקדם מעבר החום הוא: