Scribd
Cargar
143 vistas4 páginas

Taller 3-Análisis

El documento presenta 4 problemas resueltos relacionados con cálculo numérico. El primero usa el método de trapecios para aproximar el área bajo una curva. El segundo aplica la regla de Simpson 1/3 para calcular un promedio ponderado. El tercero resuelve una ecuación diferencial con el método de Runge-Kutta. Y el cuarto determina el tiempo para vaciar un tanque cilíndrico usando también Runge-Kutta.

Cargado por

Christian Hernández
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 PDF, TXT o lee en línea desde Scribd

Temas abordados

  • tanque cilíndrico,
  • análisis de rendimiento,
  • métodos de extrapolación,
  • análisis de variables,
  • intervalos de tiempo,
  • concentración de salida,
  • métodos de validación,
  • reactor químico,
  • promedio ponderado,
  • valores de función
143 vistas4 páginas

Taller 3-Análisis

El documento presenta 4 problemas resueltos relacionados con cálculo numérico. El primero usa el método de trapecios para aproximar el área bajo una curva. El segundo aplica la regla de Simpson 1/3 para calcular un promedio ponderado. El tercero resuelve una ecuación diferencial con el método de Runge-Kutta. Y el cuarto determina el tiempo para vaciar un tanque cilíndrico usando también Runge-Kutta.

Cargado por

Christian Hernández
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 PDF, TXT o lee en línea desde Scribd

Temas abordados

  • tanque cilíndrico,
  • análisis de rendimiento,
  • métodos de extrapolación,
  • análisis de variables,
  • intervalos de tiempo,
  • concentración de salida,
  • métodos de validación,
  • reactor químico,
  • promedio ponderado,
  • valores de función

ANLISIS NUMRICO

TALLER TERCER CORTE (Respuestas en rojo)


PROFESOR: GIOVANI RODRGUEZ
JOSE MIGUEL BELTRN CELIS
CHRISTIAN FERNEY HERNNDEZ TIQUE
UNIDADES TECNOLGICAS DE SANTANDER
BUCARAMANGA
1. Encontrar una aproximacin al rea bajo la curva de la siguiente funcin usando
el mtodo de trapecios. Tomar 6 sub-intervalos iguales:
() = {

( + );
+ ; <

Una aproximacin con GEOGEBRA del resultado de la Integral:

I=8.59+10 = 18.59; rea bajo la curva de f(x)


Con el mtodo de TRAPECIOS
-

Sub-intervalos

2 (-) = 3
-

3
6

Saltos entre un valor y otro

Declaracin del eje x & evaluacin de f(x) en x

f(x)

3
0

2
2
2
2
1.0078 1.5942 1.8542 1.9965 3.2533 3.7724

El resultado de la integral con el mtodo de trapecios es: 18.209

2. La concentracin de salida de un reactor se mide en distintos momentos durante


un periodo de 12 horas:
t
c

0
2
2.1 4

4
5

6
8
5.5 5

10 12
3 1.2

El caudal de salida en m3/s se puede calcular con la siguiente ecuacin:


() = 20 + 10[

( 10) ]
12

Use la regla de Simpson 1/3 para determinar el promedio ponderado () de concentracin


de salida del reactor durante el periodo de 12 horas, donde:

0 ()()

0 ()

Evaluando Q(t) en el tiempo dado, tenemos:


t=[0 2 4 6 8 10 12]
q=20+10*sin((pi/12)*(t-10))
q=[15.0000 11.3397 10.0000 11.3397 15.0000 20.0000 25.0000]
Multiplicando el vector q*c:
c=[2.1 4 5 5.5 5 3 1.2]
q=[15.0 11.3397 10.0 11.3397 15.0 20.0 25.0]
f=c*q
f= [31.5 45.359 50.0 62.3686 75.0 60.0 30.0]

Aplicando Simpson 1/3 para 0 ()()


Resultado de la integral= 654.9402

Aplicando Simpson 1/3 para 0 ()


Resultado de la integral= 173.812

Resolviendo

0 ()()

0 ()

654.9402
173.812

= .

3. Usar el mtodo de Runge Kutta para resolver la siguiente ecuacin diferencial


dy/dt = ysen3 (t) en el intervalo [0,1] con h=0,1 y condicin inicial y(0)=1.
De acuerdo al programa Runge Kutta de MATLAB declaramos las siguientes variables:
R=[0 1] % Rango de solucin
h= 0.1 % Tamaos de los sub-intervalos o pasos en el mtodo
c= 1 % Valor de la condicin inicial
fun= (y*sin(t))^3 %Definicin de la ecuacin diferencial a solucionar
Ejecutamos el programa
MetodoRK4(R,h,c,fun) Para lo cual obtenemos la tabla con sus respectivas K y las
soluciones en y:

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

x
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1

K1
0
0
0.00099504
0.0078445
0.025859
0.059413
0.11179
0.18527
0.28134
0.40097
0.54471

K2
0
0.00012484
0.0033375
0.015155
0.040449
0.083039
0.14566
0.23016
0.33773
0.46909
0.62458

K3
K4
0
0
0.00012484 0.00099502
0.0033379 0.0078442
0.015161
0.025858
0.040478
0.059409
0.083136
0.11178
0.14591
0.18525
0.23066
0.28132
0.33862
0.40095
0.47054
0.5447
0.62673
0.71258

y
1
1
1.0004
1.002
1.0061
1.0145
1.0291
1.0523
1.0862
1.1333
1.1959

4. Si se drena el agua desde un tanque cilndrico vertical por medio de una vlvula
en la base, el lquido fluir rpido cuando el tanque est lleno y despacio
conforme se drene. La tasa a la que el nivel del agua disminuye es:

donde k es una constante que depende de la forma del agujero y del rea de la seccin
transversal del tanque. La profundidad del agua y se mide en metros y el tiempo t

en minutos. Si k=0.06, determine cunto tiempo se requiere para vaciar el tanque si el


nivel de fluido se encuentra en un inicio a 3m. Utilice el mtodo RK4 con un
incremento de 0.5 segs.
De acuerdo al siguiente programa de Runge Kutta en MATLAB:
F=inline('-0.06*sqrt(y)');
h=0.5;
y(1)=3;
t(1)=0;
i=1;
while y(i)>0
k1=F(y(i));
y1(i+1)=y(i)+F(y(i))*h;
if y1(i+1)<0
break
end
t(i+1)=t(i)+h;
k2=F(y1(i+1));
yn(i+1)=y(i)+(h/2)*(k1+k2);
fprintf('%5.5f %5.5f %5.5f %5.5f %5.5f %5.5f
%5.5f\n',t(i),y(i),k1,t(i+1),y1(i+1),k2,yn(i+1))
y(i+1)=yn(i+1);
i=i+1;
end
grid on
plot(t,y)
grid on

t
56.50000

k1
0.00151

Grfica Vaciado del tanque


k2
t +0.5
k1
-0.00233 57.00000 0.00034

k2
-0.00111

y
0.00065

Para los 57 segundos el tanque se ha vaciado porque su nivel es 0.00065 metros

También podría gustarte