Trabajo Primer Corte

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

Taller 1 Métodos Numéricos

Métodos Numéricos. 452029: G-01.


Ingeniería Mecánica

𝐿𝑝𝑉 2
8.12 En una sección de tubo, la caída de presión se calcula así: ∆𝑝 = 𝑓 , donde ∆p =
2𝐷
caída de presión (Pa), f = factor de fricción, L = longitud del tubo [m], r = densidad
(kg/m3 ), V = velocidad (m/s), y D = diámetro (m). Para el flujo turbulento, la ecuación
1
de Colebrook proporciona un medio para calcular el factor de fricción, =
√𝑓
ε 2.51 ρVD
−2.0 log (3.7𝐷 + ), donde ε = rugosidad (m), y Re = número de Reynolds, 𝑅𝑒 =
𝑅ε√𝑓 µ
donde µ = viscosidad dinámica (N · s/m2 ).
a) Determine ∆p para un tramo horizontal de tubo liso de 0.2 m de longitud, dadas r =
1.23 kg/m3 , m = 1.79 × 10–5 N · s/m2 , D = 0.005 m, V = 40 m/s, y ε = 0.0015 mm.
Utilice un método numérico para determinar el factor de fricción. Obsérvese que los
tubos lisos tienen Re < 105 , un valor inicial apropiado se obtiene con el uso de la
fórmula de Blasius, f = 0.316/Re0.25.
b) Repita el cálculo pero para un tubo de acero comercial más rugoso (ε = 0.045 mm).

Solución:
Se crea un archivo llamado “caidapresion.m” en donde se escriben todos los comandos,
variables y especificaciones del ejercicio.
function [Dp, f, Re] = caidapresion(L,p,V,D,u,r,ranginf,rangsup)
% f: factor de fricción
% L: longitud del tubo (m).
% p: densidad (kg/m3).
% V: velocidad (m/s)
% D: diámetro(m)
% r: rugosidad(m)
% Re: número de reynolds
% u: viscosidad dinámica(N*m/s2)

Re = (p*V*D)/u;
% otra forma de hallar f es con la "fórmula de blasius"
% f = 0.316/(Re^0.25);
ff = ranginf;
num = 9999;

while ff <= rangsup


% evalúa la función en el rango dado, para hallar el valor de f, con la
% ecuación de colebrook.

Martínez Fernández Andrés Felipe – Negrette Avila Luis Diego


Taller 1 Métodos Numéricos
Métodos Numéricos. 452029: G-01.
Ingeniería Mecánica
ffm = ff/100000;
% busca el valor de f, con un aumento de 0,00001 en cada iteración.
dif = (1/sqrt(ffm))-(-2*log((r/(3.7*D)+(2.51/(Re*sqrt(ffm))))));
if abs(dif) <= abs(num)
num = dif;
pos = ffm;
end
ff = ff + 1;
end
f = pos;
% se halla la caída de presión con los valores encontrados
Dp = f*((L*p*V^2)/(2*D));

Y luego, aunque se puede hacer en el command Windows decidimos en lugar de


utilizarlo, crear otro archivo llamado “run_caidapresion.m” que lo ejecuta al presionar
ctrl + enter y que ya tiene los valores predeterminados para dar el resultado de los
procedimientos anteriores. A continuación, está especificado:
clear
clc
L = 0.2;
p = 1.23;
V = 40;
D = 0.005;
u = 1.79e-5;
r = 0.0015;
ranginf = 2000;
rangsup = 4500;
[Dp, f, Re] = caidapresion(L,p,V,D,u,r,ranginf,rangsup)

El archivo se ejecuta y da como resultado:

Dp =
1.5728e+03

f=
0.0400

Re =
1.3743e+04

Martínez Fernández Andrés Felipe – Negrette Avila Luis Diego

También podría gustarte