Caso de Estudio - Método de Von Karman - 4BV1 - E#1 PDF
Caso de Estudio - Método de Von Karman - 4BV1 - E#1 PDF
Caso de Estudio - Método de Von Karman - 4BV1 - E#1 PDF
NACIONAL
Unidad Profesional Interdisciplinaria de
Ingeniería
Campus Guanajuato
Ingeniería biotecnológica
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í
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
Consideraciones:
•
• (Solo existe la presión atmosférica)
• No cambia de velocidad
• Tubo liso
• Fluido Newtoniano (agua)
Siendo:
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:
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:
• 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.
• 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
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
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
ρ𝐷 ρ𝐷 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.
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( ) )
𝐷 𝐷 𝑒𝑞 𝐷 𝐷 𝑒𝑞
))
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.
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.
Re=
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
𝑠
Figura 5. Evaluación de idoneidad mediante criterio de raíz a ecuación Von Karman. Resultados a la derecha.
𝑉̇ = 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.