Dilatación y enfriamiento de barras metálicas

Mikel Solaguren Anaut

(Lucı́a Cilla)
Grupo A
11 de octubre de 2022

Resumen
Se ha realizado un ajuste lineal del tiempo en función de ln(n0 /n) para las tres diferentes barras y se han
obtenido las pendientes mp = (−265 ± 10)s, mn = (−205 ± 5)s y mb = (−209 ± 7)s a partir de las cuales se
han determinado las constantes de enfriamiento Kp = (3,77 ± 0,14) · 10−3 s−1 , Kn = (4,88 ± 0,12) · 10−3 s−1 y
Kb = (4,78 ± 0,16) · 10−3 s−1 para la barra pulida, la pintada de negro y la pintada de blanco respectivamente.
Estos valores son compatibles con la relación esperada de Kn > Kb > Kp .

1. Introducción relacionar directamente la diferencia entre la longitud

inicial y la dependiente de la temperatura con el núme-
En el presente experimento se ha analizado el com- ro de divisiones[I] .
portamiento de unas barras huecas metálicas (de latón
concretamente) al someterlas a un cambio en su tem-
peratura haciendo circular vapor en su interior. ∆L = L(T ) − L0 = L0 αL (T − TE ) = an(T ) (3)
Para poder realizar un análisis cuantitativo de un
objeto que sufre un cambio de temperatura se usa la Donde L0 es longitud inicial de la barra, αL es su coefi-
siguiente ecuación diferencial que relaciona el cambio ciente de dilatación lineal y a es el ratio de conversión
de temperatura respecto al tiempo con la diferencia de entre divisiones del dilatómetro y la variación de lon-
temperaturas entre el objeto y el entorno de la siguiente gitud de la barra.
manera: Conociendo las expresiones (2) y (3) y usando la re-
gla de la cadena es posible obtener una expresión para
dT la velocidad con la que cambian las divisiones:
= −K(T − TE )β (1)
dt
Donde K = Ah es la constante de enfriamiento deter- dn ∂ n dT L0 α dT
= = (4)
minada por el área del cuerpo en contacto con el fluido dt ∂T dt a dt
(A) y el coeficiente de transmisión superficial (h) que Con esta última expresión y usando la Ley de Newton
depende de las propiedades de cada objeto. del enfriamiento ((1) con β = 1) es posible determinar
Como en este experimento las propiedades del flui- una expresión para el número de divisiones en función
do no cambian drásticamente con la temperatura y el del tiempo:
vapor de agua es rápidamente evacuado del tubo por
estar este inclinado, se puede suponer una convección n(t) = n0 e−Kt (5)
forzada y se puede particularizar la expresión (1) para
β = 1 obteniendo ası́ la Ley de Newton del enfriamien- Esta última relación es posible reescribirla de la si-
to. Al solucionar la ecuación diferencial se obtiene la guiente manera:
siguiente expresión para la temperatura en función del  
tiempo: 1 n
t = − ln (6)
K n0
T (t) = (T0 − TE )e−Kt + TE (2)
De esta manera se obtiene una relación lineal entre el
Donde T0 es la temperatura del objeto en el instante tiempo y el logaritmo del cociente de n entre n0 . Esta
inicial de tiempo. expresión es de gran utilidad ya que permite obtener
Por otro lado, es necesario conocer una expresión K a partir de un ajuste lineal que se puede realizar con
cuantitativa para la dilatación lineal de los objetos al los datos medidos en el laboratorio.
ser sometidos a un aumento en su temperatura ya que, Usando la expresión (5), al ser una exponencial de-
para este experimento en particular, es mas sencillo creciente, es razonable preguntarse si es posible encon-
medir la dilatación lineal de las barras en lugar de su trar una estimación para K haciendo uso de las propie-
temperatura. dades de este tipo de funciones. Para esto, se analiza
Para esto se usa la siguiente expresión que relaciona primero la recta tangente en el punto (0, n0 ) para la
la longitud de un objeto en función de su temperatu- función que representa la expresión (5) que tendrá co-
ra. Además, teniendo en cuenta que en el experimento mo pendiente la derivada temporal de n(t) evaluada en
en particular se hace uso de un dilatómetro, se puede t = 0:

1
 cerrados la presión dentro del generador subirı́a rápi-
f (t) := mt + n0 ; m = −n0 K (7) damente llegando a ser peligroso.
Por otro lado, el montaje cuenta con una platafor-
A continuación se analiza la recta que pasa por el ma donde se colocan las barras. Esta plataforma ha
punto (0, n0 ) y su punto más próximo (τ, n0 − 1): de estar inclinada para conseguir que el fluido salga de
1 la barra ya que se esta tratando el sistema como un
h(t) := m̂t + n0 ; m̂ = − (8) sistema de convección forzada y para eso es necesario
τ
que tanto el vapor como el agua una vez se condense
Analizando donde estas dos rectas cortan el eje X, se evacuen rápidamente.
1
se ve que la primera lo hace para tf = K y la segunda Hablando ahora de las barras de latón, estas son tres
para th = n0 τ . Suponiendo ahora que los valores son barras huecas con una apertura por la que se introduce
suficientemente próximos como para aproximar la recta el vapor y otra apertura por la que este se evacua a un
tangente por la recta generada por los dos puntos, se recipiente mediante un tubo (es importante que la bo-
puede realizar la siguiente estimación de K: ca del tubo no esté sumergida en el agua del recipiente
para mantener la presión atmosférica dentro del tubo).
1
tf ≈ th ⇒ K ≈ (9) Lo que diferencia las tres barras entre si es que una es
n0 τ
de latón pulido, otra esta pintada de negro y otra está
Por otro lado, como al realizar el experimento se hace pintada de blanco. Las barras miden 640mm de largo
[II]
uso de un generador de vapor, es conveniente compren- . Cabe destacar que también se disponı́a de barras
der al menos de una manera general las leyes fı́sicas que de acero pero se decidió usar el latón porque se ha con-
sirven para describir el funcionamiento de este apara- siderado que los efectos del cambio de temperatura van
to. Para esto, es necesario conocer a qué temperatura a ser mas notables en este material. Al principio del ex-
se evapora el agua dada una presión. Esta relación se perimento se ha estimado que nlat ac
max ≈ 96 y nmax ≈ 48
obtiene a partir de la Ecuación de Clausius Clapeyron lo que indica que el latón se va a dilatar más. Además,
tomando el vapor de agua como un gas ideal y tenien- al tener el acero mayor densidad y capacidad calorı́fica
do en cuenta que el volumen molar del vapor de agua a volumen constate, se espera que este material tome
es mucho mayor que aquel del agua. Bajo estas condi- mas tiempo para alcanzar el equilibrio de temperatura.
ciones, resolviendo la ecuación diferencial de Clausius Por último, se dispone de un dilatómetro que se co-
Clapeyron se obtiene la siguiente relación: loca en el extremo suelto de la barra tras haber fijado
    el otro a la plataforma. Este dispositivo permite medir
P ∆hf 1 1 la variación de la longitud de la barra en función de
ln = − (10)
P0 R T0 T 100 divisiones en un cı́rculo.
Donde (T0 , P0 ) es un punto cualquiera de la curva de Previo a la realización del experimento se discute el
vapor del agua, ∆hf es la entalpı́a de vaporización del efecto que va a tener el color de la pintura o la falta
agua y R es la constante de gases. de esta en el coeficiente de enfriamiento. Por la Ley de
Por último, para finalizar este apartado, cabe men- Kirchoff de la radiación térmica, se sabe que en equi-
cionar la Ley de Kirchoff de la radiación térmica [3] que librio termodinámico si un objeto absorbe bien la ra-
cualitativamente dice que la capacidad de absorber y diación térmica también la va a emitir bien. Por esto
emitir radiación térmica de un material arbitrario en se puede esperar que la barra pulida tenga menor emi-
equilibrio termodinámico es la misma. Esto es impor- sividad ya que esta es la que se considera que mejor
tante para realizar un análisis cualitativo de como se refleja la radiación. Por otro lado, entre las dos barras
van a enfriar las diferentes barras (la pulida, la pintada restantes, se espera que aquella pintada de negro tenga
de negro y la pintada de blanco). mayor emisividad ya que esta va a estar más cerca de
comportarse como un cuerpo negro, y por tanto de ab-
sorber y emitir la mayor cantidad posible de radiación
2. Método experimental que aquella pintada de blanco. Teniendo en cuenta es-
to, se espera que la barra negra sea la que mas rápido
En este apartado del trabajo se procede a describir se enfrı́a seguida por la blanca y luego por la pulida.
el dispositivo experimentar y se comentan las estima- Con esto en mente, como la expresión (1) para β = 1 es
ciones que se han realizado previas a las medidas. una relación lineal entre la diferencia de temperaturas
El sistema que calienta las barras en el dispositivo y la velocidad a la que se enfrı́a el material, se puede
experimental es un generador de vapor. El vapor ge- deducir que Kn > Kb > Kp (Las constantes de enfria-
nerado por este puede salir por dos tubos diferentes, miento para la barra pintada de negro, la pintada de
ambos con una llave de manera que el flujo de vapor se blanco y la pulida respectivamente).
puede cortar. El primero de los tubos conduce a una en- Para realizar el experimento, primero se ha de es-
trada que lleva el vapor al interior de la barra metálica. timar la temperatura a la que sale el vapor del gene-
El segundo tubo, conectado a un recipiente de cristal, rador para lo que se utiliza la temperatura de vapor
se utiliza como salida cuando no se quiere llenar la ba- del agua a la presión de la habitación [III] . Al encon-
rra de vapor. Este segundo tubo es imprescindible si trarse Santander a h = 7m sobre el nivel del mar, se
se quiere mantener la presión atmosférica dentro del puede estimar esta altura como igual a la del nivel del
generador de vapor ya que si ambos tubos estuvieran mar y por tanto la presión sera la presión atmosféri-

2
ca Patm = 1,013 · 105 Pa. Utilizando ahora la expresión una de las barras distintas. Estas tablas son como las
(10) con (T0 = 373,15K, P0 = 1,013Pa)[1] se ve que el mencionadas en el apartado anterior. A continuación se
miembro de dentro del logaritmo se puede aproximar muestran las tablas junto con la estimación de K para
a uno, por tanto el logaritmo es cero lo que fuerza que cada una de las varillas usando el método que conduce
T ≈ T0 = 373,15K. a la expresión (9).
Una vez conocida la temperatura del vapor que se Dado que las tres tablas comparten el mismo for-
circula por dentro del tubo, como el número de divi- mato, se ha decidido mencionar qué es cada columna
siones del dilatómetro puede expresarse exclusivamente ahora en lugar de mencionarlo en cada una de ellas.
en función de la temperatura por la expresión (3), se En estas tablas se muestra el número de la división en
pueden estimar las divisiones a las que va a llegar el la que se hace cada una de las medidas (n), el tiempo
dilatómetro dada la temperatura ambiente del labo- estimado en el que se encuentra en cada división con
ratorio que se ha tomado como TE ≈ 293,15K. Esta su error (test y ∆test ) y el tiempo experimental en el
estimación ha sido de n(393,15) ≈ 96 divisiones (El que se ha medido cada división (texp ).
valor experimental en el caso de la varilla pulida ha
resultado ser de nmax = 88 divisiones).
Una vez realizadas estas estimaciones se prepara una Cuadro 1: En esta tabla se muestran las medidas
tabla en la que se van a apuntar las divisiones y sus realizadas para la barra pulida. Con los datos aquı́
tiempos correspondientes. En este punto se plantea la mostrados y usando la expresión (9), se ha obtenido
duda de cual de estas dos variables fijar como la va- un valor estimado de Kest = (4,51 ± 0,16) · 10−3 s−1
[IV]
riable independiente. Si se tomara el tiempo como va- . Con este valor estimado se ha determinado la
riable independiente, se tendrı́a que mirar en que divi- segunda columna de la tabla y haciendo uso de la
sión se encuentra la aguja cuando el cronómetro marca expresión (12) se ha determinado su error.
un tiempo especı́fico. Esta medida darı́a un valor de
n/div test /s ∆test (texp ± 0,1)/s
n(ti ) = n±1 lo que arrastrarı́a una gran incertidumbre
81 0 0 0.0
a cualquier cálculo que se quisiera hacer a partir de esos
80 2.76 0.10 2.7
valores. Sin embargo, al tomar las divisiones como va- 79 5.5 0.2 5.2
riable independiente, cuando el dilatómetro marca una 78 8.4 0.3 8.0
división ni concreta, se mide un tiempo t(ni ) = t ± δ 76 14.1 0.5 14.5
74 20.1 0.7 21.0
con δ ≈ 0,1 tomando como error el tiempo de reac-
70 32.4 1.1 34.8
ción humano en lugar de la precisión del cronómetro 67 42.1 1.5 45.3
ya que esto último no serı́a un valor realista. Es por 64 52.3 1.9 55.4
esto que se decide tomar las divisiones del dilatómetro 61 63 2 67.8
56 82 3 90.2
como variable independiente en el experimento.
51 103 4 115.3
Siguiendo con la tabla, todavı́a falta por determinar 46 126 4 142.7
cuales son las divisiones para las cuales se va a medir el 40 157 6 181.2
tiempo. Se ha decidido empezar en la división 81 para 30 220 8 263.1
29 228 8 273.7
tener ası́ tiempo para redirigir el vapor al recipiente de
cristal y cerrar la llave del tubo que comunica el gene-
rador de vapor con la barra. A partir de la división 81,
se hacen medidas de uno en uno hasta la división 78, Cuadro 2: En esta tabla se muestran las medidas rea-
luego se va ampliando la diferencia entre divisiones ya lizadas para la barra pintada de negro. Con los datos
que como se ha comentado en el desarrollo de la expre- aquı́ mostrados y usando la expresión (9), se ha obte-
sión (9) los valores más importantes son los primeros. nido un valor estimado de Kest = (5,6 ± 0,3) · 10−3 s−1
Por otro lado, también se han tomado las medidas para [IV]
. Con este valor estimado se ha determinado la
la división número 30 y la 29 ya que se espera que el segunda columna de la tabla y haciendo uso de la
valor de n(tf ) ∈ (29, 30) siendo tf el punto en el que expresión (12) se ha determinado su error.
la recta tangente al punto inicial corta el eje x ya que
k
n(tf ) = e− k = n0 e−1 ≈ 29,8. n/div test /s ∆test (texp ± 0,1)/s
Por último, se ha estimado la constante K con los dos 81 0 0 0.0
primeros valores y el intervalo de tiempo entre estos y 80 2.23 0.12 2.2
con este valor se han obtenido unos datos estimados del 79 4.5 0.2 4.3
tiempo en el que deberı́a estar sucediendo cada división 78 6.8 0.4 6.9
76 11.5 0.6 11.8
del dilatómetro. Este proceso se ha repetido para las 74 16.3 0.9 17.2
tres barras diferentes y en los siguientes apartados se 70 26.2 1.4 27.8
presentan y discuten los resultados. 67 34.1 1.8 36.9
64 42 2 46.4
61 51 3 54.7
56 66 4 71.0
3. Resultados 51 83 4 89.6
46 102 5 111.5
En este apartado se presentan los resultados obte- 40 127 7 141.3
nidos en la realización del experimento. Estos datos 30 179 10 203.4
29 185 10 211.9
se presentan en tres tablas diferentes, una para cada

3
 Conociendo estas pendientes y con la expresión (6) se
Cuadro 3: En esta tabla se muestran las medidas rea-
obtienen los valores Kp = (3,77 ± 0,14) · 10−3 s−1 , Kn =
lizadas para la barra pintada de blanco. Con los datos
(4,88 ± 0,12) · 10−3 s−1 y Kb = (4,78 ± 0,16) · 10−3 s−1 .
aquı́ mostrados y usando la expresión (9), se ha obte-
Estos valores presentan unos errores relativos con los
nido un valor estimado de Kest = (5,0 ± 0,2) · 10−3 s−1
[IV] valores estimados de δp ≈ 16 %, δn ≈ 13 %, δb ≈ 4 %
. Con este valor estimado se ha determinado la
y cumplen que Kn > Kb > Kp como se ha comentado
segunda columna de la tabla y haciendo uso de la
previamente.
expresión (12) se ha determinado su error.

n/div test /s ∆test /s (texp ± 0,1)/s

81 0 0 0.0
5. Discusión
80 2.51 0.10 2.5
79 5.0 0.2 4.6 Como se ha visto en el apartado anterior, los valo-
78 7.6 0.3 6.9 res obtenidos para la contante de enfriamiento cumplen
76 12.9 0.5 11.9 que Kn > Kb > Kp . Más concretamente, se ve que hay
74 18.2 0.7 16.5
una gran diferencia entre la barra pulida y las dos pin-
70 29.4 1.2 27.9
67 38.3 1.5 35.8 tadas. Por otro lado, entre la barra pintada de negro
64 47.5 1.9 44.6 y la pintada de blanco, a pesar de haber una ligera
61 57 2 54.8 diferencia en sus constantes de enfriamiento, esta di-
56 74 3 71.6 ferencia es muy pequeña comparada con aquella que
51 93 4 91.5
46 114 5 113.1
hay entre la pulida y las pintadas. De esto se puede
40 142 6 143.3 concluir que independientemente del color, una capa
30 200 8 207.0 de pintura mate hace aumentar considerablemente la
29 207 8 216.4 capacidad de un material de emitir radiación térmica.
Analizando ahora las dos barras pintadas, se puede
concluir que la pintada de negro tiene mayor emisivi-
4. Análisis dad y por tanto, por la Ley de Kirchoff de la radiación
térmica, mas absorbencia de radiación; lo que concuer-
En este apartado se muestra la gráfica obtenida tras da con la hipótesis de que esta barra está mas cerca de
representar el tiempo que se ha medido en el laborato- poder considerarse un cuerpo negro que la blanca.
rio en función del logaritmo del cociente n entre n0 . La En lo que a los valores estimados se trata, se puede
representación de los datos se espera que siga la expre- ver que a pesar de que estos concuerden para los va-
sión lineal (6) de manera que la pendiente del ajuste lores de mayor número de división, comienza a haber
1
lineal sea m = − K . una mayor diferencia entre los datos estimados y expe-
rimentales según aumenta el tiempo. Mirando al valor
de la función en el punto de corte de la recta tangente
del primer punto [IV] , se ve que este no se encuentra en
el intervalo esperado en ninguno de los tres casos sien-
do esto mas evidencia de que los datos experimentales
y estimados difieren para t grandes
Por último, cabe mencionar que los errores relati-
vos entre los valores experimentales y estimados de la
constante de enfriamiento son grandes en parte debido
a que la aproximación que se usa para los valores esti-
mados es solo válida para los valores que se toman al
inicio del experimento. Valores que resulta que son los
más difı́ciles de medir por la naturaleza exponencial de
los resultados.

6. Conclusiones
De este experimento se concluye que una capa mate
de pintura aplicada sobre un objeto ayuda a aumentar
Figura 1: En esta figura se muestra la dependencia su emisividad. Por otro lado, se ve que, aunque afec-
del tiempo frente a ln(n/n0 ). Las pendientes de los tando en menor medida, el color de la pintura también
diferentes ajustes lineales son mp = (−265 ± 10)s, afecta a la emisividad del material siendo el negro el
mn = (−205 ± 5)s, mb = (−209 ± 7)s siendo estas color que mejor emite y el blanco el que peor.
las pendientes de los ajustes para la barra pulida, la Resultados como este pueden ser utilizados de mane-
barra pintada de negro y la barra pintada de blanco ra práctica cuando se quiere que un cierto objeto tarde
respectivamente. Errores en las figuras 2, 3 y 4 más o menos en alcanzar el equilibrio de temperatura
[2].

4
 En lo que a la precisión del experimento respecta,
se considera que las medidas de este han tenido una
gran precisión en parte por la elección de la variable
independiente.

Apéndice
Cálculo de errores:
Error de Kest :
 
 dKest 
∆Kest =   ∆τ = 1 ∆τ (11)
dτ  n0 τ 2

Error de K: Figura 4: En esta figura se muestra el proceso del cálcu-

 
 dK  lo del error para la pendiente en el ajuste lineal de la
∆K =   ∆m = ∆m (12) barra pintada de blanco. Haciendo las respectivas res-
dm  m2
tas entre las pendientes m1 y m2 con la pendiente del
Errores de las pendientes del ajuste lineal: ajuste m se obtiene que δ1 = 7 y δ2 = 6 y por tanto
δ=7

Notas:

I Al existir una relación directa entre el número de

divisiones del dilatómetro y la variación de la lon-
gitud de la barra (∆L = an(T )) es mas sencillo
trabajar directamente con las divisiones del apa-
rato en lugar de con longitudes de la barra.
II A pesar de que las barras midan 640mm al mon-
tarlas sobre la plataforma y fijar uno de los lados
de la barra la parte que queda libre a la derecha del
tornillo con el que se fijan es de 600mm por tanto
Figura 2: En esta figura se muestra el proceso del cálcu- el efecto de dilatación que se aprecia es sobre esos
lo del error para la pendiente en el ajuste lineal de la 600mm de barra no sobre los 640mm totales.
barra pulida. Haciendo las respectivas restas entre las III El generador de vapor consta de un deposito de
pendientes m1 y m2 con la pendiente del ajuste m se agua que se hace hervir hasta evaporarla. Al tra-
obtiene que δ1 = 10 y δ2 = 9 y por tanto δ = 10 tarse de un sistema agua-vapor, la temperatura
del segundo no puede superar la temperatura de
vaporización del agua y por tanto vasta estimar
esta temperatura para la presión correspondiente.

IV Con los valores de K = 1/tf estimados se espera

que los valores tpf = 221,7s, tnf = 178,6s y tbf =
200,0s estén entre los valores medidos de t30 y t29 .

Referencias
[1] Engineering ToolBox, (2003). Liquids and
Gases - Boiling Points. [online] Available at:
[Link]
points-fluids-gases-d [Link]
Figura 3: En esta figura se muestra el proceso del cálcu-
lo del error para la pendiente en el ajuste lineal de la [2] Baruch Givoni. Passive Low Energy Cooling of
barra pintada de negro. Haciendo las respectivas res- Buildings. Pág 24-31.
tas entre las pendientes m1 y m2 con la pendiente del [3] Institution of Engineering and Technology. Ther-
ajuste m se obtiene que δ1 = 5 y δ2 = 5 y por tanto mal Microwave Radiation: Applications for Remo-
δ=5 te Sensing. Pág 5-7.