PRÁCTICA DE LABORATORIO, VOL. 1, NO.

3, Abril 2023

Variación de la resistencia de un conductor y de un semiconductor con la temperatura

Álvaro M. Rodríguez de Tembleque1 y Marcos P. Ipiña1

1 Universidad Complutense de Madrid, Laboratorio de termodinámica (LO1), Madrid, España

Realización de la práctica: 25/09/2023

Entrega del informa: 01/10/2023
Evaluación del informe: -

Resumen— En este informe de laboratorio realizaremos un estudio de las isotermas de un gas real (SF6 ), haciendo uso de dos ecuaciones
que permiten entender su comportamiento

Palabras clave—Ajuste no lineal, isoterma, gas real, Van Der Waals, Redlich-Kwong

I. I NTRODUCCIÓN
l estudio de los gases es una rama de la termodinámica
E de bastante recorrido histórico, pero cuya forma teóri-
ca (aquella en la que se analizan gases ideales) difiere de la
realidad en la gran mayoría de casos, es decir, solo se com-
portan de esa forma con unas ciertas condiciones, como que
se encuentren a altas temperaturas, presiones bajas o densi-
dades pequeñas. Es por ello que también hay una rama que
estudia el comportamiento de los gases reales, para los cua-
les se ajustan los modelos matemáticos existentes y así poder
definir el comportamiento de éstos mediante diversos pará-
metros obtenibles mediante experimentación.
El físico Johannes Diderik van der Waals en 1873 modificó
la ecuación de estado de los gases ideales para gases reales Figura 1: Montaje experimental
teniendo en cuenta el tamaño de las partículas y las fuerzas
de interacción entre las mismas que forman el conjunto de 1. Ampolla de alta presión con SF6 en su interior inmersa
un gas. Esa modificación recibe el nombre de ecuación de en un baño térmico
van der Waals y a través de dicha ecuación estudiaremos los
diferentes cambios de fase del SF6 para diferentes variables 2. Manómetro que indica la presión de dentro de la ampo-
de temperatura, volumen y presión; además estudiaremos la lla
interacción entre las fuerzas intermoleculares de dicho gas.
Esta ecuación se expresa como: 3. Termómetro que indica la temperatura del baño térmico

4. Termostato para regular la temperatura del baño térmico

n2 a
 
p + 2 (V − nb) = nRT (1)
V 5. Rueda para regular la presión del interior de la ampolla
contenedora del gas
Una segunda ecuación
b. Método

n2 a
 En este apartado definiremos el método mediante el cual
p+ √ (V − nb) = nRT (2) tomamos los datos necesarios para obtener los resultados es-
TV (V + nb)
perados en esta práctica que analizaremos y comentaremos
en la siguiente sección.
II. M ATERIALES Y MÉTODOS Para realizar las variaciones de temperatura calentamos el
a. Materiales y montaje agua con una resistencia eléctrica y la hacemos circular por
el sistema con una bomba que, a través de los tubos, va cam-
Los materiales que vamos a utilizar en esta práctica son: biando el agua para calentar y mantener la temperatura en los
 TERMÓMETROS

alrededores de la ampolla. Una vez establecida la temperatu- b. Para cada isoterma realice un ajuste no-lineal a las
ra y asegurarnos que es la que deseamos giramos la rueda de ecuaciones de van der Waals y de Redlich-Kwong.
mano para variar el volumen de la ampolla que contiene el Discuta la bondad de los ajustes en función de los
gas subiendo o bajando una barra de mercurio al mover la valores χ 2 de Pearson obtenidos. De los ajustes rea-
rueda; al modificar el volumen en el que se encuentra el gas lizados determine los valores de los parámetros a, b
la presión aumentará o disminuirá en consecuencia (a mayor
y n para cada ecuación. Interprete los resultados
volumen menor presión y a mayor temperatura más brusco
obtenidos.
es el cambio de presión).
Con efecto de mejorar la lectura de este informe se mos-
III. P RESENTACIÓN Y ANÁLISIS DE RESULTA - trará un ajuste a cada una de las ecuaciones descritas en la
DOS introducción

a. Dibuje en un único diagrama {p,V} todas las isoter- Ajuste de Van Der Waals
mas medidas, incluyendo la zona de coexistencia.
Dibuje aproximadamente la campana de coexisten-
cia de fases líquido-vapor y el punto crítico. A partir
de la gráfica, estime los valores experimentales de
las coordenadas críticas.

Figura 4: Ajuste no lineal de la isoterma de T = 293, 7 ± 0, 1 [K]

Ajuste de Redlich-Kwong

Figura 2: Isotermas medidas y la región de consistencia

A partir de los valores presentados en la figura 2 podemos

estimar las coordenadas críticas como:

P(PC ;VC ; TC ) ≈ P(37, 5; 1; 318, 8) (3)

Además de usar la gráfica, se observó que cuando se lle-

gaba al punto crítico determinado por el marco teórico y los
valores tabulados de las constantes aparecía una región de
color más blanquecino como se muestra en la figura 3.

Figura 5: Ajuste no lineal de la isoterma de T = 293, 7 ± 0, 1 [K]

1. Bondad del ajuste

A continuación comentaremos los resultados basándonos
en el χ 2 de Pearson este test estadístico nos permite observar
la correspondencia entre los datos y el modelo empleado, es
decir, nos permite ver la calidad del modelo de ajuste (junto
con el R2 ) empleado.
En la siguiente tabla podemos observar los valores de χ 2
de las diferentes isotermas según el modelo empleado para
Figura 3: Momento en el que aparece el punto de coordenadas
críticas
su ajuste:
 PRÁCTICA DE LABORATORIO, VOL. 1, NO. 3, Abril 2023

T ± ∆T [K] Ecuación de Van Der Waals Ecuación Redlich-Kwong T ± ∆T [K] (PC ± ∆PC ) · 105 [Pa] (VC ± ∆VC ) · 106 [m3 ] TC ± ∆TC [K]
293, 7 ± 0, 1 58,08 12,26 293, 7 ± 0, 1
303, 0 ± 0, 1 3,44 8,71 303, 0 ± 0, 1
308, 7 ± 0, 1 3,93 1,21 308, 7 ± 0, 1
313, 0 ± 0, 1 6,52 1,17 313, 0 ± 0, 1
318, 8 ± 0, 1 4,80 1,99 318, 8 ± 0, 1
322, 1 ± 0, 1 8,91 1,53 322, 1 ± 0, 1
Media
TABLA 1: VALORES DEL χ̄ 2 DE P EARSON DE LAS ISOTERMAS
SEGÚN LA ECUACIÓN DE AJUSTE EMPLEADA TABLA 5: VALORES OBTENIDOS PARA CADA ISOTERMA Y SU
MEDIA SEGÚN EL AJUSTE A LA ECUACIÓN DE
R EDLICH -K WONG
Como se puede observar para la temperatura más baja,
aquella que más se sale del marco de los gases ideales, el χ 2 d. A partir de las coordenadas críticas, calcule el fac-
es el más alto. Esto se debe principalmente a que las ecuacio- PCVC
tor de compresibilidad crítica ZC = nRT .Compare
nes empleadas para los modelos son correcciones a la ecua- C

ción de los gases ideales.

el valor tabulado con los obtenidos a partir de las
ecuaciones de van der Waals y de Redlich-Kwong.
Además se puede observar que los valores de χ 2 son ma-
yores en el modelo de Van Der Waals, implicando que este En función de los resultados, discuta la validez de
modelo se ajusta peor a los datos experimentales. las ecuaciones de estado propuestas.
A continuación mostramos el valor del factor de compre-
2. Determinación de los parámetros de las ecuaciones sibilidad crítica (ZC )* .
A partir de los datos obtenidos según el marco teórico re-
A continuación se muestran dos tablas con los parámetros gido por la ecuación de Van Der Waals (Tabla 4) el factor de
obtenidos a raíz de los ajustes de ambas ecuaciones: comprensibilidad crítica es:
T ± ∆T [K] (n ± ∆n) · 10−3 [mol] (a ± ∆a) · 10−1 [K] (b ± ∆b) · 10−5 [K] ZC = 0, 375 ± 0, 033 (4)
293, 7 ± 0, 1 3, 36 ± 0, 17 7, 03 ± 0, 85 7, 72 ± 0, 42
303, 0 ± 0, 1 3, 58 ± 0, 03 7, 19 ± 0, 16 7, 73 ± 0, 08
308, 7 ± 0, 1 3, 40 ± 0, 04 7, 30 ± 0, 22 7, 94 ± 0, 11 Haciendo uso de los valores obtenidos mediante el ajus-
313, 0 ± 0, 1 3, 56 ± 0, 04 7, 25 ± 0, 23 7, 83 ± 0, 12 te a la ecuación de Redlich-Kwong (Tabla 5) obtenemos el
318, 8 ± 0, 1 3, 67 ± 0, 04 7, 59 ± 0, 25 8, 22 ± 0, 13 siguiente valor para el factor de comprensibilidad crítica:
322, 1 ± 0, 1 3, 60 ± 0, 05 7, 42 ± 0, 32 8, 14 ± 0, 17
Media 3, 53 ± 0, 06 7, 30 ± 0, 34 7, 93 ± 0, 17
ZC = 0, 375 ± 0, 033 (5)
TABLA 2: VALORES OBTENIDOS PARA CADA ISOTERMA Y SU
VAN D ER WAALS
MEDIA SEGÚN EL AJUSTE A LA ECUACIÓN DE IV. C ONCLUSIÓN
R EFERENCIAS
T ± ∆T [K] (n ± ∆n) · 10−3 [mol] (a ± ∆a) · 10−1 [K] (b ± ∆b) · 10−5 [K] [1] J. R. Taylor, An introduction to error analysis: The study of uncertain-
293, 7 ± 0, 1 tiesin Physcal Measurements, 2nd ed. University Science Books, 1997.
303, 0 ± 0, 1
308, 7 ± 0, 1
313, 0 ± 0, 1 A. A NEXO
318, 8 ± 0, 1
322, 1 ± 0, 1 a. Cálculo de errores
Media
En este apartado comentaremos los diferentes errores de me-
TABLA 3: VALORES OBTENIDOS PARA CADA ISOTERMA Y SU dición que hemos tenido en cuenta y el porque de la existen-
MEDIA SEGÚN EL AJUSTE A LA ECUACIÓN DE cia de los mismos.[1]
R EDLICH -K WONG
• Errores de la instrumentación
c. Determine a partir de los parámetros del aparta- Los errores de las instrumentación son aquellos erro-
do anterior las coordenadas del punto crítico pre- res que el fabricante de los aparatos de medición esti-
dichas por las ecuaciones de van der Waals y de pula con el objetivo de fijar unos rangos de precisión en
Redlich-Kwong. Compare los resultados con la esti- cuanto a la fiabilidad de las medidas.
mación del punto crítico experimental y con el pun- Los errores de instrumentación son los siguientes:
to crítico tabulado para el gas SF6 .
1. Error del termómetro
A continuación se muestran los valores para la obtención El error del termómetro de mercurio, y por tanto
de las coordenadas del punto crítico en función de los valo- de la temperatura, es de 0, 1 [◦C]
res obtenidos haciendo uso de los valores obtenidos de los
coeficientes a y b (Véase tablas 2 t 3) : 2. Error del barómetro
El error del barómetro tomado es 0, 25 [bar]
(PC ± ∆PC ) · 105 [Pa] (VC ± ∆VC ) · 106 [m3 ] TC ± ∆TC [K] 3. Error del volumen El error del volumen es de
Media 0, 25 [ml]
* Como sugiere el guión de la práctica de laboratorio se ha usado el valor
TABLA 4: VALORES OBTENIDOS PARA CADA ISOTERMA Y SU
medio para realizar los cálculos
MEDIA SEGÚN EL AJUSTE A LA ECUACIÓN DE VAN D ER WAALS
 TERMÓMETROS

• Propagación de errores
Los errores propagados se dan debido al uso de fórmu-
las que hacen que el error se propague hasta el resul-
tado final en el tratamiento de las medidas. Dado que
las temperatura es un factor a tener en cuenta no pode-
mos asegurar la máxima precisión en todos los datos.
Por lo que se ha optado por propagar los errores siendo
consecuentes con las posibles variaciones a lo largo del
experimento; de forma que la propagación de errores se
ha llevado a cabo como se muestra a continuación:
n
∂ f (x1 , x2 , . . . , xn
∆ f (x1 , x2 , . . . , xn = ∑ ∆xi (6)
i=1 ∂ xi

• Error presión crítica Redlich-Kwong

Figura 6: Ajuste no lineal de la isoterma de T = ±0, 1 [K]
v
u Ra2 0,086645  23
u


t 2 5 · 4∆a2 b2 +25∆b2 a2
0,42748 b
a2 b2
∆pc = (7)
3

• Error volumen crítico Redlich-Kwong

q
∆b2 n2 +∆n2 b2
0,086642
∆Vc = (8)
3

• Error temperatura crítica Redlich-Kwong

s
2
a2 0,086642

3

∆a2 b2 +∆b2 a2
0,427482 R2 b2
2 a2 b2
∆Tc = (9)
3

• Error presión crítica Van Der Waals

q
∆a2 b2 +4∆b2 a2
b6
∆pc = (10)
27
Figura 7: Ajuste no lineal de la isoterma de T = ±0, 1 [K]
• Error volumen crítico Van Der Waals
p
∆Vc = 3 ∆b2 n2 + ∆n2 b2 (11)

• Error temperatura crítica Van Der Waals

q
2 2 2 a2
8 ∆a bR2+∆b
b4
∆Tc = (12)
27

• Error factor de compresibilidad crítica

s
Pc2 Tc2Vc2 ∆n2 + Pc2 Tc2 ∆Vc2 n2 + Pc2Vc2 ∆Tc2 n2 + Tc2Vc2 ∆Pc2 n2
∆Zc−=
R2 Tc4 n4

(13)

b. Ajuste de isotermas
En este apartado mostraremos el resto de ajustes de las iso-
termas medidas en el laboratorio.

Ajuste a Van Der Waals Figura 8: Ajuste no lineal de la isoterma de T = ±0, 1 [K]
 PRÁCTICA DE LABORATORIO, VOL. 1, NO. 3, Abril 2023

Figura 9: Ajuste no lineal de la isoterma de T = ±0, 1 [K] Figura 12: Ajuste no lineal de la isoterma de T = ±0, 1 [K]

Figura 10: Ajuste no lineal de la isoterma de T = ±0, 1 [K]

Figura 13: Ajuste no lineal de la isoterma de T = ±0, 1 [K]

Ajuste a Redlich-Kwong

Figura 11: Ajuste no lineal de la isoterma de T = ±0, 1 [K] Figura 14: Ajuste no lineal de la isoterma de T = ±0, 1 [K]
 TERMÓMETROS

P ± ∆P [bar] V ± ∆V [ml]
27, 50 ± 0, 25 2, 300 ± 0, 025
29, 00 ± 0, 25 2, 100 ± 0, 025
29, 50 ± 0, 25 1, 950 ± 0, 025
29, 50 ± 0, 25 1, 900 ± 0, 025
30, 00 ± 0, 25 1, 325 ± 0, 025
30, 00 ± 0, 25 1, 200 ± 0, 025
30, 00 ± 0, 25 1, 000 ± 0, 025
30, 00 ± 0, 25 0, 900 ± 0, 025
30, 00 ± 0, 25 0, 700 ± 0, 025
32, 50 ± 0, 25 0, 550 ± 0, 025
37, 00 ± 0, 25 0, 520 ± 0, 025
41, 50 ± 0, 25 0, 500 ± 0, 025

TABLA 8: DATOS DE LA ISOTERMA DE T = ±0, 1 [K]

Figura 15: Ajuste no lineal de la isoterma de T = ±0, 1 [K]

c. Datos

En este apartado mostramos todas las tablas con datos medi-

dos en el laboratorio (se excluyen aquellos datos que no se
obtengan en el proceso de medición que se expondrán en su
correspondiente apartado) con los errores, los cuales siguen P ± ∆P [bar] V ± ∆V [ml]
las fórmulas obtenidas en el Anexo a. 27, 00 ± 0, 25 2, 550 ± 0, 025
28, 00 ± 0, 25 2, 400 ± 0, 025
29, 00 ± 0, 25 2, 300 ± 0, 025
P ± ∆P [bar] V ± ∆V [ml] 29, 50 ± 0, 25 2, 200 ± 0, 025
18, 00 ± 0, 25 3, 900 ± 0, 025 30, 00 ± 0, 25 2, 100 ± 0, 025
19, 00 ± 0, 25 3, 600 ± 0, 025 30, 50 ± 0, 25 2, 000 ± 0, 025
20, 00 ± 0, 25 3, 350 ± 0, 025 31, 00 ± 0, 25 1, 900 ± 0, 025
20, 50 ± 0, 25 3, 200 ± 0, 025 33, 00 ± 0, 25 1, 550 ± 0, 025
21, 00 ± 0, 25 2, 900 ± 0, 025 33, 00 ± 0, 25 1, 400 ± 0, 025
21, 00 ± 0, 25 2, 700 ± 0, 025 33, 50 ± 0, 25 0, 800 ± 0, 025
21, 00 ± 0, 25 2, 500 ± 0, 025 34, 00 ± 0, 25 0, 600 ± 0, 025
21, 00 ± 0, 25 2, 300 ± 0, 025 37, 00 ± 0, 25 0, 550 ± 0, 025
21, 00 ± 0, 25 1, 600 ± 0, 025
21, 50 ± 0, 25 0, 700 ± 0, 025 TABLA 9: DATOS DE LA ISOTERMA DE T = ±0, 1 [K]
22, 00 ± 0, 25 0, 550 ± 0, 025
25, 00 ± 0, 25 0, 500 ± 0, 025
28, 50 ± 0, 25 0, 470 ± 0, 025

TABLA 6: DATOS DE LA ISOTERMA DE T = ±0, 1 [K]

P ± ∆P [bar] V ± ∆V [ml] P ± ∆P [bar] V ± ∆V [ml]

19, 00 ± 0, 25 3, 750 ± 0, 025 28, 00 ± 0, 25 2, 600 ± 0, 025
20, 50 ± 0, 25 3, 600 ± 0, 025 29, 00 ± 0, 25 2, 400 ± 0, 025
21, 00 ± 0, 25 3, 450 ± 0, 025 29, 50 ± 0, 25 2, 200 ± 0, 025
21, 50 ± 0, 25 3, 350 ± 0, 025 31, 00 ± 0, 25 2, 100 ± 0, 025
22, 00 ± 0, 25 3, 200 ± 0, 025 32, 00 ± 0, 25 2, 000 ± 0, 025
23, 50 ± 0, 25 2, 900 ± 0, 025 33, 50 ± 0, 25 1, 800 ± 0, 025
25, 00 ± 0, 25 2, 500 ± 0, 025 35, 00 ± 0, 25 1, 600 ± 0, 025
26, 50 ± 0, 25 1, 900 ± 0, 025 36, 00 ± 0, 25 1, 400 ± 0, 025
26, 50 ± 0, 25 1, 700 ± 0, 025 36, 50 ± 0, 25 1, 250 ± 0, 025
27, 00 ± 0, 25 0, 750 ± 0, 025 37, 00 ± 0, 25 1, 000 ± 0, 025
30, 00 ± 0, 25 0, 520 ± 0, 025 38, 00 ± 0, 25 0, 650 ± 0, 025
36, 00 ± 0, 25 0, 500 ± 0, 025 41, 50 ± 0, 25 0, 600 ± 0, 025

TABLA 7: DATOS DE LA ISOTERMA DE T = ±0, 1 [K] TABLA 10: DATOS DE LA ISOTERMA DE T = ±0, 1 [K]
 PRÁCTICA DE LABORATORIO, VOL. 1, NO. 3, Abril 2023

P ± ∆P [bar] V ± ∆V [ml]
29, 00 ± 0, 25 2, 500 ± 0, 025
30, 50 ± 0, 25 2, 300 ± 0, 025
33, 00 ± 0, 25 2, 000 ± 0, 025
34, 00 ± 0, 25 1, 900 ± 0, 025
35, 00 ± 0, 25 1, 700 ± 0, 025
37, 00 ± 0, 25 1, 500 ± 0, 025
37, 50 ± 0, 25 1, 450 ± 0, 025
38, 00 ± 0, 25 1, 350 ± 0, 025
38, 50 ± 0, 25 1, 250 ± 0, 025
39, 00 ± 0, 25 1, 100 ± 0, 025
40, 00 ± 0, 25 0, 950 ± 0, 025
41, 50 ± 0, 25 0, 700 ± 0, 025
43, 50 ± 0, 25 0, 650 ± 0, 025
46, 00 ± 0, 25 0, 600 ± 0, 025

TABLA 11: DATOS DE LA ISOTERMA DE T = ±0, 1 [K]