blob: 8dce0be28e37cfa486bcb23efffea6973dec8d34 [file] [log] [blame]
avm99963de16cf62021-01-21 13:00:04 +01001\documentclass[11pt,a4paper]{article}
2\usepackage[utf8x]{inputenc}
3\usepackage[catalan]{babel}
4\usepackage{fancyhdr}
5\usepackage{graphicx}
6\usepackage[labelfont=bf]{caption}
7\usepackage{siunitx}
8\usepackage{geometry}
9\geometry{top=25mm}
10\usepackage{amsmath}
11\usepackage{booktabs}
12\usepackage{chemformula}
13\usepackage{multicol}
14\usepackage{hyperref}
15
16\usepackage{pgfplotstable}
17\pgfplotsset{compat=1.16}
18\pgfplotstableset{
19empty cells with={--}, % replace empty cells with ’--’
20every head row/.style={before row=\toprule,after row=\midrule},
21every last row/.style={after row=\bottomrule},
22set thousands separator={\,}%,
23%every even row/.style={
24%before row={\rowcolor[gray]{0.9}}}, % Add this for stylish tables ;)
25%begin table=\begin{longtable},
26%end table=\end{longtable}
27}
28
29\setlength{\parskip}{1em}
30
31\pagestyle{fancy}
32\fancyhf{}
33\rhead{Adrià Vilanova Martínez}
34\lhead{Pràctica 9}
35\rfoot{\thepage}
36
37%%%% Title %%%%
38\title{\vspace{-2ex}Pràctica 9. Determinació de la conductivitat tèrmica d'un metall\vspace{-2ex}}
39\author{Adrià Vilanova Martínez (T1B)\vspace{-2ex} }
40\date{Tardor 2020}
41
42\begin{document}
43 \maketitle
44
45 \section{Objectiu de la pràctica}
46 L'objectiu és determinar la conductivitat tèrmica d'un metall. Això es farà escalfant dues barres de metalls diferents per un extrem. A partir de les equacions de difusió de calor i del fet que coneixem la conductivitat tèrmica d'un dels metalls, es deduirà amb les dades obtingudes per un tercer de l'experiment la conductivitat tèrmica de l'altre metall.
47
48 El tercer ha realitzat mesures de la temperatura a diferents punts de la barra en diferents moments. Després d'un temps, en un estat estacionari, ha fet mesures més precises de la temperatura a més punts de la barra. De fet, realment el que ha mesurat el tercer no és la temperatura directament, sinó el voltatge d'un termoparell, que és una propietat termomètrica a partir de la qual es pot obtenir la temperatura, a partir del calibratge del termoparell que també ha fet el tercer (la anàlisi d'aquest calibratge està tractada a l'informe de la pràctica 1).
49
50 La descripció completa i detallada de l'experiment es pot trobar al Guió de Pràctiques de Termodinàmica.
51
52 \section{Tractament de dades}
53 Les dades utilitzades són les del fitxer \texttt{Mesura 1.xslx}.
54
55 A partir del calibratge mesurat pel tercer ($\varepsilon = a T + b T^2$ on $a= (3.72 \pm 0.06) \cdot 10^{-2} \, \si{\milli\volt\per\celsius}$, $b= (2.8 \pm 0.6) \cdot 10^{-5} \, \si{\milli\volt\per\celsius\squared}$) es poden convertir els valors de voltatge en temperatura prenent l'arrel positiva de l'equació de segon ordre $- \varepsilon + a T + b T^2 = 0$. Així doncs, la temperatura es pot expressar com: \[ T = \frac{- a + \sqrt{a^2 + 4 b \varepsilon}}{2b} \]
56
57 Per tal d'obtenir les incerteses associades a la temperatura $t$ a partir de les incerteses de les constants $a$, $b$ i els voltatges $\varepsilon$, es pot fer-ho pensant $T \equiv T(a, b, \varepsilon)$, i per tant usant l'expressió de la incertesa per una funció multivariable: \[ \delta T = \left| \frac{\partial T}{\partial a} \right| \delta a + \left| \frac{\partial T}{\partial b} \right| \delta b + \left| \frac{\partial T}{\partial \varepsilon} \right| \delta \varepsilon \]
58
59 Si es calculen les derivades parcials s'arriba a les següents expressions: \[ \def\arraystretch{2.2} \left\{ \begin{array}{l}
60 \displaystyle \frac{\partial T}{\partial a} = \frac{a \Delta^{-1} - 1}{2b} \\
61 \displaystyle \frac{\partial T}{\partial b} = \frac{\varepsilon}{b \Delta} - \frac{\Delta - a}{2 b^2} \\
62 \displaystyle \frac{\partial T}{\partial \varepsilon} = \Delta^{-1}
63 \end{array} \right. \] on $\Delta := \sqrt{a^2 + 4 b \varepsilon}$ és el discriminant de l'equació de segon grau.
64
65 Això permet calcular els següents valors de la temperatura en funció del temps i la distància al forn de cada vareta, amb les seves incerteses associades:
66
67 \textbf{Barra de coure:}
68
69 \begin{center}
70 \begin{minipage}{\textwidth}
71 \begin{multicols}{2}
72 \begin{center}
73 \pgfplotstabletypeset[
74 columns/0/.style={column name=$t \, (\si{\second})$},
75 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
76 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
77 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
78 {../data/coure/10.dat}
79 \captionof{table}{Valors per la barra de coure al punt $x = \SI{10}{\centi\meter}$.}
80 \end{center}
81
82 \begin{center}
83 \pgfplotstabletypeset[
84 columns/0/.style={column name=$t \, (\si{\second})$},
85 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
86 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
87 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
88 {../data/coure/20.dat}
89 \captionof{table}{Valors per la barra de coure al punt $x = \SI{20}{\centi\meter}$.}
90 \end{center}
91 \end{multicols}
92 \end{minipage}
93 \end{center}
94
95 \begin{center}
96 \pgfplotstabletypeset[
97 columns/0/.style={column name=$t \, (\si{\second})$},
98 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
99 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
100 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
101 {../data/coure/30.dat}
102 \captionof{table}{Valors per la barra de coure al punt $x = \SI{30}{\centi\meter}$.}
103 \end{center}
104
105 \textbf{Barra de ferro:}
106
107 \begin{center}
108 \begin{minipage}{\textwidth}
109 \begin{multicols}{2}
110 \begin{center}
111 \pgfplotstabletypeset[
112 columns/0/.style={column name=$t \, (\si{\second})$},
113 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
114 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
115 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
116 {../data/ferro/10.dat}
117 \captionof{table}{Valors per la barra de ferro al punt $x = \SI{10}{\centi\meter}$.}
118 \end{center}
119
120 \begin{center}
121 \pgfplotstabletypeset[
122 columns/0/.style={column name=$t \, (\si{\second})$},
123 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
124 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
125 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
126 {../data/ferro/20.dat}
127 \captionof{table}{Valors per la barra de ferro al punt $x = \SI{20}{\centi\meter}$.}
128 \end{center}
129 \end{multicols}
130 \end{minipage}
131 \end{center}
132
133 \textbf{Règim estacionari:}
134
135 \begin{center}
136 \begin{minipage}{\textwidth}
137 \begin{multicols}{2}
138 \begin{center}
139 \pgfplotstabletypeset[
140 columns/0/.style={column name=$x \, (\si{\centi\meter})$},
141 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
142 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
143 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
144 {../data/estacionari/coure.dat}
145 \captionof{table}{Valors per la barra de \textbf{coure} un cop arribat al règim estacionari.}
146 \end{center}
147
148 \begin{center}
149 \pgfplotstabletypeset[
150 columns/0/.style={column name=$x \, (\si{\centi\meter})$},
151 columns/1/.style={column name=$\varepsilon \, (\si{\milli\volt})$, fixed, fixed zerofill, precision=2},
152 columns/2/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1},
153 columns/3/.style={column name=$\delta T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}]
154 {../data/estacionari/ferro.dat}
155 \captionof{table}{Valors per la barra de \textbf{ferro} un cop arribat al règim estacionari.}
156 \end{center}
157 \end{multicols}
158 \end{minipage}
159 \end{center}
160
161 \begin{center}
162 \centering
163 \vspace{-2em}
164 \input{../output/coure.tex}
165 \captionof{figure}{Mesures de la barra de \textbf{coure} a mesura que passa el temps i la barra es va escalfant, amb les incerteses associades. Es pot observar com s'arriba a l'estat estacionari.}
166 \end{center}
167
168 \begin{center}
169 \centering
170 \vspace{-2em}
171 \input{../output/ferro.tex}
172 \captionof{figure}{Mesures de la barra de \textbf{ferro} a mesura que passa el temps i la barra es va escalfant, amb les incerteses associades. Es pot observar com s'arriba a l'estat estacionari.}
173 \end{center}
174
175 \begin{center}
176 \centering
177 \vspace{-2em}
178 \input{../output/estacionari.tex}
179 \captionof{figure}{Logaritme de la temperatura de les dues barres en l'estat estacionari, en funció de la distància al forn, restada de la temperatura ambient $T_0 = \SI{20}{\celsius}$.}
180 \end{center}
181
182 Segons la teoria desenvolupada al Guió de Pràctiques, en el règim estacionari es té que $\log \Theta_i = C - \alpha_i x$, i que a partir dels coeficients $\alpha_i$ es pot obtenir una relació entre els coeficients de conductivitat tèrmica d'ambdues barres: \[ \frac{k_F}{k_C} = \left(\frac{\alpha_C}{\alpha_F}\right)^2 \]
183
184 Per tant, fent un ajust lineal a les dades de l'estat estacionari podem obtenir els dos coeficients $\alpha_i$ i calcular $k_F$ tenint en compte que $k_C = \SI{3.97}{\watt\per\centi\meter\per\kelvin}$.
185
186 Pels ajusts lineals mostrats a la figura 3 obtenim els següents coeficients: \[ \begin{cases}
187 \alpha_C = (222 \pm 5) \cdot 10^{-4} \, \si{\per\centi\meter} \\
188 \alpha_F = (53 \pm 2) \cdot 10^{-3} \, \si{\per\centi\meter} \\
189 \end{cases} \]
190
191 \section{Conclusió}
192
193 A partir del desenvolupament anterior s'ha deduït que el coeficient de conductivitat tèrmica del ferro és $k_F = 0.70 \pm 0.05 \, \si{\watt \per \centi\meter \per \kelvin}$. No obstant, segons (CRC) el valor de la conductivitat tèrmica del ferro és de $\SI{0.802}{\watt \per \centi\meter \per \kelvin}$. El valor de la literatura cau aproximadament a 2 marges d'incertesa del valor experimental, i per tant no es pot determinar amb certa confiança si els valors són compatibles o no.
194
195 Un dels obstacles més grans a l'hora de deduïr el valor de la conductivitat tèrmica ha sigut fitar la corba del ferro a la figura 3, degut al fet que que el comportament lineal del nostre model només es dona experimentalment en valors petits d'$x$. El fet d'escollir quins punts presenten el comportament lineal ha afegit bastanta incertesa addicional (que no s'ha tingut en compte al calcular la incertesa) a la $k_f$.
196
197 És per aquesta font d'incertesa que probablement els dos valors siguin compatibles, però en tot cas s'hauria de reptir l'experiment per poder assegurar-ho.
198
199 Un altre punt fluix és el fet que el model utilitzat per l'experiment és un model ideal i per tant no correspon totalment amb la realitat, tal com es pot veure clarament a la figura 3.
200
201 \section{Bibliografia}
202
203 (CRC): David R. Lide (ed), \textit{CRC Handbook of Chemistry and Physics, 84th Edition}. CRC Press. Boca Raton, Florida, 2003; Secció 12, Properties of Solids; Thermal and Physical Properties of Pure Metals
204
205\end{document}