Caso de Estudio - Método de Von Karman - 4BV1 - E#1 PDF

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

INSTITUTO POLITÉCNICO

NACIONAL
Unidad Profesional Interdisciplinaria de
Ingeniería

Campus Guanajuato

“La técnica al servicio de la patria”

UA: Fenómenos de Transporte| 4BV1

Ingeniería biotecnológica

CASO DE ESTUDIO EQUIPO #1: MÉTODO DE VON KARMAN

Equipo:
Cortes Gutiérrez Vanessa del Rocío
García Cervantes Fátima Jocelyn
Memije Rodríguez Fernanda Regina
Méndez Ruiz Nayeli
Sánchez González Jesús
Vega Hernández Isaí

Docente: José Alfredo Hernández Maldonado.

Marzo 23, 2021.

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 1
Grupo 4BV1
Contenido
Resumen..................................................................................................................... 3
Objetivos ..................................................................................................................... 3
Introducción ................................................................................................................ 3
Descripción y planteamiento del problema................................................................ 4
Cálculos y resultados ................................................................................................. 5
Discusión .................................................................................................................. 17
Conclusiones ............................................................................................................ 21
Bibliografía ................................................................................................................ 22

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 2
Grupo 4BV1
Resumen
En este documento se detalla paso a paso el análisis y resolución de un ejercicio
que plantea como incógnita el caudal volumétrico del agua a una temperatura de
68°F, considerando que el sistema es un tanque de carga constante y que el
diámetro interno de la tubería es de 5 in, además de tomar en cuenta los distintos
largos. La resolución de dicho problema se lleva a cabo aplicando los conocimientos
obtenidos en el primer semestre de la materia fenómenos de transporte, como lo
son: el uso de la ecuación de Bernoulli, el planteamiento de ecuaciones a tomar y
el cálculo del número de Reynolds.

Objetivos
• Identificar cada una de las operaciones básicas el mecanismo y
características de los procesos de transporte, así como establecer las
ecuaciones de conservación y de velocidad de transporte o flujos,
representativas de los mismos, a fin de determinar las propiedades de
transportado con vistas al diseño.
• Adquirir las capacidades necesarias para aplicar correctamente las leyes
básicas de conservación de la cantidad de movimiento, masa y energía a un
volumen de control, formulando las correspondientes ecuaciones de
variación de estas tres propiedades en el citado volumen.

Introducción
En el presente trabajo relacionado al estudio de un problema se abordarán algunos
de los temas vistos en el primer departamental de la materia de fenómenos de
transporte, empezando por flujo volumétrico el cual, según Cengel (2009) está
definido como el volumen de un fluido que fluye por una sección transversal por
unidad de tiempo. Después tenemos la ecuación de Bernoulli y como menciona
Calderón (2017) esta es consecuencia del principio de conservación de la energía
mecánica en un flujo ideal, el cual se supone incompresible, irrotacional, no viscoso
y estacionario. Es interesante mencionar que esta ecuación, como se conoce hoy
en día, fue derivada por Leonhard Euler, quien determinó matemáticamente que la
presión es un campo escalar variable en el espacio y el tiempo.
También, cabe mencionar que según Felder (2004) el número de Reynolds es un
grupo adimensional para un líquido que flu-ye por una tubería el cual se define con
la fórmula:
𝜌𝐷
𝑅𝑒 = 〈𝓋〉 donde 𝐷 es el diámetro de la tubería, 〈𝓋〉 es la velocidad del fluido, 𝜌
𝜇
es su densidad y 𝜇 su viscosidad. Cabe mencionar que cuando el valor del número
de Reynolds es aproximadamente menor a 2100, el flujo es laminar (es decir, el

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 3
Grupo 4BV1
líquido se mueve en líneas de flujo suave) y para números de Reynolds mayores de
2100, el flujo es turbulento, caracterizado por considerable agitación.
Es importante saber que existe el factor de fricción de Darcy, el cual es una variable
adimensional y depende tanto del número de Reynolds, como de la rugosidad
relativa de la tubería, la cual es un indicador de las imperfecciones del material de
la misma tubería (Durán, 2009).
La altura en las tuberías está compuesta por la altura geométrica total, más la altura
equivalente a todas las pérdidas por fricción. La altura geométrica total es igual a la
diferencia de altura geométrica entre el nivel de descarga del líquido y el nivel de la
superficie del mismo en la fuente de abastecimiento. La altura equivalente a las
pérdidas por fricción puede ser calculada dependiendo de si es entrada borda o
entrada ordinaria, por lo cual a continuación se realizarán los cálculos
correspondientes en cada caso.
Según Castejón (2014) el método de Von Karman tiene como objetivo calcular los
valores de espesores y esfuerzos dentro de la capa límite, partiendo de un
perfil de velocidades aproximado del se obtienen como resultado los valores
de espesor, coeficiente local de fricción, coeficiente de fricción medio y
finalmente la fuerza de fricción sobre una superficie como consecuencia de moverse
a través de un fluido viscoso incompresible. Cabe mencionar que esta ecuación
tiene solución tanto para la zona de comportamiento laminar como para la
turbulenta.

Descripción y planteamiento del problema.


Se tiene una tubería de tubo liso por la cual fluye agua, considerándose ésta un
fluido newtoniano en estado estacionario de un recipiente a otro, no hay ninguna
bomba de por medio, por lo tanto no hay trabajo. Se busca determinar la cantidad
de flujo volumétrico que pasa por esta tubería del recipiente inicial al final con una
diferencia de altura de 38 ft a una temperatura constante de 68°F. Para poder
resolver este problema es necesario tener conocimientos acerca de Fenómenos de
transporte, de manera más específica de los temas: Balance macroscópico en un
sistema isotérmico, ecuación de Bernoulli, número de Reynolds, factor de fricción,
y conceptos básicos como lo es el flujo volumétrico, densidad, gravedad y velocidad.
Es importante identificar los datos que nos proporciona para poder recurrir al método
más propicio y darle una solución.
Para este caso iniciamos con la ecuación de Bernoulli, en la cual solo tenemos como
incógnita la velocidad, por lo que se despeja, y se sustituirá en la fórmula de
Reynolds (Re), mediante un proceso se iterativo se obtendrá el Re así como el de
fricción para posteriormente, poder calcular nuestros flujo volumétrico. A
continuación, se verá de manera más detallada este proceso.

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 4
Grupo 4BV1
Cálculos y resultados
Problema: Para el sistema mostrado en la siguiente figura, calcule el flujo
volumétrico del agua a 68°F

Figura 1 Dibujo del proceso del problema.

Tomando la ecuación de Bernoulli:

Consideraciones:


• (Solo existe la presión atmosférica)
• No cambia de velocidad
• Tubo liso
• Fluido Newtoniano (agua)

Una vez tomadas estas consideraciones, la ecuación se simplifica de la siguiente


forma:

Siendo:

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 5
Grupo 4BV1

Como no se conoce la velocidad del fluido, se despeja de la ecuación anterior
quedando de la siguiente manera:

Como no se conoce el factor de fricción se procede a calcular el número de


Reynolds:

Siendo:
• ρ = 62.4 lbm/ft^3
• D = 5 in = 0.417 ft
• μ = 1cP (para el agua) = 6.72*10^-4 lbm/ft*s
• g = 32.174 ft/s^2
• z1 = 0 ft
• z2 = -38 ft (es negativo porque el punto de referencia se toma desde el borde
superior del tanque de arriba)
• L = 49 ft
Para determinar las pérdidas por componentes se utiliza la siguiente tabla:

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 6
Grupo 4BV1
Figura 2 Tabla de valores para accesorios para tubería de flujo turbulento. Tomado de Levenspiel (2014).

P1) Usando una entrada Borda: (L/D)eq = 1 entrada Borda + 2 codos de 90° = 30
+ 2(30) = 90
• Como no se puede conocer la velocidad, el número de Reynolds ni el factor
de fricción, se hace un proceso iterativo para conocer estas dos últimas
incógnitas. Para ello se hace uso de la siguiente tabla:

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 7
Grupo 4BV1
Figura 3 Gráfica para estimar el factor de fricción. Tomado de Bird (2002).

1. Suponiendo para la primera iteración un Re = 1.1*10^6, se obtiene un


ϝ=0.0028. Para confirmar que el valor es correcto se debe cumplir la siguiente
igualdad:

• Con esto se calcula el nuevo número de Reynolds y se compara con el que


se supuso.

• Después se calcula el factor de fricción con ayuda de la tabla. Con ello se


obtiene un factor de fricción de 0.0028, el cual, a comparación del factor
anterior, se tiene que es el mismo

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 8
Grupo 4BV1
• Para tener más confianza en este valor se vuelve a calcular el número de
Reynolds con el factor de fricción de 0.0028. Obteniendo lo siguiente:

• Como los valores de Reynolds coinciden, el valor del factor de fricción sería
el mismo.
• Por lo tanto, como los valores del factor de fricción y el número de fricción
coinciden, se les considera como correctos.

• Calculando el flujo volumétrico, siguiendo la siguiente fórmula:

• Por lo que el flujo volumétrico es igual a 4.40 ft^3/s usando una entrada
borda.
P2) Usando una entrada ordinaria: Si se cambia la entrada Borda por una entrada
ordinaria, el cálculo de (L/D)eq y el cálculo del número de Reynolds sería el
siguiente (los demás valores son los mismos):
• (L/D)eq = 1 entrada ordinaria + 2 codos de 90° = 16 + 2(30) = 76

• Siguiendo el mismo proceso iterativo anteriormente hecho para buscar los


valores de Reynolds y el factor de fricción, al final de este proceso se tendría

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 9
Grupo 4BV1
que los valores donde el factor de fricción y Re coinciden serían los
siguientes:


Por último, se calcula el flujo volumétrico:

• Dando un flujo volumétrico de 4.58 ft^3/s usando una entrada ordinaria.


Utilizando el método de Von Karman;
Ecuación de Von Karman para flujos turbulentos en tubos lisos, se compone de la
siguiente manera:

Siendo:
• Re = Número de Reynolds
• f(D)= factor de fricción de Darcy
Despejando en número de Reynolds (Re) quedaría como:

La ecuación de Von Karman utiliza un factor de fricción de Darcy. Este factor está
relacionado con el factor de fricción de Fanning, descrito en la siguiente ecuación:

Siendo:
• f(F) = factor de fricción de Fanning
• f(D) = factor de fricción de Darcy
Como el factor de fricción de Darcy en la ecuación de Von Karman, se encuentra de
modo implícito, no es posible poder hacer un proceso iterativo como se realizó
anteriormente, por lo que se opta por realizar una tabla tomando diversos valores
para el factor de fricción de Darcy, más precisamente a valores cercanos a el factor

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 10
Grupo 4BV1
de fricción de Fanning que se obtuvo anteriormente que quedan entre 0.00275 y
0.0028 (para el factor de fricción de Darcy sería un intervalo de 0.011 a 0.112).
Haciendo una tabla tomando un intervalo de 0.0105 a 0.012 para el f(D), con su
respectiva conversión al f(F), se calculan los Re con la ecuación de Von Karman y
también los flujos volumétricos, esto con el fin de comparar los resultados obtenidos
de flujo volumétrico con el procedimiento anterior (véase P1 y P2 de esta sección
de resultados)
Tipo de entrada usada
f(F) f(D) Re
(tanque superior)
0.002625 0.0105 1856006.2
0.00265 0.0106 1751661.161
0.002675 0.0107 1654490.500
0.0027 0.0108 1563918.392
0.002725 0.0109 1479421.407
Ordinaria 0.00275 0.011 1400523.257
0.002775 0.0111 1326790.065
Borda 0.0028 0.0112 1257826.160
0.002825 0.0113 1193270.328
0.00285 0.0114 1132792.471
0.002875 0.0115 1076090.612
0.0029 0.0116 1022888.230
0.002925 0.0117 972931.861
0.00295 0.0118 925988.957
0.002975 0.0119 881845.953
0.003 0.012 840306.541
Tabla 1. Cálculo de Re con la ecuación de Von Karman

Anteriormente se había mencionado lo de las entradas al tanque superior (véase la


figura del problema) y se había realizado el problema con ambas entradas (borda y
ordinaria). Se tenía que para una entrada borda, su factor de fricción Fanning es de
0.0028 y su Re es de 1.2*10^6, mientras que para entrada ordinaria su factor es de
0.00275 y su Re de 1.3*10^6. En la tabla anterior (tabla 1), se muestran los valores
de Re calculados, con entrada ordinaria y con de entrada borda. En la siguiente
tabla se hace una comparación y se estima el error aproximado relativo porcentual
para los valores de Re obtenidos en ambos métodos.
El error aproximado relativo porcentual se calcula de la siguiente manera:

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 11
Grupo 4BV1
Siendo (en este caso):
• V (real) = Re calculado originalmente
• V (aprox) = Re calculado con el método de Von Karman
Tipo de entrada usada
f (F) Re (original) Re (Von Karman) EARP
(tanque superior)
Borda 0.0028 1.3*10^6 1.258*10^6 3.231%
Ordinaria 0.00275 1.2*10^6 1.4*10^6 16.667%
Tabla 2. Error aproximado relativo porcentual de los Re calculados en ambos métodos

También se hace los mismo para comparar los resultados de los flujos volumétricos
obtenidos con ambos métodos (aclarando que los flujos se calcularon con la misma
fórmula del flujo volumétrico descrita anteriormente):
Flujo Volumétrico (ft^3/s)
Tipo de entrada
usada (tanque f (F) Original Con Von Karman EARP
superior)
Borda 0.0028 4.400 4.436 0.818%
Ordinaria 0.00275 4.580 4.940 7.860%
Tabla 3. Error aproximado relativo porcentual de los flujos volumétricos calculados en ambos
métodos

Utilizando la ecuación implícita de Von Karman ya que se asume rugosidad nula


(Valiente, 2016), y para el factor de Fanning (Levenspiel, 2014):
1 1.255
= −4𝑙𝑜𝑔10 ( )
√𝑓𝐹 𝑅𝑒√𝑓𝐹

Donde, de lo presentado anteriormente, se tiene

ρ𝐷 ρ𝐷 2𝑔(𝑧1 − 𝑧2 )
𝑅𝑒 = ⟨𝑣⟩ =
μ μ √ 𝐿 𝐿
1 + ((𝐷 ) + (𝐷 ) ) 𝑓𝐹
𝑒𝑞

Con lo cual puede hacerse

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 12
Grupo 4BV1
1 1.255
Φ(𝑓𝐹 ) = + 4𝑙𝑜𝑔10 =0
√𝑓𝐹 ρ𝐷 2𝑔(𝑧1 − 𝑧2 )
μ √ √𝑓𝐹
𝐿 𝐿
1 + ((𝐷) + (𝐷 ) ) 𝑓𝐹
𝑒𝑞
( )

Con lo que el problema se convierte en uno de buscar una 𝑓𝐹 tal que Φ(𝑓𝐹 ) = 0.
Puesto que esta cuestión es la de encontrar las raíces de una función trascendental,
es preciso hacer uso de un método numérico, donde el método de Newton-Raphson
es de suma utilidad (Stewart, 2016).
Este método se describe como sigue:
Sea que 𝑥𝑛 denote la 𝑛-ésima aproximación a la raíz de una función 𝑓(𝑥); siempre
que 𝑓’(𝑥𝑛 ) ≠ 0 entonces
𝑓(𝑥𝑛 )
𝑥𝑛+1 = 𝑥𝑛 −
𝑓 ′ (𝑥𝑛 )
Donde f’(x) denota la derivada de 𝑓 con respecto de 𝑥.
Otra de las ventajas de este método es que se presta a ser iterado y compilado en
software. Considerando esto, generamos el siguiente código en Matlab:

function [raiz]=niurapmn(x0,f,df,n,tol)
i=1; %Índices.
x(1)=x0; %Semilla o primer aproximación.
e(1)=inf;%El primer error no tiene sentido.
error=tol+1; %La tolerancia es la variación permitida entre las últimas dos
aproximaciones.

while((i<=n+1)&&(error>tol)) %El ciclo while corre hasta que se cumple


cualquiera de las condiciones.
x(i+1)=x(i)-f(x(i))/df(x(i));
e(i+1)=abs(x(i+1)-x(i));
error=e(i+1);
i=i+1;
end
disp("La raíz es:");
raiz=table(x',e')
end

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 13
Grupo 4BV1
Así, tenemos para nuestro caso

1 1.255
Φ= + 4𝑙𝑜𝑔10
√𝑓𝐹 𝜌𝐷 −2𝑔∆𝑧
𝜇 √ 𝐿 𝐿 √𝑓𝐹
1 + 𝑓𝐹 (𝐷 + 4 (𝐷 ) )
𝑒𝑞
( )
Y su correspondiente derivada
𝑑Φ
𝑑𝑓
1 3
= − 𝑓𝐹 −2
2
−1

−1
4 1.255 1.255 −2𝑔∆𝑧 1 3
+ (− 𝑓𝐹 −2 )
𝑙𝑛10 𝜌𝐷 √ 𝐿 𝐿 2
𝜌𝐷 −2𝑔∆𝑧 1 + 𝑓𝐹 (𝐷 + 4 (𝐷 ) )
√𝑓𝐹 𝜇 𝑒𝑞
𝜋 √ 𝐿 𝐿
1 + 𝑓𝐹 (𝐷 + 4 (𝐷 ) ) (
𝑒𝑞
( )
3
−2
1 1 −2𝑔∆𝑧
+ (𝑓𝐹 −2 ) − ( ) ∙ −2𝑔∆𝑧 ∙ −1
2 1 + 𝑓 (𝐿 + 4 (𝐿) )
𝐹 𝐷 𝐷 𝑒𝑞
(

−2
𝐿 𝐿 𝐿 𝐿
∙ (1 + 𝑓𝐹 ( + 4 ( ) )) ( +4( ) )
𝐷 𝐷 𝑒𝑞 𝐷 𝐷 𝑒𝑞

))

Luego, procedemos a introducir estos valores a la suite Matlab


format long %Indica que se generen respuestas de más de cuatro decimales.

%Estas son los valores del problema.


p=62.4;

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 14
Grupo 4BV1
D=0.417;
u=6.72*10^(-4);
g=32.174;
z1=0;
z2=-38;
L=49;
LDeq=76; %Entrada ordinaria.

Phi=@(f) 1/sqrt(f) + 4*log10(1.255/((p*(D/u)*sqrt((-2*g*(z2-z1))/(1+f*(L/D +


4*LDeq))))*sqrt(f))); %La función cuya raíz es nuestro objetivo.

dPhi=@(f) - 1/(2*f^(3/2)) -
(800*2^(1/2)*D*f^(1/2)*p*((251*2^(1/2)*u)/(800*D*f^(3/2)*p*((g*(z1 -
z2))/(f*(4*LDeq + L/D) + 1))^(1/2)) - (251*2^(1/2)*g*u*(4*LDeq + L/D)*(z1 -
z2))/(800*D*f^(1/2)*p*(f*(4*LDeq + L/D) + 1)^2*((g*(z1 - z2))/(f*(4*LDeq +
L/D) + 1))^(3/2)))*((g*(z1 - z2))/(f*(4*LDeq + L/D) +
1))^(1/2))/(251*u*log(10)); %Su correspondiente derivada.

f0=0.00265; %Aproximación inicial, proveniente del análisis previo.

niurapmn(f0,Phi,dPhi,30,1*10^(-8)) %La función programada. Indicamos una


tolerancia de ± 0.00000001 entre la última aproximación y la anterior.

La raíz es:

0.00278500189373404

𝑥𝑛 Error
1 0.00265000000000000 Inf
2 0.00278004923996585 0.000130049239965851
3 0.00278499527912451 4.94603915866202e-06
4 0.00278500189373404 6.61460952244147e-09
Tabla 4. Estimación de factores de fricción de Fanning para entrada ordinaria.

Con lo que posteriormente podemos calcular el correspondiente número de


Reynolds.
fr=0.00278500189373404;
Re= p*(D/u)*sqrt((-2*g*(z2-z1))/(1+fr*(L/D + 4*LDeq)))

Re=

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 15
Grupo 4BV1
1.298647892800684e+06

Mientras que considerando entrada borda.


La raíz es:

0.00280210055198462

𝑥𝑛 Error
1 0.00265000000000000 Inf
2 0.00279584154160369 0.000145841541603691
3 0.00280209004357017 6.24850196648220e-06
4 0.00280210055195501 1.05083848354719e-08
5 0.00280210055198462 2.96112960540551e-14
Tabla 5. Estimación de factores de fricción de Fanning para entrada borda.

Y a su vez
Re=

1.252237205102451e+06

Considerando
𝑅𝑒μ
⟨𝑣⟩ =
ρ𝐷

𝐷 2
𝐴 = π( )
2
𝑅eμDπ
𝑉̇ = ⟨𝑣⟩𝐴 =
ρ4
Entonces
Con entrada ordinaria
ft 3
𝑉̇ = 4.58038553916079
𝑠
Con entrada borda
𝑓𝑡 3
𝑉̇ = 4.57281161365895
𝑠

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 16
Grupo 4BV1
Discusión
De acuerdo con Green (2008), Crane (2010), Valiente (2016) y Levenspiel (2014)
existen varias ecuaciones que nos relacionan al factor de fricción 𝑓𝐹 y al número de
Reynolds 𝑅𝑒, así como también varios métodos para encontrar la relación entre
factor de fricción y pérdidas debido a fricción en componentes.
De acuerdo con Levenspiel (2014) y Crane (2010), la ecuación que más
exactamente describe la relación entre 𝑓𝐹 y 𝑅𝑒 es la de Colebrook-White, basada en
las ecuaciones comprobadas experimentalmente mediante los exhaustivos trabajos
de Nikuradse (Schlichting, 1979). Dicha ecuación se reduce a la propuesta por Von
Karman cuando la rugosidad es cero.
Por otra parte, existen varios métodos para calcular las pérdidas por fricción debido
a accesorios, entre los que encontramos: por longitud equivalente (Levenspiel,
2014; Perry, 20; Crane; 2010;), por cabezas de velocidad (Green, 2008), por
coeficiente de resistencia (Crane, 2010), por factores de pérdidas por fricción
(Kramers, 1958). Estos valores difieren en magnitud debido a las consideraciones
geométricas que tratan. El tipo de flujo también afecta la exactitud de los resultados
que se obtengan con un método u otro.

Figura 4. Factores de pérdida por fricción. Tomado de Bird (2002).

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 17
Grupo 4BV1
Ahora bien, y careciendo en el presente trabajo de datos experimentales para
comparar, asumimos como valor de referencia los valores obtenidos con entrada
ordinaria -puesto que de acuerdo con Levenspiel (2014) una entrada borda es la
que protruye hacia el interior del recipiente- y con el método de Newton-Raphson,
debido a buena convergencia del método y consiguientemente buena
reproducibilidad a partir de sólo datos del problema.
De igual manera, considerando de lo mencionado previamente que buenos
resultados deben cumplir la ecuación de Colebrook, se procedió a evaluar los
factores de fricción y sus 𝑅𝑒 correspondientes para corroborar que el valor de
referencia sea apropiado. Para este análisis se incluye el valor reportado por Bird
(2002) para este problema, encontrado mediante la gráfica de Moody. Se generó
lo siguiente:

Figura 5. Evaluación de idoneidad mediante criterio de raíz a ecuación Von Karman. Resultados a la derecha.

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 18
Grupo 4BV1
Puede observarse que en todos los cálculos con la ecuación de Von Karman, los
resultados aproximan adecuadamente el 0, mientras que en los cálculos con la
gráfica de Moody esta tendencia disminuye.
Luego, procedimos a encontrar la variación de los valores de la velocidad con
respecto al de referencia. Para esto, calculamos los siguientes errores relativos
absolutos aproximados porcentuales

Figura 6. Evaluación de error respecto a valor de referencia. Resultados a la derecha.

De esto se tiene que la variación no es demasiada tanto en resultados del método


de Newton-Raphson como en cálculos determinados mediante la gráfica de Moody
(<8%); no así para el reportado por Bird (2002).
Utilizando los datos en sistema internacional de unidades, así como los valores de
Kramers (1958) para encontrar las pérdidas por fricción por accesorios planteados
por Bird (2002), calculamos el factor de fricción, 𝑅𝑒 y flujo volumétrico en 𝑓𝑡 3 𝑠 −1
mediante el método de Newton-Raphson. Se obtuvo
𝑓𝐹 = 0.00276839986374562
IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato
Materia: Fenómenos de Transporte 19
Grupo 4BV1
𝑅𝑒 = 1.345777121949706𝑒 + 06

𝑉̇ = 4.758014700852419 𝑓𝑡 3 𝑠 −1
Cuya variación respecto al valor de referencia es

Figura 7. Error en flujo volumétrico mediante Newton-Raphson con datos por Bird (2002). Resultados a la
derecha.

Por lo tanto, podemos decir que en general tanto los métodos gráficos como los
computacionales llevan a una respuesta bastante similar, en este caso con no más
de 8% de error. Esto es porque, de acuerdo con Crane (2010) y Green (2008), para
los flujos en régimen turbulento con un 𝑅𝑒 alto, calcular las pérdidas por fricción por
accesorios mediante un método u otro, a pesar de considerar diferentes escenarios
geométricos, no tiene un efecto significativo, debido al carácter aproximadamente
constante de la resistencia hidráulica. Nuestros resultados verifican dicho
enunciado.
Finalmente, revisando el trabajo de Bird (2002), encontramos que la causa de su
desviación es el redondeo, ya que los autores simplificaron la expresión para
calcular el 𝑅𝑒. Esto ilustra la importancia de utilizar tantas cifras como sea posible.

Figura 8. Expresión para el Re empleada por Bird (2002).

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 20
Grupo 4BV1
Conclusiones
Se lograron cumplir los objetivos del proyecto, ya que se realizó el planteamiento
del problema al cual se dio solución identificando las características y propiedades
de este, y con ello se logró hacer uso de los conocimientos obtenidos en el primer
parcial de la materia fenómenos de transporte.
Se encontró que se puede implementar satisfactoriamente un pequeño algoritmo
computacional altamente eficiente en Matlab para resolver la ecuación de Von
Karman, así como una hoja de cálculo en Excel para estudiar su comportamiento
en intervalos definidos. Resolviendo esta ecuación para las condiciones del
𝑓𝑡 3
problema, se calculó un flujo volumétrico de aproximadamente 4.58 . Calculando
𝑠
flujos volumétricos mediante métodos gráficos (Gráfica de Moody) con diversas
correcciones a las pérdidas de por fricción (Longitud equivalente; factores de
pérdida por fricción) y mediante una hoja de cálculo, se encontró que estos diferían
en menos de 8% con respecto al obtenido mediante el algoritmo. Por lo tanto, se
confirma la observación empírica de que flujos turbulentos a 𝑅𝑒 altos en fluidos
newtonianos no dependen fuertemente de la geometría de los accesorios, así como
el que cálculos mediante método gráficos no introducen errores importantes.
Se concluye que la gran diferencia en el flujo volumétrico respecto al valor reportado
en la literatura es debida al redondeo efectuado en el procedimiento ahí referido,
pues esta desaparece al utilizar los valores originales. Indicamos que este
comportamiento debería restringirse.
Finalmente, las limitaciones a este estudio fueron la falta de mediciones
experimentales. Investigaciones posteriores podrían enfocarse a comprobar
mediante mediciones el cumplimiento de la ecuación de Von Karman.

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 21
Grupo 4BV1
Bibliografía
1. Calderón, F., Gallegos L., & Flores, F. (2017). Dificultades conceptuales para
la comprensión de la ecuación de Bernoulli. Revista Eureka sobre Enseñanza
y Divulgación de las Ciencias, 14(2),339-352.
2. Castejón, J. (2014). Resolución aproximada de la capa límite mediante la
ecuación de Kon Karman empleando una ley logarítmica. Universidad
Politécnica de Cartagena.
3. Cengel, Y. (2009). Termodinámica. (7. a ed.) Edición. McGraw Hill Education.
4. Durán, A., Segovia, G., Medrano, V. (2014) Evaluación de ecuaciones de
factor de fricción explícito para tuberías, Educación Química, Volumen 25,
Issue 2, Pág 128-134.
5. Felder, R. (2004) Principios elementales de los procesos químicos.
Elementory principies of chemicol processes. (3a. ed). México: Llmuso Wlley.
6. Levenspiel. (2014). Flujo de fluidos e intercambio de calor. Barcelona,
España. Editorial Reverté.
7. Bird. (2002). Transport phenomena. (2nd edition). USA. John Wiley & Sons
Inc. Editorial.
8. Valiente, A. (2016). Ingeniería de fluidos. CDMX, México: UNAM.
9. Stewart, J. (2016). Single Variable Calculus: Early Transcendentals.
Cengage Learning.
10. Green, D. W., Poling, B. E., Thomson, G. H., Friend, D. G., Rowley, R. L.,
Wilding, W. V., Perry, R. H. (2008). Perry’s Chemical Engineers’ Handbook.
New York, USA: McGraw-Hill.
11. Crane Co. Engineering and Research Division. (2010). Flow of fluids through
valves, fittings and pipe. Crane.
12. Schlichting, H. (1979) Boundary Layer Theory, 7th edn., transl. by J. Kestin.
New York, USA: McGraw-Hill.
13. Kramers, H. (1958). Physische transportverschijnselen. Holland: Technische
Hogeschool Delft.

IPN Unidad Profesional Interdisciplinaria de Ingeniería campus Guanajuato


Materia: Fenómenos de Transporte 22
Grupo 4BV1

También podría gustarte