משוואת האגן-פואזיי

משוואת האגן-פואזיי (Hagen–Poiseuille) או בפשטות "משוואת פואזיי" היא משוואה פיזיקלית-הנדסית במכניקת זורמים, המתארת את ירידת הלחץ בנוזל, הזורם דרך צינור גלילי אופקי. המשוואה התקבלה בנפרד עבור פואזיי בשנת 1838 ועבור האגן בשנת 1839, אך התפרסמה על ידי פואזיי בשנת 1840. המשוואה שימושית לאומדן הפרשי לחצים בצינורות בעלי קוטר קטן: כלי דם, קש שתייה, מחטי מזרקים וכדומה. ההנחות לקיום המשוואה הן: זרימה בלתי דחיסה זורם ניוטוני זרימה שכבתית הזרימה היא דרך צינור בעל אורך שגדול באופן משמעותי מקוטר הצינור. תאוצת הזורם שווה לאפס (כל אחת משכבות הזורם - זורמת במהירות קבועה). במקרים בהם ההנחות לא מתקיימות, למשל: אורך צינור קצר מידי, זרימה טורבולנטית ולא למינארית, וכדומה, ירידת הלחץ תהיה שונה משמעותית מירידת הלחץ המחושבת על ידי המשוואה. ירידת הלחץ מחושבת על ידי הנוסחה: Δ P = 8 μ L Q π r 4 = 8 π μ L Q A 2 {\displaystyle \Delta P={\frac {8\mu LQ}{\pi r^{4}}}={\frac {8\pi \mu LQ}{A^{2}}}} כאשר: Δ P {\displaystyle \Delta P} - הפרש הלחץ (ירידת הלחץ) L {\displaystyle L} - אורך הצינור μ {\displaystyle \mu } - צמיגות הזורם Q {\displaystyle Q} - הספיקה הנפחית r {\displaystyle r} - רדיוס הצינור A {\displaystyle A} - שטח חתך הזרימה המשוואה כושלת במדידת הפרש הלחץ בתנאים שבאופן טבעי יוצרים זרימה טורבולונטית. עבור זורם בעל צמיגות נמוכה, צנרת קצרה או רחבה, נדרש שימוש במודלים מורכבים יותר כמו משוואת דרסי-וויסבאך. כמו כן, לא ניתן למדוד את הפרש הלחץ בסמוך לתחילת הצינור, בשל תופעות מעבר הגורמות להפסדי לחץ נוספים. נתייחס לצינור בצורת גליל (תמונה מספר 1) ונסתכל על שכבת זורם קטנה באורך L. הכוחות הפועלים הם: כוחות לחץ וכוחות גזירה, כמסומן. לחץ: על חתך הטבעת בנקודה Z = 0 {\displaystyle Z=0} : ( 2 π r Δ r ) P 0 {\displaystyle (2\pi r\Delta r)P_{0}} ועל חתך הטבעת בנקודה Z = L {\displaystyle Z=L} : − ( 2 π r Δ r ) P L {\displaystyle -(2\pi r\Delta r)P_{L}} גזירה: כוח הגזירה הפועל על הטבעת הפנימית: ( 2 π r L τ r ) | r {\displaystyle (2\pi rL\tau _{r})|_{r}} כוח הגזירה שהטבעת החיצונית מפעילה על שכבת הנוזל הבאה איתה במגע: − ( 2 π r L τ r ) | r + Δ r {\displaystyle -(2\pi rL\tau _{r})|_{r+\Delta r}} חיבור הכוחות והשוואה לאפס: ( 2 π r Δ r ) P 0 − ( 2 π r Δ r ) P L + ( 2 π r L τ r ) | r − ( 2 π r L τ r ) | r + Δ r = 0 {\displaystyle (2\pi r\Delta r)P_{0}-(2\pi r\Delta r)P_{L}+(2\pi rL\tau _{r})|_{r}-(2\pi rL\tau _{r})|_{r+\Delta r}=0} מכאן, נקבל כי: r τ r | r + Δ r − r τ r | r = r Δ r Δ P L {\displaystyle r\tau _{r}|_{r+\Delta r}-r\tau _{r}|_{r}={\frac {r\Delta r\Delta P}{L}}} כלומר: d ( r τ ) d r Δ r = r Δ r Δ P L {\displaystyle {\frac {d(r\tau )}{dr}}\Delta r={\frac {r\Delta r\Delta P}{L}}} נבצע אינטגרל ונקבל: r τ = r 2 Δ P 2 L + C 1 {\displaystyle r\tau ={\frac {r^{2}\Delta P}{2L}}+C_{1}} ניתן להזניח את הקבוע ( C 1 = 0 {\displaystyle C_{1}=0} ), ולקבל: τ = r Δ P 2 L {\displaystyle \tau ={\frac {r\Delta P}{2L}}} ( ∗ ) {\displaystyle (*)} עבור זורם ניוטוני (הנחה מספר 2): הנוסחה עבור מאמץ הגזירה בין שני לוחות: τ = − μ d v d x {\displaystyle \tau =-\mu {\frac {dv}{dx}}} ובהתאמה, עבור צינור או גליל: τ = − μ d v d r {\displaystyle \tau =-\mu {\frac {dv}{dr}}} נשווה למשוואה ( ∗ ) {\displaystyle (*)} : τ = − μ d v d r = r Δ P 2 L {\displaystyle \tau =-\mu {\frac {dv}{dr}}={\frac {r\Delta P}{2L}}} נכפול ב- − d r μ {\displaystyle -{\frac {dr}{\mu }}} ונקבל: d v = − Δ P 2 L μ r d r {\displaystyle dv=-{\frac {\Delta P}{2L\mu }}rdr} נבצע אינטגרל ונקבל: v = − Δ P 4 L μ r 2 + C 2 {\displaystyle v=-{\frac {\Delta P}{4L\mu }}r^{2}+C_{2}} ( ∗ ∗ ) {\displaystyle (**)} כדי למצוא את הקבוע C 2 {\displaystyle C_{2}} , נשתמש בתנאי השפה: v ( r = R ) = 0 {\displaystyle v_{(r=R)}=0} - כלומר, על שפת הגליל, המהירות הזורם זהה לצינור הנייח, משמע- אפס. נציב את תנאי השפה במשוואה ( ∗ ∗ ) {\displaystyle (**)} ונחלץ את C 2 {\displaystyle C_{2}} : C 2 = Δ P R 2 4 L μ {\displaystyle C_{2}={\frac {\Delta PR^{2}}{4L\mu }}} כעת, נציב את C 2 {\displaystyle C_{2}} חזרה למשוואה ( ∗ ∗ ) {\displaystyle (**)} ונקבל: v r = Δ P 4 L μ ( R 2 − r 2 ) {\displaystyle v_{r}={\frac {\Delta P}{4L\mu }}(R^{2}-r^{2})} ( ∗ ∗ ∗ ) {\displaystyle (***)} קיבלנו ביטוי פרבולי למהירות וניתן לראות כי הוא בהחלט תואם למצופה מזרימת פואזיי (תמונה מספר 2). ניתן לראות מביטוי ( ∗ ∗ ∗ ) {\displaystyle (***)} כי המהירות המקסימלית מתקבלת במרכז, עבור r = 0 {\displaystyle r=0} . ניתן לוודא כי זוהי המהירות המקסימלית על ידי גזירת המשוואה: d v d r = − r Δ P 2 L μ {\displaystyle {\frac {dv}{dr}}=-{\frac {r\Delta P}{2L\mu }}} וכאמור, עבור r = 0 {\displaystyle r=0} , d v d r = 0 {\displaystyle {\frac {dv}{dr}}=0} - הנגזרת מתאפסת- והמשמעות היא שאכן המהירות המקסימלית נמצאת בנקודה זו. חישוב ספיקה נפחית לנוזל ניוטוני: כל אחת מטבעות הזורם בצינור נעה במהירות v r {\displaystyle v_{r}} . כדי לחשב את הספיקה עלינו לחשב את הכפל בין שטח כל טבעת ומהירות הטבעת, ולסכום עבור כל הטבעות. שטח של טבעת אינפיניטסימלית: d A = π ( r + d r ) 2 − π r 2 = 2 π r d r + ( d r ) 2 {\displaystyle dA=\pi (r+dr)^{2}-\pi r^{2}=2\pi rdr+(dr)^{2}} וכיוון שהביטוי ( d r ) 2 {\displaystyle (dr)^{2}} קטן מאוד- נוכל להזניח אותו, ולקבל: d A = 2 π r d r {\displaystyle dA=2\pi rdr} כעת נחשב את הספיקה: Q = ∫ d q = ∫ v r d A = ∫ v 2 π r d r {\displaystyle Q=\int dq=\int v_{r}dA=\int v2\pi rdr} נציב את משוואה ( ∗ ∗ ∗ ) {\displaystyle (***)} נפתור את האיטגרל בכל תחום הגליל, כלומר החל מ- r = 0 {\displaystyle r=0} ועד r = R {\displaystyle r=R} : Q = ∫ 0 R Δ P 4 L μ ( R 2 − r 2 ) 2 π r d r = Δ P 4 L μ 2 π [ R 2 r 2 2 − r 4 4 ] 0 R {\displaystyle Q=\int \limits _{0}^{R}{\frac {\Delta P}{4L\mu }}(R^{2}-r^{2})2\pi rdr={\frac {\Delta P}{4L\mu }}2\pi {\bigg [}{\frac {R^{2}r^{2}}{2}}-{\frac {r^{4}}{4}}{\bigg ]}_{0}^{R}} מכאן: Q = π Δ P R 4 8 L μ {\displaystyle Q={\frac {\pi \Delta PR^{4}}{8L\mu }}} ואם נעביר אגפים, נקבל את "משוואת האגן-פואזיי": Δ P = 8 μ L Q π r 4 {\displaystyle \Delta P={\frac {8{\boldsymbol {\mu }}LQ}{\pi r^{4}}}} בפיזיקה, משוואת וושבורן מתארת זרימה קפילרית בצינור. משוואה זו נכתבה על ידי Edward Wight Washburn. L 2 = γ D t 4 η {\displaystyle L^{2}={\frac {\gamma Dt}{4\eta }}} כאשר t מייצג את הזמן האופייני לנוזל בעל צמיגות η {\displaystyle \eta } ומתח פנים γ {\displaystyle \gamma } כדי לחדור מרחק L {\displaystyle L} לתוך קפילרה בעלת קוטר D {\displaystyle D} . במאמרו משנת 1921, משתמש וושבורן במשוואת פואזיי. שימוש בגורם אינטגרציה הנותן ביטוי לשינוי נפח אינפיניטסימלי כתלות ב- l {\displaystyle l} , אורך הנוזל בצינור, d V = π r 2 d l {\displaystyle dV=\pi r^{2}dl} נקבל: δ l δ t = ∑ P 8 r 2 η l ( r 4 + 4 ϵ r 3 ) {\displaystyle {\frac {\delta l}{\delta t}}={\frac {\sum P}{8r^{2}\eta l}}(r^{4}+4\epsilon r^{3})} כאשר ∑ P {\displaystyle \sum P} הוא סכום הלחצים במערכת כגון לחץ אטמוספירי P A {\displaystyle P_{A}} , הלחץ ההידרוסטטי P h {\displaystyle P_{h}} והכוח האקוויוילנטי הנוצר בקפילרות P c {\displaystyle P_{c}} . ϵ {\displaystyle \epsilon } הוא מקדם שיפוע המוערך כ-0 עבור נקודת מפגש בין נוזל למוצק הידוע בשם wetting. את הלחץ ניתן לבטא באמצעות: P h = h g ρ − l g ρ sin ⁡ ψ {\displaystyle P_{h}=hg\rho -lg\rho \sin \psi } P c = 2 γ r cos ⁡ ϕ {\displaystyle P_{c}={\frac {2\gamma }{r}}\cos \phi } כאשר ρ {\displaystyle \rho } מבטא את צפיפות הנוזל ו- γ {\displaystyle \gamma } את מתח הפנים שלו. ψ {\displaystyle \psi } זווית הצינור ביחס לציר האופקי ו- ϕ {\displaystyle \phi } היא זווית המגע בין הנוזל והחומר הקאפילרי. לאחר החלפת משתנים נקבל: δ l δ t = [ P A + g ρ ( h − l sin ⁡ ψ ) + 2 γ r cos ⁡ ϕ ] ( r 4 + 4 ϵ r 3 ) 8 r 2 η l {\displaystyle {\frac {\delta l}{\delta t}}={\frac {[P_{A}+g\rho (h-l\sin \psi )+{\frac {2\gamma }{r}}\cos \phi ](r^{4}+4\epsilon r^{3})}{8r^{2}\eta l}}} עבור זורם דחיס בצינור, הספיקה הנפחית ומהירות הזורם אינם אחידים לאורך הצינור. הזרימה לרוב, באה לידי ביטוי על ידי הלחץ ביציאה. כאשר הנוזל נדחס או מתרחב, מתבצעת עבודה והנוזל מתחמם או מתקרר בהתאמה. מתקבל מכך כי הספיקה תלויה במעבר החום אל ומתוך הזורם. עבור נקודת מבט איזותרמית של גז אידיאלי, כאשר טמפרטורת הזורם תלויה במאזנה עם הסביבה וכאשר הלחץ בין קצוות הצינור קטן, הספיקה הנפחית ביציאה מהצינור נתונה על ידי: כאמור, בשנת 1839 פרסם האגן מאמר בנושא זרימת מים בצינורות גליליים. תוצאות הניסויים המוזכרים במאמר דומים לאלו שהציע פואזיי, אך פחות מקיפות ופחות מדויקות מתוצאותיו של פואזיי. תוצאות הניסויים אולם, הכילו הבחנות על אפקטי קצה והבחנות על השוני שבין זרימה למינארית לטורבולנטית. התוצאה שהציע האגן להפרש הלחצים היא: P=ׂ(A*L*Qּּ+B*Q²ׁ)/Dֶ^4 כאשר A,B קבועים. האגן מצא כי הקבוע A תלוי בטמפרטורה והביע אותו באמצעות הנוסחא: A=a-b*T+c*T² כאשר a,b,c קבועים המתקבלים מניסוי. עבור ספיקות קטנות מאוד, הערך Q² ניתן להזנחה, כך שמתקבלת למעשה אותה הנוסחא שהוצעה להפרש הלחצים על ידי פואזיי. Prandtl ו Tietjens המירו את התוצאות שקיבל האגן עבור הקבוע A, על מנת לקבל גרף עבור פקטור חיכוך אל מול מספר ריינולדס. תוצאות אלו ניבאו היטב את תוצאות הקו התאורטי עבור מספרי ריינולדס בטווח 70-1000: f=64/RN כאשר f מייצג את פקטור החיכוך המקובל בצינור. בדרך כלל משוואת האגן-פואזיי מצביעה לא רק על הפרש הלחצים בזרימה בציור אלא גם על צורתה הפרבולית של הזרימה. למעשה ניתן לפתח את תוצאות הפרש הלחצים גם לזרימה טורבולנטית על ידי קביעת צמיגות טורבולנטית אפקטיבית במקרה של זרימה טורבולנטית, אפילו שצורת הזרימה הטורבולנטית היא לא לגמרי פרבולית. בשני המקרים, הירידה בלחץ בצינור קשורה לעומס על הדפנות, אשר כביכול קובעת את פקטור החיכוך. ניתן לקבוע את ירידת הלחץ בצינור כתוצאה מחיכוך על ידי משוואת דרסי-ויסבך. מדיה וקבצים בנושא משוואת האגן-פואזיי בוויקישיתוף ההיסטוריה של משוואת האגן-פואזיי, באנגלית משוואת דרסי-ויסבך זרימה דחיסה במספרי מאך נמוכים