Estudio Analítico de Flujo Monofásico en Tuberías
Estudio Analítico de Flujo Monofásico en Tuberías
Estudio Analítico de Flujo Monofásico en Tuberías
1.1 Estacionario-Transitorio
Un flujo estacionario es aquel en el que las variables que lo definen no cambian localmente con
el tiempo; uno transitorio en el que sí lo hacen.
1.2 Incomprensible-Compresible
Un fluido incompresible es aquel en el que podemos considerar su densidad como constante;
en el caso contrario el fluido se denomina compresible.
1.3 Ideal-Viscoso
Un fluido ideal es aquel en el que el efecto de la viscosidad es tan débil que puede suponerse
nula. Esto significa que podemos despreciar la fuerza de rozamiento interna entre las diferentes
capas de fluido. Por otro lado, los fluidos en los que el efecto dominante es precisamente la
viscosidad, se los denomina viscosos.
1.4 Newtonianos
Los fluidos newtonianos son aquellos en los que la viscosidad es constante, es decir, la relación
entre la fuerza por unidad de superficie que una capa de fluido ejerce sobre las vecinas y la
velocidad de deformación es constante.
𝛿
∫ 𝜌𝑑∀ + ∮ 𝜌𝑣⃗. 𝑑𝐴⃗ = 0 1.1
𝛿𝑡
𝛿𝜌
+ 𝛻. 𝜌𝑣⃗ = 0 1.2
𝛿𝑡
Luego presentamos las siguientes igualdades:
𝑑𝜌 𝛿𝜌
= + (𝑣⃗. 𝛻)𝜌 𝛻. 𝜌𝑣⃗ = (𝑣⃗. 𝛻)𝜌 + 𝜌(𝛻. 𝑣⃗)
𝑑𝑡 𝛿𝑡
Reemplazamos en la ecuación 1.2, obtenemos la ecuación de conservación en forma no-
conservativa:
𝑑𝜌
+ 𝜌(𝛻. 𝑣⃗) = 0 1.3
𝑑𝑡
𝑑𝐹𝑏𝑜𝑑𝑦 = 𝜌𝑔𝑑𝑥𝑑𝑦𝑑𝑧
−𝑝 0 0
𝑝 = [ 0 −𝑝 0]
0 0 −𝑝
𝛿𝑝
𝑑𝐹𝑥 = − 𝑑𝑥𝑑𝑦𝑑𝑧
𝛿𝑥
𝛿𝑝
𝑑𝐹 = (−𝛻𝑝)𝑑𝑥𝑑𝑦𝑑𝑧 𝑑𝐹𝑦 = − 𝑑𝑥𝑑𝑦𝑑𝑧
𝛿𝑦
𝛿𝑝
𝑑𝐹𝑧 = − 𝑑𝑥𝑑𝑦𝑑𝑧
𝛿𝑧
𝑑𝑣⃗
𝜌 = 𝜌𝑔⃗ + 𝛻. 𝜎̿𝑖𝑗 − 𝛻𝑝̿
𝑑𝑡
Se presenta la ecuación en coordenadas cartesianas:
2
𝜏𝑖𝑗 = 𝜇 [(𝛿𝑖 𝑣⃗𝑗 + 𝛿𝑗 𝑣⃗𝑖 ) − (𝛻. 𝑣⃗)𝛾𝑖𝑗 ] 𝛾𝑖𝑗 = 0 → 𝑖 ≠ 𝑗 𝛾𝑖𝑗 = 1 → 𝑖 = 𝑗
3
𝛿𝑣 𝛿𝑢 𝛿𝑢 2
𝜏𝑥𝑦 = 𝜏𝑦𝑥 = 𝜇 ( + ) 𝜎𝑥𝑥 = 𝜇 [2 𝛿𝑥 − 3 (𝛻. 𝑣⃗)]
𝛿𝑥 𝛿𝑦
𝛿𝑤 𝛿𝑣 𝛿𝑣 2
𝜏𝑦𝑧 = 𝜏𝑧𝑦 = 𝜇 ( + ) 𝜎𝑦𝑦 = 𝜇 [2 − (𝛻. 𝑣⃗)]
𝛿𝑦 𝛿𝑧 𝛿𝑦 3
𝛿𝑢 𝛿𝑤 𝛿𝑤 2
𝜏𝑧𝑥 = 𝜏𝑥𝑧 = 𝜇 ( + ) 𝜎𝑧𝑧 = 𝜇 [2 − (𝛻. 𝑣⃗)]
𝛿𝑧 𝛿𝑥 𝛿𝑧 3
Reemplazamos todas las igualdades en la ecuación general de conservación de momento:
𝑑𝑢 𝛿𝑝 𝛿 𝛿𝑢 2 𝛿 𝛿𝑣 𝛿𝑢 𝛿 𝛿𝑢 𝛿𝑤
𝜌 = 𝜌𝑔𝑥 − + {𝜇 [2 − (𝛻. 𝑣⃗)] } + {𝜇 ( + )} + {𝜇 ( + )}
𝑑𝑡 𝛿𝑥 𝛿𝑥 𝛿𝑥 3 𝛿𝑦 𝛿𝑥 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑥
𝑑𝑣 𝛿𝑝 𝛿 𝛿𝑣 𝛿𝑢 𝛿 𝛿𝑣 2 𝛿 𝛿𝑤 𝛿𝑣
𝜌 = 𝜌𝑔𝑦 − + {𝜇 ( + )} + {𝜇 [2 − (𝛻. 𝑣⃗)] } + {𝜇 ( + )}
𝑑𝑡 𝛿𝑦 𝛿𝑦 𝛿𝑥 𝛿𝑦 𝛿𝑦 𝛿𝑦 3 𝛿𝑦 𝛿𝑦 𝛿𝑧
𝑑𝑤 𝛿𝑝 𝛿 𝛿𝑢 𝛿𝑤 𝛿 𝛿𝑤 𝛿𝑣 𝛿 𝛿𝑤 2
𝜌 = 𝜌𝑔𝑧 − + {𝜇 ( + ) } + {𝜇 ( + )} + {𝜇 [2 − (𝛻. 𝑣⃗)]}
𝑑𝑡 𝛿𝑧 𝛿𝑧 𝛿𝑧 𝛿𝑥 𝛿𝑧 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑧 3
(𝛻. 𝑣⃗) = 0
𝛿𝑣 𝛿𝑢 𝛿𝑢
𝜏𝑥𝑦 = 𝜏𝑦𝑥 = 𝜇 ( + ) 𝜎𝑥𝑥 = 𝜇 (2 𝛿𝑥 )
𝛿𝑥 𝛿𝑦
𝛿𝑤 𝛿𝑣 𝛿𝑣
𝜏𝑦𝑧 = 𝜏𝑧𝑦 = 𝜇 ( + ) 𝜎𝑦𝑦 = 𝜇 (2 )
𝛿𝑦 𝛿𝑧 𝛿𝑦
𝛿𝑢 𝛿𝑤 𝛿𝑤
𝜏𝑧𝑥 = 𝜏𝑥𝑧 = 𝜇 ( + ) 𝜎𝑧𝑧 = 𝜇 (2 )
𝛿𝑧 𝛿𝑥 𝛿𝑧
Reemplazamos las variables en la ecuación general. Solo se desarrollar para la coordenada en x:
𝑑𝑢 𝛿𝑝 𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿𝑢 𝛿 𝛿𝑢 𝛿𝑤
𝜌 = 𝜌𝑔𝑥 − + {𝜇 (2 )} + {𝜇 ( + )} + {𝜇 ( + )}
𝑑𝑡 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑥
𝑑𝑢 𝛿𝑝 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑤
𝜌 = 𝜌𝑔𝑥 − +𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )
𝑑𝑡 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑦 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑧 𝛿𝑥
Ordenando los términos se tiene la siguiente ecuación:
𝑑𝑢 𝛿𝑝 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿 𝛿𝑤
𝜌 = 𝜌𝑔𝑥 − +𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )+𝜇 ( )
𝑑𝑡 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑧 𝛿𝑥
De la expresión anterior se analiza los seis últimos términos:
𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 𝛿𝑢 𝛿 2𝑢 𝛿 2𝑢 𝛿 2𝑢
𝜇 ( )+𝜇 ( )+𝜇 ( )=𝜇 2+𝜇 2+𝜇 2
𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑦 𝛿𝑧 𝛿𝑧 𝛿𝑥 𝛿𝑦 𝛿𝑧
𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿 𝛿𝑤 𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿 𝛿𝑤
𝜇 ( )+𝜇 ( )+𝜇 ( )= 𝜇 ( )+𝜇 ( )+𝜇 ( )
𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑧 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑧
Como se aprecia en el último término, el orden de la derivada es alterada en tanto las
componentes de velocidad sean funciones suaves. Desarrollamos las ecuaciones mostradas:
𝛿 2𝑢 𝛿 2𝑢 𝛿 2𝑢
𝜇 + 𝜇 + 𝜇 = 𝜇𝛻 2 𝑢
𝛿𝑥 2 𝛿𝑦 2 𝛿𝑧 2
𝛿 𝛿𝑢 𝛿 𝛿𝑣 𝛿 𝛿𝑤 𝛿 𝛿𝑢 𝛿𝑣 𝛿𝑤 𝛿
𝜇 ( )+𝜇 ( )+𝜇 ( )=𝜇 ( + + ) = 𝜇 (𝛻. 𝑣⃗)
𝛿𝑥 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑥 𝛿𝑧 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑧 𝛿𝑥
El segundo término es nulo por tratarse de un fluido incompresible, por lo tanto, la expresión
queda reducido:
𝑑𝑢 𝛿𝑝 𝛿2𝑢 𝛿 2𝑢 𝛿 2𝑢
𝜌 = 𝜌𝑔𝑥 − +𝜇 2+𝜇 2+𝜇 2
𝑑𝑡 𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑧
𝑑𝑣 𝛿𝑝 𝛿2𝑣 𝛿 2𝑣 𝛿 2𝑣
𝜌 = 𝜌𝑔𝑦 − +𝜇 2+𝜇 2+𝜇 2
𝑑𝑡 𝛿𝑦 𝛿𝑥 𝛿𝑦 𝛿𝑧
𝑑𝑤 𝛿𝑝 𝛿 2𝑤 𝛿 2𝑤 𝛿 2𝑤
𝜌 = 𝜌𝑔𝑧 − +𝜇 2 +𝜇 2 +𝜇 2
𝑑𝑡 𝛿𝑧 𝛿𝑥 𝛿𝑦 𝛿𝑧
También se escribe en forma vectorial:
⃗⃗
𝒅𝒗
𝝆 = 𝝆𝒈 ̿ + 𝝁𝜵𝟐 𝒗
⃗⃗⃗ − 𝜵𝒑 ⃗⃗
𝒅𝒕