INTRODUCCIÓN
La velocidad de gases es una magnitud importante en diversos sectores de aplicación, incluyendo el monitoreo y control de emisiones atmosféricas, la aerodinámica, los sistemas de refrigeración, la meteorología, la construcción, la evaluación del potencial eólico y distintos procesos industriales, incluyendo los biomédicos, químicos, farmacéuticos y de semiconductores [1],[2]. Existe una gran variedad de instrumentos para medir la velocidad de gases, que difieren en tamaño, geometría, nivel de precisión, intervalo de medida, condiciones de medición y principio de operación. Entre los tipos de instrumentos de medición de velocidad de aire más utilizados se encuentran los tubos Pitot de tipo estático, los anemómetros ultrasónicos, los anemómetros térmicos o de hilo caliente, los anemómetros de copas y los anemómetros de hélices [3]-[5].
Es importante garantizar la fiabilidad y precisión de las mediciones de velocidad de gas. Para ello, los equipos de medición deben someterse a diversas actividades de mantenimiento y aseguramiento metrológico, y entre estas últimas se incluye la evaluación periódica del desempeño de los anemómetros por medio de calibraciones [6]. El problema de la calibración de anemómetros ha cobrado mayor interés en la medida en que las industrias y las agencias ambientales han estado definiendo requisitos para la obtención de mediciones más precisas de velocidad de fluidos, tanto en ambientes artificiales como en naturales [7], y también como consecuencia de los avances tecnológicos en los instrumentos de medición de velocidad de aire [8]. Los métodos de calibración pueden clasificarse en dos grandes grupos: directos e indirectos [9]. El método de comparación directa es el proceso más utilizado para la calibración de velocidad de gases, donde el Master Meter (MM - instrumento de referencia o patrón) y el Meter Under Test (MUT - instrumento bajo prueba o que se está calibrando) se exponen a un mismo flujo de aire, que posea una velocidad estable y homogénea sobre el área efectiva donde se llevará a cabo la medición [1]. Una manera de lograr la estabilidad y homogeneidad requerida por la calibración por comparación directa, consiste en la utilización de un túnel de viento especializado que posea un perfil de velocidad uniforme y una limitada turbulencia, con sistemas de control muy precisos para la generación del flujo [2],[10].
Un túnel de viento es una instalación diseñada para producir una corriente estable y controlada de aire, con dos posibles propósitos: (1.) determinar experimentalmente el efecto de flujo sobre modelos de edificaciones, paisajes urbanos, autos, aviones y otros objetos, y la dispersión de sustancias, entre otras aplicaciones de modelado, y (2.) la calibración de medidores de velocidad de aire [11],[12]. Los túneles de viento para aplicaciones metrológicas (de calibración, en este caso) requieren una muy buena caracterización del perfil de velocidad en la sección de prueba. Por ejemplo, algunos túneles incluyen toberas críticas corriente abajo de la sección de prueba para lograr una alta estabilidad en el flujo, y la caracterización se basa en el conocimiento de la relación entre la velocidad promedio en la sección de prueba y una velocidad puntual en la misma región [13].
Actualmente, es posible medir velocidades en túneles de viento con incertidumbres entre 0,1 % y 10 %, dependiendo del patrón de referencia y del intervalo de medición [14]. Como patrones de calibración y de verificación, los más usados por los laboratorios son los anemómetros láser Doppler (LDA), los tubos de Pitot y los anemómetros térmicos [3]. Un LDA mide la velocidad del aire por medio de la determinación de la velocidad de una partícula en la intersección de dos rayos láser; por esto, requiere de la adición de pequeñas partículas de agua en la corriente de aire para su medición [15]. El tubo de Pitot permite determinar la velocidad del gas por medio de la medición de la presión total de la corriente y de su presión estática, y usando la ecuación de Bernoulli [16]. El anemómetro térmico determina la velocidad del fluido midiendo la transferencia de calor desde un alambre calentado eléctricamente o un sensor de película, inmerso en el fluido [3]. Independientemente de los patrones, para garantizar la comparabilidad entre las mediciones realizadas en diferentes naciones, es importante determinar el grado de equivalencia entre los patrones nacionales con que los distintos países realizan las unidades de las varias magnitudes físicas medibles [17].
Comúnmente, los procesos industriales que involucran el uso de anemómetros para medir la velocidad de distintos gases se usan para prever caudales que son determinantes en la eficiencia y seguridad de dichos procesos, por lo cual la industria requiere mediciones fiables y precisas. Esto implica que se requiere una mejora continua de los métodos de calibración en túneles de viento, los cuales son materia de estudio en distintos institutos y laboratorios alrededor del mundo [18].
A pesar de que el túnel de viento se construya y opere adecuadamente, existen factores que alteran el perfil de velocidad durante su operación, de modo que es necesario determinar el efecto de estos factores sobre el resultado del proceso de calibración. Entre los estudios aplicados a túneles de viento metrológicos se encuentran: la evaluación de la capa límite (una región cercana a la superficie en la cual el comportamiento del perfil de velocidad es diferente, con un gradiente de velocidad mayor) en el túnel; la alteración de la velocidad promedio de flujo como función de la posición en la dirección del flujo de la corriente; la disminución de la incertidumbre de calibración utilizando el túnel; la evaluación del efecto de modificaciones geométricas en el túnel; y la caracterización de perfiles de velocidad en diferentes posiciones del túnel [8],[10],[19],[20]. En general, la calibración de medidores de velocidad en túneles de viento está limitada por la ineficiencia de la ecuación de Bernoulli a bajas velocidades, lo que ocasiona un incremento notable de la incertidumbre en los rangos bajos de velocidad [2].
La anemometría en túneles de viento se enfrenta a dos desafíos principalmente: primero, es muy poco probable que la velocidad en el punto de medición sea igual a la velocidad promedio en el ducto, y segundo, el hecho de que el anemómetro bajo prueba (MUT) afecte la forma del perfil de velocidad, hace necesaria la cuantificación y posterior aplicación de una corrección por bloqueo, debido a la obstrucción [17]. El efecto de bloqueo tiene, además, un efecto significativo sobre la incertidumbre del proceso de calibración [2]. Entre las alternativas para cuantificar la corrección por bloqueo se encuentran el uso de CFD, la determinación del perfil de velocidad por medio de mediciones de velocidad en la sección circundante, y la determinación de la permeabilidad del anemómetro [17]. En los túneles de viento para metrología más avanzados, se han reportado correcciones por bloqueo del orden de 0,9 % para anemómetros de paletas, y del orden de 0,1 % para anemómetros ultrasónicos [17]. Otros túneles para calibración de menor nivel metrológico reportan relaciones de bloqueo de alrededor de 3 % para los anemómetros de mayor tamaño [2].
Los estudios sobre el impacto en la calibración de instrumentos y las distorsiones que se puedan presentar en los perfiles de velocidad generados por el tamaño de los medidores, deben ser ejecutados específicamente para cada túnel de viento en sus condiciones de operación [21]. Con base en lo expuesto, el presente trabajo de investigación se enfocó en el estudio y modelado del perfil de velocidad del aire al interior del túnel de viento WSL20, desarrollado por el Centro de Desarrollo Tecnológico del Gas (CDT de Gas), tanto en ausencia como en presencia de diferentes anemómetros en el papel de MUT. La comparación se realizó primero por medio de una simulación, y posteriormente de manera experimental; en ambos casos, el análisis se realizó considerando 4 diferentes velocidades de aire en el túnel. La experimentación se hizo midiendo la velocidad del gas en diferentes posiciones radiales de la sección de prueba (donde se ubica el MM), comparando el caso en que no se introduce un MUT en el túnel de viento contra la situación en que está presente el MUT.
1. MATERIALES Y METODOS
1.1 Materiales
El WSL20 es un túnel de viento para aplicaciones metrológicas, el primero de su tipo en el país y el único desarrollado en Colombia. Posee capacidad para generar flujos de aire entre 0 m/s y 20 m/s, aunque el rango de calibración para velocidad de aire en el cual se acreditó el servicio de acuerdo a la norma ISO 17025 está entre 2 m/s y 20 m/s. El flujo en el túnel se genera por medio de 7 ventiladores axiales, y se controla por medio de un potenciómetro lineal. La figura 1 presenta un esquema del túnel WSL20, y la figura 2 muestra una foto del mismo. Mayor información sobre el túnel puede encontrarse en [1].
El MM usado en el túnel WSL20 es un anemómetro térmico de temperatura constante (CTA). Se seleccionaron instrumentos de diferentes dimensiones como MUT, para evaluar la distorsión debida a la obstrucción del área transversal en el túnel, al incluirse en un proceso de calibración, y determinar si es posible encontrar una relación entre el área obstruida y la distorsión ocasionada. Los MUT usados fueron un tubo de Pitot tipo L, un anemómetro rotativo Kestrel Pocket Weather Meter 4000 (Kestrel), y un termohigroanemómetro (THA). La figura 3 presenta el MM usado y los 3 MUT utilizados en el estudio.
1.2 Métodos
Inicialmente, se realizó una simulación en software CFD (Computational Fluid Dynamics) para visualizar la distorsión en los perfiles de velocidad ocasionada por la introducción de los diferentes MUT usados en el estudio. El diseño CAD (Computer-Aided Design) 3D de cada uno de los instrumentos de medición y de la sección de pruebas del túnel se realizó en el software SOLIDWORKS® 2015; los diseños fueron elaborados con un alto grado de similitud y detalle. Para la simulación dinámica del aire, se utilizó el software CFD ANSYS® 15.0. Se realizó análisis en tres dimensiones, tomando como dominio computacional un volumen de aire isotérmico a 25 °C en la sección del túnel, y efectuándolo en estado estable con modelo de turbulencia simple (k-épsilon). Se impusieron condiciones de frontera simples (entrada a velocidad constante, flujo desarrollado y salida a la presión atmosférica local). Se creó una malla tetraédrica, generada por ANSYS Mesh®, con tamaño automático y enfoque lagrangiano. La malla resultante contó con 24618 nodos y 127687 elementos. Los resultados se muestran en la sección 2.1.
Para caracterizar el perfil de velocidad resultante en cada condición experimental, se utilizó el CTA para realizar mediciones en 10 posiciones diferentes a lo largo de un diámetro del área transversal del túnel. La experimentación se llevó a cabo con 4 velocidades nominales de aire en el túnel: 2 m/s, 8 m/s, 14 m/s y 18 m/s. En cada velocidad nominal, se determinó el perfil de velocidad resultante en la posición longitudinal en la cual se ubica el MM, en cuatro condiciones diferentes: (1.) ningún instrumento en la posición del MUT, (2.) el tubo de Pitot en la posición del MUT, (3) el Kestrel, y (4) el THA. De esta manera, se obtuvieron 160 datos experimentales, que configuran un diseño de experimentos de un factor (equipo en la posición del MUT) con dos bloques (velocidad nominal del fluido y posición en el perfil de velocidad), en el cual la variable de respuesta es la velocidad del aire medida por el CTA.
Para analizar los datos, inicialmente se calculó el efecto de bloqueo de cada instrumento sobre las velocidades puntuales (en cada posición) para las diferentes velocidades nominales experimentadas, calculando la desviación porcentual en cada posición de la velocidad respecto al caso base (ningún instrumento en el área de flujo), de acuerdo a (1):
Donde d i,j .. es la desviación porcentual puntual en la velocidad de fluido en la posición ; al introducir como MUT el instrumento j, v lbgse es la velocidad del fluido en la posición i; cuando no hay instrumento en la posición del MUT, y v ¡ j es la velocidad de fluido en la posición ; al introducir como MUT el instrumento j. Estos resultados se muestran en la sección 2.2.
El segundo análisis consistió en la aplicación formal de la estadística a los resultados del diseño de experimentos. Al factor estudiado se le denominó "Equipo en calibración", y sus 4 niveles fueron denotados por "ninguno", "Pitot", "Kestrel" y "THA". Al primer bloque se le llamó "Velocidad nominal", con valores posibles "2", "8", "14" y "18", en m/s. El segundo bloque se designó "Distancia al borde", con valores "0,07", "3,8", "7,5", "10,9", "14,3", "14,7", "18,1", "21,5", "25,2" y "28,93", en cm, medidos desde el mismo sitio del borde interno del túnel de viento en la sección donde se ubica el MUT. Se comenzó aplicando un Análisis de Varianza (ANOVA), con miras a determinar si se encontraba evidencia fuerte de que el factor "Equipo en calibración" tiene efecto sobre la velocidad del fluido en la sección de pruebas. A continuación, se llevó a cabo una comparación múltiple entre los niveles del factor estudiado, para concluir cuáles de los niveles ocasionan una velocidad de fluido diferente; esto se hizo aplicando el método de mínimas diferencias significativas (LSD) de Fisher. Estos análisis estadísticos se presentan en la sección 2.3.
Como tercer análisis, se determinó la corrección por bloqueo que se debe aplicar para aquellos instrumentos que se haya determinado que ocasionan una velocidad diferente a la del caso base (sin obstrucción). Dicha corrección se determina a partir de las desviaciones porcentuales di,j (para el MUT en cuestión, a las diferentes velocidades nominales) en la posición en que se ubica el MM durante las calibraciones (7,5 cm desde el borde interno del túnel). Se aplicó regresión lineal múltiple para obtener modelos de las correcciones como función de la velocidad lineal. Estos resultados se presentan en la sección 2.4.
2. RESULTADOS Y ANÁLISIS
2.1 Simulación en software CFD
La figura 4 muestra la malla resultante para uno de los casos de simulación, a modo de ejemplo. La malla presentada corresponde a la simulación del flujo en el túnel con el THA, en la posición de calibración, y sus características fueron descritas previamente en la sección 1.2.
En las figuras 5, 6 y 7 se observan los contornos de velocidad generados en la sección transversal del túnel de viento operando a diferentes velocidades nominales, mediante la simulación de la presencia de los tres medidores bajo prueba.
En la figura 5 se puede observar que no hay una variación considerable del contorno para bajas velocidades. Para velocidades superiores, es posible identificar algún nivel de distorsión en el perfil de velocidad, aunque el efecto parece ser muy bajo.
En la figura 6 se observa que la simulación indica que la presencia del medidor Kestrel afecta considerablemente el mapa de contorno del perfil de velocidad. Se presenta una zona de muy baja velocidad (región azul) alrededor del MUT, y una zona de mayor velocidad de aire (región roja y naranja rojizo) a una distancia un poco superior. Esto se debe a que el aire que debería fluir en el sector donde se encuentra el equipo choca contra éste y pierde su energía reduciendo drásticamente su velocidad, y a su vez incrementa la velocidad de las líneas de flujo cercanas, que deben modificar su trayectoria conduciendo a una mayor densidad. Un análisis similar es válido para el THA, cuya simulación se presenta en la figura 7.
2.2 Efecto de la presencia del MUT sobre la velocidad puntual en la sección de calibración
La figura 8 presenta los perfiles de velocidad obtenidos en el túnel de viento, con diámetro de 30 cm, a las velocidades de 2, 8, 14 y 18 m/s (representadas en la figura en los paneles a, b, c y d, respectivamente), comparando el perfil cuando no hay obstrucción en el túnel con el obtenido al introducir cada uno de los instrumentos estudiados (tubo Pitot, Kestrel 4000, THA).
Se puede observar que la presencia del tubo de Pitot tipo L, en el interior de la sección de pruebas del túnel de viento, no genera una variación considerable del perfil de velocidad medido por el CTA, ocasionando una distorsión que no supera el 1 % de la velocidad puntual. Por otra parte, la presencia del Kestrel genera una desviación de mayor magnitud en el perfil de velocidad del CTA, con incrementos puntuales en la velocidad que oscilan entre 2,5 % y 5,6 %. El caso del THA es similar al del Kestrel; su presencia ocasiona desviaciones puntuales en la velocidad de entre 3 % y 4,6 %.
2.3 Análisis estadísticos del diseño de experimentos
En la tabla 1 se muestra el resultado del ANOVA a los resultados del diseño con un factor y dos bloques. Se concluye que para un nivel de significancia de 0,05, la evidencia experimental apoya la hipótesis de que la velocidad del aire es afectada por la presencia de instrumentos en la posición del MUT.
Se procede a determinar cuáles de los instrumentos causan una velocidad diferente a la presentada en el caso base, que consiste en la situación en que no hay obstrucción en el túnel (es decir, no hay MUT). Para esto, se emplea el método LSD de Fisher, enfocándose en las comparaciones entre el nivel 1 del factor (ningún instrumento en la posición del MUT) y cada uno de los demás niveles; los resultados se presentan en la tabla 2. El análisis mediante el método LSD permite afirmar que no hay una diferencia significativa entre el valor medido por el CTA sin obstrucciones (nivel 1 del factor) y el valor medido por el CTA cuando se encuentra el tubo Pitot tipo L dentro del túnel (nivel 2), por lo cual se concluye que este medidor no genera un efecto de bloqueo significativo. Sin embargo, para los medidores Kestrel y THA hay evidencia fuerte de dicha diferencia, lo que indica que la velocidad en la posición donde se ubica el CTA durante una calibración es alterada significativamente por la presencia de dichos instrumentos.
2.4 Determinación de la corrección por bloqueo
Se procede a estudiar el sesgo en la velocidad de gas en la posición del instrumento patrón, para las condiciones que distorsionan la velocidad en dicha posición. La tabla 3 muestra los sesgos porcentuales en la velocidad di,j para cada una de las velocidades nominales en el túnel de viento.
Se aprecia en la tabla 3 que el error porcentual varía ostensiblemente con la velocidad nominal para el caso del Kestrel, mientras que para el THA se podría considerar que una misma corrección sería válida para las distintas velocidades estudiadas. Se decide modelar el error en la velocidad como función de la velocidad nominal en el túnel, por medio de un polinomio cuadrático para el Kestrel, y como una constante para el THA. La aplicación de regresión lineal múltiple a los datos, según los órdenes polinomiales recién mencionados, conduce a (2) y (3). De esta manera, las siguientes ecuaciones permiten determinar la corrección que debe aplicarse a la velocidad medida por el instrumento patrón en la calibración de anemómetros, aplicable a instrumentos Kestrel y THA de dimensiones similares a los usados en el estudio, en calibraciones con velocidades nominales del aire de entre 2 m/s y 18 m/s:
Vkestrel, corr: velocidad de aire, corregida por el efecto de bloqueo del Kestrel, del instrumento patrón, durante Kestrel.
VKestrel: velocidad de aire medida por el instrumento patrón durante la calibración de un Kestrel.
VTHAcorr: velocidad de aire, corregida por el efecto de bloqueo del THA, del instrumento patrón, durante la calibración de un THA.
vtha: velocidad de aire medida por el instrumento patrón, durante la calibración de un THA.
3. CONCLUSIONES
La simulación realizada en el software CFD predijo adecuadamente el comportamiento cualitativo de la velocidad de aire al interior del túnel de viento, cuando se introducen en éste diferentes tipos de instrumento de medición de velocidad de gases.
La experimentación realizada en el presente trabajo permite concluir que la existencia de efecto sobre la velocidad de aire medida por el instrumento patrón en la calibración de anemómetros en un túnel de viento depende de las dimensiones del instrumento bajo prueba. Los instrumentos de menores dimensiones (como un tubo de Pitot) no causan una distorsión significativa, mientras que instrumentos de mayor tamaño sí afectan significativamente la medición registrada por el instrumento patrón, al obstruir una porción importante del área transversal del túnel de viento. Las anteriores conclusiones están soportadas en un ANOVA y una aplicación del método LSD de Fisher a los resultados del diseño experimental. Cuando el efecto de bloqueo es significativo, es fundamental cuantificar el sesgo ocasionado en la velocidad medida por el instrumento patrón, de manera que sea posible corregir el resultado de la calibración considerando este efecto y no se presente un error sistemático en las calibraciones realizadas en el túnel de viento.
Para el caso del túnel de viento WSL20, que fue el objeto de estudio del presente trabajo, se obtuvieron ecuaciones de corrección aplicables a la velocidad de aire medida por el instrumento patrón cuando el instrumento bajo prueba es un Kestrel o un THA, a partir de la aplicación de regresión lineal múltiple a los sesgos en la velocidad identificados en la posición particular donde se encuentra el MM durante la calibración. La aplicación de estas ecuaciones de corrección impedirá que el efecto de bloqueo afecte los resultados de calibración de equipos de tamaño significativo, al ser calibrados en el túnel WSL20. La determinación de ecuaciones de corrección para la calibración de anemómetros en túneles de viento diferentes, requiere la ejecución de pruebas experimentales específicas usando el túnel de viento de interés.