Proyecto Integrador: Hidraúlica del Norte

Salma Daniela Aguirre Isais A01233498

Araceli Quezada Rodríguez A01234025
Valeria Azpilcueta Zuñiga A01235419
Samantha Atziri Chávez Urosa A01365175
José Emilio Venegas Medina A01027441

ITESM
Métodos Numéricos en Ingeniería
Prof. Alejandro Martínez Bórquez
Parcial 1
13/ Marzo/2020
Resumen
En los problemas clásicos de flujo en tuberías se presentan dos situaciones típicas en
las que el analista busca determinar el caudal (​Q) o el diámetro de diseño (​F​); estos
han sido tradicionalmente resueltos mediante procesos iterativos sobre el factor de
fricción (​D​) tratando de hallar una variable que no es propiamente la incógnita principal
del problema, con el ánimo aparente de facilitar la solución y sacar del paso los
cálculos en derivadas que implica una solución más real. En este proyecto, se presenta
una propuesta de una nueva formulación para este problema basada en una
fundamentación numérica que encuentra directamente las variables de interés,
adaptando matemática y computacionalmente los métodos de Newton Raphson y
Bisección el cual converge rápidamente a la solución deseada. Se desarrolló una
aplicación educativa en Matlab, mediante ecuaciones. Los resultados encontrados
llaman la atención de la forma como tradicionalmente se enseña a resolver los
problemas de flujo en tuberías en los textos clásicos de la mecánica de fluidos.

Introducción
En el diseño o análisis de sistemas de tuberías existen seis variables involucradas,
excluyendo la aceleración de la gravedad (g), que es considerada constante para
cualquier sistema de tuberías.
1. Pérdidas o energía del sistema de tuberías (hf).
2. Caudal (Q) o Velocidad del Flujo (V).
3. Diámetro de la tubería (D).
4. Longitud de la tubería (L).
5. Rugosidad absoluta de la pared de la tubería (e).
6. Propiedades del fluido (viscosidad cinemática (n)).
En la práctica es de interés para el analista determinar una de las tres primeras
variables, mientras que los demás se conocen o pueden ser determinadas por el
diseñador. En función de la incógnita a determinar, se definen tres clases de
problemas, estos se resumen en la siguiente tabla:
Por requerir de una explicación ampliada se describe a continuación la formulación
matemática empleada para la solución de los problemas clase II y clase III.

Marco teórico
Modelos matemático y numérico
El modelo matemático se funda en las fórmulas de Colebrook y Darcy-Weisbach
descritas en las ecuaciones (1) y (2) respectivamente.

En la ecuación (1) la variable ​Re representa el número adimensional característico de

los problemas de flujo en tuberías denominado número de Reynolds sobre el cual se
precisan detalles más adelante. La ecuación (3) que se obtiene de la ecuación (1)
permite calcular el coeficiente de fricción cuando se desconoce el caudal (Q); para esta
situación se asume turbulencia completa que conduce a un valor despreciable del
segundo término del paréntesis de (1) donde interviene el número de Reynolds ​Re​.
En el problema tipo II el método requiere de un valor semilla para inicializar el proceso
iterativo, por lo cual, utilizando la ecuación (2) y reemplazando las expresiones (4) y (5)
en dicha ecuación, se obtiene el caudal semilla descrito por la ecuación (6).

A continuación, se evalúa a partir del caudal semilla Q0 el número de Reynolds ​Re

descrito en la ecuación (7).

Si el número de Reynolds es menor de 2000 el flujo se considera laminar y es posible

calcular la fricción a partir de la ecuación (8) propuesta por Poiseuille en 1846.

Si el flujo es ​turbulento (​Re > 4000) es posible calcular el factor de fricción ​f a partir de
la ecuación explícita (9) con un error alrededor de 1% con respecto a la ecuación de
Colebrook (3).
La ecuación (9) tiene validez dentro de ciertos rangos de e/D y ​Re definidos en las
expresiones (10) y (11).

La ecuación (12) obtenida de una acomodar convenientemente la ecuación (6) es la

función utilizada para resolver el problema de caudal; esta ecuación proviene de
reemplazar las ecuaciones (4) y (5) en la ecuación (2) de Darcy-Weisbch e igualar a
cero.

Conocida la función objetivo se procede a solucionarla. Se sabe que el método de

Newton Raphson requiere de la función y su derivada, esto significa que la solución del
problema puede ser determinada iterando sobre la ecuación (13).

La derivada Z2’ de la función descrita (13) respecto a Q es altamente compleja por ser
implícita y compuesta. La expresión se obtuvo con la ayuda del CAS libre Máxima. El
número de iteraciones de cálculo es aquel que lleva la función a un valor menor o igual
a 1x10−14 .
Al igual que en el cálculo de caudal, el problema tipo III requiere de un valor semilla del
diámetro (D) para iniciar el proceso iterativo. Para establecer este valor se parte de un
diámetro arbitrario con el cual se calcula el factor de fricción ​f utilizando la ecuación (8)
o (9) según sea el caso; se despeja el diámetro ​D de la ecuación (2) previo reemplazo
de las ecuaciones (4) y (5) y se obtiene el diámetro corregido descrito en la ecuación
(14).
La ecuación (14) tienen la ventaja de que el valor inicial del diámetro supuesto
inicialmente no tiene que ser cercano al valor real ya que el valor del diámetro
corregido siempre será una muy buena aproximación a la solución, por tanto, se
requieren pocas iteraciones para obtener resultados satisfactorios. En las iteraciones
sucesivas el cálculo del factor de fricción ​f se realiza usando las ecuaciones (8) o (9)
según el caso.
Definido el valor inicial del diámetro o diámetro semilla, se plantea la solución iterativa
por medio del método Newton Raphson que requiere de la función y su derivada. La
ecuación (15) es la función utilizada. Nótese que siendo Q una constante la función
ahora es el diámetro.

La derivada de la función respecto de ​D además de implícita es compuesta por lo cual

se recurrió a la ayuda del CAS libre Máxima. La solución del problema puede ser
determinada iterando sobre la ecuación (16).

El número de iteraciones del proceso de cálculo es aquel que lleva la función a un valor
menor o igual a 1x10−14 .

Desarrollo de la solución
Primero se definen las variables y se le piden los valores al usuario. Se tienen que
definir como globales las variables D (diámetro), E (rugosidad) y Re (Reynolds) porque
se utilizan dentro de la función para calcular la ecuación de Colebrook-White.
global​ D E Re;

P=input(​'Densidad del fluido (kg/m^3): '​);

V=input(​'Viscosidad dinámica (kg/m-s): '​);
D=input(​'Diámetro interno de la tubería (mm): '​) / 1000;
E=input(​'Rugosidad De la tubería (mm): '​) / 1000;
F=input(​'Flujo volumétrico (m^3/s): '​);
Una vez definidas las variables, se puede calcular el número de Reynolds y el régimen
del fluido:
Re=P*F*D/V; ​% Número de Reynolds
regimen = ​"laminar"​;

% Regimen
if​ (Re > 2000)
​if​ (Re > 4000)
regimen = ​"turbulento"​;
​else
regimen = ​"transición"​;
​end
end

A continuación, se calcula el factor de fricción en base al régimen del fluido. Si el

régimen es de transición entre laminar y turbulento, no se calcula la fricción. Se intentó
desarrollar el código con el método de Newton, pero en muchas ocasiones la función
brincaba a valores cada vez más grandes. Esto tiene sentido cuando se toma en
cuenta que la ecuación es logarítmica y tiene la siguiente forma:

Por lo cual el método de Newton y el de la secante no son buenos candidatos para

resolver esta ecuación. Dado esto, se optó por utilizar el método de la bisección.
Además, se sabe que el factor de fricción debe ser un valor entre 0 y 1, así que
podemos utilizar esos valores como iniciales.
f = 0;

if​ (regimen == ​"laminar"​)

f = 64/Re;

elseif​ (regimen == ​"turbulento"​)

​% Ya que la ecuación es logaritmica, se optó po utilizar el

​% método de la bisección

​ Valores iniciales (La fricción suele estar entre 0 y 1)

%
x1 = 0.001;
x2 = 1;

​ Error inicial (100%)

%
e = 100;

​while​ (e > .1)

f = (x2 + x1)/2;
fr = c_w(f);

​if​ (fr > 0)

e = abs((f - x2) / f) * 100;
x1 = f;
​elseif​ (fr < 0)
e = abs((f - x1) / f) * 100;
x2 = f;
​end

​end

end

Por último, se muestran los resultados al usuario:

fprintf(​"El número de Reynolds es %f\n"​, Re);
fprintf(​"El régimen de flujo es %s\n"​, regimen);
if​ (regimen ~= ​"transición"​)
fprintf(​"El factor de fricción es %f\n"​, f);
end

En el código anterior, la función “c_w()” es la que calcula la ecuación de

Colebrook-White. Sin embargo, la ecuación que se resuelve no es la original de
Colebrook, sino la versión igualada a 0. La función se ve así:
% Ecuación de Colebrook-White
function​ fr = c_w(f)
​global​ D E Re;

​% Para prevenir errores

​if​ (f == 0)
f = 0.01;
​end

​ Ecuación original (despejada para igualar a 0)

%
syms ​fun(x)
fun(x) = 1/sqrt(x) + .869*log((E/D)/3.71 + 2.51/(Re*sqrt(x)));

fr = fun(f);
end

Resultado
Para comprobar el funcionamiento correcto del programa, se propusieron 4 situaciones
diferentes. Se corrió el programa con esos valores y se comparó el resultado con el de
la calculadora de factor de fricción LMNO Engineering (plataforma enfocada en la
ingeniería de fluidos). Sin embargo, esta calculadora requiere que se le de el número
de Reynolds, así que será calculado en cada caso.

Caso 1
Densidad del fluido (P): 700 kg/m^3
Viscosidad dinámica (V): 0.007 kg/m-s
Diámetro interno de la tubería (D): 1000 mm
Rugosidad De la tubería (E): 5 mm
Flujo volumétrico (F): 10 m^3/s
P *F *(D/1000)
Reynolds = V
= 1000000

Resultado de la calculadora: ​0.030511296

Resultado del programa: ​0.030416

Error: ​0.328%

Caso 2
Densidad del fluido (P): 250 kg/m^3
Viscosidad dinámica (V): 0.6 kg/m-s
Diámetro interno de la tubería (D): 50 mm
Rugosidad De la tubería (E): 0.3 mm
Flujo volumétrico (F): 5 m^3/s
P *F *(D/1000)
Reynolds = V
= 104.166667

Resultado de la calculadora: ​0.6144

Resultado del programa: ​0.6144

Error: ​0%

Caso 3
Densidad del fluido (P): 420 kg/m^3
Viscosidad dinámica (V): 0.04 kg/m-s
Diámetro interno de la tubería (D): 50 mm
Rugosidad De la tubería (E): .3 mm
Flujo volumétrico (F): 5 m^3/s
P *F *(D/1000)
Reynolds = V
= 2625

Resultado de la calculadora: ​N/A (régimen de transición)

Resultado del programa: ​N/A (régimen de transición)

Error: ​N/A

Caso 4
Densidad del fluido (P): 690 kg/m^3
Viscosidad dinámica (V): 0.3 kg/m-s
Diámetro interno de la tubería (D): 5000 mm
Rugosidad De la tubería (E): 7 mm
Flujo volumétrico (F): 15 m^3/s
P *F *(D/1000)
Reynolds = V
= 172500

Resultado de la calculadora: ​0.022806771

Resultado del programa: ​0.022604

Error: ​0.89%

Conclusiones
Con base a los resultados obtenidos la ecuación Colebrook-white ayuda a realizar los
cálculos necesarios para llegar a una determinación de la caídas de presión de los
diferentes fluidos pasando por las distintas tuberías analizadas, de igual manera esta
ecuación te ayuda a tener una solución rápida y económica en cuanto a situaciones
reales siempre manteniendo una precisión en cuanto a los resultados. Por otro lado con
el ejercicio resuelto se observa que el número de Reynolds al disminuir hace que de
igual manera el factor de fricción disminuye y viceversa, dicho número también
depende de la rugosidad relativa, con esto queremos decir que el factor de fricción
depende tanto de las características del tubo como de la tubería. En este trabajo
pudimos observar que el número de Reynolds contra el de fricción es muy similar al
caso de Moody solo variando un poco en el factor de fricción.

Referencias
Anaya-Durand, A. Cauich-Segovia, G. Funabazama-Bárcenas, O. Alfonso, V. (2014).
Evaluación de ecuaciones de factor de fricción explícito para tuberías.​
ScienceDirect. Recuperado de
https://www.sciencedirect.com/science/article/pii/S0187893X1470535X
Cajal, A. (2019). ​Número de Reynolds: para qué sirve, cómo se calcula, ejercicios​.
Lifeder.com. Recuperado de
https://www.lifeder.com/numero-de-reynolds/#Para_que_sirve_el_numero_de_R
eynolds
Jiménez, S. (2016). ​Factor de fricción de la ecuación de Colebrook-White, con métodos
numéricos.​ Hidráulica Fácil. Recuperado de
 https://www.hidraulicafacil.com/2016/04/FactorFriccionColebrook-White.html?m=
1
LMNO Engineering (2014). ​Moody Friction Factor Calculator.​ Recuperado de
https://www.lmnoeng.com/moody.php
Pallares, M., Pisca, R., Rodríguez, W., (2015). ​Análisis de flujo en tuberías en serie
bajo una nueva visión numérica: Desarrollo de una aplicación educativa.​ [Archivo
PDF]. Recuperado de:
file:///C:/Users/Usuario/Downloads/Dialnet-AnalisisDeFlujoEnTuberiasEnSerieBaj
oUnaNuevaVision-5432282.pdf

Liga de acceso a los videos de las sesiones

● Parte 1: ​https://youtu.be/ZPMBd-mzBK0
● Parte 2: ​https://youtu.be/j7yTRmT0XTw