Práctica Mecánica de Fluidos
Práctica Mecánica de Fluidos
Práctica Mecánica de Fluidos
Mecánica de fluidos
11/09/2020
En la instalación de transporte de agua a 20 °C (Figura 1) desde un lago (1) hasta un tanque de
almacenamiento (4) se emplea una bomba centrífuga (2) → (3) con una curva de desempeño
dada por el fabricante como se muestra (Figura 2). La elevación de la salida del tubo (4) tiene
una distancia de H = 28 [m] por arriba de la superficie del lago (1).
El tramo de tubería (1_2) y (3_4), de diámetro interno D tubo = 0.3 m, tiene una longitud total L =
150 m. La constante de pérdidas (Para cálculo con cabeza de velocidad) para todo el tramo de
tubería es de k = 85. Al encender la bomba, se establece un flujo variable que después de un
intervalo de tiempo se estabiliza.
A) Deduzca las ecuaciones diferenciales del proceso.
Mediante un modelo en simulink de MatLab, evalué para dos tipos de arranque de la bomba:
escalon y rampa (de 3 a 5 s) (Figura 3), la respuesta dinámica del sistema (Gráficas de,Hwb, Q(t),
W(t) en todo el intervalo de tiempo de variación.
SOLUCIÓN.
B) APLICACIÓN SIMULINK
DATOS
𝐷 = 0.3 𝑚
𝐿 = 150 𝑚
Δ𝐻 = 28 𝑚
𝐾 = 85
La ecuación del rendimiento de la bomba con la ayuda de la herramienta Lookup Table de EES se puede
estimar como:
𝑄 =𝑉∗𝐴
Con estos datos a la mano se puede plantear entonces la ley de Bernoulli entre los puntos (1) y (4)
obteniendo:
𝑃1 𝑉12 𝑃4 𝑉42
+ + 𝑍1 + 𝐻𝑤𝑏 = + + 𝑍4 + 𝐻𝑓
𝛾 2𝑔 𝛾 2𝑔
Como el punto (1) está ubicado sobre la superficie del lago, su presión y velocidad serán cero, del mismo
modo como el punto (4) de la tubería está abierta a la atmósfera, su presión también será cero.
Reorganizando la ecuación se obtiene:
𝑉42
𝐻𝑤𝑏 = + Δ𝑍 + 𝐻𝑓
2𝑔
Las pérdidas de la tubería se pueden calcular usando el factor de pérdidas respecto de la cabeza de
velocidad:
𝑉2
𝐻𝑓 = 𝐾 ∗
2𝑔
Resolviendo simultáneamente estas 4 ecuaciones, se puede determinar el estado estacionario de la
instalación:
𝑊𝑏 = 𝛾 ∗ 𝑄 ∗ 𝐻
El peso específico del agua a 20° es:
𝑁
𝛾 = 9793
𝑚3
D. Potencia de la bomba
𝑊𝑏 = 33640 𝑊 = 44.11𝐻𝑃
ESTADO TRANSITORIO
Cuando el sistema hasta ahora está arrancando, el líquido dentro de las tuberías se encuentra en reposo,
así que, si se hace un volumen de control en toda la tubería, la primera ley de la termodinámica dirá:
𝜋𝐷 2
𝑚 = 𝜌∗𝑉 =𝜌 𝐿 = 10584 𝐾𝑔
4
Reescribiendo la ley de conservación de energía:
𝑑𝑉
𝑚∗ = [𝑊̇𝑏 − (Δ𝐻 − 𝐻𝑓 )𝑄𝛾] ∗ Δ𝑇
𝑑𝑡
𝑑𝑉 [𝑊̇𝑏 − (Δ𝐻 + 𝐻𝑓 )𝑄𝛾]
= ∗ Δ𝑇
𝑑𝑡 𝑚
Llevando esta ecuación diferencial a Simulink se puede analizar el comportamiento transitorio del
sistema
Así se verá la implementación en Simulink. La función principal será la que tiene programada nuestra
ecuación diferencial:
Esta función tomará como entradas la potencia de la bomba, la altura de pérdidas y la altura, así como el
caudal, la masa y el peso específico.
La Potencia será en el primer caso una función escalón (STEP) donde el máximo valor será el valor de la
potencia del estado estable del sistema, por su parte, en el segundo caso su entrada será una función
lineal (RAMP) cuyo máximo valor será la potencia del estado estable del sistema y este valor se
alcanzará a los 3 segundos.
C. TIEMPO DE ESTABILIZACIÓN
Para el primer caso, donde la potencia es una función escalón (STEP), su tiempo de estabilización será:
Por su parte, el segundo caso, donde la bomba arranca con una rampa, su tiempo de estabilización es:
E. PREGUNTA ABIERTA: Análisis adicional:
Finalmente podemos deducir que, de acuerdo a los resultados para el primer caso, donde la potencia es
una función escalón, su estabilización resulta en un tiempo menor (aproximadamente entre 2 y 2.5
segundos) con respecto a donde la bomba arranca con una rampa que está aproximada entre 4 y 4.5
segundos. De este modo concluimos que resulta más optima la potencia en forma escalón.