Sesion 13 MetodosNumericosIII

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 16

MATEMATICAS AVANZADAS

780027C.
Matemáticas avanzadas

CIRCULO DE MOHR PLANAR


Matemáticas avanzadas

Ejercicio 1
Una descarga de 500 kg de latex se procesan en tanques de mezcla en serie. Los parámetros de los tanques
son:
Parámetro Tanque 1 Tanque 2
Volumen (m3) 80 x 103 150 x 103
Caudal de descarga (l/min) 1 1.5
y el modelo se define como:

𝑑𝐶1 𝑑𝐶2
= −2𝐶1 = 1, 7𝐶1 − 1, 7𝐶2
𝑑𝑡 𝑑𝑡
𝑔
Si las concentraciones iniciales en los tanques corresponden a 𝐶1 = 10 𝑙 y 𝐶2 = 0. Determine el perfil de
concentración (usando Euler, Euler modificado y Euler mejorado), de los tanques durante un día. Nota: Los
valores de 2 y 1,7 son parámetros provenientes de la relación V/Q por lo que sus unidades son (min-1) .
Realice una comparación del comportamiento de los 3 métodos con su respectivo análisis de convergencia de
h.
Matemáticas avanzadas

i=1;
% Condiciones iniciales
𝑑𝐶1 t(i)=0; C1(i)=10; C2(i)=0;
= −2𝐶1 tf=1; % dia
𝑑𝑡 h=1/(24*60); % minutos en un día
𝑑𝐶2 n=(tf-t(i))/h;
= 1, 7𝐶1 − 1, 7𝐶2 for i=1:n
𝑑𝑡 dcdt(i)=-2*C1(i);
C1(i+1)=C1(i)+h*dcdt(i);
dc2dt(i)=1.7*C1(i)-1.7*C2(i);
C2(i+1)=C2(i)+h*dc2dt(i);
𝑦𝑖+1 = 𝑦𝑖 + ℎ𝑓(𝑥𝑖 , 𝑦𝑖 ) t(i+1)=t(i)+h;
end
Respuesta=[t,C1,C2];
Matemáticas avanzadas
Matemáticas avanzadas

Ejercicio 2

Graficar los perfiles de temperatura (usando Euler, RK2 y RK4) de un líquido que fluye en un tubo
cilíndrico de 60 m en estado estacionario, que intercambia calor con una chaqueta tal y como se
presenta en la Figura 1. El modelo de transferencia se presenta a continuación:

𝑑𝑇 𝑈2𝜋𝑅
= 300 − 𝑇 Con las condiciones de contorno: x=0, T=20°C
𝑑𝑥 𝑄𝜌𝑐𝑝

𝑅 =radio=0,2 m
𝑈= Coeficiente global de transferencia de calor=60000 J/hm2°C
𝑄=caudal=4000 L/h
𝜌 =densidad=900 kg/m3
𝑐𝑝 =Calor específico del fluido=300J/Kg°C
Matemáticas avanzadas

i=1;
% Condiciones iniciales
𝑑𝑇 𝑈2𝜋𝑅 x(i)=0; T(i)=20; function dTdt =intercambiador(T)
= 300 − 𝑇 xf=60; % m
𝑑𝑥 𝑄𝜌𝑐𝑝 R=0.2; %m
h=1; % m U= 60000; %J/hm2°C
n=(xf-x(i))/h; Q= 4; %m3/h
𝑅 = 0,2 m
for i=1:n rho=900; %kg/m3
𝑈= 60000 J/hm2°C
K1 = intercambiador(T(i)); cp= 300; %J/Kg°C
𝑄= 4000 L/h = 4 m3 /h
K2 = intercambiador(T(i)+0.5*h*K1); dTdt = U*2*pi*R*(300-T)/(Q*rho*cp);
𝜌 = 900 kg/m3
K3 = intercambiador(T(i)+0.5*h*K2); end
𝑐𝑝 = 300J/Kg°C
K4 = intercambiador(T(i)+h*K3);
T(i+1)=T(i)+h/6*(K1+2*K2+2*K3+K4);
x(i+1)=x(i)+h;
end

Respuesta=[x',T'];
Matemáticas avanzadas
Matemáticas avanzadas
Matemáticas avanzadas

Discretización doble derivada centrada Discretización derivada centrada

3X3 3X1
Matemáticas avanzadas

if i==n
i=1; j=1;
M(i,end-1:end) = [A-B,-2*A-(2*h^2)];
x(i)=0; z(i)=1; xf=2;
R(i,1)=2*x(i)^2*h^2-((A+B)*z(n+1));
h=0.05;
else
n= (xf-x(i)/h;
M(i,[i-1,i,i+1]) = [A-B,-2*A-(2*h^2), A+B];
x(n+1)=2; z(n+1)=2;
R(i,1)=2*x(i)^2*h^2;
for i=2:n
end
x(i)=x(i-1)+h;
end
A= 2*(1+x(i)^2);
j=j+1;
B= h*x(i);
end
if i==2
M(1,:)=[];M(:,1)=[];
M(i,[i,i+1]) = [-2*A-(2*h^2), A+B];
R(1)=[];
R(i,1)=2*x(i)^2*h^2-((A-B)*z(i-1));
Z=inv(M)*R;
else
Z=[z(1);Z;z(n+1)];
Matemáticas avanzadas
Matemáticas avanzadas

TALLER I
• Escriba un algoritmo en MATLAB que resuelva y grafique cualquier ecuación de la forma:

Donde los datos de entrada son C(x), D(x), E(x) y F(x)


Matemáticas avanzadas

j se relaciona con el avance del tiempo.


i se relaciona con la relación espacial.
Matemáticas avanzadas
Matemáticas avanzadas

TALLER II

También podría gustarte