Gracias por visitar Nature.com. La versión del navegador que está utilizando tiene soporte limitado para CSS. Para obtener la mejor experiencia, le recomendamos que utilice un navegador actualizado (o desactive el modo de compatibilidad en Internet Explorer). Mientras tanto, para garantizar el soporte continuo, mostraremos el sitio sin estilos ni JavaScript.
Se realizaron experimentos en un canal rectangular bloqueado por líneas transversales de cuatro varillas cilíndricas inclinadas. Se midieron la presión en la superficie central de la varilla y la caída de presión a través del canal variando el ángulo de inclinación de la varilla. Se probaron tres conjuntos de varillas de diferente diámetro. Los resultados de la medición se analizaron utilizando el principio de conservación del momento y consideraciones semiempíricas. Se generaron varios conjuntos invariantes de parámetros adimensionales que relacionan la presión en ubicaciones críticas del sistema con las dimensiones características de la varilla. Se encontró que el principio de independencia se cumple para la mayoría de los números de Euler que caracterizan la presión en diferentes ubicaciones, es decir, si la presión es adimensional utilizando la proyección de la velocidad de entrada normal a la varilla, el conjunto es independiente del ángulo de inclinación. La correlación semiempírica resultante se puede utilizar para diseñar hidráulica similar.
Muchos dispositivos de transferencia de calor y masa constan de un conjunto de módulos, canales o celdas a través de los cuales pasan fluidos en estructuras internas más o menos complejas, como varillas, amortiguadores, insertos, etc. Más recientemente, ha habido un renovado interés en obtener una mejor comprensión de los mecanismos que vinculan la distribución de la presión interna y las fuerzas en los componentes internos complejos con la caída de presión general del módulo. Entre otras cosas, este interés se ha visto impulsado por las innovaciones en la ciencia de los materiales, la expansión de las capacidades computacionales para simulaciones numéricas y la creciente miniaturización de los dispositivos. Estudios experimentales recientes de la distribución interna de la presión y las pérdidas incluyen canales rugosos mediante nervaduras de varias formas 1 , celdas de reactores electroquímicos 2 , constricción capilar 3 y materiales de marco de celosía 4 .
Las estructuras internas más comunes son posiblemente varillas cilíndricas a través de módulos unitarios, ya sea agrupados o aislados. En los intercambiadores de calor, esta configuración es típica en el lado de la carcasa. La caída de presión del lado de la carcasa está relacionada con el diseño de intercambiadores de calor como generadores de vapor, condensadores y evaporadores. En un estudio reciente, Wang et al. 5 encontraron estados de flujo de reacoplamiento y codesprendimiento en una configuración en tándem de varillas. Liu et al.6 midieron la caída de presión en canales rectangulares con haces de tubos dobles en forma de U incorporados con diferentes ángulos de inclinación y calibraron un modelo numérico que simula haces de varillas con medios porosos.
Como era de esperar, hay una serie de factores de configuración que afectan el rendimiento hidráulico de un banco de cilindros: tipo de disposición (p. ej., escalonada o en línea), dimensiones relativas (p. ej., paso, diámetro, longitud) y ángulo de inclinación, entre otros. Varios autores se centraron en encontrar criterios adimensionales para guiar los diseños para capturar los efectos combinados de los parámetros geométricos. En un estudio experimental reciente, Kim et al. 7 propusieron un modelo de porosidad eficaz utilizando la longitud de la celda unitaria como parámetro de control, utilizando arreglos en tándem y escalonados y números de Reynolds entre 103 y 104. Snarski8 estudió cómo el espectro de potencia, de acelerómetros e hidrófonos unidos a un cilindro en un túnel de agua, varía con la inclinación de la dirección del flujo. Marino et al. 9 estudiaron la distribución de la presión de la pared alrededor de una varilla cilíndrica en un flujo de aire de guiñada. Mityakov et al. 10 trazaron el campo de velocidad después de un cilindro guiñado utilizando PIV estéreo. Alam et al. 11 realizaron un estudio exhaustivo de cilindros en tándem, centrándose en los efectos del número de Reynolds y la relación geométrica en el desprendimiento de vórtices. Pudieron identificar cinco estados, a saber, bloqueo, bloqueo intermitente, sin bloqueo, bloqueo subarmónico y estados de reinserción de la capa de corte. Estudios numéricos recientes han señalado la formación de estructuras de vórtice en el flujo a través de cilindros de guiñada restringida.
En general, se espera que el rendimiento hidráulico de una celda unitaria dependa de la configuración y geometría de la estructura interna, generalmente cuantificada mediante correlaciones empíricas de mediciones experimentales específicas. En muchos dispositivos compuestos por componentes periódicos, los patrones de flujo se repiten en cada celda y, por lo tanto, la información relacionada con las celdas representativas se puede utilizar para expresar el comportamiento hidráulico general de la estructura a través de modelos multiescala. En estos casos simétricos, a menudo se puede reducir el grado de especificidad con el que se aplican los principios generales de conservación. Un ejemplo típico es la ecuación de descarga para una placa de orificio 15. En el caso especial de las varillas inclinadas, ya sea en flujo confinado o abierto, un criterio interesante a menudo citado en la literatura y utilizado por los diseñadores es la magnitud hidráulica dominante (por ejemplo, caída de presión, fuerza, frecuencia de desprendimiento de vórtices, etc.) ) para contactar.) al componente de flujo perpendicular al eje del cilindro. Esto a menudo se conoce como el principio de independencia y asume que la dinámica del flujo está impulsada principalmente por el componente normal de entrada y que el efecto del componente axial alineado con el eje del cilindro es insignificante. Aunque no existe consenso en la literatura sobre el rango de validez de este criterio, en muchos casos proporciona estimaciones útiles dentro de las incertidumbres experimentales típicas de las correlaciones empíricas. Estudios recientes sobre la validez del principio independiente incluyen la vibración inducida por vórtices16 y el arrastre promedio monofásico y bifásico417.
En el presente trabajo se presentan los resultados del estudio de la presión interna y la caída de presión en un canal con una línea transversal de cuatro varillas cilíndricas inclinadas. Se miden tres conjuntos de varillas con diferentes diámetros, cambiando el ángulo de inclinación. El objetivo general es investigar el mecanismo por el cual la distribución de presión en la superficie de la varilla se relaciona con la caída de presión general en el canal. Se analizan los datos experimentales aplicando la ecuación de Bernoulli y el principio de conservación del momento para evaluar la validez del principio de independencia. Finalmente, se generan correlaciones semiempíricas adimensionales que se pueden utilizar para diseñar dispositivos hidráulicos similares.
La configuración experimental consistió en una sección de prueba rectangular que recibió un flujo de aire proporcionado por un ventilador axial. La sección de prueba contiene una unidad que consta de dos varillas centrales paralelas y dos medias varillas incrustadas en las paredes del canal, como se muestra en la Fig. 1e, todas del mismo diámetro. Las figuras 1a-e muestran la geometría y las dimensiones detalladas de cada parte de la configuración experimental. La figura 3 muestra la configuración del proceso.
a Sección de entrada (longitud en mm).Creada con Openscad 2021.01, openscad.org.Sección principal de pruebas (longitud en mm).Creada con Openscad 2021.01, openscad.org c Vista transversal de la sección principal de pruebas (longitud en mm).Creada con Openscad 2021.01, openscad.org d Sección de exportación (longitud en mm).Creada con Openscad 2021.01, vista ampliada de la sección de pruebas de openscad.org e.Creada con Openscad 2021.01, openscad.org.
Se probaron tres juegos de varillas de diferentes diámetros. La Tabla 1 enumera las características geométricas de cada caso. Las varillas están montadas en un transportador de modo que su ángulo relativo a la dirección del flujo puede variar entre 90° y 30° (Figuras 1b y 3). Todas las varillas están hechas de acero inoxidable y están centradas para mantener la misma distancia de separación entre ellas. La posición relativa de las varillas está fijada por dos espaciadores ubicados fuera de la sección de prueba.
El caudal de entrada de la sección de prueba se midió con un venturi calibrado, como se muestra en la Figura 2, y se monitoreó utilizando un DP Cell Honeywell SCX. La temperatura del fluido en la salida de la sección de prueba se midió con un termómetro PT100 y se controló a 45 ± 1 ° C. Para asegurar una distribución de velocidad planar y reducir el nivel de turbulencia en la entrada del canal, el flujo de agua entrante se fuerza a través de tres pantallas metálicas. Se utilizó una distancia de asentamiento de aproximadamente 4 diámetros hidráulicos entre la última pantalla y la varilla, y la longitud de la salida fue de 11 diámetros hidráulicos.
Diagrama esquemático del tubo Venturi utilizado para medir la velocidad del flujo de entrada (longitud en milímetros).Creado con Openscad 2021.01, openscad.org.
Monitorizar la presión en una de las caras de la varilla central mediante una toma de presión de 0,5 mm en el plano medio de la sección de prueba. El diámetro de la toma corresponde a un ángulo de 5°, por lo que la precisión angular es de aproximadamente 2°. La varilla monitorizada se puede girar sobre su eje, como se muestra en la Figura 3. La diferencia entre la presión superficial de la varilla y la presión a la entrada de la sección de prueba se mide con una celda DP diferencial Honeywell serie SCX. Esta diferencia de presión se mide para cada disposición de barra, variando la velocidad del flujo, el ángulo de inclinación \(\alpha \) y el ángulo azimutal \(\theta \).
Ajustes de flujo. Las paredes del canal se muestran en gris. El flujo fluye de izquierda a derecha y está bloqueado por la varilla. Tenga en cuenta que la vista "A" es perpendicular al eje de la varilla. Las varillas exteriores están semiincrustadas en las paredes laterales del canal. Se utiliza un transportador para medir el ángulo de inclinación \(\alpha \). Creado con Openscad 2021.01, openscad.org.
El propósito del experimento es medir e interpretar la caída de presión entre las entradas del canal y la presión en la superficie de la varilla central, \(\theta\) y \(\alpha\) para diferentes acimutes y buzamientos. Para resumir los resultados, la presión diferencial se expresará en forma adimensional como el número de Euler:
donde \(\rho \) es la densidad del fluido, \({u}_{i}\) es la velocidad media de entrada, \({p}_{i}\) es la presión de entrada y \({p }_{ w}\) es la presión en un punto dado en la pared de la varilla. La velocidad de entrada se fija dentro de tres rangos diferentes determinados por la apertura de la válvula de entrada. Las velocidades resultantes varían de 6 a 10 m/s, correspondientes al número de Reynolds del canal, \(Re\equiv {u}_{i}H/\nu \) (donde \(H\) es la altura del canal y \(\nu \) es la viscosidad cinemática) entre 40.000 y 67.000. El número de Reynolds de la varilla (\(Re\equiv {u}_{i}d/\nu \)) varía de 2500 a 6500. La intensidad de la turbulencia estimada por la desviación estándar relativa de las señales registradas en el El venturi es del 5% en promedio.
La Figura 4 muestra la correlación de \({Eu}_{w}\) con el ángulo azimutal \(\theta \), parametrizado por tres ángulos de inclinación, \(\alpha \) = 30°, 50° y 70°. Las mediciones se dividen en tres gráficos según el diámetro de la varilla. Se puede ver que dentro de la incertidumbre experimental, los números de Euler obtenidos son independientes del caudal. La dependencia general de θ sigue la tendencia habitual de la presión de la pared alrededor del perímetro de un obstáculo circular. En ángulos orientados hacia el flujo, es decir, θ de 0 a 90°, la presión de la pared de la varilla disminuye, alcanzando un mínimo a 90°, que corresponde al espacio entre las varillas donde la velocidad es mayor debido a las limitaciones del área de flujo. Posteriormente, hay una recuperación de presión de θ de 90° a 100°, después de lo cual la presión permanece uniforme debido a la separación de la capa límite trasera de la pared de la varilla. Nótese que no hay cambio en el ángulo de mínimo presión, lo que sugiere que las posibles perturbaciones de las capas de cizallamiento adyacentes, como los efectos Coanda, son secundarias.
Variación del número de Euler de la pared alrededor de la varilla para diferentes ángulos de inclinación y diámetros de varilla. Creado con Gnuplot 5.4, www.gnuplot.info.
A continuación, analizamos los resultados basándonos en el supuesto de que los números de Euler solo se pueden estimar mediante parámetros geométricos, es decir, las razones de longitud de las características \(d/g\) y \(d/H\) (donde \(H\) es la altura del canal) y la inclinación \(\alpha \). Una regla práctica popular establece que la fuerza estructural del fluido en la barra de guiñada está determinada por la proyección de la velocidad de entrada perpendicular al eje de la barra, \({u}_{n}={u}_{i}\mathrm {sin} \alpha \) . Esto a veces se denomina principio de independencia. Uno de los objetivos del siguiente análisis es examinar si este principio se aplica a nuestro caso, donde el flujo y las obstrucciones están confinados dentro de canales cerrados.
Consideremos la presión medida en la parte delantera de la superficie intermedia de la varilla, es decir θ = 0. Según la ecuación de Bernoulli, la presión en esta posición\({p}_{o}\) satisface:
Donde \({u}_{o}\) es la velocidad del fluido cerca de la pared de la varilla en θ = 0, y asumimos pérdidas irreversibles relativamente pequeñas. Nótese que la presión dinámica es independiente en el término de energía cinética. Si \({u}_{o}\) está vacío (es decir, en condición estancada), los números de Euler deberían estar unificados. Sin embargo, se puede observar en la Figura 4 que en \(\theta =0\) el \({Eu}_{w}\) resultante es cercano, pero no exactamente igual a este valor, especialmente para ángulos de inclinación mayores. Esto sugiere que la velocidad en la superficie de la varilla no se desvanece en \(\theta =0\), lo cual puede ser suprimido por la desviación ascendente de las líneas de corriente creadas por la inclinación de la varilla. Dado que el flujo está confinado a la parte superior e inferior de la sección de prueba, esta desviación debería crear una recirculación secundaria, aumentando la velocidad axial en la parte inferior y disminuyendo la velocidad en la parte superior. Suponiendo que la magnitud de la desviación anterior es la proyección de la velocidad de entrada. en el eje (es decir, \({u}_{i}\mathrm{cos}\alpha \)), el resultado del número de Euler correspondiente es:
La Figura 5 compara las ecuaciones.(3) Muestra una buena concordancia con los datos experimentales correspondientes. La desviación media fue del 25% y el nivel de confianza fue del 95%. Nótese que la ecuación.(3) De acuerdo con el principio de independencia. Asimismo, la Figura 6 muestra que el número de Euler corresponde a la presión en la superficie trasera de la varilla, \({p}_{180}\), y a la salida del segmento de prueba, \({p}_{e}\), También sigue una tendencia proporcional a \({\mathrm{sin}}^{2}\alpha \) . En ambos casos, sin embargo, el coeficiente depende del diámetro de la varilla, lo cual es razonable ya que este último determina el área obstaculizada. Esta característica es similar a la caída de presión de una placa de orificio, donde el canal de flujo se reduce parcialmente en ubicaciones específicas. En esta sección de prueba, el papel del orificio lo desempeña el espacio entre las varillas. En este caso, la presión cae sustancialmente en el estrangulamiento y se recupera parcialmente a medida que se expande. hacia atrás. Considerando la restricción como un bloqueo perpendicular al eje de la varilla, la caída de presión entre la parte delantera y trasera de la varilla se puede escribir como 18:
donde \({c}_{d}\) es un coeficiente de arrastre que explica la recuperación de presión parcial entre θ = 90° y θ = 180°, y \({A}_{m}\) y \ ({A}_{f}\) es la sección transversal libre mínima por unidad de longitud perpendicular al eje de la varilla, y su relación con el diámetro de la varilla es \({A}_{f}/{A}_{m}=\ ​​Left (g+d\right)/g\). Los números de Euler correspondientes son:
Número de Euler de Wall en \(\theta =0\) en función de la inclinación. Esta curva corresponde a la ecuación.(3). Creado con Gnuplot 5.4, www.gnuplot.info.
Los números de Euler de Wall cambian, en \(\theta =18{0}^{o}\) (signo lleno) y salida (signo vacío) con inmersión. Estas curvas corresponden al principio de independencia, es decir \(Eu\propto {\mathrm{sin}}^{2}\alpha \). Creado con Gnuplot 5.4, www.gnuplot.info.
La Figura 7 muestra la dependencia de \({Eu}_{0-180}/{\mathrm{sin}}^{2}\alpha \) en \(d/g\), mostrando la consistencia extremadamente Buena.(5).El coeficiente de arrastre obtenido es \({c}_{d}=1.28\pm 0.02\) con un nivel de confianza del 67%.Asimismo, el mismo gráfico también muestra que la caída de presión total entre la entrada y la salida de la sección de prueba sigue una tendencia similar, pero con diferentes coeficientes que tienen en cuenta la recuperación de presión en el espacio posterior entre la barra y la salida del canal.El coeficiente de arrastre correspondiente es \({c}_{d}=1.00\pm 0.05\) con un nivel de confianza del 67%.
El coeficiente de arrastre está relacionado con la caída de presión \(d/g\) delante y detrás de la varilla\(\left({Eu}_{0-180}\right)\) y la caída de presión total entre la entrada y la salida del canal. El área gris es la banda de confianza del 67 % para la correlación. Creado con Gnuplot 5.4, www.gnuplot.info.
La presión mínima \({p}_{90}\) en la superficie de la varilla en θ = 90° requiere un manejo especial. Según la ecuación de Bernoulli, a lo largo de la línea de corriente a través del espacio entre las barras, la presión en el centro\({p}_{g}\) y la velocidad\({u}_{g}\) en el espacio entre las barras (coincide con el punto medio del canal) están relacionadas con los siguientes factores:
La presión \({p}_{g}\) puede relacionarse con la presión superficial de la varilla en θ = 90° integrando la distribución de presión en el espacio que separa la varilla central entre el punto medio y la pared (véase la Figura 8). El balance de potencias da 19:
donde \(y\) es la coordenada normal a la superficie de la varilla desde el punto central del espacio entre las varillas centrales, y \(K\) es la curvatura de la línea de corriente en la posición \(y\). Para la evaluación analítica de la presión sobre la superficie de la varilla, suponemos que \({u}_{g}\) es uniforme y \(K\left(y\right)\) es lineal. Estas suposiciones se han verificado mediante cálculos numéricos. En la pared de la varilla, la curvatura está determinada por la sección de elipse de la varilla en el ángulo \(\alpha \), es decir, \(K\left(g/2\right)=\left(2/d\right){\ mathrm{sin} }^{2}\alpha \) (véase la Figura 8). Entonces, con respecto a la curvatura de la línea de corriente que se desvanece en \(y=0\) debido a la simetría, la curvatura en la coordenada universal \(y\) viene dada por:
Vista transversal característica, frontal (izquierda) y superior (abajo).Creado con Microsoft Word 2019.
Por otra parte, por conservación de la masa, la velocidad promedio en un plano perpendicular al flujo en el lugar de medición \(\langle {u}_{g}\rangle \) está relacionada con la velocidad de entrada:
donde \({A}_{i}\) es el área de flujo transversal en la entrada del canal y \({A}_{g}\) es el área de flujo transversal en la ubicación de medición (ver Fig. 8) respectivamente por:
Nótese que \({u}_{g}\) no es igual a \(\langle {u}_{g}\rangle \). De hecho, la Figura 9 representa la relación de velocidad \({u}_{g}/\langle {u}_{g}\rangle \), calculada por la ecuación.(10)–(14), graficada de acuerdo con la relación \(d/g\). A pesar de cierta discreción, se puede identificar una tendencia, que se aproxima mediante un polinomio de segundo orden:
La relación de las velocidades máxima\({u}_{g}\) y promedio\(\langle {u}_{g}\rangle \) de la sección transversal del centro del canal\(.\) Las curvas sólidas y discontinuas corresponden a las ecuaciones.(5) y el rango de variación de los coeficientes correspondientes\(\pm 25\%\).Creado con Gnuplot 5.4, www.gnuplot.info.
La figura 10 compara \({Eu}_{90}\) con los resultados experimentales de la ecuación.(16). La desviación relativa media fue del 25% y el nivel de confianza fue del 95%.
El número de Euler de Wall en \(\theta ={90}^{o}\).Esta curva corresponde a la ecuación.(16).Creado con Gnuplot 5.4, www.gnuplot.info.
La fuerza neta \({f}_{n}\) que actúa sobre la barra central perpendicular a su eje se puede calcular integrando la presión sobre la superficie de la barra de la siguiente manera:
donde el primer coeficiente es la longitud de la varilla dentro del canal, y la integración se realiza entre 0 y 2π.
La proyección de \({f}_{n}\) en la dirección del flujo de agua debe coincidir con la presión entre la entrada y la salida del canal, a menos que la fricción sea paralela a la varilla y menor debido al desarrollo incompleto de la sección posterior. El flujo de momento está desequilibrado. Por lo tanto,
La figura 11 muestra un gráfico de las ecuaciones.(20) mostró un buen acuerdo para todas las condiciones experimentales. Sin embargo, hay una ligera desviación del 8% hacia la derecha, que se puede atribuir y utilizar como una estimación del desequilibrio de momento entre la entrada y la salida del canal.
Balance de potencia del canal.La línea corresponde a la ecuación.(20).El coeficiente de correlación de Pearson fue 0,97.Creado con Gnuplot 5.4, www.gnuplot.info.
Variando el ángulo de inclinación de la varilla, se midieron la presión en la pared de la superficie de la varilla y la caída de presión en el canal con las líneas transversales de las cuatro varillas cilíndricas inclinadas. Se probaron tres conjuntos de varillas de diferente diámetro. En el rango de números de Reynolds probado, entre 2500 y 6500, el número de Euler es independiente del caudal. La presión en la superficie central de la varilla sigue la tendencia habitual observada en los cilindros, siendo máxima en la parte delantera y mínima en el espacio lateral entre las varillas, recuperándose en la parte trasera debido a la separación de la capa límite.
Se analizan datos experimentales utilizando consideraciones de conservación de momento y evaluaciones semiempíricas para encontrar números adimensionales invariantes que relacionan los números de Euler con las dimensiones características de los canales y las varillas. Todas las características geométricas del bloqueo están completamente representadas por la relación entre el diámetro de la varilla y el espacio entre las varillas (lateralmente) y la altura del canal (vertical).
Se ha comprobado que el principio de independencia se cumple para la mayoría de los números de Euler que caracterizan la presión en diferentes puntos; es decir, si la presión es adimensional utilizando la proyección de la velocidad de entrada normal a la varilla, el conjunto es independiente del ángulo de inclinación. Además, la característica está relacionada con la masa y el momento del flujo. Las ecuaciones de conservación son consistentes y respaldan el principio empírico mencionado. Solo la presión superficial de la varilla en el espacio entre ellas se desvía ligeramente de este principio. Se generan correlaciones semiempíricas adimensionales que pueden utilizarse para diseñar dispositivos hidráulicos similares. Este enfoque clásico es consistente con aplicaciones similares de la ecuación de Bernoulli en hidráulica y hemodinámica, reportadas recientemente (20,21,22,23,24).
Un resultado particularmente interesante surge del análisis de la caída de presión entre la entrada y la salida de la sección de prueba. Dentro de la incertidumbre experimental, el coeficiente de arrastre resultante es igual a la unidad, lo que indica la existencia de los siguientes parámetros invariantes:
Nótese el tamaño \(\left(d/g+2\right)d/g\) en el denominador de la ecuación.(23) es la magnitud entre paréntesis en la ecuación.(4), de lo contrario se puede calcular con la sección transversal mínima y libre perpendicular a la varilla, \({A}_{m}\) y \({A}_{f}\). Esto sugiere que se supone que los números de Reynolds permanecen dentro del rango del estudio actual (40 000-67 000 para canales y 2500-6500 para varillas). Es importante tener en cuenta que si hay una diferencia de temperatura dentro del canal, puede afectar la densidad del fluido. En este caso, el cambio relativo en el número de Euler se puede estimar multiplicando el coeficiente de expansión térmica por la diferencia de temperatura máxima esperada.
Ruck, S., Köhler, S., Schlindwein, G. y Arbeiter, F. Medidas de transferencia de calor y caída de presión en un canal rugoso mediante nervaduras de diferentes formas en la pared. Expert. Heat Transfer 31, 334–354 (2017).
Wu, L., Arenas, L., Graves, J. y Walsh, F. Caracterización de celdas de flujo: visualización del flujo, caída de presión y transporte de masa en electrodos bidimensionales en canales rectangulares. J. Electrochemistry. Socialist Party. 167, 043505 (2020).
Liu, S., Dou, X., Zeng, Q. y Liu, J. Parámetros clave del efecto Jamin en capilares con secciones transversales constreñidas. J. Gasoline.science.Britain.196, 107635 (2021).