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

מתוך ויקיפדיה, האנציקלופדיה החופשית
קפיצה אל: ניווט, חיפוש

בהידרודינמיקה, משוואת דרסי-ויסבךאנגלית 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), שהוא רבע מערכו של מקדם דרסי.

המשוואה[עריכת קוד מקור | עריכה]

מפל לחץ[עריכת קוד מקור | עריכה]

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

\Delta P = f_D \cdot \frac{L}{D} \cdot \frac{\rho V^2}{2}

כאשר:

  • \Delta P - מפל הלחץ בצינור כתוצאה מחיכוך בדפנות (ביחידות של Pa)
  • f_D - מקדם החיכוך של דרסי – מקדם אל-ממדי (עבור זרימה למינרית או טורבולנטית)
  • \frac{L}{D} -יחס בין אורך הצינור לקוטרו ההידראולי
  • \rho - צפיפות הנוזל (ביחידות של \frac{kg}{m^3})
  • V - המהירות הממוצעת של הזורם בצינור (ביחידות של \frac{m}{s})

את המקדם f_D ניתן למצוא מדיאגרמות מודי או על ידי פתרון ה- Modified Colebrook Equation.

ניתן לעבור מתצוגה זו לתצוגת מפל העומד על ידי חלוקה ב \rho \cdot g.

מפל עומד[עריכת קוד מקור | עריכה]

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

\Delta P = \rho \cdot g \cdot h_f או h_f = \frac{\Delta P}{\rho \cdot g}

ונקבל:

h_f = f \cdot \frac{L}{D} \cdot \frac{V^2}{2g}

כאשר:

  • h_f – מפל העומד בצינור כתוצאה מחיכוך בדפנות (ביחידות של m)
  • f_D – מקדם החיכוך של דרסי (עבור זרימה למינרית או טורבולנטית)
  • L – אורך הצינור (ביחידות של m)
  • D – הקוטר ההידראולי של הצינור (עבור צינור בעל חתך עגול זהו קוטרו הפנימי, עבור צורות אחרות ניתן להשתמש בקוטר אפקטיבי) (ביחידות של m)
  • gתאוצת הכובד (ביחידות של \frac{m}{s^2})
  • V - המהירות הממוצעת של הזורם בצינור (ביחידות של \frac{m}{s})

מקדם החיכוך של דרסי[עריכת קוד מקור | עריכה]

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

עבור זרימות איטיות ולמינריות, הגודל של f_D נקבע לפי חוק פואזיל (Poiseuille's Law), ומקבלים:

f_D = \frac{64}{Re}

כאשר:

  • f_D – מקדם החיכוך של דרסי
  • Reמספר ריינולדס של הזורם, המחושב כאשר האורך האופייני הינו קוטר הצינור ההידראולי

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

הבדלה בין מקדם פאנינג (Fanning) למקדם דרסי (Darcy)[עריכת קוד מקור | עריכה]

חשוב לשים לב כי מקדם החיכוך של דרסי גדול פי 4 ממקדם החיכוך של Fanning, כלומר: f_D = 4 \cdot f כאשר אנו מעיינים בדיאגרמות מודי למציאת גודל החיכוך, חשוב שנבדוק אם הדיאגרמה הינה עבור המקדם של Darcy או Fanning. מקדם Darcy בדרך כלל משומש יותר בתחום הנדסת המכונות וההנדסה האזרחית, לעומת מקדם Fanning המשמש יותר להנדסה כימית.

נשים לב: \Delta P = f_D \cdot \frac{L}{D} \cdot 2 \cdot \rho \cdot V^2

עבור זרימה למינרית,ראינו כי f_D = \frac{64}{Re} ולכן נקבל כי מקדם Fanning הינו f = \frac{16}{Re}.

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

  1. מצא את ערך מקדם החיכוך עבור זרימה למינרית בעלת מספר ריינולדס של 1000.
  2. אם ערך המקדם הינו 0.064, ניתן להסיק שמדובר במקדם דרסי - f_D.
  3. אם ערך המקדם הינו 0.016, ניתן להסיק שמדובר במקדם Fanning – f.

היסטוריה[עריכת קוד מקור | עריכה]

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

h_f = \frac{L}{D} \cdot (aV+bV^2)

כאשר:

  • h_f – מפל העומד בצינור
  • \frac{L}{D} – יחס אורך הצינור לקוטרו
  • V – המהירות הממוצעת של הזורם בצינור
  • a,b – קבועים הנקבעים אמפירית

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

פיתוח המשוואה[עריכת קוד מקור | עריכה]

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

ללחץ יש יחידות של אנרגיה ליחידת נפח:

[\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 נקודות לאורך הצינור צריך להיות פרופורציונלי ל- \frac{1}{2} \cdot \rho \cdot V^2, ביחידות של אנרגיה קינטית ליחידת נפח. בנוסף אנו יודעים כי מפל הלחץ פרופורציונלי לאורך הצינור, כאשר מפל הלחץ ליחידת אורך הינו קבוע. על מנת להפוך קשר זה לקשר אל-ממדי, ננרמל את L – אורך הצינור, לפי D – קוטר הצינור, כאשר גם הוא קבוע לאורך הצינור.

נקבל:

\Delta P \propto \frac{L}{D} \cdot \frac {1}{2} \cdot \rho \cdot V^2

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

חשוב לציין כי הגודל \frac{1}{2} \cdot \rho \cdot V^2 אינו האנרגיה הקינטית ליחידת נפח של הזורם בצינור, עקב מספר סיבות:

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

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

יישום פרקטי[עריכת קוד מקור | עריכה]

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

V^2 = \frac{Q^2}{A_w^2} \cdot \rho \cdot V^2

כאשר:

  • V – המהירות הממוצעת של הזורם לאורך הצינור (ביחידות של \frac{m}{s})
  • Q – הספיקה הנפחית של הזורם בצינור (ביחידות של \frac{m^3}{s})
  • A_w – השטח הרטוב בצינור (ביחידות של m^2)

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

A_w^2 = ( \frac{\pi D^2}{4} )^2 = \frac{\pi^2 D^4}{16}

כאשר D הינו קוטר הצינור.

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

h_f = \frac{8 f_D L Q^2}{g \pi^2 D^5}

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

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

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

לקריאה נוספת[עריכת קוד מקור | עריכה]

  • 1982 Handbook of Hydraulics by Ernest F. Brater and Horace Williams King.Sixsth Edition. McGraw-Hill Book Company. ISBN 0-07-007243-4 פרק 6.

קישורים חיצוניים[עריכת קוד מקור | עריכה]


חישוב מפל לחץ בזרימה

משוואת דרסי ויסבךמשוואת היזן ויליאמסמשוואת מנינגמשוואת סווימי ג'ייןמשוואת קולברוקדיאגרמת מודיהלם מים