lunes, 8 de mayo de 2017

Termodinámica

Ecuaciones de la termodinámica


La ecuación de Van der Waals es una ecuación de estado de un fluido compuesto de partículas con un tamaño no despreciable y con fuerzas intermoleculares, como las fuerzas de Van der Waals. La ecuación, cuyo origen se remonta a 1873, debe su nombre a Johannes Diderik van der Waals, quien recibió el premio Nobel en 1910 por su trabajo en la ecuación de estado para gases y líquidos, la cual está basada en una modificación de la ley de los gases ideales para que se aproxime de manera más precisa al comportamiento de los gases reales al tener en cuenta su tamaño no nulo y la atracción entre sus partículas.


Ecuación

Las isotermas de Van der Waals: el modelo predice correctamente la fase de líquido prácticamente incompresible, pero las oscilaciones en la fase de transición no se corresponden con los resultados experimentales.
Una forma de esta ecuación es:
donde:
p es la presión del fluido, medido en atmósferas,
v es el volumen en el que se encuentran las partículas dividido por el número de partículas (en litros),
k es la constante de Boltzmann,
T es la temperatura, en kelvin,
a' es un término que tiene que ver con la atracción entre partículas,
b' es el volumen medio excluido de v por cada partícula.
Si se introducen el número de AvogadroNA, el número de moles n y, consecuentemente, el número total de partículas n•NA, la ecuación queda en la forma siguiente:
donde:
p es la presión del fluido,
V es el volumen total del recipiente en que se encuentra el fluido,
a mide la atracción entre las partículas ,
b es el volumen disponible de un mol de partículas ,
n es el número de moles,
R es la constante universal de los gases ideales,
T es la temperatura, en kelvin.
Debe hacerse entre una distinción cuidadosa entre el volumen disponible para una partícula y el volumen de una partícula misma. En particular, en la primera ecuación  se refiere al espacio vacío disponible por partícula. Es decir que , es el volumen  del recipiente dividido por el número total de  de partículas. El parámetro , por el contrario, es proporcional al volumen ocupado de una partícula —únicamente delimitado por el radio radio atómico. Este es el volumen que se restará de  debido al espacio ocupado por una partícula. En la derivación original de Van der Waals, que figura a continuación,  es cuatro veces el volumen disponible de la partícula. Observe además que la presión  tiende a infinito cuando el contenedor está completamente lleno de partículas de modo que no hay espacio vacío dejado por las partículas a moverse. Esto ocurre cuando .

Validez

Por encima de la temperatura crítica la ecuación de Van der Waals es una mejora de la ley del gas ideal, y para temperaturas más bajas la ecuación es también cualitativamente razonable para el estado líquido y estado gaseoso a baja presión. Sin embargo, el modelo Van der Waals no es adecuado para los cálculos cuantitativos rigurosos, útil restante sólo con fines educativos y de calidad.
En la transición de fase de primer orden, el rango de (P, V, T), donde la fase líquida y la fase gaseosa se encuentran en equilibrio, no lo exhibe el hecho empírico de que p es constante en función de V a una temperatura dada aunque este comportamiento se pueda insertar fácilmente en el modelo la Van der Waals, el resultado ya no es un modelo analítico simple, y otros (como los basados en el principio de estados correspondientes) logran un mejor ajuste con más o menos el mismo trabajo.

Derivación

La mayoría de los libros de texto dan dos distintas derivaciones. Una de ellas es la derivación convencional que se remonta a Van der Waals y la otra es una derivación de la mecánica estadística. Este último tiene la gran ventaja de que se hace explícito el potencial intermolecular, que se descuida en la primera derivación.
El volumen excluido por partícula es , que hay que dividir por dos para no contar dos veces la misma interacción, por lo que el volumen excluido del  es , que es cuatro veces el volumen adecuado de la partícula. Fue un punto de interés para Van der Waals, que el factor de cuatro es en realidad un rendimiento superior de los valores empíricos límites,  son generalmente más bajos. Por supuesto, las moléculas no son totalmente rígidas como Van der Waals pensaba, sino con frecuencia bastante suaves.
A continuación, se introduce una fuerza atractiva por parejas entre las partículas. Van der Waals supone que, a pesar de la existencia de esta fuerza, de la densidad del fluido es homogéneo. Además se asume que el rango de la fuerza de atracción es tan pequeño que la gran mayoría de las partículas no sienten que el contenedor es de tamaño finito. Es decir, la mayor parte de ellas tienen más atracción a las partículas a su derecha que a su izquierda cuando están relativamente cerca de la pared de la izquierda del contenedor (y viceversa). Teniendo en cuenta la homogeneidad del líquido, la mayor parte de las partículas no experimentan una fuerza neta que tire de ellas hacia la derecha o hacia la izquierda. Esto es diferente para las partículas en las capas superficiales directamente adyacentes a los muros. Se sienten una fuerza neta de las partículas de mayor tirando de ellos hacia el recipiente, ya que esta fuerza no es compensado por las partículas en el lado donde la está pared (otro supuesto es que no hay interacción entre las paredes y las partículas, lo cual no es cierto como puede verse en el fenómeno de la formación de gotas; la mayoría líquidos muestran adhesión). Esta fuerza neta disminuye la fuerza ejercida sobre la pared por las partículas en la capa superficial. La fuerza neta sobre una partícula de la superficie, tirando de ella hacia el recipiente, es proporcional a la densidad numérica  (número de partículas por unidad de volumen). El número de partículas en las capas superficiales es, de nuevo, asumiendo homogeneidad del fluido, también proporcional a la densidad. En total, la fuerza sobre las paredes se reduce por un factor proporcional al cuadrado de la densidad y la presión (fuerza por unidad de superficie) se reduce en
de modo que:
Es de algún interés histórico señalar que Van der Waals en su conferencia del premio Nobel le dio crédito a Laplace argumentando que la presión se reduce de forma proporcional al cuadrado de la densidad.
Esto hace que el promedio de energía Helmholtz por partícula sea reducida en una cantidad proporcional a la densidad del fluido. Sin embargo, la presión obedece a la relación termodinámica:
donde A* es la energía Helmholtz del sistema por partícula. La atracción, por lo tanto, reduce la presión en una cantidad proporcional a . Denota la constante de proporcionalidad por a, se obtiene:
que es la ecuación de Van der Waals.

Derivaciones convencionales

Se considera en primer lugar un mol de gas que se compone de partículas puntuales sin interacción que satisfacen la ley de los gases ideales.
A continuación se asume que todas las partículas son esferas duras del mismo radio finito r (el radio de Van der Waals). El efecto del volumen finito de las partículas es disminuir el espacio vacío disponible en el cual se mueven libremente las partículas. Se debe reemplazar V por V − b, donde b se le llama el volumen excluido. La ecuación corregida se convierte en:
El volumen excluido  no es exactamente igual al volumen ocupado por el tamaño de sólidos, de tamaño finito, partículas, pero en realidad cuatro veces ese volumen. Para ver esto hay que darse cuenta de que una partícula está rodeada por una esfera de radio r = 2r (dos veces el radio original) que está prohibido para los centros de las otras partículas. Si la distancia entre dos centros de las partículas es más pequeño que 2r, lo que significaría que las dos partículas penetran entre sí, que, por definición, las esferas duras no son capaces de hacer.

Derivación termodinámica estadística

La función de partición canónica Q de un gas ideal consistente en N = nNA partículas idénticas, es
donde  es la longitud de onda térmica de De Broglie,
con las definiciones usuales: H es la constante de Planckm, la masa de una partícula, k, la constante de Boltzmann y T la temperatura absoluta. En un gas ideal  es la función de partición de una partícula en un recipiente de volumen V. Con el fin de obtener la ecuación de Van der Waals se supone ahora que cada partícula se mueve de forma independiente en un campo de potencial promedio ofrecido por las otras partículas. El promedio más las partículas es fácil, porque se supone que la densidad de las partículas del fluido de Van der Waals es homogénea. La interacción entre un par de partículas, que son esferas rígidas, se considera
r es la distancia entre los centros de las esferas y D es la distancia que las esferas duras se tocan (el doble del radio de Van der Waals). La profundidad de la de Van der Waals es .
Debido a que las partículas son independientes, la función de partición total todavía se factoriza, , pero el potencial intermolecular necesita dos modificaciones a . En primer lugar, debido al tamaño finito de las partículas, no todos los V están disponible, pero sólo , donde (como en la derivación convencional anterior). En segundo lugar, introducir un factor de Boltzmann  para atender el potencial intermolecular medio. Se divide aquí el potencial de dos, porque esta energía de interacción es compartida entre dos partículas. Así:
Toda la atracción que es ejercida sobre una partícula es:
donde se asume que en una capa de espesor dr hay N/V 4π r2dr partículas. Esta es una aproximación del campo medio, la posición de las partículas es un promedio. En realidad, la densidad cercana de la partícula es diferente que la que están lejanas, como pueden ser descritas por una función de correlación par. Además, se descuida que el fluido está encerrado entre las paredes. Realizando la integral, se obtiene:
Por lo tanto, se obtiene:
De la termodinámica estadística, se sabe que:
de modo que sólo se tienen que diferenciar los términos que contienen V. Se obtiene:

Otros parámetros termodinámicos

Se reitera que el volumen extenso V se relaciona con el volumen por partícula v=V/N  donde N = nNA  es el número de partículas en el sistema. La ecuación de estado no proporciona todos los parámetros termodinámicos del sistema. Se puede tomar la ecuación de la energía Helmholtz A:
A partir de la ecuación deducida anteriormente para lnQ, se tiene:
Esta ecuación expresa A  en términos de su variables naturales V  y T , y por lo tanto nos da toda la información sobre el sistema termodinámico. La ecuación mecánica de estado ya se deriva por encima
La ecuación de la entropía del estado de los rendimientos de la entropía (S )
de la cual se puede calcular la energía interna
Se pueden escribir ecuaciones similares para los otros potenciales termodinámicos y químicos, pero expresando cualquier potencial en función de la presión p  requerirá la solución de un polinomio de tercer orden, que produce una expresión complicada. Por lo tanto, va a ser complicado expresar la entalpía y la energía libre de Gibbs en función de sus variables naturales.

Forma reducida

A pesar de las constantes de material a y b en la forma usual de la ecuación de Van der Waals es diferente para cada líquido de una sola cuenta, la ecuación puede ser refundida en una forma invariante aplicable a todos los líquidos.
Definir las variables de reducción siguientes ( es la versión de reducción y variables críticas de , respectivamente),
,
donde:
como muestra Salzman.1
La primera forma de la ecuación de Van der Waals, mostrada más arriba, puede ser una refundición de la siguiente forma reducida:
Esta ecuación es invariante para todos los líquidos, es decir, se aplica la misma ecuación en forma reducida de estado, no importa a cual  y  puede ser que el fluido particular.
Esta invariancia también puede ser entendida en términos del principio de estados correspondientes. Si dos fluidos tienen la misma presión, volumen y temperatura reducida, se puede decir que sus estados son correspondientes. Los estados de dos fluidos puede ser correspondientes, incluso si su presión medida, volumen y temperatura sean muy diferentes. Si los estados de los dos fluidos son correspondientes, existen en el mismo régimen de la ecuación de forma reducida del Estado. Por lo tanto, van a responder a los cambios en más o menos de la misma manera, a pesar de sus características físicas mensurables puedan diferir significativamente.

Ecuación cúbica

La ecuación de Van der Waals es una ecuación cúbica de estado. Es decir, se puede escribir la ecuación en una forma cúbica del volumen. En la formulación de la ecuación cúbica reducida es la siguiente:
A la temperatura crítica, donde  se tiene, como era de esperar,
⇔ 
Para , hay 3 valores de . Para , existe un valor real para .

Aplicación a fluidos compresibles

La ecuación también se puede utilizar como una ecuación dependiente de P, V, T para fluidos compresibles, ya que, en este caso, los cambios en el volumen específico son pequeños, y se puede escribir de la siguiente manera:
donde:
p es la presión
V es el volumen específico
T es la temperatura
AB and C son parámetros.

Regla de área igual de Maxwell

Por debajo de la temperatura crítica (T’ < 1) como se muestra, oscila una isoterma de la ecuación de Van der Waals.
La regla de Maxwell elimina el comportamiento oscilante de la isoterma en la zona de cambio de estado definiéndolo como una cierta isobara en esa zona.
A lo largo de la porción roja de la isoterma  que es inestable, la ecuación de Van der Waals falla para describir las sustancias reales en esta región debido a que dicha ecuación siempre asume que el fluido es uniforme, mientras que entre a y c en la isoterma se vuelve más estable que un coexistencia de dos fases, una fase más densa que normalmente se llama líquida y una fase más dispersa que normalmente se llama gaseosa. Para solucionar este problema James Clerk Maxwell (1875) sustituyó a la isoterma de entre a y c, con una línea horizontal de modo que las áreas de las dos regiones establecidas sean iguales. La parte de la línea fija de la isoterma de ahora corresponde al equilibrio líquido-vapor. Las partes ad y ce se interpretan como estados metaestables de líquidos super-calientes y super-refrigerados por vapor, respectivamente.2
Maxwell justifica la regla, diciendo que el trabajo realizado en el sistema al pasar desde c a b debe ser igual a trabajo liberado al pasar de a a b. (El área en el diagrama PV se corresponde a un trabajo mecánico). Eso es porque el cambio en la función energía libre A (TV) es igual al trabajo realizado durante un proceso reversible de la función energía libre siendo una variable de estado debe tomar un valor único sin tener en cuenta de la ruta. En particular, el valor de A en el punto b debe calcular el mismo, independientemente de si la ruta procede de la izquierda o de la derecha, o se dirigen directamente a través de la isoterma horizontal o alrededor de la isoterma original de Van der Waals. El Argumento de Maxwell no es totalmente convincente, ya que requiere de un proceso reversible a través de una región de inestabilidad termodinámica. Sin embargo, los argumentos más sutiles sobre la base de las modernas teorías de equilibrio de fases parecen confirmar la construcción del área igualitaria de Maxwell y sigue siendo válida una modificación de la ecuación de Van der Waals.3
La regla de igual área de Maxwell se pueden derivar del supuesto de la igualdada del potencial químico, μ, en la coexistencia de las fase líquida y de vapor.4

Obtención de parámetros a partir del punto crítico

Los parámetros de la ecuación de Van der Waals,  y  son obtenidos en general a partir de las condiciones de punto crítico de un compuesto, es decir, el punto en que las fases en equilibrio líquida y gaseosa se hacen idénticas. En dicho punto se verifica
Donde el subíndice 'c' implica que la propiedad es la correspondiente al punto crítico. Aplicando estas condiciones a la ecuación de Van der Waals se obtiene un sistema de 3 ecuaciones, por lo que una de las propiedades deberá variar para tener un sistema de ecuaciones definido. Si bien posible usar como parámetro ajustable la constante [cita requerida], la misma suele dejarse fija, obteniendo el límite del gas ideal en el caso de que  y  sean cero. El parámetro que se permite variar es entonces el volumen crítico.[cita requerida]
La resolución del sistema entonces conlleva a las siguientes expresiones:
Es aquí donde se ve que fijando temperatura y presión en su punto crítico, la ecuación de Van der Waals no da buenas predicciones de los volúmenes molares de compuestos puros. Siendo el factor de compresibilidad:
Para la ecuación de Van der Waals se obtiene el valor fijo de:
Como esta ecuación está derivada para moléculas con simetría esférica, si se compara este valor con el de compuestos como los gases nobles, o incluso metanooxígeno y nitrógeno moleculares (que no son exactamente esféricos, pero se les aproximan), se encuentra que su valor experimental se haya alrededor de 0,29 aproximadamente.5 Es entonces que la ecuación de Van der Waals original, en general, predice valores de volúmenes molares mayores (o equivalentemente, densidades menores) a las experimentales.

La ecuación de van der Waals

La ecuación de estado de un gas ideal es PV=nRT
P es la presión del gas, V es el volumen, T es la temperatura, es el número de moles y R la constante de los gases
Describe aproximadamente la conducta de los gases reales a muy bajas presiones.
La ecuación de van der Waals tiene en cuenta el volumen finito de las moléculas y las fuerzas atractivas que una molécula ejerce sobre otra a distancias muy cercanas entre ellas.
(P+n2aV2)(Vnb)=nRT
Las constantes a y b son característicos de cada gas y se obtienen a partir de los datos de la presión, Pc, volumen Vc y la temperatura Tc crítica. El punto crítico es un punto de inflexión de la isoterma Tc en el diagrama P-V de modo que se cumple que
(PV)T=0nRT(Vnb)2+2n2aV3=0(2PV2)T=02nRT(Vnb)36n2aV4=0
De estas dos ecuaciones obtenemos, el volumen Vc y la temperatura Tc crítica.
Vc=3nbTc=8a27Rb
Sustituyendo Vc y Tc en la ecuación de van der Waals obtenemos la presión crítica, Pc,
Pc=a27b2
Escribimos la ecuación de van der Waals universal más simple en términos de nuevas variables: p=P/Pcv=V/Vc y t=T/Tc
(p+3v2)(3v1)=8t
En el punto crítico pc=1, tc=1 y vc=1.
Como vamos a dibujar las isotermas en un diagrama P-V, despejamos la presión p.
p=8t3v13v2
Creamos un script para representar las isotermas de temperaturas t=0.8,0.9,1.0,1.1 y 1.2. Señalamos con una marca el punto crítico tal como se aprecia en la figura
f=@(t,x) 8*t./(3*x-1)-3./(x.^2);
hold on
axis([0.5 4 0 2])
for t=0.8:0.1:1.2;
    v=linspace(0.5,4,100);
    p=f(t,v);
    plot(v,p)
end
plot(1,1,'ro')
xlabel('V_R')
ylabel('P_R')
grid on
hold off
Examinamos las isotermas por debajo de la temperatura crítica t<1 .="" p="">
En la figura vemos que para una determinada presión p, la recta p=cte corta a la curva en tres puntos A, B y C. Los volúmenes v1v2 y v3 que corresponden a estos tres puntos se obtienen resolviendo la ecuación cúbica
v38t+p3pv2+3pv1p=0
Para obtener las raíces podemos utilizar la función roots de MATLAB o la función raices_3 que nos calcula las tres raíces exactas de una ecuación de tercer grado.
function x = raices_3(p)
    Q=(p(2)*p(2)-3*p(3))/9;
    R=(2*p(2)^3-9*p(2)*p(3)+27*p(4))/54;
    x=zeros(3,1);
    if (R*R)<(Q^3)
        tetha=acos(R/sqrt(Q^3));
        x(1)=-2*sqrt(Q)*cos(tetha/3)-p(2)/3;
        x(2)=-2*sqrt(Q)*cos((tetha+2*pi)/3)-p(2)/3;
        x(3)=-2*sqrt(Q)*cos((tetha-2*pi)/3)-p(2)/3;
    else
        A=-sign(R)*nthroot(abs(R)+sqrt(R*R-Q^3),3);
        if A==0
            B=0;
        else
            B=Q/A;
        end
        x(1)=(A+B)-p(2)/3;
        x(2)=-(A+B)/2-p(2)/3+(sqrt(3)*(A-B)/2)*sqrt(-1); 
        x(3)=-(A+B)/2-p(2)/3-(sqrt(3)*(A-B)/2)*sqrt(-1);
    end
end
Ordenamos los valores de las raíces de menor a mayor mediante la función sort de MATLAB.
La isoterma de temperatura t, tiene un mínimo y un máximo local. Determinamos sus coordenadas (vm1pm1) y (vm2pm2) calculando la derivada primera e igualando a cero.
p=8t3v13v2(pv)t=024t(3v1)26v3=0v394tv2+32tv14t=0
Tenemos de nuevo, que calcular las raíces de una ecuación cúbica.
En la zona de transición entre la fase gaseosa y la fase líquida la presión no oscila, se mantiene constante y el potencial químico se mantiene constante. La regla de Maxwell elimina el comportamiento oscilante de la isoterma de la ecuación de van der Waals y la sustituye por un segmento horizontal de presión pr tal que el área comprendida entre el segmento horizontal AB y la isoterma es igual al área comprendida entre la isoterma y el segmento horizontal BC, tal como se ve en la figura más abajo.
En la figura podemos ver la isoterma de temperatura t=0.9,
  • los extremos locales identificados con el símbolo '*',
  • la presión pr para la cual las áreas sombreadas de color amarillo son iguales.
  • los puntos de corte de la isoterma con la recta pr =cte
El el área A1comprendida entre el segmento horizontal AB y la isoterma se calcula integrando
p(v2v1)v1v2p·dv=p(v2v1)v1v2(8t3v13v2)·dv=p(v2v1)83tln3v213v113(1v21v1)
El el área A2 comprendida entre el segmento horizontal BC y la isoterma se calcula integrando
v2v3p·dvp(v3v2)=83tln3v313v21+3(1v31v2)p(v3v2)
Creamos una función denominada igualArea que calcule la diferencia entre las dos áreas A1-A2.
function res=igualArea(p,t)
    pol=[1 -(p+8*t)/(3*p) 3/p -1/p];
    v = sort(raices_3(pol));
    area_1 = p*(v(2)-v(1))-8/3*t*log((3*v(2)-1)
/(3*v(1)-1))-3*(v(1)-v(2))/(v(1)*v(2));
    area_2 =8/3*t*log((3*v(3)-1)/(3*v(2)-1))+3*(v(2)-v(3))
/(v(2)*v(3))-p*(v(3)-v(2));
    res=area_1-area_2;
end
Creamos un script para realizar las siguientes tareas:.
  1. Establezca el valor de la temperatura t reducida, menor que la unidad
  2. Represente la isoterma de color rojo en el diagrama P-V.
  3. Calcule el mínimo pm1 y máximo pm2 local resolviendo la ecuación de tercer grado en v.
  4. La presión pr buscada estará en el intervalo ( pm1pm2 ), o bien en el intervalo ( 0, pm2 ) si pm1 fuese negativa (la presión es siempre positiva). El manejador (handle) de esta función igualArea se lo pasaremos a la función fzero de MATLAB para que calcule el valor pr de la presión que hace la que diferencia de áreas sea nula.
  5. Una vez obtenido el valor de la presión pr, se calcula las raíces de la ecuación cúbica para determinar los volúmenes v1v2 y v3. El primero corresponde al volumen de la fase líquida y el tercero al del gas, para la presión pr y la temperatura t
  6. Rellene de color amarillo las dos regiones, para mostrar visualmente que las áreas son iguales.
t=input('temperatura reducida (0.8-1): ');
%gráfica de la isoterma
f=@(x) 8*t./(3*x-1)-3./(x.^2);
Vr=linspace(0.5,4,100);
Pr=f(Vr);
hold on
axis([0 4 0 1.4])
plot(Vr,Pr,'r')
xlabel('V_R')
ylabel('P_R')

%máximo y mínimo
pol=[1 -9/(4*t) 3/(2*t) -1/(4*t)];
z=sort(raices_3(pol));
plot(z(1),f(z()),z(2),f(z(2)),'k*',z(3),f(z(3)),'k*')
if f(z(2))>0
    a=f(z(2))+0.01;
else
    a=0.01;
end
b=f(z(3))-0.01;
%calcula p para que las áreas sean iguales
f1=@(p)igualArea(p,t);
p=fzero(f1,[a b]);
%igualArea(p,t) %para comprobar que las areas son iguales
plot([0 4],[p p],'k--') 
pol=[1 -(p+8*t)/(3*p) 3/p -1/p];
v = sort(raices_3(pol));
plot(v(1),p,'bo',v(2),p,'bo',v(3),p,'bo')

%rellenar la areas
xx = [v(1) Vr(Vr >= v(1) & Vr <= v(2))  v(2)];
yy = [f(v(1)) Pr(Vr >= v(1) & Vr <= v(2)) f(v(2))];
fill(xx,yy,'y')

xx = [v(2) Vr(Vr >= v(2) & Vr <= v(3)) v(3)];
yy = [f(v(2)) Pr(Vr >= v(2) & Vr <= v(3)) f(v(3))];
fill(xx,yy,'y')
plot(Vr,Pr,'r')
%imprimir datos
text(1,1.35,sprintf('presión %1.3f',p))
text(1,1.25,sprintf('volumen gas %1.3f',v(3)))
text(1,1.15,sprintf('volumen líquido %1.3f',v(1)))

hold off
Ahora vamos a dibujar varias isotermas por encima y por debajo del punto crítico t=1. Sustituiremos la región ondulada entre vliquido y vgas a la presión pr por un segmento y uniremos los puntos que marcan la región de coexistencia de las dos fases, tal como se muestra en la figura.
Creamos un script para que realice las siguientes tareas.
  1. Dibuje las isotermas t≥1 (para t=1,1.05,1.1,1.15 y 1.2)
  2. Dibuje las isotermas para t<1 em="" nbsp="" para="">t
=0.8, 0.85, 0.9, 0.95).
  • Para cada isoterma se calcula el volumen vliquido=v(1) y vgas=v(3) a la presión prresolviendo la ecuación trascedente A1-A2=0, mediante la función fzero, guardando los datos en una matriz bidimiensional (volumen, presión).
  • Se sustituye la región ondulada por el segmento que une v(1) y v(3) a la presión pr,
  • Se añade a la matriz bidimensional las coordendas del punto crítico v=1 y p=1. Ordenamos la matriz de datos (v,p) utilizando la función MATLAB sortrows, y unimos los puntos (v,p) mediante líneas de color azul tal como se muestra en la figura.
  • Vr=linspace(0.1,4,100);
    hold on
    i=1;
    PV=zeros(9,2);
    axis([0.4 4 0.2 1.4])
    for t=0.8:0.05:1.2;
        f=@(x) 8*t./(3*x-1)-3./(x.^2);
        if t<1 -1="" -9="" 3="" calcula="" de="" f="" if="" local="" m="" nimo="" p="" pol="[1" t="" ximo="" y="" z="">0
                a=f(z(2))+0.01;
            else
                a=0.01;
            end
            b=f(z(3))-0.01;
            %calcula p para que las áreas sean iguales
            f1=@(p)igualArea(p,t);
            p=fzero(f1,[a b]);
            pol=[1 -(p+8*t)/(3*p) 3/p -1/p];
            v = sort(raices_3(pol));
            %guarda Vlíquido, Vgas y la presión constante p 
            PV(i,1)=v(1);PV(i,2)=p;
           PV(i+1,1)=v(3);PV(i+1,2)=p;
           i=i+2;
     %gráfica
            Pr=f(Vr);
            Vf = [Vr(Vr <= v(1)) v(1) v(3) Vr(Vr>= v(3))];
            Pf = [Pr(Vr <= v(1)) f(v(1)) f(v(3)) Pr(Vr>= v(3))];
            plot(Vf,Pf,'r')
        else %en el punto crítico y por encima
            Pr=f(Vr);
            plot(Vr,Pr,'r')
        end 
    end
    %une los puntos (p Vlíquido), (p Vgas) cambio de fase
    PV(i,1)=1; PV(i,2)=1; %punto crítico
    PVs=sortrows(PV,1);
    plot(PVs(:,1),PVs(:,2),'b')
    plot(1,1,'bo') %punto crítico
    xlabel('V_R')
    ylabel('P_R')
    hold off

    No hay comentarios:

    Publicar un comentario