Resumen Técnico: Advección de Elementos Espectrales Miméticos

Planteamiento del Problema

El artículo aborda la discretización numérica de la ecuación de advección lineal para formas diferenciales en dominios acotados. El problema clásico de advección escalar, $\partial_t \rho + \nabla \cdot (v\rho) = 0$, se generaliza dentro de un marco de geometría diferencial como:
$$\frac{\partial \alpha^{(k)}}{\partial t} + \mathcal{L}_v \alpha^{(k)} = 0$$
donde $\alpha^{(k)}$ es una $k$-forma diferencial (que representa escalares, vectores o cantidades de mayor dimensión dependiendo de $k$), $v$ es un campo de velocidad Lipschitz continuo prescrito, y $\mathcal{L}_v$ es la derivada de Lie. El desafío radica en desarrollar una discretización que preserve las estructuras geométricas subyacentes de las leyes físicas, garantice la conservación de masa local y logre una precisión espectral, distinguiendo entre las cantidades evaluadas en instantes de tiempo específicos y aquellas integradas sobre intervalos de tiempo.

Metodología

El método propuesto, Advección de Elementos Espectrales Miméticos, extiende el marco mimético establecido en el trabajo previo [4] incorporando la derivada de Lie mediante la fórmula de homotopía de Cartan. El enfoque se basa en tres pilos fundamentales:

1. Fundamento de Geometría Diferencial
El método utiliza el lenguaje de las formas diferenciales. La derivada de Lie $\mathcal{L}_v$ se expresa mediante la fórmula de homotopía de Cartan:
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
donde $d$ es la derivada exterior e $\iota_v$ es el producto interior. Un aspecto crucial de la formulación es la relación de adjunción entre el producto interior y el producto exterior, definida mediante el producto interno $L^2$:
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
Esta dualidad permite que las cantidades físicas representadas por productos interiores se mapeen a sus formas diferenciales duales de orden 1.

2. Discretización Espacial (Elementos Espectrales Miméticos)
El dominio espacial se discretiza utilizando un complejo de celdas compuesto por puntos, segmentos de línea, superficies y volúmenes. El espacio de $k$-formas suaves se aproxima mediante un espacio de dimensión finita $\Lambda_h^k$ generado por formas base $\epsilon_i^{(k)}$.

• Construcción de la Base: Las funciones de base se construyen utilizando expansiones polinómicas por partes (productos tensoriales) en elementos cuadriláteros. En 1D, las 0-formas se representan mediante polinomios de Lagrange evaluados en nodos de Gauss-Lobatto, mientras que las 1-formas se representan mediante "polinomios de arista" derivados para satisfacer propiedades integrales específicas.
• Proyección Conmutativa: Se define un operador de proyección $\pi_h$ tal que conmuta con la derivada exterior ($\pi_h d = d \pi_h$). Esto asegura que la derivada exterior discreta $d_h$ esté representada por matrices de incidencia que contienen únicamente valores $\{0, 1, -1\}$, preservando la estructura topológica de los operadores continuos.

3. Discretización Temporal (Integrador Mimético Canónico)
La integración temporal emplea un integrador de orden arbitrario simpléctico derivado de integradores de colocación de Gauss canónicos [5].

• Nodos Temporales Escalonados: El método distingue entre dos conjuntos de nodos temporales:
  • Nodos de Gauss-Lobatto ($t_k$): Donde se definen las variables de la solución (por ejemplo, $\rho$).
  • Nodos de Gauss ($\tilde{t}_q$): Donde se evalúan las derivadas temporales (flujos).
• Evolución Discreta: Este escalonamiento resulta en un integrador discreto donde el cambio en la solución sobre un intervalo se iguala al flujo evaluado en los nodos de Gauss internos. Esta estructura emula los métodos de tipo leap-frog y la regla del punto medio implícito, asegurando propiedades simplécticas.

4. Producto Interior Discreto
El producto interior $\iota_v$ se discretiza imponiendo la dualidad de emparejamiento (Ec. 16) en el entorno discreto. Esto conduce a un sistema donde los flujos discretos se computan resolviendo un sistema lineal que involucra los productos internos de las formas base y el campo de velocidad.

Contribuciones Clave

• Consistencia Geométrica: El esquema incorpora explícitamente la derivada de Lie usando la fórmula de Cartan, asegurando que la discretización respete la representación libre de métrica de los operadores diferenciales y sus adjuntos de Hilbert.
• Precisión Espectral: El método utiliza funciones de base polinómicas de alto orden, permitiendo la convergencia espectral en el espacio.
• Conservación de Masa Local: El uso de matrices de incidencia y la estructura específica de la derivada exterior discreta garantiza la conservación de la masa local.
• Distinción de Paso de Tiempo: El marco separa rigurosamente las cantidades evaluadas en instantes de tiempo de aquellas integradas sobre intervalos, utilizando una rejilla temporal escalonada para lograr una integración simpléctica de alto orden.

Resultados Numéricos

Los autores presentan experimentos numéricos en dominios 2D con condiciones de contorno periódicas, probando la advección de ondas senoidales y campanas senoidales en campos de velocidad constantes y de vórtice de Rudman.

• Convergencia: El método demuestra una convergencia algebraica $h$ de orden $(p+1)$ y una convergencia espectral $p$, siempre que el error de integración temporal no domine el error espacial.
• Precisión de la Integración Temporal: Se muestra que el error en la solución depende del orden del esquema de integración temporal ($p_t$). Si el esquema temporal es suficientemente preciso, se preserva el error de discretización inicial; de lo contrario, la dispersión artificial aumenta con el tiempo.
• Conservación de Masa: El error de masa total permanece en el nivel de la precisión de máquina durante los primeros $10^3$ pasos de tiempo y se mantiene por debajo de $10^{-12}$ incluso después de $2 \times 10^4$ pasos, demostrando excelentes propiedades de conservación.
• Reversibilidad: Una prueba que implica la advección de una onda senoidal en un vórtice de Rudman, seguida de la reversión de la dirección del flujo, muestra que el método puede recuperar la solución inicial, probando la reversibilidad del método de integración.
• Dispersión Artificial: El artículo señala que la dispersión artificial es una función del orden de la integración temporal; los esquemas de tiempo de menor orden introducen errores de dispersión que dependen de la frecuencia de la onda advectada.

Significado y Reivindicaciones

El artículo afirma que el esquema derivado combina con éxito la precisión espectral con la conservación de masa local dentro de un marco compatible con la física (mimético). Al extender el marco mimético para incluir la derivada de Lie mediante la fórmula de homotopía de Cartan, los autores proporcionan una discretización que clarifica la estructura geométrica de la ecuación de advección. El método se presenta como un enfoque robusto para resolver problemas de advección donde preservar las propiedades geométricas y topológicas subyacentes de las leyes físicas es esencial. Los autores señalan modestamente que, si bien la discretización espacial es altamente precisa, la precisión general y las características de dispersión dependen del orden del esquema de integración temporal utilizado.