• Author / Uploaded
• tbautistadiego

• Categories
• Viscosidad
• Líquidos
• Presión
• Gases
• Densidad

Citation preview

ESCURRIMIENTO DE FLUIDOS Aplicaciones

ALEJANDRO REYES SALINAS

ESCURRIMIENTO DE FLUIDOS Aplicaciones

© Editorial Universidad de Santiago de Chile Av. Libertador Bernardo O`Higgins #2229 Santiago de Chile Tel.: 56-2-7180080 www.editorial.usach.cl [email protected] © Alejandro Reyes Salinas Inscripción Nº: 194.968 I.S.B.N.: 978-956-303-100-3 Portada y diseño: Andrea Meza Vergara Diagramación: Andrea Meza Vergara Primera edición, septiembre de 2010 Impreso en Gráfica LOM Ninguna parte de esta publicación puede ser reproducida, almacenada o transmitida en manera alguna ni por ningún medio, ya sea eléctrico, químico o mecánico, óptico, de grabación o de fotocopia, sin permiso previo de la editorial. Impreso en Chile.

ÍNDICE Prefacio

11

CAPÍTULO 1 CONCEPTOS GENERALES DE FLUIDOS 1.1. Concepto de fluido 1.2. Propiedades de los fluidos 1.3. Hidrostática Ejercicios

15 15 17 28 33

CAPÍTULO 2 BALANCES MACROSCÓPICOS 2.1. Ecuación de balance 2.2. Balance de masa 2.3. Balance de energía Ejercicios

37 37 38 42 48

CAPÍTULO 3 FLUJO DE FLUIDOS 3.1. Naturaleza del flujo de fluidos 3.2. Aplicaciones de la ecuación de Bernoulli, en situaciones que la fricción es despreciable 3.3. Definición de factor de fricción 3.4. Descripción de cañerías, válvulas y accesorios 3.5. Evaluación de pérdidas en válvulas y accesorios 3.6. Aplicaciones que involucran la evaluación de factores de fricción 3.7. Evaluación de pérdidas en ductos no circulares 3.8. Diámetro óptimo económico (DOE) 3.9. Situaciones complejas 3.10. Fluidos no newtonianos Ejercicios

51 51 54 63 68 71 74 82 84 88 99 104

CAPÍTULO 4 MEDIDORES DE FLUJO 4.1. Medidores que funcionan en base a principios fluido-dinámicos Tubo de Pitot Placa orificio, tobera y Venturi Ejercicios

116 116 120 148

CAPÍTULO 5 TRANSPORTE DE LÍQUIDOS A.1 Bombas alternativas o recíprocas A.2 Bombas rotatorias B.1 Bombas centrífugas Teoría de bombas centrífugas Funcionamiento real de una bomba centrífuga Determinación del caudal de operación Selección de bombas centrífugas B.2 Bombas con efectos especiales Ejercicios

151 151 156 162 170 177 186 194 195 205

CAPÍTULO 6 FLUJO DE GASES 6.1. Flujo isotérmico 6.2. Flujo adiabático 6.3. Equipos impulsores de gases Ejercicios

209 210 215 219 227

CAPÍTULO 7 CARACTERIZACIÓN DE PARTÍCULAS SÓLIDAS 7.1. Tamaño de partículas 7.2. Forma de las partículas Ejercicio

229 229 235 241

CAPÍTULO 8 FACTORES DE FRICCIÓN EN MEDIOS POROSOS 8.1. Teoría del conjunto de tubos 8.2. Teoría del medio continuo Ejercicios

243 243 248 255

111

CAPÍTULO 9 FILTRACIÓN 9.1. Filtros de lecho profundo 9.2. Filtros sobre superficies 9.3. Teoría de filtración de lechos profundos 9.4. Teoría de filtración sobre superficies 9.5. Filtración de tortas compresibles 9.6. Auxiliares filtrantes Ejercicios

257 257 263 275 279 307 309 313

CAPÍTULO 10 FLUJO ALREDEDOR DE OBJETOS SUMERGIDOS 10.1. Velocidad Terminal 10.2. Coeficientes de arrastre Ejercicios

317 319 320 323

CAPÍTULO 11 FLUIDIZACIÓN 11.1. Introducción a la fluidización 11.2. Aplicaciones industriales de la fluidización 11.3. Calidad de la fluidización 11.4. Equipos de fluidización 11.5. Parámetros de diseño y de operación de lechos fluidizados Ejercicios

325 325 330 333 335 336 343

CAPÍTULO 12 SEDIMENTACIÓN Sedimentación libre Sedimentación retardada Sedimentación floculada 12.1. Análisis de sedimentadores con sedimentación retardada 12.2. Dimensionamiento de espesadores continuos Ejercicios

345 347 348 349 350 360 376

CAPÍTULO 13 SEPARACIONES GAS-PARTÍCULA 13.1. Separadores gravitatorios 13.2. Separadores centrífugos (ciclones) Ejercicios

377 381 385 392

CAPÍTULO 14 TRANSPORTE HIDRÁULICO DE SÓLIDOS EN TUBERÍAS 14.1. Introducción 14.2. Modelos para predecir pérdidas de carga Modelo de Durand Modelo de Newitt Ejercicios

395 395 398 398 401 411

Apéndice

413

Capítulo 1: Conceptos generales de fluidos

PREFACIO Con el propósito de reunir los aspectos más importantes, originalmente con diversas nomenclaturas propias de cada autor, se presenta este texto que pretende entregar en forma coherente y sencilla los temas más relevantes del Escurrimiento de Fluidos y sus aplicaciones. Este material ha sido utilizado en cursos de Mecánica de Fluidos y de Separaciones Fluidopartícula, de la carrera de Ingeniería Civil Química y de Fluidodinámica de la carrera de Ingeniería en Biotecnología de la Universidad de Santiago de Chile, USACH, y en programas de Postítulo y cursos cerrados de formación profesional ofrecidos por el DIQ/USACH. El autor agradece a los colegas del área de Operaciones Unitarias del Departamento de Ingeniería Química de la USACH por haberme facilitado diversos apuntes de clases, en especial al Profesor Rolando Vega. También se agradece a los alumnos ayudantes de diversas promociones, que han colaborado revisando los ejercicios propuestos.

11

A mi esposa Elizabeth que me incentivó y apoyó incondicionalmente en la preparación de este texto, y a nuestros hijos Esteban, Alejandro y Daniela.

Capítulo 1: Conceptos generales de fluidos

CAPÍTULO 1 CONCEPTOS GENERALES DE FLUIDOS

1.1.

Concepto de fluido

Previo a la definición del término “fluido”, se aclarará el concepto de esfuerzo de corte, el cual es fácil de comprender analizando las 3 situaciones físicas siguientes:

Figura 1.1. Tipos de fuerza ejercida por el peso P.

En la Figura 1.1.a, la cuerda está manteniendo al peso P, el que ejerce una fuerza que tiende a estirar la cuerda. Por otro lado, una tensión es el cuociente entre una fuerza y el área sobre la cual se ejerce esta fuerza. La fuerza que trata de estirar la cuerda se llama fuerza de tracción y la tensión producida se llama tensión de tracción. En Figura 1.1.b, una columna mantiene un peso P, el que ejerce una fuerza que tiende a comprimir la columna. Esta clase de fuerza se llama fuerza compresiva y la fuerza en la columna se llama tensión compresiva. En Figura 1.1.c, un pegamento mantiene adherido a las paredes el peso P, el que ejerce una fuerza que tiende a deslizarlo, bajando por las paredes y así tensiona el pegamento. Esta fuerza, que tiende a deslizar una

15

Escurrimiento de fluidos

USACH

superficie paralela a una superficie adyacente se denomina fuerza de corte y la tensión en el pegamento, vale decir la fuerza dividida por el área de pegamento se llama tensión de cortadura o esfuerzo de corte (τ). Luego, (1.1)

Con el propósito de diferenciar fluidos y sólidos, es posible establecer que los sólidos son sustancias que pueden resistir permanentemente grandes esfuerzos de corte. Al ser sometidos a una fuerza de corte pueden moverse solamente una pequeña distancia (por deformación elástica) y entonces dejarán de moverse. Si la fuerza de corte es aún más grande se cortarán, es decir se rompen (Figura 1.2).

Figura 1.2. Aplicación de una fuerza de corte sobre un sólido.

Los materiales que son fluidos no pueden resistir permanentemente una fuerza de corte, no importando cuan pequeña sea ésta. Cuando son sometidos a dichas fuerzas ellos comienzan a moverse y continúan en movimiento mientras se aplica la fuerza (es decir fluyen). Dentro de este contexto el término fluido es general, incluyendo tanto a líquidos como a gases.

16

Capítulo 1: Conceptos generales de fluidos

1.2.

Propiedades de los fluidos

Entre las propiedades que caracterizan a los fluidos se tiene la densidad, gravedad específica, viscosidad, tensión superficial, presión de vapor, conductividad térmica, etc. Definiremos a continuación algunas de éstas propiedades físicas y conceptos fundamentales en mecánica de fluidos: A) Densidad (ρ): Se define como el cuociente entre masa y volumen.

ρ=

masa volumen

(1.2)

Existen numerosas técnicas experimentales e instrumentos que permiten evaluar la densidad de gases, líquidos y sólidos, existiendo los densímetros, picnómetros de líquidos y gases. Frecuentemente los textos de mecánica de fluidos y manuales en general, incluyen tablas con valores de densidad para diferentes líquidos y sólidos, referidos a 1 atmósfera y temperatura de 4 ºC o 20 ºC. Para el caso de gases, dado que su densidad depende fuertemente de la presión y temperatura, es práctica común emplear relaciones P-V-T (ecuaciones de estado) para tal efecto. Si las condiciones de presión no son severas puede usarse la ecuación de gas ideal:

ρ=

p⋅M R ⋅T

en que: p = presión absoluta M = peso molecular T = temperatura absoluta R = constante de los gases

17

(1.3)

Escurrimiento de fluidos

USACH

Ejemplo 1.1 En un ducto de ventilación circula aire a 20 ºC y P = 1.1 at. ¿Cuál es la densidad del aire? Solución: Considerando que el aire se comporta como gas ideal, usaremos la ecuación de gases ideales, ecuación 1.3. Introduciendo los valores respectivos, se obtiene ρ = 1,33 (kg/m3). Note que en esta situación, si la presión aumenta al doble, la densidad del aire también se duplica. Ejemplo 1.2 Determine la densidad de metano a 50 ºC y 2.5 atmósferas. Solución:

B) Gravedad específica (G.E.): Gravedad específica se define como el cuociente entre la densidad del fluido y la densidad de un fluido de referencia, en condiciones establecidas. Vale decir: (1.4) Esta definición tiene el mérito de ser una razón, por lo tanto, es un número puro, sin dimensión. Debe tenerse en cuenta el fluido de referencia y las condiciones de presión y temperatura de éste. Normalmente el fluido de referencia es agua y la presión de referencia es 1 atmósfera. Como temperatura de referencia para sólidos y líquidos se utiliza 4 ºC (39 ºF) o 21.1 ºC (70 ºF).

18

Capítulo 1: Conceptos generales de fluidos

Para gases suele utilizarse como fluido de referencia aire a 1 atmósfera, a las temperaturas ya señaladas. Algunos procesos industriales utilizan escalas especiales de densidades de fluidos, las que usualmente se refieren como “gravedades”. Algunas de ellas son: “gravedad API” para aceites y petróleo, “gravedad Brix” para industria del azúcar y “gravedad Baumé” para ácido sulfúrico. Cada una de estas escalas es convertible en densidad, existiendo tablas y fórmulas para este efecto. C) Viscosidad: La viscosidad es una medida de la resistencia a fluir. Por ejemplo el agua es menos viscosa que la miel. Una definición más precisa para viscosidad es posible en términos del siguiente experimento:

Consideremos dos placas largas, separadas por un pequeño film de líquido. Si se mueve la placa superior en la dirección X con velocidad Vo, se requerirá una fuerza para vencer la fricción en el fluido entre las placas. Esta fuerza será diferente para diferentes velocidades, diferentes fluidos y diferentes separaciones entre placas. Midiendo la fuerza por unidad de área de la placa se obtiene el esfuerzo de corte (μ).

Figura 1.3. Experimento de la placa deslizante.

19

Escurrimiento de fluidos

USACH

Experimentalmente se ha demostrado que a bajos valores de Vo el perfil de velocidades en el fluido entre placas es lineal, es decir:

También se ha demostrado experimentalmente que para la mayoría de los fluidos los resultados de este experimento pueden mostrarse en forma más conveniente graficando τ vs dV/dy. En la Figura 1.4 se describen los comportamientos más frecuentes:

Figura 1.4. Comportamiento reológico de diferentes fluidos, en función del esfuerzo de corte (τ) vs dV/dy.

El comportamiento más común es el representado por la línea recta que parte del origen, llamada newtoniana porque describe la ley de Newton de la viscosidad: τ = −µ

dV dy

Aquí μ corresponde a la viscosidad o coeficiente de viscosidad.

20

(1.5)

Capítulo 1: Conceptos generales de fluidos

Los gases como el aire, presentan bajos valores de la viscosidad, quedando la recta μ vs dV/dy muy cerca del eje horizontal. Para fluidos tales como jarabe de maíz el valor de μ es muy grande, pasando la línea recta cercana al eje τ. Los fluidos que siguen el comportamiento ya descrito son llamados fluidos newtonianos. Todos los otros son llamados fluidos no-newtonianos. ¿Qué fluidos son no-newtonianos? Aquellos formados de moléculas o partículas mucho mayores que las moléculas de agua, tales como suspensiones concentradas y pastas en general. Dimensiones y unidades de la viscosidad (μ). Las dimensiones de la viscosidad se pueden deducir de la ley de Newton de la viscosidad:

Utilizando dimensiones de masa, longitud, tiempo (es decir, un sistema absoluto), tenemos:

En el sistema centímetro-gramo-segundo (CGS), la unidad de viscosidad es el poise:

21

Escurrimiento de fluidos

USACH

Por comodidad se emplea frecuentemente el centipoise. 1 centipoise = 0.01 poise = 0.001 [kg /(m s)] Variación de la viscosidad con la temperatura y presión La variación de la viscosidad con la temperatura es significativa, por lo que se debe considerar este efecto en el diseño de equipos. La viscosidad de líquidos disminuye al aumentar la temperatura, mientras que para gases el efecto es inverso. El efecto de la presión en la viscosidad de líquidos normalmente es despreciable. Para gases a presiones elevadas, la viscosidad aumenta con la presión. En tablas y nomogramas se dispone de valores de la viscosidad de líquidos y gases, normalmente referidos a 1 atmósfera. En literatura especializada se dispone de ecuaciones para determinar la viscosidad en función de la presión. Fluidos no newtonianos Existen numerosos fluidos que no siguen el comportamiento newtoniano, tales como suspensiones concentradas y pastas en general. Para estos fluidos se han propuesto ecuaciones que representan su comportamiento específico en términos de la relación entre esfuerzo de corte (τ) y el gradiente de velocidad. A continuación se muestran las ecuaciones de dos modelos: Modelo de Bingham Toda sustancia que sigue este modelo de dos parámetros se denomina plástico de Bingham; permanece rígida mientras el esfuerzo cortante es menor de un determinado τo, comportándose como fluido newtoniano para τ > τo. 22

Capítulo 1: Conceptos generales de fluidos

(1.6) y (1.7)

Tabla 1.1. Parámetros para fluidos de Bingham (Levenspiel, 1993). Fluido

Tensión de fluencia, τo (Pascales)

Ketchup (30 ºC) Mostaza (30 ºC) Oleomargarina (30 ºC) Mayonesa

14 38 51 85

Viscosidad plástica, (kg/m s) 0.08 0.25 0.72 0.63

Modelo de Ostwald de Waele (o ley de la potencia) En este modelo, K es el índice de consistencia de flujo y n es el índice de comportamiento de flujo. Si n1, el fluido es dilatante. En Tabla 1.2. se presentan parámetros para algunos fluidos que siguen la ley de la potencia (1.8)

23

Escurrimiento de fluidos

USACH

Tabla 1.2. Parámetros de fluidos que siguen la ley de la potencia (Levenspiel, 1993). Fluido Compota de manzana (24 ºC) Papilla de plátanos (24 ºC) Sangre humana Sopas y salsas Jugo de tomate (5.8% de sólidos, 32 ºC) Jugo de tomate (30% de sólidos, 32 ºC) Pasta de papel en agua (4%) Cal en agua (33%) Carboximetilcelulosa en agua (15%)

n (-) 0.41 0.46 0.89 0.51 0.59 0.40 0.575 0.171 0.554

K (kg/ms2-n) 0.66 6.5 0.00384 3.6 - 5.6 0.22 18.7 20.7 7.16 3.13

D) Presión: La fuerza normal que actúa sobre un área plana dividida por el área es la presión media. Si el área tiende a cero, se habla de presión puntual, la cual es la misma en todas las direcciones (x, y, z). En situaciones de mecánica de fluidos es frecuente trabajar con escalas de presión absolutas y escalas relativas como por ejemplo presiones manométricas y presiones de vacío. Las relaciones entre ellas son directas: Pmanométrica = P de vacío

=

P absoluta

-

Patmosférica local

(1.9)

Pabsoluta

(1.10)

Patmosférica local -

En la Figura 1.5 se muestran gráficamente estas relaciones. Debe observarse que las presiones absolutas siempre serán positivas, no así las manométricas que podrían ser negativas.

24

Capítulo 1: Conceptos generales de fluidos

Figura 1.5. Escalas de presión (Absoluta y relativas).

Instrumentos para medir presión. El manómetro Bourdon (Figura 1.6) es uno de los aparatos típicos que se usan para medir presiones manométricas. El elemento que soporta la presión es un tubo metálico curvado, cerrado por un extremo, y que por el otro se conecta al recipiente que contiene el fluido cuya presión va a medirse. Cuando la presión interna aumenta, el tubo tiende a enderezarse, tirando de un eslabón que actúa sobre la aguja obligándola a moverse. En la esfera se lee cero cuando en el interior y en el exterior del tubo reina la misma presión, cualesquiera que sean sus valores particulares. La esfera puede ser graduada con las unidades que se prefieran, tales como mm de mercurio o metros de agua. Por su construcción, este manómetro sirve para medir presiones relativas a la presión del medio que rodea al tubo, que suele ser la presión atmosférica local.

25

Escurrimiento de fluidos

USACH

Figura 1.6. Instrumentos de medición. a) manómetro de Bourdon, b) Vista interna del manómetro de Bourdon, c) Calibrador de peso muerto.

El manómetro de Bourdon, por ser un dispositivo de tipo mecánico debería ser calibrado periódicamente utilizando, por ejemplo, un calibrador de peso muerto (Figura 1.6.c). Otros instrumentos para medir presión funcionan en base a diafragmas (Figura 1.7), transductores piezoresistivos, transductores piezoeléctricos, columnas de líquidos manométricos. Estos últimos serán descritos más adelante.

Figura 1.7. Instrumentos para medir presión en base a diafragmas.

E) Presión de vapor: Los líquidos se evaporan porque las moléculas se escapan de su superficie. Cuando el espacio por encima del líquido está limitado, las moléculas de vapor ejercen una presión parcial en dicho espacio, llamada

26

Capítulo 1: Conceptos generales de fluidos

presión de vapor. Después de un tiempo suficiente, el número de moléculas de vapor que chocan contra la superficie del líquido y de nuevo se condensan es justamente igual al número de las que escapan en un intervalo de tiempo, existiendo un equilibrio. Como este fenómeno depende únicamente de la actividad molecular, la cual es función de la temperatura, la presión de vapor de un fluido dado depende de la temperatura y aumenta con ella. Cuando la presión encima del líquido se iguala a la presión del vapor del líquido, éste hierve. La ebullición del agua, por ejemplo, puede ocurrir a la temperatura ambiente si la presión se reduce suficientemente. Así a 20 ºC el agua tiene una presión de vapor de 2337 (Pa), equivalente a 0.00000176 kgf/cm2. En algunas situaciones que implican el movimiento de líquidos es posible que se produzcan presiones muy bajas en determinados sectores del sistema. Bajo tales circunstancias la presión puede llegar a ser igual o menor que la presión del vapor. Cuando ocurre esto, el líquido se transforma en vapor. Este es el fenómeno de cavitación, en el cual se forman bolsas o cavidades de vapor, que normalmente son sacadas de su punto hacia zonas donde la presión es mayor que la presión del vapor, produciéndose el colapso de estas cavidades o burbujas. Este crecimiento y decaimiento de las burbujas de vapor afecta al rendimiento de funcionamiento de las bombas y turbinas hidráulicas y puede dar como resultado erosiones en las partes metálicas de la región de cavitación (Ver capítulo 5). F) Tensión superficial: En la superficie de contacto entre líquido y gas, parece formarse en el líquido una película o capa especial, debida en apariencia a la atracción de las moléculas del líquido situadas por debajo de la superficie. Esto se comprueba fácilmente al colocar una pequeña aguja en la superficie del agua en reposo y observando como es soportada por la película. Esta propiedad de la película superficial de ejercer una tensión se llama tensión superficial y es la fuerza necesaria para mantener la unidad de longitud de la película en equilibrio. La tensión superficial agua/aire dismi-

27

Escurrimiento de fluidos

USACH

nuye desde 72.8 mili Newton/m (0,00745 kgf/cm) a 20 ºC hasta 58.9 mili Newton/m (0,00596 kgf/m) a 100 ºC.

1.3.

Hidrostática



Un fluido se considera estático si todos los elementos constituyentes

del fluido se encuentran en reposo o se mueven con velocidad uniforme, con respecto a un sistema de referencia. Esto se cumple cuando existe un equilibrio de las fuerzas que actúan sobre el fluido.

En un fluido en reposo sólo actúan esfuerzos normales (presión) y necesariamente los esfuerzos de corte (τ) deben ser cero.

Para un fluido en reposo la presión es la misma en todas las direccio-

nes, pero varía con la altura, según la ecuación fundamental de la hidrostática (conocida también como ecuación de Torricelli), es decir: (1.11) p = presión ρ = densidad del fluido h = altura g = aceleración de gravedad De acuerdo con la Figura 1.8, pa = pb = pc, siempre que entre a, b, y c exista un medio continuo y homogéneo. Conociendo la presión en la aparte inferior del estanque, se puede evaluar la presión en h = h1, integrando la ecuación fundamental de la hidrostática:

28

Capítulo 1: Conceptos generales de fluidos

d p = - ρ⋅g⋅dh

Considerando que el fluido es incompresible (ρ = constante) se obtiene: pa = pd – ρ g h1

(1.12)

Esta última expresión se utiliza para evaluar cambios de presión en estanques, columnas de líquido, etc.

Figura 1.8. Presiones en un estanque de líquido, de densidad ρ.

29

Escurrimiento de fluidos

USACH

Ejemplo 1.3 Deduzca una expresión que permita evaluar la presión en el fondo de un tambor (px) que contiene en su parte inferior agua y en la parte superior aceite: Solución:

p2 = patmosférica



p1 = p2+ ρaceite· g·h2



px = p1 + ρagua· g·h1

Luego: px = p2 + g [ρaceite· h2 + ρagua · h1] Ejemplo 1.4 Un manómetro en U consiste en una columna de vidrio curvada, que contiene en su interior un líquido, inmiscible con el fluido del estanque al cual se conecta. Al respecto, un manómetro en U se encuentra conectado a un gasómetro que contiene 02. Si la presión atmosférica es de 760 mm de Hg, determine la presión absoluta y la presión manométrica en el interior del estanque, considerando que el líquido manométrico es mercurio y que la rama de la derecha está abierta a la atmósfera. Solución: p1 = p2; (puntos a igual nivel) p3 ≠ p4; (están a un mismo nivel, pero entre ellos no existe un medio continuo y homogéneo)

30

Capítulo 1: Conceptos generales de fluidos

Aplicando la ecuación de la hidrostática a la rama izquierda y luego a la rama derecha del manómetro: p1 = pO2 + ρO2 · g · 0.5(m) + ρO2 · g · 0.2(m) p2 = patmosférica + ρHg · g · 0.5(m) + ρaire · g · 0.2(m) Igualando estas expresiones: pO2 = patmosférica + ρHg · g · 0.5(m) + ρaire · g · 0.2(m) – ρO2 · g · 0.5(m) - ρO2 · g · 0.2(m) pO2 = Patmosférica + g [ρHg · 0.5 + ρaire · 0.2 – ρO2 · 0.7] pO2 = 760 mmHg + 9.8 m/s2 [13.600 kg/m3 · 0.5 m + 1.1 kg/m3 · 0.2 m – 1.3 kg/m3 · 0.7] pO2 = 760 mmHg + 66633.2 kg/s2m; pero: 1 kg/s2m = 0.0075 mmHg, luego: pO2 = 1259.75 [mmHg] (presión absoluta) Para obtener la presión manométrica se debe restar a este valor la presión atmosférica local: pmanométrica de O2 = 499.75 mm Hg

Ejemplo 1.5 Un manómetro Bourdon se encuentra conectado a una tubería por la cual circula agua, según se muestra en la figura. Si el manómetro indica una presión de 60 psig, ¿Cuál es la presión absoluta en la cañería si en el lugar la presión atmosférica es de 14 psia (724 mm Hg)?

31

Escurrimiento de fluidos

USACH

Solución:

p1 (manométrica) = 60 psig = 3102.9 mm Hg p1 (abs) = 60 psi + 14 psia = 74 psia = 3826.9 mm Hg p1 = p2 + ρagua· g · 0.61 m

p2 = p1 – ρagua· g · 0.61 = 3826,9 mm Hg –1000 kg/m3 · 9.8 m/s2 · 0.6 m p2 = 3826,9 mm Hg – 5978 kg/ms2 · 0.0075 mm Hg/(kg/ms2) = 3782 mmHg p2 = 73.13 psia (absoluta) = 58.4 psi (manométrica)

32

Capítulo 1: Conceptos generales de fluidos

EJERCICIOS 1.1. A nivel del mar se tiene una presión de 14,7 lbf/pulg2 y una temperatura de 20 ºC. Determine la presión atmosférica a alturas de 300; 3.000 y 30.000 m. sobre el nivel del mar, considerando una atmósfera isotérmica. 1.2. Grafique la viscosidad de aire, agua y ácido sulfúrico al 60%, a una temperatura de 0; 50 y 75 ºC, a la presión atmosférica. ¿Qué conclusiones deduce de estos resultados? 1.3. Utilizando las propiedades críticas del N2, determine su viscosidad a 1; 10; 50 y 100 at. y a una temperatura de 293 K.1.4. Determine la densidad de flujo de cantidad de movimiento (esfuerzo de corte) en estado estacionario, expresado en kgf/m2 y kg/m.s2, cuando la velocidad V de la lámina superior es de 0,3 m/s. La distancia entre las láminas es de 0,3 mm. y la viscosidad del fluido es 0,7 c.p.

1.4. Exprese una presión manométrica de 0,6 kgf/cm2, cuando la presión atmosférica es de 750 mm de Hg, en: a) metros de columna de agua. b) psia. c) kgf/cm2 (abs). d) Newton/m2 (abs). e) Kilo-Pascales (abs). f) torr (abs).

33

Escurrimiento de fluidos

USACH

1.5. El agua que se encuentra en un depósito cerrado está sometida a una presión manométrica de 0,3 kgf/cm2, ejercida por aire comprimido introducido en la parte superior del depósito. En la pared lateral del mismo hay un pequeño orificio, situado 5 m. por debajo del nivel del agua. Calcular la velocidad con la cual sale el agua por este orificio. 1.6. El tubo de la figura está lleno de aceite. Determine la presión en A y B en metros de columna de agua.

1.7. Una manera de determinar la densidad de un líquido, ρx, (con densidad cercana pero menor que la del agua) es el sistema que se muestra en la figura siguiente:

34

Capítulo 1: Conceptos generales de fluidos

1.8. La pérdida de presión a través del accesorio X se mide utilizando un manómetro de la forma indicada, con un aceite cuya densidad relativa es 0,75. El líquido que fluye tiene una densidad relativa de 1,5. a) Utilizando las medidas que se indican, calcular la diferencia de presión entre las dos conexiones piezométricas en unidades MKS y en columna de fluido circulante. b) Determine qué densidad debe tener el líquido manométrico para que una diferencia de presión de 0,1 atmósferas produzca una diferencia de 100 cm.

35

Capítulo 2: Balances macroscópicos

CAPÍTULO 2 BALANCES MACROSCÓPICOS

2.1.

Ecuación de balance



La expresión general para un balance tiene la forma siguiente: (2.1)

Esta ecuación de balance puede ser aplicada a cualquier propiedad extensiva X, entendiéndose por propiedad extensiva aquella que se duplica al duplicarse la masa del sistema, por ejemplo, energía, entalpía, masa, cantidad de movimiento, etc. También es aplicable a unidades contables: dinero, personas, árboles, etc. Al aplicar la ecuación de balance deben especificarse claramente los límites del sistema sobre el cual se aplica el balance. Los tipos de sistema son: • Un sistema abierto, aquel cuyos límites permiten la entrada y salida de energía y masa. • Un sistema cerrado, aquel cuyos límites sólo permiten la entrada y salida de energía. • Un sistema aislado, aquel que no permite el intercambio de energía ni de masa. Si se elige como sistema alguna región arbitraria del espacio en la cual puedan existir flujos de entrada y/o salida, entonces este sistema será el volumen de control.

37

Escurrimiento de fluidos

USACH

La ecuación de balance no es aplicable a unidades incontables ni a propiedades intensivas (temperaturas, presión, viscosidad, dureza, color, densidad, etc.).

2.2.

Balance de masa



Como ya se señaló, a la masa, propiedad extensiva, se le puede apli-

car la ecuación de balance. Los términos de creación y destrucción de masa son cero, por lo que el balance de masa será:

(2.2)

La ecuación de balance de masa se conoce también con los nombres

de: ecuación de conservación de masa, ecuación de continuidad o principio de continuidad.

Si el término de acumulación de masa se hace despreciable, se habla

de balance de masa en estado estacionario. En este caso se obtiene: {Flujo másico de entrada} = {Flujo másico de salida}

(2.3)

Esta expresión se conoce como ecuación de continuidad en estado estacionario.

38

Capítulo 2: Balances macroscópicos

Ejemplo 2.1 Por una tubería circula agua en estado estacionario. En un sector de la tubería de diámetro 0,1 m la velocidad es de 3 m/seg. Si luego el diámetro de la tubería aumenta a 0,2 m, determine la nueva velocidad del agua. ρagua = 1000 (kg/m3)

V1 = 3 (m/s) D1 = 0.1 (m) D2 = 0.2 (m) V2 = ¿? Solución: El sistema será el fluido comprendido entre 1 y 2. Flujo másico = w1 = ρ1 (kg/m3) · A1 (m2) · V1 (m/s)

Pero, de acuerdo con el balance de masa: w1 = w2 = ρ2 A2 V2, de aquí:

Note que:

w1 = w2 = ρ1 A1 V1 = ρ2 A2 V2

Si el fluido es un líquido como en este caso (fluido incompresible): D  V2 = V1  1   D2 

2

Ejemplo 2.2 Un estanque cilíndrico de diámetro 3 m posee una tubería de entrada de diámetro 10 cm y una tubería de salida de diámetro 15 cm. En la primera entra

39

Escurrimiento de fluidos

USACH

agua a una velocidad de 3 (m/seg) y en la segunda sale agua a una velocidad de 1 (m/seg) ¿El nivel de agua en el estanque sube o baja? Esta situación corresponde a un balance de masa en estado no estacionario. D1 = 0.1 m V1 = 3 m/s D2 = 0.2 m V2 = 1 m/s

Dado que el fluido es incompresible ⇒ ρ ρ= =ρ1ρ= ρ2ρ 2 1 =

El valor positivo de dv/dt nos indica que el nivel del estanque sube. Ejemplo 2.3 Se desea saber en cuánto tiempo se llenará el estanque del problema anterior, considerando que cuando se abren las válvulas de entrada y salida contiene 5 m3 de agua. La altura del estanque es de 2 metros. Solución: El planteamiento del problema es similar al anterior, vale decir:

40

; luego

Capítulo 2: Balances macroscópicos

vinicial = 5 m3, para t = 0 vfinal = (π/4) D2 L = (π/4) · 32 · 2 = 14.137 m3 vfinal - vinicial = 0.00589 · t ⇒ t t==1551.3(s) 1551.3 ( s )==25.9 25 .9(min) (min) Ejemplo 2.4 El estanque de la figura contiene aire que puede ser considerado como un gas ideal. La bomba de vacío está extrayendo aire del estanque. Durante el proceso los serpentines mantienen la temperatura constante en el interior del estanque en 70 ºF. El flujo volumétrico de salida de aire es de 1 pie3/min, independiente de la presión. ¿Cuánto tiempo emplea la bomba en disminuir la presión desde 1 atm a 0.001 atm? Solución: El estanque es el sistema en estudio:

Pi = 1 at Pf = 0.001 at Qsalida = 1 (pie3/min)

La ecuación general es: = wentrada – wsalida

;

pero: wentrada = 0

41

Escurrimiento de fluidos

= – wsalida

USACH

;

wsalida = Q · ρ

d ( ρ v) = −–QQ⋅ ρ · ρ ; ;v = cons v = constante tan te dt

2.3.

Balance de energía



Un sistema puede poseer diversas clases de energía: interna, cinética,

potencial, electrostática, magnética, etc. Definiremos las tres primeras por ser las más importantes en el flujo de fluidos: •

• •

Energía interna (U). Es una propiedad intrínseca del fluido. El sistema a estudiar estará integrado por moléculas orientadas según una geometría particular, en el caso de sólidos, o bien siguiendo un movimiento errático en el de los fluidos. Energía potencial (Φ). Se debe a la posición del fluido con respecto a un plano de referencia arbitrario. Energía cinética (K). Está asociada al movimiento del fluido.



Para estas dos últimas existen expresiones que permiten su evaluación. El sistema a considerar será abierto, siendo el sistema el fluido que pasa por el circuito general mostrado a continuación:

42

Capítulo 2: Balances macroscópicos

Figura 2.1. Circuito sobre el cual se realizará el balance de energía.

Consideremos que el sistema puede almacenar energía interna, potencial y cinética y que pueden existir además energía “en tránsito” como calor (Q) y trabajo (W).

Usaremos la siguiente nomenclatura:

Energía (Energía/tiempo) (Energía/masa) uˆ Energía interna U u ˆ Energía potencial Φ φ φ kˆ Energía cinética K k Aceptaremos además la siguiente convención de signos:

Figura 2.2. Convención de signos.

43

Escurrimiento de fluidos

USACH

En el sistema elegido puede entrar o salir cualquiera de los tipos acumulables de energía (interna, potencial, cinética), al igual que las energías en tránsito (calor y trabajo). Usando la ecuación general de balance, aplicada al sistema señalado en Figura 2.1, se obtiene:

(2.4)



En este balance no se han considerado términos de “generación” o de

“destrucción” de energía, ya que normalmente no se producen en situaciones corrientes de escurrimiento de fluidos. A continuación, plantearemos expresiones para cada uno de los términos que aparecen en la ecuación de balance, y luego los introduciremos en ella:

φˆ = g⋅Z g = aceleración de gravedad Z = altura respecto a un nivel de referencia uˆ = (energía interna/masa) kˆ = ½ ⋅V2/α (energía cinética/masa) V = velocidad media en la tubería α = parámetro de corrección de la energía cinética. Se evalúa con Figura 2.3 para fluidos newtonianos. Con respecto al término de trabajo, en un sistema como el mostrado en Figura 2.1 existen dos formas de trabajo: a) Un trabajo que debe realizarse para que la masa de fluido entre y salga del sistema. Se le llama trabajo de flujo. En unidades de energía/masa, el trabajo de flujo = p/ρ. b) Un trabajo relacionado con partes móviles (bomba, compresor, turbina, etc.). Se denomina W en unidades de energía/tiempo y se llama W cuando está en unidades de energía/masa.

44

Capítulo 2: Balances macroscópicos

Figura 2.3. Factor de corrección de energía cinética (α) vs Reynolds (Adaptado de Hicks T. Pump, selection and application, 1957).

Para fluidos no-newtonianos que siguen la ecuación 2.5, α puede ser evaluado con expresión 2.6, válida en régimen laminar. (2.5) (2.6) Introduciendo cada uno de los términos en la expresión general de balance de energía, para un fluido incompresible, se obtiene: (2.7) Si el en sistema no hay acumulación de masa, podemos considerar w1 = w2 = w. Además si consideramos que el balance se realiza en estado estacionario, el término de acumulación de energía será cero. Luego, dividiendo por w, se obtiene: (2.8)

45

Escurrimiento de fluidos

USACH

Reordenando esta expresión y multiplicando por -1 se llega a: (2.9) (2.10)

Introduciendo las expresiones para la energía cinética y potencial: (2.11) Reescribiendo esta ecuación se obtiene: (2.12)

En la mayoría de las aplicaciones de interés práctico es posible con-

siderar α = 1 [Régimen turbulento].

De acuerdo con lo especificado anteriormente, el significado de cada

uno de estos términos es: p/ρ = energía de flujo por unidad de masa del fluido gZ = representa la energía potencial sobre un nivel de referencia por unidad de masa del fluido. 1V2 = energía cinética por unidad de masa del fluido. 2 α = es el trabajo realizado por el fluido o el trabajo realizado sobre el fluido. Es distinto al trabajo entregado por los equipos mecánicos, aunque se relacionan a través de la eficiencia de las unidades mecánicas.

46

Capítulo 2: Balances macroscópicos

(∆ uˆ − Qˆ ) = corresponde a las pérdidas por fricción al ir el fluido desde 1 a 2. Llamaremos a (∆ uˆ − Qˆ ) = Eˆ v Siempre se cumplirá que Eˆ v > 0 (de lo contrario significaría que es posible enfriar un fluido mediante fricción).



De acuerdo con esto, la ecuación de balance de energía será: (2.13) (2.14) Pero p/ρ = vˆ p ; siendo vˆ = volumen específico (2.15) Por definición: entalpía: Hˆ = uˆ + pvˆ Luego: ∆Hˆ = ∆ uˆ + ∆( pvˆ)

(2.16)

Introduciendo el término de entalpía en la ecuación de balance: (2.17) Esta ecuación es de gran utilidad para evaluar cambios energéticos en sistemas en que la temperatura varía considerablemente. Se conoce como balance de energía total.

47

Escurrimiento de fluidos

USACH

EJERCICIOS 2.1. Un estanque de agua tiene una tubería de entrada de 30 cm. de diámetro y dos tuberías de salida de 15 cm. y 10 cm., respectivamente. La velocidad en la tubería de entrada es 1,5 m/s y en la salida de 15 cm. de diámetro es 2 m/s. La masa del estanque aumenta a razón de 23,6 kg/s. ¿Cuál es la velocidad, el flujo volumétrico y el flujo másico en la tubería de 10 cm.? ¿Si esta tubería se cierra cuál será la variación de la masa del estanque en el tiempo? 2.2. Gas Metano (M=16 kg/kgmol) circula por una cañería en estado estacionario a 25 ºC. En la sección (1) de la cañería, la velocidad es 2,2 (m/s) y la presión es 1,2 (at). ¿En qué proporción se debe disminuir el diámetro de la cañería para que en el punto (2) la velocidad del gas aumente a 8,8 (m/s), considerando que la presión en esa sección es de 1,1 (atm)? [Solución: D1/D2 = 1.9]

48

Capítulo 2: Balances macroscópicos

2.3. Dos estanques de grandes dimensiones conteniendo agua se encuentran conectados con una tubería. El estanque A tiene una presión sobre el agua de 1,2 (at), mientras que el estanque B se encuentra abierto a la atmósfera. ¿Es posible que físicamente se dé esta situación? Justifíquelo.

49

Capítulo 3: Flujo de fluidos

CAPÍTULO 3 FLUJO DE FLUIDOS

3.1.

Naturaleza del flujo de fluidos

Antes de iniciar las aplicaciones de los balances de masa y energía al flujo de fluidos, es fundamental considerar algunos aspectos del comportamiento de los fluidos. En 1883, Osborne Reynolds (1842-1912), un físico británico, concluyó que a bajas velocidades de agua, esta fluía en láminas o capas paralelas. En efecto, Reynolds en 1883 realizó el siguiente experimento: conectó un depósito de agua a un tubo de vidrio horizontal. En el extremo por donde ingresa la corriente de agua instaló una boquilla por la que se inyecta agua coloreada, tal como se esquematiza en la Figura 3.1. A bajas velocidades, a lo largo del tubo permanece el filamento de tinta, ya que las partículas de tinta difunden lentamente y no tienen tiempo de diseminarse. A este flujo se le llama “laminar”. Reynolds probó disminuir y aumentar la viscosidad del fluido, calentando y enfriando el agua respectivamente. El experimento mostró que en todos los casos existe una velocidad crítica que varía en proporción directa con la viscosidad del flujo. Al aumentar la velocidad del agua encontró que a una velocidad determinada, velocidad crítica, desaparecía el chorro coloreado y la masa global de fluido se coloreaba uniformemente, concluyendo que sobre la velocidad crítica las partículas dejan de moverse en forma ordenada y paralela, moviéndose en forma caótica, mezclándose completamente.

51

Escurrimiento de fluidos

USACH

Figura 3.1. Experimento de Reynolds.

Reynolds determinó que la velocidad crítica dependía del diámetro del tubo y de las propiedades físicas del fluido: densidad (ρ) y viscosidad (µ). Posteriormente, se concluyó que estos parámetros se pueden agrupar en un número adimensional, número de Reynolds, Re, definido de la forma siguiente:

(3.1)

ρ = densidad del fluido. μ = viscosidad del fluido. V = velocidad media del fluido. L = longitud que caracteriza al ducto. Para tuberías circulares L = D (diámetro). Experimentalmente se determinó que para el flujo en el interior de tuberías, cuando Re ≤ 2100 el flujo es laminar y que cuando Re > 4000, el régimen de flujo es turbulento. En la zona comprendida entre 2100 y 4000 el flujo puede ser laminar o turbulento, dependiendo de las características superficiales de la tubería, llamándose esta zona de régimen de transición. En régimen laminar es posible describir totalmente el comportamiento del fluido en el interior de ductos, a partir de ecuaciones teóricas. Por ejemplo, para un fluido newtoniano, se ha demostrado que el perfil de velocidades en

52

Capítulo 3: Flujo de fluidos

el interior de un ducto circular es parabólico, con la velocidad máxima en el centro del ducto, cuando el flujo esta establecido (sin perturbaciones). La expresión para el perfil de velocidades en esta situación es:

(3.2)

en que: v = velocidad puntual o local para cualquier r v máx = velocidad en el centro del ducto

r R

= posición dentro del ducto (radial) = radio del ducto

En problemas de escurrimiento de fluidos, normalmente interesa conocer la velocidad media del fluido, V. Esta puede evaluarse utilizando la definición de velocidad media:

(3.3)

Al introducir en esta definición general de velocidad media, el perfil para régimen laminar de un fluido newtoniano e integrar, se obtiene:

(3.4)

Por otro lado, existen instrumentos que permiten medir la velocidad puntual o local v máx en función del radio. Si este instrumento, por ejemplo un tubo de Pitot, se coloca en el centro del ducto y el régimen es laminar, es posible conocer la velocidad media usando la expresión V = v máx /2. Dado que en general, en régimen turbulento no es posible describir teóricamente el comportamiento del fluido, se han establecido relacio-

53

Escurrimiento de fluidos

USACH

nes empíricas, observándose que en régimen altamente turbulento V/ v máx ~ 0.81. En la Figura 4.6 se presenta un gráfico que relaciona el cuociente V/ v máx con el número de Reynolds.

3.2.

Aplicaciones de la ecuación de Bernoulli, en situaciones que la fricción es despreciable



En la sección 2.3. se obtuvo la ecuación de Bernoulli, expresión que

corresponde a un balance de energía mecánica en estado estacionario:

(2.13)



Considerando que en algunas aplicaciones el término de pérdidas por fricción (Êv) es pequeño frente a los otros términos, se obtiene una expresión simplificada de la ecuación 2.13:

(3.5)

Esta última expresión es aplicable sólo en algunas situaciones particulares: Ejemplo 3.1 Un estanque se encuentra lleno de agua y abierto en el tope. Posee un pequeño orificio cerca del fondo, cuyo diámetro es pequeño comparado con el diámetro del estanque, ¿Cuál es la velocidad del agua a la salida del orificio? Solución: Aplicaremos el balance de energía mecánica (Bernoulli) entre la superficie del estanque (Punto 1) y el fluido que sale por el orificio (Punto 2).

54

Capítulo 3: Flujo de fluidos

Debe observarse que las presiones Consideraciones: absolutas siempre serán positivas, no así las manométricas que podrían ser a. v1 = 0 (Velocidad de descenso del negativas. agua en el estanque). b. P1 = P2 Presión atmosférica. c. Fricción despreciable

∧

⇒

Ev

≈0

∧

d. No hay trabajo externo ⇒ W = 0 e. El flujo es estacionario, es decir el nivel del agua no baja. f. Z2 = 0 (nivel de referencia). g. α2 =1

Luego:

Simplificando:

Introduciendo valores numéricos:

Nota: Si se consideran las pérdidas por fricción debido a la expansión, dada por ÊV = K · ν2/2, con K=1, se obtiene:

Este último valor es más cercano a la realidad.

55

Escurrimiento de fluidos

USACH

Ejemplo 3.2 El mismo estanque del ejemplo anterior se encuentra ahora sumergido en otro estanque de mayor diámetro, que contiene gasolina. (G.E. = 0,72) ¿Cuál será el valor de la velocidad de salida del agua? Solución:

Datos: ρH20 = 1000(kg/m3); ρgasolina = 720(kg/m3) • De acuerdo con lo visto en hidrostática ⇒ p1 = p1* . • Punto 2 se encuentra ubicado a la salida del estanque de agua. p1 = p0 + ρgasolina · g · x p2 = p1 + ρgasolina · g · h • Nótese que p 2 ≠ p 2* , ya que * si p 2 = p 2 , no se produciría la descarga de agua. Además

p 2* > p 2 .

• Z2 = 0 (nivel de referencia)

Aplicando un balance de energía mecánica (Bernoulli) entre 1* y 2 (el sistema es el agua) se obtiene:

56

Capítulo 3: Flujo de fluidos

Introduciendo los valores numéricos:

v2 = 2,34(m/s) Nota: Si se evalúan las pérdidas por fricción debido a la expansión, con la expresión ÊV = K · ν2/2, con K=1, se obtiene:

Este último valor es razonable, ya que la pérdida por expansión reduce la energía disponible para la salida del agua. Ejemplo 3.3 En la figura que se entrega a continuación, se muestran dos estanques conectados a través de una tubería. ¿El agua circula de A a B o a la inversa? Solución: En esta situación debe escribirse la ecuación de Bernoulli con todos sus términos:

Debe suponerse una dirección de flujo y aplicar Bernoulli en esa dirección. Si el valor que se obtiene para ÊV es positivo, la dirección del flujo supuesto es correcta, en caso contrario, la dirección será opuesta.

57

Escurrimiento de fluidos

USACH

Supondremos flujo de 1 a 2:

Luego:

∧

–101,318 (m2/s2) + 9,8 (m/s2) · 8(m) = – EV = 22,9 ( m 2 s 2 ) ∧

2 2 2 EV == 22,9(m 22,9 ( m2/s s) )

Dado que ÊV > 0, la dirección de flujo supuesta es correcta y el fluido va desde el estanque A al B. Ejemplo 3.4 Un flujo volumétrico de 0,25(m3/s) de agua circula a través de una turbina. ¿Cuál es la potencia entregada por el fluido a la turbina, si a la entrada el diámetro es de 30(cm) y la presión de 2,5(atm), mientras que en la salida el diámetro es de 60(cm)y la presión de 0,8(atm)?

58

Capítulo 3: Flujo de fluidos

Solución: Cálculos de las velocidades:

Aplicando un balance de energía mecánica (Bernoulli) entre los puntos 1 y 2:

(m ) Nivel de referencia: z 22 == 0 ⇒ z11 = 11(m)

Introduciendo estos valores en Bernoulli:

Pero: Potencia = Ŵ · w = Ŵ · Q · ρ Potencia = 187,9(m2/s2) · 0,25(m3/s) · 1000(kg/m3) Potencia =

59

Escurrimiento de fluidos

USACH

Ejemplo 3.5 Una bomba centrífuga se utiliza para elevar petróleo Nº 6, desde un depósito subterráneo hasta un estanque ubicado sobre la bomba. Manómetros ubicados a la entrada y salida de la bomba señalan presiones de –10(psig) y 20(psig), respectivamente. Si los diámetros de succión y descarga son iguales y el flujo volumétrico transportado es de 300(l/min), ¿Cuál es la potencia que entrega la bomba al fluido? Solución:

ρpetróleo = 910(kg/m3) patmosférica = 14,696(psi) pabs 1 = (14,696 – 10)(psi) = 4,696(psi) pabs 1 = 4,696(psi) = 0,32(atm) pabs 2 = (14,696 + 20)(psi) = 34,696(psi) pabs 2 = 34,696(psi) = 2,36(atm)

Aplicando Bernoulli entre la entrada y salida de la bomba se obtiene:

Se ha considerado despreciable la diferencia de altura entre succión y descarga de la bomba. Luego:

60

Capítulo 3: Flujo de fluidos

Luego: 227,13(m2/s2) = −Ŵ Ŵ = – 227,13(m2/s2) Potencia = Potencia = Nota 1: El signo negativo indica que la energía es entregada por la bomba al fluido. En el ejemplo 3.4 la potencia era positiva ya que la energía la suministraba el fluido a la turbina. Nota 2: La potencia así calculada corresponde a la energía requerida por el fluido. La potencia “consumida” por el grupo moto-bomba debe ser mayor, a fin de compensar las diferentes pérdidas que se producen: fricción del fluido en el interior de la bomba, roce de las partes mecánicas de la bomba, flujo circulatorio, etc. Estas pérdidas se consideran en la eficiencia o rendimiento del grupo moto-bomba, definido como:

La eficiencia es un parámetro característico de cada bomba, es función del caudal y debería ser entregado por el fabricante del grupo motobomba. Valores corrientes de eficiencia oscilan entre 0,55 y 0,75. Ejemplo 3.6 Freón 12 circula en estado estacionario a través de una válvula reductora de presión. A la entrada de la válvula de presión es de 100(psia) y la tempera-

61

Escurrimiento de fluidos

USACH

tura de 100 ºF. Si la presión a la salida de la válvula es de 20(psia). ¿Cuál es la temperatura a la salida de la válvula? Solución:

El balance de energía mecánica (Bernoulli) entre los puntos 1 y 2 es:

∧  Considerando sistema adiabático  Q = 0  , flujo horizontal y despre  ciando la variación de energía cinética, se obtiene: ∧

∧

H1 = H 2

Desde un diagrama presión v/s entalpía para Freón 12, se obtiene: ∧

(100ººF, ( psia )) ==888,5(Btu/lb) H 1 (100 F ,100 100(psia)) ,5(Btu lb ) Pero, por lo anterior:

∧

(T2, 2020(psia)) ( psia )) = 888,5(Btu/lb) H 2 (T ,5(Btu lb )

Del mismo diagrama citado se obtiene que la temperatura es aproximadamente 75 ºF.

62

Capítulo 3: Flujo de fluidos

3.3.

Definición de factor de fricción

Para evaluar el término de pérdidas por fricción (Êv) que aparece en la ecuación de Bernoulli, es necesario definir y luego evaluar factores de fricción. Se considerará el flujo estacionario de un fluido incompresible (ρ = cte.) que circula por un ducto recto de sección uniforme. El fluido ejerce sobre la superficie interna del ducto una fuerza F, que puede descomponerse en dos: Fs y Fk, definidas de la siguiente forma: Fs = fuerza estática. Es la fuerza que ejerce el fluido, aunque esté en reposo. Fk = fuerza dinámica. Es la fuerza relacionada con el comportamiento cinético del fluido. Tiene la misma dirección que la velocidad media V en el ducto. El valor de la fuerza Fk puede expresarse como el producto entre un área característica A, una energía cinética característica por unidad de volumen K y un número adimensional f, denominado factor de fricción:

FK = A · K · f

(3.6)

Esta expresión no es una ley de mecánica de fluidos, sino una definición de f. Es evidente, que para un determinado sistema de flujo, f no está definido mientras no se especifiquen A y K. Esta definición es general, válida incluso para la situación en que un fluido circula alrededor de un objeto sumergido. Para el flujo en ductos A corresponde a la superficie mojada y K representa la energía cinética por unidad de volumen, dada por ½ ·ρ · V2. Para tubos circulares de radio R y longitud L, f está definido por:

FK = (2πRL) · (½ · ρ · V2) · f 63

(3.7)

Escurrimiento de fluidos

USACH

Dado que generalmente lo que se mide no es FK sino que la caída de presión y la diferencia de altura, se buscará una expresión para FK aplicando un balance de fuerzas al fluido que circula por la tubería mostrada en la Figura 3.2, entre 1 y 2:

Figura 3.2. Esquema para la aplicación de fuerzas.

En estado estacionario se tiene que:





π R2 · p1 – π R2 · p2 – FK – π R2 · L · ρ · g · senα = 0

Pero: sen α =

(3.8)

(3.9)

(3.10)

Z 2 − Z1 L



π R2 · p1 – π R2 · p2 – FK – π R2 · ρ · g · (Z2– Z1) = 0

(3.11)



FK = π R2 · [(p1 – p2) – ρ · g · (Z2– Z1)]

(3.12)

64

Capítulo 3: Flujo de fluidos

Utilizando el concepto de potencial fluido-dinámico: P = p + ρ · g · Z, se obtiene: FK = π R2 (–ΔP) (3.13) Reemplazando esta última expresión en la definición de factor de fricción, se obtiene: (3.14)

Despejando f:



(3.15)

Este factor así definido, se denomina factor de fricción de Fanning. En algunos textos de mecánica de fluidos en vez de utilizar el factor de Fanning, definen un factor f* que equivale a 4f, denominado factor de Darcy. Un análisis del factor de fricción muestra que éste depende de características, tanto del fluido como del sistema de escurrimiento. Planteando un balance de energía mecánica al sistema mostrado en Figura 3.2, se obtiene:

(3.16)

Introduciendo la definición de potencial fluido-dinámico:

(3.17)

Reemplazando (–ΔP/ρ) por Êv, en la definición de factor de fricción, se obtiene: (3.18)

65

Escurrimiento de fluidos

USACH

Esta última expresión es de gran utilidad ya que permite evaluar las pérdidas por fricción, Êv, si se conoce el valor de f. Si se utiliza el factor de Darcy (f*), entonces:

(3.19)

La técnica de análisis dimensional, en unión con lo observado experimentalmente, permite establecer de cuántos y cuáles grupos adimensionales depende cualquier parámetro de naturaleza física. Para el factor de fricción se encontró que f = f(Re, L/D, ε/D), donde ε es la rugosidad del material, definida como el promedio de altura de las irregularidades de la superficie del material del ducto. En los sistemas de flujo, habitualmente se considera que ε es un parámetro constante. Sin embargo, en ductos de grandes longitudes (L >> 1000), pequeñas variaciones en ε debido a corrosión y/o depósito de incrustaciones a lo largo del tiempo, pueden afectar significativamente los requerimientos de potencia para una determinada exigencia de flujo, por su efecto en el factor de fricción.

ε/D se conoce con el nombre de rugosidad relativa. En el apéndice

se muestra un gráfico de ε/D vs D para diferentes materiales de ductos. Para ductos en que L/D >> 1 (como generalmente ocurre en las situaciones prácticas), entonces el factor de fricción se hace independiente de L/D, luego: f = f (Re, ε/D). Valores experimentales de f vs Re para diversos tipos de materiales y por ende diferentes ε/D, fueron graficados por Moody en un gráfico que lleva su nombre, obteniéndose resultados satisfactorios (Figura 3.3).

66

Capítulo 3: Flujo de fluidos

Figura 3.3. Gráfico de Moody. Factor de fricción de Fanning, f, versus rugosidad relativa, ε/D.

Con respecto a este gráfico hay que señalar que aparece una recta para f (independiente de la rugosidad relativa) para la zona de flujo laminar. La ecuación de f en esta zona es:

f = 16/Re

(3.20)

Pasada una cierta zona de número de Reynolds, las curvas se hacen paralelas al eje de Reynolds, independizándose de éste, dependiendo sólo de ε/D. Se habla entonces de régimen altamente turbulento. Para las zonas de transición y de régimen turbulento, en la literatura se encuentran varias ecuaciones, que ajustan en diferentes rangos de Reynolds y condiciones de ε/D, los factores de fricción del gráfico de Moody. En general, a medida que se complican estas ecuaciones, aumenta la exactitud de sus predicciones.

67

Escurrimiento de fluidos



USACH

Una de las ecuaciones más sencillas indica que:



(3.21)

la cual tiene un error máximo de 0,75% para Re > 30000 y ε/D > 0,004.

Una ecuación más general es la ecuación de Shacham (Ecuación 3.22), la que presenta un error máximo de 1%, respecto a las curvas del gráfico de Moody. Esta ecuación es válida para Re > 4000 y 0.005 < 4f < 0.08

3.4.

(3.22)

Descripción de cañerías, válvulas y accesorios



Las tuberías pueden fabricarse con cualquier material de construcción disponible, dependiendo de las propiedades corrosivas del fluido que se maneja, su temperatura y presión. Entre los materiales se incluyen aceros, cobre, polímeros, vidrio, concreto, etc. Sin embargo, los materiales de tubería más comunes en la industria son el acero (varias calidades), el cobre y el bronce, seleccionándose el más adecuado de acuerdo con la aplicación específica y los costos involucrados. Dado que los tubos se fabrican de diversos diámetros y espesores de pared, existe una normalización establecida por la American Standards Associations (ASA), la que establece las características de las dimensiones de los tubos.

En el caso específico de los tubos de acero (común o comercial), se

ha establecido que el tamaño de los tubos y de las conexiones asociadas se realice en función del diámetro nominal y del espesor de pared. Por lo tanto, el diámetro nominal no corresponde para los tubos de acero, ni al diámetro

68

Capítulo 3: Flujo de fluidos

exterior ni interior de los tubos. El espesor de pared se expresa por el número de cédula, el cual es un cociente entre la presión interna y la tensión permisible (# cédula ∼ 1000 ∙ presión interna/tensión permisible). Se utilizan diez números de cédula, a saber: 10, 20, 30, 40, 60, 80, 100, 140 y 160. El espesor de pared aumenta con el número de cédula. Para tubos de acero comercial, la cédula 40 corresponde al tubo “normal”, para emplear en aplicaciones sin mayores exigencias. En las instalaciones es imprescindible el uso de diversas conexiones para trasladar el fluido de un sector a otro. Una conexión cumple el papel de: a) Juntar dos tuberías [coplas, unión americana] b) Cambiar la dirección de la tubería [codos] c) Cambiar la sección de flujo [reducciones] d) Terminar la tubería [tapones] e) Unir o diversificar una corriente [tees, cruces e yes] f) Control de flujo [válvulas]

Figura 3.4. Conexiones.

69

Escurrimiento de fluidos

USACH

Las válvulas son conexiones que cumplen diversas funciones en un circuito de escurrimiento de fluidos. Las válvulas se utilizan para regular el flujo o bien cerrar el paso completamente. Entre las válvulas de mayor uso se tienen las de compuerta, de globo y en los últimos años, la de bola. La válvula de compuerta consiste en un disco que se desliza perpendicularmente al flujo. Su uso principal es para sellar o detener el flujo en forma rápida, ya que pequeñas variaciones en la altura del disco se traducen en grandes cambios en el área disponible al flujo. La válvula de bola, ampliamente utilizada tanto a nivel doméstico como industrial, sirve para los mismos fines que la válvula de compuerta. Su principal característica es que cuando está 100% abierta, prácticamente no produce ninguna obstrucción al paso del fluido.

La válvula de globo por su diseño es más adecuada para regular el paso de fluido. En esta válvula, el fluido pasa a través de una abertura cuya área se controla mediante un disco colocado en forma casi paralela a la dirección del flujo.

Figura 3.5. Válvulas. a) Compuerta, b y c) Bola.

70

Capítulo 3: Flujo de fluidos

Figura 3.6. Válvula de globo: 1- Disco que interrumpe el paso del fluido. 2- Eje o husillo (conduce y fija el obturador). 3- Asiento: Sector de la válvula donde se realiza el cierre con el disco. 4- Empaquetadura del eje. 5- Juntas de cierre. 6- Cuerpo de la válvula. 7- Extremos de la válvula que permiten la conexión a la tubería. 8- Pernos de unión. 9- Manilla de accionamiento.

3.5.

Evaluación de pérdidas en válvulas y accesorios

Las pérdidas por fricción, provocadas por conexiones (válvulas y accesorios) pueden ser determinadas inicialmente en forma experimental, para obtener un valor característico para cada conexión en particular. A continuación se describe el procedimiento para determinar la pérdida provocada por una válvula.

Figura 3.7. Esquema de un sistema para determinar experimentalmente pérdidas en una válvula.

71

Escurrimiento de fluidos

USACH

De acuerdo con lo presentado en la sección anterior, las pérdidas en los tramos de cañería recta están dados por:

(3.23)



(3.24)

Si D1 = D2, entonces V1= V2 y f1 = f2, luego:

(3.25)



(3.26)

ÊV total se obtiene al aplicar Bernouilli entre 1 y 4:

(3.27)

Las presiones p4 y p1 se leen en los manómetros ubicados en 1 y 4. Luego, las pérdidas por fricción en la válvula quedan expresadas por:

(3.28)

De esta forma es posible evaluar, experimentalmente, las pérdidas provocadas por cada uno de los accesorios utilizados en las diferentes instalaciones de redes de flujo de fluidos. Los valores de ÊV de conexiones, obtenidos experimentalmente, pueden ser presentados en dos formas:

72

Capítulo 3: Flujo de fluidos

a) Método del coeficiente de resistencia, K En este método las pérdidas son presentadas en función de un parámetro K, característico de cada accesorio, el cual se considera constante, independiente del régimen de flujo. Entonces:

(3.29)

El método del coeficiente de resistencia K, se utiliza también para evaluar las pérdidas por expansiones o contracciones producidas al cambiar de diámetro una tubería o en entradas y salidas de estanques. En este caso, la velocidad V corresponde a la velocidad en la sección de menor área. b) Método de longitud equivalente En este método se caracteriza la fricción de un accesorio por una longitud de tubería ficticia, la que produciría la misma fricción que el accesorio. La fricción provocada por esta longitud de tubería ficticia se evalúa como ya se presento anteriormente, es decir:

(3.30)

Aunque el término (L/D) es adimensional, es costumbre denominarlo longitud equivalente.

73

Escurrimiento de fluidos

USACH

Figura 3.8. Equipo para determinar pérdidas en accesorios (LOPU/DIQ/USACH).

3.6.

Aplicaciones que involucran la evaluación de factores de fricción

Una vez que se han planteado las ecuaciones válidas para un determinado sistema de flujo, los problemas pueden ser resueltos en forma directa (tipo 1), o se requerirá iterar (tipo 2), aunque esta iteración puede ser obviada si se utilizan los llamados gráficos de Von Kárman. Finalmente, si no se puede evitar la iteración, se tiene un problema tipo 3. La tabla siguiente muestra esta clasificación: Tipo de problema

Datos

Incógnita

1

D, ε, μ, ρ, Q

Êv, ó Ŵ, ó potencia

2

D, ε, μ, ρ, Êv

Q, ó w, ó V

3

ε, μ, ρ, Q, Êv

D

Problemas tipo 1: La solución de este tipo de problema es directa, no se requiere iterar, como se muestra en el siguiente ejemplo. 74

Capítulo 3: Flujo de fluidos

Ejemplo 3.7 Dos estanques de agua son conectados mediante 200(m)de cañería de 3" de acero comercial tipo 40. Deben transportarse 12,5(l/s) desde uno a otro estanque. Los niveles en los estanques son los mismos y ambos están abiertos a la atmósfera. En las pérdidas deben incluirse un codo de 90º, una válvula de globo totalmente abierta, una expansión y una contracción. ¿Cuál es la potencia de la bomba requerida? Solución:



Aplicando un balance de energía mecánica (Bernoulli) entre 1 y 2:

Considerando que: p2= p1 = presión atmosférica z2= z1 v2= v1 = 0

75

Escurrimiento de fluidos

USACH

Se obtiene que:

∧

∧

W = − EV ∧

Determinar el trabajo, W y luego la potencia de la bomba, debe evaluarse en primer lugar el término de pérdidas por fricción, ÊV:

De las tablas del apéndice se obtienen los siguientes valores: (L/D)válvula = 340; (L/D)codo = 30; Kexpansión = 1,0; Kcontracción = 0,5 ρagua = 1000(kg/m3) μagua = 1(cp) = 0,001(kg/m · s) di = 3,068(pulg) = 0,078(m); ε/D = 0,0006 Para evaluar el factor de fricción f ó f * = 4 f, debe evaluarse previamente el número de Reynolds:

76

; Pero:

Capítulo 3: Flujo de fluidos

Luego: Re = 2,04 ·105 del gráfico de Moody se lee f = 0,0049 Introduciendo estos valores en (*) se obtiene:

Entonces: Potencia = Potencia = –202(m2/s2) · 0,0125(m3/s) · 1000(kg/3) Potencia = Potencia = 3,4(hp) Problemas tipo 2: En este tipo de problema, en que se conocen las pérdidas por fricción, pero se desconoce la velocidad del fluido, la solución al problema se obtiene mediante una iteración. Esta última puede evitarse si se utiliza el gráfico de Von Kárman, en el cual se representa (4 f)-0,5 v/s Re · (4 f)0,5, donde:

(A)



(B)

El parámetro Re · (4 f)0,5 se conoce como el número de Kárman y puede evaluarse sin conocer v. Con este gráfico se evita la iteración que sería necesaria si se emplea el gráfico de Moody. El procedimiento es el siguiente:

77

Escurrimiento de fluidos

USACH

De la ecuación (B) se calcula Re · (4 f)-0,5, con ε/D se lee 1 de Von Kárman y de la ecuación (A) se despeja v.

4 f en el gráfico

Ejemplo 3.8 Petróleo a 70 ºF es transportado desde un lugar A al otro B, a través de 4000(ft) de tubería de diámetro interno de 6" y ε = 0,0002(ft). El punto B está ubicado a 50,5(ft)sobre el punto A y la presión en A y B son de 123(psi) y 48,6(psi) respectivamente. ¿Cuál es el flujo volumétrico, Q, del petróleo? Solución:

Datos: (μ/ρ) = υ = 4,12 · 10-5 (ft2/s) ρ = 854(kg/m3) = 53,3(lb/ft3) di = 6(in) = 0,5(ft) ε/D = 0,0004 Aplicando Bernoulli entre (1) y (2):

78

(*)

Capítulo 3: Flujo de fluidos

Considerando que: v1 = v2 z1 = 0; z2 = 50,5(ft)

Reemplazando estos valores en (*):

4f · v2 = 1,211 Una forma de resolver esta ecuación sería suponerse un valor de f, calcular v y con esta velocidad comprobar el valor de f supuesto, previo cálculo de Re. Otra forma de resolver es utilizar el método de Von Kárman:

Leyendo en el gráfico de Von Kárman, con ε/D = 0,0004 se obtiene que 1 4 f = 7 , pero:

Problemas tipo 3: Como se muestra en el ejemplo 3.9, en este tipo de problemas, es inevitable la iteración.

79

Escurrimiento de fluidos

USACH

Ejemplo 3.9 Debe transportarse agua desde un estanque abierto a la atmósfera a través de 200(ft) de longitud equivalente (cañería recta más codos, expansión y contracción), para ser descargada a la atmósfera en un punto 12(ft) sobre el estanque. ¿Cuál es el diámetro mínimo de cañería requerido para asegurar un flujo de 200(gal/min), si la bomba es de 2(hp), con una eficiencia de 60,7%? Solución:

Aplicando un balance de energía mecánica entre 1 y 2:

Reemplazando valores:

80

Capítulo 3: Flujo de fluidos

Pero:

Luego:

(*)



Como puede observarse, la ecuación (*) debe resolverse en forma iterativa: supondremos un diámetro, calcularemos Re, se leerá f y se comprobará la igualdad en la ecuación (*). A continuación se muestran algunos valores de las iteraciones: Dno min al

di (ft)

Re

ε/D

4f

Igualdad en (*)

2"

0,1723

3 · 105

0,0009

0,021

4284

3"

0,256

2 · 105

0,0006

0,019

520,9

3½"

0,296

1,8 · 105

0,0005

0,0185

242

El valor del diámetro que “hace cumplir” la igualdad en (*) está entre valores de diámetro interno de 0,256(ft) y 0,296(ft), vale decir entre diámetros nominales de 3" y 3½". Dado que no existen cañerías entre 3" y 3½" se seleccionará la tubería de 3½" para la instalación propuesta.

81

Escurrimiento de fluidos

3.7.

USACH

Evaluación de pérdidas en ductos no circulares

Cuando la sección del ducto no es circular, o cuando el fluido no llena totalmente la tubería, se debe utilizar un parámetro empírico, que se ha verificado entrega buenas predicciones cuando el régimen es turbulento. Es llamado radio hidráulico, RH. Este RH debe relacionarse con el diámetro de un ducto circular, a fin de emplear las fórmulas habituales de ductos efectivamente circulares. Para este fin se calculará el RH de un ducto circular, obteniéndose:

Esto último es lo que se emplea como diámetro equivalente, Deq, es decir: Deq = 4 ∙ RH. Luego, una vez calculado Deq, el problema se resuelve igual que como si fuera un ducto circular, reemplazando en las ecuaciones que se requiera el diámetro del ducto D por Deq. Ejemplo 3.10 Determine el diámetro equivalente de un canal de regadío, abierto a la atmósfera, con una geometría y dimensiones señaladas en la figura:

82

Capítulo 3: Flujo de fluidos

Solución:

Ejemplo 3.11 Vapor saturado circula por el ducto anular de un intercambiador de calor concéntrico. Determine el Deq, que se empleará para evaluar la caída de presión del fluido que circula por la sección anular.

83

Escurrimiento de fluidos

3.8.

USACH

Diámetro óptimo económico (DOE)

Con las ecuaciones anteriormente planteadas (balances de masa y energía), es posible resolver numerosas situaciones de flujo de fluidos, aunque en aplicaciones de diseño se requiere una ecuación o información adicional. Por ejemplo, consideremos la necesidad de especificar el diámetro de la tubería (D), junto con la potencia de la bomba necesaria para transportar un flujo Q de agua desde un estanque a otro, según muestra la Figura 3.9:

Figura 3.9. Transporte de agua desde un estanque a otro.

Para esta situación, al plantear Bernoulli entre 1 y 2, se observa que la potencia de la bomba requerida depende del diámetro de la tubería, existiendo infinitas soluciones. A medida que se disminuye D, aumenta el consumo de energía (potencia de “bombeo”), para un mismo flujo Q:

84

Capítulo 3: Flujo de fluidos

En situaciones de este tipo la elección del diámetro se decide por alguna de las consideraciones siguientes: a) Stock en bodega de cañerías de un determinado tamaño. b) Consideraciones de espacio disponible para el paso de la tubería. c) Consideraciones económicas. En esta última consideración, los costos a tomar en cuenta son: i) Costos de energía para el transporte del fluido. ii) Costos de mantención de bomba, tubería y conexiones. iii) Costos de inversión e instalación de bomba y tubería.

Al disponer de ecuaciones para cada uno de estos costos, el diámetro

óptimo económico (DOE) corresponde al diámetro que se obtiene al derivar la expresión de costos totales, con respecto al diámetro de la tubería, despejándose de la ecuación resultante el diámetro:



Si bien es cierto en algunos textos se presentan expresiones analíticas para cada uno de estos costos, el procedimiento aceptado es trabajar con valores de DOE recomendados, en función del caudal y densidad del fluido.

85

Escurrimiento de fluidos

USACH

Figura 3.10. Diámetro económico en función del flujo volumétrico y de la densidad del fluido, válido para tuberías de acero comercial, cédula 40 (Adaptado de Perry y Chilton, Manual del Ingeniero Químico,1973).

86

Capítulo 3: Flujo de fluidos

En vez de trabajar con el concepto de diámetro económico, se puede emplear el concepto de velocidad económica, Veco, dada la relación entre ambos términos:

En la tabla siguiente se entregan algunos valores de velocidad económica en función de la densidad del fluido: Tabla 3.1. Velocidades económicas en función de la densidad del fluido.

Veco (m/s)

1.7

1.9

3.1

6.0

11.9

24.0

ρ (kg/m3)

1600

800

160

16

1.6

0.16

Ejemplo 3.12 Se desea transportar 200(gal/min) de agua a 60 ºF, a través de una tubería de acero comercial tipo 40 de 5000(ft) de longitud. ¿Qué diámetro de tubería especifica? ¿Cuál es la potencia que debería entregar la bomba al fluido? Solución:

Q = 200(gal/min) = 0,446(ft3/s) ρH2O = 62,4(lb/ft3) μH2O = 7,4 · 10-4(lb/ft · s)

87

Escurrimiento de fluidos

USACH

Aplicando Bernoulli entre 1 y 2:

De la Figura 3.8 con Q = 200(gpm) y ρ = 62,4(lb/ft3) se lee: di = 3,5"

ε/D = 0,0006, usando gráfico de Moody ⇒ f ==0,0047 0,0047 Introduciendo valores se obtiene ÊV = 7169(ft2/s2) ∧

(

2 2 2 Luego W == –7169(ft −7169 ft2/s )s

3.9.

)

Situaciones complejas

En esta sección se presentarán, a través de ejemplos, diversas situaciones que involucran la utilización del balance de energía mecánica (Bernoulli).

88

Capítulo 3: Flujo de fluidos

Ejemplo 3.13 En la figura se muestra un sifón por el cual circula agua. La cañería tiene un diámetro interno de 3(cm) y su longitud es de 16 m. Los codos instalados son estándar de 90º. La rugosidad relativa puede considerarse en 0,0005. Suponiendo que el nivel del estanque permanece constante y que la temperatura del agua es de 20 ºC, ¿Cuál es el flujo volumétrico Q, que circula por la tubería? Solución: Q=? di = 0,03(m) L = 15(m)

ε/D = 0,0005 f = ≈0,35 0047 Kcontracción = 0,78 ⇒ L/D f = ≈0,45,6 0047 Kexp ansión = 1,0 ⇒ L/D

Aplicando un balance de energía mecánica entre (1) y (2) se obtiene:

Considerando: p2 = p1 =1(atm) ,0047 z2 = 0 ⇒ zf1 = 01,25(m) v1 = 0

89

Escurrimiento de fluidos

USACH

Introduciendo los valores en el balance de energía se obtiene:

Dado que f es función de v2, esta última ecuación debe resolverse en forma iterativa: Suponerse un f → calcular v2 → calcular Re y verificar el valor supuesto de f. Considerando régimen turbulento ⇒ ff==0,00425, 0,00425con este valor se ob4 tiene v2 = 1,46(m/s) y Re = 4,4 · 10 . Del gráfico de Moody se lee f = 0,006. Con este nuevo valor de f se obtiene v2 = 1,19(m/s) y Re = 3,6 · 104. Del gráfico de Moody se obtiene f = 0,006, el cual coincide con el valor supuesto. Para este valor:

Ejemplo 3.14 Un depósito cilíndrico de 1(m) de diámetro y 4(m) de altura está lleno de agua a 20 ºC. El fondo del depósito está conectado a un tubo de 1,5" y 5" de longitud, a través del cual se vacía. ¿Cuál es el tiempo que tarda en descender 1(m) el nivel del agua en el depósito?

90

Capítulo 3: Flujo de fluidos

Solución:

p1 = p2 = 1(atm) v1 = 0 ; (Incluye el Kcontracción= 0,6(m))

Dado que el nivel del estanque disminuye, disminuye también la presión a la entrada del tubo de descarga, por lo que la velocidad del agua a través del tubo varía en función del tiempo. Considerando un punto del depósito a una altura h, al descender el nivel dh en el tiempo dt, el caudal estará dado por:

(1)

En este instante, a través del tubo de sección A2, circulará el mismo caudal: Q = A2 · v2 (2) Tomando como nivel de referencia z2 = 0 y aplicando un balance de energía mecánica (Bernoulli) entre (1) y (2), se obtiene:

91

Escurrimiento de fluidos

USACH

z1 = h = variable



(3)

Igualando (1) y (2), habiendo introducido v2 en (2), se obtiene:

(*) Para integrar la ecuación * supondremos un valor promedio de f, entre las condiciones iniciales y finales:

h = Tiempo de descarga del estanque. El valor de finicial debe evaluarse en forma iterativa, vale decir: Suponerse una vinicial → Calcular Re → Leer f y calcular v2 de la ecuación (3). Si v2 calculado coincide con el supuesto, se tendrá el ffinal. Para el valor de ffinal, el procedimiento iterativo es similar, solo que en la ecuación (3) el valor de h será menor en 1(m) con respecto al considerado en el cálculo de finicial. Con la rutina de cálculo señalada se obtiene 4 finicial = 0,0212 y 4 ffinal = 0,021: ∴ 4 fpromedio 0,0211 promedio==0,0211 Con este valor se obtiene que el tiempo de descarga es de 91(s).

92

Capítulo 3: Flujo de fluidos

Ejemplo 3.15 Por una tubería de 25(cm) de diámetro interno se transporta petróleo a lo largo de una longitud total de 30(km), con un caudal de 1000(m3/día). Con el objeto de aumentar el caudal, conservando las mismas presiones de entrada y salida, se conecta a la tubería primitiva, 5(km) antes del lugar de descarga, otra tubería del mismo diámetro y paralela a la primitiva. Si en las condiciones de transporte la densidad del petróleo es 920(kg/m3) y su viscosidad es de 5(poises), ¿Cuál es el aumento de caudal? Solución:

Calcularemos la fricción en la tubería antes de hacer la conexión adicional:

Una vez realizada la conexión, dado que se conservan las presiones de entrada y salida, “la carga de fricción total ha de ser la misma” que en la situación original.

93

Escurrimiento de fluidos



USACH

La fricción antes de hacer la conexión será:

En que v1 es la nueva velocidad en la tubería, en el tramo 1-2. Si v2 es la velocidad en el tramo final de 5(km), en cualquiera de las 2 tuberías, la fracción en este tramo será:



Por otro lado, el caudal antes de la ramificación ha de ser igual a la

suma de los caudales a lo largo de las tuberías paralelas, es decir:

Luego:

Pero:

Despejando se obtiene v1 =0,276 (m/s) y el caudal será 1171(m3/día). En consecuencia se logra aumentar en un 17,1% al transporte de petróleo. Nota: En tuberías en paralelo las pérdidas de energía mecánica (fricción) son las mismas en cualquiera de las ramas, y no son acumulativas.

94

Capítulo 3: Flujo de fluidos

Ejemplo 3.16 Como se muestra en la figura, desde un recipiente, fluye agua a través de una cañería a un punto de bifurcación, desde donde circula a otros recipientes mediante cañerías separadas. Las cañerías son de acero comercial, catálogo 40. Calcular el flujo en (gal/min) de agua que llega a cada recipiente, suponiendo flujo estacionario. Solución:



Sección 1 = 2000', cañería de 6" Sección 2 = 2000', cañería de 3" Sección 3 = 1000', cañería de 4"



Para grandes líneas, el término cinético en la ecuación de energía mecánica (Bernoulli) se puede despreciar (comprobando al final su efecto), y las presiones manométricas en las superficies libres son cero. Por lo tanto, las expresiones para la ecuación de Bernoulli, para las tres secciones de tubería son:

95

Escurrimiento de fluidos

USACH



(1)



(2)



(3)



Sumando las ecuaciones (1) y (2) y luego las ecuaciones (1) y (3) se

obtienen:

Por otro lado, de un balance de masa se obtiene que:

v1 = v2 ·(D2/D1)2 + v3 ·(D3/D1)2

(4) (5)

(6)

Despejando v2 de la ecuación (4) y v3 de la ecuación (5) se obtiene:



(7)



(8)

Para evaluar los factores de fricción se puede emplear el gráfico de Moody, o bien alguna correlación empírica, como la de Shacham:

96

Capítulo 3: Flujo de fluidos

En el siguiente diagrama de bloques se presenta un esquema de resolución para evaluar Q1, Q2 y Q3. Los valores iniciales necesarios (v2; v3), pueden estimarse empleando el concepto de diámetro óptimo económico (DOE).

Los valores finales son: Q1 = 285(gal/min) Q2 = 60(gal/min) Q3 = 225(gal/min)

97

Escurrimiento de fluidos

USACH

Diagrama de bloques para resolver ejemplo 3.10

98

Capítulo 3: Flujo de fluidos

3.10. Fluidos no newtonianos Para fluidos no newtonianos sigue siendo válida la ecuación 2.13 (Ecuación de Bernoulli), con la salvedad de que α debe ser evaluado con la expresión correspondiente. Para fluidos que siguen la ley de la potencia (o modelo de Ostwald de Waele), se dispone de la ecuación 2.6, válida en régimen laminar:

(2.6)

Si el régimen de flujo es turbulento, se considera α = 1.

El factor de fricción, f, necesario para evaluar Êv, debe leerse de figuras de f vs Reynolds. Por ejemplo, para fluidos cuyo comportamiento se puede representar por la ley de la potencia, debe usarse la Figura 3.11. Para fluidos que siguen el comportamiento de Bingham (ec. 1.6) se debe utilizar la Figura 3.12, la cual lleva como parámetro el número adimensional de Hedstrom (He).

99

Escurrimiento de fluidos

USACH

Figura 3.11. Gráfico de f vs Re’, para líquidos que siguen la ley de la potencia (Adaptado de Levenspiel, Flujo de Fluidos e Intercambio de Calor, 1993).

Si el fluido sigue el comportamiento de Bingham, entonces debe usarse la Figura 3.12:

Figura 3.12. Gráfico de f vs Re’, para líquidos de Bingham [ [

],

] (Adaptado de Levenspiel, Flujo de Fluidos e Intercambio de Calor, 1993).

100

Capítulo 3: Flujo de fluidos

Ejemplo 3.17 En la figura se muestra un sistema de tuberías, a través del cual se transporta, desde un estanque cocedor de diámetro 0,5(m) compota de manzana hasta el sector de envases, la que sigue la ley de la potencia (n = 0,4; K = 0,6(kg/m · s2-n)). Determine el tiempo que demora disminuir desde un nivel inicial de 1(m) a un nivel final de 0,1(m). La longitud de la tubería de acero inoxidable es de 2,75(m) y su diámetro interno es de 25 (mm). El codo indicado es de radio largo. La rugosidad del acero es ε = 0,001(mm).

Considerando la ecuación de Bernoulli:

Igualando los caudales de (1) y (2) tenemos:

101

Escurrimiento de fluidos

USACH

Reemplazando el factor de fricción en la ecuación de Bernoulli se tiene:

Ordenando y despejando en función de v2 se tiene:

Para Regen se tiene que:

FINAL

INICIAL

Iterando: f =gen0,=00425 Si v2 = 1(m/s) ⇒ Re 1168 ⇒ f =gen0,=00425 Si v2 = 2,4(m/s) ⇒ Re 4740 f =gen0,=00425 Si v2 = 3,2(m/s) ⇒ Re 7510 f =gen0,=00425 Si v2 = 2,5(m/s) ⇒ Re 5060 Si v2 = 0,5(m/s) Si v2 = 0,9(m/s) Si v2 = 1,5(m/s) Si v2 = 1,9(m/s)

⇒ ⇒ ⇒ ⇒

0,00425= 0,0137 ⇒ vf2==1,89(m/s) 0,00425 f ==16/1168 ⇒ f = 0,0054 0,00425⇒ vf2 == 3,2(m/s) 0,00425 ⇒ f ==0,0046 0,00425⇒ vf2 == 2,6(m/s) 0,00425 ⇒ f ==0,0050 0,00425⇒ vf2 == 2,56(m/s) 0,00425

f =gen0,=00425 0,00425⇒ vf2 == 0,94(m/s) 0,00425 Re 385 ⇒ f ==0,0415 f =gen0,=00425 0,00425 ⇒ vf2 == 1,38(m/s) 0,00425 Re 986,8 ⇒ f ==0,016 f =gen0,=00425 0,00425 ⇒ vf2 == 1,82(m/s) 0,00425 Re 2235 ⇒ f ==0,007 f =gen0,=00425 0,00425⇒ vf2 == 1,9(m/s) 0,00425 Re 3262 ⇒ f ==0,0060



;

Por otro lado, de un balance de masa (continuidad):

102

Capítulo 3: Flujo de fluidos

103

Escurrimiento de fluidos

USACH

EJERCICIOS 3.1. Determinar la potencia entregada por la bomba que descarga 28 lt/s de agua en el sistema mostrado en la figura. Las pérdidas del sistema son equivalentes a 10 · v2/2 y el diámetro interior de la cañería es de 15 cm.

3.2. Un líquido circula con flujo estacionario, a través de una tubería de 3" de diámetro interno. La densidad es 1050 kg/m3 y la viscosidad de 2 cp. Empleando un tubo de Pitot, se obtuvieron los siguientes valores de velocidades puntuales:

V(m/s ) (Pulg.) -------------------------2,28 0 2,26 0,15 2,23 0,30 2,18 0,45 2,14 0,60 2,09 0,75 2,01 0,90 1,89 1,05 1,77 1,20 1,53 1,35 1,11 1,425 0,00 1,59

Calcular: a) V empleando el gráfico V/Vmáx vs Re máx. b) α empleando gráfico α vs Re. c) V empleando la definición e integrando numéricamente. d) α empleando la definición e integrando numéricamente.

104

Capítulo 3: Flujo de fluidos

3.3. Para transportar aceite desde un depósito A a otro B con un caudal de 200 lt/min., es necesario instalar un grupo motobomba, cuya potencia se desea determinar. La tubería es de 3" de diámetro y las pérdidas a lo largo del sistema son de 3,5 kgf · m/kg. El nivel del aceite en el depósito B se mantiene 12 m. por encima del nivel del estanque A. En las condiciones de transporte, la densidad del aceite es de 840 kg/m3. Obtenga la potencia de la bomba y compare sus resultados al realizar balances entre: a) 1 – 4. b) 2 – 3. c) 2 – 4.

3.4. Un motor suministra 20 HP a una bomba por la que circula agua a razón de 500 gal/min. 25% de la potencia se consume elevando la energía interna del agua y para vencer el roce en las partes mecánicas de la bomba. El diámetro de la cañería de entrada a la bomba es de 6" y el de descarga de la bomba es de 4". Ambas tuberías están a la misma altura. Calcular el aumento de presión del agua.

105

Escurrimiento de fluidos

USACH

3.5. Circulan 3600 ton/h de agua desde un estanque, elevado a una turbina ubicada a menor nivel, a través de un conducto circular de diámetro uniforme. En un punto del conducto ubicado a 300 pies sobre la turbina la presión es 30 psia. A 10 pies bajo la turbina la presión es 18 psia. La potencia de la turbina a la salida del eje es de 1000 Hp. Si la eficiencia de la turbina es de 90%, a) Calcule las pérdidas por fricción en el conducto. b) Si no se produce transferencia de calor al ambiente, ¿En cuánto podría subir la temperatura del agua? 3.6. En el sistema de la figura circula agua. Calcular la potencia entregada por la bomba si el caudal es de 80 lt/s. El diámetro de succión es de 6" y el de descarga es de 4". Considere despreciables las pérdidas.

3.7. A través de una cañería horizontal de acero comercial, tipo 40 de D.N. 10", circulan 2 m3/min. de agua a 20 ºC. La cañería de 20 m. de longitud tiene además un codo standard de 90º, una válvula de retención (“check”), una válvula de compuerta (“gate”) y una válvula de globo (“globe”), todas totalmente abiertas ¿Cuál es la caída de presión? 3.8. En la figura se muestra un sifón, el cual es utilizado para extraer agua desde un estanque. El sifón está construido con cañería de 10" y su

106

Capítulo 3: Flujo de fluidos



longitud total es de 18 metros. Cuando el agua se encuentra en su mínimo nivel, calcule: a) El flujo. b) La presión en el punto A. La curvatura del tope del sifón es equivalente a dos codos de 90º standard.

3.9. Dos grandes recipientes de agua están conectados por una cañería de acero comercial de 8 pulgadas de diámetro nominal y 5000 pies de largo. El nivel de uno de los recipientes está 200 pies por sobre el nivel del otro, ambos abiertos a la atmósfera. ¿Cuál es el flujo de agua que circula? 3.10. Se desea transportar 100 gal/min. de agua, desde un lago a un estanque ubicado a una altura de 100 metros. La longitud de la cañería debe ser de 400 m. y contar con 10 codos de 90º y 2 válvulas de compuerta. La temperatura media es de 15 ºC. a) Determine el diámetro óptimo económico. b) ¿Cuál es la potencia de la bomba en HP? c) ¿Qué sucede con sus respuestas a) y b) si el caudal a transportar aumenta al doble?

107

Escurrimiento de fluidos

USACH

3.11. Para concentrar una solución de NaCl, ésta se bombea desde un depósito hacia un evaporador, a través de una tubería lisa de P.V.C. de diámetro interno 3 cm, a razón de 8 m3/h. A la temperatura de bombeo, la solución tiene una densidad de 1150 kg/m3 y una viscosidad de 2,3, cp. Si el trayecto a recorrer es de 350 m, ¿Cuál es la potencia de la bomba requerida? 3.12. Se retira agua de un depósito y se bombea una longitud equivalente de 2 millas, a través de un ducto horizontal, circular, de concreto de 10 pulgadas de diámetro interno. Al final de este ducto el flujo se divide en 2 cañerías de acero comercial, cédula 40, una de 4" (D.N.) y otra de 3" (D.N.). La línea de 4" tiene una longitud equivalente de 200 pies y se eleva hasta un punto a 50 pies sobre la superficie del agua del depósito, donde el flujo se descarga a la atmósfera. Este flujo debe mantenerse en 100 gal/min. La línea de 3" (horizontal) descarga también a la atmósfera en un punto a 700 pies, desde la bifurcación y al nivel del agua en el depósito. Determinar los HP entregados a la bomba, la que tiene una eficiencia de 70%. 3.13. Una bomba de 5 C.V. con una eficiencia de 70% toma amoniaco al 20%, desde un depósito y lo transporta a lo largo de una tubería de 100 m. de longitud total, hasta el lugar de descarga situado a 15 m. por encima del lugar de succión. Determine el diámetro de tubería a emplear si el caudal que circula es de 10 m3/h. 3.14. Debe bombearse una pasta de dientes (no-newtoniana), a través de una tubería de acero inoxidable de 50 mm d.i., desde la máquina de mezclado de los ingredientes hasta la máquina de llenado de los tubos de pasta. La longitud equivalente de la línea, incluyendo las pérdidas en codos, uniones, entrada y salida, es 10 m y la velocidad media del flujo es de 1 m/s.

108

Capítulo 3: Flujo de fluidos



a) ¿Qué diferencia de presión (en at.) dará este caudal? b) ¿Que tamaño de motor hará el trabajo para una eficacia de la moto bomba de un 30%?

3.15. Como se muestra en la figura, desde un recipiente fluye agua a través de una cañería a un punto de bifurcación, desde donde fluye a otros recipientes mediante cañerías separadas. Las cañerías son de acero comercial, tipo 40. Calcule el flujo en (gal/min) de agua que llega a cada recipiente, suponiendo flujo estacionario.



Sección 1 = 2000', cañería de 3" Sección 2 = 2000', cañería de 3" Sección 3 = 1000', cañería de 3"

109

Capítulo 4: Medidores de flujo

CAPÍTULO 4 MEDIDORES DE FLUJO

Para cuantificar el flujo de fluidos existen métodos directos e indirectos. En los primeros se determina el peso o volumen de fluido, que pasa a través de una sección en un determinado tiempo, registrado con un cronómetro. En los métodos indirectos, se determina la altura, la diferencia de presiones o de velocidades en varios puntos de una sección, calculándose el flujo con esta información. Adicionalmente, existen métodos basados en conceptos físicos de electrónica, electromagnetismo, óptica, ultrasonido, entre otros. La medición puede obedecer a distintos fines: información, control automatizado de un proceso, registro de variaciones de caudal en el tiempo, etc. Para estos efectos existen en el mercado una variedad de dispositivos cuyo funcionamiento obedece a distintos principios. Cada uno de ellos tiene una serie de ventajas y limitaciones, que lo hacen más o menos recomendable para una aplicación determinada, en función del fin que se persiga y del medio en el cual se utilizará (gas, líquido, fluido corrosivo o no, sólidos, etc.). A continuación se dará una breve descripción de medidores que funcionan en base a principios electrónicos, magnéticos, radioactivos, etc. Al final, se describirán instrumentos que funcionan en base a principios hidrodinámicos, indicándose además sus ecuaciones de diseño. A) Anemómetro de hilo caliente: Este medidor se utiliza para determinar velocidades puntuales o locales, en ductos por los que circulan gases. El sensor es un hilo de platino (a veces revestido en cuarzo), que es calentado por una corriente eléctrica. La resistencia eléctrica del hilo es función de su temperatura. El flujo de gas alrededor del hilo caliente lo enfría y de este modo varía su resistencia eléctrica. Dejando constante la tensión o la intensidad de

111

Escurrimiento de fluidos

USACH

la corriente en el hilo mediante un circuito conveniente, la variación en la intensidad o tensión, respectivamente, son función de la velocidad del flujo gaseoso que rodea al hilo caliente. Este hilo debe ser calibrado previamente, colocándolo en una corriente gaseosa de velocidad conocida. B) Medidores electromagnéticos: Estos se basan en la ley de Faraday, según la cual, cuando un conductor eléctrico se mueve dentro de un campo magnético, se induce en él una fuerza electromotriz que es proporcional a su velocidad. Un medidor electromagnético consiste esencialmente en un carrete de tubería de material no magnético, que lleva adosado una serie de bobinas las cuales, una vez conectadas a un circuito eléctrico, crean un campo magnético transversal al tubo. Al atravesar este campo el fluido, hace el papel del conductor eléctrico; la fuerza electromotriz inducida en el fluido, cuyo valor es proporcional a la velocidad media del flujo, es detectada por dos sensores instalados en paredes opuestas al tubo. Un circuito electrónico auxiliar recibe las señales de los sensores y, después de procesar esta información determina la velocidad o caudal que está circulando.

Figura 4.1. Medidor electromagnético.

112

Capítulo 4: Medidores de flujo

La gran ventaja de este medidor es que no presenta ninguna obstrucción al paso del fluido, e incluso su pared interna se puede recubrir con revestimientos adecuados, para evita el ataque de fluidos corrosivos. Puede ser utilizado para la medición de fluidos pastosos. La principal desventaja es su alto precio, que es en función del diámetro de la tubería. C) Medición radioactiva: La medición radioactiva se puede emplear en conducciones abiertas o cerradas. El método consiste en inyectar, en un punto del circuito, un trazador radioactivo y contar el tiempo que transcurre, hasta ser detectado por un contador Geiger, situado aguas abajo del anterior. Conocida la distancia entre ambos puntos, se determina la velocidad del fluido. Este método tiene el inconveniente, además del derivado de las precauciones que hay que observar en el manejo del material radioactivo, de que no se tiene una medición continua, sino únicamente la correspondiente al momento de introducir el trazador. D) Medición por ultrasonidos: Existen distintos dispositivos que se basan en el efecto Doppler o en la recepción de ecos. En ambos métodos se utiliza un elemento emisor, que envía impulsos de ultrasonido a través del fluido en distintas direcciones. Los impulsos que viajan en el mismo sentido que el fluido, lo hacen a mayor velocidad que aquellos otros que viajan en sentido contrario. Una serie de sensores recogen estos impulsos y un circuito electrónico auxiliar procesa la información, calculando el tiempo transcurrido entre la emisión y la recepción del impulso. Como las distancias entre el emisor y los sensores son fijos y conocidos, se puede obtener la velocidad a la que viajó el impulso que, a su vez, es función de la velocidad del fluido. E) Medición por impacto: La medición se efectúa introduciendo en el seno del fluido a medir una paleta, que se desplaza en la dirección de la corriente. La fuerza que ejerce el fluido en movimiento sobre esta paleta será proporcional al caudal que circula. Aunque este dispositivo se puede emplear con

113

Escurrimiento de fluidos

USACH

cualquier fluido, su utilización está en general limitada a la medición de líquidos extremadamente viscosos o de alto punto de fusión: pastas, asfaltos, resinas, azufre fundido, etc.

Figura 4.2. Medidor de impacto.

F) Anemómetro de aire: En este medidor, como se muestra en Figura 4.3, las aspas actúan sobre contadores que indican el número de revoluciones, las que se relacionan con la velocidad del gas.

Figura 4.3. Anemómetro de aire.

G) Medidores de desplazamiento positivo (volumétricos): Estos medidores constan de émbolos o tabiques, que son desplazados por la corriente fluida y de un mecanismo contador, que registra el número de desplazamientos en una unidad conveniente, como litros o metro cúbicos. Un medidor volumétrico es el utilizado en la mayoría de sistemas domésticos de distribución de agua. El caudal se obtiene determinando con un cronómetro el tiempo que tarda en pasar un volumen dado de fluido.

114

Capítulo 4: Medidores de flujo

H) Medidor de turbina: El principio de funcionamiento del medidor de turbina es similar al del anemómetro de aire, con la diferencia de que se hacen impermeables y se colocan en el interior de una tubería, para medir caudales de fluidos. El núcleo del rodete puede contener un magneto, que al girar el rotor, produce un campo magnético variable. El campo se detecta en una bobina eléctrica que va montada en la armadura externa de la unidad. La frecuencia de los impulsos magnéticos indican el caudal en un equipo de lectura. Debido a que la construcción y control de estos medidores es complicado, estos medidores son de elevado precio, aunque su gran exactitud los hace recomendables en operaciones automatizadas. I) Medición por producción de torbellinos: Cuando un fluido en movimiento se encuentra en su camino con un obstáculo, las líneas de flujo intentan seguir el contorno de éste. Cuando el cuerpo no tiene un contorno aerodinámico, el fluido no se puede adaptar a su forma y, entonces, se crean torbellinos aguas abajo del cuerpo en la zona de menor presión. Estos torbellinos se desprenden alternativamente de ambos lados; la frecuencia con que se desprenden de un lado u otro resulta ser proporcional a la velocidad del fluido. La velocidad en ambos lados del obstáculo nunca es la misma. En efecto, cuando se desprende un torbellino, la velocidad del fluido aumenta en ese punto, y como consecuencia, disminuye la presión. En el lado opuesto ocurre lo contrario, por lo que en todo momento existe una diferencia de presiones a ambos lados del obstáculo, cuyo valor cambia con una frecuencia idéntica a la del momento del desprendimiento de los torbellinos. Midiendo las variaciones de presión a ambos lados del obstáculo se puede determinar el caudal. J) Codos: Un codo de una tubería puede utilizarse en ciertas condiciones, para estimar el caudal (Figura 4.4). En efecto, la fuerza centrífuga que actúa sobre el fluido al girar 90º, produce que la presión en la pared externa de la curva tenga un valor superior al que existe en la pared interna. Esta diferencia se relaciona con el caudal. Este método de medición, en determinadas

115

Escurrimiento de fluidos

USACH

condiciones tiene la suficiente reproducibilidad como para permitir su uso con fines de control.

Figura 4.4. Medidor de codo.

4.1.

Medidores que funcionan en base a principios fluidodinámicos

K) Medición por diferencia de presiones: La obtención de velocidades y/o caudales, a partir de la información de la diferencia de presión entre dos puntos de la conducción es el principio en el cual se basa la mayoría de las mediciones de fluidos, en conducciones cerradas. Dentro de los instrumentos que funcionan en base a este método se tiene a: tubos de Pitot, placa orificio, Venturi y tobera (o boquilla). K.1) Tubo de Pitot. Este dispositivo se caracteriza por medir velocidades puntuales ( V ) de fluidos. Existen diversas configuraciones, pero en esencia consisten en dos tubos abiertos: uno está orientado de forma que su abertura se enfrente al movimiento del fluido (tubo de impacto). El otro, tubo estático, está orientado de forma que no le influya la presión debida a la velocidad del fluido. El tubo de impacto detecta tanto la presión estática del fluido como la equivalente a su energía cinética. La diferencia de ambas mediciones dará por resultado la presión de impacto, a partir de la cual se puede obtener la velocidad del fluido. En la Figura 4.5 se muestran algunas configuraciones de tubos de Pitot.

116

Capítulo 4: Medidores de flujo

A continuación se hará la deducción de la ecuación del tubo de Pitot simple. Como se muestra en la Figura 4.5(a), la abertura del tubo horizontal está dirigida agua arriba, de modo que el fluido impacte la abertura.

Figura 4.5. a) Tubo de Pitot simple, b) Tubo de Pitot compuesto, c) Tubo de Pitot comercial, d) Tubo de Pitot traverso.

El punto 2 se llama punto de estancamiento y en éste, la velocidad del fluido es cero. La toma piezométrica bajo el punto 1 detecta la presión estática del fluido. Aplicando la ecuación de Bernoulli a la línea de corriente entre 1 y 2 se obtiene: (4.1)

117

Escurrimiento de fluidos

USACH

Considerando que el fluido es incompresible, que Z2 = Z1 y que  V2 = 0 se obtiene: (4.2)

Además entre 1 y 2, Eˆ v ≈ 0 , luego: (4.3)



A fin de compensar las aproximaciones, se introduce un factor experimental Cp, el cual se conoce como coeficiente de descarga del tubo de Pitot. Para tubos bien diseñadosCCpp ≅ 1. 1 Utilizando ecuación de la hidrostática se obtiene: Δp = g ΔH (ρm – ρ), luego: (4.4) Todas las versiones de tubos de Pitot tienen la misma ecuación, que relaciona la velocidad puntual con la diferencia de presiones que señala algún manómetro. Dado que este medidor determina velocidad local, lo usual es situarlo en el centro de la tubería, a fin de calcular con la ecuación 4.4 la velocidad máxima. Obtenido este valor, es posible obtener la velocidad media utilizando Figura 4.6. Para obtener buenas medidas, es necesario que el Pitot esté ubicado lejos de cualquier perturbación: a lo menos uno 100 diámetros de distancia de cualquier accesorio. El medidor en sí no produce alteraciones en el flujo, es de bajo costo y ofrece posibilidades de determinar velocidades puntuales en cualquier ducto y posición. Si es inevitable la presencia de accesorios cerca del medidor, estos provocarán distorsión en el perfil de velocidades y ya no podría usarse

118

Capítulo 4: Medidores de flujo

la Figura 4.6 para obtener la velocidad media. En este caso deben hacerse lecturas en diferentes posiciones (radiales) del ducto y obtener la velocidad media, usando la definición de ésta, vale decir:



(4.5)

La selección de los puntos adecuados de medición, es decir, los ra-

dios en el caso de ductos circulares, o la posición x,y en ductos de sección rectangular, implica dividir el área de flujo en N sub-áreas iguales, colocando el medidor en el centro de cada una de estas sub-áreas. Luego, estos . valores ( V , r) se introducen en la integral anterior para así obtener V.

Figura 4.6. Gráfico para obtener velocidad media a partir de la velocidad máxima, obtenida en el centro de la tubería (Adaptado de Foust y col., 1961).

119

Escurrimiento de fluidos

USACH

K.2) Placas orificio, tobera y Venturi. Estos medidores de ΔP variable con el flujo, tienen ciertas características similares, siendo posible deducir una expresión general para todos ellos.

Figura 4.7. Placas orificio.

Una corriente forzada aumenta su velocidad, al pasar a través de una estrangulación a costa de perder presión. Este fenómeno se representa en Figura 4.8, para el flujo a través de una placa orificio.

Figura 4.8. Flujo a través de una placa orificio.

Al pasar el fluido por un orificio de menor diámetro que el conducto, la vena fluida experimenta una contracción considerable. El punto en que la sección transversal es mínima no se encuentra en el orificio mismo, sino aguas abajo de él. Esta zona se llama “vena contracta” y es donde tienen lugar las máximas diferencias de presión (caída de presión temporal). Aunque

120

Capítulo 4: Medidores de flujo

la presión vuelve a aumentar una vez pasada esta zona como consecuencias de la disminución de velocidad, nunca llegará al valor original, produciéndose una “pérdida permanente” de presión, la cual depende de las características de la placa orificio, Venturi o tobera. Consideremos el esquema siguiente, válido para los medidores recién descritos:

Figura 4.9. Esquema general para medidores de ΔP variable con el flujo.

En esta situación los puntos 1 y 2 se encuentran muy próximos, separados una distancia despreciable comparada con la longitud total del sistema. El punto 3, denominado “aguas abajo” no es afectado por la presencia del medidor. Aplicando Bernoulli entre 1 y 2, considerando fluido incompresible y que Z1 = Z2, se obtiene: (4.6) Luego: (4.7)

121

Escurrimiento de fluidos

USACH

De un balance de masa entre 1 y 2 (ecuación de continuidad) se obtiene: (4.8) Reemplazando en la ecuación de Bernoulli, se obtiene: (4.9)

Esta es una expresión general para fluidos incompresibles, para medidores de Δp variable. El término (–Δp/ρ –Êv) representa el efecto en el fluido del cambio de energía cinética. Es decir, a causa del estrechamiento, el fluido aumenta su velocidad, lo cual se traduce en una disminución de presión, asociado con una disipación de energía por fricción. Placa orificio. Este medidor consiste básicamente de una placa plana, con un orificio concéntrico, aunque existen diseños con algunas variaciones, tal como muestra la Figura 4.7. Aplicando la ecuación general 4.9 a este medidor (Figura 4.8) y considerando que el término Êv es proporcional a –Δp (caída de presión temporal, que comenzaremos a llamar –Δpt), se tendrá: (4.10) Con lo cual se obtiene: (4.11)

Definiendo un coeficiente de contracción Cc como el cuociente entre la sección de flujo de la vena contracta A2 y la sección del orificio se obtiene: (4.12)

122

Capítulo 4: Medidores de flujo



Entonces se tiene que: (4.13)

Usando la ecuación de continuidad es posible cambiar V2 por Vo (velocidad en el orificio): V2 A2 = Vo Ao, luego: (4.14)



Llamando β = D0/D1, entonces: (A0/A1)2 = β4. Luego:



(4.15)

Multiplicando arriba y abajo por 1 − β 4 , con el objeto de definir un coeficiente que caracterice al medidor, de modo que:

(4.16)

Se define Co, coeficiente de descarga del medidor de orificio como: ,

luego:

123

(4.17)

Escurrimiento de fluidos

USACH

(4.18) El coeficiente de descarga Co depende de la geometría del orificio, de la relación de diámetros β y también del número de Reynolds. En principio, Co se debería determinar experimentalmente para cada caso. Sin embargo, existen gráficos de Co vs Re con β como parámetro que permiten evaluar Co. En Figura 4.10 se muestra este gráfico para una placa orificio centrada, con bordes afilados, con toma de presión en la vena contracta. Para estas condiciones se cumple que con Re > 30.000, Co se aproxima a 0.61 y se hace independiente de Re y β.

Figura 4.10. Coeficiente de descarga para placa orificio y rotámetro en función del número de Reynolds (Adaptada de Brown G., Operaciones básicas de la Ingeniería, 1963).

124

Capítulo 4: Medidores de flujo

Para ubicar el punto de la vena contracta (punto de mínima presión) se puede utilizar la Figura 4.11, donde se presenta la distancia desde el plato orificio (como fracción del diámetro interno de la tubería) vs β. La primera toma piezométrica se ubica un diámetro interno de tubería antes de la paca orificio.

Figura 4.11. Ubicación de la vena contracta (Adaptado de Spink, L. Principles and Practice of Flow Meter Engineering, 1958).

En la Figura 4.12 se sugieren las mínimas longitudes de tubería recta, que deberían respetarse antes y después del plato orificio, dependiendo de los accesorios presentes en la tubería.

125

Escurrimiento de fluidos

USACH

Figura 4.12. Distancias recomendadas de tubería recta antes y después de la placa orificio.

Para evaluar la caída de presión permanente provocada por la presencia de medidores de placa orificio, con orificio centrado, se dispone de la expresión: (4.19) Para la medición de líquidos o gases limpios y no corrosivos, se utiliza la placa de orificio concéntrico y canto vivo. Si el gas arrastra ligeras cantidades de condensado o el líquido está próximo a las condiciones de evaporación, se pueden practicar, si la placa va a ser instalada en un tramo horizontal, orificios de drenaje o venteo debajo o encima del orificio principal. Si las cantidades de condensado son significativas, se utiliza la placa de orificio excéntrico. Si el fluido arrastra sedimentos, se puede utilizar la placa de orificio segmentado, en la cual el radio de este orificio es el 98 por 100 del radio interno de la tubería y la cuerda del círculo también tiene el canto vivo, en su cara de entrada.

126

Capítulo 4: Medidores de flujo

Es fundamental que el canto de entrada en los orificios de las placas sea afilado, a fin de asegurar la exactitud de la medida. Un desgaste incipiente de este canto puede introducir errores de hasta un 20 por 100 del fondo de escala. Por otro lado, la placa orificio de canto afilado puede introducir grandes errores para bajos caudales y elevadas viscosidades. Para estas aplicaciones se puede utiliza la placa de orificio con canto en cuarto de círculo. El radio de curvatura del borde de entrada es función del diámetro del orificio. Venturi Una de las desventajas de la placa de orificio es que produce una gran pérdida irrecuperable de presión. Otra desventaja es su poca precisión cuando se trata de medir fluidos con sólidos en suspensión y el inconveniente que supone colocar barreras al paso de estos fluidos. Tales inconvenientes se soslayan en el tubo de Venturi. El tubo de Venturi es, en esencia, un tubo con una garganta que presenta una forma tronco-cónica aguas arriba y después del estrechamiento. El ángulo de entrada suele ser de 21º ya que ángulos mayores producirían cavitación en los líquidos al atravesar la garganta, y menores supondrían un tubo demasiado largo. El ángulo de salida está entre los 5 y 15º con el fin de recuperar el máximo de presión. (Figura 4.13). Entre las desventajas del tubo Venturi está su construcción delicada, que se traduce en un elevado precio, y la necesidad de disponer de largos tramos rectos, aguas arriba y abajo, en su instalación.

127

Escurrimiento de fluidos

USACH

Figura 4.13. Dimensiones críticas para el medidor Venturi (Adaptado de Spink, L. Principles and Practice of Flow Meter Engineering, 1958).

La ecuación general de medidores de ΔP variable con el flujo, ecuación 4.9, es aplicable al Venturi, obteniéndose:

(4.20) Donde Vv es la velocidad en la garganta del Venturi, Cv es el coeficiente del Venturi. En la Figura 4.14 se presenta un gráfico de Cv vs Re en la tubería. De allí se observa que para Re > 200.000, Cv = 0.985. En esta figura la línea continua representa el promedio de la data experimental disponible, las líneas punteadas indican el rango de dispersión de la información experimental.

128

Capítulo 4: Medidores de flujo

Figura 4.14. Coeficiente de descarga para medidores Venturi. El punto 1 se refiere a la velocidad y diámetro en la tubería (Adaptado de Foust y col., 1961).

La caída de presión permanente en un Venturi es entre 10 y 15% de la caída de presión temporal, cuando el ángulo en la salida está entre 5 y 7% y de 10 y 30% para ángulos del orden de 15º, siendo mayores las pérdidas para bajos valores de β. Toberas (boquillas) Las toberas son elementos derivados del Venturi, que se insertan en los conductos para producir la estrangulación del fluido. En ellos se ha eliminado la zona de salida existente en el tubo Venturi y la zona convergente o de entrada queda convertida en una forma más redondeada. Se utilizan formas muy distintas para la boca de entrada de las boquillas, pero las más utilizada es la de cuadrante de elipse. El costo de estos elementos, su tamaño y la exigencia de tramos rectos a la entrada son muy inferiores a los del tubo Venturi. Producen una pérdida irrecuperable de presión ligeramente superior a las del tubo Ventu-

129

Escurrimiento de fluidos

USACH

ri, pero muy inferior a la producida por la placa de orificio. En Figura 4.15 se muestran algunas toberas.

Figura 4.15. a) Dimensiones para la construcción de una tobera, b) Vista en corte de una instalación que incluye una tobera, c) Esquema de una tobera.

También para las toberas es aplicable la ecuación 4.9, obteniéndose: (4.21) En Figura 4.16, se presenta un gráfico de C vs Re (basado en el diámetro de la tubería), donde: (4.22)

130

Capítulo 4: Medidores de flujo

Figura 4.16. a) Vista en corte de una tobera, d1 es el diámetro de la cañería y do es el diámetro interno de la tobera. b) Gráfico de coeficiente de descarga, C, versus Re1 (Adaptado de Spink, L. Principles and Practice of Flow Meter Engineering, 1958).

(4.23) La expresión recién escrita permite evaluar la caída de presión permanente, provocada por la presencia de la tobera. Para el flujo de gases, las ecuaciones para el plato orificio, Venturi y tobera, son modificadas por el factor Y, el cual considera que el gas se expande adiabáticamente desde P1 a P2. Para el plato orificio, la ecuación es: (4.24) Donde p1: presión antes del plato orificio. p2: presión en la vena contracta. En la Figura 4.17 se entrega el factor Y para orificio, Venturi y tobera:

131

Escurrimiento de fluidos

USACH

Figura 4.17. Factor de expansión Y para orificios, Venturi y tobera. r = P2/P1 y k = Cp/Cv.

Ejemplo 4.1 Para medir petróleo que circula a 15 ºC a través de una cañería de acero comercial tipo 40 de 2", se instalará una placa orificio. Se desea que para un caudal de 3(l/s) el manómetro marque una diferencia de nivel de 40(cm). El líquido manométrico es mercurio. Calcule: a) Diámetro del orificio (D0). b) La pérdida de potencia provocada por el orificio, para un caudal de 3(l/s). Solución: ρpetróleo = 965(kg/m3) μpetróleo = 3,5(kg/m · s) ρmercurio = 13600(kg/m3) Dc = 0,0525(m) Q = 0,003(m3/s)

132

Capítulo 4: Medidores de flujo



La ecuación para el plato orificio es:



(a)

Del enunciado ⇒ –∆p f == 0,00425 ∆H · g · (ρmercurio – ρpetróleo) = 0,4(m)·9,8(m/s2)·(13600 – 965)(kg/m3) f = 0,= 00425 51,33(m2/s2) –∆p = 49529,2(kg/m · s2) ⇒ –∆p/ρ Inicialmente consideraremos: β4 ≈ 0, luego:

Supondremos C0 ≈ 0,61 ∴D D00 == 0,0249(m) 0,0249(m ) 0,00425 ⇒ βf4 = 0,05 0,00425 Con este valor de D0 ⇒ βf ==0,473 Introduciendo este valor en la ecuación (a) y considerando C0 ≈ 0,61

0,00425 ⇒ βf4 = 0,047 0,00425 Con este valor de D0 ⇒ βf ==0,4667 Introduciendo este valor en la ecuación (a) se obtiene D0 = 0,0246(m), eRe0 0==4242,8 ,8 ⇒ 0,600425 4 Con este valor de C0 se obtiene final⇒CCf00 ≅ = 0,64. para el cualR mente D0 = 0,024(m).

133

Escurrimiento de fluidos

USACH

De la ecuación correspondiente se obtiene:

Luego La potencia perdida en el medidor es: Potencia = 40(m/s)2 · 0,003(m3/s) · 965(kg/m3) = 116(kg · m/s3) Potencia =

Ejemplo 4.2 Agua asciende a través de un medidor Venturi de coeficiente 0,98 y β = 0,5. Un manómetro en U indica una diferencia de nivel de 3,88(ft), siendo el líquido manométrico un fluido de densidad 78(lb/ft3). ¿Cuál es el flujo de agua en (ft3/s)?

134

Capítulo 4: Medidores de flujo

Solución: Debe aplicarse ecuaciones de la hidrostática al manómetro, para poder determinar Δp: pa = p1 + ρ · g ·(Y +3,88) pb = p2 + ρ · g ·(1,5+ Y) + ρM · g · 3,88 pa = pb = p1 + ρ · g · Y + ρ · g · 3,88 = p2 + ρ · g ·1,5 + ρ · g · Y + ρM · g · 3,88 p2 – p1 = ρ · g · 3,88 – ρ · g ·1,5 – ρM · g · 3,88 p2 – p1 = g · (62,4 · 2,38 – 78 · 3,88) p2 – p1 = – 4962,92(lb/ft · s2)

La ecuación del venturi es:

(b)

Pero, esta caída de presión es la suma de la caída de presión provocada por el estrangulamiento, más el efecto de la columna de líquido entre (1) y (2). Para “aislar” el efecto de caída de presión por el estrangulamiento, debe restarse el peso de la columna de líquido: – (Δp)temporal = 4962,9 – 62,4 · 32,2 · 1,5 = 1949(lb/ft · s2) Introduciendo los valores en (b):

135

Escurrimiento de fluidos

USACH

Nota: Para este mismo medidor (pero invertido), considerando que el flujo es de “bajada”, manteniendo los valores señalados por el manómetro, se obtiene: Q = 3,18(ft3/s) Ejemplo 4.3 Un tubo de Pitot se introduce en el centro de una cañería de 3" tipo 40 de acero comercial que conduce nitrógeno, registrándose una lectura de 35(mm) en un manómetro inclinado (1/10) en conexión con el Pitot. La temperatura del nitrógeno en la tubería es de 15 ºC y la presión en el lugar en que se introduce el tubo de Pitot es 850(mmHg). Determine el caudal, referido a condiciones normales de presión y temperatura. Solución: La ecuación del Pilot es:

R eConsiderando 2 ,8 ⇒ CC0p ≅ 10,64 0 = 4

p1 = 850(mmHg) = 1,118(atm) DC = 0,0779(m)

En la rama inclinada se lee 35(mm), pero la variación vertical es de 3,5(mm) (de acuerdo con la inclinación del manómetro). Por lo tanto, usando las ecuaciones de la hidrostática se obtiene: (p2 – p1) = ΔH · g · (ρM – ρN2)

136

Capítulo 4: Medidores de flujo



Considerando comportamiento de gas ideal para el nitrógeno:

(p2 – p1) = 0,0035(m) · 9,8(m/s2) · (1000 – 1,33) (kg/m3) (p2 – p1) = 34,25(kg/m · s2) Introduciendo valores en las ecuaciones del Pitot se obtiene:



Esta es la velocidad máxima, ya que el Pitot se ubicó en el centro ⋅ de la tubería. Usando la figura se obtiene una razón V V max , entrando con Remax:

Este es el caudal a 850(mmHg) y 15 ºC. En las condiciones normales de (0 ºC y 760(mmHg)) será:

137

Escurrimiento de fluidos

USACH

Ejemplo 4.4 A través de una cañería de 2" de acero comercial, tipo 40, se conduce hidrógeno a 30 ºC. Para medir el flujo volumétrico (Q) se instala una placa orificio de 2(cm) de diámetro. La lectura obtenida en un manómetro diferencial de mercurio, conectado a ambos lados del plato, es de 5(cm), y la presión de hidrógeno en las proximidades de la placa es 1,5(atm). ¿Cuál es el caudal? Solución:



La ecuación para la placa orificio, para un fluido compresible es:

(–Δp) = ΔH · g · (ρM – ρ) = 0,05(m) · 9,8(m/s2) · (13600 – 0,121)(kg/m3) (–Δp) = 6664(kg/m · s2) Para el hidrógeno: Cp / Cv = K = 1,405 f =r)/K 0,00425 r = p2 / p1 = 0,956 ⇒ (1– = 0,031 Con estos valores en la figura se lee Y = 0.988. Considerando inicialmente C0 = 0,61, se obtiene V0 = 202(m/s). Con este valor se calculará el Re0, para verificar el calor de C0.

138

Capítulo 4: Medidores de flujo

Con este valor, de la Figura 4.9 se lee C0 = 0,61, luego el valor supuesto es correcto. , en las condiciones del problema. L) Medidores de área variable con el flujo: La utilización de placas, boquillas, etc., supone la creación de una caída de presión variable, a través de una estrangulación cuya área permanece constante. Para caudales pequeños hemos visto que dichos dispositivos resultan inexactos, por la relación cuadrática entre el caudal y la diferencia de presión. Se puede realizar la medición de caudales a la inversa: manteniendo constante la diferencia de presión, a través de la estrangulación y haciendo variaFigura 4.18. Rotámetros. ble el área de la misma. La cantidad del fluido que circula en la unidad de tiempo, es en este caso proporcional a dicha superficie. El dispositivo más utilizado que sigue este principio es el rotámetro.

139

Escurrimiento de fluidos

USACH

El rotámetro consiste, en síntesis, en un tubo colocado verticalmente que tiene forma tronco-cónica ensanchado hacia la parte superior. En su interior lleva un flotador. Al fluido se le hace circular hacia arriba y al pasar entre el flotador y el tubo, se crea una diferencia constante de presión, cuyo valor se puede modificar variando la masa y forma del flotador. La conicidad del tubo determina que el espacio anular por donde debe pasar el fluido aumenta de sección, conforme el caudal aumenta.

Al circular fluido, la corriente “trata” de llevar el flotador hacia arri-

ba. Para un caudal dado, existe una posición del flotador en la cual queda en equilibrio y donde se cumple: Σ Fi = 0

(4.24)

Las fuerzas que actúan sobre el flotador son:

FE = fuerza de empuje

FG = fuerza de gravedad FD = fuerza de arrastre En el equilibrio = FG = FE + FD

140

Capítulo 4: Medidores de flujo

FG = m g = ρ flotador Vflotador g

→

F G = ρF V F g

FE = ρ VF g FD = (– ΔPF) SF = (– ΔP12) CF2 SF (– ΔPF): caída de presión a través del flotador. (– ΔP12): caída de presión entre la entrada del rotámetro y la posición del flotador. SF: área de la sección transversal máxima del flotador.



Reemplazando en el balance de fuerzas:

ρF VF g = VF ρ g + (– ΔP12) SF CF2

(4.25)

De aquí se obtiene:

(4.26)

Por otro lado, en forma similar a lo demostrado para medidores de ΔP variable con el flujo, es posible plantear que: (4.27) en que:

141

Escurrimiento de fluidos

USACH

Introduciendo (– ΔP12) del balance de fuerzas, en esta última ecuación se obtiene: (4.28) (4.29) Pero, de la ecuación de continuidad: W = VR S2 ρ = Q ρ = Flujo másico S2 es la sección anular, en la posición del flotador Luego:

w = CR S2



(4.30)

CR puede ser evaluado con Figura 4.9, usando: ReR =

(4.31)

En que: D2 = diámetro del tubo en la posición del flotador. DF = diámetro del flotador en su parte más ancha. M) Medidores en canales abiertos: Existen dos tipos de medidores: rebosaderos y canales de medida. M.1) Rebosaderos. En un canal se coloca una represa cuyo rebosadero puede adoptar distintas formas. El líquido alcanzará diferentes alturas en función del caudal: a mayor caudal, mayor altura. Esta altura se mide en un tubo tranquilizador por cualquier método establecido para la medición de niveles, como se muestra en Figura 4.19.

142

Capítulo 4: Medidores de flujo

Figura 4.19. Rebosadero.

La altura está relacionada con el caudal por ecuaciones semi-empíricas que dependen del tipo de rebosadero, como se muestra en la Tabla 4.1. Tabla 4.1. Forma rebosadero Ecuación para Q de descarga en pie3/seg. Q = 1,43 x H2,5 Q = 2,49 x H2,48 Q = 3.33 H3/2 (L - 0.2 H) Q = (pie3/seg), H = pie, L = pie Estos son los rebosaderos más sencillos, existiendo otros de forma más complicada. M.2) Canales. Existen distintos tipos de canales que se diferencian por sus dimensiones y ángulos de entrada y salida. En la Figura 4.20 se muestra uno de ellos. Como se ve, consisten en esencia en intercalar una garganta, para conseguir que las variaciones de caudal se reflejen en una variación de nivel en el tubo de medida.

143

Escurrimiento de fluidos

USACH

La relación existente en los canales entre el caudal y la altura del nivel adopta la siguiente forma:

Figura 4.20. Canal Parshall.

En la Tabla 4.2 se presenta una comparación de los distintos medidores, con sus ventajas y desventajas.

144

Capítulo 4: Medidores de flujo

Tabla 4.2. Características generales y limitaciones de medidores (Revista Española de Ingeniería Química, 1978). Características generales

A favor

En contra

- Fácil instalación. - Bajo costo. - Precio independiente del tamaño tubería. - Estático. - Aceptación universal. - No necesita calibrar.

- Necesidad de tramos rectos de tubería. - Obstruye el paso. - La exactitud depende del grado de desgaste del canto. - Influyen la viscosidad y la densidad. - Alta pérdida irrecuperable. - Necesita otro elemento auxiliar para determinar la medición.

- Baja pérdida permanente. - Elemento estático. - No necesita calibrar. - Auto-limpiable. - Bajo mantenimiento.

- Necesidad de tramos rectos de tubería. - Precio alto. - Tamaño grande. - Influye la viscosidad. - Necesita otro elemento auxiliar para determinar la medición.

Boquillas y tubos

- Baratos. - Pérdida permanente aceptable. - Estáticos. - Bajo mantenimiento. - Auto-limpiables.

- Datos experimentales son escasos. - Necesita calibración. - Necesidad de tramos. - Influenciados por condiciones operación.

Tubo Pitot

- Muy barato. - Precio independiente tamaño conducto. - Fácil instalación. - Pérdida permanente mínima.

- Medición puntual. - Sujeto a ensuciamiento. - Necesidad de tramos. - Necesita otro elemento auxiliar para determinar la medición.

Placa de orificio

Servicio: Líquidos + gases Rango: 3:1 Exactitud: Buena Señal: Cuadrática Error: 3,4% Pérdida permanente: 50-90% Tramos rectos: 10-30 diámetros

Tubo Venturi

Servicio: Líquidos + gases + vapor agua Rango: 3:1 Exactitud: Buena Señal: Cuadrática Error: 1% Pérdida permanente: 10-20% Tramos rectos: 5-10 diámetros

Servicio: Líquidos + gases + vapor agua Rango: 3:1 Exactitud: Buena Señal: Cuadrática Error: 1,5% Pérdida permanente: 30-70% Tramos rectos : 10-30 diámetro

Servicio: Líquidos + gases + vapor agua Rango: 3:1 Exactitud: Mediana Señal: Cuadrática Error: 1% Pérdida permanente: Despreciable Tramos rectos: 20-30 diámetros

145

Escurrimiento de fluidos

USACH

Codo

- Muy barato. - Puede ser bi-direccional. - Auto limpiable. - Pérdida permanente casi nula.

- No recomendable para medir. - Necesita velocidad min. - Necesita calibración. - Necesita de tramos.

Rotámetro

- Barato en general. - Autosuficiente si es indicación local. - Casi auto-limpiable - Puede leer másico. - Ideal caudales bajos. - Independiente de viscosidad. - Pérdida permanente constante. - Casi estático.

- Instalación vertical solamente. - En gas necesita contrapresión. - Tamaño limitado por peso y precio.

Volumétricos

- Ideal para líquidos viscosos. - Ídem para mezclas y oper. de venta. - Pueden ser autónomos. - No necesita tramos. - Posibilidad de lectura másica. - Gran variedad de rangos.

- Desgaste mecánico. - Sólo fluidos limpios. - Limitación de peso, tamaño y precio. - Mantenimiento. - Sobre velocidad. - Propenso a daño por fluidos abrasivos.

- Exactitud excelente. - Gran rango de caudales. - Mantenimiento bajo. - Ligero. - Fácil instalación. - Casi auto-limpiable.

- Necesita conversión de señal para determinar medición. - Necesita calibrar. - Influye viscosidad. - Sólo fluidos limpios. - Necesita de tramos. - Caro. - Propenso a sobre velocidad. - Necesita contrapresión. - Sujeto a desgaste. - Re. Min. 10.000.

- Ideal para baja velocidad gas. - Posibilidad lectura másica. - Nula pérdida permanente.

- Delicado. - Sujeto a ensuciamiento por fluido. - Necesita tramos. - Precio independiente tamaño conducto. - Necesita calibrar.

Servicio: Líquidos Rango: 3:1 Exactitud: Mala Señal: Cuadrática Error: Mayor 1% Pérdida permanente: Despreciable Tramos: 20-30 diámetros

Servicio: Líquidos + gases Rango: 10:1 Exactitud: Mediana Señal: Lineal logarítmica Error: 2% Pérdida permanente pequeña Tramos rectos: No

Servicio: Líquido + gas Rango: 10:1 Exactitud: Muy buena Señal: Lineal Error: 1/2% Pérdida permanente: Baja Tramos rectos: Despreciables

Medidor de turbina

Servicio: Líquido + gas + vapor agua Rango: 14:1 Exactitud: Excelente Señal: Lineal Error: 1/2% Pérdida permanente: Despreciable Tramos rectos: 5-10 diámetros

Hilo caliente

Servicio: Líquido + gas Rango: 10:1 Exactitud: Buena Señal: Exponencial Error: 1% Pérdida permanente: Nula Tramos rectos: 5-10 diámetros

146

Capítulo 4: Medidores de flujo

Electromagnético

- No influenciado por densidad ni viscosidad. - Pérdida permanente nula. - Bi-direccional. - Estático. - Sanitizado. - Auto-limpiable.

- Caro. - Voluminoso y pesado en tamaños grandes. - No funciona con gas ni fluidos no conductores. - Limitación por temperatura. - Necesita calibrar. - Necesita potencia eléctrica.

Medidor de impacto

- Barato. - Casi estático. - Ideal para fluidos muy viscosos. - Idem caudales bajos. - Instalación fácil.

- Necesita de tramos. - Necesita calibrar.

Medidor ultrasonidos

- Independiente tamaño del tubo. - Paso diáfano. - Sanitizado. - Bi-direccional. - Estático. - Pérdida permanente nula.

- En fase de desarrollo. - Necesita calibrar. - Necesita de tramos. - Únicamente líquidos.

- Precio razonable. - Mantenimiento bajo. - No necesita calibrar. - No influyen la viscosidad ni densidad. - Estático. - Instalación sencilla. - Independiente del tamaño del tubo.

- No adecuado para líquidos sucios ni abrasivos. - No adecuado para líquidos viscosos. - Necesita de tramos. - Limitación de velocidad en líquidos por cavitación. - Necesidad de velocidad mínima.

Servicio: Líquido Rango: 30:1 Exactitud: Muy buena Señal: Lineal Error: 1% Pérdida permanente: Nula Tramos rectos: No necesarios

Servicio: Líquido + gas + vapor agua Rango: 10:1 Exactitud: Mediana Señal: Cuadrática lineal Error: 2% Pérdida permanente: Pequeña Tramos rectos: 5-10 diámetros

Servicio: Líquido Rango: ± 40 pies/seg.

Exactitud: Buena Señal: Lineal Error: 1% Pérdida permanente: Despreciable Tramos rectos: 5-10 diámetros

Medidores de torbellinos

Servicio: Líquidos + gas + vapor agua Rango: 20:1 Exactitud: Buena Señal: Lineal Error: 1% Pérdida permanente: Baja Tramos rectos: 5-10 diámetros

Nota: El error se refiere al fondo de escala. Son valores típicos que pueden variar según el tipo de instalación y de fluido. Los porcentajes de pérdida permanente se refieren al valor de la diferencia de presiones que se haya elegido en cada caso.

147

Escurrimiento de fluidos

USACH

EJERCICIOS 4.1. A través de un medidor Venturi de coeficiente 0,98 y β = 0,5 circula agua. Un manómetro en U indica una diferencia de nivel de 1m, siendo el líquido manométrico un fluido de densidad 1400 kg/m3.



a) ¿Cuál es el flujo si el agua asciende en el Venturi? b) ¿Cuál es el flujo si el agua desciende en el Venturi, manteniéndose los valores de la figura?

4.2. Nitrógeno a 20 ºC y 15 cm de agua de sobre presión fluye a través de una tubería de acero comercial de 4". Para determinar su flujo se dispone de una placa orificio (d0 =30 mm). Si el manómetro conectado a la placa indica una diferencia de presión equivalente a 1,2 m de agua, determine el flujo si la presión atmosférica es de 710 mm Hg. 4.3. Para medir un flujo de petróleo a 15 ºC (ρ = 950 kg/m3; μ = 35 cp) se desea utilizar una placa orificio, que indique una diferencia de nivel de 40 cm de Hg. Determine el diámetro de la placa orificio y la pérdida de potencia producida por el medidor de orificio.

148

Capítulo 4: Medidores de flujo

4.4. Un tubo de Pitot instalado en el centro de una tubería de 3" por la que circula oxígeno (15 ºC y 820 mm de Hg), indica una lectura de 50 mm de agua en el manómetro inclinado (1/10) conectado al Pitot. Determine la velocidad máxima, la velocidad media y el flujo másico. 4.5. En una tubería de 2" se dispone de una placa orificio de diámetro 3 cm, acoplado a un manómetro en U con mercurio como fluido manométrico. Determine el diámetro del orificio si se desea que la lectura disminuya a la mitad, cuando circula el mismo flujo de agua. 4.6. Para medir petróleo que circula a 15 ºC a través de una cañería de acero comercial tipo 40 de 2", se instalará una placa orificio. Se desea que para un caudal de 5 (l/s) el manómetro marque una diferencia de nivel de 50 cm. El líquido manométrico es mercurio. Calcule: c) Diámetro del orificio (D0). d) La pérdida de potencia provocada por el orificio, para un caudal de 5 l/s. Datos: ρpetróleo = 965(kg/m3); μpetróleo = 3,5(kg/m. s); ρmercurio = 13600(kg/m3)

149

Capítulo 5: Transporte de líquidos

CAPÍTULO 5 TRANSPORTE DE LÍQUIDOS

Los equipos impulsores de líquido reciben la denominación general de bombas. Definición de bomba La bomba es el aparato destinado a extraer, elevar, e impulsar un fluido en una dirección determinada. Convierte la energía suministrada por un elemento motriz en energía de fluido, ya sea energía cinética, potencial o de presión. Las bombas se clasifican en dos grandes grupos: de desplazamiento positivo y dinámicas. Las bombas de desplazamiento positivo se subdividen en las categorías de alternativas y rotatorias, mientras que las dinámicas comprenden las centrífugas y las bombas de efectos especiales. A.1. Bombas alternativas o recíprocas. Pertenecen a este grupo las bombas de pistón o émbolo y las bombas de diafragma. (Figura 5.1).

Figura 5.1. Clasificación de bombas de desplazamiento positivo.

151

Escurrimiento de fluidos

USACH

Figura 5.2. Clasificación de bombas dinámicas.

Figura 5.3. Esquema de una bomba de pistón de doble efecto.

152

Capítulo 5: Transporte de líquidos

Figura 5.4. Bomba de pistón de simple efecto.

Bombas de un pistón o émbolo: En éstas, un pistón o émbolo desplaza en cada carrera un volumen dado de líquido, generando un flujo pulsante. Las pulsaciones pueden ser disminuidas, usando una bomba de doble acción o aumentando el número de cilindros, como se muestra en la Figura 5.5.

Figura 5.5. Curvas de descarga para bombas de pistón. a) Acción sencilla de un pistón, b) Doble acción de un pistón, c) Acción doble de dos pistones.

La mayoría de las bombas de pistón son de doble efecto, es decir, el líquido puede admitirse a cada lado del pistón, de manera que mientras se llena una parte del cilindro, la otra se vacía. El movimiento del pistón o émbolo puede realizarse por la acción de un motor eléctrico, o mediante un ci-

153

Escurrimiento de fluidos

USACH

lindro de vapor que acciona directamente el vástago del émbolo. La presión máxima de descarga en las bombas comerciales de pistón es de aproximadamente 50 at. Para presiones superiores se emplean bombas de émbolo, en las cuales el cilindro es de pequeño diámetro y de pared gruesa. Las bombas de émbolo son de efecto sencillo, generalmente accionadas por un motor, y pueden alcanzar presiones de 1400 at, e incluso mayores. Tanto en las bombas de pistón como de émbolo, debe instalarse una válvula de seguridad y una línea de recirculación (by-pass) para evitar daños en la bomba, en el caso de un cierre inadvertido de la válvula en la descarga (Figura 5.6). Se utilizan para impulsar fluidos viscosos. Los líquidos que contienen sólidos abrasivos no pueden ser bombeados con bombas de pistón o émbolo, ya que dañan la superficie de las piezas. Tienen la desventaja de su gran tamaño y alto costo inicial.

Figura 5.6. Válvula de alivio (LOPU/USACH).

Bomba de diafragma: Éstas utilizan una membrana flexible como elemento de desplazamiento del líquido. Pueden ser móviles mecánicamente por una excéntrica o hidráulicamente por un líquido de bombeo secundario. Se usan cuando no se pueden permitir fugas del fluido o bien éste es corrosivo. Si son movidas mecánicamente, su presión de bombeo máxima es de 9-11 atm. Si la bomba es accionada hidráulicamente se pueden obtener presiones de hasta 340 atm. Los diafragmas se construyen de elastómetros, plásticos y metales (Figura 5.7).

154

Capítulo 5: Transporte de líquidos

Figura 5.7. Bomba de diafragmas.

Parámetros de diseño de una bomba alternativa Capacidad de la bomba. Es el volumen total de fluido realmente descargado por unidad de tiempo. Incluye tanto el líquido como los gases en él disueltos en las condiciones de operación (temperatura y presión de descarga). Está dada por: Q = D (1 – S), donde: Q = capacidad (L3/θ). D = desplazamiento de la bomba (L). S = fuga de líquido, en porcentaje. Desplazamiento de la bomba (D). Es el volumen desplazado por cilindro en una carrera del pistón, émbolo o diafragma sin pérdidas por fugas o por compresibilidad del fluido. D depende del desplazamiento de cada cilindro, del número de pistones, émbolos o diafragmas y del número de ciclos de bombeo.

155

Escurrimiento de fluidos

USACH

Rendimiento o eficiencia global. Es el cuociente entre la potencia de salida de la bomba y la potencia de entrada total (potencia requerida para manejar la bomba: puede tratarse de potencia eléctrica o la entregada por el vapor en las bombas de accionamiento directo). Rendimiento o eficiencia volumétrica. Es el cuociente entre la capacidad y el desplazamiento de la bomba. Condiciones de succión. El sistema de bombeo debe satisfacer ciertas condiciones en la succión, requeridas por el equipo impulsor, y dadas por el parámetro NPSH (net positive suction head). Este parámetro, similar para bombas alternativas, rotatorias y centrífugas, será definido más adelante. A.2. Bombas rotatorias. La acción de bombeo se produce por el movimiento relativo entre los elementos rotantes y los estacionarios de la bomba. Existe una gran variedad de bombas rotatorias, algunas de las cuales se muestran en la Figura 5.8. Constan de las siguientes partes: Cámara de bombeo: Es el espacio interno que se llena de líquido mientras funciona la bomba. Cuerpo o carcasa: Es la parte de la bomba que rodea la cámara y cuando permanece estacionario se llama estator. Rotor: Es la parte que rota cuando la bomba funciona y recibe un nombre específico de acuerdo al tipo de bomba (engranaje, lóbulo, tornillo, etc.). Sellos: Los hay estáticos y móviles. Los primeros dan un cierre hermético de líquido entre las partes estacionarias desmontables de la bomba.

156

Capítulo 5: Transporte de líquidos

Figura 5.8. Bombas rotatorias. a) Engranajes internos, b) Trilobular, c) Engranajes externos.

En las bombas rotatorias, la presión de descarga puede aumentar peligrosamente si se obstruye la descarga, por lo que se instalan válvulas de seguridad (alivio), que se abren a una presión determinada. Funcionamiento. En todos los tipos de bombas rotatorias se distinguen tres etapas: 1) Se abre la aspiración, aumentando el volumen de líquido a medida que rota la bomba. 2) El volumen permanece constante, se cierra la aspiración y la descarga. 3) Se abre la descarga, disminuyendo el volumen de líquido. Características de operación de las bombas rotatorias. Las bombas rotatorias son capaces de entregar una capacidad aproximadamente constante, contra cualquier presión dentro de los límites del diseño de la bomba. El flujo de descarga proveniente de una bomba rotatoria, varía directamente con la velocidad. La descarga está casi libre de pulsaciones, particularmente para las bombas de engranajes. Las características de capacidad típicas para una bomba de engranajes externos se ilustran en la Figura 5.9.

157

Escurrimiento de fluidos

USACH

Figura 5.9. Características de capacidad de bomba de engranajes rotatorios.

Bombas de engranajes. Las bombas de engranajes son el tipo más simple de bombas rotatorias. Como se muestra en la Figura 5.10, las hay de engranajes internos y externos. Un caso especial de bomba se muestra en la Figura 5.10.c, la que puede ser utilizada para bombear suspensiones, ya que sus engranajes son de polímeros.

Figura 5.10. a) Bomba de engranajes externos, b) Bomba de engranajes internos, c) Bomba de neopreno (engranajes flexibles).

En la bomba de engranajes externos, el centro de rotación de cada engranaje es externo al diámetro mayor del engranaje vecino y todos los

158

Capítulo 5: Transporte de líquidos

engranajes tienen dientes externos. Al girar los engranajes, una cantidad de líquido queda atrapado entre los dientes y es devuelto al lado de aspiración. El caudal desplazado por revolución (capacidad) depende del diámetro y ancho de los engranajes. En la bomba de engranajes internos, el líquido es introducido en el cuerpo de la bomba y queda atrapado entre los dientes del rotor y la corona dentada. La forma creciente del cabezal de la bomba, divide el líquido y sirve como un sello entre las compuertas de entrada y descarga. En la Figura 5.11 se muestra el ciclo de funcionamiento de una bomba de engranajes internos. En éstos, cualquiera de las dos ruedas de engranajes puede accionarse por el motor, siendo la otra “loca”.

Figura 5.11. Ciclo de funcionamiento de una bomba de engranajes internos.

Bombas lobulares. Son similares a la bomba de engranajes, excepto en que los engranajes son reemplazados por rotores que tienen dos o más lóbulos. Ambos rotores están accionados externamente. En la Figura 5.12 se muestran alguna vistas de bombas lobulares.

159

Escurrimiento de fluidos

USACH

Figura 5.12. Esquemas de bombas lobulares. a) Bomba bilobular, b) Bomba trilobular.

Bomba de tornillos. Existen varios tipos de bombas de tornillos, mostrándose en la Figura 5.13 una de ellas. El rotor central va accionado mecánicamente y engrana con los rotores locos, formando cavidades llenas de fluido. El fluido penetra por cada uno de los extremos de los rotores, queda aprisionado entre éstos y la caja, y es entonces impulsado suavemente hacia la cámara central de descarga. Otras bombas pueden contar con dos o más tornillos, manteniéndose el principio de funcionamiento.

Figura 5.13. Sección de una bomba de dos tornillos.

Las bombas de tornillos son eficientes, silenciosas, el caudal es continuo con poca fluctuación. Entregan presiones en la descarga de hasta 170 at. y caudales de 450 lt/min.

160

Capítulo 5: Transporte de líquidos

Bombas de aletas. Un disco provisto de ranuras, en las cuales se insertan unas aletas deslizantes, constituye el mecanismo esencial de este tipo de bomba. Las aletas mantienen contacto y se deslizan por la superficie interna de una caja elíptica, tal como se representa en la Figura 5.14. La fuerza centrífuga lanza las aletas hacia fuera y, por el giro del rotor, el espacio posterior a cada una de ellas, primero aumenta y aspira fluido, más luego disminuye y lo impele hacia afuera. Prácticamente, son las aletas las que soportan todo el desgaste y pueden reemplazarse fácilmente. Una variante de este tipo lo constituye la forma especial de aleta (tipo martillo) y la disposición mostrada por la Figura 5.14.b. Los “martillos”, moviéndose hacia afuera por la fuerza centrífuga, atrapan al fluido. Como estas aletas se desgastan, la variación de su tamaño está compensada automáticamente, hasta que llega un momento en que el cierre deja de ser perfecto, obligando a substituirlas por otras nuevas, lo que se efectúa sin necesidad de herramientas especiales.

Figura 5.14. Bombas de aletas deslizantes.

Bomba rotatoria de pistón. En lugar de engranajes, este tipo de bomba rotatoria consta de un rotor circular, montado excéntricamente en el centro del cuerpo de la bomba. La Figura 5.15 es un corte de una bomba rotatoria de pistón. El movimiento en el ciclo de operación origina un espacio para el fluido en la cámara de bombeo, mientras simultáneamente se descarga fluido a través de la válvula de salida. Esta bomba da excelentes resultados para bombear gases.

161

Escurrimiento de fluidos

USACH

Figura 5.15. Corte parcial de una bomba rotatoria de pistón.

Bomba peristáltica. Este tipo de bomba es muy utilizada en instalaciones de laboratorio, ya que permite impulsar líquidos corrosivos y/o conteniendo partículas, al no estar en contacto con la superficie de la bomba en ningún momento, debido a que lo hace por el interior de una manguera flexible. Se utiliza industrialmente para dosificar reactivos, por lo que se conoce también como bomba dosificadora.

Figura 5.16. Bomba peristáltica.

B.1. Bombas centrífugas. En su forma más simple, la bomba centrífuga consiste en un impulsor (rodete) que gira dentro de una carcasa. El fluido entra a la bomba, cerca del centro del impulsor rotatorio y es llevado hacia arriba por acción centrífuga. La energía cinética del fluido aumenta, desde el

162

Capítulo 5: Transporte de líquidos

centro del impulsor hasta los extremos de las aletas impulsoras. Esta carga de velocidad se convierte en carga de presión cuando el fluido sale de la bomba. Las bombas centrífugas se usan ampliamente en los procesos industriales, debido a la simplicidad de su diseño, bajo costo inicial, bajo mantenimiento y flexibilidad de aplicación. Se han construido bombas centrífugas para bombear cantidades tan pequeñas como unos cuantos litros/minuto, con una pequeña altura de carga y también para bombear cantidades como 2000m3/ minuto, contra alturas de carga de 100 metros. Sin embargo, el uso frecuente se restringe a volúmenes del orden de 35 m3/minuto, con presiones de descarga moderadas (inferiores a 20 at.). Existe una gran variedad de bombas centrífugas que difieren en la forma de la carcasa, tipo de rodete, succión sencilla o doble, etc. (Figura 5.17).

Figura 5.17. Bombas centrífugas.

Dependiendo de la forma de la carcasa se tienen las bombas de voluta y las bombas con difusores, como se muestra en la Figura 5.18. Si el área disponible para el flujo aumenta gradualmente entre el borde externo del impulsor y la carcasa, la bomba se llama de voluta. Este aumento en el área de flujo origina que la velocidad del fluido disminuya gradualmente, reduciendo por tanto la formación de remolinos. La carcasa con difusores tiene guías estacionarias que ofrecen al líquido una trayectoria estrecha del impulsor a la carcasa. Los difusores sirven para el mismo propósito que la voluta y ambos tipos de bomba tienen aproximadamente la misma eficiencia.

163

Escurrimiento de fluidos

USACH

Figura 5.18. a) bomba de voluta, b) bomba con difusores.

Impulsores El impulsor o rodete consiste en una serie de aletas o alabes, curvados o rectos, para que el flujo sea lo más liso posible dentro de la bomba. Al aumentar el número de alabes, se controla mejor la dirección de movimiento del fluido y se disminuyen las pérdidas por turbulencia entre aletas. De acuerdo a su construcción mecánica, los impulsores se subdividen en: Abiertos: Las aletas están fijadas a un eje; su principal desventaja es la debilidad estructural. Resultan aptos para manipular líquidos con sólidos pulposos y sólidos abrasivos. Semi-cerrados: Las aletas están apoyadas en un plato y pueden tener o no aletas suplementarias, colocadas en la cara posterior del plato. Se utilizan para transportar aguas servidas (alcantarillado). Cerrados: Son casi universalmente usados para líquidos limpios o suspensiones de pequeñas partículas. En ellos, los platos cierran completamente la trayectoria del fluido, desde la aspiración a la periferia del impulsor.

164

Capítulo 5: Transporte de líquidos

Por otro lado, de acuerdo a la dirección principal del flujo, con referencia al eje de rotación, los impulsores se clasifican en: Impulsores con flujo radial: Las superficies de las aletas son generadas por líneas rectas, perpendiculares al eje de rotación. Impulsores con flujo axial: Si la aleta forma un ángulo con el eje, el flujo es estrictamente paralelo al eje de rotación. En otras palabras, se mueve el fluido solo axialmente. Impulsores con flujo mixto: Los bordes de las aletas de un impulsor de flujo axial se curvan para dirigir el flujo perpendicular al eje; así se agrega la fuerza centrífuga a la fuerza de empuje, dada al fluido por las aletas. En la Figura 5.19 se muestran diferentes tipos de impulsores.

165

Escurrimiento de fluidos

USACH

Figura 5.19. Diseños de impulsores o rodetes (Adaptado de Foust y col., 1961).

Las bombas centrífugas pueden tener: Impulsores con aspiración única, donde el fluido ingresa por una sola entrada. Figura 5.20.a. Impulsores con doble aspiración, con el fluido entrando al impulsor simétricamente a cada lado. Figura 5.20.b.

166

Capítulo 5: Transporte de líquidos

Figura 5.20. a) Bomba centrífuga con aspiración unilateral, b) Bomba centrífuga de doble aspiración.

Bombas centrífugas de auto-cebado. Una bomba centrífuga común es incapaz de auto-cebarse e iniciar su trabajo, si en su interior contiene sólo aire o vapor. En las bombas auto-cebantes se ha modificado la ubicación de la succión, como se muestra en la Figura 5.21, con las conexiones de succión y de descarga en la parte superior de una carcasa, con dos cámaras encima del rotor, lo que permite mantener en la bomba una cantidad de líquido, aún cuando la bomba no esté en funcionamiento. Cuando se reinicia el funcionamiento de la bomba, el rodete lanza al líquido hacia los lados de la caja y comienza a bombear aire, tomándolo de la tubería de aspiración. En las cercanías de los extremos de las aletas se mezclan el aire y el líquido, formando una espuma. Esta espuma sigue alrededor del rodete hasta el borde, desviándose la espuma fuera de la zona del rodete, hacia la cámara superior. Allí el aire se eleva y escapa, agotándose de esta forma el aire en la tubería de aspiración, efectuado lo cual, la acción se normaliza y sólo se aspira líquido.

167

Escurrimiento de fluidos

USACH

Figura 5.21. Bomba auto-cebante.

Bombas centrífugas de achique, o bombas sumergidas. La bomba centrífuga se puede construir en tamaños reducidos, con dimensiones tales que puede proyectarse una unidad de múltiple efecto en una carcasa de sólo 4 pulgadas (10 cm) de diámetro. La bomba puede trabajar sumergida en el agua, al igual que el motor eléctrico (Figura 5.22).

Figura 5.22. Bombas sumergibles.

Bombas de turbina. Antiguamente, la bomba con difusores recibía el nombre de turbina por su similitud con las turbinas generadoras de energía. Actualmente, el nombre de turbina se reserva a la bomba cuyo rotor es como el mostrado en Figura 5.23. La turbina está construida de una pieza que gira y actúa como una centrífuga de rodete semiabierto.

168

Capítulo 5: Transporte de líquidos

En esta bomba la recirculación está muy favorecida. El fluido que abandona una aleta es arrastrado alrededor del canal por la propias aletas, vuelve a penetrar por la zona de admisión de las aletas y recibe así uno o más impulsos, antes de que recorra una vuelta completa por la periferia, desde el lugar de succión hasta el de descarga o impulsión. Aunque estas bombas poseen una holgura pequeña, no pertenecen al tipo de desplazamiento positivo. No es posible en ellas el “taponamiento” por aire y no necesitan ser cebadas.

Figura 5.23. a) Vista en corte de bomba de turbina, b) Funcionamiento de una bomba de turbina.

Bombas de múltiple efecto (multi-etapas). La carga proporcionada por un rotor simple viene limitada por las propias restricciones prácticas de su diámetro y velocidad de giro del rodete. Para conseguir elevadas presiones o cargas, se recurre a dos o más rodetes montados sobre un eje común que actúan en serie (Figura 5.24).

169

Escurrimiento de fluidos

USACH

Figura 5.24. Bomba múltiple etapas.

Teoría de bombas centrífugas Las ecuaciones fundamentales que relacionan la potencia, carga desarrollada y capacidad de una bomba centrífuga, se deducen para una bomba ideal a partir de los principios fundamentales de la mecánica de fluidos. Dado que el funcionamiento de una bomba real difiere considerablemente de una ideal, las bombas reales se diseñan corrigiendo la situación ideal con parámetros experimentales.

Figura 5.25. Esquemas de bombas centrífugas.

170

Capítulo 5: Transporte de líquidos

En la Figura 5.25, se representa esquemáticamente como fluye el líquido a través de una bomba centrífuga. El líquido ingresa axialmente por la toma de succión, punto a. En el orificio giratorio del rodete, el líquido se proyecta radialmente y penetra en los canales que existen entre los alabes, en el punto 1. Circula a través del rodete, y abandona la periferia del mismo en el punto 2, es recogido en la voluta y finalmente sale por el conducto de descarga de la bomba, en el punto b.

El funcionamiento de la bomba, se analiza considerando separadamente las tres partes del recorrido total: primero, el flujo desde el punto a al 1; segundo, el flujo a través del rodete desde el punto 1 al punto 2; y tercero, el flujo a través de la voluta desde el punto 2 al b. La parte esencial de una bomba es el rodete, y por esta razón se realizará en primer lugar el análisis teórico de esta etapa. En la Figura 5.26 se representa uno de los alabes de un rodete. Los vectores representan las distintas velocidades en los puntos 1 y 2, a la entrada y a la salida del alabe, respectivamente. Consideremos en primer lugar los vectores en el punto 2. Como consecuencia del diseño de la bomba, la tangente al rodete en este punto, forma un ángulo β2 con la tangente a la circunferencia descrita por el extremo del rodete.

171

Escurrimiento de fluidos

USACH

Figura 5.26. a) Velocidades a la entrada y en la descarga de un alabe de una bomba centrífuga, b) Diagrama vectorial en un extremo de un alabe.

El vector v2* es la velocidad del fluido, que vería un observador situado sobre el extremo del alabe, punto 2, y moviéndose con él, siendo por lo tanto una velocidad relativa. Supongamos ahora: primero, que todo el líquido situado sobre la circunferencia descrita por el rodete, se mueve con la misma velocidad, es decir, que el valor numérico (aunque no la dirección del vector) es v2* en todos los puntos; segundo, que el ángulo formado por el vector v2* con la tangente, es el ángulo real del alabe β2. Esta suposición, es equivalente a su vez a admitir que hay un número infinito de alabes, sin espesor y separados por una distancia infinitesimal. Este estado ideal se conoce como orientación perfecta. El punto 2 está en movimiento de rotación alrededor del eje de la bomba y se mueve con una velocidad periférica u2 con respecto al eje. El vector V2 es la velocidad resultante de la corriente de fluido, que abandona el rodete, tal como la vería un observador en reposo y recibe el nombre de velocidad absoluta del fluido. Esta velocidad, por la ley del paralelogramo, es igual al vector suma de la velocidad relativa v2* y la velocidad periférica u2. El ángulo formado por los vectores V2 y u2 se representa por α2. A la entrada de los alabes, en el punto 1, se tiene un conjunto semejante de vectores, tal como se representa en la Figura 5.26.a. En

172

Capítulo 5: Transporte de líquidos

los diseños corrientes, α1 es aproximadamente igual a 90º y el vector V1 es prácticamente radial. La Figura 5.26.b, que es el diagrama vectorial para el punto 2, indica la relaciones más útiles que existen entre los diversos vectores. Indica igualmente como se descompone el vector velocidad absoluta V2 en dos componentes: una radial que se representa por Vr2 y otra periférica Vu2. La potencia que se suministra al rodete, y por consiguiente, la potencia que necesita la bomba, se obtiene utilizando principios de la mecánica de cuerpos en rotación: Uno de ellos establece que la potencia necesaria para mantener un cuerpo en movimiento circular es igual al producto de la velocidad angular por el par, es decir: P = ω (ΔT)

(5.1)

T = par ω = velocidad angular (=) rad/tiempo

El otro principio establece que, en el movimiento de rotación de una partícula alrededor de un eje, el par es igual a la velocidad de aumento del momento de la cantidad de movimiento, estando el momento de la cantidad de movimiento definido como el producto de la cantidad de movimiento por la distancia normal, medida desde el eje a la línea de acción del vector cantidad de movimiento. Con el fin de establecer este concepto, consideremos una pequeña masa m, moviéndose en la dirección del vector V2 según muestra la Figura 5.26.b. La cantidad de movimiento en la dirección V2, es igual a m ⋅ V2. La descomposición del vector cantidad de movimiento en dos componentes proporciona el componente radial m ⋅ Vr2 y el componente tangencial m ⋅ Vu2. La línea de acción del componente radial pasa por el centro de rotación y el momento de este componente de la cantidad de movimiento,

173

Escurrimiento de fluidos

USACH

es igual a cero. La distancia desde el componente tangencial Vu2 es igual al radio r2, y el momento de la cantidad de movimiento de este componente con respecto al eje es igual a m ⋅ r2 ⋅ Vu2. Luego:

(5.2)

Al aplicar esta ecuación al flujo estacionario, a través de un volumen de control, se obtiene el balance de momento de la cantidad de movimiento, similar a un balance de cantidad de movimiento. Considerando como volumen de control el volumen total, limitado por los alabes del rodete, el par es igual a la diferencia que existe entre la velocidad de flujo de cantidad de movimiento que abandona el rodete en el punto 2, y la que entra al rodete por el punto 1. Si la velocidad de flujo de masa del líquido a través del rodete • es igual a m , entonces: •

•

Δ T = m (r2 Vu2 – r1 Vu1) = m (r2 V2 cos α2 – r1 V1 cos α1)

(5.3)

En flujo radial, α1 = 90º y por lo tanto Vu1 = 0, luego: •

T = m r2 Vu2



(5.4)

De las ecuaciones 5.1 y 5.4 se obtiene la expresión que permite determinar la potencia de una bomba ideal, sin fricción: •

Pideal = Pi = m ω r2 Vu2



(5.5)

Por otro lado, el trabajo realizado por unidad de masa de fluido es: ∧

•

W = Pi / m = ω r2 Vu2



174

(5.6)

Capítulo 5: Transporte de líquidos

Aplicando la ecuación de Bernoulli a los puntos 1 y 2, suponiendo que no existe fricción, despreciando Za – Zb, y suponiendo que la orientación es perfecta, se obtiene: ω r2 Vu2 = p2/ρ – p1/ρ + V22/2 – V12/2



(5.7)



Por otra parte, la ecuación de Bernoulli aplicada entre los puntos a -1, y b -2, proporcionan: pa/ρ + Va2/2 = p1/ρ + V12/2



(5.8)

p2/ρ + V22/2 = pb/ρ + Vb2/2



(5.9)



Sumando ecuaciones (5.7), (5.8) y (5.9) se obtiene:

pa/ρ + Va2/2 = pb/ρ + Vb2/2 – ω r2 Vu2



(5.10)

Por otro lado, definiremos altura H (dimensiones de longitud) a la expresión que incluye los términos de presión, energía potencial y cinética, divididos por la aceleración de gravedad. (5.11)

Luego, la ecuación 5.10 puede ser reescrita usando esta definición:

g Ha = g Hb – ω r2 Vu2

(5.12) ∧

g (Hb – Ha) = ω r2 Vu2 = g ΔH = W i



∧

(5.13)

Donde W i es el trabajo ideal, sin fricción, realizado por la bomba.

175

Escurrimiento de fluidos

USACH

Relación entre la carga y la velocidad de flujo de una bomba ideal. El componente radial Vr2, es la velocidad absoluta del líquido que abandona el rodete, medida perpendicularmente a la periferia del mismo. El producto de Vr2 y el área total de las secciones transversales de los canales que determinan los alabes en la periferia del rodete, Ap, da la velocidad volumétrica del flujo a través de la bomba, en metro cúbicos por segundo. Q = Vr2 ⋅ Ap

(5.14)

De la Figura 5.26 b: (5.15)

Y utilizando la definición de velocidad angular: ω =

u2 r2

Sustituyendo los valores de ω y Vr2, se obtiene:



(5.16)

De aquí se observa que la carga virtual (ideal) o variación de altura es una función lineal de la velocidad volumétrica de flujo. La pendiente de esta recta depende del signo tg β2. Si β2 es menor de 90º, la pendiente es negativa, si es igual a 90º la recta es horizontal y si es mayor de 90º, la pendiente es positiva. El primer caso corresponde a que los alabes del rodete son convexos en el sentido de giro, el segundo a alabes rectos, y el tercero a alabes cóncavos. El primer tipo de diseño es el que se elige casi siempre, ya que si la curva de carga frente a velocidad volumétrica de flujo es horizontal o de pendiente positiva, puede dar lugar a inestabilidad en el sistema de flujo.

176

Capítulo 5: Transporte de líquidos

Funcionamiento real de una bomba centrífuga. La carga desarrollada por una bomba real, es considerablemente menor que la calculada por la ecuación 5.16 para una bomba ideal. Por consiguiente, el rendimiento es menor que la unidad y la potencia utilizada por el fluido será menor que la potencia al freno. Curvas características I) Curvas ΔH vs Q: La curva real ΔH vs Q, difiere de la curva teórica, por las siguientes pérdidas: a) Flujo circulatorio: Una suposición en la teoría de la bomba ideal ha sido la de orientación perfecta, de forma que el ángulo real formado por los vectores V2 y u2, es igual al ángulo del alabe β2. En realidad, la orientación no es perfecta en una bomba real y la corriente de fluido sale con un ángulo considerablemente menor que β2. Esto produce un flujo circulatorio de líquido dentro del rodete que se superpone al flujo neto. b) Fricción del fluido: De acuerdo con la ecuación de Bernoulli, la fricción del fluido disminuye la carga. El flujo a través de los pasos y canales de la bomba, va acompañado de una fricción, que es aproximadamente proporcional al cuadrado de la velocidad volumétrica de flujo. La pérdida de carga por fricción es por lo tanto mayor para flujos elevados, tal como se indica en la Figura 5.27. c) Pérdida de choque: En una bomba de voluta, el líquido que sale del rodete en la dirección del vector V2 se introduce repentinamente en la corriente de líquido que se mueve alrededor de la carcasa. El cambio brusco de dirección provoca turbulencia y conduce a pérdidas, tanto de carga como de potencia. Las pérdidas de este tipo se denominan pérdidas

177

Escurrimiento de fluidos

USACH

de choque. Estas pérdidas son menores en una bomba con difusores, ya que éstos varían gradualmente el cambio de dirección.

Figura 5.27. Curvas de ∆H vs Q.

II) Curva de potencia vs Q. Existen varias causas que producen pérdidas de potencia. a) Fugas de líquido: La barrera al flujo de retroceso que opone el anillo de cierre no es perfecta, originándose una inevitable fuga interior, desde la descarga del rodete hacia el orificio de succión. El efecto de una fuga es el de reducir el volumen de descarga de la bomba, por unidad de potencia que se consume. b) Fricción de disco: Es la fricción que tiene lugar entre la superficie exterior del rodete y el líquido situado entre el rodete y el interior de la carcasa. Corresponde a la mayor fricción que se origina por la acción de bombeo del líquido, que tiene lugar en este mismo espacio. El líquido que está en contacto con el rodete giratorio es recogido y lanzado hacia la voluta. El líquido fluye hacia atrás a lo largo de la pared interior de

178

Capítulo 5: Transporte de líquidos

la carcasa, es tomado de nuevo por el rodete y vuelto a bombear. Esta acción secundaria inútil consume potencia. c) Pérdidas en los cojinetes: Corresponden a la potencia que se requiere para vencer la fricción mecánica que tiene lugar en los cojinetes y en las cajas prensa estopas de la bomba. d) Pérdidas por choque y fricción del fluido, ya fueron descritas para el caso de la curva ΔH vs Q. En la Figura 5.28 se representa cómo disminuye la potencia al freno consumida por una bomba, debido a las diferentes pérdidas de potencia, para quedar finalmente la potencia comunicada al fluido. El punto b corresponde a las condiciones de diseño de la bomba.

Figura 5.28. Curvas de potencia vs Q.

III) Curva de rendimiento vs Q. El rendimiento es la relación que existe entre la potencia, absorbida por el fluido (potencia hidráulica), y la potencia entregada por el motor a la bomba (potencia al freno). El rendimiento es máximo para las condiciones de diseño, y disminuye cuando las velocidades de flujo son diferentes de aquella para la que ha

179

Escurrimiento de fluidos

USACH

sido proyectada, que recibe el nombre de velocidad de flujo característica. A medida que la velocidad de flujo se aparta de las condiciones de diseño, aumentan las pérdidas por choque, disminuyendo la eficiencia de la bomba. IV) Curva de NPSH vs Q. Para asegurar el buen funcionamiento de la bomba, se deben verificar ciertas condiciones mínimas en la succión de la bomba. El parámetro NPSH es una característica de la bomba y depende de una serie de factores intrínsecos a su diseño y fabricación. El valor de NPSHr aumenta con el acrecentamiento del caudal. Cavitación. La potencia hidráulica depende de la diferencia de presiones existente entre la descarga y la succión, por lo tanto, es independiente de la presión absoluta y por consideraciones estrictamente energéticas, las presiones de succión y descarga podrían estar por debajo o sobre la atmosférica. Sin embargo, en la práctica, el límite inferior de la presión de succión está fijado por la presión de vapor del fluido a la temperatura de la succión, ya que si la presión sobre el líquido alcanza su presión de vapor, parte de él sufre una vaporización súbita, proceso que recibe el nombre de cavitación. La cavitación no tendrá lugar, si la suma de las cargas de velocidad y presión en la succión de la bomba, son sensiblemente mayores que la presión de vapor de líquido. Esta diferencia recibe el nombre de carga neta de succión positiva (NPSH), que es igual a: En la práctica, para evitar la cavitación se debe verificar que: NPSH > NPSHr

180

(5.17)

Capítulo 5: Transporte de líquidos

Si no se cumple esta desigualdad, se produce la gasificación total o parcial de líquido en el interior de la bomba, siendo arrastradas las burbujas de gas por el líquido y al llegar a puntos de mayor presión, la fase de vapor condensa nuevamente de forma brusca, implosionando y produciéndose una repentina e importante disminución de volumen, que se traduce en una disminución de caudal, acompañada de ruidos, vibraciones y una gran erosión en las zonas de implosión, principalmente en los bordes de los alabes; todo esto bajo el punto de vista mecánico. Hidráulicamente se observa una importante disminución del rendimiento, al perder capacidad y presión.

Figura 5.29. Fenómeno de cavitación. a) Burbujas implosionando, b) Daño en rodete.

Consideraciones para tener una correcta aspiración Para lograr que en un determinado sistema NPSH sea mayor que NPSHr, debe analizarse la posibilidad de modificar algunos de los siguientes factores: Desnivel de aspiración: Si éste es positivo, al aumentarlo elevamos la carga de aspiración y si es negativo, al disminuirlo reducimos la elevación de aspiración. Diámetro de la tubería de aspiración: Con diámetros mayores se reducen las pérdidas por fricción.

181

Escurrimiento de fluidos

USACH

Longitud de la tubería de aspiración: Al reducirla, se disminuyen las pérdidas por fricción, aumentando el NPSH disponible. Estado superficial de la tubería de aspiración: Las tuberías nuevas o de materiales poco rugosos disminuyen las pérdidas por fricción. Accesorios en la aspiración: La mayor simplicidad en el recorrido disminuye las pérdidas por fricción. Viscosidad cinemática del líquido: Su variación modifica la pérdida de presión por fricción, de tal suerte que al aumentar aquél la disminuye el NPSH disponible. Temperatura del líquido: La presión de vapor aumenta al elevar la temperatura, disminuyendo el NPSH disponible. Flujo: Influye principalmente en la pérdida de carga por fricción, en proporción al cuadrado de las variaciones del mismo. De acuerdo con lo recién descrito, son varios los parámetros que afectan el valor del NPSH disponible, aunque no todos son susceptibles de modificar en una instalación dada. De los factores dependientes del sistema sólo será posible modificarlos en caso de nuevas instalaciones; de los que afectan al tipo de fluido sólo es factible modificar su temperatura, instalando enfriadores previos al bombeo, aunque esto aumenta la pérdida de carga en la instalación y por otro lado, modifica su viscosidad, por lo que deberá meditarse su viabilidad técnica y económica. El flujo se puede reducir recurriendo a la instalación de varias bombas que trabajen en paralelo. También será muy importante tener en cuenta al realizar la instalación que ésta impida la entrada de aire en el sistema, para lo cual se procurará que

182

Capítulo 5: Transporte de líquidos

todos los accesorios y tramos de tubería ajusten perfectamente a fin de que en la tubería no se formen bolsas de aire. La Figura 5.30, muestra algunos detalles de la correcta instalación de una bomba centrífuga.

Figura 5.30. Instalación de bomba centrífuga.

En Figuras 5.31, 5.32 y 5.33 se muestran las curvas características para tres modelos de bombas centrífugas. Dependiendo del catálogo de bombas que se use, puede variar ligeramente la forma de entregar dicha información.

183

Escurrimiento de fluidos

USACH

Figura 5.31. Curvas características para una “familia” de bombas centrífugas (AURORA).

184

Capítulo 5: Transporte de líquidos

Figura 5.32. Curvas características para una “familia” de bombas centrífugas (Wernert Pumpen).

185

Escurrimiento de fluidos

USACH

Figura 5.33. Curvas características para una “familia” de bombas centrífugas (Wernert Pumpen).

Determinación del caudal de operación Para determinar el caudal que entregará la bomba para un determinado sistema de flujo (caudal de operación), debe graficarse la curva de trabajo ∧ requerido del sistema (ΔHr = W /g) vs Q. La intersección de esta curva con la curva ΔH vs Q de la bomba, entregará el caudal de operación como se muestra en Figura 5.34.

186

Capítulo 5: Transporte de líquidos

Figura 5.34. Determinación del caudal de operación (Q1).

Modificación del punto de trabajo El punto de trabajo de la instalación es fijo para unas condiciones dadas y si éste no corresponde al requerido por el proceso debería ser modificado, existiendo tres opciones: 1. Modificar la curva de la bomba Esto se puede conseguir con la misma bomba, variando el diámetro de su impulsor o su velocidad de giro, dentro de ciertos límites. También se puede cambiar la bomba o instalar dos o más bombas en serie o en paralelo. Debe tenerse presente que si se colocan dos bombas en paralelo, el nuevo flujo total descargado no corresponde al doble del entregado por una bomba. En este caso, el caudal de operación será el que corresponda al punto de trabajo con la nueva curva equivalente de las dos bombas. Para obtener la curva equivalente a 2 bombas en paralelo (ambas idénticas), a iguales alturas le corresponde la suma de caudales de las bombas. Así, en el ejemplo que ilustra la Figura 5.35.a, con dos bombas en paralelo (idénticas) el caudal de operación es Q2 y no el doble de Q1 que correspondería a Q3. Es fundamental 187

Escurrimiento de fluidos

USACH

tener presente que para instalar bombas a trabajar en paralelo, sus curvas características han de ser idénticas, pues de lo contrario se produciría un desequilibrio de carga para cada una de ellas. Si el problema es de alturas y no de caudales, se puede recurrir al acoplamiento de bombas en serie, en las cuales la descarga de la primera bomba alimenta la aspiración de la segunda y así sucesivamente. Por tanto, el caudal será el mismo, pero las alturas se suman sucesivamente. En la Figura 5.35.b podemos ver un sistema resuelto por este procedimiento.

Figura 5.35.a. Dos bombas idénticas en paralelo.

Figura 5.35.b. Bombas en serie.

2. Modificar la curva del sistema Consiste en modificar las pérdidas por fricción entre la succión y descarga, lo cual se puede realizar por varios procedimientos: - Instalando accesorios de pérdida de carga variable, como válvulas de regulación. - Cambiando el diámetro de una parte del recorrido (líneas en serie). - Cambiando el diámetro en todo el recorrido. - Colocando una tubería en paralelo a la original, parcialmente o en todo el recorrido.

188

Capítulo 5: Transporte de líquidos

Figura 5.36. Modificación de la curva del sistema al aumentar el diámetro desde φ1 hasta φ2.

En la Figura 5.36 se muestra el efecto de un cambio de diámetro en la tubería. En estos casos hay que pensar en la economía, pues aumentar el diámetro de la tubería conservando la misma bomba puede ser claramente antieconómico, lo mismo que cambiar la bomba conservando la tubería existente. Hay que hacer un balance de costos para analizar que es lo más conveniente. Existe un diámetro económico en cada caso, el que pone en equilibrio el costo de la instalación y el de la energía. La curva de un sistema en paralelo se obtiene aplicando el procedimiento descrito en Capítulo 4, disminuyendo las pérdidas por fricción, al aumentar el área de flujo utilizando tuberías en paralelo, mientras que para un sistema en serie, que consiste en un tendido único, pero formado por varios diámetros de tubería, en éstos se cumplirá que para un determinado caudal la pérdida de carga será la suma de las correspondientes a cada uno de los ramales.

189

Escurrimiento de fluidos

USACH

3. Modificación en bomba y sistema El procedimiento más radical es la unión de los dos ítems anteriores. En efecto, cuando se modifica el punto de trabajo en una bomba hay que tener presente que la potencia absorbida es variable, aumentando con el caudal. En algún caso puede que el motor instalado sea de suficiente potencia para la que se requiere en el punto de funcionamiento, pero que al variar éste puede ser insuficiente. Por este motivo, será necesario verificar cada vez si el nuevo requerimiento de potencia es entregado por el motor disponible, especialmente cuando el nuevo punto está hacia la derecha de la curva, como se indica en la Figura 5.37. Esto mismo ocurre si se cambia el diámetro del impulsor o si por cambio de servicio se ha de bombear líquido de mayor densidad.

Figura 5.37. Curvas carga y potencia vs caudal.

Modificación de las curvas características de la bomba Cuando la bomba es capaz de ser operada a velocidades variables, se obtienen curvas características como las mostradas en la Figura 5.38.a. Por otro lado, lo usual es que las curvas características estén referidas al agua. Las curvas se modifican cuando se bombean fluidos con una viscosidad dis-

190

Capítulo 5: Transporte de líquidos

tinta, aumentando el consumo de potencia y reduciéndose la altura, flujo bombeado y eficiencia (Figura 5.38.b).

Figura 5.38. Curvas características a) Efecto de la velocidad de giro del rodete, b) Efecto del aumento de la viscosidad del fluido.

Para bombas con dimensiones similares pero diferentes diámetro de impulsor, se presentará un cambio en el área de descarga, lo que origina directamente un cambio en la capacidad. El cambio del diámetro también tiene influencia directamente sobre la presión total de descarga, la que aumenta en relación al cuadrado de los diámetros de los impulsores y en la potencia, que aumenta directamente proporcional al cubo de los diámetros. El efecto del diámetro del impulsor o rodete en las características de una bomba centrífuga, se ilustra en la Figura 5.39.

191

Escurrimiento de fluidos

USACH

Figura 5.39. Efecto del diámetro del rodete en las curvas características.

Leyes de semejanza de bombas centrífugas Si llamamos n1 y n2 a las revoluciones por minuto que giran los rodetes de 2 bombas centrífugas iguales, se tienen las siguientes relaciones:

Dichas relaciones pueden ser utilizadas para cambios moderados en la velocidad. Velocidad específica (Ns) Velocidad específica es un concepto que une los tres parámetros más importantes de una bomba centrífuga: capacidad (Q), carga (H) y velocidad de giro del rodete (n), en un solo término. Ns está dado por: 192

Capítulo 5: Transporte de líquidos

H (=) pie Si bien es cierto que Ns podría evaluarse para cualquier condición de operación, la definición de Ns asume que la carga (H) y capacidad (Q), usadas en la ecuación deben corresponder a las condiciones de máxima eficiencia de la bomba. Esto es porque Ns se utiliza para comparar diferentes tipos de bombas centrífugas, especialmente disímiles en los rodetes. Usando leyes de semejanza se puede demostrar que bombas geométricas y dinámicamente semejantes tienen el mismo rendimiento. En la Figura 5.40 se muestran los rangos de velocidades específicas para cada tipo de rodete.

Figura 5.40. Ns para diferentes rodetes (Adaptado de Foust y col., 1961).

En la tabla siguiente se entrega la zona de eficiencia máxima para diferentes tipos de bombas centrífugas.

193

Escurrimiento de fluidos

USACH

Tabla 5.1. Zonas de máxima eficiencia para cada tipo de bomba. Tipo Flujo radial “ “ “ “ “ Flujo mixto “ “ “ Flujo axial “ “

Velocidad específica 500 750 1000 2000 3500 4500 4000 5000 7000 10000 9000 12000 15000

Zona de eficiencia máx. en % 45 – 70 55 – 75 65 – 85 70 – 90 82 – 92 85 – 90 85 – 92 80 – 90 80 – 88 80 – 85 82 – 88 80 – 85 78 – 82

Selección de bombas centrífugas Para elegir una bomba para una determinada aplicación, se deben considerar los siguientes aspectos (normalmente solicitados por los proveedores de bombas): 1. Capacidad requerida (Q). 2. Carga que requiere el fluido, normalmente expresada como “altura” de fluido (ΔH). 3. Potencia requerida por el fluido. 4. Eficiencia del grupo moto-bomba (se debe elegir un grupo moto-bomba que funcione en la zona de mayor eficiencia). 5. NPSH del sistema. Debe ser mayor que el requerido por la bomba.

194

Capítulo 5: Transporte de líquidos

6. Naturaleza del fluido (corrosivo, temperatura, presión de vapor, viscosidad, suspensión, abrasivo). 7. Presión de descarga requerida. 8. Localización de la bomba (espacio disponible, aire libre o no, etc.). 9. Naturaleza de la energía disponible para accionar la bomba.

En la Figura 5.41 se muestra la altura que alcanza una columna de tres líquidos de diferente densidad, siendo la presión de descarga de la bomba de 1 bar, en los tres casos.

Figura 5.41. Efecto de la densidad del líquido, para una misma presión de descarga de la bomba.

B.2. Bombas con efectos especiales Electromagnéticas. Una bomba electromagnética no tiene partes móviles, por lo que no se requiere usar sellos de ningún tipo. Dado que la bomba electromagnética trabaja según el mismo principio que el motor de inducción, sólo metales líquidos, con una conductividad eléctrica alta, pueden ser bombeados con este dispositivo. El metal líquido fluye a través del campo de un electroimán. Una corriente externa pasa a través del metal, en una dirección perpendicular al campo magnético. La fuerza ejercida como resultado de estos campos origina el flujo.

195

Escurrimiento de fluidos

USACH

Empuje a gas. Una forma sencilla y económica para transportes líquidos y que no incluye partes móviles, se realiza mediante el uso de un gas a presión, tal como se muestra en la Figura 5.42. En la ilustración, el líquido se alimenta dentro de un estanque, por gravedad y es fugado hacia afuera por medio de aire comprimido. Este sistema es muy útil para transportar fluidos corrosivos. Figura 5.42. Empuje a gas.

Eyectores. Un eyector es un dispositivo, que usando un fluido impulsor provoca un vacío que se puede utilizar para transportar otro fluido. En Figura 5.43 se muestra un esquema de un eyector.

Figura 5.43. Eyector.

En la literatura existe poca uniformidad con respecto al empleo del término eyector e inyector. Usaremos el término eyector cuando se trate de un dispositivo que utiliza tanto un gas como un líquido como agente motor (fluido impulsor), y que descarga a una presión intermedia entre la presión 196

Capítulo 5: Transporte de líquidos

motriz y la presión de succión. En cambio un inyector es un aparato que utiliza un gas condensable (vapor) para arrastrar a un líquido, al que descarga a una presión mayo que la inicial. Un eductor es un eyector que utiliza un líquido como fluido motor para arrastrar a otro líquido. Las partes esenciales de un eyector son la boquilla (tobera) y el difusor. La boquilla es un dispositivo que transforma la presión elevada del fluido motor en una gran velocidad. El difusor es una cámara de mezcla, y es la zona donde se convierte de nuevo la velocidad final de la mezcla en presión.

Usando la nomenclatura de la Figura 5.43, consideremos que por el tubo de diámetro D circula un fluido, por ejemplo, aire comprimido. Su presión se controla por una válvula de estrangulamiento, no indicada en la figura. Gracias a la depresión que se crea en d el agua sube por la tubería de diámetro D’. Este inyector es, pues, una bomba cuya gran ventaja consiste en carecer de partes móviles.

Despreciando las pérdidas escribamos la ecuación de Bernoulli entre

las secciones 1 y 2: si los puntos 1 y 2 están en el mismo plano horizontal Z1 = Z2, y por lo tanto p1/ρ + V12/2 = p2/ρ + V22/2 ; luego p2/ρ = P1/ρ + V12/2 – V22/2

Usando continuidad (balance de masa):

197

Escurrimiento de fluidos

USACH



Donde Q – caudal de aire que pasa por la tubería D y, por tanto



Introduciendo estos valores en la ecuación de Bernoulli, se obtiene:

Con una válvula (no indicada en la figura) se puede regular P2. Así por ejemplo, al abrir la válvula aumenta Q, con lo que disminuye P2. Aplicaciones de los eyectores. Los eyectores son frecuentemente utilizados para provocar vacío, aunque también son usados como compresores. Dependiendo del vacío requerido pueden ser instalados en serie o en paralelo (si los flujos son muy altos). La Tabla 5.2 muestra los límites de aplicación, de acuerdo con el número de eyectores. Tabla 5.2. Presión absoluta en función del número de eyectores. Presión absoluta (mm Hg)

Número de eyectores en serie

760 – 76

uno

130 – 13

dos

50 – 1

tres

5 – 0.05

Cuatro o cinco

198

Capítulo 5: Transporte de líquidos

Ejemplo 5.1 Una bomba se instala en un sistema para bombear 240 m3/h de agua a 20 ºC, con p1 = 0.5 bar y p2 = 1.1 bar. El diámetro de succión es 150 mm y el de descarga es de 125 mm. La diferencia de nivel entre los puntos 1 y 2 es de 355 mm. Determine la altura de descarga de la bomba. Solución:

La altura es dada por:

Ejemplo 5.2 Una industria química obtiene agua desde un río cercano, debiendo ser tratada y purificada es estanques abiertos, para su posterior distribución a los servicios de la planta. La conducción está construida en cañería de acero comercial de 6(in) de diámetro nominal, cédula 40. El tramo de succión tiene una longitud equivalente de 10(m); la descarga tiene 310(m) de cañería y los accesorios incluidos en este tramo equivalen a 30(m) de cañería de igual diámetro. Si la bomba instalada tiene una curva característica representada por:

199

Escurrimiento de fluidos

ΔH = 200 – 1585 · Q2

USACH

Con: H[=](m) y Q[=](m3/s)

a) ¿Cuál es el caudal que circula por el sistema? b) ¿Cuál es la presión de descarga de la bomba? c) Si el NPSH requerido por la bomba es de 2(m), ¿Cavitará?

Solución: Aplicando un balance de energía mecánica (Bernoulli) entre (1) y (2) se tiene:



200

(a)

Capítulo 5: Transporte de líquidos

Para determinar el caudal de operación de la bomba, debe igualarse ΔHsistema = ΔHbomba. 30 + 4f · 3,34 · 105 · Q2 = 200 – 1585 · Q2

(b)

ε/D = 0,0003 y considerando régimen altamente turbulento a 4f = 0,015 tenemos: 6595 · Q2 = 170 Q = 0,161(m3/s)

4f = 0,0155 Usando la ecuación (b) ⇒ Q = 00,159(m ,159(m 33/s) s) Re = 1,31 · 106 ⇒ 4f Q = 00,0155 ,159(m 3 s ) ∴Q Q == 0,159(m 0,159(m33/s)s ) Por otro lado Pd = Ps + 1000(kg/m3) · 9,8(m/s2) · (200 – 1585 · 0,1592)(m) Pd = Ps + 1,57 · 106(kg/m·s2)

201

(c)

Escurrimiento de fluidos

USACH

Pero, de un balance entre el punto (1) y la succión de la bomba:

Reemplazando este valor en la ecuación(c) se obtiene: Pd = 1576030(kg/m · s2) = 15,55(atm) Para determinar si la bomba cavita o no, debe calcularse el NPSH del sistema:

Dado que NPSHsistema > NPSHrequerido, la bomba no cavitará.

Ejemplo 5.3 Se necesita especificar una bomba para impulsar agua de mar a 20 ºC, la que se empleará como refrigerante en un intercambiador de calor. El flujo de agua necesario es de 300(gpm), con una presión de 30(psia) a la entrada del

202

Capítulo 5: Transporte de líquidos

intercambiador de calor. La cañería es de 4" de acero comercial, tipo 40, con una longitud de 200(m) entre la bomba y el intercambiador. Debe incluirse en ese trayecto 8 codos de 90º y una válvula de compuerta. En el tramo de succión de 12(m), se debe considerar un codo de radio largo de 90º.

Aplicando un balance de energía mecánica (Bernoulli) entre (1) y (2):

(a)



(b)

f ==00,0185 ,00425 ε/D = 0,00045 ⇒ 4f

203

Escurrimiento de fluidos

USACH

Reemplazando los valores en (b), se obtiene:

Introduciendo estos valores en la ecuación (a): ∧

337,6(m/s) ·⋅ gg ⇒ ΔH f H=sistema 0,00425 ==34,5(m) W == −–337 ,6(m s ) 2 = ΔH ∆Hsistema ∆ 34 ,5(m ) üüü üüü 2

Esta carga o energía, requerida por el sistema. ∧ Potencia Potencia= W ⋅·QQ⋅ ρ · ρ==−–8,8,6(hp) 6(hp ) (Potencia requerida por el fluido)

(c)

La presión de succión se obtiene de un balance entre el punto (1) y la succión de la bomba:

(d)

Introduciendo estos valores en la ecuación (c):

De acuerdo con la información para bombas centrífugas (curvas características), se debe seleccionar la bomba centrífuga más adecuada.

204

Capítulo 5: Transporte de líquidos

EJERCICIOS 5.1. Se necesita especificar una bomba para impulsar agua de mar a 20 ºC, la que se empleará como refrigerante, en un intercambiador de calor. El flujo de agua necesario es de 500 gpm, con una presión de 50 psia a la entrada del intercambiador de calor. La cañería es de 4" de acero comercial, tipo 40, con una longitud de 150 m entre la bomba y el intercambiador. Deben incluirse en ese trayecto 6 codos de 90º, una válvula de compuerta 50% abierta y una placa orificio con β = 0.5. En el tramo de succión de 12(m), se debe considerar un codo de radio largo de 90º.

5.2. Hasta el circuito de refrigeración de una serie de compresores, que se emplea para obtener nitrógeno líquido a partir de aire, se bombeará agua desde un pozo. Por consideraciones técnicas, especificadas por los representantes del equipo de compresión, la presión a la entrada del circuito de refrigeración debe ser de 2(kgf/cm2), con un caudal mínimo de agua de 0,002(m3/s). Para bombear el agua se dispone de varias bombas, todas ellas iguales y con las mismas curvas características. Estas son:

205

Escurrimiento de fluidos





USACH

Q(m3/s)

H(m)

NHPSrequerido(m)

ηGMB%

0

45

-

50

0,0002

44

1,2

55

0,0004

43

1,3

60

0,0006

42

1,4

60

0,0008

40

1,6

65

0,001

38,5

1,8

70

0,0012

35

2,1

75

0,0014

31

2,4

75

0,0016

25

2,8

67

0,0018

19

3,3

60

0,002

12

3,8

55

El diámetro de la tubería en el sistema es de acero comercial tipo 40 de 1" y el plato orificio representado en la figura, tiene un diámetro de 1.5 cm. Considerando que la temperatura del agua es de 15 ºC, responda: a) ¿Cuántas bombas serán necesarias? ¿En serie o en paralelo? (R: 2 en serie). b) ¿Cuál será el caudal de operación? (R: Qoperación = 0,0014(m3/s); ΔHoperación = 58(m)). c) ¿Cuál será la potencia hidráulica en HP, en las condiciones de operación? d) ¿Cuál será la potencia eléctrica en (KW) en las condiciones de operación? e) ¿Cuál es la desventaja de utilizar varias bombas en vez de una sola? f) ¿Se producirá cavitación? g) ¿Qué lectura en (cmHg) se leerá en el manómetro, en las condiciones de operación?

206

Capítulo 5: Transporte de líquidos



h) Explique cómo determinaría la fracción a cerrar en la válvula de globo, si se desea trabajar con el mínimo de agua requerida.

5.3. La curva característica de una bomba centrífuga está dada por: H= 35 – 0.1 · Q2 , con: H[≡](m) y Q[≡](l/s) La bomba está instalada en una red de cañerías de 2" y 30(m) de longitud. Los accesorios de la red significan además una longitud equivalente de 30(m). La bomba impulsa agua a 20 ºC desde un estanque a otro, ambos abiertos a la atmósfera, estando el estanque receptor a tres metro sobre el primero. a) Graficar las curvas características de la bomba y del sistema. b) ¿Qué caudal circula por el sistema? (R: 9.6 l/s). c) Si se instalara una placa orificio que produce una pérdida de carga adicional equivalente a 10(m). ¿Cuál será el nuevo caudal? p1 = p2 dN = 2(in) ⇒ dfi ==0,0525(m) 0,00425 z2 – z1= 3(m) Lequibalente(accesorios) = 30(m) Lcañería = 30(m) Lequibalente(placa orificio) = 10(m) ρH2O = 999(kg/m3) µH2O = 1 · 10–3(kg/m · s)

207

Capítulo 6: Flujo de gases

CAPÍTULO 6 FLUJO DE GASES

En el flujo de gases por tuberías, si la variación de presión es significativa, se produce una variación de la densidad del gas, que inhabilita el uso de la ecuación de Bernoulli. Una situación especial son los ductos de acondicionamiento de aire, donde los bajos flujos y caídas de presión permiten seguir considerando al aire como fluido incompresible, y por lo tanto emplear la ecuación de Bernoulli. En el flujo de fluidos compresibles debe considerarse que a lo largo de la tubería la densidad del gas disminuye (proporcional a su presión), con lo cual su velocidad aumenta a expensas de su energía interna (o sea baja su temperatura), todo lo cual dificulta la integración del Balance de Energía Mecánica. Para el estudio del escurrimiento de fluidos compresibles suelen plantearse las siguientes consideraciones: 1) Estado estacionario. 2) Flujo unidimensional. 3) Sin perfiles radiales de presión ni velocidad. 4) Fricción solo en las paredes del ducto. 5) Trabajo mecánico nulo (Para su evaluación debe usarse la ecuación correspondiente aplicada a la entrada y salida del equipo impulsor). 6) Energía potencial despreciable. En el estudio del flujo de gases, además del Balance de Energía Mecánica, se utiliza el Balance de Energía Total (si el flujo es adiabático), la ecuación de continuidad y una ecuación de estado. No obstante lo anterior,

209

Escurrimiento de fluidos

USACH

es posible aplicar el Balance de Energía Mecánica a un ducto de diámetro D constante y de longitud dL, tal como muestra la figura siguiente:

Figura 6.1. Sección de tubería de longitud dL.



(6.1)



(6.2)

Introduciendo la densidad de flujo másico, G = w/A = ρV, la cual permanece constante en la tubería, e integrando, se obtiene:

6.1.

(6.3)

Flujo isotérmico

La integración del término que incorpora la densidad del gas requiere el empleo de una ecuación de estado. Considerando gas ideal (ρ = p · M/RT), y que la temperatura permanece constante [Flujo Isotérmico], se obtiene:

210

(6.4)

Capítulo 6: Flujo de gases

Además, usando la definición de G y la ecuación de gas ideal, se tiene que V2/V1 = p1/p2, con lo cual:

(6.5)

Si Δp/p1 ≤ 0.20 es posible eliminar el término ln (p1/p2), obteniéndose la ecuación de Waymouth: (6.6) Despejando G desde la ecuación 6.5, es posible analizar qué sucede con la densidad de flujo másico, cuando se conectan dos gasómetros a diferentes presiones. En este ejercicio virtual consideremos que p1 permanece constante y que gradualmente es posible disminuir p2, desde un valor p1 hasta cero (Figura 6.2).

Figura 6.2. Representación esquemática de la conexión de dos gasómetros.

211

(6.7)

Escurrimiento de fluidos

USACH

En primer lugar es posible observar que la disminución de p2 afecta tanto el numerador como el denominador de la ecuación 6.7. Por un lado, tiende a disminuir G, mientras que por otro, su tendencia es incrementarlo. Esto nos indica matemáticamente que dicha expresión presenta un valor máximo. Este máximo se obtiene derivando G con respecto a p2, manteniendo los demás parámetros constantes, obteniéndose la expresión siguiente, donde p2C, es la presión que produce el flujo máximo, Gmáx.

(6.8)

Figura 6.3. Representación esquemática de la variación del flujo másico versus la presión de descarga (p2).

La comparación de la predicción de la ecuación con valores experimentales muestra, que si bien es cierto sus predicciones son razonables para presiones entre p2C y p1, bajo p2C, el flujo no disminuye, como lo predice dicha ecuación, sino que éste permanece constante. La existencia de Gmáx puede explicarse usando conceptos de física. Consideremos una tubería conteniendo un gas y que repentinamente un

212

Capítulo 6: Flujo de gases

pistón se mueve a la izquierda, generando una subida de presión en dicho sector. Esta región de mayor presión, tal como una onda sonora, viaja a la velocidad del sonido (VS). Por delante del frente de onda, la velocidad del fluido es cero y detrás de él hay una velocidad V+dV, provocada por el incremento de presión dp. Un observador en posición estacionaria observaría que se le aproxima una onda de presión a una velocidad V y que se aleja a V+dV, siendo dV 1000), esta diferencia no supera el 5%. Por esto, es corriente hacer estimaciones utilizando las ecuaciones para flujo isotérmico, dada su facilidad de manejo Cálculo de trabajo mecánico y potencia

Entre los equipos impulsores de gases se tienen sopladores, compre-

sores y ventiladores. El trabajo necesario se obtiene aplicando un balance de energía mecánica a la entrada y salida del equipo impulsor.

Figura 6.8. Balance de energía mecánica a equipo impulsor.

En esta situación son despreciables las variaciones de energía cinética y potencial. Por otro lado las pérdidas por fricción en el paso del gas por el equipo impulsor, serán consideradas en la eficiencia del equipo y por lo tanto no evaluadas en forma separada. De acuerdo con esto, el balance de energía mecánica se reduce a: (6.19) Suponiendo un comportamiento ideal del gas y que la compresión es reversible y adiabática, el trabajo ideal efectuado sobre el gas esta dado por: (6.20)

218

Capítulo 6: Flujo de gases



Luego, la potencia teórica requerida se calcula como: (6.21)

Donde n corresponde a moles/tiempo, k = Cp/Cv Durante la pasada del gas por el equipo impulsor, la fricción produce el calentamiento del gas, debido a lo cual la potencia real consumida es mayor que la potencia ideal, relacionándose ambas a través de la eficiencia del equipo impulsor. (6.22)

6.3.

Equipos impulsores de gases

Existe una similitud entre los equipos impulsores de líquidos descritos en el capítulo 5 y los equipos utilizados para impulsar gases. Estos son clasificados, de acuerdo con el nivel de sobrepresión que entregan, en compresores, sopladores y ventiladores. Los compresores de gases son subdivididos en compresores de desplazamiento positivo y compresores centrífugos. A su vez, entre los de desplazamiento positivo se distinguen los compresores recíprocos (con pistones) y los compresores o sopladores rotatorios. Compresores recíprocos. Estos pueden entregar el gas con sobrepresiones desde algunas psi hasta presiones mayores a 35000 psi. En estos equipos el gas que va a ser comprimido, entra y sale del cilindro, a través de válvulas diseñadas para actuar cuando la diferencia de presión entre el contenido del cilindro y las condiciones exteriores es la deseada. Si se usa un compresor con varias etapas, generalmente se debe enfriar el gas entre ellas (hasta 12 etapas). Con una etapa se consiguen presiones de descarga de 150 psi, su-

219

Escurrimiento de fluidos

USACH

biendo a 750 psi con dos etapas y 3000 psi con tres etapas. Estos equipos presentan eficiencias entre 65 y 80%. Compresores o sopladores rotatorios. Al igual que los equipos impulsores de líquidos, existe una gran variedad de elementos rotantes, que entregan características específicas a cada equipo: de engranajes, de tornillo, de lóbulos, de aletas deslizantes, etc. Se caracterizan por entregar una descarga prácticamente continua.

La selección de los compresores se efectúa generalmente mediante catálogos, en los cuales se indican para diversas capacidades y presiones de descarga, los tipos más adecuados, junto con datos de tamaño y potencia de los motores. Compresores centrífugos. Consisten en uno o más rodetes montados en un mismo eje, dentro de una carcasa. En un compresor centrífugo se aumenta la presión del gas, cuando éste circula a través del impulsor o rodete, en forma análoga a lo que sucede en una bomba centrífuga (Figura 6.8). Son utilizados tanto para impulsar un gas como para ejercer vacío. Normalmente entregan sobrepresiones inferiores a 9 psi, aunque diseños especiales entregan presiones de salida hasta de 800 psi.

Figura 6.9. Compresores centrífugos.

220

Capítulo 6: Flujo de gases

Ventiladores. Los ventiladores se utilizan para transportar gases a bajas presiones, generalmente bajo 40 cm col de agua (~0.6 psi). Su principio de funcionamiento es semejante a un compresor centrífugo, aunque al operar a bajas presiones, los volúmenes de succión y descarga son similares, siendo posible despreciar el efecto de compresibilidad sobre el gas. De acuerdo con la dirección del flujo de aire, los ventiladores se clasifican en ventiladores de flujo radial y ventiladores de flujo axial. Ambos tipos se utilizan para trabajos de ventilación, suministro de aire a quemadores y calderas, suministro de aire para secado, transporte neumático de sólidos, eliminación de humos, etc. En términos generales, la eficiencia de funcionamiento de los ventiladores se encuentra entre 40 a 70%. La presión operacional es la suma de la presión estática y la carga de velocidad del aire que sale del ventilador. Normalmente se expresa en cm de columna de agua.

El volumen de gas desplazado varía directamente con la velocidad de giro del impulsor (rodete). La presión estática cambia en relación con el cuadrado de la velocidad del rodete y la potencia varía en relación al cubo de la velocidad. En la práctica, la eficiencia final de un ventilador (depende fuertemente del tipo de álabe), debe ser determinada experimentalmente. El tamaño de los ventiladores va desde pequeñas unidades, para circular aire en vehículos, hasta enormes ventiladores industriales, diseñados para mover hasta 60 m3/s a 15 cm col de agua. Los ventiladores de aleta radial son recomendados para el transporte de gases que contienen sólidos en suspensión, debido a que la acción centrífuga tiende a mantener limpias las aletas. Su eficiencia es moderada. El ventilador curvado hacia delante es apropiado para baja velocidad y grandes volúmenes. Su eficiencia es media y es adecuado para gases limpios. Los ventiladores con las aletas curvadas hacia atrás se caracterizan por una mayor eficiencia mecánica (cercana al 90%) y son adecuados parta gases limpios.

221

Escurrimiento de fluidos

USACH

Los ventiladores axiales son utilizados en aplicaciones de ventilación, donde se requieren bajas presiones. Se instalan en el interior del mismo ducto o en ventanas y/o paredes. Determinación del caudal de operación. Para determinar el caudal de operación debe trazarse la curva de resistencia del sistema en función del flujo de aire. La intersección de ambas curvas define el caudal de operación, tal como se esquematiza en la Figura 6.10.

Figura 6.10. Determinación del punto de operación y efecto de la velocidad de giro en las curvas características.

Ejemplo 6.1 Se debe transportar hidrógeno (M=2 kg/kmol) desde un estanque a 20(atm) hasta otro estanque a igual presión. Para realizar la operación es necesario elevar la presión hasta 25(atm) a la salida del primer estanque, mediante una bomba. La tubería de conducción es de acero de 2" y su longitud es de 500(m). Calcular el flujo de gas y la potencia de bombeo.

222

Capítulo 6: Flujo de gases

Solución: A continuación se entregan una serie de datos: di = 52,5(mm) ε/D = 0,0009 T = 20 ºC ⇒ μf H2= =0,900425 · 10–5 (poise) Como Δp/p1 = (25 – 20)/25 = 20%

Se usará la ecuación de Weymouth donde:

(1) El método de cálculo consiste en suponer un valor de f *, con ε/D determinar el Re del gráfico f * v/s Re. Con la ecuación (1) calcular G y el Re = G · D/μ hasta obtener Re aproximadamente constante. f* 0,0250 0,0220 0,0200 0,0195

Re 3,5 ·104 9,0 ·104 3,5 ·105 6,0 ·105

G, (kg/cm2 · s) 8,975 ·10–3 9,567 ·10–3 10,030 ·10–3 10,160 ·10–3

Re 5,24 ·105 5,58 ·105 5,85 ·105 5,93 ·105

Se considerará el valor real G = 101,6 (kg/cm2 · s) w = G · A = 101,6 · 21,6 · 10–4 = 0,219(kg/s)

Para calcular la potencia real se requiere disponer de la eficiencia del equipo impulsor.

223

Escurrimiento de fluidos

USACH

Ejemplo 6.2 Un gasómetro contiene metano (M =16 kg/kmol)) a 5(atm). El gas debe llevarse hasta su lugar de uso a través de una tubería recta de 250(m)de acero comercial, tipo 40 de diámetro nominal 2" (Di = 2,067(in)). Determinar la presión en el lugar de llegada, para que el flujo sea de 4/5 del máximo. Considerar temperatura constante de 20 ºC. Solución:

μ = 0,017(cp) (ε/D) = 8,71 · 10–4

Suponiendo flujo muy turbulento, f = 0,0047 De ln(p1/p2C)2 + 1 – (p1/p2C)2 + 4f · (L/D) = 0 L/D = 4761,73 ln(X2) – X2 = – 4 · 0,0047 · 4761 – 1 = – 90,52 X 10 9 9,7 9,75

ln(X2) – X2 – 97,69 – 76,60 – 89,55 – 90,51

X = 9,75 = p1/p2C p2C = 5/9,75 = 0,513(atm)

Gmax = 133,13(kg/m2 · s)

Re = G · D/μ = 6,5 · 105 (Altamente turbulento) f = 0,0047 (Suposición correcta) G = 4/5 · Gmax = 106,5(kg/m2 · s) Re = 3,28 · 105 Del gráfico f v/s Re, f = 0,0049

224

Capítulo 6: Flujo de gases



Reemplazando en:

p2(atm) φ(p2)

1 4,57

4 47,71

3 27,22

3,3 32,74

3,29 32,546

p2 = 2,39(atm) Ejemplo 6.3 En un estanque se mantiene aire a 20 ºC y 5 kgf/cm2. Determine el flujo de aire que saldrá a la atmósfera, cuando se descarga a través de una tubería de longitud 7,9 m y de DN 2" con 3 codos de 90º. El extremo del tubo penetra ligeramente en el depósito. Solución: Є/D = 0.0004; M = 29 kg/kg-mol ; γ = k= Cp/Cv = 1.4 Suponiendo f* = 4 f = 0,016 Resistencias

L/Dequivalente

N = f * L/D

Contracción y expansión en el tubo Tubo recto 3 Codos Total:

~72 150 90 354

1.15 2.41 + 1,44 5,00

225

Escurrimiento de fluidos

USACH

Usando la Figura 6.6 se obtiene: G/Gnz*=0.5  G = 712,9 (kg/m2s) Luego, el flujo másico de aire será:

Ahora debe verificarse si efectivamente nos encontramos en la zona alta turbulencia:

En este caso la suposición considerada es correcta. Si el valor de f * es inadecuado, deberá repetirse el cálculo.

226

Capítulo 6: Flujo de gases

EJERCICIOS 6.1. Una tubería que conduce gas natural tiene un diámetro interior de 36". Las estaciones compresoras están distanciadas 60 millas, siendo la presión de salida de la primera 750 psia y la presión de entrada de la segunda 500 psia. El gas puede ser considerado como gas ideal de peso molecular 16. Si la temperatura se mantiene constante en 70 ºF. a) ¿Cuál es flujo másico? b) ¿Hasta qué valor se puede bajar la presión a la entrada de la segunda estación? c) Manteniendo el valor de 750 psia a la salida de la primera estación, ¿Hasta qué valor puede aumentarse como máximo el flujo gaseoso? 6.2. Se bombea hidrógeno a 20 ºC, desde un estanque a 20 atm de presión (absolutas) a través de 1600 pies de cañería horizontal de 2" de diámetro nominal. En la descarga, la presión también es de 20 atm. La presión es elevada a 25 atm. por una bomba ubicada a la salida del estanque. Si en las condiciones del transporte la viscosidad es de 0.009 cp. ¿Cuál es el flujo transportado? ¿Cuál es la potencia de la bomba? 6.3. Nitrógeno fluye desde una tubería de 2" en la cual la temperatura y presión son de 40 ºF y 40 psig, hacia una cañería de 1" en la cual la presión es de 21,3 psig. Calcule la velocidad en cada tubería, suponiendo condiciones isotérmicas y que no se producen pérdidas. 6.4. Un gasómetro contiene SO2 a la presión de 5 atm. El gas ha de llevarse hasta el lugar de aplicación a través de una tubería de 2" y longitud total de 250 m. Determine la presión en el lugar de aplicación, para que el flujo sea 2/3 del máximo, considerando flujo isotérmico.

227

Escurrimiento de fluidos

USACH

6.5. El hidrógeno empleado en una planta de síntesis de amoniaco ha de entrar en los convertidores a 75 atm. Si el gasómetro dispone de hidrógeno a 90 atm. y la línea de conducción tiene una longitud de 220 m. Determine el diámetro de la tubería, si el flujo másico debe ser de 60 kg/min, en condiciones isotérmicas a 27 ºC.

228

Capítulo 7: Caracterización de partículas sólidas

CAPÍTULO 7 CARACTERIZACIÓN DE PARTÍCULAS SÓLIDAS

Las propiedades físicas y fuerzas de atracción que existen entre partículas sólidas, caracterizan el escurrimiento de fluidos a través de ellas. Entre las propiedades físicas, deben mencionarse la densidad, el tamaño y forma de las partículas. Con respecto a la densidad, debe diferenciarse entre la densidad de la partícula (llamada también densidad real), ρs y la densidad aparente (“bulk”), ρb. Esta última corresponde a la densidad global del material disperso, incluyendo los espacios libres entre partículas. Su evaluación experimental es simple, utilizándose por ejemplo con un picnómetro o una probeta graduada. La relación entre ambos parámetros es a través de la porosidad, ε, la cual se define como: (7.1) De esta forma:

ρb = ρs (1 - ε)

(7.2)

Valores frecuentes de porosidad en lechos fijos constituidos por partículas se encuentran entre 0.35 y 0.5. En aplicaciones en lechos fluidizados, estos valores suben a porosidades sobre 0.6 y son superiores a 0.85 en transporte neumático.

7.1.

Tamaño de partículas

El tamaño de la partícula es una descripción de su extensión en el espacio, la cual puede describirse en forma lineal (diámetro, ancho, alto, espesor), superficial o de volumen. El concepto de diámetro, originalmente definido para describir el tamaño de partículas esféricas, también puede

229

Escurrimiento de fluidos

USACH

ser empleado para caracterizar partículas irregulares. El diámetro de éstas últimas puede ser definido en términos de la geometría de partículas individuales, existiendo varias alternativas al respecto. Una de ellas considera el diámetro de la partícula como un promedio aritmético de sus dimensiones lineales características (di = (l + b + h)/3). Otra definición utiliza el promedio geométrico de las dimensiones lineales (di = (l b h)1/3). Un diámetro recomendado para emplear con partículas de pequeño tamaño (< 40µm) es el diámetro de Stokes, Dst, el cual esta basado en la ley de Stokes: (7.3) donde µ es la viscosidad del fluido, V es la velocidad de descenso de la partícula de densidad ρs, en el fluido de densidad ρf. Aunque teóricamente la Ley de Stokes es válida solo para partículas esféricas, la diferencia entre el volumen de una partícula irregular y el volumen de una esfera equivalente, no representa un error, sino que por el contrario, da información útil sobre la forma de la partícula. Para la determinación experimental del tamaño de mezclas de partículas irregulares, entre 7.6 mm y 0.038 mm, lo usual es emplear una serie de tamices, cuyas mallas de alambres y aberturas están cuidadosamente estandarizadas. Los tamices se identifican por el número de mallas por pulgada, debiendo considerarse que las aberturas son menores a causa del espesor de los alambres. Las series de tamices más usadas son la Sieves y la Tyler. Esta última se basa en la abertura del tamiz de 200 mallas, la que se establece en 0.074 mm. El área de las aberturas, en cualquier tamiz de la serie, es exactamente el doble de las aberturas del tamiz inmediato inferior. La relación de la dimensión real de las mallas de cualquier tamiz a la del inmediato inferior es √ 2. Para análisis especiales es posible disponer de tamices intermedios.

230

Capítulo 7: Caracterización de partículas sólidas

Tabla 7.1. Abertura de tamices Tyler. # malla 3 4 6 8 10 14 20 28

Abertura, µm 6680 4699 3327 2362 1651 1168 833 589

# malla 35 48 65 100 150 200 325 400

Abertura, µm 417 295 208 147 104 74 53 38

El análisis de una muestra de sólidos se realiza colocando la muestra (alrededor de 200 gramos, aunque varía según la densidad de las partículas) en el tamiz superior, y agitando la serie de tamices seleccionados durante un tiempo definido (entre 10-25 minutos). Luego se pesan las fracciones retenidas en cada tamiz y se informan los % o fracciones respecto a la masa total inicial. Las partículas que pasan por el tamiz más fino se recogen sobre un colector colocado en el fondo de la serie de tamices. Los resultados de un análisis por tamizado se tabulan o grafican, para indicar la fracción másica sobre cada tamiz, en función del intervalo del tamaño de malla entre dos tamices. Puesto que las partículas retenidas en cualquier tamiz provienen del tamiz inmediatamente superior, se necesitan dos números para especificar el tamaño de la fracción retenida entre dos tamices consecutivos, uno para el tamiz a través del cual pasan las partículas y otro para el tamiz donde son retenidas. Así la notación 14/20 significa que las partículas pasaron el tamiz de 14 mallas y quedaron retenidas en el de 20 mallas. También se usa la notación – 14 + 20, con igual interpretación. La dispersión de una mezcla de partículas es frecuentemente representada por una curva de distribución de tamaño de partículas, en la cual se grafica el diámetro de las partículas vs la fracción (o porcentaje) de partículas de tamaño i. Esta representación corresponde a la curva de distribu-

231

Escurrimiento de fluidos

USACH

ción diferencial (se conoce también como análisis de frecuencia). También se emplea la distribución de partículas en forma acumulativa (integral), en cuyo caso se grafica di vs fracción (o porcentaje). En el análisis acumulativo (Figura 7.1) se suman los incrementos diferenciales individuales, comenzando por el retenido en el tamiz superior y tabulando o representando las sumas acumuladas frente a la dimensión de malla del tamiz que retiene la última fracción acumulada o añadida. En cálculos basados en el análisis diferencial (Figura 7.2), se considera que todas las partículas de una fracción tienen igual tamaño, y que éste es la media aritmética de las dimensiones de malla de los tamices que definen la fracción. En algunas aplicaciones se emplea la media geométrica para obtener el tamaño representativo.

Figura 7.1. Análisis acumulativo por tamizado.

232

Capítulo 7: Caracterización de partículas sólidas

Figura 7.2. Análisis por tamizado.

Funciones empíricas de distribución. Frecuentemente es necesario representar las distribuciones de tamaño de partículas a través de diversos modelos empíricos, con el fin de facilitar el manejo y/o análisis de la información. Entre los modelos más usados están los de Gates-Gaudin-Schumann y el de Rosin-Rammler-Bennett. El primer modelo corresponde a una representación del tipo: (7.4) donde K y m son parámetros del modelo. Yi representa la fracción acumulada pasante en el tamiz i, es decir la fracción de partículas con tamaño menor a dpi. El modelo (o función) de Rosin-Rammler-Bennett, tiene la forma: (7.5)

donde la representación adecuada de los valores experimentales permite determinar los parámetros xo y m. Llamando Ri a la fracción acumulada rete233

Escurrimiento de fluidos

USACH

nida, la cual representa la fracción de partículas con tamaño mayor a dpi, definida como: Ri = 1 – Yi, se tiene: (7.6)



Aplicando logaritmo dos veces a esta ecuación, y representando el resultado en un gráfico log-log vs log (papel Rosin-Rammler) se obtiene una recta de pendiente m. Si se introduce dpi = xo, entonces Ri = 1/e. El tamaño representativo de mezclas de sólidos irregulares, debe ser caracterizado a través de un diámetro promedio: Dprom = d1G1/G + d2G2/G + ···· dn Gn/G = Σ di Gi/G (7.7) donde di corresponde a la partícula de diámetro di, mientras que Gi representa el peso de cada fracción i. Otras formas para expresar el diámetro de mezclas, consisten en el empleo de promedios geométricos. Un tamaño representativo muy utilizado lo constituye el diámetro superficie/volumen, conocido como diámetro de Sauter, dps: (7.8)

Otra alternativa es el empleo de un diámetro medio volumétrico, dv, cuya expresión para una mezcla constituida por partículas de igual factor de forma (volumétrico) está dada por:

(7.9)

234

Capítulo 7: Caracterización de partículas sólidas

7.2.

Forma de las partículas

Generalmente la forma de partículas irregulares se expresa como una aproximación a una esfera, utilizándose el concepto de esfericidad, φs, definido como: (7.10) Si dp representa el diámetro de una esfera de igual volumen que la partícula, Vp, entonces (7.11) Si Sp representa la superficie de la partícula, entonces:

(7.12) Por otro lado, la superficie específica, av, se define como: av = Sp/Vp. Si la partícula es de forma esférica, entonces, se obtiene: av = 6/dp, siendo dp el diámetro de una esfera de igual volumen que la partícula. Esta ecuación permite definir un diámetro equivalente, Dp, para partículas no esféricas. Por lo tanto, en general: (7.13)

235

Escurrimiento de fluidos

USACH

Tabla 7.2. Esfericidad de algunas partículas. Tipo de partícula Esfera Cubo Cilindro (h = d) Cilindro ( h = 10 d) Carbón activado y sílica gel Carbón calidad antracita Carbón pulverizado Láminas de mica Arena redondeada Arena de playa Arena de río Trigo

Esfericidad (φ) 1.0 0.81 0.87 0.47 0.70 – 0.90 0.63 0.73 0.28 0.86 < 0.86 > 0.53 0.85

El número de partículas de una muestra de masa m, está dado por: (7.14) Si llamamos S a la superficie total de las partículas, entonces, empleando la definición de esfericidad, se tiene: (7.15)

Ahora bien, si se dispone de una mezcla de partículas, constituida por n fracciones de masa total m y cada una de las fracciones con masas m1, m2, mn y densidades ρ1, ρ2 ρn, esfericidades φ1, φ2, φn, y tamaño medio en cada fracción dp1, dp2, dpn, y siendo xi =mi/m, la fracción másica de la

236

Capítulo 7: Caracterización de partículas sólidas

fracción i de la mezcla, la superficie específica de la mezcla por unidad de masa, Sm, será: (7.16)

Ejemplo 7.1 La tabla siguiente resume valores de un análisis por tamizado de 117.431 gramos de partículas de aserrín de madera, obtenidos con tamices con aberturas entre 4.00 mm y 0.090 mm. Abertura tamiz [mm] 4.000 3.360 2.830 2.360 2.000 1.400 1.000 0.710 0.600 0.500 0.356 0.217 0.150 0.090 0.000 TOTAL

mret [g] 1.5598 2.2786 3.4165 5.0329 3.2621 17.2481 46.9982 25.5679 4.2913 3.9841 2.6114 0.8997 0.1569 0.0859 0.0376 117.431

a) Represente gráficamente la distribución en forma diferencial y acumulativa. b) ¿Qué función representa mejor la distribución de la galena? ¿Cuáles son sus parámetros? c) Determine el diámetro medio de Sauter, promedio y volumétrico de la mezcla.

237

Escurrimiento de fluidos

USACH

Solución: a) La siguiente tabla muestra los valores necesarios para graficar la distribución, tanto diferencial como acumulativa. A continuación se muestra un ejemplo de cálculo para la primera corrida. dp = Abi = 4.000(mm), para el resto de la tabla se utiliza la fórmula

dpi [mm] 4,000 3,680 3,095 2,595 2,180 1,700 1,200 0,855 0,655 0,550 0,428 0,287 0,184 0,120 0,045

xi 0,013 0,019 0,029 0,043 0,028 0,147 0,400 0,218 0,037 0,034 0,022 0,008 0,001 0,001 0,000

238

Yi 0,013 0,033 0,062 0,105 0,132 0,279 0,680 0,897 0,934 0,968 0,990 0,998 0,999 1,000 1,000

Capítulo 7: Caracterización de partículas sólidas

239

Escurrimiento de fluidos

USACH

b) Modelo Gates-Gaudin-Schuman (G-G-S):

Aplicando logaritmos: ln(1 – yi) = – m · ln(k) + m · ln(dpi)

ln(1 – yi) = 2,4138 – 2,4182 · ln(dpi) → R2 = 0,9531 m = 2,4182 – m ∙ ln(k) = 2,4138 → k = 2,713 Modelo Rosin-Rammler-Bennett (R-R-B)



→ R2 = 0,9844

m = 2,8411 m ∙ ln(x0) = 1,9268 → x0 = 1,970 El modelo de R-R-B representa mejor la distribución de galena, ya que presenta un mayor valor del parámetro de ajuste R2. c)

240

Capítulo 7: Caracterización de partículas sólidas

EJERCICIO 1.

Galena triturada proveniente de un molino de bolas, presenta la siguiente distribución granulométrica: Mallas Tyler -28+35 -35+48 -48+65 -65+100 -100+150 -150+200 -200



Esfericidad, φ 0.85 0.85 0.86 0.87 0.88 0.90 0.93

Fracción en peso 0.150 0.200 0.171 0.134 0.104 0.080 0.161

a) Represente gráficamente la distribución en forma diferencial y acumulativa. b) ¿Qué función representa mejor la distribución de la galena? ¿Cuáles son sus parámetros? c) Determine la superficie específica y el diámetro medio de Sauter de la mezcla.

241

Capítulo 8: Factores de fricción en medios porosos

CAPÍTULO 8 FACTORES DE FRICCIÓN EN MEDIOS POROSOS

En diversas aplicaciones industriales se produce el flujo de un fluido a través de un medio poroso (Figura 8.1): filtración, reacciones químicas con uso de catalizadores, columnas de relleno, secado de cereales, etc. Con respecto a la determinación de la caída de presión a través del lecho de partículas, se distinguen tres métodos: en uno de ellos se considera la columna de relleno como un manojo de tubos enmarañados de sección irregular, a los cuales se adapta el desarrollo de la teoría aplicada al flujo en el interior de un ducto. En otro método, se representa la columna de relleno como un conjunto de objetos sumergidos, calculándose la pérdida de presión como la suma de las resistencias de las partículas sumergidas. Finalmente, existe también un método teórico que aplica la mecánica del medio continuo (ecuación de continuidad y la ecuación del movimiento) a la matriz porosa.

8.1.

Teoría del conjunto de tubos

En la teoría del conjunto de tubos, el relleno puede estar constituido por partículas de diferentes formas (Figura 8.2), debiendo aceptarse las siguientes consideraciones: a) Relleno uniforme y sin canalizaciones. b) Diámetro de la columna, Dc >> D partículas. c) Diámetro constante de la columna. d) No existe efecto de pared.

243

Escurrimiento de fluidos

USACH

Figura 8.1. Columna de relleno.

Figura 8.2. Materiales típicos de columnas de relleno.

Para el flujo por el interior de tuberías se definió el factor de fricción (Ecuación de Fanning). En forma análoga para columnas de relleno se define: (8.1) donde el Dp equiv corresponde al diámetro equivalente de las partículas (Dp = φs dp), siendo dp el diámetro de una esfera de igual volumen que la partícula. vo es la velocidad superficial, es decir, la velocidad lineal que tendría el flui244

Capítulo 8: Factores de fricción en medios porosos

do en la columna sin relleno. P corresponde al potencial fluido-dinámico y L es la longitud de la columna de relleno (vo = Q/Acolumna). La velocidad real del fluido, v, se relaciona con vo a través de la porosidad del lecho, ε: v = vo/ε (8.2) Para describir el flujo a través de medios porosos, es necesario definir un diámetro equivalente de la columna de relleno, De: (8.3)

Por otro lado, para el escurrimiento en régimen laminar de un fluido

newtoniano, por el interior de un ducto circular, es utilizable la ecuación de Hagen-Poiseuille: (8.4)

Aplicando esta ecuación para la columna de relleno e incorporando la definición de De =Dc = 4 RH se obtiene: (8.5) Re-ordenando, utilizando la definición de radio hidráulico y de Dp, se obtiene: (8.6) En esta ecuación, el empleo del concepto de radio hidráulico y el hecho de que la longitud real que recorre el fluido en el interior del lecho sea mayor que L, se traduce en que sus predicciones deban ser corregidas, a fin de representar adecuadamente la situación real. En efecto, sus predicciones

245

Escurrimiento de fluidos

USACH

mejoran bastante si en vez del factor 72, se introduce el factor empírico 150, de tal modo que: (8.7) Esta ecuación se conoce como ecuación de Blake-Koseny (ó ecuación de Koseny-Carman) y es válida para: (8.8) con

Go = ρ vo

Introduciendo en la ecuación de Koseny-Carman la definición del factor de fricción, se obtiene: (8.9) Para flujo altamente turbulento en columnas de relleno (Rep > 1000), con un desarrollo análogo al recién planteado, pero considerando un factor de fricción independiente del flujo, se obtiene: (8.10) Estas ecuaciones se conocen con el nombre de ecuaciones de BurkePlummer.

246

Capítulo 8: Factores de fricción en medios porosos

Sumando las expresiones de Koseny-Carman (válida en régimen laminar), y la de Burke-Plummer (válida en régimen turbulento), se obtiene la siguiente ecuación general, deducida por ERGUN: (8.11)

En esta ecuación Dp = dp ∙ φs, siendo dp el diámetro de la esfera equivalente de la partícula.

Figura 8.3. Representación gráfica del factor de fricción vs número de Reynolds (Adaptado de Levenspiel, Flujo de Fluidos e Intercambio de Calor, 1993).

247

Escurrimiento de fluidos

USACH

Figura 8.4. Relaciones gráficas para estimar la porosidad de columnas de relleno (Adaptadas de Levenspiel, 1993).

8.2.

Teoría del medio continuo

Otro enfoque teórico para el análisis del flujo en medios porosos es la aplicación de la mecánica del continuo a mezclas fluido-partícula. En este método se aplica la ecuación de continuidad y del movimiento a una matriz porosa indeformable: (8.12) (8.13) donde m es la fuerza resistiva, para la cual existen varios modelos. τ es el tensor tensión extra, que depende de la velocidad superficial. Existe muy poca información para este parámetro, aunque afortunadamente es un término de menor importancia frente a los otros.

248

Capítulo 8: Factores de fricción en medios porosos

Henry Darcy (1856) propuso para la fuerza resistiva, en cualquier régimen de flujo, la siguiente expresión: (8.14) k y c son parámetros que dependen solamente de la estructura de la matriz porosa, excepto cuando se producen interacciones físico-químicas entre la matriz porosa y el fluido. k representa la permeabilidad del medio poroso, con dimensión L2, mientras que c es adimensional.

El flujo lento (escurrimiento Darcyniano), se observa cuando . En esta situación, m = µ vo/k, con lo cual se obtiene la cono-

cida “Ley de Darcy”: -grad p = µ vo/k + ρ g, o bien, utilizando el concepto de altura piezométrica (P = p + ρg h): -grad P = µ vo/k (8.15) Los parámetros k y c son determinados experimentalmente en un permeametro:

Si el escurrimiento es incompresible, para cualquier régimen de flujo: (8.16)

249

Escurrimiento de fluidos

USACH

Si se trata de un gas ideal: (8.17) Relacionando el modelo capilar de ERGUN con el modelo recién expuesto se observa que: (8.18) Con

el diámetro medio de Sauter y β con valores entre 4 y 5.

Estas expresiones reproducen los valores experimentales de lechos con porosidades alrededor de 0.4 y valores de 10-5 < k < 10-4 cm2. Churchill y Usagi (1972) proponen una correlación para estimar c, cuyo rango de validez es: 0.15< ε Remf 0.001

> 0.3 - 0.5

Filtros de papel de alta eficiencia

> 0.5 – 1.0

Filtros de manga

Precipitadores electrostáticos

> 0.3 – 1.0

Venturi scrubbers

Torres spray

Ciclones

Cámaras de sedimentación

Colector

0.5 2.0

0.25 -- 0.5

1 - 10

15 –30

500 50

Sobre el 99.9% de todas las partículas bajo 5 µm

200-250

95-99 % de partículas bajo 5 µm 80% a más de 99% de todas las partículas

200-250

200-250

No sensible Poco efecto

No sensible Poco efecto

Seco

Muy crítico

No afecta la eficiencia

Seco o No afecta la Muy crítico eficiencia húmedo Poco efecto Seco o Crítico en la eficienhúmedo cia

Húmedo

Húmedo

Máxima Condición Sensibilidad Efecto de la temperatura del gas de al punto de densidad de de trabajo salida rocío las partículas (°C) Aumenta la Seco o No muy e f i c i e n c i a 500 húmedo sensible con la densidad Aumenta la Seco o eficiencia 500 Crítico húmedo con la densidad

90-99% de partículas bajo 5 µm

Rango de Caída de Grado de limpieza tamaño de presión, cm partículas, esperado agua µm 50% de partículas bajo 50 µm y µm > 150 10 1 - 3 de 95% de partículas sobre 50 µm 98 % de partículas >3 2 - 7 sobre 5 µm y 50% de partículas bajo 3 µm

Tabla 3.4. Carta de selección de equipos (Chemical Engineering, Julio 1993).

Escurrimiento de fluidos USACH

Capítulo 13: Separaciones gas-partícula

13.1. Separadores gravitatorios Los sedimentadores gravitatorios o cámaras de sedimentación, han sido tradicionalmente utilizados para remover partículas desde corrientes gaseosas. Entre sus ventajas debe considerarse su sencillez, bajo costo inicial y de inversión, y sus bajas caídas de presión. El diseño básico de una cámara de sedimentación consiste de una cámara de expansión, en la cual la velocidad de la partícula se reduce a tal grado, que ella puede sedimentar bajo la acción de la gravedad. Industrialmente se utiliza para remover partículas relativamente grandes (>40 µm), con velocidad del gas inferior a 1 m/s. Aún así, la eficiencia de separación es inferior al 50%. Normalmente este equipo cumple el rol de pre-separador o separador primario, ya que para conseguir mejores eficiencias de separación, el tamaño de equipo requerido sería demasiado grande.

Figura 13.1. Separador gravitatorio.

Para determinar las dimensiones básicas de un sedimentador gravitatorio, es imprescindible conocer las características de las partículas, específicamente su velocidad terminal, vt. Para que se produzca la retención de la partícula, esta debe llegar a la superficie horizontal inferior (piso), en el tiempo que ésta permanece en la cámara de separación (tiempo de residencia, tr). Considerando flujo pistón: tr = L/V ; V = Q/BH ⇒ tr = LBH/Q, (13.1)

381

Escurrimiento de fluidos

USACH

donde V es la velocidad del fluido, Q el flujo volumétrico de gas, H el alto de la cámara, L su largo y B su ancho. El tiempo requerido para que la partícula sedimente, ts, en la distancia H, mientras desciende a su velocidad terminal vt, es: ts = H/vt (13.2) Para que la captura de la partícula efectivamente se produzca, se debe verificar que: ts ≤ tr. En el límite: ts = tr. Por lo tanto: H/vt = LBH/Q, o bien: vt = Q/LB. Utilizando expresiones adecuadas para las velocidades de sedimentación (zona de Stokes, zona de transición o zona de Newton), es posible obtener ecuaciones para el tamaño de partícula (dp), a partir del cual son retenidas en la cámara de sedimentación. Estas ecuaciones se resumen en la tabla 13.5. Tabla 13.5. Expresiones para determinar tamaños de partícula. Régimen de flujo Transición

Stokes {18µQ/(gρpBL)}

0.5

Q

0.88

ρ

0.254

µ

/{0.193(g ρp)

0.377

Newton (LB)

0.623

}

0.88

0.33 (ρ/g ρp)(Q/LB)2

Los diámetros de partícula calculados con las ecuaciones anteriores representan valores límites, dp mín = dp*, puesto que partículas con tamaños iguales o mayores llegarán a la superficie inferior y partículas con diámetros menores saldrán de la cámara de sedimentación. En las ecuaciones anteriores no aparece la altura de la cámara de sedimentación, H, y por lo tanto puede concluirse que la eficiencia de la captura es independiente de dicho valor. Sin embargo, existe un valor mínimo de H, usualmente establecido para asegurar que la velocidad del gas sea lo suficientemente baja y que no produzca el arrastre de las partículas ya depositadas.

382

Capítulo 13: Separaciones gas-partícula

Si en el diseño de una determinada cámara de sedimentación se utiliza dp mín, la eficiencia de recolección de partículas, η, para partículas mayores o iguales a dicho tamaño será 100%. Puede demostrarse que: η = vt B L/Q 100% = (vt/V)(L/H) = (H*/H)100% (13.3) donde H* es la mínima altura que asegura la captura de las partículas.

Si en la corriente gaseosa existen partículas con una distribución de

tamaño, existirá una curva de eficiencia de recolección de partículas para cada tamaño. Si una partícula de tamaño dp está a H*, entonces (H*/H) representa la fracción de partículas de este tamaño que será recolectado. Si H* es igual o mayor que H, todas las partículas de ese tamaño o mayores serán recolectadas en la cámara de sedimentación. Para evitar el re-arrastre de las partículas depositadas, la velocidad del gas no debe ser superior a la velocidad de arrastre, vp. Una expresión semi-empírica para estimar vp es: (13.4) donde (W/G) es la velocidad de los sólidos arrastrados (masa de sólidos/ masa de gas); Ps es la presión estática del gas.

Figura 13.2. Cámara de Howard.

383

Escurrimiento de fluidos

USACH

Para aumentar la eficiencia y disminuir el tamaño del equipo, se introducen planchas horizontales dentro del separador, constituyéndose la llamada “cámara de Howard”. De esta forma se consigue separar partículas de hasta 15 µm. La distancia entre las placas puede ser tan pequeña como una pulgada. Una de las desventajas de este diseño es la dificultad para retirar los sólidos que se depositan en la cámara. Por este motivo no se recomienda su uso para concentraciones de partículas superiores a 1 gr/pie3. En esta familia de separadores, dadas las bajas velocidades del gas (1 a 10 pie/s), cabe esperar bajas caídas de presión, menores a 0.2 pulgadas de agua. El conjunto de dos o más tubos verticales, con diámetros crecientes por los cuales pasa el gas, constituyen los llamados elutriadores. En estos, las partículas mayores sedimentan en el fondo del tubo, mientras que las partículas menores son arrastradas hacia el tope.

Figura 13.3. Elutriadores.

Existe otro tipo de cámara de sedimentación, en la cual, junto con el efecto de la gravedad se inducen fuertes cambios de dirección del gas, de forma tal que las partículas impactan en placas internas. Estos equipos se

384

Capítulo 13: Separaciones gas-partícula

conocen con el nombre de separadores de momentum. Su eficiencia puede superar el 70%, para partículas mayores a 50 µm.

Figura 13.4. Separadores de momentum.

13.2. Separadores centrífugos (ciclones) Los ciclones son ampliamente utilizados en diversos procesos industriales: minería, metalurgia, procesos farmacéuticos, operaciones de combustión, etc. En los ciclones, el gas ingresa tangencialmente en la parte superior, adquiriendo un movimiento centrífugo, siendo las partículas desplazadas hacia la periferia, continuando con su movimiento espiral hacia abajo debido a la fuerza de gravedad, mientras tanto el gas libre de partículas (o con partículas de menor tamaño), asciende por la parte central. Cuando el gas ingresa al ciclón, su velocidad se distribuye de tal manera que su componente tangencial (Vct) aumenta al reducirse el radio, mientras

385

Escurrimiento de fluidos

USACH

la velocidad radial (Vcr) se dirige hacia el centro en casi todas las partes del ciclón. La separación de las partículas es influenciada por el campo centrífugo, velocidad radial, tiempo de residencia y la turbulencia dentro del ciclón. El tiempo de residencia depende del número de orbitales que se forman en el interior del ciclón, el que a su vez depende de factores geométricos y del flujo de gas. La caída de presión en el ciclón es provocada por: - Expansión a la entrada del ciclón. - Pérdida de energía cinética en la trayectoria espiral del gas, dentro

del ciclón. - Fricción en las paredes internas y en el ducto de salida.

Figura 13.5. Dimensiones básicas de un ciclón Lapple (Adaptado de Perry, Manual del Ingeniero Químico).

386

Capítulo 13: Separaciones gas-partícula

Obedeciendo en general a una misma forma, es posible agrupar los ciclones en tres categorías: a) Ciclones convencionales (eficiencia media). Con eficiencias entre 80 a 95%, son capaces de manejar flujos con caídas de presión entre 5 y 13 cm col agua. b) Ciclones de alta eficiencia (95 a 99%). En general son de menor diámetro y mayor longitud. Sus rangos de caídas de presión típicos están entre 5 y 16 cm col de agua. c) Sistemas multi-ciclones. Están constituidos por un conjunto de ciclones de pequeño diámetro que funcionan en paralelo. Sistemas bien diseñados pueden presentar eficiencias tan altas como 90% para partículas en el rango de 5 a 10 µm.

La eficiencia de un ciclón es usualmente especificada en términos de

un tamaño de corte, el cual corresponde al tamaño de partícula recolectada con una eficiencia del 50%, es decir: dp50. Este tamaño de corte depende del gas, de las propiedades de las partículas, de las dimensiones del ciclón (normalmente el tamaño de un ciclón se representa en función del diámetro de la sección cilíndrica, Dc) y de las condiciones de operación. Puede calcularse con la expresión: (13.5) µ es la viscosidad del gas en lb/pie-s, Bc es el ancho de la entrada del gas, pie, Nt es el número efectivo de espiras que el gas desarrolla en el ciclón, vi es la velocidad del gas en la entrada, pie/s, ρs es la densidad de las partículas, lb/pie3 y ρ es la densidad del gas, lb/pie3. Lapple entrega las dimensiones básicas de un ciclón convencional, demostrando que para esta configuración Nt = 5.44. Partículas mayores a dp50 serán recolectadas con eficiencias superiores al 50%, mientras que las de 387

Escurrimiento de fluidos

USACH

menor tamaño serán recolectadas con una eficiencia menor. Para determinar la eficiencia de colección puede utilizarse la ecuación siguiente: (13.6)



En general, la eficiencia de recolección aumenta con el aumento de

la velocidad de entrada, encontrándose ésta en el rango de 20 a 70 pie/s. Velocidades mayores aumentan la eficiencia, aunque valores superiores a 100 pie/s no son recomendables, ya que se produce el arrastre de partículas separadas, bajando la eficiencia de colección.

Instalaciones en serie se justifican solo si:

• El polvo presenta una amplia distribución de tamaños, incluyendo partículas bajo 15 µm. En esta situación un primer gran ciclón permite separar las partículas gruesas, luego varios ciclones menores en paralelo permiten separa las partículas menores. • El polvo es fino, pero con partículas que forman flóculos.

Los ciclones operan de igual modo a vacío o a presión.

Ejemplo 13.1 Debido a su cargo de Inspector Controlador de la Polución Atmosférica, usted ha sido comisionado para evaluar el funcionamiento de un ciclón usado a la salida de un secador de áridos de una planta de asfaltos. La norma vigente prohíbe la emisión de más de 229 [mg/m3] de partículas menores de 10 [μm]. Un análisis de partículas a la entrada del ciclón indica que el 10% de la emisión másica particulada está entre 5 y 10 [μm], siendo despreciable la concentración de partículas menores de 5 [μm]. La carga total de entrada al ciclón es de 11450 [mg/m3]. El ciclón disponible es uno convencional, de

388

Capítulo 13: Separaciones gas-partícula

diámetro 0,61 [m], con una velocidad a la entrada de 15,2 [m/s]. La densidad del polvo es 2,75 [g/cm3]. El número de giros del gas se estima en 4,5 y opera a una temperatura ambiente de 21 [ºC]. ¿Daría usted el permiso de operación para el ciclón? Solución: Considerando que el ciclón es del tipo Lapple, se tiene: Total de partículas: 11450 [mg/m3] Partículas ≤10 [μm]: 11450 · 0.1 = 1145 Se pide: Tamaño promedio de partículas ≤10 [μm]: Para el cálculo del diámetro de corte se utiliza la ecuación:

donde

;luego:

El cálculo de la eficiencia, Ei, se realiza mediante la ecuación:

El ciclón no cumple con la Emin requerida, por lo tanto no se da el permiso.

389

Escurrimiento de fluidos

USACH

Ejemplo 13.2 Una dada muestra de una corriente de aire tiene la siguiente distribución en % en peso con cada tamaño medio de partículas: Tamaño (μm)

1

10

50

100

200

% en peso

10

20

40

20

10

Se usará un separador ciclónico convencional, basado en lo siguiente: ancho de entrada, 12 pulg; Nº efectivo de giros o vueltas, 5; velocidad del gas de entrada, 60 (pie/s); densidad de la partícula, 1,6 (g/cm3); densidad del gas, 0,074 (lb/ft3); viscosidad del gas, 0,045 (lb/pie·h). a) Use el método de Lapple para estimar el % del peso total que sería removido en el ciclón. b) ¿Qué flujo de aire contaminado podría tratarse en el ciclón? Solución: Primero se determinará dp50. Para esto se transformaran los datos a SI Bc = 0,3048 (m) vi = 18,288 (m/s) ρs = 1600 (kg/m3) ρ = 1,185 (kg/m3) μ = 1,9x10-5 (kg/ m*s)

390

Capítulo 13: Separaciones gas-partícula

A continuación se muestra un ejemplo para el primer rango de tamaño. Los cálculos siguientes son análogos.

El cálculo de Ei se puede realizar mediante la ecuación:

y por último: x · Ei = 0,1 · 1,74 = 0,174 dp (μm) 1 10 50 100 200

x

dp/dp50

Ei

x · Ei

0,1 0,2 0,4 0,2 0,1

0,133 1,328 6,640 13,280 26,560

1,74 63,82 97,78 99,44 99,86

0,17 12,76 39,11 19,89 9,99

del peso total será removido en el ciclón b) Q = vi · Bc · Hc

391

Escurrimiento de fluidos

USACH

EJERCICIOS 13.1. En una de las etapas del proceso de fabricación de cemento se utiliza un ciclón Lapple, de diámetro 8 [ft]. Medidas de la velocidad del aire a la entrada del ciclón señalan que ésta es de 60 [ft/s] y que el número efectivo de espiras dentro del ciclón, Nt, es 5,4. En las condiciones del proceso la viscosidad del gas es 0,02 [cp] y la densidad de las partículas es 181 [lb/ft3]. El análisis granulométrico de las partículas es: Dp [μm] 2 8 15 20 35 45 55 >55

% peso 5 15 20 28 15 10 5 2

Determine el diámetro de corte y la eficiencia global del ciclón. [Solución: Eglobal = 78,9%].

13.2. La distribución de tamaños de partículas emitidas por un incinerador es: Rango de tamaños (μm)

d p (μm) 1 3 7 12,5 17,5 25 > 30

0–2 2–4 4 – 10 10 – 15 15 – 20 20 – 30 30 - +

392

% en peso 17 5 8 3 3 5 59

Capítulo 13: Separaciones gas-partícula



a) ¿Cuál es la eficiencia global de un ciclón convencional si se usa un dp50 de 10 (μm)? b) Si la caída de presión máxima permisible es de 5 (pulg de agua), ¿cuáles serían las dimensiones del ciclón para tratar 1500 (m3/h)? Solución [Et = Σx · Ei = 70,59%, Dc = 47,1(cm)].

393

Capítulo 14: Transporte hidráulico de sólidos en tuberías

CAPÍTULO 14 TRANSPORTE HIDRÁULICO DE SÓLIDOS EN TUBERÍAS

14.1. Introducción Existen numerosas aplicaciones de transporte hidráulico de sólidos a través de tuberías, especialmente en el área de la minería. Cabe destacar el transporte de suspensiones de cobre en el norte de Chile, de mineral de fierro en Brasil y de carbón mineral en EE.UU., lugares donde los minerales recorren distancias superiores a cien kilómetros. Prácticamente en todas las operaciones conocidas el fluido de transporte es el agua. Desafortunadamente, el avance en las aplicaciones no significa que se tenga el completo conocimiento teórico del efecto de las variables involucradas en el transporte hidráulico. Para que el transporte de mezclas líquido-sólido a través de cañerías sea técnicamente factible, se deben cumplir las siguientes condiciones: 1) El sólido debe poder mezclarse y separarse fácilmente. 2) No debe existir el riesgo de taponamiento de la cañería. 3) El desgaste y ruptura que sufren las partículas durante el transporte, no debe tener efectos adversos para el proceso posterior de ellas. 4) La cantidad del fluido empleado para el transporte debe ser adecuada. Variables del sistema El flujo de mezclas sólido-líquido por cañerías depende de diversos parámetros, no conociéndose con exactitud la influencia de ellas. Estos parámetros se sintetizan en la Figura 14.1:

395

Escurrimiento de fluidos

USACH

Figura 14.1. Variables a considerar en el diseño de un sistema hidráulico de material particulado.

Aspectos Básicos: Al aplicar un balance de energía mecánica a una tubería horizontal, por la cual circula un fluido, debido a un gradiente de presión, se obtiene: (14.1) donde (–∆p) representa la caída de presión a través del sistema, ρ la densidad del fluido y Eˆ v , las pérdidas por fricción. Si en el mismo sistema, en vez de un fluido puro, circula una suspensión se obtiene: (14.2) En esta situación, las pérdidas por disipación de la mezcla, pueden expresarse en altura de suspensión al dividir la última expresión por g. Si se divide además por L y denominando Jm al término 2fv2/gD, entonces, las pérdidas

396

Capítulo 14: Transporte hidráulico de sólidos en tuberías

por disipación viscosa estarán expresadas en altura de suspensión por unidad de longitud de tubería. (14.3) En forma equivalente, si el fluido que circula es agua, llamaremos en este caso a las pérdidas de disipación viscosa, expresadas en altura de agua por unidad de tubería, como Jw. La densidad de la mezcla, ρm, si se comporta como un fluido homogéneo está dada por: (14.4) donde Vw y Vs corresponden a los volúmenes de agua y sólido, respectivamente y V representa el volumen de la suspensión. Llamando C a la concentración en volumen de los sólidos en la mezcla, se tiene, que C = Vs/V y por lo tanto, se puede escribir también: 1 – C = Vw/V. Introduciendo, estas expresiones en la ecuación que representa la densidad de la mezcla, se obtiene: ρm = ρw + C(ρs – ρw) (14.5) Newitt, relacionó las pérdidas por fricción de la mezcla en altura de agua, J, con las pérdidas por fricción del agua sin partículas, Jw, en términos de la función φ, característica de cada sistema: J = Jw + Jw C φ Entonces la función φ se define como:

(14.6)

(14.7)

397

Escurrimiento de fluidos

USACH

14.2. Modelos para predecir pérdidas de carga La modelación del transporte hidráulico de sólidos está siendo perfeccionada continuamente por el aporte de diversos investigadores. Los modelos tradicionales corresponden a los formulados por Durand y también el modelo de Newitt. Modelo de DURAND (Durand, 1953) Basándose en numerosos valores experimentales, Durand propone la siguiente clasificación para las mezclas: A. Mezclas homogéneas. Válido para sólidos con tamaño inferior a 20 ó 30 micrones. Las mezclas homogéneas presentan dos tipos de regímenes de flujo: laminar y turbulento. A.1 Régimen laminar. Las mezclas se comportan como fluidos no newtonianos, con formación de depósito en la sección inferior de la tubería. Este tipo de comportamiento debe ser evitado en las instalaciones. A.2 Régimen turbulento. En este régimen, las pérdidas por disipación viscosa son iguales a las del agua pura, por lo tanto, J = Jw, y: (14.8) donde J representa las pérdidas por fricción de la suspensión, por unidad de longitud de cañería, Jw las pérdidas por fricción del agua pura por unidad de longitud de cañería y (Ev/g)m las pérdidas por disipación viscosa de la mezcla, expresada como altura de mezcla.

398

Capítulo 14: Transporte hidráulico de sólidos en tuberías

B. Mezclas intermedias. Válido para sólidos con tamaño entre 25 y 50 micrones. Durand no presenta expresiones para este tipo de régimen. C. Mezclas heterogéneas. Están constituidas por sólidos que exceden los 50 micrones de diámetro. Es posible encontrar la siguiente subdivisión, válida en rigor para sólidos con una estrecha distribución de tamaños. Para mezclas con una amplia granulometría, la clasificación puede variar, dependiendo de la concentración de los sólidos. C.1. Transportadas en suspensión (diámetro de sólidos entre 50 micrones y 0.2 mm). C.2. Categoría de transición (diámetro de sólidos entre 0.2 y 2 mm). C.3. Transportadas por saltos (diámetro de sólidos superior a 2 mm).

Los movimientos del agua son insuficientes para mantener a los sóli-

dos y el agua homogeneizados y la gravedad tiende a sedimentar a las partículas. Durand encontró que existe una velocidad límite del agua, vL, bajo la cual se inicia el depósito de partículas en la tubería. Encontró que vL depende del diámetro de la cañería D, de la densidad de los sólidos, ρs de la densidad del líquido, ρw, de la concentración volumétrica de sólidos C y del diámetro de las partículas Dp. Entrega la siguiente relación para estimar vL: (14.9) donde FL es una variable que depende únicamente de la concentración en volumen (C) y del diámetro de los sólidos, tal como se muestra en la figura siguiente:

399

Escurrimiento de fluidos

USACH

Figura 14.2. FL vs Dp para diversas concentraciones en volumen (Adaptado de Durand, 1953).

Por lo tanto, dependiendo de la velocidad límite de depósito, se desprende que en mezclas heterogéneas existen dos tipos de regímenes: régimen sin depósito de sólidos (v > vL) y régimen con depósito de sólidos (v < vL). Régimen sin depósito de sólidos en la tubería. Para el parámetro φ, informan de la siguiente expresión: (14.10) donde CD es el coeficiente de fricción representativo de las partículas. S es el cuociente ρs/ρw.

400

Capítulo 14: Transporte hidráulico de sólidos en tuberías

Régimen con depósito de sólidos en la tubería. Si la velocidad de operación es menor que la velocidad límite de depósito vL, en la parte inferior de la cañería comenzará a formarse un lecho, que con el tiempo puede obstruir totalmente la cañería. Sin embargo, experiencias visuales muestran que la cima de los depósitos es aproximadamente plana, y que por este motivo existe una altura h y una sección de flujo libre, siendo aplicable también la ecuación para φ, introduciendo ve en vez de v, siendo ve la velocidad de flujo sobre el depósito y cambiando D por 4 RH, donde RH es el radio hidráulico. En el diseño de sistemas de transporte hidráulicos, se recomienda emplear v > vL. Modelo de NEWITT (Newitt, 1955) Newitt clasifica las suspensiones en las siguientes categorías: a. Mezclas homogéneas. b. Mezclas heterogéneas en suspensión. c. Mezclas heterogéneas por saltos o lecho móvil. d. Mezclas heterogéneas con lecho estacionario. Esta clasificación puede visualizarse en la figura siguiente, en la que se grafica el diámetro de la partícula en función de la velocidad media, para una densidad de sólidos y un diámetro de cañería dados. Newitt propone además, ecuaciones que permiten calcular la transición entre cada uno de los patrones de flujo indicados en la Figura 14.3.

401

Escurrimiento de fluidos

USACH

Figura 14.3. Diámetro de partícula vs velocidad media.

Para mezclas homogéneas, se propone la siguiente correlación: φ = 0.6 (S – 1)

(14.11)

Para suspensiones heterogéneas se propone: (14.12)

402

Capítulo 14: Transporte hidráulico de sólidos en tuberías

Para suspensiones con lecho móvil: (14.13) Velocidades de transición

El cálculo de las velocidades de transición permite decidir que mo-

delo debe aplicarse a una determinada suspensión. Es decir, dichas velocidades permiten clasificar las mezclas. Las velocidades de transición vH y vB, se obtienen igualando las correlaciones correspondientes a dos zonas sucesivas y reemplazando v por vH o vB, según corresponda. De esta forma se obtiene para la transición entre mezcla homogénea y heterogénea: vH = (1800 g D vt)1/3 (14.14) Análogamente, para la transición entre mezcla heterogénea y flujo con lecho móvil, se tiene: vB = 17 vt (14.15) El criterio que permite decidir establece que si la velocidad media es mayor que la velocidad de transición, debe considerarse la mezcla en la categoría inferior, entendiéndose como orden creciente: homogéneas, heterogéneas, lecho móvil. Procedimiento de diseño Para diseñar un sistema de transporte hidráulico de sólidos por tuberías bajo presión, se recomienda que la mezcla fluido partícula sea una suspensión heterogénea; además de un completo conocimiento de las características de la mezcla y del sistema (topografía). Los modelos para evaluar las pérdidas de disipación viscosa están formulados para mezclas con estrechas granulometrías de tamaño, situación que en las aplicaciones industriales no

403

Escurrimiento de fluidos

USACH

se presenta. Con el propósito de incorporar en las ecuaciones el efecto de la distribución de tamaño, se recomienda utilizar un tamaño ponderado vía un coeficiente de arrastre representativo de la mezcla, el que se determina con la expresión siguiente: (14.16) donde Xi corresponde a la fracción en peso de las partículas de tamaño Dpi y CDi su correspondiente coeficiente de arrastre. Luego, con el valor del coeficiente de arrastre representativo de la mezcla, CD, usando las ecuaciones que representan este coeficiente, se obtiene el tamaño representativo de la mezcla, Dp y la velocidad de arrastre (llamada también velocidad terminal) de la mezcla, vt. Consideraciones particulares: i) Amplia granulometría y C ≤ 15% Para esta situación se recomienda el siguiente procedimiento de diseño: 1. Calcular el diámetro promedio de la mezcla de partículas, utilizando los coeficientes de arrastre. 2. Clasificar la suspensión según Durand. 3. Seleccionar un diámetro de cañería. 4. Calcular las pérdidas por disipación viscosa de la mezcla, Ev, usando preferentemente el modelo de Durand. 5. Si corresponde, calcular los requerimientos de potencia de bombeo. La selección del diámetro de la cañería, está directamente relacionada con la velocidad de la suspensión, considerándose que un diámetro inadecuado puede producir depósitos estacionarios de sólidos en la tubería.

404

Capítulo 14: Transporte hidráulico de sólidos en tuberías

La velocidad de operación óptima, vópt, es aquella que produce las pérdidas por disipación mínimas y puede obtenerse derivando la ecuación que relaciona velocidad con diámetro de tubería. Al usar la expresión correspondiente a Durand, se obtiene una expresión para la velocidad óptima, a la cual le corresponde un diámetro óptimo de cañería: (14.17)

(14.18)

Calculada la velocidad y el diámetro óptimos, se procede a calcular la velocidad límite de depósito, vL. Si es necesario, debe rediseñarse la velocidad de operación a fin de evitar la formación de depósitos, entrando en un ciclo iterativo que incluye el diámetro de la cañería. Una vez conseguidas las dimensiones con un régimen sin formación de depósitos, se procede a calcular las pérdidas por disipación viscosa. ii) Amplia granulometría y C > 15% Cuando se presenta este caso, no es posible calcular vL, ya que la ecuación planteada se aplica a concentraciones bajo 15%. En esta situación, las etapas básicas del procedimiento de diseño son las siguientes: 1. Calcular el diámetro promedio de la mezcla, usando los coeficientes de arrastre. 2. Seleccionar un diámetro mínimo de cañería, D, verificando que se cumpla la relación: vB < v > v,v se forma depósito.  Se tiene que disminuir el diámetro de la tubería. D = 8”  Di = 7,981”  Di = 0,203 (m) v = 1,72 (m/s) vL = 1,71 (m/s) ∴vvLL v,v si bien se cumple la condición para que no haya deposito, está muy al límite por lo que se opta por un diámetro aun menor. D = 6”  Di = 6,065”  Di = 0,154 (m) v = 2,98 (m/s) vL = 1,49 (m/s) ∴vvLL vv. Este será el diámetro de tubería a usar. Resolviendo el problema por el método Durand queda:

409

Escurrimiento de fluidos

USACH

Utilizando la expresión de Shacham para determinar el factor de fricción: fW = 0,0041

Realizando un Bernoulli entre los puntos (1) y (2)

entonces: , por lo tanto, la potencia es: Potenciateórica= Potenciareal =

410

Capítulo 14: Transporte hidráulico de sólidos en tuberías

EJERCICIOS 14.1. Determine el flujo de agua y la potencia de bombeamiento requeridas para el transporte hidráulico de 40 ton/h de arena. La tubería es de acero comercial de diámetro nominal de 5", tipo 40. El sistema debe operar con una velocidad de mezcla 20% mayor que la velocidad de depósito. La pérdida de carga total de los accesorios puede ser estimada en 25% de la proporcionada por las tuberías. Temperatura del sistema: 20 ºC. Densidad y esfericidad de la arena: 2600 kg/m3 y 0,78 respectivamente. La distribución granulométrica de la arena se muestra en la tabla siguiente: # Tyler -35+48 -48+65 -65+100



Fracción retenida 0,3 0,4 0,3

Solución: Flujo = 128 m3/h; Potencia teórica = 6.5 KW.

14.2. En una planta de fabricación de cemento se desea implementar un sistema para transportar hidráulicamente 140 m3/h de una suspensión de concentración 3.16% en peso y de 1.22% en volumen. En la succión se dispone de una presión de 2.76 · 105 Pascales (manométrica), estando la descarga a presión atmosférica, 6 metros bajo la succión. En la tabla siguiente se muestra la distribución granulométrica de las partículas, de densidad, ρs = 2600 kg/m3. 411

Escurrimiento de fluidos

USACH

Análisis por tamizado Diámetro (mm) 0,037 0,053 0,074 0,105 0,149 0,21 0,297 0,420 0,595

% peso (pi) 11,50 10,50 12,60 20,70 20,45 19,95 3,80 0,35 0,15



Utilizando los métodos de Newitt y Durand, determine las pérdidas por fricción, considerando que se utilizará una tubería de 220 metros de longitud y diámetro nominal 4".



Solución: Por ambos métodos se obtiene:

412

Eˆ v

g

≈ 3 ,5 (m

mezcla) .

Apéndice

APÉNDICE

413

Apéndice

Factores de conversión

415

Escurrimiento de fluidos

USACH

416

Apéndice

Temperatura t(ºF) = 1.8 t (ºC) + 32 t(ºC) = [ t (ºF) – 32]/1.8 t(ºK) = t (ºC) + 273,15 T(ºR) = t (ºC) + 459,67 T(ºR) = 1.8 T (ºK) T(ºR) =1.8 t (ºC) + 459,67 T(ºK) = t (ºF)/1.8 + 255,4 t(ºF) =1.8 T (ºK) – 459,67 t(ºC) = T (ºR)/1.8 – 273.15 Constante Universal de los Gases (R) R = 1,987 [cal/gramo mol ºK] R = 10,731 [psi pie3 /lbmol ºR] R = 0,08206 [lt at/g mol ºK] R = 0,7302 [at pie3 /lbmol ºR] R = 82,06 [cm3 at/g mol ºK] R = 62,36 [lt mm Hg/g mol ºK] R = 8,309 [Joule /gmol ºK] R = 8308,7 [m2 kg/kgmol seg2 ºK]

417

Escurrimiento de fluidos

USACH

TABLAS A.1. Dimensiones de tubos de acero normal (ASA Standards). Tamaño nominal (plg)

Diámetro exterior (plg)

1/8

0.405

1/4

0.540

3/8

0.675

1/2

0.840

3/4

1.050

1

1.315

1 4 1 1 2

1

1.660 1.900

2

2.375

1 2

2.875

3

3.500

1 2

4.000

4

4.500

5

5.563

6

6.625

8

8.625

10

10.75

12

12.75

2

3

Cédula Nº 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80 40 80

Espesor de pared (plg) 0.068 0.095 0.088 0.119 0.091 0.126 0.109 0.147 0.113 0.154 0.133 0.179 0.140 0.191 0.145 0.200 0.154 0.218 0.203 0.276 0.216 0.300 0.226 0.318 0.237 0.337 0.258 0.375 0.280 0.432 0.322 0.50 0.365 0.593 0.406 0.687

418

Diámetro interno (plg) 0.269 0.215 0.364 0.302 0.493 0.423 0.622 0.546 0.824 0.742 1.049 0.957 1.380 1.278 1.610 1.500 2.067 1.939 2.469 2.323 3.068 2.900 3.548 3.364 4.026 3.826 5.047 4.813 6.065 5.761 7.981 7.625 10.020 9.564 11.938 11.376

Diámetro interno (m) 0.00683 0.00546 0.00925 0.00767 0.01252 0.10744 0.01580 0.01387 0.02030 0.01885 0.02665 0.02431 0.03505 0.03246 0.04089 0.03810 0.05250 0.04925 0.06271 0.05900 0.07793 0.07366 0.09012 0.08545 0.10226 0.09718 0.12819 0.12225 0.15405 0.14633 0.20272 0.19368 0.25451 0.24293 0.30323 0.28895

Apéndice

A.2. Abertura de tamices Tyler. Nº de malla 3 4 5 6 7 8 9 10 12 14 16 20 24 28 32 35 42 48 60 65 80 100 115 150 170 200 230 270 325 400

Abertura (mm) 6.680 4.699 3.962 3.327 2.794 2.362 1.981 1.651 1.397 1.168 0.991 0.833 0.701 0.589 0.495 0.417 0.351 0.295 0.248 0.208 0.175 0.147 0.124 0.104 0.088 0.074 0.061 0.053 0.043 0.038

419

Escurrimiento de fluidos

USACH

A.3. Largo equivalente (L/D) en accesorios [Adaptado de Foust y col., 1961].

Accesorio

L/D

Convencional, 100% abierta Con disco de chaveta, 100% abierta Convencional, 100% abierta Válvula angular Con disco de chaveta, 100% abierta 100 % abierta Válvula de 75 % abierta 50% abierta compuerta 25% abierta 100 % abierta 75 % abierta Válvulas de lodo 50% abierta 25% abierta Giro convencional, 100% abierta Válvulas de retención Giro de despeje, 100% abierta Válvulas de aspira- Con disco de alza vertical, 100% abierta Con disco articulado, 100% abierta ción con cedazo

340 450 145 200 13 35 160 900 17 50 260 1200 135 50 420 75

Válvula de globo

Válvula de tres pasadas Codos T normal

Flujo directo Flujo a través de la bifurcación De 90º, radio normal De 90º, radio largo De 45º, radio normal Con flujo a lo largo Con flujo a través de la ramificación

420

44 140 30 20 16 20 50

Apéndice

A.4. Aspereza relativa en función del diámetro de la tubería [Adaptado de Foust y col., 1961].

421

Escurrimiento de fluidos

USACH

A.5. Factor de fricción en función del número de Karman [Adaptado de Foust y col., 1961].

422

Apéndice

A.6. Coeficiente de resistencia (K) para expansiones y contracciones. [Adaptado de Foust y col., 1961].

423

Escurrimiento de fluidos

USACH

A.7. Densidad y valores x, y para leer viscosidades de líquidos [Adaptado de Ocon y Tojo 1979].

424

Apéndice

425

Escurrimiento de fluidos

USACH

A.8. Viscosidad de gases y vapores a 1 atmósfera [Adaptado de Ocon y Tojo 1979].

426

Apéndice

427

Escurrimiento de fluidos

USACH

A.9. Presión de vapor y densidad de agua a diferentes temperaturas.

428

Apéndice

A.10. Curvas características de bombas centrífugas.

429