A Generalised Curvilinear Coordinate system-based Patch Dynamics Scheme in Equation-free Multiscale Modelling

Resumen Técnico: Un esquema de dinámica de parches basado en un sistema de coordenadas curvilíneas generalizadas en la modelización multiescala libre de ecuaciones

Planteamiento del problema
El artículo aborda los desafíos computacionales inherentes a la modelización multiescala de sistemas complejos, específicamente aquellos gobernados por ecuaciones de convección-difusión-reacción (CDR) lineales, heterogéneas y no estacionarias. Los esquemas numéricos tradicionales suelen requerir simulaciones a nivel de escala fina sobre dominios espaciales extensos e intervalos de tiempo largos, lo que conduce a costes computacionales y requisitos de memoria prohibitivos. Si bien se han desarrollado métodos multiescala "libres de ecuaciones" (equation-free), como el esquema de dinámica de parches, para cerrar la brecha entre las escalas microscópica y macroscópica mediante la simulación de solo fracciones del dominio espacio-temporal, las implementaciones existentes están mayoritariamente restringidas a dominios rectangulares con mallas uniformes y parches rectangulares. Esta limitación dificulta su aplicación a problemas de la vida real caracterizados por geometrías irregulares, requisitos de malla no uniformes (por ejemplo, para capturar capas límite) y dominios físicos complejos.

Metodología
Los autores proponen una novedosa extensión del esquema de dinámica de parches que integra coordenadas curvilíneas ortogonales generalizadas. La metodología central consiste en:

1. Transformación de Coordenadas: El dominio físico ($\Omega_p$), que puede ser no rectangular o requerir un espaciado de malla no uniforme, se mapea a un dominio computacional rectangular ($\Omega_c$) con una malla uniforme mediante una transformación suave y biunívoca. Esta transformación convierte la ecuación CDR heterogénea original en una forma donde el término de derivada mixta desaparece (condición de ortogonalidad), facilitando el cálculo en una malla estructurada.
2. Marco de Dinámica de Parches: El esquema opera en el dominio computacional utilizando una combinación del esquema de diente de hueco (gap-tooth scheme) y la integración proyectiva de escala gruesa.
   • Microsimulación: El problema microscópico se resuelve solo dentro de parches pequeños y dispersos centrados en los puntos de la malla macroscópica.
   • Levantamiento y Restricción: Un operador de levantamiento reconstruye las condiciones iniciales microscópicas dentro de los parches a partir del estado macroscópico utilizando interpolación polinómica bi-cuadrática. Inversamente, un operador de restricción promedia la solución microscópica sobre cada parche para actualizar el estado macroscópico.
   • Condiciones de Borde: Se imponen condiciones de contorno de Neumann en los bordes de los parches utilizando interpolación polinómica de Lagrange basada en los puntos de la malla macroscópica vecinos (un estenco de 9 puntos).
   • Integración Proyectiva: La derivada temporal macroscópica se estima a partir de microsimulaciones cortas, y la solución se extrapola hacia adelante sobre un gran paso de tiempo macroscópico utilizando el método de Euler hacia adelante.
3. Gestión de la Heterogeneidad: El esquema acomoda tensores de difusión dependientes del espacio, velocidades convectivas y términos fuente dependientes del espacio-tiempo. Está diseñado para problemas con coeficientes suaves y no oscilatorios.

Contribuciones Clave

• Marco Generalizado: Esta es la primera aplicación de un sistema de coordenadas curvilíneas ortogonales generalizadas dentro del esquema de dinámica de parches, lo que permite el manejo tanto de mallas no uniformes en dominios rectangulares como de dominios físicos no rectangulares.
• Parches No Uniformes y No Rectangulares: El artículo introduce el concepto de configuraciones de parches no uniformes y no rectangulares en el dominio físico, que se mapean a parches rectangulares uniformes en el dominio computacional.
• Formulación Explícita: Se proporciona una representación explícita del esquema para problemas CDR lineales, no estacionarios y heterogéneos en dos dimensiones, incluyendo la derivación de los coeficientes transformados para las ecuaciones gobernantes.
• Implementación Algorítmica: Los autores detallan el algoritmo completo, incluyendo el manejo específico de las condiciones de borde, el levantamiento y la extrapolación proyectiva en el contexto de las coordenadas curvilíneas.

Resultados
El esquema propuesto fue validado mediante dos casos de prueba numéricos distintos:

1. CDR Heterogénea Dominada por Convección en una Malla Estirada:

   • Se resolvió un problema con un dominio rectangular que presenta altos gradientes (capas límite) cerca de los bordes superior y derecho.
   • Mallas Estiradas: Utilizando un parámetro de estiramiento $\lambda = 0.1$, el esquema logró una precisión significativamente mayor en comparación con las mallas uniformes ($\lambda = 0$), particularmente en regiones de alto gradiente.
   • Rendimiento: El esquema de dinámica de parches demostró una precisión superior (errores hasta 13 veces menores en casos específicos) en comparación con una simulación de dominio completo (FDS) utilizando el método de líneas. También requirió menos tiempo de computación (reducción de aproximadamente el 40%) y memoria (reducción de aproximadamente el 10%) utilizando menos puntos de red.
   • Independencia de la Malla: Se demostró que las soluciones son independientes de la malla con una resolución de $20 \times 20$ mallas macro con 8,500 pasos de tiempo.

2. Difusión No Aximétrica en un Anillo:

   • Se resolvió un problema de difusión en un dominio anular (geometría no rectangular) utilizando un sistema de coordenadas curvilíneas ajustadas al cuerpo (body-fitted).
   • Precisión: Los resultados numéricos mostraron una excelente concordancia con la solución analítica.
   • Convergencia: El esquema demostró convergencia a medida que se refinaban los niveles de malla y tiempo macroscópicos.
   • Condiciones de Contorno: El método manejó con éxito tanto las condiciones de contorno periódicas (en la dirección angular) como las condiciones de Dirichlet (en la dirección radial).

Significancia y Reivindicaciones
El artículo afirma que el esquema propuesto ofrece una alternativa robusta y eficiente a las simulaciones de dominio completo para problemas multiescala complejos. Su significancia radica en:

• Flexibilidad Geométrica: Supera la limitación de los dominios rectangulares, permitiendo la simulación precisa de sistemas con geometrías complejas y regiones localizadas de alto gradiente sin la necesidad de una resolución global de malla fina.
• Eficiencia Computacional: Al restringir las simulaciones microscópicas a pequeños parches y utilizar la integración proyectiva, el esquema reduce significativamente el tiempo de computación y el uso de memoria manteniendo una alta precisión.
• Versatilidad: El método es aplicable tanto a problemas de "Tipo-A" (donde los detalles de la microescala se necesitan solo en regiones específicas como las capas límite) como a problemas de "Tipo-B" (donde se requieren modelos de microescala en todo el dominio debido a la ausencia de una descripción macroscópica).

Los autores señalan que la implementación actual es un desarrollo inicial centrado en ecuaciones CDR lineales, no estacionarias y heterogéneas con coeficientes suaves. Reconocen que se necesita más trabajo para extender el marco a problemas con coeficientes altamente oscilatorios, donde los errores de resonancia de las condiciones de borde artificiales podrían requerir regiones de amortiguación o núcleos de promediado.