Prawo stygnięcia

Prawo stygnięcia (prawo stygnięcia Newtona) – w fizyce prawo określające z jaką szybkością ciała przekazują sobie energię cieplną w wyniku przewodnictwa ciepła. Prawo zostało sformułowane przez Izaaka Newtona.

Prawo nie obowiązuje jeżeli przekazywanie energii cieplnej odbywa się przez promieniowanie cieplne, konwekcję lub przewodzeniu towarzyszy zmiana stanu skupienia (np. parowanie).

Sformułowanie prawa

Prawo stygnięcia (prawo stygnięcia Newtona) mówi, że:

„Szybkość z jaką układ stygnie jest proporcjonalna do różnicy temperatur pomiędzy układem a otoczeniem”.

Matematycznie można to wyrazić jako:

d T d t = k ( T T R ) = k Δ T , {\displaystyle {\frac {\mathrm {d} T}{\mathrm {d} t}}=-k\left(T-T_{\mathrm {R} }\right)=-k\Delta T,}

gdzie:

T {\displaystyle T} – temperatura ciała,
T R {\displaystyle T_{\mathrm {R} }} – temperatura otoczenia,
Δ T {\displaystyle \Delta T} – różnica temperatur układu i otoczenia,
t {\displaystyle t} – czas,
k {\displaystyle k} – stała dla danego układu (zależna m.in. od fizycznej wielkości układu, jego pojemności cieplnej i jego wewnętrznej struktury, przenikalności cieplnej ścianek układu, rodzaju otoczenia).

Stygnięcie przy stałej temperaturze otoczenia

Z powyższego, przy założeniu stałości temperatury otoczenia, otrzymujemy eksponencjalną zależność temperatury stygnącego układu od czasu stygnięcia:

T ( t ) T R = Δ T ( t ) = Δ T ( 0 )   e k t , {\displaystyle T(t)-T_{\mathrm {R} }=\Delta T(t)=\Delta T(0)\ \mathrm {e} ^{-kt},}

gdzie:

Δ T ( 0 ) {\displaystyle \Delta T(0)} – początkowa różnica temperatur.

Wyprowadzenie

Krzywa ostygania ilustrująca spadek temperatury stygnącego ciała w czasie

Prawo ostygania zapisane w postaci

d T d t = k ( T T R ) , {\displaystyle {\frac {\mathrm {d} T}{\mathrm {d} t}}=-k\left(T-T_{\mathrm {R} }\right),}

gdzie:

T R {\displaystyle T_{\mathrm {R} }} – temperatura otoczenia,
T {\displaystyle T} – aktualna temperatura układu,

jest równaniem różniczkowym, w którym można rozdzielić zmienne

d T T T R = k d t , {\displaystyle {\frac {\mathrm {d} T}{T-T_{\mathrm {R} }}}=-k\mathrm {d} t,}
T 0 T d T T T R = 0 t k d t , {\displaystyle {}\,\int \limits _{T_{0}}^{T}{\frac {\mathrm {d} T}{T-T_{\mathrm {R} }}}=-\int \limits _{0}^{t}{k\mathrm {d} }t,}

gdzie:

T 0 {\displaystyle T_{0}} – temperatura początkowa układu.

Po wycałkowaniu

ln ( T T R ) | T 0 T = k t , {\displaystyle \left.\ln \left(T-T_{\mathrm {R} }\right)\right|_{T_{0}}^{T}=-kt,}
ln T T R T 0 T R = k t , {\displaystyle \ln {\frac {T-T_{\mathrm {R} }}{T_{0}-T_{\mathrm {R} }}}=-kt,}
T T R = ( T 0 T R ) exp ( k t ) {\displaystyle T-T_{\mathrm {R} }=\left(T_{0}-T_{\mathrm {R} }\right)\exp \left(-kt\right)}

i ostatecznie:

T = T R + ( T 0 T R ) exp ( k t ) . {\displaystyle T=T_{\mathrm {R} }+\left(T_{0}-T_{\mathrm {R} }\right)\exp \left(-kt\right).}

Temperaturę stygnącego ciała w funkcji czasu ilustruje krzywa ostygania.

Realność warunku stałości temperatury otoczenia

Sprzątanie Wikipedii
Tę sekcję należy dopracować:
Należy poprawić punkt 2 dotyczący przemiany fazowej..
Po wyeliminowaniu niedoskonałości należy usunąć szablon {{Dopracować}} z tej sekcji.

Gdy ciało stygnie, wówczas temperatura otoczenia może się podnosić. Warunek stałości temperatury otoczenia może być jednak utrzymany, gdy otoczenie ma dużą pojemność cieplną w porównaniu ze stygnącym przedmiotem (np. szklanka z herbatą na dworcu);

W praktyce stałość temperatury otoczenia można uzyskać przez użycie takich warunków eksperymentalnych jak:

Stygnięcie przy zmiennej temperaturze otoczenia

Założenia

Jeżeli stygnący układ i bezpośrednie otoczenie układu są odizolowane od otoczenia termodynamicznego, prawo stygnięcia Newtona pozostaje słuszne, pomimo tego że temperatura otoczenia układu nie jest stała.

Najprościej można sobie wyobrazić 2 układy odizolowane termicznie od otoczenia, a w kontakcie ze sobą poprzez przegrodę, przy czym wnętrza obu układów mają jednorodny rozkład temperatury (uzyskuje się np. poprzez mieszanie lub gdy szybkość przepływu ciepła przez przegrodę jest dużo mniejsza niż przepływ wewnątrz obu układów). Konieczne jest też założenie o (przynajmniej w przybliżeniu) stałości pojemności cieplnych obu układów (stałości ciepeł właściwych).

Przepływ ciepła przez przegrodę zależy od różnicy temperatur obu układów:

d Q 1 d t = k q ( T 1 T 2 ) . {\displaystyle {\frac {\mathrm {d} Q_{1}}{\mathrm {d} t}}=-k_{q}(T_{1}-T_{2}).}

Oba układy są izolowane od otoczenia, a więc:

d Q 2 d t = d Q 1 d t . {\displaystyle {\frac {\mathrm {d} Q_{2}}{\mathrm {d} t}}=-{\frac {\mathrm {d} Q_{1}}{\mathrm {d} t}}.}

Rozwiązanie

Różnice pojemności cieplnej obu układów (inna masa m , {\displaystyle m,} i inne ciepło właściwe C {\displaystyle C} ), powodują, że ta sama ilość ciepła (energii) zmienia temperaturę w różny sposób:

d Q 2 = d Q 1 {\displaystyle \mathrm {d} Q_{2}=-\mathrm {d} Q_{1}}   i   Δ Q 2 = Δ Q 1 , {\displaystyle \Delta Q_{2}=-\Delta Q_{1},}
d Q 1 = m 1 C 1 d T 1 {\displaystyle \mathrm {d} Q_{1}=m_{1}C_{1}\mathrm {d} T_{1}}   i   Δ Q 1 = m 1 C 1 Δ T 1 , {\displaystyle \Delta Q_{1}=m_{1}C_{1}\Delta T_{1},}
d Q 2 = m 2 C 2 d T 2 {\displaystyle \mathrm {d} Q_{2}=m_{2}C_{2}\mathrm {d} T_{2}}   i   Δ Q 2 = m 2 C 2 Δ T 2 , {\displaystyle \Delta Q_{2}=m_{2}C_{2}\Delta T_{2},}

a także:

d T 1 d T 2 = m 2 C 2 m 1 C 1 = c o n s t , {\displaystyle {\frac {\mathrm {d} T_{1}}{\mathrm {d} T_{2}}}=-{\frac {m_{2}C_{2}}{m_{1}C_{1}}}=\mathrm {const} ,}
( T 1 T e q ) = ( T 2 T e q ) m 2 C 2 m 1 C 1 , {\displaystyle (T_{1}-T_{\mathrm {eq} })=-(T_{2}-T_{\mathrm {eq} }){\frac {m_{2}C_{2}}{m_{1}C_{1}}},}

gdzie temperatura końcowa T e q {\displaystyle T_{\mathrm {eq} }} jest funkcją temperatur początkowych T 1 , 0 {\displaystyle T_{1,0}} i T 2 , 0 {\displaystyle T_{2,0}} oraz pojemności cieplnych układów:

T e q = T 1 , 0 + T 2 , 0 m 2 C 2 m 1 C 1 1 + m 2 C 2 m 1 C 1 . {\displaystyle T_{\mathrm {eq} }={\frac {T_{1,0}+T_{2,0}{\frac {m_{2}C_{2}}{m_{1}C_{1}}}}{1+{\frac {m_{2}C_{2}}{m_{1}C_{1}}}}}.}

Stąd:

d T 1 d t = k q m 1 C 1 ( T 1 T 2 ) = k q m 1 C 1 Δ T , {\displaystyle {\frac {\mathrm {d} T_{1}}{\mathrm {d} t}}=-{\frac {k_{q}}{m_{1}C_{1}}}(T_{1}-T_{2})=-{\frac {k_{q}}{m_{1}C_{1}}}\Delta T,}
d T 2 d t = k q m 2 C 2 ( T 2 T 1 ) = k m 2 C 2 Δ T , {\displaystyle {\frac {dT_{2}}{dt}}=-{\frac {k_{q}}{m_{2}C_{2}}}(T_{2}-T_{1})={\frac {k}{m_{2}C_{2}}}\Delta T,}

gdzie Δ T {\displaystyle \Delta T} jest różnicą temperatur układów „1” i „2”:

Δ T = T 1 T 2 . {\displaystyle \Delta T=T_{1}-T_{2}.}

Skąd wynika:

d Δ T d t = k T , 12 Δ T , {\displaystyle {\frac {\mathrm {d} \Delta T}{\mathrm {d} t}}=-k_{T,12}\Delta T,}

gdzie:

k T , 12 = k q ( 1 m 1 C 1 + 1 m 2 C 2 ) . {\displaystyle k_{T,12}=k_{q}\left({\frac {1}{m_{1}C_{1}}}+{\frac {1}{m_{2}C_{2}}}\right).}

I ostatecznie:

T 1 ( t ) T 2 ( t ) = Δ T ( t ) = Δ T 0 exp ( k T , 12 t ) , {\displaystyle T_{1}(t)-T_{2}(t)=\Delta T(t)=\Delta T_{0}\exp(-k_{T,12}t),}

gdzie Δ T 0 {\displaystyle \Delta T_{0}} jest początkową różnicą temperatur:

Δ T 0 = Δ T ( 0 ) = T 1 , 0 T 2 , 0 {\displaystyle \Delta T_{0}=\Delta T(0)=T_{1,0}-T_{2,0}}

oraz

T 1 ( t ) T e q = m 2 C 2 m 1 C 1 1 + m 2 C 2 m 1 C 1 Δ T o exp ( k T , 12 t ) {\displaystyle T_{1}(t)-T_{\mathrm {eq} }={\frac {\frac {m_{2}C_{2}}{m_{1}C_{1}}}{1+{\frac {m_{2}C_{2}}{m_{1}C_{1}}}}}\Delta T_{o}\exp(-k_{T,12}t)}

lub

T 1 ( t ) = T e q + ( T 1 , 0 T 2 , 0 ) m 2 C 2 m 1 C 1 + m 2 C 2 exp ( k T , 12 t ) . {\displaystyle T_{1}(t)=T_{\mathrm {eq} }+(T_{1,0}-T_{2,0}){\frac {m_{2}C_{2}}{m_{1}C_{1}+m_{2}C_{2}}}\exp(-k_{T,12}t).}

Wynik końcowy zgodny jest więc (co do charakteru przebiegu eksponencjalnego) z prawem stygnięcia Newtona dla stygnącego układu w kontakcie z otoczeniem o stałej temperaturze. To tłumaczy również sukces tego prostego prawa nawet gdy jego podstawowe założenia nie są spełnione.

Przypadek graniczny – stała temperatura otoczenia

Gdy pojemność układu „2” traktowanego tutaj jako „bezpośrednie otoczenie” jest dużo większa niż pojemność cieplna układu stygnącego:

m 2 C 2 m 1 C 1 , {\displaystyle m_{2}C_{2}\gg m_{1}C_{1},}

wówczas temperatura układu „2” (bezpośredniego otoczenia stygnącego układu) pozostaje stała:

T e q T 2 T 2 , 0 {\displaystyle T_{\mathrm {eq} }\approx T_{2}\approx T_{2,0}}

oraz

Δ T = T 1 ( t ) T 2 , 0 = Δ T 0 exp ( k T , 1 t ) , {\displaystyle \Delta T=T_{1}(t)-T_{2,0}=\Delta T_{0}\exp(-k_{T,1}t),}

gdzie współczynnik k T , 1 {\displaystyle k_{T,1}} w równaniu jest tożsamy z wartością k {\displaystyle k} w oryginalnym równaniu eksponencjalnym:

k T , 1 = k q ( 1 m 1 C 1 ) = k . {\displaystyle k_{T,1}=k_{q}\left({\frac {1}{m_{1}C_{1}}}\right)=k.}
Encyklopedia internetowa (prawo fizyki):