הצגת המקור של הדף "זרימה בצנרת"

[[קובץ:Naware pipes.JPG|ממוזער|צנרת]]
'''זרימה בצנרת'''{{הערה|שם=:1|{{צ-ספר|מחבר = Fox, Robert W., and Alan T. Mcdonald|שם = Introduction to Fluid Mechanics. 6th ed.:|מו"ל = John Wilet & Sons|שנת הוצאה = 2004}}}} הוא ענף ב[[מכניקת זורמים]] אשר דן במשוואות התנועה, פרופילי המהירות והכוחות הפועלים על [[זרימה בלתי דחיסה|זורם בלתי דחיס]] כאשר הוא מבצע תנועה בצינור סגור (לא בתעלה פתוחה).

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

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

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

== משוואת ברנולי – חוק שימור האנרגיה{{הערה|שם=:1}} ==
{{ערך מורחב|משוואת ברנולי}}

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

בהזנחת החיכוך ניתן להניח שהאנרגיה של הזורם נשמרת בזמן תנועתו בצינור. על פי משוואת [[יאקוב ברנולי|ברנולי]]: <math>
P+\rho gh+\frac{1}{2}\rho {{v}^{2}}=const
</math>

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

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

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

== משטר הזרימה בצנרת:{{הערה|שם=:1}} ==
למעט זורמים צמיגים מאוד, לרוב הזרימה בצנרת היא טורבולנטית לחלוטין. המשטר ייקבע מתוך ערכו של מספר ריינולדס (מסומן Re, [[גודל חסר ממד]] אשר מייצג את היחס בין אינרציית הזורם וצמיגותו). מספר ריינולדס גבוה מעיד על זרימה בה האינרצייה שולטת בבעיה ומספר ריינולדס קטן מעיד על זרימה בה הצמיגות שולטת בבעיה.

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

== זרימה למינרית מפותחת בצינור עגול{{הערה|שם=:1}} ==
זרימה למינרית בצינור תתקיים עבור: <math>\operatorname{Re}\le 2300</math> והיא זרימה המונעת על ידי מפל לחצים מסוג [[:en:Hagen–Poiseuille_flow_from_the_Navier–Stokes_equations|זרימת פואסיל]].

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

[[קובץ:Infinitesmal_volume_flow.jpg|447x447px]]

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

מאמץ הגזירה מוגדר להיות:  <math>
{{\tau }_{rx}}=\mu \frac{\partial u}{\partial r}</math>

ופתרון המשוואה הדפרנציאלית הואו: <math>
\mu \frac{du}{dr}={{\tau }_{rx}}=\frac{r}{2}\left( \frac{\partial p}{\partial x} \right)+\frac{{{c}_{1}}}{r}</math>

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

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

== זרימה טורבולנטית מפותחת בצינור עגול ==
{{ערך מורחב|זרימה טורבולנטית}}
{{הערה|שם=:1}}

[[זרימה טורבולנטית]] תתקיים עבור: <math>\operatorname{Re}\ge 4000</math>

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

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

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

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

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

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

עבור כל צינור ניתן לחשב [[קוטר הידראולי]]: <math>{{D}_{h}}=\frac{4\cdot Area}{Perimeter}</math>. קוטר זה ישמש בתפקיד הקוטר בכל החישובים בהמשך.

== החיכוך בצנרת ==
{{הערה|שם=:1}} 

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

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

בעזרת [[משוואת ברנולי]] ניתן לתאר את הפסדי החיכוך כך:         <math>
\left( \frac{{{p}_{1}}}{\rho }+{{\alpha }_{1}}\frac{{{{\bar{V}}}_{1}}^{2}}{2}+g{{z}_{1}} \right)-\left( \frac{{{p}_{2}}}{\rho }+{{\alpha }_{2}}\frac{{{{\bar{V}}}_{2}}^{2}}{2}+g{{z}_{2}} \right)={{h}_{{{l}_{T}}}}</math>

כאשר המהירויות שוות ותכונות הזורם נשארות קבועות ניתן לקבל את התוצאה המופשטת: <math>
\frac{{{p}_{1}}-{{p}_{2}}}{\rho }=g({{z}_{2}}-{{z}_{1}})+{{h}_{l}}</math>

בנוסף, אם הצינור אופקי, אזי הגובה איננו משתנה ונקבל: <math>
\frac{{{p}_{1}}-{{p}_{2}}}{\rho }=\frac{\Delta p}{\rho }={{h}_{l}}</math>

== הפסדים ראשיים – הנובעים מזרימה מתמדת בצינור אופקי ללא שינוי שטח החתך ==
==== הפסדי לחץ הנובעים מחיכוך בזרימה למינרית ====
{{הערה|שם=:2|{{קישור כללי|כתובת = http://www.codecogs.com/library/engineering/fluid_mechanics/pipes/head_loss/index.php|כותרת = code coges}}}}{{הערה|שם=:1}}

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

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

\Delta p=32\frac{L}{D}\cdot \frac{\mu \bar{V}}{D}
</math>

ולכן: <math>

{{h}_{l}}=\frac{\Delta p}{\rho }=\left( \frac{64}{\operatorname{Re}} \right)\frac{L}{D}\frac{{{{\bar{V}}}^{2}}}{2}
</math>.

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

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

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

מאנליזה ממדית על פי [[:en:Buckingham_π_theorem|עקרון <math>
\pi</math> של בקינגהם]] נקבל:<math>

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

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

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

\frac{{{h}_{l}}}{\frac{1}{2}{{{\bar{V}}}^{2}}}=\frac{L}{D}\phi \left( \operatorname{Re},\frac{e}{D} \right)</math>

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

ולסיכום נרשום ביטוי פשוט יחסית: <math>
{{h}_{l}}=f\cdot \frac{L}{D}\cdot \frac{{{{\bar{V}}}^{2}}}{2}</math>

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

[[קובץ:Moody EN.svg|ממוזער]]לכן, על מנת להשתמש בה עלינו לחשב את ערכו של <math>
e</math> שאיננו ידוע עדיין.

את ערך <math>
e</math> נחשב מתוך טבלה מספר 1:
{| class="wikitable"
|+טבלה מספר 1 - Roughness{{הערה|שם=:1}}
!
{| class="MsoTableGrid"

 |
 <math>
e [mm]</math>

|}
!
{| class="MsoTableGrid"

 |
'''Pipe Type'''

|}
|-
|0.9 - 9
|Riveted Steel
|-
|0.3 - 3
|Concrete
|-
|0.2 - 0.9
|Wood Stave
|-
|0.26
|Cast Iron
|-
|0.15
|Galvanized Iron
|-
|0.046
|Commercial Steel or Wrought Iron
|-
|0.0015
|Drawn Tubing
|}

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

\operatorname{Re}=\frac{\rho UD}{\mu }=\frac{UD}{\nu }</math>.

[[קובץ:Moody EN.svg|alt=איור 3 - דיאגרמת "מודי"|centre|ממוזער|797x797px|איור 3 - דיאגרמת "מודי"]]

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

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

השימוש בדיאגרמה פשוט ובכל זאת בוצעו ניסיונות רבים לתרגם את הדיאגרמה למשוואה שתייתר את השימוש בדיאגרמה.
* '''Colebrook''' - משוואה סתומה (לא מפורשת ב-<math>
f</math>) : <math>
\frac{1}{{{f}^{0.5}}}=-2.01\cdot \log \left( \frac{e}{3.7\cdot D}+\frac{2.51}{\operatorname{Re}{{f}^{0.5}}} \right)</math> . המשוואה ניתנת לפתרון בפשטות בכל תוכנת חישוב בסיסית או על ידי איטרציות או על ידי [[פונקציית W של למברט]]{{אנ|Lambert W function}}.
* '''Miller - '''הציע את הניחוש ההתחלתי למשוואה האיטרטיבית בעזרתו נקבל דיוק של אחוז בודד אחרי איטרציה בודדת בלבד. <math>
{{f}_{o}}=0.25{{\left[ \log \left( \frac{e}{3.7\cdot D}+\frac{5.74}{{{\operatorname{Re}}^{0.9}}} \right) \right]}^{-2}}</math>
* '''Blausius '''- הציע ביטוי למקדם החיכוך לזרימות טורבולנטיות בצינורות חלקים תחת התנאי <math>
\operatorname{Re}\le {{10}^{5}}</math>: <math>

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

{{\tau }_{w}}=0.0332\rho {{\bar{V}}^{2}}{{\left( \frac{\nu }{R\bar{V}} \right)}^{0.25}}</math>

== הפסדים משניים – נובעים משינוי שטח החתך, פניות, שינויי גובה וכדומה{{הערה|שם=:2}}{{הערה|שם=:1}} ==
באופן כללי מחושבים ההפסדים בצורה הבאה, כאשר הקבוע <math>
\Kappa</math> נקבע בצורה שונה עבור כל סיטואציה: <math>

{{h}_{{{l}_{m}}}}=\kappa \frac{{{{\bar{V}}}^{2}}}{2}</math>

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

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

ערכו של <math>
\Kappa</math> נקבע בניסויים עבור כל כניסה אך כ[[כלל אצבע]] נשתמש במקדמים הבאים: (באיור 4{{הערה|שם=:1}}):
[[קובץ:PIPE ENTRANCE FRICTION LOSS COEF.JPG|alt=איור 4 - כניסה לצנרת|centre|ממוזער|365x365px|איור 4 - כניסה לצנרת]]

=== שינויים חדים בשטח החתך ===
שינויים בשטח החתך מאלצים את הזרימה לשנות את מהירותה על מנת לקיים שימור מסה. הדבר גורם להפסדי לחץ.
נגדיר יחס שטחי חתך כך: <math>
AR=\frac{{{A}_{small}}}{{{A}_{big}}}</math>. מקדם זה נע בתחום <math>0\le AR\le 1</math>.

בעזרת הגרף (איור 5{{הערה|שם=:1}}) נחשב את <math>
\Kappa</math>  בצורה דומה הן להצרה והן להרחבת שטח החתך.
[[קובץ:AREA CHANGE FRICTION LOSS COEF.png|alt=איור 5 - שינוי חד בשטח החתך|centre|ממוזער|370x370px|איור 5 - שינוי חד בשטח החתך]]

=== שינויים מתונים בשטח החתך ===
נשתמש בטבלה 2 אשר תגדיר את <math>
\Kappa</math> עבור ערכי AR וזווית פתיחה שונים. הגאומטריה הרלוונטית מופיעה באיור מספר 6.[[קובץ:ENTRANCE TO PIPE WITH OPEN ANGLE.png|ממוזער|alt=איור 6 - גאומטריה של שינוי מתון בשטח חתך|centre|379x379px|איור 6 - גאומטריה של שינוי מתון בשטח חתך]]
{| class="wikitable"
|+טבלה 2 - מקדם <math>
\Kappa</math> לשינוי מתון בשטח החתך{{הערה|שם=:1}}
! colspan="7" |Included Angle <math>
\theta
</math>, Degrees
|-
|150
|120
|90
|50-60
|15-40
|10
|<math>
A2/A1
</math>
|-
|0.24
|0.18
|0.12
|0.06
|0.05
|0.05
|0.50
|-
|0.35
|0.27
|0.17
|0.07
|0.04
|0.05
|0.25
|-
|0.37
|0.29
|0.19
|0.08
|0.05
|0.05
|0.10
|}

=== כיפופים בצינור ===
על פניו, נרצה לומר שהפסדי חיכוך בצינור מעוקל דומים
לאלו של צינור ישר כל עוד נלקח אורכו המלא של העיקול. בפועל, נוצרת בעיקול זרימה משנית אשר גורמת להפסדים נוספים. לכן נמצא אורך שקול מנורמל <math>L_e/D</math> מתוך דיאגרמות מבוססות ניסויים ונציבו לנוסחה שקיבלנו להפסדים ראשיים: <math>
{{h}_{l}}=f\cdot \frac{{{L}_{eq}}}{D}\cdot \frac{{{{\bar{V}}}^{2}}}{2}</math>

עבור כיפוף רציף של 90 מעלות ברדיוסים שונים האורך השקול המנורמל ייקבע על ידי איור 6{{הערה|שם=:1}}:[[קובץ:TURNING PIPE FRICTION LOSS COEF.png|ממוזער|centre|291x291px|איור 6- פנייה רציפה של 90 בצינור]]עבור פניות חדות בזוויות שונות האורך השקול המנורמל <math>L_e/D</math> ייקבע על ידי איור 7{{הערה|שם=:1}}:

[[קובץ:PIPE WITH SHARP ANGLE FRICTION LOSS COEF.png|ממוזער|centre|293x293px|איור 7 - פנייה חדה בזווית כלשהי בצינור]]

=== שסתומים ואביזרים ===
שסתומים סגורים מהווים מחסום עבור הזורם. בבואנו לחשב הפסדים הנובעים מחיכוך נחשב אותם עם שסתום פתוח לרווחה. שוב נשתמש בנוסחה שקיבלנו להפסדים ראשיים : <math>
{{h}_{l}}=f\cdot \frac{{{L}_{eq}}}{D}\cdot \frac{{{{\bar{V}}}^{2}}}{2}</math>

כאשר את האורך השקול המנורמל <math>L_e/D</math> נוציא מטבלה מספר 3 או מקטלוג יצרן.
{| class="wikitable"
|+{{הערה|שם=:1}}'''טבלה 3 - אורך שקול מנורמל לחישוב הפסדי לחץ באביזרים שונים'''
!<math>L_e/D</math> Equivalent Legth
! colspan="2" |Fitting Type
|-
|
| colspan="2" |'''Valves'''
|-
|8
|Gate Valve
|-
|340
|Globe Valve
|-
|150
|Angle Valve
|-
|3
|Ball Valve
|-
|600
|Globe Lift Check Valve
|-
|55
|Angle Lift Check Valve
|-
|30
| colspan="2" |'''Standart Elbow 90°'''
|-
|16
| colspan="2" |'''Standart Elbow 45°'''
|-
|50
| colspan="2" |'''Return Bend, Close Pattern'''
|}

== מקורות ==
<div class="mw-content-ltr">
*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
</div>

==הערות שוליים==
{{הערות שוליים|יישור=ימין}}


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