>	restart:#"m02_p12"

Dentro de un cilindro vertical de 0,01 m2 de sección hay 0,01 kg de nitrógeno encerrado con un émbolo superior de acero de baja fricción cuyo peso da lugar a una sobrepresión de 5 kPa, estando el ambiente a 20 ºC y 93 kPa. Considérese la siguiente evolución: 1) mediante las fuerzas apropiadas se obliga al pistón a reducir lentamente en un 10% el volumen ocupado por el gas, y 2) se libera el anclaje y se permite el libre movimiento del émbolo. Se pide:
a) Esquema de la evolución en un diagrama altura-tiempo.
b) Determinar el espesor del émbolo y su altura inicial.
c) Valores p-V-T en los estados de equilibrio considerados y diagrama p-V.
d) Variación de energía entre los estados antedichos y para todo el ciclo, así como calor y trabajo transferidos entre los sistemas involucrados.
e) Diagrama T-s de la evolución.
f) Variación de entropía y generación de entropía para todos los sistemas entre los estados antedichos y para todo el ciclo.

Datos:

>	read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):unprotect(gamma):

>	su1:="N2":su2:="Acero_inox":dat:=[A=0.01*m_^2,m=0.01*kg_,DpE=5e3*Pa_,T0=293*K_,p0=93e3*Pa_,V2_V1=0.9];

Esquema:

>

>

Ecs. const.:

>	eqET:=subs(eq1_11,eq1_12);eqEE:=eq1_16;gdat:=get_gas_data(su1):sdat:=get_sol_data(su2):rho[acero]:=subs(sdat,rho);dat:=op(dat),Const,gdat,sdat,SI2,SI1:

a) Determinar el espesor del émbolo y su altura inicial.

>

eqBF:=mE*diff(z(t),t,t)=(p-p0)*A-mE*g-Ff;eqBFint:=int(mE*diff(z(t),t,t)*diff(z(t),t),t)=Int((p-p0)*A-mE*g-Ff,z);T1:=T0:T2:=T0:p1:=p0+DpE;p1_:=subs(dat,p1):'p1'=evalf(%/(1000*Pa_/kPa_),2);V1:=solve(subs(p=p1,T=T1,eqET),V);V1_:=subs(dat,subs(p1=p1_,dat,V1)):'V1'=evalf(V1_,2);z1_:=subs(dat,V1_/A):'z1'=evalf(z1_,3);e_:=solve(DpE=rho*e*A*g/A,e);e__:=subs(dat,e_):'e'=evalf(e__,2);mE_:=subs(e=e_,dat,rho*e*A):'mE'=evalf(mE_,2);

b) Valores p-V-T en los estados de equilibrio considerados, numerándolos convenientemente.

>

V2_:=subs(dat,dat,V1*V2_V1):'V2'=evalf(V2_,2);p2:=solve(subs(T=T2,V=V2_,eqET),p):p2_:=subs(dat,dat,p2):'p2'=evalf(p2_/(1000*Pa_/kPa_));z2_:=subs(dat,V2_/A):'z2'=evalf(z2_,3);DU12:=0;W12:=Int(-p,V='V1..V2');W12:='Int(-m*R*T0/V,V=V1_..V2_)';W12_:=evalf(subs(dat,SI0,value(W12)))*J_:'W12'=evalf(W12_,2);Q12:='DU12-W12';Q12_:=DU12-W12_:'Q12'=evalf(Q12_,2);Fdedo:='(m*R*T0/V-p0)*A-mE*g';Wdedo:='Int(-Fdedo/A,V=V1..V2)';Wdedo_:=subs(dat,int(expand(subs(mE=mE_,dat,V=V*m_^3,-Fdedo/A*m_^3)),V=V1_/m_^3..V2_/m_^3)):'Wdedo'=evalf(Wdedo_,2);Watm_:=subs(mE=mE_,dat,p0*A*(z1_-z2_)):'Watm'=evalf(Watm_,2);WE_:=subs(mE=mE_,dat,mE*g*(z1_-z2_)):'WE'=evalf(WE_,2);Wneto_:=W12_-Watm_-WE_:'Wneto'=evalf(Wneto_,2);

b) cont.) Suelta desde 2 hasta 3=max, 4=equil.mec. y 5=reposo

>

eqBE23:=0=rhs(eqBFint);eqBE23:=m*c[v]*(T3-T2)=-p0*A*(z3-z2)-mE*g*(z3-z2);eqBE23_:=subs(T3=T2*(z2/z3)^(gamma-1),mE=mE_,z2=z2_,dat,SI0,eqBE23):z3_:=fsolve(eqBE23_,z3=z2_/m_+1e-6..1)*m_:'z3'=evalf(z3_,3);T3_:=subs(dat,T2*(z2_/z3_)^(gamma-1));V3_:=subs(dat,z3_*A);p3_:=subs(dat,m*R*T3_/V3_):'p3'=evalf(p3_/(1000*Pa_/kPa_));Q23_:=0;W23_:=subs(dat,m*c[v]*(T3_-T2));

>	p4_:=subs(dat,p1):'p4'=evalf(%/(1000*Pa_/kPa_));z4_:=subs(dat,dat,z2_*(p2/p4_)^(1/gamma)):'z4'=evalf(z4_,3);T4_:=subs(dat,T2*(z2_/z4_)^(gamma-1));V4_:=subs(dat,z4_*A);Q24_:=0;W24_:=subs(dat,m*c[v]*(T4_-T2));

p5_:=subs(dat,p1):'p5'=evalf(%/(1000*Pa_/kPa_));z5_:=z1_;T5_:=subs(dat,T1);V5_:=subs(dat,z5_*A);W45_:=subs(dat,-p5_*(V5_-V4_));Q45_:=subs(dat,m*c[v]*(T5_-T4_)-W45_);

c) Representación de la evolución en el diagrama p-V.

>	plot({[[V1_/m_^3,p1_/Pa_],[V2_/m_^3,p2_/Pa_]],p1_/Pa_*(V1_/m_^3/V_),[[V4_/m_^3,p4_/Pa_],[V5_/m_^3,p5_/Pa_]],[[.9*V2_/m_^3,p3_/Pa_],[V3_/m_^3,p3_/Pa_]],p2_/Pa_*(V2_/m_^3/V_)^subs(dat,gamma)},V_=.9*V2_/m_^3..1.1*V3_/m_^3,p=0.9*p3_/Pa_..1.1*p2_/Pa_,colour=black);

e) Diagrama T-s de la evolución.

Bastaría un esquema, pero se va a hacer a escala.

>

eqDS:=DS12='m*(c[p]*ln(T2/T1)-R*ln(p2/p1))';DS12_:=subs(dat,m*(c[p]*ln(T2/T1)-R*ln(p2_/p1_))):DS12=evalf(%,3);DS24_:=subs(dat,dat,evalf(subs(dat,m*(c[p]*ln(T4_/T2)-R*ln(p4_/p2_))))):DS24=evalf(%,3);DS45_:=subs(dat,dat,evalf(subs(dat,m*(c[p]*ln(T5_/T4_)-R*ln(p5_/p4_))))):DS45=evalf(%,3);S0:=0.5*J_/K_:plot(subs(dat,SI0,{[S0+m*(c[p]*ln(T/T0)),T,T=270..300],[[S0,T0],[S0+DS12_,T2],[S0+DS12_,T3_]]}),S=0..1,T=250..300);

f) Variación de entropía y generación de entropía para todos los sistemas entre los estados antedichos y para todo el ciclo..

De 1 a 2 hay un flujo de entropía isotermo con el ambiente, así que si no hay fricción no habrá generación de entropía.

De 2 a 4 no hay flujo de entropía (se ha supuesto isoentrópica la evolución del gas), pero sí hay generación de entropía for fricción entre el émbolo y el cilindro (puede despreciarse la disipación viscosa a ambos lados del émbolo).

De 4 a 5 hay un flujo isobárico de entropía.

>	DS12amb:=Q12amb/Tamb;DS12amb_:=subs(dat,-Q12_/T0):'DS12amb'=evalf(%,3);DS24amb:=Q24amb/Tamb;DS24amb:=-W24/Tamb;DS24amb_:=subs(dat,-W24_/T0):'DS24amb'=evalf(%,3);DS45amb:=Q45amb/Tamb;DS45amb_:=subs(dat,-Q45_/T0):'DS45amb'=evalf(%,3);

>

Sgen:='DS+DSamb';Sgen12_:=DS12_+DS12amb_:'Sgen12'=evalf(%,2);Sgen24_:=DS24_+DS24amb_:'Sgen24'=evalf(%,3);Sgen45_:=DS45_+DS45amb_:'Sgen45'=evalf(%,2);DS15:='DS12+DS24+DS45';DS15_:=DS12_+DS24_+DS45_:'DS15'=evalf(%,3);DS15amb:='DS12amb+DS24amb+DS45amb';DS15amb_:=DS12amb_+DS24amb_+DS45amb_:'DS15amb'=evalf(%,3);Sgen15:='Sgen12+Sgen24+Sgen45';Sgen15_:=Sgen12_+Sgen24_+Sgen45_:'Sgen15'=evalf(%,3);

viéndose que casi toda la entropía se genera de 2 a 4 (tan sólo un 2% se genera de 4 a 5)

>