משוואת דרסי-ויסבך

משוואת דרסי-ויסבך (באנגלית :Darcy–Weisbach equation) משמשת בהידרודינמיקה לכמת את מפל העומד (Head Loss), או מפל הלחץ (Pressure Loss) בצינור, כתוצאה מחיכוך בדפנות, עבור מהירות ממוצעת של הזורם בצינור. משוואה זו נקראת על שמם של אנרי דרסי (Henry Darcy) ויוליוס ויסבך (Julius Weisbach). משוואת דרסי-ויסבך מכילה מקדם חיכוך אל-ממדי, לו קוראים מקדם החיכוך של דרסי (Darcy Friction Factor), ובנוסף ידוע גם כפקטור חיכוך דרסי-ויסבך (Darcy-Weisbach friction factor), או מקדם חיכוך מודי (Moody Friction Factor). בנוסף למקדם זה קוראים לפעמים ה-Blasius Friction Factor. אין להתבלבל בין מקדם דרסי למקדם חיכוך של פאנינג (Fanning friction factor), שהוא רבע מערכו של מקדם דרסי.

המשוואה

מקדם החיכוך של דרסי

מקדם החיכוך ${\displaystyle f_{D}}$ הוא מקדם אל-ממדי התלוי בגאומטריית הצינור, ובמהירות הזורם הנע בתוכו. גודלו ידוע בדיוק גבוה עבור משטרי זרימה מסוימים. ניתן להעריך את גודל מקדם החיכוך בעזרת מספר שיטות אמפיריות או קשרים תאורטיים, או בעזרת טבלאות וגרפים שפורסמו. הצורה הנפוצה ביותר למציאת מקדם זה הוא בעזרת דיאגרת מודי (Moody Diagram), על שמו של Lewis Ferry Moody, ומסיבה זו לפעמים מקדם החיכוך נקרא על שמו של מודי ולא של דרסי-ויסבך.

עבור זרימות איטיות ולמינריות, הגודל של ${\displaystyle f_{D}}$ נקבע לפי חוק פואזיל (Poiseuille's Law), ומקבלים:

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

כאשר:

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

עבור זרימה טורבולנטית, ניתן לחשב את ${\displaystyle f_{D}}$ מתוך דיאגרמות מודי, או מפתרונות משוואות כמו משוואת קולברוק (Colebrook-White Equation),  או משוואת סוואמי-ג'יין (Swamee-Jain Equation).

הבדלה בין מקדם פאנינג (Fanning) למקדם דרסי (Darcy)

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

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

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

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

1. מצא את ערך מקדם החיכוך עבור זרימה למינרית בעלת מספר ריינולדס של 1000.
2. אם ערך המקדם הוא ${\displaystyle 0.064}$, ניתן להסיק שמדובר במקדם דרסי - ${\displaystyle f_{D}}$.
3. אם ערך המקדם הוא עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 0.016} , ניתן להסיק שמדובר במקדם Fanning – עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle f} .

היסטוריה

משוואת דרסי-ויסבך היא וריאציה על משוואת פרוני (Prony Equation):

${\displaystyle h_{f}={\frac {L}{D}}\cdot (aV+bV^{2})}$

כאשר:

• עיבוד הנוסחה נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle h_{f}} – מפל העומד בצינור
• עיבוד הנוסחה נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle {\frac {L}{D}}} – יחס אורך הצינור לקוטרו
• ${\displaystyle V}$ – המהירות הממוצעת של הזורם בצינור
• ${\displaystyle a,b}$ – קבועים הנקבעים אמפירית

וריאציה ראשונית ומשופרת פותחה על ידי אנרי דרסי, מהנדס צרפתי שחי במאה ה-19, הידוע גם בזכות גילוי חוק דרסי עבור זרימת נוזלים במצע נקבובי, ולאחר מכן שופרה שוב לצורתה הסופית על ידי יוליוס ויסבך בשנת 1845. במקור, מידע על השתנות ${\displaystyle f_{D}}$ כתלות במהירות הייתה חסר, ולכן בשנים הראשונות שאחרי הצעתו של ויסבך לתיקון המשוואה, תוצאות של Prony עדיין סיפקו דיוק יותר טוב. בשנים שאחרי, עקב קשיים בחישוב, היה נוח לרוב להשתמש בתוצאות אמפיריות התקפות למשטרי זרימה ספציפיים, כמו משוואת היזן-ויליאמס (Hazen-Williams Equation) או משאוות מנינג (Manning Equation), שהקלו על החישוב המתמטי. אך כאשר מחשבונים נהיו נגישים לציבור הרחב, משוואת דרסי-ויסבך נהפכה לדומיננטית, כתוצאה מהטווח הרחב שעבורו היא תקפה.

פיתוח המשוואה

משוואה זו פותחה בעזרת אנליזת סדרי-גודל. באזורים רחוקים מקצוות הצינור, מתקבלת זרימה מפותחת לאורך הצינור, כלומר פרופיל הזרימה ותכונותיו לא משתנות כתלות במיקום לאורך הצינור. כתוצאה מכך, הפרמטרים החשובים לחישוב הזרימה הם מפל הלחץ ליחידת אורך בצינור - ${\displaystyle {\frac {\Delta P}{L}}}$, והספיקה הנפחית של הזורם בצינור. ניתן להמיר את הספיקה הנפחית למהירות ממוצעת על ידי חלוקתה בשטח הרטוב של הצינור, השקול לשטח החתך שבצינור עבור צינור מלא בנוזל.

ללחץ יש יחידות של אנרגיה ליחידת נפח: עיבוד הנוסחה נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle [\Delta P]={\frac {N}{m^{2}}}={\frac {N}{m^{2}}}\cdot {\frac {m}{m}}={\frac {N\cdot m}{m^{3}}}={\frac {J}{m^{3}}}} לכן, מפל המתח בין 2 נקודות לאורך הצינור צריך להיות פרופורציונלי ל- עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{1}{2} \cdot \rho \cdot V^2} , ביחידות של אנרגיה קינטית ליחידת נפח. בנוסף אנו יודעים כי מפל הלחץ פרופורציונלי לאורך הצינור, כאשר מפל הלחץ ליחידת אורך הוא קבוע. על מנת להפוך קשר זה לקשר אל-ממדי, ננרמל את עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle L}  – אורך הצינור, לפי ${\displaystyle D}$ – קוטר הצינור, כאשר גם הוא קבוע לאורך הצינור.

נקבל: עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \Delta P \propto \frac{L}{D} \cdot \frac {1}{2} \cdot \rho \cdot V^2} את קבוע הפרופורציה נגדיר כמקדם החיכוך של דרסי אותו חיפשנו - ${\displaystyle f_{D}}$. מספר אל-ממדי זה הוא פונקציה של מספר גדלים גאומטריים וקבועים נוספים, כגון ${\displaystyle \pi }$, מספר ריינולדס, ו- עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e}  - החספוס היחסי של הצינור (גובה החספוס חלקי הקוטר).

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

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

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

יישום פרקטי

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

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

כאשר:

• ${\displaystyle V}$ – המהירות הממוצעת של הזורם לאורך הצינור (ביחידות של עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{m}{s}} )
• עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle Q} – הספיקה הנפחית של הזורם בצינור (ביחידות של עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac{m^3}{s}} )
• עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle A_w} – השטח הרטוב בצינור (ביחידות של ${\displaystyle m^{2}}$)

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

${\displaystyle A_{w}^{2}=({\frac {\pi D^{2}}{4}})^{2}={\frac {\pi ^{2}D^{4}}{16}}}$

כאשר ${\displaystyle D}$ הוא קוטר הצינור.

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

עיבוד הנוסחה נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle h_f = \frac{8 f_D L Q^2}{g \pi^2 D^5}}

כאשר משמעות כל אחד הפרמטרים הוגדרה לעיל.

ניתן גם לקבל את מפל הלחץ לפי:

${\displaystyle \Delta P={\frac {8}{\pi ^{2}}}\cdot f_{D}\cdot {\frac {L}{D}}\cdot \rho \cdot {\frac {Q^{2}}{D^{4}}}}$

לקריאה נוספת

קישורים חיצוניים

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