Problema 8.

13
En una sección de tubo, la caída de presión se calcula así:
2
Lρ v
Δ P =f .
2D

donde Δp=¿ 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
√f
=−2.0 log
ε
+
(2.51
3.7 D ℜ √ f )
donde ε = rugosidad (m) y Re = número de Reynolds
ρVD
Re =
u

donde m = viscosidad dinámica ( N · s /m2).

a) Determine Δp para un tramo horizontal de tubo liso de 0.2m de longitud, dadas
R=1.23 kg/m3 , m=1.79 ×10 – 5 N · s /m2 , D=0.005 m ,V =40 ms y ε 0.0015 mm. 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
f=
( ℜ) 0.25
b) Repita el cálculo, pero para un tubo de acero comercial más rugoso (e = 0.045 mm).
SOLUCIÓN:

Datos:

Longitud del tubo=0.2 m

kg
Densidad=1.23
m3
m
velocidad=40
s
Diámetro=0.005 m
−6
Rugosidad=1.5 ×10 m
−5 2
Viscosidad dinámica=1.79× 10 N . s/m
El número de Reynolds se obtiene mediante:
ρ× V × D
ℜ=
μ
1.23 × 40 ×0.005
ℜ=
1.79× 10−5
ℜ=13743.01676
El valor inicial ( x 0 ) se obtiene por la fórmula de Blasius:
0.25
f =0.316 /ℜ
 f =0.0291854614

a) La función tomada es:

1
√f
=−2.0 × log ⁡
ε
+
2.51
3.7 × D ℜ× √ f ( )
F (f )=0=2.0 × log
( 3.7ε× D + ℜ×2.51√ f )+ √1f
( )
−6
1.5× 10 m 2.51 1
F (f )=2.0 × log + +
3.7 ×0.005 m 13743.01676 × √ f √ f

La primera derivada es:

−44929 1
F '( f )= −
( )
3 3
44929 3
246000000× ln ( 10 ) × + ×f 2
2f 2
246000000 √ f 37000

Método de Newton-Raphson

F(f)
x i+1=x i− '
F (f )
2−5 −4
ε s=0.5 ×10 =5 ×10
i=0
F ( 0.029185) −0.0249
x 1=0.029185− =0.029185−
F
'
( 0.029185)
−114.115

x 1=0.02893652

i=1
F (0.02893652 ) −0.00361
x 2=0.02893652− =0.02893652−
'
F ( 0.02893652 ) 115.5329

x 2=0.02897208

( 0.02897208−0.02893652
|ε a|(1) = 0.02897208 )× 100=1.22 ×10 −1
<5 ×10−4 , no cumple , seguir

i=2
F ( 0.02897208 ) 0.000492
x 3=0.02897208− =0.02897208−
'
F (0.02897208 ) 115.328

x 3=0.02896722
 |ε a|(2 )=( 0.02896722−0.02897208
0.02896722 )× 100=1.67 × 10 −2
<5 ×10−4 ,no cumple , seguir .

i=3
F( 0.02896722 ) −0.000068
x 4 =0.02896722− =0.02896722−
'
F (0.02896722 ) 115.356

x 4 =0.02896789

|ε a|(3 )=( 0.02896789−0.02896722

0.02896789 )× 100=2.31 ×10 −3 −4
< 5× 10 , no cumple , seguir .

i=4
F (0.02896789 ) 0.0000092
x 4 =0.02896789− =0.02896789−
'
F (0.02896789 ) 115.352

x 4 =0.02896780

|ε a|(4)= ( 0.02896780−0.02896789
0.02896780 )×100=3.17 ×10 −4 −4
<5 ×10 , si cumple , fin.

El error aproximado de la iteración n°4 ya se cumple. Entonces el valor de f es f =0.028967 .

Para hallar la caída de presión en el tubo usamos:

2
Lρ v
∆ p=f
2D
3 2
0.2 m ×1.23 kg / m ×(40 m / s)
∆ p=0.028967 ×
2 ×0.005 m
∆ p=1140.18 Pa
b) La función tomada es:

1
√f
=−2.0 × log ⁡
ε
(+
2.51
3.7 × D ℜ× √ f )
F (f )=0=2.0 × log
( 3.7ε× D + ℜ×2.51√ f )+ √1f
( )
−5
4.5 ×10 m 2.51 1
F (f )=2.0 × log + +
3.7 ×0.005 m 13743.01676 × √ f √ f

La primera derivada es:

−44929 1
F '( f )= −
( )
3 3
44929 9
246000000× ln ( 10 ) × + ×f 2
2f 2
246000000 √ f 37000
MÉTODO DE NEWTON-RAPHSON

F(f)
x i+1=x i− '
F (f )
2−5 −4
ε s=0.5 ×10 =5 ×10
i=0
F ( 0.029185) 0.942073
x 1=0.029185− =0.029185−
F
'
( 0.029185)
−104.827

x 1=0.03857908
i=1
F (0.03857908 ) 0.144516
x 2=0.03857908− =0.03857908−
F
'
( 0.03857908 )
−69.0978

x 2=0.04076923

|ε a|(1) = ( 0.04076923−0.03857908
0.04076923 ) ×100=5.37 × 10 <5 ×10 0 −4
, no cumple , seguir

i=2
F ( 0.04076923 ) −0.00069
x 3=0.04076923− =0.04076923−
'
F (0.04076923 ) −63.627

x 3=0.04075786

|ε a|(2 )= ( 0.4075786−0.04076923
0.4075786 ) ×100=2.78× 10 −2 −4
<5 ×10 , no cumple , seguir .

i=3
F (0.04075786 ) −0.000033
x 4 =0.04075786− =0.04075786−
'
F ( 0.04075786 ) −63.653

x 4 =0.04075840

|ε a|(3 )= ( 0.04075840−0.04075786
0.04075840 ) ×100=1.32× 10 −3
<5 ×10−4 , no cumple , seguir .

i=4
F (0.04075840 ) −0.0000012
x 4 =0.04075840− =0.04075840−
'
F (0.04075840 ) −63.652

x 4 =0.04075838

|ε a|(4)= ( 0.04075838−0.04075840
0.04075838 )× 100=4.90× 10 −5
<5 ×10−4 , si cumple , fin.
El error aproximado de la iteración n°4 ya se cumple. Entonces el valor de f es f =0.040758 .

Para hallar la caída de presión en el tubo usamos:

2
Lρ v
∆ p=f
2D
3 2
0.2m ×1.23 kg / m ×(40 m / s )
∆ p=0.040758 ×
2 ×0.005 m
∆ p=1604.23 Pa

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 x1 +20 x 2+ 25 x 3=4800

30 x 1+50 x 2+ 20 x 3 =5800

18 x 1+30 x 2+ 55 x 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
x 1= = =40.157
∆ 86000
∆2 6135900
x 2= = =71.348
∆ 86000
∆3 4420000
x 3= = =51.395
∆ 86000
En conclusión, de la cantera 1 se debe extraer 40.157 m3 , en la cantera 2 se debe extraer
3 3
71.347 m y por último, de la cantera 3 se debe extraer 51.395 m .
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