זרימה בצנרת

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

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

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

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

## משוואת ברנולי – חוק שימור האנרגיה

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

בהזנחת החיכוך ניתן להניח שהאנרגיה של הזורם נשמרת בזמן תנועתו בצינור. על פי משוואת ברנולי: 
P+\rho gh+\frac{1}{2}\rho =const


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

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

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

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

בזרימה בצנרת:
* \operatorname{Re}\le 2300 - זרימה למינארית
* 2300 - אזור מעבר בין זרימה למינארית לטורבולנטית.
* \operatorname{Re}\ge 4000 - זרימה טורבולנטית
כאשר: 
\operatorname{Re}=\frac{\rho UD}{\mu }=\frac{UD}{\nu }

## זרימה למינרית מפותחת בצינור עגול
זרימה למינרית בצינור תתקיים עבור: \operatorname{Re}\le 2300 והיא זרימה המונעת על ידי מפל לחצים מסוג זרימת פואסיל.

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

447x447px

לאחר חלוקה ב 2\pi r מתקבל: \frac{\partial p}{\partial x}=\frac}{r}+\frac{d}{dr}=\frac{1}{r}\cdot \frac{d\left( r \right)}{dr}

מאמץ הגזירה מוגדר להיות:  
=\mu \frac{\partial u}{\partial r}

ופתרון המשוואה הדפרנציאלית הואו: 
\mu \frac{du}{dr}==\frac{r}{2}\left( \frac{\partial p}{\partial x} \right)+\frac}{r}

מכיוון שמאמץ גזירה איסופי לא ייתכן, אזי C_1=0 ומכאן נקבל ביטוי למאמץ הגזירה: 
=\frac{r}{2}\left( \frac{\partial p}{\partial x} \right)

ולכן, על פי החוק השלישי של ניוטון, **הגזירה שמפעילה הדופן על הזורם היא**:  
=-\frac{R}{2}\left[ \frac{\partial p}{\partial x} \right]

## זרימה טורבולנטית מפותחת בצינור עגול

זרימה טורבולנטית תתקיים עבור: \operatorname{Re}\ge 4000

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

ניתן לחלק את המערבולות הללו למהירויות **ממוצעות **(לכן בנוסחה מופיע קו מעליהן) בזמן בכיוון x ובכיוון r .הפחתה זו של מאמץ הגזירה קיבלה את השם "מאמצי ריינולדס".

לכן, מאמצי הגזירה בזרימה טורבולנטית מפותחת הם: 
=\mu \frac{\partial u}{\partial x}-\rho \overline}

מתוך ניסויים שנעשו במספרי ריינולדס שונים, התגלה כי תופעה זו מתאפסת בסמוך לדפנות הצינור. לכן על הדופן, על אף שהזרימה טורבולנטית: 
=\mu \frac{\partial u}{\partial r}.

התופעה אינה זניחה בשכבות פנימיות בצינור. 

## קוטר אפקטיבי (הקוטר ההידראולי)
מכיוון שלרוב היישומים ההנדסיים משתמשים בצינורות עגולים, הניתוח יעשה בקואורדינטות גליליות ובחתך עגול. הפתרונות המתקבלים נכונים גם לצינורות מרובעים ומשולשים תחת הדרישה הבאה:  \frac{height}{width}\le 4. כלומר, יחס גובה לרוחב לא גדול מ-4.

עבור כל צינור ניתן לחשב קוטר הידראולי: =\frac{4\cdot Area}{Perimeter}. קוטר זה ישמש בתפקיד הקוטר בכל החישובים בהמשך.

## החיכוך בצנרת
 

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

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

בעזרת משוואת ברנולי ניתן לתאר את הפסדי החיכוך כך:         
\left( \frac}{\rho }+\frac^{2}}{2}+g \right)-\left( \frac}{\rho }+\frac^{2}}{2}+g \right)=

כאשר המהירויות שוות ותכונות הזורם נשארות קבועות ניתן לקבל את התוצאה המופשטת: 
\frac-}{\rho }=g(-)+

בנוסף, אם הצינור אופקי, אזי הגובה איננו משתנה ונקבל: 
\frac-}{\rho }=\frac{\Delta p}{\rho }=

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

בזרימה למינרית, מכיוון ש \frac{\partial p}{\partial x}=const=\frac{\Delta p}{\Delta x}.

ניתן לתאר את הפסדי הלחץ בצורה פשוטה יחסית: 

\Delta p=32\frac{L}{D}\cdot \frac{\mu \bar{V}}{D}


ולכן: 

=\frac{\Delta p}{\rho }=\left( \frac{64}{\operatorname{Re}} \right)\frac{L}{D}\frac}{2}
.

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

נתחיל בהגדרת הפרמטרים בהם תלוי הפסד הלחץ:


\Delta p=f\left( L(length),D(diameter),e(roughness),\mu (viscosity),\bar{V}(average\text{ }velocity),\rho (density) \right)

מאנליזה ממדית על פי עקרון <math>
\pi</math> של בקינגהם נקבל:

\left. \begin{align}
 & \frac{\Delta p}{\rho }^{2}}}=f\left( \frac{\mu }{\rho \bar{V}D},\frac{L}{D},\frac{e}{D} \right) \\
 & \frac{\Delta p}{\rho }^{2}}D}=\phi \left( \operatorname{Re},\frac{L}{D},\frac{e}{D} \right) \\
 & \frac}}=\phi \left( \operatorname{Re},\frac{L}{D},\frac{e}{D} \right) \\
 & \frac}{\frac{1}{2}}^{2}}}=\phi \left( \operatorname{Re},\frac{L}{D},\frac{e}{D} \right) \\
\end{align} \right\}


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

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

\frac}{\frac{1}{2}}^{2}}}=\frac{L}{D}\phi \left( \operatorname{Re},\frac{e}{D} \right)

הפונקציה הנעלמת כאן מוגדרת כ"מקדם החיכוך" ותסומן כך:  
f=\phi \left( \operatorname{Re},\frac{e}{D} \right).

ולסיכום נרשום ביטוי פשוט יחסית: 
=f\cdot \frac{L}{D}\cdot \frac}{2}

את המקדם 
f  נמצא מתוך דיאגרמה אמפירית בשם: "דיאגרמת מודי". דיאגרמה זו (המופיעה באיור 3) היא בעצם, שרטוט הפונקציה 
f עבור ערכי Re שונים כאשר 
e/D ידוע.

ממוזערלכן, על מנת להשתמש בה עלינו לחשב את ערכו של 
e שאיננו ידוע עדיין.

את ערך 
e נחשב מתוך טבלה מספר 1:

בשימוש ב"דיאגרמת Moody" יש צורך במספר ריינולדס אותו נחשב, כאמור, ע"י: 

\operatorname{Re}=\frac{\rho UD}{\mu }=\frac{UD}{\nu }.

alt=איור 3 - דיאגרמת "מודי"|centre|ממוזער|797x797px|איור 3 - דיאגרמת "מודי"

באזור הלמינארי ההפסדים הם: 

=\left( \frac{64}{\operatorname{Re}} \right)\frac{L}{D}\frac}{2} . לכן, לפי הגדרת מקדם החיכוך נקבל כי כאשר הזרימה למינארית ערכו של 
f  פרופורציוני הפוך ל 
Re בלבד וערכו: 
f=\frac{64}{\operatorname{Re}} .

השימוש בדיאגרמה פשוט ובכל זאת בוצעו ניסיונות רבים לתרגם את הדיאגרמה למשוואה שתייתר את השימוש בדיאגרמה.
* **Colebrook** - משוואה סתומה (לא מפורשת ב-
f) : 
\frac{1}}=-2.01\cdot \log \left( \frac{e}{3.7\cdot D}+\frac{2.51}{\operatorname{Re}} \right) . המשוואה ניתנת לפתרון בפשטות בכל תוכנת חישוב בסיסית או על ידי איטרציות או על ידי פונקציית W של למברט.
* **Miller - **הציע את הניחוש ההתחלתי למשוואה האיטרטיבית בעזרתו נקבל דיוק של אחוז בודד אחרי איטרציה בודדת בלבד. 
=0.25} \right) \right]}^{-2}}
* **Blausius **- הציע ביטוי למקדם החיכוך לזרימות טורבולנטיות בצינורות חלקים תחת התנאי 
\operatorname{Re}\le : 

f=\frac{0.316}^{0.25}}}
על ידי שימוש ב: ביטוי זה של Blausius, הגדרת מקדם החיכוך ובביטוי לגזירה בדופן נוכל לקבל ביטוי פשוט למאמץ הגזירה בדופן הצינור:  

=0.0332\rho ^{2}} \right)}^{0.25}}

## הפסדים משניים – נובעים משינוי שטח החתך, פניות, שינויי גובה וכדומה
באופן כללי מחושבים ההפסדים בצורה הבאה, כאשר הקבוע 
\Kappa נקבע בצורה שונה עבור כל סיטואציה: 

=\kappa \frac}{2}

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

ברור כי ככל שהפינות תהיינה מעוגלות יותר כך ההפסדים יקטנו. במידה ובכניסה הפינות תהיינה מעוגלת היטב 
\left( \frac{r}{D}\ge 0.15 \right) ההפסד בכניסה יהיה זניח לחלוטין.

ערכו של 
\Kappa נקבע בניסויים עבור כל כניסה אך ככלל אצבע נשתמש במקדמים הבאים: (באיור 4):
alt=איור 4 - כניסה לצנרת|centre|ממוזער|365x365px|איור 4 - כניסה לצנרת

### שינויים חדים בשטח החתך
שינויים בשטח החתך מאלצים את הזרימה לשנות את מהירותה על מנת לקיים שימור מסה. הדבר גורם להפסדי לחץ.
נגדיר יחס שטחי חתך כך: 
AR=\frac}}. מקדם זה נע בתחום 0\le AR\le 1.

בעזרת הגרף (איור 5) נחשב את 
\Kappa  בצורה דומה הן להצרה והן להרחבת שטח החתך.
alt=איור 5 - שינוי חד בשטח החתך|centre|ממוזער|370x370px|איור 5 - שינוי חד בשטח החתך

### שינויים מתונים בשטח החתך
נשתמש בטבלה 2 אשר תגדיר את 
\Kappa עבור ערכי AR וזווית פתיחה שונים. הגאומטריה הרלוונטית מופיעה באיור מספר 6.ממוזער|alt=איור 6 - גאומטריה של שינוי מתון בשטח חתך|centre|379x379px|איור 6 - גאומטריה של שינוי מתון בשטח חתך

### כיפופים בצינור
על פניו, נרצה לומר שהפסדי חיכוך בצינור מעוקל דומים
לאלו של צינור ישר כל עוד נלקח אורכו המלא של העיקול. בפועל, נוצרת בעיקול זרימה משנית אשר גורמת להפסדים נוספים. לכן נמצא אורך שקול מנורמל L_e/D מתוך דיאגרמות מבוססות ניסויים ונציבו לנוסחה שקיבלנו להפסדים ראשיים: 
=f\cdot \frac}{D}\cdot \frac}{2}

עבור כיפוף רציף של 90 מעלות ברדיוסים שונים האורך השקול המנורמל ייקבע על ידי איור 6:ממוזער|centre|291x291px|איור 6- פנייה רציפה של 90 בצינורעבור פניות חדות בזוויות שונות האורך השקול המנורמל L_e/D ייקבע על ידי איור 7:

ממוזער|centre|293x293px|איור 7 - פנייה חדה בזווית כלשהי בצינור

### שסתומים ואביזרים
שסתומים סגורים מהווים מחסום עבור הזורם. בבואנו לחשב הפסדים הנובעים מחיכוך נחשב אותם עם שסתום פתוח לרווחה. שוב נשתמש בנוסחה שקיבלנו להפסדים ראשיים : 
=f\cdot \frac}{D}\cdot \frac}{2}

כאשר את האורך השקול המנורמל L_e/D נוציא מטבלה מספר 3 או מקטלוג יצרן.

## מקורות

*Fox, Robert W., and Alan T. Mcdonald. *Introduction to Fluid Mechanics*. 6th ed. N.p.: John Wilet & Sons, 2004
*https://web.archive.org/web/20150218120914/http://www.efm.leeds.ac.uk/CIVE/CIVE1400/PDF/Notes/section_all2.pdf - An Introduction to Fluid Mechanics School of Civil Engineering, University of Leeds
*http://www.codecogs.com/library/engineering/fluid_mechanics/pipes/head_loss/index.php - calculator for head loss


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

קטגוריה:מכניקת הזורמים