Aplicación de Las Transformadas de Laplace

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 6

Nombre: Segoviano Villegas Diego Aaron

Número de Matrícula: 311198535

Aplicación de las transformadas de Laplace en Ingeniería


Química.

En este trabajo se tiene la solución de una reacción dada por dos métodos, lo
primero fue proponer la cinética de dicha reacción y posteriormente resolverla por
medio de una simulación en el programa Scilab con la herramienta Xcos y
aplicando las Transformadas de Laplace simultáneamente, teniendo que
observarse que en ambos casos los resultados son iguales.

Cinética asignada:

Considerando las condiciones iniciales siguientes.

t=0 min , C A =C A 0=10 mol / LC B =0 CC =0 C D=0

Las ecuaciones de velocidad de reacción para cada componente y las constantes


k que se han propuesto son las siguientes:

dC A
=−k 1 C A
dt
k 1=0.98
dC B
=k 1 C A −k 2 C B −k 3 C B k 2=0.53
dt
dC C k 3=0.50
=k 3 C B −k 4 C C +k 5 C D
dt k 4=0.31
dC D k 5=0.1 1
=k 2 C B + k 4 CC −k 5 C D
dt
Simulando en Scilab con la herramienta Xcos:

Obtenemos el siguiente gráfico:


Comportamiento de las variables (CA, CB, CC y CD), al resolver el sistema de
ecuaciones diferenciales propuesto en Xcos.
Ahora, se procedió a resolver el sistema de ecuaciones con transformadas de La
Place y posteriormente se graficarán las funciones resultantes en Excel.
Aplicando la Transformada de La Place:

s C A ( s ) −C A ( 0 )=−k 1 C A ( s )

s C B ( s )−C B ( 0 )=k 1 C A ( s )−k 2 C B ( s )−k 3 C B ( s )

s CC ( s )−C C ( 0 )=k 3 C B ( s )−k 4 C C ( s ) + k 5 C D ( s )

s C D ( s )−C D ( 0 )=k 2 C B ( s )+ k 4 C C ( s ) −k 5 C D ( s )

Reordenando y asignando valores iniciales:

s C A ( s ) + k 1 C A ( s )=10

s C B ( s )−k 1 C A ( s ) +k 2 C B ( s )+ k 3 CB ( s ) =0

s CC ( s )−k 3 C B ( s ) +k 4 CC ( s )−k 5 C D ( s ) =0

s C D ( s )−k 2 C B ( s )−k 4 CC ( s ) +k 5 C D ( s )=0

Factorizando los términos posibles, se tiene:

−k 1 C A ( s ) +¿

−k 3 C B ( s ) + [ s+ k 4 ] C C ( s ) −k 5 C D ( s )=0
−k 2 C B ( s )−k 4 CC ( s ) +[ s+k 5]C D ( s ) =0

Resolviendo las ecuaciones se encuentra que:

10
C A ( s )=
s +k 1

10 k 1
C B ( s )=
( s+ k 1)(s+ k 2+ k 3)

10 k 1[k 3 s 2+ k 3 k 5 s+ k 4 k 3 s+ k 2 k 5 s+k 2 k 4 k 5+k 3 k 4 k 5]


C C ( s) =
s (s +k 1)( s+ k 4)( s+ k 2+ k 3)(s +k 5+k 4)

10 k 1 [k 2 s +k 2 k 4 +k 3 k 4 ]
C D ( s )=
s(s+ k 4+ k 5)( s+ k 1)(s+ k 2+ k 3)

Aplicando la anti transformada de La Place a las cuatro ecuaciones, se obtienen


las siguientes expresiones que calculan las concentraciones de los componentes
con respecto al tiempo:

C A=10 e−k 1 t

10 k 1 10 k 1
C B= e−k 1 t + e−(k 2+k 3)t
k 2+ k 3−k 1 k 1−k 2−k 3

10 k 5 10[k 3 k 12−k 3 k 4 k 1−k 3 k 5 k 1−k 2 k 5 k 1+ k 2 k 4 k 5+k 3 k 4 k 5]e−k 1 t


C C= +
k 4+ k 5 −(k 4−k 1)( k 2+ k 3−k 1)(k 4+k 5−k 1)

+ 10 k 1 [k 3 k 22−2k 2 k 32 +k 33 −k 2 k 3 k 4−2 k 2 k 3 k 5−k 4 k 32−k 5 k 32−k 5 k 22+ k 2 k 4 k 5+k 3 k 4 k 5]e−(k 2+


(k 2−k 3)( k 1−k 2−k 3)( k 4+ k 5−k 2−k 3)(k 4−k 2−k 3)

+10 k 1[k 3 k 4−k 2 k 5]e−(k 4+k 5)t


−(−k 4−k 5)(k 1−k 4−k 5)(k 2+k 3−k 4−k 5)
10 k 4 10 k 1[k 3 k 4−k 2 k 5]e−(k 4+ k5 )t 10 [k 1 k 2−k 2 k 4−k 3 k 4 ]e−k 1 t
CD= + +
k 4+ k 5 (−k 4−k 5)(k 1−k 4−k 5)(k 2+k 3−k 4−k 5) (k 4 +k 5−k 1)( k 2+ k 3−k 1)

+ 10 k 1 [k 3 k 4−k 2 k 4−k 2 k 3−k 22]e−(k 2+k 3)t


(−k 2−k 3)( k 4+ k 5−k 2−k 3)(k 1−k 2−k 3)

Asignando los mismos valores de las constantes de reacción propuestos:

k 1=0.98 k 2=0.53 k 3=0.50 k 4=0.31 k 5=0.11

Haciendo el cálculo correspondiente, se simplifican las ecuaciones anteriores:

C A=10 e−0.98 t

C B=196 e−0.98t −196 e−1.03 t

C C =2.61904−134.53571 e−0.98 t + 125.31148 e−1.03 t +6.60519 e−0.42 t

C D =7.38095−71.46429 e−0.42 t−6.60519 e−0.98 t + 70.68852e−1.03 t

En la siguiente tabla se muestran los resultados de las concentraciones de los


cuatro componentes involucrados en la reacción, durante el transcurso de 10
minutos.

t (min) Ca Cb Cc Cd Ca + Cb+Cc+Cd
0 10 0 2.13163E-14 0 10
1 3.75311099 3.5876111 1.203284528 1.45599338 10
2 1.40858421 2.6272724 2.491527964 3.47261543 10
3 0.52865729 1.44329977 3.082217638 4.94582531 10
4 0.19841095 0.70492974 3.216363901 5.88029541 10
5 0.07446583 0.32284696 3.152788624 6.44989859 10
6 0.02794785 0.14197405 3.033940947 6.79613715 10
7 0.01048914 0.06071233 2.919736679 7.00906185 10
8 0.00393669 0.02543782 2.828578167 7.14204732 10
9 0.00147748 0.01049381 2.761716191 7.22631252 10
10 0.00055452 0.00427643 2.714843027 7.28032603 10
Obsérvese que se cumple el criterio establecido, que en cada intervalo la suma de
las concentraciones es igual a la inicial. Es decir:

mol
CA =10 =CA +CB +CC +CD
L

Al graficar las concentraciones con respecto al tiempo, se obtiene el siguiente


comportamiento.

10

7
Concentraciones

5 CA
CB
4 CC
CD
3

0
0 1 2 3 4 5 6 7 8 9 10
Tiempo

Con esto podemos comparar ambos gráficos obtenidos y observar que el


comportamiento es el mismo y que ambos métodos son efectivos para la
resolución de problemas de la Ingeniería Química.

También podría gustarte