Proceso politrópico

Representamos los cuatro procesos termodinámicos

hold on
line([0,2],[1,1], 'color', 'g')
line([1,1],[0,2.5],'color','m')
fplot(@(V) 1./V,[0.25,2], 'color','b')
fplot(@(V) 1./V.^1.4,[0.25,2], 'color','r')
hold off
grid on
ylim([0,2.5])
xlabel('V')
legend('p=cte','V=cte','isotermo','adiabático')
ylabel ('p')
title('Procesos termodinámicos')

Conocida la relación entre presión y volumen en un proceso politrópico, pVβ=cte, para un gas ideal cuya ecuación de estado es pV=nRT

Representamos procesos politrópicos de índices: 0.2, 0.5, 3, 10, 100, -10, -2, -1, -0.5, 0

hold on
for n=[0.2,0.5,3,10,100, -10,-2,-1,-0.5,0]
    fplot(@(V) 1./V.^n,[0.25,2])
end
n=1.4;
fplot(@(V) 1./V.^n,[0.25,2],'color', 'r','lineWidth',1.2)
n=1;
fplot(@(V) 1./V.^n,[0.25,2],'color', 'b','lineWidth',1.2)
hold off
grid on
ylim([0,2.5])
xlabel('V')
ylabel ('p')
title('Procesos politrópicos')

Los procesos con β>0 se representan en los cuadrantes 2° y 4°. Los procesos con β<0 se representan en los cuadrantes 1°y 3°. β=0 es un proceso a presión constante (línea horizontal), y β=100 es un proceso a volumen casi constante (línea vertical)

Las dos líneas más gruesas: la azul, es un proceso con β=1 (isotermo), la roja, es un proceso con β=1.4 (adiabático)

Trabajo

Se cumple que

p V β = p 1 V 1 β = p 2 V 2 β =cte

El trabajo vale

W= V 1 V 2 p·dV = V 1 V 2 cte V β dV= cte β+1 ( V 2 β+1 V 1 β+1 )= 1 β+1 ( p 2 V 2 β V 2 β+1 p 1 V 1 β V 1 β+1 ) W= p 2 V 2 p 1 V 1 β+1

Otra forma de expresar el trabajo

W= p 1 V 1 β β+1 ( V 2 β+1 V 1 β+1 )= p 1 V 1 β V 1 β+1 β+1 ( ( V 2 V 1 ) β+1 1 )= p 1 V 1 β+1 ( ( V 2 V 1 ) β+1 1 )

Teniendo en cuenta la relación entre los calores específicos a presión constante cp y a volumen constante cv

cp=cv+R, γ=cp/cv, cv=R/(γ-1)

T 1 V 1 β1 = T 2 V 2 β1 , T 1 V 1 β+1 = T 2 V 2 β+1 , ( V 2 V 1 ) β+1 = T 2 T 1 W= nR T 1 β+1 ( T 2 T 1 1 )= n( c p c v ) β+1 ( T 2 T 1 )=n c v γ1 β+1 ( T 2 T 1 )

Variación de energía interna

La variación de energía interna entre los estados inicial y final es

ΔU=n c v ( T 2 T 1 )

Calor

El primer principio de la Termodinámica, ΔU=Q-W, Q=ΔU+W

Q=n c v ( T 2 T 1 )+n c v γ1 β+1 ( T 2 T 1 )=n c v βγ β1 ( T 2 T 1 )=n c β ( T 2 T 1 ) c β = c v βγ β1 = β c v c p β1 = R γ1 βγ β1

Es la definición de calor específico cβ en un proceso politrópico

Representamos cβ/R en función de β, tomando el índice addibático γ=1.4

gamma=1.4;
f=@(x) (x-gamma)./((gamma-1)*(x-1));
hold on
fplot(f, [-1,3])
plot(0,7/2,'ro','markersize',3,'markerfacecolor','r')
line([-1,0],[7/2,7/2],'lineStyle','--')
line([0,0],[0,7/2],'lineStyle','--')
line([-1,3],[5/2,5/2],'lineStyle','--')
plot(gamma,0,'ro','markersize',3,'markerfacecolor','r')
hold off
grid on
ylim([-2,10])
xlabel('\beta')
ylabel ('c_\beta/R')
title('Calor específico c_\beta')

El calor específico cβ es negativo en el intervalo (1,γ)

Deducción alternativa de proceso politrópico

Hemos partido de la definición de proceso politrópico como pVβ=cte. Partimos ahora de, dQ=ncβdT y llegamos a pVβ=cte, del siguiente modo

dU=dQpdV n c v dT=n c β dTpdV n( c β c v )dT= nRT V dV dT T = R c β c v dV V dT T c p c v c β c v dV V =0 T V c p c v c β c v =cte p V 1 c p c v c β c v =cte p V c β c p c β c v =cte β= c β c p c β c v , c β = β c v c p β1

Relación entre trabajo y calor

Relacionamos el calor Q y el trabajo W

Q=n c v βγ β1 ( T 2 T 1 ) W=n c v γ1 β+1 ( T 2 T 1 ) W Q = γ1 γβ

Representamos el cociente W/Q en función de β tomando γ=1.4 (aire)

gamma=1.4;
f=@(x) (gamma-1)./(gamma-x);
fplot(f, [-5,5])
line([1,1],[0,f(1)],'lineStyle','--')
line([-5,1],[f(1),f(1)],'lineStyle','--')
line([0,0],[0,f(0)],'lineStyle','--')
line([-5,0],[f(0),f(0)],'lineStyle','--')
grid on
ylim([-2.5,2.5])
xlabel('\beta')
ylabel ('W/Q')
title('Procesos politrópicos')

Para un proceso termodinámico caracterizado por un índice β<γ, el resultado es W/Q>0, parte del calor absorbido Q>0 se convierte en trabajo W>0 (el gas se expande). Si el gas se comprime W<0 el gas cede calor Q<0 a los alrededores del sistema.

Indice politrópico

Cuando se comprime el aire muy rápidamente, no hay tiempo de que el sistema transfiera calor a los alrededores, el aire incrementa su temperatura, el proceso es adiabático β=γ

Cuando se comprime aire muy lentamente, se transfiere calor a los alrededores, la temperatura del aire permanece constante, el proceso es isotermo, β=1

Un proceso real no será ni infinitamente rápido, ni infinitamente lento, el proceso politrópico tendrá un índice β comprendido entre 1 y γ.

Dados el estado inicial y final, el índice del proceso politrópico se calcula despejando β.

Cuando 1<β<γ el calor específico cβ<0. La idea de calor específico negativo parace algo extraña, ya que la tempartaura del gas ΔT>0 se incrementa a la vez que el gas cede calor Q<0 o viceversa

Procesos isotermo, adiabático y politrópico

Estado inicial del gas

p1=17.5 atm, V1=0.2124 l, T1=2000 K

Estado final

V1=1 l

Representamos los tres procesos: isotermo, adiabático y politrópico entre el estado incial y el estado final. Calculamos la presión p2 y la temperatura T2 del estado final

Calculamos el calor Q, trabajo W y variación de energía interna, ΔU en cada uno de los procesos

gamma=1.4;
beta=1.1689;
p1=17.5; %estado inicial
V1=0.2124;
T1=2000;
V2=1; %estado final
hold on
fplot(@(x) p1*V1^gamma./x.^gamma, [V1,V2])
fplot(@(x) p1*V1^beta./x.^beta, [V1,V2])
fplot(@(x) p1*V1./x, [V1,V2])
hold off
grid on
xlabel('V (l)')
legend('adiabático','politrópico','isotermo','location','best')
ylabel ('p(atm)')
title('Adiabático, isotermo, politrópico')

%isotermo
p2=p1*V1/V2;
T2=2000;
Du=0;
W=p1*V1*log(V2/V1);
Q=W;
fprintf('isotermo: p2=%1.3f, T2=%4.1f, Du=%1.3f, Q=%1.3f, W=%1.3f\n',
 p2, T2, Du, Q, W)
%adiabático
p2=p1*V1^gamma/V2^gamma;
T2=T1*V1^(gamma-1)/V2^(gamma-1);
Du=5*p1*V1*(T2-T1)/(2*T1);
Q=0;
W=(p2*V2-p1*V1)/(-gamma+1);
fprintf('adiabático: p2=%1.3f, T2=%4.1f, Du=%1.3f, Q=%1.3f, W=%1.3f\n', 
p2, T2, Du, Q, W)
%politrópico
p2=p1*V1^beta/V2^beta;
T2=T1*V1^(beta-1)/V2^(beta-1);
Du=5*p1*V1*(T2-T1)/(2*T1);
Q=p1*V1*(beta-gamma)*(T2-T1)/(T1*(gamma-1)*(beta-1));
W=(p2*V2-p1*V1)/(-beta+1);
fprintf('politrópico: p2=%1.3f, T2=%4.1f, Du=%1.3f, Q=%1.3f, W=%1.3f\n', 
p2, T2, Du, Q, W)

isotermo: p2=3.717, T2=2000.0, Du=0.000, Q=5.759, W=5.759
adiabático: p2=2.000, T2=1076.2, Du=-4.292, Q=0.000, W=4.292
politrópico: p2=2.861, T2=1539.5, Du=-2.139, Q=2.927, W=5.067

Expansión rápida de una mezcla de aire con vapor saturado

Un recipiente de volumen V0 contiene aire y vapor saturado a la temperatura T0. El émbolo ejerce una presión p0.

El número de moles de vapor saturado es

n v0 = P v0 V 0 R T 0

Donde MH2O=18 g/mol es la masa molecular del agua y R=8.3143 J/(K·mol) es la constante de los gases.

Supondremos que el vapor se comporta como un gas ideal cuya presión parcial Pv0 se puede consultar en las tablas o deducir de fórmulas empíricas de ajuste de los datos experimentales. Por ejemplo, para T0=19 °C (292 K) la presión de vapor saturado es Pv0=2 198 Pa

Cuando el émbolo se desplaza rápidamente. El volumen de aire aumenta a V y la temperatura disminuye de T0 a T=T0-ΔT, con ΔT<<T0

Al disminuir la temperatura de la mezcla, la presión de vapor Pv se reduce, una parte del vapor inicial nv se condensa.

n v = P v V RT

El número de moles de vapor que se ha condensado es la diferencia

n= P v0 V 0 R T 0 P v V RT

El calor generado en la condensación será

Q=nL= L R ( P v0 V 0 T 0 P v V T )

Donde L=43 509 J/mol es el valor de la entalpía de vaporización

Como el proceso de expansión es rápido, la mezcla aire-vapor no tiene tiempo de intercambiar calor con el medio ambiente a través de las paredes del recipiente

La presión de vapor viene dada por la fórmula de Clausius-Clapeyron

P v = P v 0 exp( L R ( 1 T 0 1 T ) )

Sustituimos la presión parcial de vapor Pv en la expresión del calor generado en la condensación

Q= L R P v0 ( V 0 T 0 exp( L R ( 1 T 0 1 T ) ) V T )

Proceso politrópico

Entre el estado inicial y final se produce un proceso politrópico de índice β que vamos a determinar

p V β = p 0 V 0 β p ( nRT p ) β = p 0 ( nR T 0 p 0 ) β p 1β T β = p 0 1β T 0 β p p 0 = ( T 0 T ) β/(1β) V T = nR p = nR p 0 ( T 0 T ) β/(β1)

El número n de moles de la mezcla (aire+vapor) no permence constante ya que una parte del vapor inicial se condensa. Si esta parte es pequeña, en comparación con el número total de moles, se podría considerar que n es casi constante.

En un proceso politrópico, el calor Q es

Q=n c β ( T T 0 )=n c β ΔT c β = R γ1 βγ β1

Determinación del índice β

Introduciendo V/T, la expresión del calor generado en la condensación es

Q= L R P v0 ( nR p 0 exp( L R ( 1 T 0 1 T ) ) nR p 0 ( T 0 T ) β/(β1) ) Q=n P v0 p 0 L( 1exp( L R ( 1 T 0 1 T ) ) ( T 0 T ) β/(β1) ) Q=n P v0 p 0 L( 1exp( L R T 0 ( 1 T 0 T 0 ΔT ) ) ( T 0 T 0 ΔT ) β/(β1) ) Q=n P v0 p 0 L( 1exp( L R T 0 ( x 1x ) ) ( 1x ) β/(β1) ),x= ΔT T 0

Dado que las temperaturas inicial y final son próximas, T=T0-ΔT, con ΔT<<T0, desarrollamos en serie tomando el término en ΔT

Utilizamos la función taylor de MATLAB, se omite la constante delante del paréntesis, las variables x es ΔT/T0, b es β y a es L/(RT0)

>> syms x b a;
>> y=1-exp(-a*x/(1-x))*(1-x)^(-b/(b-1));
>> taylor(y,x,0,'order',2)
ans =x*(a - b/(b - 1))

Qn{ P v0 p 0 L T 0 ( L R T 0 β β1 ) }ΔT

Comparamos con el calor del proceso politrópico

Q=n{ R γ1 βγ β1 }ΔT

Despejamos el índice β

P v0 p 0 L T 0 ( L R T 0 β β1 )= R γ1 βγ β1 β= ( L R T 0 ) 2 P v0 p 0 ( γ1 )+γ ( L R T 0 ) P v0 p 0 ( γ1 )( L R T 0 1 )+1

Si no hay vapor o la presión es muy pequeña, Pv0≈0, entonces β=γ el proceso es adiabático

Creamos un script para calcular el índice β dada la presión inicial p0 en atmósferas. Datos

gamma=1.4; %índice adiabático
L=43509; %calor de vaporización J/mol
T0=292; %temperatura inicial
Pv=2.1983e3;%presión de vapor
R=8.3143; %constante de los gases 
for p=[4.8,4.4,4.0,3.6,3.2] %presión en atmósferas
    P0=p*1.013e5;
    beta=(Pv*L^2*(gamma-1)/(R^2*T0^2*P0)+gamma)/(1+Pv*L*(gamma-1)*
(L/(R*T0)-1)/(R*P0*T0));
    disp([p,beta])
end
    4.8000    1.2793
    4.4000    1.2724
    4.0000    1.2647
    3.6000    1.2560
    3.2000    1.2461
Presión, p0 (atm)índice, β
4.81.2793
4.41.2724
4.01.2647
3.61.2560
3.21.2461

Son similares a la tercera columna de la tabla I del artículo de Romanelli, pág. 764

Fórmula barométrica

En la página titulada 'Modelos simples de atmósfera', calculamos la variación de la presión del aire en función de la altura por encima del nivel del mar. Suponíamos que la atmósfera está en equilibrio y que su temperatura se mantenía constante, pV=cte.

En este caso, la relación entre presión y volumen responde al la fórmula pVβ=cte

Las fuerzas que actúan sobre el elemento de fluido (aire) de densidad ρ, un cilindro de base A y altura dz

Las fuerzas que ejerce la presión sobre la superficie lateral se anulan

En el equilibrio

(p+dp)A+ρ(A·dz)g=pA dp=ρg·dz

Supongamos ahora que la relación entre presión y volumen es de la forma

p V β = p 0 V 0 β p ρ β = p 0 ρ 0 β ρ= ρ 0 ( p p 0 ) 1/β

Donde p0 y ρ0 son la presión y densidad a la altura z=0

Integramos la ecuación diferencial

p= ρ 0 ( p p 0 ) 1/β g·dz p 0 p dp p 1/β = ρ 0 p 0 1/β g 0 z dz 1 1 β +1 ( p 1/β+1 p 0 1/β+1 )= ρ 0 p 0 1/β gz z= β β1 1 ρ 0 g ( p 0 p β1 β p 0 1/β )= β β1 p 0 ρ 0 g ( 1 ( p p 0 ) β1 β )

Utilizando la ecuación de estado del un gas ideal

pV= m M RT p= ρ M RT

El resultado es

z= β β1 ρ 0 M R T 0 ρ 0 g ( 1 ( p p 0 ) β1 β )= β β1 R T 0 Mg ( 1 ( p p 0 ) β1 β )

Despejamos la presión

p= p 0 ( 1 β1 β Mgz R T 0 ) β β1

Llamamos

x= β1 β Mgz R T 0 p= p 0 { ( 1+x ) 1 x } Mgz R T 0

Cuando β→1 (proceso isotermo), x→0, (1+x)1/x→e (número e)

>> syms x;
>> limit((1+x)^(1/x),x,0)
ans =exp(1)

La variación de presión p con la altura z que obtuvimos para una atmósfera isoterma

p= p 0 exp( Mgz R T 0 )

Datos

R=8.3143; %constante de los gases
T0=300; %temperatura a nivel del mar
M=28.9/1000; %masa molecular del aire en kg/mol
hold on
f=@(z) exp(-M*9.8*1000*z/(R*T0)); %n=1
fplot(f,[0,40])
for n=1.1:0.1:1.4
    f=@(z) (1-(n-1)*M*9.8*z*1000/(n*R*T0)).^(n/(n-1));
    fplot(f,[0,40])
end
grid on
xlabel('z (km)')
ylabel('p/p_0')
title('Fórmula barométrica')

Referencias

Polytropic process in a closed system. tec-science (07/11/2021)

Alejandro Romanelli, Italo Bove, Federico González Madina. Air expansion in a water rocket. Am. J. Phys. 81 (10), October 2013, pp. 762-766

W EI Guo-zhu,GU Yong-wei,DU An. The improvement on the barometric formula in isothermal process of the gas. College Physics, Volume 25, Issue 11, 20 November 2006