Lattice Boltzmann modeling for conduction heat transfer based on the Qadyan method

Resumen Técnico: Modelado de Lattice Boltzmann para la Conducción de Transferencia de Calor Basado en el Método de Qadyan

Planteamiento del Problema
El artículo aborda la simulación de la transferencia de calor por conducción utilizando el Método de Lattice Boltzmann (LBM). Aunque el LBM es eficaz para la dinámica de fluidos, su aplicación a la conducción de calor enfrenta limitaciones respecto a la uniformidad de la malla y la estabilidad numérica. Los esquemas estándar de LBM, particularmente aquellos basados en la aproximación de Bhatnagar–Gross–Krook (BGK), suelen estar restringidos a mallas uniformes con pasos espaciales iguales. Además, los métodos de discretización espacial estándar, como los esquemas de diferencias finitas de segundo orden, pueden inducir oscilaciones numéricas en la solución. El estudio busca evaluar estrategias alternativas de discretización espacial para mejorar la precisión y la estabilidad sin introducir complejidad adicional al marco estándar de LBM o recurrir a mallas no uniformes.

Metodología
La investigación implementa y compara tres enfoques distintos para resolver la ecuación de lattice Boltzmann (LBE) para la conducción de calor:

1. LBM Estándar: Utiliza modelos de red convencionales (D2Q9, D2Q8 y D2Q5) con procesos estándar de colisión-transporte (collision-streaming). El estudio señala que para problemas puramente difusivos, los modelos con menos conjuntos de velocidades (como D2Q8 y D2Q5) son suficientes, ya que la difusión de calor carece de preferencia direccional.
2. Método de Lattice Boltzmann de Diferencias Finitas (FDLBM): Este enfoque discretiza el término convectivo de la LBE utilizando esquemas de diferencias finitas. Específicamente, emplea una diferencia hacia adelante de primer orden para las derivadas temporales y un esquema de upwind (corriente arriba) de primer orden para las derivadas espaciales para evitar oscilaciones, aunque esto limita la precisión espacial al primer orden cerca de las fronteras.
3. El Método de Qadyan: Un nuevo enfoque propuesto por Ghadyani et al. que combina esquemas semidiscretos con el Método de Cuadratura Diferencial (DQM). En este método, las derivadas espaciales en la LBE se discretizan utilizando coeficientes de ponderación de DQM. Para replicar el efecto de upwind inherente a los problemas de advección-difusión, el método modifica los coeficientes de ponderación para asignar un mayor peso a los nodos de corriente arriba que a los de corriente abajo. Esto permite la resolución de la LBE como un sistema de ecuaciones diferenciales ordinarias sin estar ligado a un tipo específico de modelo de velocidad discreta.

El estudio valida estos métodos utilizando la conducción de calor transitoria unidimensional en una placa con condiciones de contorno tanto de Dirichlet (temperatura constante) como de Neumann (aislamiento/flujo cero). Las soluciones analíticas derivadas de la ecuación de calor sirven como referentes. Se calculan las normas de error ($L_1$, $L_2$ y $L_\infty$) para evaluar la precisión, y se resuelve un problema de difusión en una placa 2D para demostrar la aplicabilidad del método en casos multidimensionales.

Contribuciones Clave y Resultados

• Precisión vs. Estabilidad: Los resultados indican que, si bien el LBM tradicional ofrece una alta precisión, el método Qadyan propuesto proporciona una alternativa robusta que mejora la precisión sobre el enfoque FDLBM, particularmente para condiciones de contorno de Neumann. Para condiciones de Dirichlet, el LBM D2Q5 estándar mostró las normas de error más bajas, pero el método Qadyan superó significativamente a las formulaciones FDLBM.
• Independencia de la Malla: El método Qadyan logra ampliar el rango de precisión y estabilidad de las soluciones, aumentando la exactitud sin introducir complejidad adicional al método de lattice Boltzmann estándar ni recurrir a mallas no uniformes.
• Costo Computacional: Un hallazgo crítico es la relación entre precisión y eficiencia. El método Qadyan conlleva un mayor costo computacional y un tiempo de simulación más largo en comparación con tanto el LBM estándar como el FDLBM, debido al incremento en el número de operaciones de cálculo requeridas para los coeficientes de ponderación de DQM.
• Sensibilidad de Parámetros: El estudio establece que, para el LBM estándar y el FDLBM, los valores de referencia para la velocidad de la red dependen del tipo de condición de contorno (temperatura vs. flujo). En contraste, el método Qadyan utiliza un valor de velocidad de red consistente ($c=0.015$) independientemente del tipo de condición de contorno, simplificando la selección de parámetros.
• Validación 2D: En una simulación de placa cuadrada bidimensional, el método Qadyan produjo contornos de isotermas y distribuciones de temperatura de la línea central que mostraron una buena concordancia con los resultados del método de diferencias finitas, demostrando su capacidad para problemas de transferencia de calor multidimensionales.

Significancia y Reivindicaciones
El artículo afirma que el método Qadyan es un enfoque adecuado para una amplia gama de problemas de transferencia de calor por conducción. Su principal significancia radica en su capacidad para aumentar la precisión numérica sin introducir la complejidad de las mallas no uniformes o abandonar el marco de LBM. Los autores sostienen que el método es robusto y sencillo de implementar. Sin embargo, el documento mantiene una postura modesta respecto a la eficiencia, reconociendo que el método es computacionalmente más costoso que los enfoques tradicionales. El estudio concluye que, si bien el método Qadyan requiere un mayor tiempo de cómputo, su mejora en la precisión y estabilidad para condiciones de contorno específicas justifica su aplicación, ofreciendo una alternativa competitiva a las formulaciones estándar de LBM para simulaciones de transferencia de calor por conducción.