APLICACIÓN DEL MODELO DE 3 CONSTANTES PARA ALGUNAS FUNCIONES CRECIENTES Y

DECRECIENTES

Luis Andy Espino Martinez, Lima-Perú, [email protected]

Sea el siguiente modelo matemático empleado en el Trabajo Investigativo denominado: “MODELO

MATEMÁTICO EMPÍRICO QUE RELACIONA LA VISCOSIDAD DE UN FLUIDO CON LA
TEMPERATURA”

Y = α. (EXP (β.Xγ)) … (1)

DONDE ‘Y’ ES LA VARIABLE DEPENDIENTE Y ‘X’ ES LA VARIABLE INDEPENDIENTE.

Podemos usar el mismo procedimiento de desarrollo y determinación de las constantes

mostradas, en el mencionado trabajo investigativo, para funciones cuadráticas, lineales, etc.
Procedamos a verificar.

ECUACIONES DE MOVIMIENTO:

1. MOVIMIENTO RECTILÍNEO ACELERADO (MODELO CRECIENTE CUADRÁTICO):

X (t) = X0 + V0t + at2/2 … (2)

POR EJEMPLO: X = 3 + 4t + 5t2/2 … (3)

TOMEMOS 9 PUNTOS INCREMENTADOS DESDE 0,1 A 0,9. LOS PUNTOS DE LA SOLUCIÓN SERÍAN:

T X
0,1 3,425
0,2 3,9
0,3 4,425
0,4 5
0,5 5,625
0,6 6,3
0,7 7,025
0,8 7,8
0,9 8,625

EMPLEANDO EL MODELO DE 3 CONSTANTES, usando un programa de Búsqueda de Intervalos y

cálculo de las raíces por el Método de Newton-Raphson, puedo simplificar el trabajo de calcular:
HALLO LAS CONSTANTES EN UN INTERVALO DE ‘C’ EVALUADO DE -10 A 10:

HALLO LOS PUNTOS PARA EL TIEMPO QUE VA DE 0,1 A 0,9 INCREMENTADO EN 0,1:
TABULAMOS LOS RESULTADOS:

t X
0,1 3,425
0,2 3,9
0,3 4,409
0,4 4,961
0,5 5,565
0,6 6,226
0,7 6,952
0,8 7,749
0,9 8,625

INCLUSO PARA t=0s, X=2,941; LO CUAL CONCUERDA CON: X (0) = 3 + 4(0) + 5(0)2/2 = 3

*PRIMERA RESTRICCIÓN: NO SE PUEDE USAR EL VALOR DE VARIABLE INDEPENDIENTE IGUAL A

‘CERO’ PARA HALLAR LAS CONSTANTES, PORQUE AL EVALUAR LA VARIABLE INDEPENDIENTE EN
‘CERO’, PARA UNA CONSTANTE ‘GAMMA’ IGUAL A ‘CERO’ NO TENDRÍA SOLUCIÓN:

[Ln(Y1 / Y2)/ Ln(Y2 / Y3) + 1].X2γ – [Ln(Y1 / Y2)/ Ln(Y2 / Y3)].X3γ - X1γ = 0 …(4)
*SEGUNDA RESTRICCIÓN: LOS PUNTOS QUE SE EMPLEAN PARA HALLAR LAS CONSTANTES DEL
MODELO DEBEN SER: EL PUNTO DE PARTIDA, UN PUNTO MEDIO (COMO EL SEGUNDO PUNTO POR
EJEMPLO) Y EL ÚLTIMO PUNTO DE UN GRUPO DE DATOS EN EL QUE SE DESEA USAR EL MODELO.

*TERCERA RESTRICCIÓN: EL MODELO SOLO PUEDE SER CRECIENTE O DECRECIENTE, ENTONCES

PARA X1<X2<X3, CUMPLE QUE Y1<Y2<Y3, O SINO Y1>Y2>Y3.

AHORA COMPARAMOS LOS RESULTADOS OBTENIDOS CON LA SOLUCIÓN ANALÍTICA Y EL

MODELO DE 3 CONSTANTES:

SOLUCIÓN ANALÍTICA MODELO DE 3 CONSTANTES

2. MOVIMIENTO ARMÓNICO SOBREAMORTIGUADO (MODELO DECRECIENTE

EXPONENCIAL):

X (t) = A0.e-γt. Cos(wt) … (5)

POR EJEMPLO PARA EL CASO EN QUE:

X (t) = (0,5).e-2t. Cos(√3t) … (6)

TOMEMOS 6 PUNTOS INCREMENTADOS DESDE 0,1 A 0,6. LOS PUNTOS DE LA SOLUCIÓN SERÍAN:

t X
0,1 0,409363506
0,2 0,3351538973
0,3 0,2743945336
0,4 0,2246480574
0,5 0,1839187093
0,6 0,1505723344

EMPLEANDO EL MODELO DE 3 CONSTANTES:

HALLO LAS CONSTANTES EN UN INTERVALO DE ‘C’ EVALUADO DE -10 A 10:

HALLO LOS PUNTOS PARA EL TIEMPO QUE VA DE 0,1 A 0,6 INCREMENTADO EN 0,1:
TABULAMOS LOS RESULTADOS:

t X
0,1 0,409363506
0,2 0,3351538973
0,3 0,274393576778
0,4 0,224646701718
0,5 0,18391769029127
0,6 0,15057221928202

INCLUSO PARA t=0s, X=0,499987846; LO CUAL CONCUERDA CON: X(0)=(0,5).e-2(0). Cos(√3(0))= 0,5

AHORA COMPARAMOS LOS RESULTADOS OBTENIDOS CON LA SOLUCIÓN ANALÍTICA Y EL MODELO

DE 3 CONSTANTES:
 SOLUCIÓN ANALÍTICA MODELO DE 3 CONSTANTES

3. MOVIMIENTO RECTILÍNEO UNIFORME (MODELO CRECIENTE LINEAL):

X (t) = X0 + Vt … (2)

POR EJEMPLO: X = 2 + 4t … (3)

TOMEMOS 9 PUNTOS INCREMENTADOS DESDE 0,1 A 0,9. LOS PUNTOS DE LA SOLUCIÓN SERÍAN:

t X
0,1 2,4
0,2 2,8
0,3 3,2
0,4 3,6
0,5 4,0
0,6 4,4
0,7 4,8
0,8 5,2

EMPLEANDO EL MODELO DE 3 CONSTANTES:

HALLO LAS CONSTANTES EN UN INTERVALO DE ‘C’ EVALUADO DE -10 A 10:

HALLO LOS PUNTOS PARA EL TIEMPO QUE VA DE 0,1 A 0,8 INCREMENTADO EN 0,1:
TABULAMOS LOS RESULTADOS:

t X
0,1 2,40
0,2 2,80
0,3 3,18
0,4 3,56
0,5 3,95
0,6 4,35
0,7 4,77
0,8 5,20

INCLUSO PARA t=0s, X=1,83; LO CUAL NO ESTÁ TAN LEJOS DE: X(0)=2+4(0)= 2

AHORA COMPARAMOS LOS RESULTADOS OBTENIDOS CON LA SOLUCIÓN ANALÍTICA Y EL MODELO

DE 3 CONSTANTES:

SOLUCIÓN ANALÍTICA MODELO DE 3 CONSTANTES

*CUARTA RESTRICCIÓN: PARA EL CASO DE UNA FUNCIÓN LINEAL, NO SE PUEDE APLICAR ESTE
MODELO A FUNCIONES LINEALES DEL TIPO: Y=mX, YA QUE NECESARIAMENTE CUMPLEN PARA EL
PUNTO (0,0) Y NO ES POSIBLE QUE ESTO CUMPLA, YA QUE NECESARIAMENTE PARA EL MODELO
DE TRES CONSTANTES, PARA QUE LA VARIABLE DEPENDIENTE SEA CERO (Y=0), EL RESULTADO DEL
EXPONENTE AL QUE DEBE QUEDAR ELEVADO EL EXPONENCIAL, DEBE SER -∞, DE ESTE MODO,
PARA X=-∞, Y=0, POR LO CUAL NO HAY SOLUCIÓN FINITA Y NO SE PUEDE AJUSTAR A DICHA
FUNCIÓN ESTE MODELO TAMPOCO.
DETERMINACIÓN DE LA RELACIÓN PRESIÓN-TEMPERATURA Y VOLUMEN-TEMPERATURA DE UN
GAS

TOMANDO DE REFERENCIA EL ARTÍCULO DE LA WEB RESPECTO DEL AIRE:

https://espanol.libretexts.org/Quimica/Libro%3A_Quimica_General_(OpenSTAX)/09%3A_Gases/
9.2%3A_La_presi%C3%B3n%2C_el_volumen%2C_la_cantidad_y_la_temperatura_relacionados
%3A_la_ley_del_gas_ideal

ENCONTRAMOS LA SIGUIENTE TABLA DE DATOS EXPERIMENTALES PARA LA RELACIÓN P-T:

PROBEMOS POR LA ECUACIÓN DEL GAS IDEAL PѴ=RT, SIENDO R = 8,314472

KPa.L/mol°K, HALLANDO Ѵ:

Ѵ1 = 39,95 L/mol

Ѵ2 = 39,96 L/mol

Ѵ3 = 40,03 L/mol

Ѵ4 = 40,02 L/mol

Ѵ5 = 40,02 L/mol

Ѵ6 = 39,97 L/mol

HALLAMOS EL PROMEDIO: ѴPROMEDIO = 39,99 L/mol

CON NUESTRO MODELO HALLAMOS LA RELACIÓN PRESIÓN-TEMPERATURA:

EL CÁLCULO DE LAS CONSTANTES:
Y LA COMPROBACIÓN:

ENTONCES:

P = (3,491022217967968X10-23). (EXP (50,41126744536285.T0,01793565)) …(4)

SI ѵ = 39,99 L/mol …(4.1)

MULTIPLICO (4) POR (4.1) Y OBTENGO:

Pѵ = (1,396059785X10-21). (EXP (50,41126744536285.T0,01793565)) …(4.2)

ENCONTRAMOS TAMBIÉN LA SIGUIENTE TABLA DE DATOS EXPERIMENTALES PARA LA RELACIÓN
V-T (PARA EL VOLUMEN MOLAR):
PROBEMOS POR LA ECUACIÓN DEL GAS IDEAL PѴ=RT, SIENDO R =
8,314472 KPa.L/mol°K, HALLANDO P:

P1 = 102,04 KPa

P2 = 98,98 KPa

P3 = 101,62 KPa

P4 = 101,62 KPa

HALLAMOS EL PROMEDIO: PPROMEDIO = 101,07 KPa

CON NUESTRO MODELO HALLAMOS LA RELACIÓN VOLUMEN-TEMPERATURA:
EL CÁLCULO DE LAS CONSTANTES:

Y LA COMPROBACIÓN:
ENTONCES:

Ѵ = (9216050,255125787). (EXP (-19,71237238756384.T-0,07511131)) …(5)

SI P = 101,07 KPa …(5.1)

MULTIPLICO (5) POR (5.1) Y OBTENGO:

PѴ = (9,314661993X108). (EXP (-19,71237238756384.T-0,07511131)) …(5.2)

MULTIPLICAMOS (4.2) POR (5.2) Y ELEVAMOS AL EXPONENTE ½:

Pѵ = (1,140343151X10-6).(EXP(25,20563372T0,01793565-9,856186194T-0,07511131))
…(6)
CALCULEMOS PARA LOS DATOS QUE TENEMOS A PARTIR DE (6):

T=111°K, Ѵ=9 L/mol, P= 102,1 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=0,99)

T=173°K, Ѵ=39,9 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0,25)

T=220°K, Ѵ=18 L/mol, P=101,3 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=0,20)

T=223°K, Ѵ=40,7 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=1,75)

T=250°K, Ѵ=21 L/mol, P=98,6 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=2,47)

T=270°K, Ѵ=22 L/mol, P=101,5 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=0,40)

T=273°K, Ѵ=39,8 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0,50)

T=323°K, Ѵ=39,7 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0,75)

T=373°K, Ѵ=39,6 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=1,00)

T=423°K, Ѵ=39,4 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=1,50)

%EPROMEDIO=0,98

PODEMOS NOTAR QUE, LOS RESULTADOS OBTENIDOS PARA LOS DATOS DE LAS TABLAS, SON MUY
CERCANOS AL VALOR PROMEDIO (MOSTRADO ENTRE PARÉNTESIS). POR LO TANTO, LA ECUACIÓN TIENE
COHERENCIA Y PODEMOS COMPARARLO CON LA ECUACIÓN DEL GAS IDEAL:

T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=1,38)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=0,49)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=2,08)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (PPROMEDIO = 101,1 KPa) (%E=0,89)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (ѴPROMEDIO = 40,0 L/mol) (%E=0)

%EPROMEDIO=0,48
AHORA PROBEMOS APLICANDO LA ECUACIÓN DE VAN DER WAALS DE LA PÁGINA WEB:

https://edico3004.blogspot.com/2017/11/ejercicios-con-la-ecuacion-de-estado.html

…(7)

SI 1atm = 0,101325 x 106 = 101,325 KPa, ENTONCES: a = 134,76 KPa(L/mol)2,

ADEMÁS SABEMOS QUE R = 8,314472 KPa.L/mol°K. HALLAMOS P COMPARANDO
CON EL MODELO DEL GAS IDEAL:

T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (PVDV = 101,3 KPa) (%E=1,17)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (PVDV = 35,9 KPa) (%E=0,28)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (PVDV = 101,4 KPa) (%E=0,20)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (PVDV = 46,3 KPa) (%E=0,22)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (PVDV = 98,8 KPa) (%E=0,20)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (PVDV = 101,9 KPa) (%E=0,10)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (PVDV = 56,7 KPa) (%E=0)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (PVDV = 67,1 KPa) (%E=0)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (PVDV = 77,5 KPa) (%E=0)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (PVDV = 87,9 KPa) (%E=0,11)

%EPROMEDIO=0,23
 COMPAREMOS LOS RESULTADOS DE LA ECUACIÓN DEL GAS IDEAL CON LOS
RESULTADOS DE NUESTRO MODELO HALLANDO LA PRESIÓN PNM:

T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 102,1 KPa) (%E=0,39)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (PNM = 35,9 KPa) (%E=0,28)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 101,3 KPa) (%E=0,30)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (PNM = 46,2 KPa) (%E=0,43)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 98,6 KPa) (%E=0,40)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 101,5 KPa) (%E=0,49)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (PNM = 56,4 KPa) (%E=0,53)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (PNM = 66,6 KPa) (%E=0,75)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (PNM = 76,7 KPa) (%E=1,03)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (PNM = 86,7 KPa) (%E=1,48)

%EPROMEDIO=0,61
COMO PODEMOS NOTAR, HAY COHERENCIA EN LOS RESULTADOS, POR LO CUAL,
PODEMOS HACER UNA MEJOR APROXIMACIÓN A PARTIR DE LAS ECUACIONES (4.2) Y
(5.2). PRIMERO PARA EL CASO DE (4.2) A VOLUMEN CONSTANTE:

Pѵ = (1,396059785X10-21). (EXP (50,41126744536285.T0,01793565)) …(4.2)

Y PARA (5.2) A PRESIÓN CONSTANTE:

PѴ = (9,314661993X108). (EXP (-19,71237238756384.T-0,07511131)) …(5.2)

PODEMOS DERIVAR (4.2) RESPECTO A LA TEMPERATURA A VOLUMEN CONSTANTE:

ѵdP = (1,262259808X10-21 T0,01793565-1).(EXP (50,41126744536285.T0,01793565))dT

…(8)
DE (4.2) DESPEJO EL EXPONENCIAL:

Pѵ/(1,396059785X10-21) = EXP (50,41126744536285.T0,01793565) …(9)

REEMPLAZO (9) EN (8) Y OBTENGO:
 ѵdP = (1,262259808X10-21 T0,01793565-1).(Pѵ/1,396059785X10-21)dT …(10)
AHORA, PODEMOS DERIVAR (5.2) RESPECTO A LA TEMPERATURA A PRESIÓN CONSTANTE:

PdѴ = (1,379149452X109 T-0,07511131-1). (EXP (-19,71237238756384.T-0,07511131))dT …

(11)

DE (5.2) DESPEJO EL EXPONENCIAL:

PѴ/(9,314661993X108) = (EXP(-19,71237238756384.T-0,07511131)) …(12)

REEMPLAZO (12) EN (11) Y OBTENGO:

PdѴ = (1,379149452X109 T-0,07511131-1). (PѴ/9,314661993X108)dT …(13)

SUMAMOS (13) Y (10):

PdѴ + ѵdP = d(PV) = (PV/T)(1,480622112T-0,07511131+0,904158848T0,01793565)dT …(14)

SEPARO A UN LADO TODO PѴ Y A OTRO LADO T:

d(PѴ)/(PѴ) =(dT/T)(1,480622112T-0,07511131+0,904158848T0,01793565) …(15)

INTEGRANDO TENEMOS PARA UNA CONSTANTE K:

ln(PѴ) = (0,904158848/0,01793565)T0,01793565- (1,480622112/0,07511131)T-0,07511131 + K …

(16)

RESOLVIENDO LAS DIVISIONES DE LAS CONSTANTES:

ln(PѴ) = (50,411267392)T0,01793565 - (19,712372371)T-0,07511131 + K …(17)

DE LA FORMA EXPONENCIAL:

EXP(ln(PѴ)) = PѴ = (EXP(K)).(EXP((50,411267392)T0,01793565-(19,712372371)T-0,07511131)

…(18)

DESIGNAMOS LA CONSTANTE Z= EXP(K), ENTONCES:

PѴ = Z (EXP((50,411267392)T0,01793565-(19,712372371)T-0,07511131) …(19)

HALLEMOS Z PARA LOS SIGUIENTES DATOS OBTENIDOS PARA EL GAS IDEAL PARA EL CASO
DEL AIRE:
 T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (Z1 = 1,419994242X10-15)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (Z2 = 0,9084164139X10-15)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (Z3 = 0,7150796411X10-15)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (Z4 = 0,7064828133X10-15)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (Z5 = 0,6311428216X10-15)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (Z6 = 0,5852117307X10-15)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (Z7 = 0,5787236687X10-15)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (Z8 = 0,4919896656X10-15)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (Z9 = 0,4286725946X10-15)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (Z10 = 0,3808126896X10-15)

SE OBSERVA CLARAMENTE QUE Z ES INVERSAMENTE PROPORCIONAL A T, ENTONCES

EXISTE UN θ = ZT, QUE LO PODEMOS HALLAR COMO UN PROMEDIO DE:

θ1 = Z1T1 = 1,576193609X10-13

θ2 = Z2T2 = 1,571560396X10-13

θ3 = Z3T3 = 1,57317521X10-13

θ4 = Z4T4 = 1,575456674X10-13

θ5 = Z5T5 = 1,577857054X10-13

θ6 = Z6T6 = 1,580071673X10-13

θ7 = Z7T7 = 1,579915616X10-13

θ8 = Z8T8 = 1,58912662X10-13

θ9 = Z9T9 = 1,598948778X10-13

θ10 = Z10T10 = 1,610837677X10-13

PROMEDIANDO LOS DATOS 1, 2, 3, 4, 5,6 Y 7 HALLAMOS θ = ZT = 1,576318605X10-13, POR

LO CUAL Z = θ/T, LO QUE PODEMOS REEMPLAZAR EN (19):
 PѴ = (1,576318605X10-13/T) (EXP((50,411267392)T0,01793565-(19,712372371)T-0,07511131) …
(20)

PROBEMOS PARA EL CASO DEL GAS IDEAL:

T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 102,5 KPa) (%E=0)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (PNM = 36,1 KPa) (%E=0,28)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 101,8 KPa) (%E=0,20)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (PNM = 46,4 KPa) (%E=0)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 98,9 KPa) (%E=0,10)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (PNM = 101,8 KPa) (%E=0,20)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (PNM = 56,6 KPa) (%E=0,18)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (PNM = 66,6 KPa) (%E=0,75)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (PNM = 76,4 KPa) (%E=1,42)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (PNM = 86,1 KPa) (%E=2,16)

%EPROMEDIO= 0,53
ESTOS RESULTADOS TIENEN COHERENCIA, POR LO CUAL:

RT2 ≈ (1,576318605X10-13/T).EXP ((50,411267392)T0,01793565-(19,712372371)T-0,07511131 )…(21)

AHORA PROBEMOS CON DATOS EXPERIMENTALES DE VAPOR DE ÁCIDO ACÉTICO DE LA

PÁGINA WEB:
https://www.modeladoeningenieria.edu.ar/images/IntegracionIII/TP/02_Problemas_presi
%C3%B3n_de_vapor.pdf
 SI 1mmHg = 0,133322 KPa, Y °K = °C +273, CAMBIAMOS DE UNIDADES:

T (°K) P (KPa)
290,5 1,33322
302,9 2,66644
316 5,33288
324,7 7,99932
336 13,33220
353 26,66440
372 53,32880
391 101,32472
416,5 202,64944

CON NUESTRO MODELO HALLAMOS LA RELACIÓN PRESIÓN-TEMPERATURA:

HALLAMOS LAS CONSTANTES:

COMPROBAMOS CON LOS DATOS DE LA TABLA (APROXIMADOS A 2 CIFRAS DECIMALES):

ENTONCES:

P = (1,362556601X106). (EXP (-1,678961605X104.T-1,252053)) …(22)

SE PUEDE APROXIMAR EL VOLUMEN APLICANDO LA ECUACIÓN DEL GAS IDEAL.
PROBEMOS POR LA ECUACIÓN DEL GAS IDEAL PѴ=RT, SIENDO R = 8,314472 KPa.L/mol°K,
HALLANDO Ѵ:

Ѵ1 = 1816,06 L/mol

Ѵ2 = 943,24 L/mol

Ѵ3 = 492,94 L/mol

Ѵ4 = 337,46 L/mol

Ѵ5 = 209,58 L/mol

Ѵ6 = 110,09 L/mol

Ѵ7 = 58,00 L/mol

Ѵ8 = 32,09 L/mol

Ѵ9 = 17,09 L/mol

SI SEGÚN LA PUBLICACIÓN:
http://erecursos.uacj.mx/bitstream/handle/20.500.11961/3856/Formulario%20de
%20Fisicoqu%C3%ADmica-M%C3%B3nica%20Galicia-14agosto.pdf?
sequence=1&isAllowed=y

LAS CONSTANTES DE LA ECUACIÓN DE VAN DER WAALS PARA EL ÁCIDO ACÉTICO SON:

a = 17,82 bar.L2/mol2 = 17,82x102 KPa.L2/mol2, b = 0,1068 L/mol y R = 8,314472

KPa.L/mol°K, HALLAMOS EL VOLUMEN MOLAR RESOLVIENDO DE FORMA NUMÉRICA LA
ECUACIÓN DE VDW Y NOTAMOS QUE LOS RESULTADOS OBTENIDOS POR LA ECUACIÓN
DEL GAS IDEAL SON SIMILARES:
 T (°K) Ѵ (L/mol)
290,5 1811,039
302,9 943,899
316 492,102
324,7 336,938
336 209,010
353 109,570
372 57,525
391 31,637
416,5 16,671

HALLAMOS CON NUESTRO MODELO LA RELACIÓN VOLUMEN-TEMPERATURA:

CALCULAMOS LAS CONSTANTES:
VERIFICAMOS LOS RESULTADOS:

ENTONCES:
 Ѵ = (6,240296208x10-3). (EXP (1,940015296x104.T-1,294353)) …(23)
MULTIPLICAMOS (22) POR (23):

PѴ = (8,502756790x103).(EXP(1,940015296x104.T-1,294353-1,678961605X104.T-1,252053)) …
(24)

COMPARAMOS CON LAS OTRAS ECUACIONES OBTENIDAS PARA EL CASO DEL AIRE:

Pѵ = (1,396059785X10-21). (EXP (50,41126744536285.T0,01793565)) …(4.2)

PѴ = (9,314661993X108). (EXP (-19,71237238756384.T-0,07511131)) …(5.2)

Pѵ = (1,140343151X10-6).(EXP(25,20563372T0,01793565-9,856186194T-0,07511131))
…(6)

PѴ = (1,576318605X10-13/T) (EXP((50,411267392)T0,01793565-(19,712372371)T-0,07511131) …
(20)

SI PUDIÉRAMOS EXPRESARLO COMO PѴ = RT(F(T)), DONDE F(T) ES UNA FUNCIÓN

EN LA TEMPERATURA Y R = 8,314472; ENTONCES:

LA ECUACIÓN (24) QUEDARÍA COMO:

PѴ = (RT).[(1022,64543/T)EXP(1,940015296x104.T-1,294353-1,678961605X104.T-1,252053)] …
(25)

LA ECUACIÓN (4.2) QUEDARÍA COMO:

Pѵ = RT[(1,679072087X10-22/T)EXP(50,41126744536285.T0,01793565)] …(26)
LA ECUACIÓN (5.2) QUEDARÍA COMO:

PѴ = RT[(1,12029507X108/T)EXP(-19,71237238756384.T-0,07511131)] …(27)
LA ECUACIÓN (20) QUEDARÍA COMO:

PѴ = RT[(1,89587337X10-14/T2)EXP(50,411267392T0,01793565-19,712372371T-0,07511131)] …(28)

COMPARAMOS PRIMERO (26) Y (27), APLICAMOS EL LOGARITMO NATURAL A LA PARTE

ENTRE CORCHETES Y TENEMOS:
ln(1,679072087X10-22)-lnT+50,41126744536285.T0,01793565=

ln(1,12029507X108)-lnT-19,71237238756384.T-0,07511131
ENTONCES:

50,41126744536285T0,01793565+19,71237238756384T-0,07511131=68,67290358 …
(29)

PROBANDO VALORES PARA T, LA ECUACIÓN (29) ES APROXIMADAMENTE

CORRECTA.
COMPARANDO (25) Y (28), APLICAMOS EL LOGARITMO NATURAL A LA PARTE ENTRE
CORCHETES Y TENEMOS:

ln(1022,64543)-lnT+1,940015296x104.T-1,294353-1,678961605X104.T-1,252053

= ln(1,89587337X10-14)-2lnT+50,411267392T0,01793565-19,712372371T-0,07511131 …(30)

ENTONCES:

lnT+38.52665983+1,940015296x104.T-1,294353-1,678961605X104.T-1,252053

=50,411267392T0,01793565-19,712372371T-0,07511131 …(31)

SUMANDO (29) MÁS (31) TENEMOS:

lnT+38.52665983+50,41126744536285T0,01793565+19,71237238756384T-0,07511131

+1,940015296x104.T-1,294353-1,678961605X104.T-1,252053

= 68,67290358+50,411267392T0,01793565-19,712372371T-0,07511131 …(32)

ENTONCES:

lnT=30,14624375-1,940015296x104.T-1,294353+1,678961605X104.T-1,252053

-39,42474474T-0,07511131 …(33)

PROBANDO VALORES PARA T, LA ECUACIÓN (33) TAMBIÉN ES

APROXIMADAMENTE CORRECTA.
POR LO CUAL, PODEMOS DEDUCIR ENTONCES QUE LA ECUACIÓN DE LOS GASES SE PUEDE
EXPRESAR SOLO COMO UNA FUNCIÓN EN T DEL PRODUCTO PV QUE REPRESENTA AL
TRABAJO TERMODINÁMICO W (W=PV):

Pѵ = RT (1 + f(T))
LO CUAL CONCORDARÍA CON NUESTRA EXPRESIÓN GENERALIZADA APROXIMADA (DEL
MODELO DE LA CONSERVACIÓN DE LA ENERGÍA-TIEMPO):

Pѵ ≈ RT (1 – β(T,t).(1/T))

LO EXPRESAMOS DE FORMA DETALLADA:

PV ≈ η(T,t).RT (1 – (nђ/(η(T,t).Rt))(1/T)) = η(T,t).RT (1 – β(T,t).(1/T))

PARA EL CASO DE UN SISTEMA CERRADO EN EL QUE EL NÚMERO DE MOLES η SE DEBE

MANTENER CONSTANTE Y EL TIEMPO t CONTADO DESDE UN TIEMPO t=0 DESDE EL INICIO
DEL TIEMPO DE VIDA DEL UNIVERSO HASTA EL TIEMPO ACTUAL, QUE ES UNA CANTIDAD
MUCHO MAYOR QUE CERO Y POR LO CUAL SE PUEDE ASUMIR COMO CONSTANTE,
ENTONCES:

PV ≈ ηRT (1 – (nђ/(ηRt))(1/T))

LO CUAL LO PODEMOS EXPRESAR EN FUNCIÓN AL VOLUMEN MOLAR ѵ = V/η:

Pѵ ≈ RT (1 – (nђ/(ηRt))(1/T))

SI DEL MODELO DE CONSERVACIÓN DE LA ENERGÍA-TIEMPO:

n= (N°partículasi X N°microestadoj)
SI CADA MICROESTADO SE PUEDE REPRESENTAR POR LA DISTRIBUCIÓN DE BOLTZMANN:

ni = n0EXP(-Ei/KT)
SI LA SUMA DE CADA UNO DE ESTOS MICROESTADOS ES IGUAL A n Y LA SUMATORIA DE
LOS ELEMENTOS EXPONENCIALES EN FUNCIÓN DE LA TEMPERATURA LO EXPRESAMOS
POR MEDIO DE NUESTRO MODELO DE 3 CONSTANTES, ENTONCES:

n = ∑ni = n0 ∑EXP(-Ei/KT) = n0(EXP(βTγ))

 ENTONCES EL VALOR DE n SE PUEDE REPRESENTAR POR EL MODELO ESTADÍSTICO DE
BOLTZMANN, AL CUAL LO PODEMOS GENERALIZAR COMO:

n = n0EXP(βTγ), SIENDO n0 EL ESTADO MÍNIMO CONSTANTE.

POR LO CUAL PUEDO REPRESENTAR LA APROXIMACIÓN DE LA ECUACIÓN DE ESTADO
GENERALIZADA COMO:

Pѵ ≈ RT (1 – (n0ђ/(ηRt))(EXP(βTγ))(1/T))

SI COMO YA MENCIONAMOS n0, η y t SE PUEDEN ASUMIR PARA UN SISTEMA

CERRADO COMO MUY APROXIMADAMENTE CONSTANTES, Y ADEMÁS SABEMOS
QUE ђ y R SON CONSTANTES UNIVERSALES, ENTONCES PODEMOS DEFINIR UNA
CONSTANTE α:

α = n0ђ/(ηRt) ≈ CONSTANTE

ENTONCES PODEMOS EXPRESAR LA ECUACIÓN DE ESTADO GENERALIZADA DE LA

SIGUIENTE FORMA:

Pѵ ≈ RT (1 – α(EXP(βTγ))/T)

LO CUAL PODEMOS RESOLVERLO APLICANDO NUESTRO MODELO DE 3 CONSTANTES.

EXPRESÉMOSLO COMO UNA APROXIMADA IGUALDAD Y EN FUNCIÓN A UNA FUNCIÓN EN
T QUE DENOMINAREMOS G(T). SI DESPEJAMOS TENEMOS:

G(T) = [(Pѵ/RT) – 1]T = -α(EXP(βTγ))

PROBEMOS PARA EL CASO DEL AIRE TABULANDO PARA CADA T EL VALOR DE G(T):
T (°K) G(T)
111 -0,04888
173 0,19199
220 -0,04616
223 0,22523
250 0,04594
270 -0,10913
273 -0,22261
323 -0,18936
373 -0,1561
423 0,35821

PROBEMOS PARA EL CASO DEL VAPOR DE ÁCIDO ACÉTICO TABULANDO PARA CADA T EL
VALOR DE G(T):
 T (°K) G(T)
290,5 0,000684
302,9 -0,000333
316 -0,000355
324,7 -0,003494
336 0,00467
353 -0,001108
372 0,018812
391 0,048138
416.5 0,037394

PARA EL CASO DEL AIRE HALLAMOS UNA APROXIMACIÓN SIMPLE DE FORMA

COMPUTACIONAL USANDO LOS DATOS SOMBREADOS EN VERDE:

α= 5,081105352x10-3, β= 0,453078018, γ=0,3415953

PROBEMOS LA ECUACIÓN PARA EL CASO DEL AIRE (CON LOS DATOS DEL MODELO DEL
GAS IDEAL):

Pѵ = RT (1 – α(EXP(βTγ))/T)

T=111°K, Ѵ=9 L/mol, P= 102,5 KPa (DE LA GRÁFICA Ѵ-T) (PCALCULADO = 102,5 KPa) (%E=0)

T=173°K, Ѵ=40,0 L/mol, P=36 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 39,9 L/mol) (%E=0,25)

T=220°K, Ѵ=18 L/mol, P=101,6 KPa (DE LA GRÁFICA Ѵ-T) (PCALCULADO = 101,6 KPa) (%E=0)

T=223°K, Ѵ=40,0 L/mol, P=46,4 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 39,9 L/mol) (%E=0,25)

T=250°K, Ѵ=21 L/mol, P=99,0 KPa (DE LA GRÁFICA Ѵ-T) (PCALCULADO = 98,9 KPa) (%E=0,10)

T=270°K, Ѵ=22 L/mol, P=102,0 KPa (DE LA GRÁFICA Ѵ-T) (PCALCULADO = 102,0 KPa) (%E=0)

T=273°K, Ѵ=40,0 L/mol, P=56,7 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 40,0 L/mol) (%E=0)

T=323°K, Ѵ=40,0 L/mol, P=67,1 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 40,0 L/mol) (%E=0)

T=373°K, Ѵ=40,0 L/mol, P=77,5 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 40,0 L/mol) (%E=0)

T=423°K, Ѵ=40,0 L/mol, P=88 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 39,9 L/mol) (%E=0,25)

%EPROMEDIO=0,09

PROBEMOS LA ECUACIÓN PARA EL CASO DEL ÁCIDO ACÉTICO (CON LOS DATOS DEL
MODELO DEL GAS IDEAL):
 Pѵ = RT (1 – α(EXP(βTγ))/T)

T=290,5°K, Ѵ=1816,06 L/mol, P= 1,33 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 1815,3 KPa) (%E=0,04)

T=302,9°K, Ѵ=943,24 L/mol, P=2,67 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 942,9 L/mol) (%E=0,04)

T=316°K, Ѵ=492,94 L/mol, P=5,33 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 492,7 KPa) (%E=0,05)

T=324,7°K, Ѵ=337,46 L/mol, P=8,00 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 337,3 L/mol) (%E=0,05)

T=336°K, Ѵ=209,58 L/mol, P=13,33 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 209,5 KPa) (%E=0,04)

T=353°K, Ѵ=110,09 L/mol, P=26,66 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 110,0 KPa) (%E=0,08)

T=372°K, Ѵ=58,00 L/mol, P=53,33 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 58,0 L/mol) (%E=0)

T=391°K, Ѵ=32,09 L/mol, P=101,32 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 32,1 L/mol) (%E=0,03)

T=416,5°K, Ѵ=17,09 L/mol, P=202,65 KPa (DE LA GRÁFICA P-T) (ѴCALCULADO = 17,1 L/mol) (%E=0,06)

%EPROMEDIO=0,04

COMPARAMOS CON LOS RESULTADOS OBTENIDOS DE LA EC. DE VDW:

T (°K) ѴVDW (L/mol) ѴCORRECCIÓN (L/mol) %E

290,5 1811,039 1815,3 0,24
302,9 943,899 942,9 0,11
316 492,102 492,7 0,12
324,7 336,938 337,3 0,11
336 209,010 209,5 0,19
353 109,570 110,0 0,39
372 57,525 58,0 0,83
391 31,637 32,1 1,46
416,5 16,671 17,1 2,57

%EPROMEDIO=0,67

SE PUEDE NOTAR QUE HAY UNA TENDENCIA A ACERCARSE A LOS DATOS OBTENIDOS POR
EL MODELO DE VDW PARA VALORES BAJOS DE PRESIÓN. POR LO TANTO SÍ HAY UNA
TENDENCIA AL MODELO EXPONENCIAL DE 3 CONSTANTES.

AHORA, SI PUDIÉRAMOS ASUMIR EL CASO PARTICULAR DE QUE EL NÚMERO n

DEPENDIERA DE LA DISTRIBUCIÓN DE BOLTZMANN EN FUNCIÓN DE LAS ENERGÍAS DE
ESTADO, ES DECIR:

n = n0EXP((-E/K)T-1)
 QUE ES LO MISMO QUE DECIR:

n = n0EXP(-E/(KT))
POR LO CUAL ENTONCES DEBERÍA CUMPLIR QUE:

Pѵ = RT (1 – α[EXP(-E/(KT))]/T)

PARA PODER VERIFICAR TAL SITUACIÓN, REQUERIMOS DE DATOS EXPERIMENTALES QUE

PUEDAN RESPALDAR TAL SUPOSICIÓN. SIN EMBARGO AL PROBAR PARA EL CASO DE LA
RELACIÓN VISCOSIDAD-TEMPERATURA DE LÍQUIDOS POR EL MODELO DE ANDRADE
(ECUACIÓN CON UNA DISTRIBUCIÓN MUY SIMILAR A LA DE BOLTZMANN EN FUNCIÓN DE
LAS ENERGÍAS DE ESTADO), SE PUDO VERIFICAR QUE LA TENDENCIA AL MODELO DE 3
CONSTANTES MOSTRÓ RESULTADOS MÁS EXACTOS, Y POR LO TANTO, ESTE CASO
PARTICULAR ES PROBABLE QUE NO CUMPLA DE FORMA GENERALIZADA PARA TODO GAS,
TAL COMO TAMPOCO OCURRE PARA EL CASO DE LA RELACIÓN VISCOSIDAD-
TEMPERATURA DE LOS LÍQUIDOS. SE PUEDE SUPONER TAMBIÉN, QUE LA DISTRIBUCIÓN
DE BOLTZMANN ES SOLO UN CASO PARTICULAR DEL MODELO DE 3 CONSTANTES.

AHORA, SI EXPRESAMOS TODA LA EXPRESIÓN ENTRE PARÉNTESIS COMO UNA FUNCIÓN

Z(T) TENDRÍAMOS LA MISMA EXPRESIÓN DE LA ECUACIÓN VIRIAL DE ESTADO:

Pѵ = Z(T) RT

DONDE Z(T) REPRESENTA EL FACTOR DE COMPRESIBILIDAD QUE PUEDE SER

EXPRESADO COMO:

Z(T) = 1 + A/T + B/T2 + … = 1 – α(EXP(βTγ))/T

CON DATOS EXPERIMENTALES DE ESTE MODELO VIRIAL, PODEMOS HALLAR UNA

SOLUCIÓN PARA EL CÁLCULO DE NUESTRAS CONSTANTES PARA NUESTRO
MODELO.

PARTAMOS DEL TEXTO OBTENIDO DE LA SIGUIENTE PÁGINA:

http://webdelprofesor.ula.ve/ingenieria/josedel/josedel/asignatura
s-pregrado_files/Ecuaciones%20de%20Estado_may_2017.pdf

DONDE LA ECUACIÓN VIRIAL SE PUEDE EXPRESAR DE LA SIGUIENTE

FORMA:
SIN EMBARGO, A BAJAS PRESIONES (ALREDEDOR DE 100 bar =
10000 KPa), SE PUEDE USAR LA APROXIMACIÓN TRUNCADA:

DONDE B SE PUEDE HALLAR COMO:

DONDE: Tr = T/TC Y, PC Y TC SON LA PRESIÓN Y LA TEMPERATURA

CRÍTICA Y R = 8,314472 KPa.L/mol°K.

EN EL TEXTO MENCIONADO, SE ENCUENTRAN LAS SIGUIENTES

CONSTANTES PARA EL CASO DEL ETILENO (A BAJAS PRESIONES):

PC = 50,41 bar = 5041 KPa; TC =283,34°K; ω = 0,087

1. CALCULEMOS PRIMERO EL VOLUMEN MOLAR PARA UNA

PRESIÓN CONSTANTE DE 90 bar = 9000KPa, PARA LAS
TEMPERATURAS DE 310-320°K POR EL MODELO VIRIAL:

Ѵ (L/mol) Z (1-Z) X T T (°K) Tr

0,1552479772 0,5420908674 141,951831 310 1,094091904
1
0,1571140033 0,5468427040 140,931919 311 1,097621232
1
0,1589716540 0,5515348165 139,921137 312 1,101150561
3
0,1608210339 0,5561684472 138,919276 313 1,104679890
0
0,1626622463 0,5607444085 137,926255 314 1,108209219
7
0,1644953926 0,5652635982 136,941966 315 1,111738547
6
0,1663205734 0,5697268998 135,966299 316 1,115267876
7
0,1681378867 0,5741351769 134,999148 317 1,118797205
9
0,1699474291 0,5784892774 134,040409 318 1,122326533
8
0,1717492964 0,5827900358 133,089978 319 1,125855862
6
0,1735435820 0,5870382682 132,147754 320 1,129385191
2

CON 3 DE LOS DATOS DE LA TABLA HALLAMOS LAS CONSTANTES DE FORMA

COMPUTACIONAL PARA EXPRESARLO EN FUNCIÓN A NUESTRO MODELO.
INGRESAMOS PRIMERO LOS DATOS (DE LA RELACIÓN (1-Z)T VS. T):
LUEGO HALLAMOS LAS CONSTANTES:

Y VERIFICAMOS LOS RESULTADOS:

CON ESTOS DATOS VERIFICAMOS LOS RESULTADOS DEL VOLUMEN MOLAR
HALLADO CON LA APROXIMACIÓN VIRIAL PARA BAJAS PRESIONES DE LA
PÁGINA WEB DE REFERENCIA Y NUESTRO MODELO EXPRESADO COMO:

Z(T) = 1 – α(EXP(βTγ))/T

DONDE: α = 14644,39776480088; β = -0,2911594071929635; γ = 0,4824844

Y CUMPLE QUE:

Pѵ = Z(T) RT

COMPAREMOS LOS RESULTADOS DEL VOLUMEN MOLAR OBTENIDO POR

NUESTRO MODELO Y EL MODELO VIRIAL, PARA LA PRESIÓN CONSTANTE DE
9000 KPa:

Ѵ (L/mol) MODELO Ѵ (L/mol) MODELO T (°K) (PRESIÓN

 VIRIAL A BAJAS VIRIAL DE 3 CONSTANTE =
PRESIONES CONSTANTES 9000 KPa)
0,155 0,156 310
0,157 0,158 311
0,159 0,160 312
0,161 0,162 313
0,163 0,164 314
0,164 0,165 315
0,166 0,167 316
0,168 0,169 317
0,170 0,171 318
0,172 0,173 319
0,174 0,174 320

DE ESTA FORMA, PODEMOS COMPROBAR QUE NUESTRO MODELO SE

AJUSTA MUY BIEN A LOS RESULTADOS DEL MODELO VIRIAL PARA EL CASO
DE LA COMPARACIÓN DE LOS VOLÚMENES PARA UNA PRESIÓN
CONSTANTE.

2. CALCULEMOS AHORA LA PRESIÓN PARA UN VOLUMEN CONSTANTE DE

0,16266 L/mol, PARA LAS TEMPERATURAS DE 310-320°K. PRIMERO
DEBEMOS DESPEJAR LA PRESIÓN A PARTIR DE LA EXPRESIÓN:

PѴ = RT (1 + BP/(RT))

AGRUPAMOS LOS ELEMENTOS QUE CONTIENEN A P:

PѴ – BP = RT = P (Ѵ – B)

DESPEJANDO P NOS QUEDA:

P = RT/(Ѵ – B)

HALLAMOS LOS RESULTADOS DE P PARA AMBOS MODELOS:

P (KPa) MODELO P (KPa) T (°K) Tr

 VIRIAL A BAJAS MODELO DE 3
PRESIONES CONSTANTES
8772,95 8589,90 310 1,094091904
8829,56 8693,14 311 1,097621232
8886,29 8795,92 312 1,101150561
8943,12 8898,25 313 1,104679890
9000,07 9000,13 314 1,108209219
9057,12 9101,56 315 1,111738547
9114,29 9202,55 316 1,115267876
9171,56 9303,06 317 1,118797205
9228,93 9403,23 318 1,122326533
9286,41 9502,92 319 1,125855862
9344,00 9602,20 320 1,129385191

POR LO QUE SE PUEDE APRECIAR, EXISTE CIERTA CERCANÍA ENTRE LOS

DATOS DE LA PRESIÓN CALCULADA POR AMBOS MODELOS.
COMPARAREMOS ENTONCES CON LOS RESULTADOS OBTENIDOS POR EL
MODELO DE VDW, Y SI NO HAY CONCORDANCIA, PROCEDEREMOS A
CALCULAR LAS CONSTANTES DE NUESTRO MODELO PARA UN SEGUNDO
VALOR DE PRESIÓN (CONSTANTE).
DE LA PÁGINA WEB PODEMOS HALLAR LAS CONSTANTES DE VDW PARA EL
ETILENO:
file:///C:/Users/A/Downloads/SEMINARIOS%20PRIMERA%20PARTE-
2015%20(4).pdf

PARA LA ECUACIÓN DE VDW:

SI 1atm = 101,325 KPa:

a = 453,0241 KPa.L2.mol-2; b = 0,05714 L.mol-1
 COMPARAMOS LOS RESULTADOS OBTENIDOS PARA LA PRESIÓN POR EL
MODELO VIRIAL, EL MODELO DE VDW Y NUESTRO MODELO PARA UN
VOLUMEN MOLAR CONSTANTE DE 0,16266 L/mol:

P (KPa) MODELO VIRIAL P (KPa) P (KPa) MODELO DE 3 T (°K)

A BAJAS PRESIONES MODELO DE CONSTANTES
VAN DER
WAALS
8772,95 7304,31 8589,90 310
8829,56 7383,11 8693,14 311
8886,29 7461,90 8795,92 312
8943,12 7540,70 8898,25 313
9000,07 7619,49 9000,13 314
9057,12 7698,29 9101,56 315
9114,29 7777,08 9202,55 316
9171,56 7855,88 9303,06 317
9228,93 7934,67 9403,23 318
9286,41 8013,46 9502,92 319
9344,00 8092,26 9602,20 320
EN ESTA OTRA PUBLICACIÓN HALLAMOS LOS SIGUIENTES DATOS:
http://depa.fquim.unam.mx/termofisica/silvia/GR2.pdf

POR LO CUAL SI 1bar = 100 KPa, ENTONCES: a = 461,12 KPa.L2.mol-2; b = 0,058199 L.mol-1.
CALCULAMOS CON ESTOS NUEVOS DATOS:

P (KPa) MODELO VIRIAL P (KPa) P (KPa) MODELO DE 3 T (°K)

A BAJAS PRESIONES MODELO DE CONSTANTES
 VAN DER
WAALS
8772,95 7245,95 8589,90 310
8829,56 7325,55 8693,14 311
8886,29 7405,14 8795,92 312
8943,12 7484,74 8898,25 313
9000,07 7564,33 9000,13 314
9057,12 7643,92 9101,56 315
9114,29 7723,52 9202,55 316
9171,56 7803,11 9303,06 317
9228,93 7882,71 9403,23 318
9286,41 7962,30 9502,92 319
9344,00 8041,89 9602,20 320

PODEMOS NOTAR QUE LOS RESULTADOS OBTENIDOS CON ESTOS ÚLTIMOS DATOS DE LAS
CONSTANTES, QUE SON MUY SIMILARES E INCLUSO, EN EL SEGUNDO CASO, UN POCO
MENORES A LOS OBTENIDOS ANTERIORMENTE POR EL MODELO DE VDW. SIN EMBARGO,
NUESTROS RESULTADOS TIENEN UNA TENDENCIA A PROPORCIONAR VALORES MÁS
CERCANOS A LOS DEL MODELO VIRIAL. POR LO CUAL PODEMOS MENCIONAR QUE
NUESTRO MODELO PROPORCIONA RESULTADOS COHERENTES. VERIFIQUEMOS ENTONCES
AHORA PARA UNA MAYOR PRESIÓN (CONSTANTE) PARA EL MODELO VIRIAL Y PODREMOS
COMPARAR LUEGO CON NUESTRO MODELO.
3. CALCULEMOS AHORA EL VOLUMEN MOLAR PARA UNA
PRESIÓN CONSTANTE DE 100 bar = 10000KPa, PARA LAS
TEMPERATURAS DE 310-320°K POR EL MODELO VIRIAL:

Ѵ (L/mol) Z (1-Z) X T T (°K) Tr

0,1266092403 0,4912120747 157,724256 310 1,094091904
9
0,1283828834 0,4964917786 156,591056 311 1,097621232
9
0,1301481511 0,5017053516 155,467930 312 1,101150561
3
0,1319051479 0.5068538301 154,354751 313 1,104679890
2
0,1336539773 0,5119382317 153,251395 314 1,108209219
2
0,1353947406 0,5169595538 152,157740 315 1,111738547
 6
0,1371275384 0,5219187776 151,073666 316 1,115267876
3
0,1388524687 0,5268168633 149,999054 317 1,118797205
3
0,1405696281 0,5316547527 148,933788 318 1,122326533
6
0,1422791123 0,536433373 147,877754 319 1,125855862
0
0,1439810149 0,5411536314 146,830838 320 1,129385191
0

CON 3 DE LOS DATOS DE LA TABLA HALLAMOS LAS CONSTANTES DE FORMA

COMPUTACIONAL PARA EXPRESARLO EN FUNCIÓN A NUESTRO MODELO.
INGRESAMOS PRIMERO LOS DATOS (DE LA RELACIÓN (1-Z)T VS. T):

LUEGO HALLAMOS LAS CONSTANTES:

Y VERIFICAMOS LOS RESULTADOS:
CON ESTOS DATOS VERIFICAMOS LOS RESULTADOS DEL VOLUMEN MOLAR
HALLADO CON LA APROXIMACIÓN VIRIAL PARA BAJAS PRESIONES DE LA
PÁGINA WEB DE REFERENCIA Y NUESTRO MODELO EXPRESADO COMO:

Z(T) = 1 – α(EXP(βTγ))/T

DONDE: α = 16526,72605140074; β = -0,294833089125136; γ = 0,4808827

Y CUMPLE QUE:

Pѵ = Z(T) RT

COMPAREMOS LOS RESULTADOS DEL VOLUMEN MOLAR OBTENIDO POR

NUESTRO MODELO Y EL MODELO VIRIAL, PARA LA PRESIÓN CONSTANTE DE
10000 KPa:

Ѵ (L/mol) MODELO Ѵ (L/mol) MODELO T (°K) (PRESIÓN

VIRIAL A BAJAS VIRIAL DE 3 CONSTANTE =
PRESIONES CONSTANTES 10000 KPa)
0,1266092 0,1266091 310
0,1283829 0,1283827 311
0,1301481 0,1301480 312
0,1319051 0,1319050 313
0,1336540 0,1336539 314
0,1353947 0,1353947 315
0,1371275 0,1371275 316
0,1388525 0,1388525 317
0,1405696 0,1405697 318
0,1422791 0,1422791 319
0,1439810 0,1439809 320

DE ESTA FORMA, PODEMOS COMPROBAR NUEVAMENTE QUE NUESTRO

MODELO SE AJUSTA MUY BIEN A LOS RESULTADOS DEL MODELO VIRIAL
PARA EL CASO DE LA COMPARACIÓN DE LOS VOLÚMENES PARA UNA
PRESIÓN CONSTANTE. AHORA COMPAREMOS LAS CONSTANTES OBTENIDAS
A ESTA PRESIÓN CONSTANTE, CON LAS CONSTANTES OBTENIDAS
ANTERIORMENTE. PODEMOS NOTAR QUE LAS CONSTANTES β Y γ SE
MANTIENEN PRÁCTICAMENTE CONSTANTES, POR LO TANTO, PODEMOS
PROMEDIARLAS APROXIMÁNDOLAS A 3 CIFRAS DECIMALES:

β = (-0,2912 + -0,2948)/2 = -0,293

γ = (0,4825 + 0,4809)/2 = 0,482

EXPRESAMOS α COMO UNA FUNCIÓN DE LA PRESIÓN:

α1 = X + 9000Y = 14644,39776480088
α2 = X + 10000Y = 16526,72605140074

HALLAMOS LAS CONSTANTES X Y Y RESOLVIENDO EL SISTEMA DE 2

ECUACIONES:

X = -2296,557
Y = 1,882

SIENDO:

α = 1,882P – 2296,557

POR LO CUAL PODEMOS EXPRESAR NUESTRO MODELO ADAPTADO A LOS

DATOS DEL MODELO VIRIAL, COMO:

γ
PѴ = RT (1 – (1,882P – 2296,557)(EXP(βT ))/T)

ENTONCES:

PѴ + (1,882R)(EXP(βTγ))P = RT (1 + (2296,557)(EXP(βTγ))/T)

DESPEJANDO P TENEMOS:

P = RT [1 + (2296,557)(EXP(βTγ))/T]/(Ѵ + (1,882R)(EXP(βTγ)))
 COMPARAMOS LOS RESULTADOS OBTENIDOS PARA LA PRESIÓN POR EL
MODELO VIRIAL, EL MODELO DE VDW Y NUESTRO MODELO ADAPTADO AL
MODELO VIRIAL PARA UN VOLUMEN MOLAR CONSTANTE DE 0,16266
L/mol:

P (KPa) MODELO VIRIAL A BAJAS P (KPa) MODELO DE VAN P (KPa) MODELO DE 3 CONSTANTES T (°K)
PRESIONES DER WAALS ADAPTADO A LOS DATOS VIRIALES

8772,95 7304,31 8848,28 310

8829,56 7383,11 8901,40 311
8886,29 7461,90 8954,66 312
8943,12 7540,70 9007,93 313
9000,07 7619,49 9061,57 314
9057,12 7698,29 9115,22 315
9114,29 7777,08 9168,86 316
9171,56 7855,88 9222,89 317
9228,93 7934,67 9276,79 318
9286,41 8013,46 9331,06 319
9344,00 8092,26 9385,34 320

SE PUEDE APRECIAR QUE LOS RESULTADOS DE NUESTRO MODELO, ADAPTADO A LOS

DATOS DEL MODELO VIRIAL, ESTÁN MUY CERCANOS A LOS OBTENIDOS POR EL MODELO
VIRIAL TRUNCADO. LA VENTAJA DE NUESTRO MODELO ADAPTADO ES QUE SOLO
REQUIERE DE 4 CONSTANTES PARA HACER EL CÁLCULO, A DIFERENCIA DEL MODELO
VIRIAL TRUNCADO QUE REQUIERE DE 10 CONSTANTES. NUESTRO MODELO DE FORMA
RESUMIDA SE EXPRESA COMO:

PѴ = RT (1 – (1,882P – 2296,557)(EXP(-0,293T0,482))/T)

CON ESTA ECUACIÓN CALCULEMOS LOS VOLÚMENES A LAS PRESIONES CONSTANTES DE

9000 KPa Y 10000 KPa:

Ѵ (L/mol) MODELO VIRIAL A Ѵ (L/mol) MODELO VIRIAL DE 3 T (°K) (PRESIÓN

BAJAS PRESIONES CONSTANTES CONSTANTE = 9000 KPa)
0,155 0,157 310
0,157 0,159 311
0,159 0,161 312
0,161 0,163 313
0,163 0,165 314
0,164 0,167 315
0,166 0,168 316

0,168 0,170 317

0,170 0,172 318
0,172 0,174 319
 0,174 0,176 320

Ѵ (L/mol) MODELO Ѵ (L/mol) MODELO T (°K) (PRESIÓN

VIRIAL A BAJAS VIRIAL DE 3 CONSTANTE =
PRESIONES CONSTANTES 10000 KPa)
0,1266 0,1267 310
0,1284 0,1285 311
0,1301 0,1303 312
0,1319 0,1320 313
0,1337 0,1338 314
0,1354 0,1351 315
0,1371 0,1373 316
0,1389 0,1390 317
0,1406 0,1407 318
0,1423 0,1424 319
0,1440 0,1441 320

CON ESTA ECUACIÓN VERIFIQUEMOS LOS VOLÚMENES A LA PRESIÓN CONSTANTE DE

8000 KPa:

Ѵ (L/mol) Ѵ (L/mol) T (°K) (PRESIÓN

MODELO VIRIAL A MODELO VIRIAL CONSTANTE =
BAJAS PRESIONES DE 3 8000 KPa)
CONSTANTES
0,191 0,196 310
0,193 0,198 311
0,195 0,200 312
0,197 0,202 313
0,199 0,203 314
0,201 0,205 315
0,203 0,207 316
0,205 0,209 317
0,207 0,211 318
0,209 0,213 319
0,210 0,215 320

PODEMOS COMPROBAR QUE LOS RESULTADOS TIENEN COHERENCIA CON EL AJUSTE

HECHO A α COMO UNA FUNCIÓN LINEAL DE LA PRESIÓN. COMPROBEMOS AHORA PARA
EL CASO EN EL QUE LA PRESIÓN CONSTANTE FUERA DE 11000 KPa:
 Ѵ (L/mol) Ѵ (L/mol) T (°K) (PRESIÓN
MODELO VIRIAL A MODELO VIRIAL CONSTANTE =
BAJAS PRESIONES DE 3 11000 KPa)
CONSTANTES
0,103 0,102 310
0,105 0,103 311
0,107 0,105 312
0,108 0,107 313
0,110 0,108 314
0,112 0,110 315
0,113 0,112 316
0,115 0,113 317
0,117 0,115 318
0,118 0,117 319
0,120 0,118 320

PROBEMOS UN PUNTO INTERMEDIO PARA EL CASO DE LA PRESIÓN

CONSTANTE DE 9500 KPa:

Ѵ (L/mol) Ѵ (L/mol) T (°K) (PRESIÓN

MODELO VIRIAL A MODELO VIRIAL CONSTANTE =
BAJAS PRESIONES DE 3 9500 KPa)
CONSTANTES
0,140 0,141 310
0,142 0,143 311
0,144 0,145 312
0,146 0,147 313
0,148 0,148 314
0,149 0,150 315
0,151 0,152 316
0,153 0,154 317
0,154 0,156 318
0,156 0,157 319
0,158 0,159 320
PROBEMOS POR ÚLTIMO PARA LAS PRESIONES CONSTANTES DE 7000 KPa Y 12000 KPa:

Ѵ (L/mol) MODELO VIRIAL A Ѵ (L/mol) MODELO VIRIAL T (°K) (PRESIÓN

BAJAS PRESIONES DE 3 CONSTANTES CONSTANTE = 7000 KPa)
0,237 0,245 310
0,239 0,247 311
0,241 0,249 312
0,243 0,251 313
0,246 0,253 314
0,248 0,255 315
0,250 0,257 316
0,252 0,259 317
0,254 0,261 318
0,256 0,263 319
0,258 0,265 320

Ѵ (L/mol) MODELO VIRIAL A Ѵ (L/mol) MODELO VIRIAL T (°K) (PRESIÓN

BAJAS PRESIONES DE 3 CONSTANTES CONSTANTE = 12000
KPa)
0,084 0,081 310
0,085 0,082 311
0,087 0,084 312
0,089 0,086 313
0,090 0,087 314
0,092 0,089 315
0,093 0,091 316
0,095 0,092 317
0,097 0,094 318
0,098 0,095 319
0,100 0,097 320

DE ESTA FORMA, PODEMOS COMPROBAR QUE HEMOS PODIDO AJUSTAR LOS DATOS DEL
MODELO VIRIAL A NUESTRO MODELO CON GRAN PROXIMIDAD DE LOS RESULTADOS
OBTENIDOS.

MODELO DE 3 CONSTANTES PARA LA RELACIÓN TEMPERATURA-TIEMPO

AHORA, SABIENDO QUE LA APROXIMACIÓN GENERALIZADA ERA DE LA FORMA:

Pѵ = RT (1 – (nђ/(ηRt))(1/T))
 SI YA HEMOS COMPROBADO QUE CUMPLE LA EXPRESIÓN:

Pѵ = RT (1 – α(EXP(βTγ))/T)

COMPARANDO AMBAS EXPRESIONES TENEMOS QUE:

1 – (nђ/(ηRt))(1/T) = 1 – α(EXP(βTγ))/T

LO CUAL NOS DEMUESTRA QUE TEMPERATURA Y TIEMPO ESTÁN RELACIONADAS, DE

MODO QUE SI LOGRÁRAMOS DESPEJAR T EN FUNCIÓN DE t PODRÍAMOS REPRESENTARLA
COMO UNA APROXIMACIÓN POR NUESTRO MODELO DE 3 CONSTANTES DE MODO QUE:

T = T0 (EXP (ktn))

PROBEMOS EXPERIMENTALMENTE, PARA UN SISTEMA CERRADO Y AISLADO APROXIMADAMENTE,

INTRODUCIENDO UN VASO CON AGUA DESTILADA EN UN PLANCHA DE CALENTAMIENTO CUYA
TEMPERATURA SE MANTIENE CASI CONSTANTE EN EL TIEMPO. LA TEMPERATURA DEL AGUA
VARÍA EN EL TIEMPO EXPERIMENTALMENTE DE LA SIGUIENTE FORMA:

t (segundos) T (°C) T (°K)

0 27 300
140 42 315
168 46 319
193 50 323
208 52 325
229 56 329

CALCULAMOS DE FORMA NUMÉRICA PARA 3 DATOS:

CALCULAMOS LAS CONSTANTES:

VERIFICAMOS LOS RESULTADOS:

DE ESTA FORMA COMPROBAMOS LA RELACIÓN ENTRE LA TEMPERATURA Y EL TIEMPO, LA CUAL
SE PUEDE EXPRESAR POR NUESTRO MODELO EXPONENCIAL. LAS CONSTANTES OBTENIDAS SON:

T0 = 303,752°K; k = 1,132X10-5; n = 1,634

SI EXPRESAMOS A: k = 1/tEn, DESPEJAMOS: tE = (1/k)(1/n) = 1064,196s

DE ESTE MODO PODEMOS EXPRESAR LA TEMPERATURA EN FUNCIÓN DEL TIEMPO COMO:

T = T0 (EXP (t/tE)n)

DONDE: T0 = 303,752°K; tE = 1064,196s; n = 1,634

PARA PODER EXPRESARLO DE FORMA DIFERENCIAL, DE LA EXPRESIÓN INICIAL:

T = T0 (EXP (ktn))

MULTIPLICAMOS POR 1/T0:

T/T0 = EXP (ktn)

SACAMOS LOGARITMO NATURAL:

 ln(T/T0) = ln(EXP (ktn))

LO QUE ES LO MISMO QUE:

lnT -lnT0 = ktn – k(0)n

QUE SE PUEDE EXPRESAR COMO:

lnT = ktn + cte.

DE LA FORMA DIFERENCIAL:

dT/T = nktn-1dt

SI ASUMIMOS A Kα = nk Y u = nk:

dT = Kα T tu dt

SIENDO ESTA ÚLTIMA EXPRESIÓN, LA FORMA DIFERENCIAL DE LA RELACIÓN

TEMPERATURA-TIEMPO.

POR ÚLTIMO, PODEMOS MOSTRAR OTRA FORMA DE PODER HALLAR LAS

CONSTANTES DE NUESTRO MODELO PARA ESTE CASO. PARA ESTO, APLICAMOS
PARA 3 DATOS NUESTRO MODELO:

t=0, ENTONCES: T=T0 EXP(k(0)n) = T0 = 300°K (HALLAMOS UNA DE LAS CONSTANTES) …(i)

t=140s, ENTONCES: T=(300°K) EXP(k(140)n) = 315°K, POR LO CUAL: (315/300) = EXP(k(140)n) …(ii)

t=168s, ENTONCES: T=(300°K) EXP(k(168)n) = 319°K, POR LO CUAL: (319/300) = EXP(k(168)n) …(ii)

APLICAMOS LOGARITMO NATURAL A LOS RESULTADOS OBTENIDOS EN (ii) Y (iii):

ln(315/300) = k(140)n …(iv)

ln(319/300) = k(168)n …(v)

DIVIDIMOS LO RESULTADOS OBTENIDOS EN (iv) Y (v):

(ln(315/300))/(ln(319/300)) = (140/168)n …(vi)

APLICO EL LOGARITMO NATURAL Y DESPEJO n:

n = ln[(ln(315/300))/(ln(319/300))]/ln[(140/168)] = 1,262

DE LA ECUACIÓN (iv) DESPEJO k:

 k = (ln(315/300))/(140)n = 9,566 x 10-6

COMPROBAMOS PARA LOS DEMÁS DATOS APLICANDO ESTAS CONSTANTES:

t (segundos) T (°K) T (°K)

CALCULADA EXPERIMENTAL
0 300,0 300
140 315,0 315
168 319,0 319
193 322,8 323
208 325,2 325
229 328,6 329

DE ESTA FORMA, PODEMOS COMPROBAR QUE LA RELACIÓN TEMPERATURA-TIEMPO SE AJUSTA A

NUESTRO MODELO DE 3 CONSTANTES.