Effects of Thermal Boundary Conditions on Natural Convection and Entropy Generation in Non-Newtonian Power-Law Fluids

Este estudio utiliza simulaciones de elementos finitos para demostrar que, en fluidos no newtonianos de ley de potencia, el comportamiento de adelgazamiento por cizallamiento mejora la transferencia de calor, mientras que las condiciones de frontera térmicas uniformes promueven una convección más fuerte y una mayor generación de entropía en comparación con el calentamiento no uniforme, ofreciendo conocimientos clave para optimizar el diseño de sistemas térmicos.

Lo esencial

Imagina que intentas comprender cómo se mueve el calor a través de un líquido espeso, como la miel o la pintura, dentro de un recipiente. Este artículo es como una receta detallada y un conjunto de experimentos para determinar exactamente cómo se comporta ese líquido cuando se calienta y cuánta "energía desperdiciada" (entropía) se genera en el proceso.

Aquí está el desglose de lo que hicieron y descubrieron los investigadores, utilizando analogías sencillas:

El Escenario: Dos Recipientes Diferentes

Los científicos examinaron dos formas específicas para observar cómo se mueve el líquido:

1. Una Caja Cuadrada: Imagina un marco de cuadro cuadrado. El fondo está caliente, los lados están fríos y la parte superior está cubierta (para que el calor no pueda escapar).
2. Una Forma de Rosquilla (Anillo): Imagina una tubería grande con una tubería más pequeña dentro de ella. La tubería interior está caliente y la tubería exterior está fría.

En ambos casos, la gravedad tira del líquido hacia abajo. Cuando el líquido cerca de la pared caliente se calienta, se vuelve más ligero e intenta flotar hacia arriba (como un globo de aire caliente), mientras que el líquido frío y pesado se hunde. Esto crea un bucle de circulación natural sin necesidad de una bomba ni un ventilador.

El Ingrediente Especial: Líquidos "Inteligentes"

La mayoría de los líquidos (como el agua) tienen un espesor constante, o viscosidad. Pero los líquidos en este estudio son No Newtonianos. Esto significa que su espesor cambia dependiendo de la velocidad a la que se mueven.

• Fluidos que se adelgazan con el cizallamiento (Shear-Thinning) (El Fluido "Fluido"): Imagina el ketchup. Cuanto más lo agitas o lo empujas, más fino y fluido se vuelve. En el artículo, estos son fluidos donde el índice de ley de potencias es menor que 1.
• Fluidos que se espesan con el cizallamiento (Shear-Thickening) (El Fluido "Rígido"): Imagina una mezcla de maicena y agua. Si la golpeas o la empujas con fuerza, instantáneamente se convierte en un bloque sólido. En el artículo, estos son fluidos donde el índice es mayor que 1.
• Newtoniano (El Fluido "Normal"): Este es el punto medio, como el agua o el aceite, donde el espesor permanece igual sin importar la velocidad a la que se mueva.

El Experimento: Cambiando la Fuente de Calor

Los investigadores no solo calentaron los recipientes uniformemente. Probaron dos formas de aplicar calor:

1. Calentamiento Uniforme: Imagina encender un calentador que calienta toda la pared inferior (o la tubería interior) por igual.
2. Calentamiento No Uniforme (Sinusoidal): Imagina un calentador que está más caliente en el medio y se vuelve más frío hacia los bordes, como una onda suave de calor.

Lo Que Descubrieron: El Baile del Calor y el Flujo

1. Cómo se Mueve el Líquido (El Flujo)

• El Fluido "Fluido" (Shear-Thinning): Cuando este fluido se calienta, se vuelve más fino y se mueve mucho más rápido. Crea bucles de remolinos fuertes y vigorosos (vórtices) que transportan el calor con gran eficiencia. Es como una licuadora de alta velocidad.
• El Fluido "Rígido" (Shear-Thickening): Cuando este fluido intenta moverse, se vuelve más grueso y resiste el movimiento. Los bucles de remolinos se vuelven débiles y lentos. El calor se mueve principalmente filtrándose lentamente a través del líquido (conducción) en lugar de fluir. Es como intentar caminar por un barro profundo.
• El Patrón de Calentamiento: Cuando el calor se aplicó uniformemente (Uniforme), el líquido creó bucles grandes y fuertes que llenaron todo el recipiente. Cuando el calor se aplicó en forma de onda (No uniforme), el líquido solo giró con fuerza justo donde el calor era más intenso, creando un "penacho" localizado de líquido caliente ascendente, mientras que el resto del recipiente permaneció relativamente tranquilo.

2. Cuánto Calor se Transfiere

• Los fluidos "Fluidos" transfirieron el calor mejor porque se movían tan rápido.
• Los fluidos "Rígidos" transfirieron el calor peor porque apenas se movían.
• Curiosamente, los fluidos "Fluidos" fueron incluso más sensibles al patrón de calentamiento. Cuando el calor era ondulado, la diferencia en el rendimiento entre los fluidos "Fluidos" y "Rígidos" se volvió aún más dramática.

3. La "Energía Desperdiciada" (Generación de Entropía)
Los investigadores también calcularon la "entropía", que es una medida de cuánta energía se desperdicia o pierde como desorden durante el proceso. Piénsalo como el "costo de fricción" de mover el calor.

• La Gran Sorpresa: Para los fluidos "Fluidos", el mayor desperdicio de energía provino del líquido frotándose contra sí mismo (disipación viscosa) mientras giraba rápidamente. Era como el motor de un coche acelerando demasiado y quemando combustible solo para girar las ruedas.
• El Cambio: A medida que el fluido se volvía "más rígido" (moviéndose hacia el lado Newtoniano o de espesamiento por cizallamiento), el desperdicio por fricción disminuyó drásticamente. Eventualmente, la fuente principal de desperdicio se convirtió en el calor mismo intentando moverse de las zonas calientes a las frías.
• El Efecto del Patrón de Calentamiento: El calentamiento "Ondulado" (No uniforme) siempre resultó en menos energía desperdiciada en total que el calentamiento "Uniforme" (Uniforme). Al concentrar el calor en un solo punto, el sistema no tuvo que trabajar tan duro para mover todo, haciéndolo ligeramente más "eficiente termodinámicamente".

La Conclusión

El estudio muestra que si quieres controlar cómo se mueve el calor a través de fluidos especiales (como pinturas, polímeros o fluidos biológicos), tienes dos palancas para accionar:

1. El Tipo de Fluido: Elegir un fluido que se vuelve más fino cuando se mueve (shear-thinning) hará que la transferencia de calor sea más rápida, pero podría crear más desperdicio por fricción.
2. El Diseño de Calentamiento: Calentar una superficie uniformemente crea corrientes fuertes y generalizadas. Calentarla en un patrón específico (como una onda) crea corrientes enfocadas y generalmente desperdicia menos energía en total.

Los investigadores construyeron una potente simulación por computadora (usando una herramienta llamada Gridap.jl) para probar estos puntos, y pusieron su código a disposición para que otros pudieran verificar su trabajo. Confirmaron que la forma en que calientas un recipiente es tan importante como el tipo de líquido que hay dentro al diseñar sistemas térmicos eficientes.

Resumen técnico

Resumen Técnico: Efectos de las Condiciones de Contorno Térmicas sobre la Convección Natural y la Generación de Entropía en Fluidos No Newtonianos de Ley de Potencia

Planteamiento del Problema
Este estudio investiga los efectos acoplados de las condiciones de contorno térmicas y la reología del fluido sobre la convección natural y la irreversibilidad termodinámica en fluidos no newtonianos de ley de potencia. La investigación se centra en dos geometrías bidimensionales distintas: una cavidad cuadrada y un anillo cilíndrico concéntrico. Si bien la convección natural en fluidos newtonianos está bien documentada, la interacción entre las condiciones de contorno térmicas que varían espacialmente (calentamiento uniforme frente a calentamiento sinusoidal no uniforme) y la viscosidad dependiente del cizallamiento de los fluidos de ley de potencia sigue siendo un área que requiere un análisis sistemático. Específicamente, el trabajo aborda la brecha en la literatura respecto a un marco unificado que evalúe simultáneamente la estructura del flujo, el rendimiento de la transferencia de calor y la generación de entropía en los regímenes de adelgazamiento por cizallamiento ($n < 1$), newtoniano ($n = 1$) y espesamiento por cizallamiento ($n > 1$).

Metodología
Los autores desarrollaron un solver numérico robusto utilizando el marco de elementos finitos Gridap.jl dentro del lenguaje de programación Julia. Las ecuaciones gobernantes para el flujo laminar bidimensional, incompresible y en estado estacionario se formularon bajo la aproximación de Boussinesq, incorporando el modelo constitutivo de ley de potencia (Ostwald–de Waele).

• Esquema Numérico: El estudio emplea el método de elementos finitos de Galerkin (MEF) con elementos de Taylor–Hood (cuadráticos para velocidad y temperatura, lineales para presión) para garantizar la estabilidad inf-sup.
• Estrategia de Solución: Debido a la fuerte no linealidad introducida por el término de viscosidad de ley de potencia, particularmente en índices extremos de ley de potencia, las ecuaciones residuales no lineales discretas se resuelven utilizando tanto el método de Newton–Raphson como un método de región de confianza (mediante el paquete NLsolve.jl) para garantizar la convergencia en casos mal condicionados.
• Validación: El solver fue validado rigurosamente contra soluciones de referencia establecidas para la convección natural tanto newtoniana como no newtoniana en geometrías cuadradas y cilíndricas. Las comparaciones incluyeron campos de isotermas, líneas de corriente, distribuciones del número de Nusselt local y contornos de generación de entropía, demostrando un buen acuerdo con los datos de la literatura (por ejemplo, Basak et al., Turan et al., Matin y Khan).
• Post-procesamiento: Se evaluaron cantidades físicas clave, incluyendo la función de corriente, la función de calor, los números de Nusselt local y promedio, y las tasas de generación de entropía. La generación de entropía se descompuso en contribuciones provenientes de la transferencia de calor ($S_\theta$) y la fricción del fluido ($S_\psi$), utilizando el número de Bejan promedio ($Be_{av}$) para cuantificar la dominancia de los mecanismos de irreversibilidad.

Resultados Clave
Las simulaciones se realizaron con un número de Prandtl de $Pr = 100$ y un número de Rayleigh de $Ra = 10^4$ para índices de ley de potencia $n = 0.6, 1.0,$ y $1.4$.

1. Influencia Reológica:

   • Adelgazamiento por cizallamiento ($n=0.6$): La reducción de la viscosidad aparente mejora la circulación impulsada por flotabilidad, lo que conduce a vórtices más fuertes, gradientes térmicos más pronunciados y tasas de transferencia de calor más altas. Sin embargo, esto también resulta en una disipación viscosa significativamente elevada.
   • Espesamiento por cizallamiento ($n=1.4$): El aumento de la viscosidad aparente suprime el movimiento del fluido, desplazando el mecanismo de transporte hacia la dominancia de la conducción, con circulación más débil y tasas de transferencia de calor más bajas.
   • Newtoniano ($n=1.0$): Exhibe un comportamiento intermedio, sirviendo como línea base para la comparación.

2. Efectos de las Condiciones de Contorno Térmicas:

   • Calentamiento Uniforme: Produce estructuras convectivas más fuertes y espacialmente distribuidas. En la cavidad cuadrada, esto conduce a células de circulación simétricas; en el anillo, genera vórtices intensos y asimétricos que giran en sentido contrario.
   • Calentamiento No Uniforme (Sinusoidal): Localiza la fuerza térmica en regiones específicas (por ejemplo, el centro de la cavidad o el sector superior del anillo). Esta localización reduce consistentemente la generación total de entropía en comparación con el calentamiento uniforme al confinar las regiones de alto gradiente y reducir la intensidad general de la fricción del fluido.

3. Generación de Entropía e Irreversibilidad:

   • Mecanismos de Dominancia: En fluidos que adelgazan por cizallamiento, la disipación viscosa es la fuente principal de irreversibilidad, particularmente en la geometría anular donde la circulación es más vigorosa. A medida que aumenta el índice de ley de potencia, la contribución de la irreversibilidad por fricción del fluido disminuye bruscamente, y la irreversibilidad por transferencia de calor se vuelve dominante.
   • Tendencias del Número de Bejan: El número de Bejan promedio ($Be_{av}$) aumenta monótonamente con $n$. Para fluidos que adelgazan por cizallamiento bajo calentamiento uniforme, $Be_{av}$ es muy bajo (por ejemplo, $\approx 0.034$ en el anillo), indicando dominancia viscosa. A medida que $n$ aumenta hasta 1.6, $Be_{av}$ se acerca a la unidad, confirmando una transición hacia una irreversibilidad dominada por la transferencia de calor.
   • Impacto de la Condición de Contorno: Si bien las condiciones de contorno térmicas alteran significativamente la magnitud y la distribución espacial de la generación de entropía (el calentamiento no uniforme reduce la entropía total), la naturaleza del mecanismo de irreversibilidad dominante está gobernada principalmente por el índice de ley de potencia.

Significado y Afirmaciones
El artículo afirma que el diseño adecuado de las condiciones de contorno térmicas, cuando se considera junto con la reología del fluido, proporciona una vía efectiva para controlar el rendimiento de la transferencia de calor y minimizar las pérdidas termodinámicas en sistemas de convección no newtonianos. El estudio demuestra que el calentamiento sinusoidal no uniforme puede utilizarse para localizar la fuerza térmica y reducir la generación total de entropía sin alterar fundamentalmente las características reológicas del régimen de flujo. Además, el trabajo establece que el índice de ley de potencia es el parámetro crítico que determina si un sistema está limitado por la disipación viscosa o por la irreversibilidad de la transferencia de calor. Los autores concluyen que sus hallazgos ofrecen directrices fiables para la optimización de sistemas térmicos eficientes energéticamente que involucran fluidos complejos, como el procesamiento de polímeros y aplicaciones de enfriamiento industrial, al proporcionar una comprensión integral de la interacción entre las condiciones de contorno y la reología no newtoniana. El código fuente y los metadatos se ponen a disposición del público para garantizar la reproducibilidad.