Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

Resumen Técnico: Redistribución de Estado por Magnitud de Actualización (UM-SRD)

Enunciado del Problema
Los métodos de celdas cortadas permiten que los esquemas de volumen finito manejen geometrías complejas sin mallas adaptadas al cuerpo, pero introducen desafíos severos de estabilidad debido a celdas cortadas arbitrariamente pequeñas. Los métodos de Redistribución de Estado (SRD), específicamente el marco SRD ponderado desarrollado por Berger y Giuliani (2024), abordan estos problemas de estabilidad fusionando celdas pequeñas en vecindades, evolucionando promedios en la malla fusionada y redistribuyendo los resultados de vuelta a la malla original.

Sin embargo, una limitación crítica de la formulación SRD ponderada existente es la falta de un mecanismo de "desactivación". A diferencia de los métodos clásicos de redistribución de flujos, donde el posprocesamiento desaparece naturalmente cuando los desequilibrios de flujo son cero, el SRD ponderado aplica un operador de redistribución no trivial en cada paso de tiempo, independientemente de la magnitud de la actualización del volumen finito. Esta activación continua conduce a dos deficiencias prácticas:

1. Fallo en preservar estados estacionarios: Incluso cuando el balance de flujo es exactamente cero, el SRD reemplaza los valores puntuales con promedios ponderados por volumen, impidiendo la convergencia a precisión de máquina e introduciendo errores de $O(1)$ en regiones estáticas.
2. Disipación innecesaria: En regiones suaves con actualizaciones pequeñas pero no nulas, la aplicación continua de la redistribución introduce difusión numérica proporcional al gradiente de la vecindad, independientemente de la dinámica física real.

Metodología
Este artículo introduce la Redistribución de Estado por Magnitud de Actualización (UM-SRD), una extensión de desactivación del SRD ponderado diseñada para mezclar el operador SRD con el operador identidad basándose en la magnitud local de la actualización del volumen finito.

El método opera dentro del marco estándar de pasos de tiempo de tres etapas:

1. Actualización de Volumen Finito: Calcular el estado pre-redistribución $U^*$ utilizando el esquema base.
2. Agregación de Vecindad: Calcular promedios ponderados de vecindad $\hat{Q}^*_j$ en vecindades fusionadas $M_j$.
3. Redistribución Mezclada: En lugar de aplicar el operador SRD $S_j$ directamente, el UM-SRD aplica una combinación convexa:
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   donde $I$ es el operador identidad y $s_j \in [0, 1]$ es un parámetro de mezcla.

El parámetro de mezcla $s_j$ se determina mediante una función de desactivación suave y localmente definida $h(\eta_j)$:
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
Aquí, $\eta_j$ es un indicador normalizado de la magnitud máxima de actualización absoluta dentro de la vecindad $M_j$ (específicamente $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$).

• Régimen de Actualización Cero: Si todas las actualizaciones en la vecindad están a precisión de máquina ($\Delta U_i \approx 0$), entonces $\eta_j = 0$, $s_j = 0$, y el método se reduce exactamente a la identidad ($U^{n+1} = U^*$).
• Régimen Activo: Si ocurren actualizaciones significativas (por ejemplo, cerca de discontinuidades o cuando la rigidez de celdas pequeñas está activa), $\eta_j \approx 1$, $s_j \approx 1$, y el método se comporta como un SRD ponderado completo.

Contribuciones Analíticas Clave
El artículo proporciona un análisis riguroso del UM-SRD para un modelo de advección lineal unidimensional con una única celda cortada pequeña:

• Conservación: Se demuestra que el operador mezclado es conservativo por construcción para cualquier $s_j \in [0, 1]$.
• Disminución de la Variación Total (TVD): Bajo la misma condición CFL que el esquema de viento ascendente de primer orden base, se demuestra que el UM-SRD es TVD en la malla de celdas cortadas original. Esto depende de la combinación convexa de la identidad y el operador SRD TVD.
• Preservación de Estados Estacionarios: Se demuestra que el UM-SRD preserva exactamente los estados estacionarios del esquema base. Si la actualización del volumen finito es cero, $s_j=0$, y no ocurre redistribución, evitando la deriva espuria inherente al SRD estándar.
• Precisión: En regiones suaves, se muestra que el UM-SRD es una perturbación del esquema base. Si el parámetro de mezcla escala apropiadamente con el tamaño de la malla (o mediante el indicador normalizado utilizado en la práctica), la contribución adicional al error de truncamiento local es de orden superior al esquema de primer orden base. En consecuencia, el UM-SRD conserva la tasa de convergencia de primer orden del método de viento ascendente subyacente.

Resultados Numéricos
Los experimentos numéricos en problemas modelo unidimensionales y bidimensionales confirman las predicciones teóricas:

• Convergencia: El UM-SRD exhibe convergencia de primer orden en las normas $L_1$ y $L_\infty$ para advección suave, con errores indistinguibles del esquema base y del SRD ponderado estándar (diferencias < 0.01%).
• Comportamiento de Desactivación: En pruebas con pulsos compactos, el parámetro de mezcla $s_j$ cae exactamente a cero cuando el pulso está ausente, recuperando el esquema base, y sube a $\approx 0.99$ cuando el pulso interactúa con la celda pequeña.
• Estabilidad: El UM-SRD estabiliza el esquema base bajo pasos de tiempo de malla completa (donde el CFL local en la celda pequeña es 10 veces el límite de estabilidad), previniendo la divergencia exponencial observada en el esquema de viento ascendente base.
• Propiedades No Oscilatorias: Las pruebas con datos discontinuos muestran que el UM-SRD permanece no oscilatorio, coincidiendo con el comportamiento del SRD estándar cerca de la celda pequeña.
• Preservación de Estados Estacionarios: Los experimentos confirman que, mientras que el SRD estándar introduce errores de $O(1)$ en perfiles estáticos, el UM-SRD mantiene el perfil inicial con precisión de máquina.
• Extensión 2D: El algoritmo se aplicó con éxito a una malla cartesiana 2D con una frontera embebida, demostrando convergencia de primer orden sin modificación a la lógica de mezcla.

Significado
El artículo afirma que el UM-SRD ofrece una modificación práctica y mínima a la infraestructura existente de SRD ponderado. Al reutilizar las estructuras de vecindad y pesos establecidas y agregar solo un mecanismo de mezcla local, dota al SRD de una propiedad de actualización que desaparece. Esto combina los beneficios de estabilidad de la redistribución de estado para celdas cortadas arbitrariamente pequeñas con las propiedades de preservación de precisión de los métodos de redistribución de flujos, permitiendo específicamente la preservación exacta de estados estacionarios y reduciendo la disipación innecesaria en regiones suaves. Los autores señalan que, aunque el análisis actual se restringe a la advección lineal escalar unidimensional, la estructura algorítmica es independiente de la dimensión y está lista para su integración en códigos de celdas cortadas de producción.