Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation

Resumen Técnico: Singularidades en modelos de separación de fases: un enfoque de elementos espectrales para la ecuación de Cahn–Hilliard no local

1. Planteamiento del Problema

El artículo aborda el desafío numérico de resolver la ecuación de Cahn–Hilliard no local (NLCH), una extensión del modelo clásico para la separación de fases en aleaciones binarias. A diferencia del modelo local, la versión no local incorpora interacciones de largo alcance mediante un núcleo de convolución singular (en este caso, el potencial newtoniano, que es logarítmico en 2D).

El sistema se describe mediante:
$$\frac{d\rho}{dt} = \nabla \cdot \left( m(\rho) \nabla \mu \right), \quad \mu = F'(\rho) - K * \rho$$
donde:

• $\rho$ es el parámetro de orden (fracción de masa).
• $m(\rho)$ es la movilidad (constante en este trabajo).
• $F(\rho)$ es un potencial logarítmico (singular), que asegura que $\rho$ permanezca en el intervalo físico $[-1, 1]$.
• $K$ es un núcleo de interacción singular en el origen (tipo Newtoniano).

Desafíos principales:

1. Singularidad del núcleo: La convolución $K * \rho$ involucra una singularidad integrable en el origen, lo que dificulta la discretización directa en esquemas de colocación estándar.
2. Potencial singular: El potencial logarítmico $F$ tiene derivadas que tienden a infinito en los estados puros ($\pm 1$), lo que genera rigidez numérica.
3. Operadores densos: La naturaleza no local de la convolución genera matrices densas, aumentando el costo computacional.
4. Dominios complejos: La teoría de existencia suele requerir dominios lisos, pero las simulaciones numéricas a menudo utilizan mallas poligonales o poliédricas.

2. Metodología Propuesta

Los autores desarrollan un marco numérico eficiente basado en un método de elementos espectrales (SEM) combinado con una técnica de descomposición para manejar la singularidad.

A. Análisis Teórico y Marco Matemático:

• Se extiende la teoría de existencia y unicidad de soluciones débiles a dominios de Lipschitz (que incluyen dominios con esquinas, típicos en mallas computacionales), bajo la condición de cono.
• Se demuestra la propiedad de separación estricta instantánea: para $t > 0$, la solución $\rho$ se mantiene uniformemente alejada de los estados puros $\pm 1$, lo que garantiza la regularidad de la solución.
• Se establece la convergencia de las soluciones hacia un único equilibrio con una tasa de decaimiento polinomial, utilizando la desigualdad de Łojasiewicz–Simon.

B. Estrategia Numérica: Descomposición de Corte (Cutoff Decomposition):
Para aproximar la integral de convolución $(K * \rho)(x)$, se divide el dominio de integración en dos regiones:

1. Región externa ($\Omega \setminus B_\epsilon(x)$): Donde el núcleo es suave. Se aproxima utilizando un método de elementos espectrales sobre una "multiforma" (multishape).
2. Región interna ($B_\epsilon(x)$): Donde ocurre la singularidad. Se utiliza una aproximación de Taylor de orden cero (asumiendo que $\rho$ es constante en la vecindad pequeña), lo que permite calcular la integral del núcleo singular analíticamente.

La aproximación total es:
$$(K * \rho)(x) \approx \int_{\Omega \setminus B_\epsilon} K(x-y)\rho(y)dy + \rho(x) \int_{B_\epsilon \cap \Omega} K(u)du$$

C. Implementación con Elementos Espectrales (Multishape):

• Se utiliza el paquete de software MultiShape (basado en ChebClass).
• Para cada punto de colocación $x_i$, se define una caja flotante alrededor de la singularidad. El dominio restante se descompone en un número mínimo de cuadriláteros (partición máxima) para formar una "multiforma".
• Se precalcula una matriz de convolución para la parte no singular, aprovechando que el núcleo es constante en el tiempo.
• Se emplea un solucionador de ecuaciones diferenciales-algebraicas (DAE) rígido (ode15s en MATLAB) para la evolución temporal, manejando las condiciones de frontera de flujo nulo.

3. Contribuciones Clave

1. Extensión Teórica: Demostración de la existencia, unicidad y convergencia a equilibrio para el sistema NLCH con potencial singular en dominios de Lipschitz (incluyendo 3D), un resultado no disponible previamente en la literatura para este contexto.
2. Nuevo Esquema Numérico: Desarrollo de un método de elementos espectrales específico para convoluciones con núcleos singulares, utilizando una descomposición de corte y aproximación de Taylor de orden cero. Este enfoque es flexible para diferentes núcleos y dominios.
3. Estimación de Error Rigurosa: Se proporcionan estimaciones de error analíticas que muestran una convergencia de orden $O(\epsilon^2 \log \epsilon)$, verificadas computacionalmente.
4. Código Abierto: Se libera el código fuente para reproducir los experimentos y generar soluciones de alta resolución.

4. Resultados Numéricos

Los autores realizaron pruebas exhaustivas en dominios 2D (cuadrado unitario) y 3D (caja unitaria), utilizando el potencial logarítmico y el núcleo newtoniano.

• Dinámica de Separación de Fases:
  • Se observó que el parámetro de escalado del núcleo ($\eta$) controla la competencia entre difusión y atracción no local.
  • Para $\eta > 0$ (dominio difusivo), el sistema converge a un estado homogéneo.
  • Para $\eta < 0$ (dominio atractivo), se produce una separación de fases rápida, formando estructuras estacionarias no triviales (similares a funciones escalón).
• Validación de la Aproximación:
  • Los errores numéricos para funciones constantes y soluciones de prueba se mantuvieron por debajo de los umbrales de precisión de máquina, confirmando la eficacia de la descomposición de corte.
  • Se demostró que la precomputación de la matriz de convolución hace que el método sea computacionalmente eficiente.
• Potencial Regularizado:
  • Se introdujo una versión regularizada del potencial logarítmico ($F_\omega$) para acelerar las simulaciones de largo plazo.
  • Los resultados mostraron que la solución regularizada converge a la solución original con alta precisión (error $\sim 10^{-14}$) pero con un costo computacional significativamente menor, ya que elimina la rigidez asociada a las singularidades en los estados puros.
• Núcleos Mixtos:
  • Se exploró una combinación lineal del núcleo newtoniano con un suavizador (mollifier), lo que permitió generar dinámicas ricas con estructuras de fase separada más complejas y patrones distintivos.

5. Significado e Impacto

Este trabajo es pionero en proporcionar soluciones numéricas de alta resolución para el sistema de Cahn–Hilliard no local con potencial logarítmico y núcleo singular, logrando un equilibrio entre precisión y costo computacional limitado.

• Avance Metodológico: La técnica de "multiforma" combinada con la descomposición de corte permite tratar singularidades en dominios complejos sin perder la alta precisión espectral.
• Aplicabilidad: El marco desarrollado es generalizable a otros sistemas de dinámica de fluidos y modelos de densidad funcional (DDFT).
• Futuro: Los autores planean utilizar este enfoque para resolver sistemas acoplados, como el modelo de crecimiento tumoral de Hawkins-Daruud o el sistema Cahn–Hilliard-Navier-Stokes, para aplicaciones en biología y ciencia de materiales.

En resumen, el artículo presenta una solución robusta y eficiente a un problema numérico difícil, combinando análisis matemático riguroso con técnicas computacionales avanzadas para modelar fenómenos de separación de fases con interacciones no locales.