Problema 8.
13
En una sección de tubo, la caída de presión se calcula así:
𝐿𝜌𝑣 2
𝛥𝑃 = 𝑓.
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
de Colebrook proporciona un medio para calcular el factor de fricción.
1 𝜀 2.51
= −2.0 log ( + )
√𝑓 3.7𝐷 Re √𝑓
donde 𝜀 = rugosidad (m) y Re = número de Reynolds
𝜌𝑉𝐷
𝑅𝑒 =
𝑢
donde m = viscosidad dinámica (𝑁 · 𝑠/𝑚2 ).
a) Determine 𝛥𝑝 para un tramo horizontal de tubo liso de 0.2m de longitud, dadas
𝑅 = 1.23 𝑘𝑔/𝑚3 , 𝑚 = 1.79 × 10– 5 𝑁 · 𝑠/𝑚2, 𝐷 = 0.005 𝑚, 𝑉 =
40 𝑚𝑠 𝑦 𝜀0.0015 𝑚𝑚. Utilice un método numérico para determinar el factor de fricción.
Obsérvese que los tubos lisos tienen 𝑅𝑒 < 105 , un valor inicial apropiado se obtiene con el
uso de la fórmula de Blasius, 𝑓 = 0.316⁄(Re)0.25
b) Repita el cálculo, pero para un tubo de acero comercial más rugoso (e = 0.045 mm).
SOLUCIÓN:
Datos:
𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒𝑙 𝑡𝑢𝑏𝑜 = 0.2 𝑚
𝑘𝑔
𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 = 1.23
𝑚3
𝑚
𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 = 40
𝑠
𝐷𝑖á𝑚𝑒𝑡𝑟𝑜 = 0.005 𝑚
𝑅𝑢𝑔𝑜𝑠𝑖𝑑𝑎𝑑 = 1.5 × 10−6 𝑚
𝑉𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑 𝑑𝑖𝑛á𝑚𝑖𝑐𝑎 = 1.79 × 10−5 𝑁. 𝑠/𝑚2
El número de Reynolds se obtiene mediante:
𝜌×𝑉×𝐷
𝑅𝑒 =
𝜇
1.23 × 40 × 0.005
𝑅𝑒 =
1.79 × 10−5
𝑅𝑒 = 13743.01676
El valor inicial (𝑥0 ) se obtiene por la fórmula de Blasius:
𝑓 = 0.316/𝑅𝑒 0.25
𝑓 = 0.0291854614
a) La función tomada es:
1 𝜀 2.51
= −2.0 × log ( + )
√𝑓 3.7 × 𝐷 𝑅𝑒 × √𝑓
𝜀 2.51 1
𝐹(𝑓) = 0 = 2.0 × log ( + )+
3.7 × 𝐷 𝑅𝑒 × √𝑓 √𝑓
1.5 × 10−6 𝑚 2.51 1
𝐹(𝑓) = 2.0 × log ( + )+
3.7 × 0.005 𝑚 13743.01676 × √𝑓 √𝑓
La primera derivada es:
44929 1
𝐹 ′ (𝑓) = − − 3
44929 3 3
246000000 × ln(10) × ( + 37000) × 𝑓 2 2𝑓 2
246000000√𝑓
Método de Newton-Raphson
𝐹(𝑓)
𝑥𝑖+1 = 𝑥𝑖 −
𝐹 ′ (𝑓)
𝜀𝑠 = 0.5 × 102−5 = 5 × 10−4
𝑖=0
𝐹(0.029185) −0.0249
𝑥1 = 0.029185 − ′
= 0.029185 −
𝐹 (0.029185) −114.115
𝑥1 = 0.02893652
𝑖=1
𝐹(0.02893652) −0.00361
𝑥2 = 0.02893652 − ′
= 0.02893652 −
𝐹 (0.02893652) 115.5329
𝑥2 = 0.02897208
0.02897208 − 0.02893652
|𝜀𝑎 |(1) = ( ) × 100 = 1.22 × 10−1 < 5 × 10−4 , 𝑛𝑜 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑠𝑒𝑔𝑢𝑖𝑟
0.02897208
𝑖=2
𝐹(0.02897208) 0.000492
𝑥3 = 0.02897208 − ′
= 0.02897208 −
𝐹 (0.02897208) 115.328
𝑥3 = 0.02896722
0.02896722 − 0.02897208
|𝜀𝑎 |(2) = ( ) × 100 = 1.67 × 10−2 < 5 × 10−4 , 𝑛𝑜 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑠𝑒𝑔𝑢𝑖𝑟.
0.02896722
𝑖=3
𝐹(0.02896722) −0.000068
𝑥4 = 0.02896722 − ′
= 0.02896722 −
𝐹 (0.02896722) 115.356
𝑥4 = 0.02896789
0.02896789 − 0.02896722
|𝜀𝑎 |(3) = ( ) × 100 = 2.31 × 10−3 < 5 × 10−4 , 𝒏𝒐 𝒄𝒖𝒎𝒑𝒍𝒆, 𝒔𝒆𝒈𝒖𝒊𝒓.
0.02896789
𝑖=4
𝐹(0.02896789) 0.0000092
𝑥4 = 0.02896789 − ′
= 0.02896789 −
𝐹 (0.02896789) 115.352
𝑥4 = 0.02896780
0.02896780 − 0.02896789
|𝜀𝑎 |(4) = ( ) × 100 = 3.17 × 10−4 < 5 × 10−4 , 𝒔𝒊 𝒄𝒖𝒎𝒑𝒍𝒆, 𝒇𝒊𝒏.
0.02896780
El error aproximado de la iteración n°4 ya se cumple. Entonces el valor de f es 𝑓 = 0.028967.
Para hallar la caída de presión en el tubo usamos:
𝐿𝜌𝑣 2
∆𝑝 = 𝑓
2𝐷
0.2 𝑚 × 1.23 𝑘𝑔/𝑚3 × (40 𝑚/𝑠)2
∆𝑝 = 0.028967 ×
2 × 0.005 𝑚
∆𝑝 = 1140.18 𝑃𝑎
b) La función tomada es:
1 𝜀 2.51
= −2.0 × log ( + )
√𝑓 3.7 × 𝐷 𝑅𝑒 × √𝑓
𝜀 2.51 1
𝐹(𝑓) = 0 = 2.0 × log ( + )+
3.7 × 𝐷 𝑅𝑒 × √𝑓 √𝑓
4.5 × 10−5 𝑚 2.51 1
𝐹(𝑓) = 2.0 × log ( + )+
3.7 × 0.005 𝑚 13743.01676 × √𝑓 √𝑓
La primera derivada es:
44929 1
𝐹 ′ (𝑓) = − − 3
44929 9 3
246000000 × ln(10) × ( + 37000) × 𝑓 2 2𝑓 2
246000000√𝑓
MÉTODO DE NEWTON-RAPHSON
𝐹(𝑓)
𝑥𝑖+1 = 𝑥𝑖 −
𝐹 ′ (𝑓)
𝜀𝑠 = 0.5 × 102−5 = 5 × 10−4
𝑖=0
𝐹(0.029185) 0.942073
𝑥1 = 0.029185 − = 0.029185 −
𝐹 ′ (0.029185) −104.827
𝑥1 = 0.03857908
𝑖=1
𝐹(0.03857908) 0.144516
𝑥2 = 0.03857908 − ′
= 0.03857908 −
𝐹 (0.03857908) −69.0978
𝑥2 = 0.04076923
0.04076923 − 0.03857908
|𝜀𝑎 |(1) = ( ) × 100 = 5.37 × 100 < 5 × 10−4 , 𝑛𝑜 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑠𝑒𝑔𝑢𝑖𝑟
0.04076923
𝑖=2
𝐹(0.04076923) −0.00069
𝑥3 = 0.04076923 − ′
= 0.04076923 −
𝐹 (0.04076923) −63.627
𝑥3 = 0.04075786
0.4075786 − 0.04076923
|𝜀𝑎 |(2) = ( ) × 100 = 2.78 × 10−2 < 5 × 10−4 , 𝑛𝑜 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑠𝑒𝑔𝑢𝑖𝑟.
0.4075786
𝑖=3
𝐹(0.04075786) −0.000033
𝑥4 = 0.04075786 − ′
= 0.04075786 −
𝐹 (0.04075786) −63.653
𝑥4 = 0.04075840
0.04075840 − 0.04075786
|𝜀𝑎 |(3) = ( ) × 100 = 1.32 × 10−3 < 5 × 10−4 , 𝑛𝑜 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑠𝑒𝑔𝑢𝑖𝑟.
0.04075840
𝑖=4
𝐹(0.04075840) −0.0000012
𝑥4 = 0.04075840 − ′
= 0.04075840 −
𝐹 (0.04075840) −63.652
𝑥4 = 0.04075838
0.04075838 − 0.04075840
|𝜀𝑎 |(4) = ( ) × 100 = 4.90 × 10−5 < 5 × 10−4 , 𝑠𝑖 𝑐𝑢𝑚𝑝𝑙𝑒, 𝑓𝑖𝑛.
0.04075838
El error aproximado de la iteración n°4 ya se cumple. Entonces el valor de f es 𝑓 = 0.040758.
Para hallar la caída de presión en el tubo usamos:
𝐿𝜌𝑣 2
∆𝑝 = 𝑓
2𝐷
0.2 𝑚 × 1.23 𝑘𝑔/𝑚3 × (40 𝑚/𝑠)2
∆𝑝 = 0.040758 ×
2 × 0.005 𝑚
∆𝑝 = 1604.23 𝑃𝑎
function prob_8_13()
resp='s'
while resp=='s' | resp=='S'
L=input('Ingresar la longitud del tubo: ')
d=input('Ingresar la densidad del tubo: ')
v=input('Ingresar la velocidad del flujo: ')
D=input('Ingresar el diámetro del tubo: ')
E=input('Ingresar la rugosidad del tubo: ')
u=input('Ingresar la viscosidad dinámica: ')
disp('El número de reinold es: ')
Re=(d*v*D)/u;
disp(Re);
disp('La funcion es : ');
disp('f(x) = 2.0 * log10(((E)/(3.7 * D))+(2.51/(Re * x^(1/2)))) + (x^(-1/2)))')
deff('y=f(x)','y=2.0*log10(((E)/(3.7*D))+(2.51/(Re*x^(1/2))))+(x^(-1/2))');
disp('El valor inicial tomado es: ');
Vi=0.316/(Re^0.25);
disp(Vi)
n=input('Ingresar número de cifras significativos a considerar: ')
Es=0.5*10^(2-n)
printf('\n\t%s\t\t%s\t\t%s\n','Iteración','Xi','Ea(%)')
printf('\n\t%d\t\t%.2f\n',0,Vi)
i=1
x(i)=Vi-(f(Vi)/numderivative(f,Vi))
i=i+1
x(i)=x(i-1)-(f(x(i-1))/numderivative(f,x(i-1)))
Ea(i-1)=100*(x(i)-x(i-1))/x(i)
printf('\t%d\t\t%.6f\t\t%e\n',i-1,x(i-1),Ea(i-1))
while abs(Ea)>=Es
x(i)=x(i-1)-(f(x(i-1))/numderivative(f,x(i-1)))
i=i+1
x(i)=x(i-1)-(f(x(i-1))/numderivative(f,x(i-1)))
Ea(i-1)=100*(x(i)-x(i-1))/x(i)
printf('\t%d\t\t%.6f\t\t%e\n',i-1,x(i-1),Ea(i-1))
end
resp=input('Si desea repetir el programa presione s/S: ','s')
end
endfunction
Problema 12.13
Un ingeniero civil que trabaja en la construcción requiere 4800, 5800 y 5 690 m3 de arena,
grava fina y grava gruesa, respectivamente, para cierto proyecto constructivo. Hay tres
canteras de las que puede obtenerse dichos materiales. La composición de dichas canteras es
la que sigue:
Arena% Grava fina% Grava gruesa%
Cantera 1 52 30 18
Cantera 2 20 50 30
Cantera 3 25 20 55
¿Cuántos metros cúbicos deben extraerse de cada cantera a fin de satisfacer las necesidades
del ingeniero?
Solución:
Generamos el sistema de ecuación para que satisfacer las necesidades del ingeniero.
52𝑥1 + 20𝑥2 + 25𝑥3 = 4800
30𝑥1 + 50𝑥2 + 20𝑥3 = 5800
18𝑥1 + 30𝑥2 + 55𝑥3 = 5690
Si se trabaja por la regla de Cramer:
Primero hallamos la determinante principal.
52 20 25
∆= (30 50 20) = 86000
18 30 55
Determinante para cada incógnita.
4800 20 25
∆1 = (5800 50 20) = 3453500
5690 30 55
52 4800 25
∆2 = (30 5800 20) = 6135900
18 5690 55
52 20 4800
∆3 = (30 50 5800) = 4420000
18 30 5690
Ahora hallamos el valor de cada incógnita:
∆1 3453500
𝑥1 = = = 40.157
∆ 86000
∆2 6135900
𝑥2 = = = 71.348
∆ 86000
∆3 4420000
𝑥3 = = = 51.395
∆ 86000
En conclusión, de la cantera 1 se debe extraer 40.157 𝑚3 , en la cantera 2 se debe extraer
71.347 𝑚3 y por último, de la cantera 3 se debe extraer 51.395 𝑚3 .
function cramer()
printf('Válido para un sistema de 3 ecuaciones con 3 variables\n')
n=input('Ingresar al cantidad de ecuaciones del sistema : ')
for i=1:n
for j=1:n
princ(i,j)=input('Ingrese los coeficientes de cada ecuación : ')
end
end
for i=1:n
indep(i)=input('Ingrese los términos independientes : ')
end
DP=det(princ)
A1=princ
A1(:,1)=indep
D1=det(A1)
A2=princ
A2(:,2)=indep
D2=det(A2)
A3=princ
A3(:,3)=indep
D3=det(A3)
X1=D1/DP
X2=D2/DP
X3=D3/DP
printf('Las soluciones son x1 = %.3f, x2 = %.3f y x3 = %.3f ',X1,X2,X3)
endfunction
