Scribd
Cargar
401 vistas6 páginas

Cálculo de disociación del agua usando Newton-Raphson

El documento describe un problema de ingeniería química sobre la disociación del agua a altas temperaturas en hidrógeno y oxígeno. Se utiliza el método de Newton-Raphson para calcular la fracción molar x de H2O que se disocia basándose en la ecuación de equilibrio químico dada una constante de equilibrio K = 0.05 y una presión total de 3 atmósferas. El método iterativo arroja como solución una fracción molar x = 0.028249.

Cargado por

Ximena Ac
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
401 vistas6 páginas

Cálculo de disociación del agua usando Newton-Raphson

El documento describe un problema de ingeniería química sobre la disociación del agua a altas temperaturas en hidrógeno y oxígeno. Se utiliza el método de Newton-Raphson para calcular la fracción molar x de H2O que se disocia basándose en la ecuación de equilibrio químico dada una constante de equilibrio K = 0.05 y una presión total de 3 atmósferas. El método iterativo arroja como solución una fracción molar x = 0.028249.

Cargado por

Ximena Ac
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

PROBLEMA 8.

3
En un proceso de ingeniería química el vapor de agua (H2O) se calienta a temperaturas lo
suficientemente altas para que una porción significativa del agua se disocie, o se rompa, para formar
oxígeno (O2) e hidrógeno (H2):

H2O ↔ H2 + ½ O2

Si se supone que ésta es la única reacción que se lleva a cabo, la fracción molar x de H2O que se
disocia se representa por

x 2 pt
K=
1−x √ 2+ x

donde K = la constante de equilibrio de la reacción y pt = la presión total de la mezcla. Si pt = 3 atm


y K = 0.05, determine el valor de x que satisfaga la ecuación.

SOLUCIÓN
*Reemplazar los valores dados

x 6
0.05=
1−x 2+ x
x 6
√
f ( x )=
1−x 2+ x √
−0.05=0
*Método Newton Raphson
El cálculo se hará con dos cifras significativas
Es=0.5(10)2−5=0.0005 %

− √6 (0.5 x 2−0.5 x −2)


f ' ( x )=
( 1−x )2 (2+ x) √ 2+ x

f ( x i)
x i+1=x i−
f ' ( x i)
Valor inicial de x será 0, puesto que, al principio de la disociación, es decir en el tiempo 0, la
fracción molar de H2 será 0.
a) Para i=0
f ( x0 )
x 1=x 0−
f ' ( x0 )

0 6
x 1=0−
1−0 2+ 0
−0.05
√
−√ 6(0.5 ( 0 )2−0.5(0)−2)
( 1−0 )2 (2+0) √ 2+0

−0.05
x 1=0− =0.028868
1.73
0.02887−0
Ea= .100=100 %
0.02887
b) Para i=1

0.028868 6
x 2=0.028868−
√
1−0.028868 2+0.028868
2
−0.05
=0.028250
−√ 6 ( 0.5 ( 0.028868 ) −0.5 ( 0.028868 )−2 )
( 1−0.028868 )2 ( 2+0.028868 ) √ 2+0.028868
0.028250−0.028868
Ea= .100=−2.18682%
0.028250

c) Para i=2
0.02825 6
x 3=0.02825−
√
1−0.02825 2+ 0.02825
2
−0.05
=0.028249
−√ 6 ( 0.5 ( 0.02825 ) −0.5 ( 0.02825 ) −2 )
( 1−0.02825 )2 ( 2+0.02825 ) √ 2+ 0.02825

0.028249−0.0 2825
Ea= .100=−0.0035399 %
0.028249
d) Para i=3
0.028249 6
x 4 =0.028249−
√
1−0.028249 2+0.028249
2
−0.05
=0.028249
−√ 6 ( 0.5 ( 0.028249 ) −0.5 ( 0.028249 )−2 )
( 1−0.028249 )2 ( 2+0.028249 ) √ 2+0.028249

0.028249−0.028249
Ea= .100=0 %
0.028249
CÓDIGO

function newt()
n=0; Eap=100
vo=input('Ingrese el valor inicial: ');
E=input('Ingrese el error preespecificado: ');
imax=input('Ingrese el número máximo de iteraciones: ');
deff('y=f(x)','y=(x/(1-x))*sqrt(6/(2+x))-0.05');
xi=vo-(f(vo)/numderivative(f,vo));
Eap=abs(((xi-vo)/xi)*100);
printf('i\t\t xo\t\t\t xi\t\t\t Error\n');
printf('%d\t\t %3.6f\t\t %3.6f\t\t %3.6f\n',n,vo,xi,Eap);
while (Eap>E && n<imax)
vo=xi;
xi=vo-(f(vo)/numderivative(f,vo));
Eap=abs(((xi-vo)/xi)*100);
n=n+1
printf('%d\t\t %3.6f\t\t %3.6f\t\t %3.6f\n',n,vo,xi,Eap);
if (Eap<=E | n==imax)
break
end
end
printf('El valor de la fracción molar es : %3.6f\n',xi);
x=-3:0.2:10;
plot(x,f);
title('Disociacion del agua','fontsize',6);
xlabel('x','interpreter','latex','fontsize',3)
ylabel('y','interpreter','latex','fontsize',3)
xgrid(2)

endfunction
EJECUCIÓN
CONCEPTOS
Reacción de descomposición:
Las reacciones de descomposición son lo opuesto de las reacciones de combinación. Concretamente,
una reacción de descomposición es la ruptura de un compuesto en dos o más componentes.

Método de Newton-Raphson:
Es un método ampliamente usado, se basa en el concepto de la derivada de una función. Si el valor
inicial para la raíz es xi, entonces se puede trazar una tangente desde el punto [xi, f(xi)]. La primera
derivada será:

f ( x i)
f ' ( x i )=
x i−x i+1

Despejando:

f ( x)
x i+1=x i−
f ' (x)

La cual se conoce como fórmula de Newton-


Raphson

PROBLEMA 8.3 En un proceso de ingeniería química el vapor de agua (H2O) se calienta a temperaturas lo suficientemente altas
b) Para i=1 x2=0.028868− 0.028868 1−0.028868√ 6 2+0.028868−0.05 −√6(0.5 (0.028868) 2−0.5 (0.028868)−2) (1−0.028868) 2 (2+0.02
CÓDIGO function newt() n=0; Eap=100 vo=input('Ingrese el valor inicial: '); E=input('Ingrese el error preespecifi
EJECUCIÓN
CONCEPTOS Reacción de descomposición: Las reacciones de descomposición son lo opuesto de las reacciones de combinación. Conc

También podría gustarte