Exponentes de escala de las funciones de estructura longitudinales y transversales en turbulencia homogénea e isotrópica generada en un túnel de viento
Scale exponents of
longitudinal and transverse structure functions in homogeneous and isotropic turbulence
generated in a wind tunnel
1Departamento de Ingeniería Mecánica, Facultad de Ingeniería, Universidad EAFIT, Medellín, Colombia. Correo electrónico: lvelas30@eafit.edu.co
2 Grupo Energía Alternativa, Departamento de Ingeniería Mecánica, Facultad de Ingeniería, Universidad de Antioquia, Medellín, Colombia. Correo electrónico: lisabel.velasquez@udea.edu.co
Usando una malla activa y una malla pasiva
se caracterizó un flujo turbulento homogéneo e isotrópico estudiando los
exponentes de las funciones de estructura longitudinal y transversales y , donde u es la fluctuación de la velocidad longitudinal, y v es la
fluctuación de la velocidad transversal, en función del número de
Reynolds para órdenes, n, entre 2 y 8. La turbulencia se generó en un
túnel de viento horizontal de 91 × 91 cm2 de sección transversal, y 9,14 m de
longitud. El campo de velocidades se midió con anemometría de hilos calientes
con sensores TSI 1241 y TSI 1210 fabricados con una aleación de platino y rodio
de 2,54 de diámetro. Los resultados indicaron que los exponentes de
escala varían en función de , cuanto más aumentaba
, mayor era la diferencia entre los exponentes de las funciones de
estructura y los exponentes predichos por la teoría de Kolmogorov
sobre turbulencia.
Palabras clave: anemometría de hilos calientes; funciones de estructura; malla activa; malla pasiva; túnel de viento; turbulencia homogénea; turbulencia isotrópica.
By using an active
grid and a passive grid, a homogenous and isotropic turbulent flow was
characterized by studying its longitudinal structure functions and transverse where
u is the longitudinal velocity fluctuation, and v is the
fluctuation of the transverse velocity, as a function of the Reynolds number
for orders, n, between 2 and 8. The
turbulence was generated in a wind tunnel with 91 × 91 cm2 in
cross-section and 9,14 m in length. The velocity field was measured with
hot-wire anemometry with TSI 1241 and TSI 1210 sensors manufactured with an
alloy of platinum and rhodium of 2,54 diameter.
The results indicated that the exponents of scale vary according to
; if increases,
the difference between the exponents of the structure functions and the
exponents predicted by the K41 theory also increases.
Keywords: active grid; homogenous turbulence; isotropic
turbulence; structure function; passive grid; hot-wire anemometry; wind tunnel.
La turbulencia generada en túneles de viento a través del paso de un flujo de aire por mallas ha sido extensamente estudiada por más de 50 años. Estos túneles son comúnmente usados para crear turbulencia homogénea e isotrópica, gracias a que es posible controlar las condiciones de operación. La turbulencia isotrópica es el tipo de turbulencia donde el valor de cualquier función de las componentes de la velocidad y sus derivadas en un punto, definido en relación a un conjunto particular de ejes, son inalterados si los ejes de referencia son rotados en cualquier dirección, o si el sistema coordenado es reflejado en cualquier plano a través del origen. La turbulencia homogénea describe un tipo de turbulencia cuyas propiedades estadísticas son independientes de la posición en el fluido. Por tanto, la turbulencia isotrópica y homogénea representa la turbulencia, que es la misma en todas las direcciones y en todas las posiciones.
Entre los diferentes análisis que pueden realizarse sobre este fenómeno, las propiedades estadísticas de la turbulencia homogénea e isotrópica a menudo se investigan mediante el estudio de las funciones de estructura. Las funciones de estructura de orden n para la velocidad longitudinal u se definen así:
Donde son velocidades a lo largo del eje x, eje paralelo al flujo, para dos puntos separados por una escala o distancia, r, y 〈-〉 indica valores promedio.
Para la turbulencia homogénea e isotrópica, Kolmogorov estableció en su teoría K41 [1] que la relación entre las funciones de estructura y r, que corresponde a una distancia dentro del rango inercial, podía representarse siguiendo una ley de potencias:
Donde es el exponente de escala de la función de estructura longitudinal ∆ de orden n. El rango inercial corresponde a longitudes donde la viscosidad y la tasa de difusividad son despreciables, 𝐿 >> 𝑟 >> 𝜂, donde es la escala disipativa, y 𝐿 es la escala integral (escala donde la energía es inyectada). Kolmogorov supuso además que, como la tasa de disipación, ϵ, permanecía constante en este rango, existía una dependencia lineal entre los exponentes de escala y el orden de la función de estructura, . Sin embargo, existe una diferencia significativa entre su predicción y los resultados experimentales. Estas diferencias o desviaciones respecto de esa forma lineal se conocen como intermitencia de la turbulencia a pequeña escala.
Es posible también, con la descripción de Kolmogorov, el cálculo de las funciones de estructura transversal, de la que existen dos tipos: la función de estructura para la velocidad transversal, v, medida en la dirección del flujo, x, que se calcula mediante la ecuación 〈∆v(x)n〉 =
n
〈(v(x + r) − v(x)) 〉, y habrá una serie de valores diferentes para sus exponentes de
escala, ; y la función de estructura para la velocidad longitudinal,
u, medida con separaciones en dirección transversal al flujo, y, que se
calcula mediante la ecuación 〈∆u(y)n〉 =
n
〈(u(y + r) − u(y)) 〉, y que también tendrá una serie de valores propios para los
exponentes de escala, .
En este reporte se documentan los resultados de una investigación experimental de flujos turbulentos. Estos experimentos estuvieron enfocados en la medición de los exponentes de escala longitudinal y transversal para turbulencia generada en un túnel de viento por una malla pasiva y una malla activa. Un aspecto importante de estos experimentos fue la posibilidad de medir simultáneamente las funciones de estructura 〈∆u(x)n〉 y 〈∆v(x)n〉 usando sensores bidimensionales.
La turbulencia fue generada en un túnel de viento horizontal, de succión y de circuito abierto de 91 ×91 cm2 de sección transversal y 9,14 m de longitud ubicado en la Facultad de Ingeniería Mecánica y Aeroespacial de la Universidad de Cornell, en Ithaca, Estados Unidos. Un túnel de viento de succión y de circuito abierto es aquel en que el aire se aspira directamente de la atmósfera, y después de hacerlo pasar por la sección de ensayo o prueba, se devuelve nuevamente a ella. El túnel fue diseñado inicialmente para el estudio de flujos estratificados que involucraban no solo perfiles de velocidad, sino también de temperatura [4]. Un esquema del túnel con los diferentes elementos que lo componen se proporciona en la figura 1.
El túnel fue fabricado con madera contrachapada de 0,19 m de espesor, y estaba cubierto con hojas aislantes de fibra de vidrio de 0,05 m de espesor. La sección de ensayo permitía el acceso óptico hacia el interior del túnel, a través de 34 agujeros ubicados en las paredes laterales y 18 agujeros en las paredes superior e inferior. Una de las paredes laterales del túnel era ligeramente divergente, con ángulo de paso de 0,3°, para asegurar una velocidad media constante en la línea central del túnel a pesar del desarrollo de las capas límite; la superficie interna de la sección de pruebas no tenía una rugosidad significativa, para que las variaciones de velocidad y presión dinámica del flujo no superaran el 25 % y 50 % de sus valores medios, respectivamente.
El flujo fue acelerado antes de entrar en la sección de ensayo en una contracción con una relación de áreas de 9:1 de 1,83 m longitud. En la unión entre la contracción y la sección de ensayo fue instalada la malla encargada de la generación del flujo turbulento. El final de la sección de ensayo se conectó a un difusor, en el que se hacía una transición suave de la sección cuadrada a una sección circular que albergaba un ventilador axial 1,24 m de diámetro. El ventilador era alimentado por un motor de 60 HP, que giraba a una velocidad máxima de 1070 rpm. El motor que impulsaba este ventilador estaba situado en el exterior del túnel de viento, aislado del resto del túnel para reducir la transmisión de vibraciones mecánicas. El motor del ventilador estaba conectado a un variador de frecuencia que permitía ajustar la velocidad de giro del motor con incrementos de 0,5 Hz, para generar diferentes velocidades dentro del túnel. La velocidad máxima resultante de la instalación fue de 20 m/s. Para las pruebas se seleccionaron velocidad 𝑈0 del aire dentro del túnel de 3, 6, 9 y 12 m/s, lo que generó variaciones para el número de Reynolds definido en función de la escala de Taylor λ de 100 hasta 1200.
Los generadores de turbulencia pasivos son el método más simple para producir turbulencia aguas abajo en túneles de viento. Una malla pasiva consiste en un arreglo de barras cilíndricas o planas que forman una sección cuadrada. La malla pasiva usada en las mediciones estaba compuesta por 8 barras planas horizontales y 8 barras planas verticales de aluminio, con un espesor de 19 mm unidas a un marco, figura 2a. Las barras están separadas por una distancia de 8,26 cm, horizontal y verticalmente, lo que produce una separación entre los puntos medios de cada barra, M, de 10,16 cm. Una malla activa es un generador de turbulencia que utiliza elementos activos para producir turbulencia. La malla activa estilo Makita [2] es un conjunto de 100 alas cuadradas acopladas a una malla de barras redondas. Las barras cilíndricas son accionadas por motores paso a paso (200 pasos por revolución) en una secuencia aleatoria para crear una distribución de obstrucción transitoria en la sección de entrada de prueba del túnel de viento. La malla activa está compuesta por 8 barras horizontales y 7 barras verticales de aluminio con un diámetro de 12,7 mm unidas a un marco a través de rodamientos de baja fricción, figura 2b.
Las barras están separadas horizontal y verticalmente por una distancia de 82,6 mm, lo que produce una separación M de 12,1 cm. Ambas mallas fueron diseñadas y fabricadas para experimentos previos realizados en el mismo túnel de viento [3][4] [5].
El campo de velocidades fue medido con anemometría de hilos calientes en 17 diferentes posiciones, o puntos de monitoreo, a lo largo de túnel usando anemómetros de temperatura constante Dantec 55M01. Un anemómetro de hilo caliente consiste en un filamento calentado que queda expuesto al paso de un flujo. Este filamento está conectado a un circuito eléctrico que es capaz de monitorear las variaciones de resistencia eléctrica por la acción del flujo. Entonces, se puede establecer una relación de pertinencia entre la velocidad del flujo y la resistencia observada en el filamento calentado a determinada temperatura.
Para la adquisición de datos fueron usados dos tipos de hilos calientes: sensor unidimensional TSI 1210 y sensor bidimensional TSI 1241. El sensor TSI 1241, figura 3, fue usado para medir las componentes axial y tangencial de la velocidad, componentes necesarios para el cálculo de las funciones de estructuras 〈∆u(x)n〉 y 〈∆v(x)n〉. Para la medición de 〈∆u(y)n〉 fueron usados dos sensores TSI 1210, uno de los cuales estaba ubicado en una posición fija dentro del túnel, mientras que el otro se movía de forma perpendicular al flujo, para que la función de estructura en dirección y pudiera ser determinada.
Todos los sensores fueron fabricados de manera manual, con hilos de aleación de platino y rodio (al 10 %) de 2,54 μm de diámetro. Los hilos fueron calibrados usando el enfoque de ángulo efectivo, antes y después de cada medición, para encontrar las constantes de la ley de King [7], que luego permitirían relacionar la velocidad del flujo con el voltaje necesario para mantener el hilo a la temperatura fijada dentro del flujo.
Para la adquisición de los datos se utilizó
una tarjeta de adquisición de datos NI USB-6221 de ocho canales; por cada
velocidad y cada localización medida dentro del túnel el tamaño de la muestra
fue de 1638400 datos; la adquisición de la señal fue controlada con el software
comercial LabView. Debido al elevado volumen de
datos experimentales recogidos, se procedió a su reducción mediante técnicas
estadísticas para ser analizada en el dominio del tiempo y de la frecuencia,
este proceso de reducción condujo a obtener en cada punto de medición un solo
valor característico por cada variable de interés; el posprocesamiento
de los datos fue llevado a cabo con el software comercial Matlab.
Hay varias características clave que deben buscarse para verificar que la turbulencia producida por una malla sea homogénea e isotrópica. Por lo general, la velocidad media es homogénea dentro de un umbral de ±1% [7]. En segundo lugar, la relación de isotropía, que es la relación de las raíces medias cuadráticas de la velocidad longitudinal y transversal, deberá ser la unidad, en un caso ideal, aunque normalmente toma valores entre 1,2 y 1,4 [8]. Además, la distribución estadística de la turbulencia debe ser aproximadamente Gaussiana, y la gráfica de la distribución de probabilidad debería colapsar para diferentes posiciones en el flujo [2] [7] [9] [10].
La homogeneidad e isotropía del campo de velocidad fue analizada tomando mediciones en dirección transversal al flujo en diferentes posiciones x para ambas mallas. La distancia x corresponde a la separación entre el sensor de medición y la malla; por tanto, una distancia 𝑥 = 0 corresponde a la ubicación de la malla. Durante la primera parte del estudio, se usó la malla pasiva y se midieron velocidades 𝑈0 de 3, 6 y 12 m/s, se tomaron mediciones para distancias x entre 482 mm y 6070 mm; para la segunda parte, se usó la malla activa, con
La figura 4 muestra las variaciones transversales de la velocidad media 〈𝑈〉 normalizada por la velocidad del viento 𝑈0, como función de la posición transversal y para diferentes posiciones x para malla pasiva y activa, para solo una velocidad de operación por cada malla. Las distancias x y y fueron normalizadas por el tamaño de malla M. La coordenada 𝑦 = 0 corresponde al centro del túnel. La sección de prueba para la malla pasiva comenzó en 𝑥/𝑀=4,75, mientras que para la malla activa la sección inició en 𝑥/𝑀=29,18.
La figura 4 indica que los perfiles transversales de la velocidad media normalizada en la malla pasiva son no homogéneos cerca a la malla, pero esta falta de homogeneidad disminuyó con el aumento de x, debido a la acción de la difusividad, consistente con los resultados presentados por [11]. Las variaciones más pronunciadas de 〈U〉/𝑈0, que implican mayor falta de homogeneidad, se generaron en la malla pasiva para una velocidad 𝑈0=6 m/s, en donde la relación 〈U〉/𝑈0 tomó valores de 1,32 en 𝑥/𝑀=4.75. Esta falta de homogeneidad cerca a la malla es causada por la formación de estelas al paso del aire, a través de las barras que conforman la malla.
Para las mediciones con malla activa no hubo problemas de homogeneidad en los perfiles transversales de la velocidad media, porque las mediciones se realizaron a mayores distancias de la malla. La relación 〈U〉/𝑈0 fue muy similar en las cinco posiciones 𝑥/𝑀 medidas en cada velocidad. Para la posición 𝑥/𝑀=45,97 esta relación tendió a 1,26, 1,08 y 1,01 para las velocidades 𝑈0 de 3, 6 y 9 m/s, respectivamente, lo que indicó que cuanto mayor fue 𝑈0 más cercano fue 〈U〉/𝑈0, a la unidad, que es el valor esperado para la turbulencia homogénea. Para la malla activa se puede asegurar que por los valores de 〈U〉/𝑈0 se está siempre ante un flujo homogéneo, mientras que para la malla pasiva se requiere que el valor de 𝑥/𝑀 sea mayor de 19,75.
Para evaluar la distancia desde la malla donde el flujo es isotrópico, se introduce el factor de isotropía, que es la relación de las raíces medias cuadráticas de las fluctuaciones de la velocidad longitudinal y transversal,
𝑢′/𝑣′, donde ; el factor de isotropía deberá ser 1,0 cuando el flujo es isotrópico. En experimentos previos con mallas activas se reportaron valores de 𝑢′/𝑣′cercanos a 1,21 [5] y para mallas pasivas valores de 1,1 o ligeramente menores [8]. Valente y Vassilicos [12] han obtenido en sus experimentos con mallas pasivas valores de anisotropía en el rango de 1,10 y 1,15.
La figura 5 muestra el factor de isotropía como función de la posición transversal 𝒚 para ambos tipos de malla. El factor de isotropía fue de aproximadamente 1,05 para la malla pasiva y de 1,11 para la malla activa, que son valores muy cercanos a los valores de referencia. Este factor disminuyó al aumentar 𝑈0 y fue constante en todo el rango de medición, lo que indica que es independiente de la posición 𝑥/𝑀.
Con el objetivo de perfeccionar la
estimación del inicio de la región homogénea e isotrópica, se usa la evolución
longitudinal de la asimetría y la curtosis. La
asimetría o skewness está definida
como 𝑆(𝑢) = 〈𝑢3 〉/(〈𝑢2 〉)3/2 y es una medida de la
falta de simetría estadística en el flujo; la curtosis
o flatness está definida como 𝐾(𝑢) = 〈𝑢4 〉/(〈𝑢2 〉)2
y es la medida de concentración de los datos alrededor del valor medio. En turbulencia
homogénea e isotrópica las tres componentes de las fluctuaciones de la
velocidad tienen una distribución de densidad de probabilidad muy cercana a una
función normal; una distribución de probabilidad normal o Gaussiana tendrá un
valor de asimetría de cero y una curtosis de 3.0 [13]
[14].
La asimetría puede tomar valores positivos y negativos, pero, de acuerdo con los experimentos observados en turbulencia, la asimetría debe ser negativa [13] [14], porque representa un proceso de transferencia de energía. La figura 6 muestra la evolución longitudinal de la asimetría y la curtosis de la velocidad longitudinal en el centro del túnel (𝑦/𝑀 = 0).
El flujo se vuelve aproximadamente
Gaussiano alrededor de 𝑥/𝑀 = 20 para la malla pasiva donde 𝑆(𝑢) ≈ 0,05 y 𝐾(𝑢) ≈ 2,95. Estas estadísticas toman valores más altos para la malla
activa, pero también se vuelven constantes alrededor de 𝑥/𝑀 = 30, donde 𝑆(𝑢) ≈ 0,3 y 𝐾(𝑢) ≈ 3,2.
Estos resultados para la asimetría y la curtosis están en buen acuerdo con los estudios anteriores de [9] y [7]. Las gráficas para las variaciones de 𝑆(𝑣) y 𝐾(𝑣) no son mostradas, pero su comportamiento fue consistente con los de 𝑆(𝑢) y 𝐾(𝑢).
La figura 7 muestra la distribución de
probabilidad (o PDF, por sus siglas en inglés Probability
Density Function) de la
velocidad longitudinal 𝑢, medida para una posición en el centro del túnel para las
seis configuraciones o velocidades estudiadas en 𝑥/𝑀=44,75 en la malla pasiva y 𝑥/𝑀=45,97 en la malla activa, normalizada por la raíz de la media
cuadrática, PDF, y han sido denotados
como PDF (z).
Las distribuciones de probabilidad son claramente cercanas a la distribución Gaussiana, debido a que los valores asimetría y curtosis son muy cercanos a cero y tres, respectivamente.
De acuerdo con los perfiles de velocidad media, del factor de isotropía, de la asimetría y la curtosis, el flujo fue considerado razonablemente homogéneo e isotrópico para posiciones 𝑥/𝑀 >19,75 en malla pasiva y 𝑥/ 𝑀>29,18 en malla activa. La tabla 1 muestra los principales parámetros del flujo para diferentes posiciones 𝑥/𝑀 y diferentes velocidades 𝑈0 para ambos tipos de malla.
Una vez que se comprobó que el flujo fuera homogéneo e isotrópico, se procedió con el cálculo de los exponentes de las funciones de estructura. Los exponentes de escala fueron calculados usando el método ESS (por sus siglas en inglés Extended Self Similarity), propuesto por Benzi et al. [15], que permite el cálculo de los exponentes graficando la función de estructura de tercer orden en la abscisa frente a la función de estructura de cualquier otro orden en la ordenada. Benzi et al. [15] formularon lo siguiente:
La convergencia estadística de la función de estructura de orden impar es significativamente más lenta que el de las órdenes pares. Para superar esto, han sustituido las funciones de estructura de la ecuación 3 por su módulo:
Y establecieron además lo siguiente:
La ecuación 5 también es válida para el cálculo de los exponentes transversales.
Si se representa como función de 𝑟/𝜂, existe una zona plana que indica la existencia de un exponente
constante. En la figura 8 se han representado los valores de
para tres diferentes números de Reynolds y para valores
de 𝑛 entre 2 y 8; se ha delimitado el inicio y el fin de esta zona plana
con líneas punteadas verticales; las líneas horizontales punteadas corresponden
al valor promedio del exponente
en esta zona
delimitada. La determinación de los exponentes para órdenes más grandes que 6
es difícil para números de Reynolds altos, porque, por lo general, no es
apreciable una zona plana donde el exponente sea constante, ese fenómeno indica
que las funciones de estructura de alto orden son dominadas por eventos
inusuales, que provocan que algunas estadísticas no converjan de manera
satisfactoria [16].
Las tablas 2 y 3 muestran los resultados obtenidos mediante el método ESS para , para diferentes posiciones 𝑥/𝑀 y diferentes velocidades 𝑈0
para ambos tipos de malla. Los exponentes solo fueron medidos para la configuración de malla activa. En las
tablas no se incluye ningún valor para los exponentes de escala de orden 3,
porque para la metodología ESS el exponente para este orden es igual a 1,0.
Los valores obtenidos para son consistentes con los resultados de Toshiyuki et al. [16] y Shen et al. [17], quienes también utilizaron la metodología ESS para determinar los exponentes de escalamiento de las funciones de estructura longitudinales y transversales en sus experimentos.
La figura 9 muestra los exponentes de escalamiento graficados en función del orden n,
junto con los exponentes de escala de la teoría K41. Los símbolos ●, ▼ y ■ representan los exponentes de escala para la malla pasiva, mientras los símbolos ○, ▽ y □ representan los exponentes para la malla activa.
Se observa claramente que el alejamiento
respecto a la línea K41 es variable con el orden n considerado, cuanto más
aumenta el orden, más se aleja el exponente de escalamiento de la predicción de
Kolmogorov y es aún mayor la diferencia con respecto
a K41 para el exponente de escalamiento transversal , que toma valores significativamente menores que los exponentes
y
. Los exponentes de escala orden 2 y 4 parecen ser independientes
del número de Reynolds y del tipo de malla, porque los valores obtenidos para
estos órdenes no tienen diferencias significativas apreciables en los seis
experimentos reportados; todos los puntos de las gráficas para estos ordenes
convergen en el mismo valor. A medida que 𝑛
aumenta, se comienzan a notar diferencias cada vez más significativas entre
exponentes; para 𝑛 > 5 se evidencia una clara dependencia de los exponentes de
escala del tipo de malla y del número de Reynolds. Para la malla activa, los
exponentes de escala, longitudinales y transversales, alcanzan valores más pequeños
que los hallados para la malla pasiva.
En las figura 10 se han graficado los
exponentes ,
de orden 4, 6 y 8, en función de 𝑅𝜆 mientras se mantiene constante el grado 𝑛 de
la función de estructura, para observar más detalladamente la variación de los
exponente con el número de Reynolds y con el orden.
Los exponentes de escala longitudinales y
transversales resultaron ser mayores que la predicción de Kolmogorov
para la función de estructura de orden 2 en todas las velocidades estudiadas;
para los órdenes de 4 a 8, ,
fueron menores a
K41, siempre permaneciendo el exponente de
escala
aún más alejado de la predicción.
no mostró una fuerte
dependencia de 𝑅𝜆;
razón por la cual los resultados no fueron graficados.
Se ha estudiado el movimiento turbulento en turbulencia generada por una malla pasiva y una malla activa en un túnel de viento. La turbulencia fue estudiada para un intervalo de números de Reynolds basados en la escala de Taylor entre 100 y 1200. En general, la discusión se enmarcó en el contexto del estudio de los perfiles de velocidad, la evolución del factor de isotropía, asimetría y curtosis, las distribuciones de probabilidad y las funciones de estructura y sus exponentes de escala. Se demostró que para posiciones cercanas a la malla, 𝑥/𝑀 < 19,75, el flujo medio normalizado se caracteriza por una pobre homogeneidad transversal. La falta de homogeneidad cercana a la malla comenzó a disminuir conforme aumentaba la distancia x/M.
De acuerdo con los perfiles de velocidad media, del factor de isotropía, de la asimetría y la curtosis, el flujo fue considerado razonablemente homogéneo e isotrópico para posiciones 𝑥/𝑀 >19,75 en malla pasiva y 𝑥/ 𝑀 >29,18 en malla activa. La malla activa, debido a su relación de bloqueo constante al paso del aire, generó los números de Reynolds más altos.
Los resultados indicaron que los exponentes de escala de las funciones de estructura siguiente la metodología EES varían en función de Rλ; cuanto más aumentaba Rλ, mayor era la diferencia entre los exponentes de las funciones de estructura y los exponentes predichos por Kolmogorov en su teoría. Los exponentes de escala orden 2 y 4 fueron independientes del número de Reynolds y del tipo de malla, porque los valores obtenidos para estos órdenes no tienen diferencias significativas apreciables en los experimentos reportados.
Las líneas de tendencia de la figura 10
indicaron que tiene una fuerte dependencia del número de Reynolds, pues los
valores obtenidos para los exponentes varían significativamente entre el número
de Reynolds más alto y más bajo. Se espera que cuando Rλ
→ ∞, el exponente
sea el exponente
que más aleje de la predicción K41.
Estos resultados sugieren entonces que la postulación de Kolmogorov en su teoría K41 para los exponentes de escala es incorrecta; se propone entonces la realización de más experimentos con número de Reynolds más altos y con configuraciones de malla diferentes para validar o rechazar esta hipótesis.
El autor está agradecido con Zellman Warhaft, por su acompañamiento y guía para la realización de la fase experimental; gracias también a Juan Camilo Isaza, quien, como asesor, compartió su conocimiento y experiencia para llevar cabo esta investigación. Además, quiere dar las gracias a la Facultad de Ingeniería Mecánica y Aeroespacial de la Universidad de Cornell, por proveer los recursos e instalaciones necesarias para llevar a cabo cada uno de los experimentos; a la Universidad EAFIT, por su constante interés y apoyo en el desarrollo de esta investigación, y a la Universidad de Antioquia, por la motivación para el desarrollo de este estudio.
[1] A. N. Kolmogorov, “Dissipation of Energy in Locally Isotropic
Turbulence,” Akad. Nauk
SSSR Dokl., vol. 32, p. 16+, 1941.
[2] H. Makita and K. Sassa, “Active Turbulence
[3] J. C. Isaza, R. Salazar, and Z. Warhaft, “On gridgenerated
turbulence in the near- and far field regions,” Journal of Fluid Mechanics,
vol. 753, pp. 402– 426, 2014. doi:
10.1017/jfm.2014.375
[4] K. Yoon and Z. Warhaft, “The evolution of gridgenerated turbulence under conditions of stable thermal
stratification,” J. Fluid Mech., vol. 215, pp. 601–638,
[5]
L. Mydlarski
and Z. Warhaft, “On the onset of highReynolds-number
grid-generated wind tunnel turbulence,” J. Fluid Mech., vol. 320, pp.
331–368, 1996.
[6]
H. H. Bruun,
Hot wire anemometry : principles and signal analysis.
Oxford University Press, 1995.
[7]
P.-Å. KROGSTAD and P. A. DAVIDSON,
“Freely decaying, homogeneous turbulence generated by multiscale grids,” J.
Fluid Mech., vol. 680, pp. 417–434, 2011.
[8]
A. Sirivat
and Z. Warhaft, “The effect of a passive cross-stream
temperature gradient on the evolution of temperature variance and heat flux in
grid turbulence,” J. Fluid Mech., vol. 128, no. 1, p. 323, Mar. 1983. doi:
[9]
P.-Å. KROGSTAD and P. A. DAVIDSON,
“Is grid turbulence Saffman turbulence?,”
J. Fluid Mech., vol.
[10]
R. J. Hearst and P. Lavoie, “Decay
of turbulence generated by a square-fractal-element grid,” J. Fluid
[11]
Ö. ERTUNÇ, N. ÖZYILMAZ, H.
LIENHART, F. DURST, and K. BERONOV, “Homogeneity of turbulence generated by
static-grid structures,” J. Fluid
[12]
P. C. Valente and J. C. Vassilicos, “The decay of turbulence generated by a class
of multiscale grids,” J. Fluid Mech., vol. 687, pp. 300–340, 2011. doi:
[13]
A. A. Townsend and G. K. Batchelor, “Decay of isotropic turbulence in the initial
period,” Proc. R. Soc. London A Math. Phys. Eng. Sci., vol. 193, no.
1035, pp.
[14]
L. F. G. Simmons and C. Salter,
“Experimental investigation and analysis of the velocity variations in
turbulent flow,” Proc. R. Soc. London A Math. Phys. Eng. Sci., vol. 145,
no. 854, pp. 212–234, 1934. doi:
[15]
R. Benzi,
S. Ciliberto, R. Tripiccione,
C. Baudet, F. Massaioli,
and S. Succi, “Extended self-similarity in turbulent
flows,” Phys. Rev. E, vol. 48, no. 1, pp. R29-R32, Jul. 1993. doi: 10.1103/PhysRevE.48.R29
[16]
T. Gotoh,
D. Fukayama, and T. Nakano, “Velocity field
statistics in homogeneous steady turbulence obtained using a high-resolution
direct numerical