הנוסחה הברומטרית

בתרמודינמיקה, הנוסחה הברומטרית (Barometric formula) היא נוסחה פיזיקלית המתארת את הלחץ האטמוספירי כתלות בגובה בשדה כבידה, תחת הנחות המודל של גז אידיאלי.

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

רקע היסטורי


מימין פייר סימון לפלס ומשמאל אדמונד האלי

העובדה שלחץ האוויר דועך ככל שעולים בגובה הייתה כבר ידועה היטב לפיזיקאים במאה ה-17, אשר ידעו גם להסביר כי הסיבה לכך היא שככל שמגביהים יש פחות אוויר מעליך הנלחץ מטה על ידי כוח הכובד. אך תיאור מתמטי של תלות הלחץ בגובה התאפשר לראשונה רק לאחר פרסומו בשנת 1662 של חוק בויל הקובע יחס הפוך בין הלחץ לבין הנפח של גז אידיאלי.^[1]

הנוסחה הברומטרית פורסמה לראשונה בשנת 1686 על ידי הפיזיקאי אדמונד האלי, אולם הנוסחה פורסמה ללא פיתוח מתמטי מסודר. הראשון לפרסם פיתוח מסודר של הנוסחה הברומטרית היה המתמטיקאי והפיזיקאי פייר-סימון לפלס בשנת 1805, בכרך הרביעי של עבודתו בת חמשת הכרכים מכניקה שמיימית (אנ') למעלה ממאה שנים לאחר שהנוסחה פורסמה לראשונה על ידי האלי.

הנוסחה הברומטרית

מודל הגז האידיאלי מתאר גזים טוב יותר ככל שהטמפרטורה גבוהה יותר, והלחץ נמוך יותר, כך שהאנרגיה הקינטית של החלקיקים משמעותית יותר מהאינטראקציה הבין חלקיקית^[2]. אטמוספירת כדור הארץ מורכבת בעיקר ממולקולות דו-אטומיות של חנקן וחמצן, ובטמפרטורה ולחץ סטנדרטים ( $T=300K,p=1atm$ ) מתוארת היטב בתור גז אידיאלי. לעומת זאת, האטמוספירה של נוגה, כוכב הלכת השני במרחקו מהשמש, דחוסה מאוד. הלחץ על פני הכוכב הוא בערך פי 90 מהלחץ על פני כדור הארץ, ולכן לא ניתן להשתמש במודל הגז האידיאלי ובנוסחה הברומטרית על פני כוכב נוגה. ההנחה שהאוויר מתנהג כמו גז אידיאלי טובה בעיקר עבור אוויר יבש, כאשר ישנם אדי מים באוויר ההנחה נפגעת עקב אינטראקציה חזקה יחסית של מולקולות המים.

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

הנוסחה הברומטרית האיזותרמית

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

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

$p=p_{0}e^{-{\frac {mgh}{k_{B}T}}}$

כאשר

$h$	הגובה מעל נקודת היחוס שהוגדרה בתור $h=0$ (לרוב זהו גובה פני הים)
$p$	הלחץ בגובה $h$ מעל נקודת היחוס
$p_{0}$	הלחץ בגובה נקודת היחוס ( $h=0$ )
$m$	מסתה הממוצעת של מולקולה באטמוספירה
$k_{B}$	קבוע בולצמן
$T$	הטמפרטורה ביחידות קלווין
$g$	תאוצת הכובד

כאשר $p_{0}$ , $p$ צריכים להמדד באותן יחידות של לחץ (לדוגמה פסקל), ויש להתאים את קבוע בולצמן ותאוצת הכובד ליחידות בהן מודדים את הגובה ומסת החלקיקים.

הנוסחה הברומטרית בטמפרטורה משתנה

גרף המתאר את הטמפרטורה, הלחץ, הצפיפות, ומהירות הקול כתלות בגובה לפי ה- U.S Standard Atmosphere.

אם מניחים שישנה תלות ליניארית של הטמפרטורה בגובה מהצורה:

$T=T_{0}-ah$

כאשר $a$ קבוע ביחידות של קלווין ליחידת אורך, ו- $T_{0}$ הטמפרטורה בגובה שהוגדר להיות $h=0$

מתקבלת הנוסחה הברומטרית:

$p=p_{0}\left(1-{\frac {a}{T_{0}}}h\right)^{\frac {mg}{ak_{B}}}$

כאשר כל הקבועים הוגדרו בטבלה לעיל, למעט $a,T_{0}$ שהוגדרו בפסקה הקודמת.

עבור אטמוספירת כדור הארץ ההנחה שהטמפרטורה ליניארית בגובה שימושית במקטעים רחבים יחסית, אך לא לאורך כל האטמוספירה ^[3]. מדידות נסיוניות מראות כי התלות הכללית של הטמפרטורה בגובה מאוד לא טריוויאלית, שינוי הטמפרטורה אפילו אינו מונוטוני כתלות בגובה. למרות זאת, ניתן לאמר בקירוב טוב שהתלות היא ליניארית למקוטעין.

באופן סכמטי, הטמפרטורות יורדות ככל שעולים בגובה באזור הטרופוספירה (0-10 ק"מ מעל פני הים). לאחר מכן, באזור הסטרטוספירה (10-50 ק"מ גובה), הטמפרטורות עולות ככל שמגביהים. תופעה זו נובעת בין היתר בגלל בליעה של קרני השמש בטווח העל-סגול בשכבת האוזון הנמצאת בסטרטוספירה. לאחר מכן בגבהים של מעל חמישים קילומטרים מגובה פני הים, במזוספירה, הטמפרטורות יורדות שוב עם הגובה וכן הלאה. למרות התלות הכללית המסובכת, כפי שנכתב לעיל, ממדידות גם ידוע כי בכל אזור בנפרד הטמפרטורה בקירוב טוב ליניארית בגובה. כך למשל מפני הים עד לגובה של 10 ק"מ הטמפרטורה צונחת בערך ב-6.5 מעלות קלווין לכל קילומטר, ולאחר מכן נשארת בערך קבועה עד גובה של 20 ק"מ.

פיתוח הנוסחאות

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

פיתוח משיקולים הידרוסטטים

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

נתבונן בעמוד גז (למשל האוויר באטמוספירה) הנמצא במנוחה, גובהו של עמוד הגז קטן מאוד, ונסמנו $dh$ ושטח החתך שלו $A$ .לכל חלקיק גז מסה $m$ , וצפיפות החלקיקים $n(h)$ וגם היא יכולה להיות תלויה בגובה, אך מכוון שגובה העמוד, $dh$ , קטן מאוד הצפיפות כמעט ולא משתנה לאורכו ולכן מספר חלקיקי הגז הכלואים בו הוא הצפיפות מוכפלת בנפח העמוד $N=n(h)Adh$

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

$-Nmg$	כוח המשיכה הפועל על כל החלקיקים כלפי מטה
$A*p(h)$	הלחץ בתחתית העמוד המפעיל עליו כוח כלפי מעלה
$-A*p(h+dh)$	הלחץ בחלקו העליון של העמוד שמפעיל עליו כוח כלפי מטה

נזכור שמספר החלקיקים בעמוד הגז הוא $N=n(h)Adh$ ולכן משוואת הכוחות היא:

$-n(h)Amgdh+p(h)A-p(h+dh)A=0$

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

${\frac {dp}{dh}}=-mgn(h)$

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

$p(h)=n(h)k_{B}T(h)$

נציב את הצפיפות ממשוואת המצב במשוואה ההידרוסטטית ונקבל משוואה דיפרנציאלית על הלחץ

${\frac {dp}{dh}}=-{\frac {mg}{k_{b}T(h)}}p(h)$

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

פיתוח הנוסחה האיזותרמית

הגענו למשוואה דיפרנציאלית על הלחץ התלויה בטמפרטורה. אם נניח שהטמפרטורה קבועה כשעולים בגובה $T(h)=T_{0}$ ניתן מיד לפתור את המשוואה באמצעות הפרדת משתנים:

${\frac {dp}{p}}=-{\frac {mg}{k_{b}T_{0}}}dh$

אינטגרציה מ- $h=0$ עד $h$ תניב: (בלוגריתם אין צורך בערך מוחלט מכוון שהלחץ תמיד חיובי)

$ln\left({\frac {p(h)}{p_{0}}}\right)=-{\frac {mgh}{k_{B}T_{0}}}$

ניקח אקספוננט של שני האגפים ונקבל:

$p=p_{0}e^{-{\frac {mgh}{k_{B}T_{0}}}}$

וזו הנוסחה הברומטרית האיזותרמית.

הנוסחה עבור טמפרטורה משתנה

נחזור חזרה לנוסחה שפיתחנו משיקולים הידרוסטטים:

${\frac {dp}{dh}}=-{\frac {mg}{k_{b}T(h)}}p(h)$

כעת נניח את התלות $T(h)=T_{0}-ah$ , כאשר $a$ קבוע ביחידות של קלווין ליחידת אורך.

הפרדת משתנים תניב את המשוואה:

${\frac {dp}{p}}=-{\frac {mg}{k_{B}(T_{0}-ah)}}dh$

אינטגרציה מ- $h=0$ עד $h$ תניב: (בלוגריתם אין צורך בערך מוחלט מכוון שהלחץ תמיד חיובי)

$ln\left({\frac {p(h)}{p_{0}}}\right)={\frac {mg}{k_{B}a}}ln\left({\frac {T_{0}-ah}{T_{0}}}\right)$

בעזרת חוקי הלוגריתם ניתן לרשום:

$ln\left({\frac {p(h)}{p_{0}}}\right)=ln\left[\left(1-{\frac {a}{T_{0}}}h\right)^{\frac {mg}{k_{B}a}}\right]$

ניקח אקספוננט של שני האגפים ונקבל:

$p(h)=p_{0}\left(1-{\frac {a}{T_{0}}}h\right)^{\frac {mg}{ak_{B}}}$

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

פיתוח באמצעות דרישה לשיווי משקל

על מנת להבין הוכחה זו דרוש ידע בסיסי בתרמודינמיקה.

ראשית נוכיח טענת עזר כללית על שתי מערכות הנמצאות בשיווי משקל במגע דיפוזיבי, כלומר חלקיקים יכולים לעבור בין המערכות. ידוע ששיווי משקל עבור מערכת בטמפרטורה קבועה משמעותו מינימום של האנרגיה החופשית של הלמהולץ. נסתכל על שתי מערכות A,B אשר נמצאות באותה טמפרטורה. מכוון שהטמפרטורה והנפח לא משתנים ניתן לרשום את הדיפרנציאל של האנרגיה החופשית הכוללת של שתי המערכות $F_{tot}=F_{A}+F_{B}$ באופן הבא:

$dF_{tot}=\left({\frac {\partial F_{A}}{\partial N_{A}}}\right)_{T,V_{A}}dN_{A}+\left({\frac {\partial F_{B}}{\partial N_{B}}}\right)_{T,V_{B}}dN_{B}$

אך מכוון שחלקיקים היוצאים ממערכת A יכולים לעבור אך ורק למערכת B בהכרח מתקיים $dN_{A}=-dN_{B}$ . נציב זאת, נשווה לאפס מהדרישה לשיווי משקל ונקבל:

$\left[\left({\frac {\partial F_{A}}{\partial N_{A}}}\right)_{T,V_{A}}-\left({\frac {\partial F_{B}}{\partial N_{B}}}\right)_{T,V_{B}}\right]dN_{A}=0$

מכוון שנרצה שהשוויון יהיה נכון לכל $dN_{A}$ בהכרח נקבל:

$\left({\frac {\partial F_{A}}{\partial N_{A}}}\right)_{T,V_{A}}=\left({\frac {\partial F_{B}}{\partial N_{B}}}\right)_{T,V_{B}}$

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

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

כעת נתייחס לאטמוספירה בצורה מקורבת כאל עמוד גז הנמצא בשיווי משקל ובטמפרטורה קבועה $T_{0}$ לכל אורכו. נחלק את עמוד הגז לשכבות שונות לכל גובהו. השכבות השונות נמצאות במגע דיפוזיבי זו עם זו ובשיווי משקל.

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

$\ F(T,N,V)=-Nk_{B}T\left(\log {\left({\frac {Vn_{Q}(T)}{N}}\right)}+1\right)$

כאשר $V$ נפח השכבה, $N$ מספר החלקיקים הכולל בשכבה, $n_{Q}(T)$ הצפיפות הקוונטית אשר נתונה על ידי $n_{\rm {Q}}=\left({\frac {mk_{b}T}{2\pi \hbar ^{2}}}\right)^{3/2}$ , $m$ המסה הממוצעת של חלקיק בגז, $k_{B}$ קבוע בולצמן, ו- $\hbar$ קבוע פלאנק המצומצם.

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

$\ F=-Nk_{B}T_{0}\left(\log {\left({\frac {Vn_{Q}(T_{0})}{N}}\right)}+1\right)+Nmgh$

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

$\left({\frac {\partial F}{\partial N}}\right)_{T,V}=k_{B}Tln\left({\frac {n}{n_{Q}}}\right)$

ולכן נגזור את האנרגיה החופשית הכוללת של כל אחת מהשכבות לפי מספר החלקיקים ונשווה, נקבל את המשוואה:

$k_{B}T_{0}ln\left({\frac {n(0)}{n_{Q}}}\right)=k_{B}T_{0}ln\left({\frac {n(h)}{n_{Q}}}\right)+mgh$

נשתמש בחוקי הלוגריתם כדי לרשום:

$ln\left({\frac {n(h)}{n(0)}}\right)=-{\frac {mgh}{k_{B}T_{0}}}$

ניקח אקספוננט של שני האגפים ונקבל:

$n(h)=n_{0}e^{-{\frac {mgh}{k_{B}T_{0}}}}$

לבסוף נשתמש במשוואת המצב של הגז האידיאלי $p=nk_{B}T$ כדי לקבל את הנוסחא הברומטרית:

$p=p_{0}e^{-{\frac {mgh}{k_{B}T_{0}}}}$

ראו גם

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

לקריאה נוספת

David Goodstein, Thermal Physics, Cambridge University Press, 2015
Charles Kittel and Herbert Kroemer, Thermal Physics, University of California, W.H. Freeman and Company, 1980

הערות שוליים והפניות

^ Mario N. Berberan-Santos, Evgeny N. Bodunov and Lionello Pogliani, On the Barometric formula, 1996
^ התנאי המדויק עבור גז חלקיקים ללא אינטראקציה להחשב גז אידיאלי הוא שצפיפות הגז תהיה קטנה בהרבה מהצפיפות הקוונטית, כלומר $n<<n_{Q}$ כאשר הצפיפות הקוונטית נתונה על ידי $n_{\rm {Q}}=\left({\frac {mk_{b}T}{2\pi \hbar ^{2}}}\right)^{3/2}$ ו- $\hbar ,k_{b},T,m$ הם מסת החלקיקים, הטמפרטורה, קבוע בולצמן וקבוע פלאנק המצומצם בהתאמה
^ לקוח מתוך U.S Standard Atmosphere 1976 (קישור). הנתונים ב-U.S standard atmosphere נמדדים ביחידות גובה גאופוטנציאלי (Geopotential height), שזוהי סקלת גובה מקובלת בגאופיזיקה הלוקחת בחשבון את תלות תאוצת הכובד בקו הגובה, ובגובה עצמו. מכוון שהתיקון קטן יחסית, ואנו מעוניינים בהתנהגות כללית, בערך זה הנחנו תאוצת כובד קבועה ולכן ניתן להתייחס לנתוני הגובה כאל גובה אבסולוטי מעל פני הים.

[1] Mario N. Berberan-Santos, Evgeny N. Bodunov and Lionello Pogliani, On the Barometric formula, 1996

[2] התנאי המדויק עבור גז חלקיקים ללא אינטראקציה להחשב גז אידיאלי הוא שצפיפות הגז תהיה קטנה בהרבה מהצפיפות הקוונטית, כלומר $n<<n_{Q}$ כאשר הצפיפות הקוונטית נתונה על ידי $n_{\rm {Q}}=\left({\frac {mk_{b}T}{2\pi \hbar ^{2}}}\right)^{3/2}$ ו- $\hbar ,k_{b},T,m$ הם מסת החלקיקים, הטמפרטורה, קבוע בולצמן וקבוע פלאנק המצומצם בהתאמה

[3] לקוח מתוך U.S Standard Atmosphere 1976 (קישור). הנתונים ב-U.S standard atmosphere נמדדים ביחידות גובה גאופוטנציאלי (Geopotential height), שזוהי סקלת גובה מקובלת בגאופיזיקה הלוקחת בחשבון את תלות תאוצת הכובד בקו הגובה, ובגובה עצמו. מכוון שהתיקון קטן יחסית, ואנו מעוניינים בהתנהגות כללית, בערך זה הנחנו תאוצת כובד קבועה ולכן ניתן להתייחס לנתוני הגובה כאל גובה אבסולוטי מעל פני הים.

[1]

[2]

[3]