INSTITUTO POLITÉCNICO NACIONAL

“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”

DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela
Tarea 1. Ecuaciones no Lineales
Problema 1.- Calcule el volumen específico molar del metano a 200ºC para los
valores de presión dados en la siguiente tabla:

P 5 5 5 6 6 7 7 8 8 9 9 9 10 12
0 5 8 2 5 0 7 3 7 4 6 9 0 0
V

La presión está dada en atm y el volumen en lt. Utilice la ecuación de estado de

Gas Ideal como valor inicial y luego encuentre el volumen utilizando para ello la
ecuación de estado de Van der Waals:

( P+ v́a ) ( v́ −b)=RT
2

atm∗¿ ¿2∗atm
Con los valores de n = 1mol,R=0.08205 con a=2.2732 y
gmol∗K mol 2
b=0.04306 ¿ con el Método de Newton Raphson con un número máximo de
mol
iteraciones de 10 y  = 1.0E-4.

SOLUCION:

1) La ecuación de gas ideal dará el valor inicial de volumen para cada presión, de
esa manera podremos hacer la iteración del volumen por medio de un valor
conocido
PV =nRT ( Ecuación de gas ideal)

2) Obtendremos la función , despejando la ecuación de estado de Van der Waals:

a
( )
f ( x )= P+ 2 ( v́−b )−RT
v́

3) Para derivar con respecto a la variable “v” del volumen molar, se necesita despejar
para dejar la ecuación en términos de “v”

( P+ v́a ) ( v́ −b)=RT
2

V́ a ab
P V́ + 2
−Pb− 2 =RT
V́ V́
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela

Para quitar denominadores , multiplicamos ambos lados por V́ 2

V́ a ab
( P V́ +
V́ 2
V́ )
−Pb− 2 =RT ∗V́ 2

P V́ 3+ V́ a−Pb V́ 2−ab=RT V́ 2

P V́ 3−RT V́ 2 + V́ a−Pb V́ 2−ab=0

P V́ 3−V́ 2 (Pb+ RT )+ V́ a−ab=0

( P V́ 3−V́ 2 ( Pb+ RT )+ V́ a−ab=0 )∗1

P

V́ 2 ( Pb+ RT ) V́ a ab
V́ 3− + − =0
P P P

RT V́ a ab
V́ 3−V́ 2 b+
P (+
P
− =0
P )
Simplificando la ecuación:

RT
A=− b+ ( P )
a
B=
P
−ab
C= =−Bb
P

Sustituyendo los términos en la ecuación, para obtener la función:

f ( v )=V́ 3− A V́ 2 + B V́ +C

Derivando la ecuación:
dP
=3 V́ 2−2 A V́ + B
dV

 Sustituimos para el primer valor de presión

Datos:
P=50 atm T =475.15 K
2
¿ ∗atm atm∗¿
a=2.2732 R=0.08205
mol2
gmol∗K
b=0.04306 ¿ n=1 mol
mol
4) Obtener el volumen molar de la presión a 50 atm, por la ecuación de gas ideal:
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela
atm∗¿
( 1mol)(0.08205 )(475.15 K )
gmol∗K =0.77643915lt
V=
50 atm

Sustituir A , B y C :

atm∗¿
A=− 0.04306( ¿
mol
+
0.08205
gmol∗K

2
50 atm
∗( 475.15 K )
)
=-0.81949915

¿ ∗atm
2.2732
mol2 =0.045464
B=
50 atm
C=−0.045464∗0.04306 ¿ =0.00195768
mol
5) Obtener f ( v ) , con el valor inicial del volumen:

f ( v o ) =V́ 3− A V́ 2 +B V́ +C
f ( v o ) =¿´¿ ¿0.00195768
f ( v o ) =¿0.011298654

dP
6) Obtener la derivada ( ¿ del volumen molar dado
dV
dP
=3(0.77643915)2−2(−0.81949915) ( 0.77643915 ) +0.045464
dV

dP
=¿0.58145481
dV
7) Aplicando el método de Newton Raphson
f ( v́ )
v1 =v́ 0−
dp
dv
0.011298654
v1 =−0.77643915 0
´ − =0.75700745
0.58145481

8) Cálculo de la norma:
2
norma= ( V́ ¿ ¿ k +1− V́ k ) = √(0.75700745−0.77643915)2=¿ ¿ 0.0194317
√
9) TABLA DE RESULTADOS
P 50 55 58 62 65 70 77 83 87 94 96 99 100 120
V 0.75 0.68 0.64 0.60 0.57 0.53 0.48 0.44 0.42 0.38 0.37 0.36 0.360 0.292
59 48 79 42 49 16 02 30 08 65 75 48 7 1

10) GRAFICA DE REUSLTADOS:

 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela

PV
0.8
0.7
0.6
Volumen molar

0.5
0.4
0.3
0.2
0.1
0
40 50 60 70 80 90 100 110 120 130
Presion (atm)

Con los resultados obtenidos observamos que se cumple la norma ya que van
disminuyendo los valores de la norma conforme se va iterando, hasta obtener un
resultado correcto de volumen para cada presión en el sistema.
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela

Problema 2.- Para calcular el radio hidráulico r en un canal abierto, Francis

propuso la siguiente ecuación:
v=c √ ℜ---(1)
75
c=
m ---(2)
0.45+
√r
Donde m = 1.05, e = 0.002 y v = 5.5 Realice un máximo de 10 iteraciones con el
método de Newton-Raphson para aproximar una raíz de la ecuación con el valor
inicial r(0) = 4.1 y ε = 0.1

SOLUCION:
De la ecuación 1 se despeja a c, para sustituirlo en la ecuación 2 y dejar todo en términos
conocidos dejando solo una variable sin conocer “v” e igualar la ecuación 2 a cero.
Despeje de ecuación 1 en termino de c
v
c= ----(3)
√ℜ

Sustituyendo (3) en (2)

v 75
=
√ℜ m
0.45+
√r
m
(
v 0.45+ )
√r
−75 √ ℜ=0

m
(
f ( v )=v 0.45+
√r )
−75 √ ℜ

Derivar
d −1
= ¿
dv 2
Simplificando:
d
=−¿
dv
d −2.8875 37.5 √ 0.002
= −
dv √r3 √r
3) El valor inicial es dado r=4.1

4) sustituimos el valor de r para la función

1.05
(
f ( v )=( 5.5 ) 0.45+
√ 4.1 )
−75 √ ( 4.1 ) ( 0.002 )=¿-1.46446964

5) sustituimos para l derivada

d −2.8875 37.5 √ 0.002
= − =−1.17604976
dv √ 4.13 √ 4.1
6) aplicando el método:
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela
f ( v́ )
r 1=ŕ 0−
d
dr
−1.46446964
v1 =4.1− =2.85475538
−1.17604976
7) Cálculo de la norma:
2
norma= (ŕ ¿ ¿ k +1−ŕ k ) =√(2.85475538−4.1)2 =1.24524462¿
√

TABLA DE RESULTADOS:

i r f(r) df/dr r^k+1 norma

0 4.1 -1.46446964 -1.17604976 2.85475538 1.24524462
1 2.85475538 0.22586658 -1.59121555 2.99670132 0.14194594
2 2.99670132 0.00475219 -1.52539603 2.9998167 0.00311538
3 2.9998167 2.1349E-06 -1.52402599 2.9998181 1.4008E-06
4 2.9998181 4.3077E-13 -1.52402537 2.9998181 2.8244E-13
5 2.9998181 1.7764E-15 -1.52402537 2.9998181 1.3323E-15
6 2.9998181 -1.7764E-15 -1.52402537 2.9998181 1.3323E-15
7 2.9998181 1.7764E-15 -1.52402537 2.9998181 1.3323E-15
8 2.9998181 -1.7764E-15 -1.52402537 2.9998181 1.3323E-15
9 2.9998181 1.7764E-15 -1.52402537 2.9998181 1.3323E-15
10 2.9998181 -1.7764E-15 -1.52402537 2.9998181 1.3323E-15

En la tabla de resultados, observamos que el valor de error solicitado en el problema era

de 0.1 lo cual se llegó desde la segunda iteración, sin embargo, el valor ya no cambia
desde la cuarta iteración con lo que se selecciona como resultado final, el radio es de
2.9998181.
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela
Problema 3.- Una mezcla equimolar de Monóxido de carbono y Oxígeno alcanza el
equilibrio a 300°K y una presión de 5 atm. La reacción teórica es:
CO + 0.5O2  CO2
La reacción química real se escribe como:

CO + O2  xCO 1
+ (1 + x )O2 + (1 – x )CO2
2

La ecuación de equilibrio químico para determinar la fracción del CO restante x

se describe como:

(1−x)(3+ x )1/ 2
Kp= −−Ec(1)0< X <1
x ( x+ 1)1 /2 P 1/ 2

Donde Kp = 3.06 es la constante de equilibrio para CO + 0.5O2  CO2 a 300°K

y P = 5 atm es la presión. Calcule el valor de x por medio del Método de Newton-
Raphson hasta obtener un error relativo de 0.0001.
Solución : 0.1930
SOLUCION:
1) Despejar la ecuación (1) a cero, para obtener la función con la que se trabajará
1 1 1
( )
Kp x ( x +1 ) 2 ( P ) 2 =( 1−x )( 3+ X ) 2
1 1 1
Kp ( x ( x +1 ) ( P ) )−( 1−x )( 3+ X ) =0
2 2 2

1 1 1
f ( x )=Kp ( x ( x +1 ) ( P ) )−( 1−x )( 3+ X )
2 2 2

2) Derivar la función obtenida en función a la variable x

1 1 1
( )
Kp x ( x +1 ) 2 ( P ) 2 −( 1−x )( 3+ X ) 2 =0

1
1 1
( P )2 ∗d d
f ( x )=Kp
dx
( dx
)
x ( x +1 ) 2 − ( 1−x )( 3+ X ) 2 ----Derivada de la regla de la cadena

−1
1 1
1 ( 1−x ) ( 3+ X ) 2
−1
2
(
f ´ ( x )=Kp ( P ) ∗ x ( x+1 ) + ( x +1 ) −
2
2
2
2
)
1 1 1
x ( 1−x )
f ´ ( x )=Kp ( P ) ∗ 2

( 2 ( x+ 1 )
1
2
2
+ ( x+1 ) −
) 2 ( 3+ X )
1
2
−( 3+ X ) 2

Kp √ P(x ) ( 1−x )
f ´ ( x )=Kp √ P √ x +1+ ( 2 √ x+1
− )
2 √ 3+ x
−√ 3+ x
 INSTITUTO POLITÉCNICO NACIONAL
“ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS”
DEPARTAMENTO DE INGEMNIERÍA QUÍMICA PETROLERA
Alumna: Hernández Quintana Ana Gabriela

TABLA DE RESULTDOS:
Para estos resultados la variable fue cambiando su valor entre 0 y 1 según lo
escrito en el problema para la constante de equilibrio en función a la presión
i X fx d/dx   x^k+1 norma
0 0.2 0.06800422 4.23373766 6.10759732 0.18886564 0.01113436
1 0.18886564 -0.03942679 4.26252625 6.04033174 0.19539289 0.00652725
2 0.19539289 0.02345284 4.24560108 6.07981975 0.1915354 0.00385749
3 0.1915354 -0.01374192 4.25558682 6.05650204 0.19380435 0.00226895
4 0.19380435 0.00812392 4.24970744 6.07022401 0.19246603 0.00133832
5 0.19246603 -0.00477758 4.25317332 6.06213251 0.19325413 0.0007881
6 0.19325413 0.00281833 4.25113165 6.06689817 0.19278959 0.00046454
7 0.19278959 -0.00165953 4.25233486 6.06408936 0.19306326 0.00027367
8 0.19306326 0.00097824 4.25162596 6.06574415 0.19290198 0.00016127
9 0.19290198 -0.00057627 4.25204369 6.064769 0.192997 9.502E-05
10 0.192997 0.00033961 4.25179756 6.06534356 0.19294101 5.5991E-05

En este caso no teníamos ningún valor fijo para el valor x, así que se fue
calculando respecto al valor en el que tenia que estar, sin embargo, este resultado
igual se podía hacer con una grafica para llegar a la x correcta para calcular, pero
cualquiera de los dos métodos resulta de la misma manera.