Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Gutarra Rios, Juan Peceros Lecca, Brenda

Descargar como pptx, pdf o txt
Descargar como pptx, pdf o txt
Está en la página 1de 14

GUTARRA RIOS, JUAN

PECEROS LECCA,
BRENDA
Las ecuaciones diferenciales ordinarias que
modelan una realidad especca, aumentan su
complejidad en la medida que se aproximen cada
vez ms al comportamiento del objeto o fenmeno
en estudio, razn por la cual en la mayoria de
los casos hallar su solucin por mtodos
analticos es imposible lo que nos lleva a
utilizar los mtodos numricos.
El trabajo de investigacin tiene
como objetivo; la resolucin de
ecuaciones Diferenciales ordinarias
con Matlab, en el que se desarrollan
los mtodos numricos: De Euler, de
Taylor de orden superior y Runge-
Kutta.
El mtodo de Runge-Kutta es un mtodo genrico de
resolucin numrica de ecuaciones diferenciales.
El mtodo de Runge-Kutta no es slo un nico mtodo,
sino una importante familia de mtodos iterativos,
tanto implcitos como explcitos, para aproximar las
soluciones de ecuaciones diferenciales ordinarias;
estas tcnicas fueron desarrolladas alrededor de 1900
por los matemticos alemanes Carl David Tolm Runge y
Martin Wilhelm Kutta
Mtodo de Runge-Kutta
de Cuarto Orden


+ = + + + + ,

= , , , . . .

Dnde:
= ( , )

= + , +


= + , +


= + , +

ALGORITMO

Este algoritmo calcula la solucin del problema de valor inicial en puntos
equidistantes 1 = 0 + , 2 = 0 + 2, 3 = 0 + 3, , = 0 + , aqu
tiene una solucin nica en [0 , ].

1. Entrada: Valores iniciales 0 , 0 , tamao de paso y nmero de pasos N


2. Para n=0,...,N-1, hacer
1 = ( , )
1 1
2 = + , + 1
2 2
1 1
3 = + , + 2
2 2
1
4 = + , + 3
2
1
+1 = + 1 + 22 + 23 + 4
6
Salida +1 , +1
1. Parar
PROGRAMACIN EN
MATLAB

APLICACIONES

El comportamiento de muchos procesos fisicoqumicos
en ingeniera qumica, particularmente aquellos que
experimenten cambios dependientes del tiempo, se
modelan a travs de un sistema de ecuaciones
deferenciales ordinarias, de aqu resulta la
importancia que adquiere para ingenieros y
cientficos el conocimiento de mtodos de
resolucin de este tipo de ecuaciones. Si bien
muchas ecuaciones diferenciales ordinarias
importantes se pueden resolver a travs de tcnicas
analticas bien conocidas, un gran numero de
ecuaciones diferenciales fsicamente significativas
no pueden resolverse de esta forma.
PROBLEMAS DE
APLICACION

1. Se tiene que enfriar 4536 L. de cido sulfrico (C.E 0.36 kcal(kg.^o C))
y densidad relativa (1.85 KgL) en un CSTR, el cido a 174oC se introduce
en el tanque donde es bien agitado con un serpentn refrigerante de rea 8
m2 y que se mantiene constante a la misma temperatura de 20oC. La capacidad
del tanque es de 4536 L. y el coeficiente de transmisin de calor entre el
serpentn y el cido es 635 kcal((h.m^(2 o) C)) y puede suponerse
constante, suponiendo que el caudal de salida es el mismo que el entrada A
qu temperatura sale el cido en cualquier instante de tiempo?
La ecuacin es:

= + + ( + )

Dnde:

q=4356 V=4356 L. = 635 A=8 m2


(. 2 )

p=(1.85 )
c=0.36 . Te=174oC Ta=20oC
Con la condicin inicial T(0)=174 OC; h=0.5 para hallar T(1)

Solucin:

Reemplazando los valores en la funcin queda:


= 2.681575756 + 207.6315151

Aplicando mtodo de Runge-Kutta de orden 4:

+1 = +

1era iteracin:
1
0.5 = 174 + + 22 + 23 + 4 0.5
6 1
1 = 0,174 = 258.9626664
2 = 0.25,109.2593334 = 85.35566446
3 = 0.25,152.6610839 = 201.7407463
4 = 0.5,73.12962685 = 11.52888069
Reemplazando en la ecuacin

1
0.5 = 174 + + 22 + 23 + 4 (0.5)
6 1
0.5 = 105.5311163

2da iteracin:

1 = 0.5 +
1
1 = 105.5311163 + + 22 + 23 + 4 (0.5)
6 1
1 = 75.35816778
2 = 24.83850894
3 = 58.70658203
4 = 3.354905671

Reemplazando en la ecuacin:
1
1 = 105.5311163 + + 22 + 23 + 4 (0.5)
6 1
1 = 85.60666263
2. Un Balance de masa para un qumico en un reactor completamente mezclado puede
escribirse como:

2
=

Donde:
V=volumen(10 3 ) c=concentracin F=flujo msico(200 )
3 3
Q=flujo volumtrico(1 ) k=flujo de reaccin de 2do orden(0.1 .)
Condicin inicial C(0)=0, hallar C(0.5) para un paso h de 0.5
Solucin:
Reemplazando los valores en la funcin queda:

200 1 (0.1)(10) 2
=
10

Iterando:
0.5 = 0 +
1
0.5 = 0 + 1 + 22 + 23 + 4 (0.5)
6
1 = 20
2 = 17
3 = 17.76875
4 = 11.21835059
Reemplazando:
1
0.5 = 0 + + 22 + 23 + 4 (0.5)
6 1
0.5 = 8.396320883

También podría gustarte