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

Riepilogo Tecnico: Ridistribuzione dello Stato basata sull'Entità dell'Aggiornamento (UM-SRD)

Enunciato del Problema
I metodi a celle tagliate permettono agli schemi a volumi finiti di gestire geometrie complesse senza mesh adattate al corpo, ma introducono sfide di stabilità severe a causa di celle tagliate arbitrariamente piccole. I metodi di Ridistribuzione dello Stato (SRD), in particolare il framework SRD pesato sviluppato da Berger & Giuliani (2024), affrontano questi problemi di stabilità fondendo le celle piccole in vicinati, evolvendo le medie sulla mesh fusa e ridistribuendo i risultati sulla mesh originale.

Tuttavia, una limitazione critica della formulazione SRD pesata esistente è la mancanza di un meccanismo di "spegnimento". A differenza dei metodi classici di ridistribuzione del flusso, dove il post-processing svanisce naturalmente quando gli squilibri di flusso sono nulli, la SRD pesata applica un operatore di ridistribuzione non banale ad ogni passo temporale, indipendentemente dall'entità dell'aggiornamento a volumi finiti. Questa attivazione continua porta a due carenze pratiche:

1. Mancata conservazione degli stati stazionari: Anche quando il bilancio di flusso è esattamente zero, la SRD sostituisce i valori puntuali con medie pesate per volume, impedendo la convergenza alla precisione di macchina e introducendo errori $O(1)$ nelle regioni statiche.
2. Dissipazione non necessaria: Nelle regioni lisce con aggiornamenti piccoli ma non nulli, l'applicazione continua della ridistribuzione introduce diffusione numerica proporzionale al gradiente del vicinato, indipendentemente dalla dinamica fisica effettiva.

Metodologia
Questo articolo introduce la Ridistribuzione dello Stato basata sull'Entità dell'Aggiornamento (UM-SRD), un'estensione di spegnimento della SRD pesata progettata per fondere l'operatore SRD con l'operatore identità basandosi sull'entità locale dell'aggiornamento a volumi finiti.

Il metodo opera all'interno del framework standard di integrazione temporale a tre stadi:

1. Aggiornamento a Volumi Finiti: Calcolare lo stato pre-ridistribuzione $U^*$ utilizzando lo schema di base.
2. Aggregazione del Vicinato: Calcolare le medie pesate del vicinato $\hat{Q}^*_j$ sui vicinati fusi $M_j$.
3. Ridistribuzione Fusa: Invece di applicare direttamente l'operatore SRD $S_j$, la UM-SRD applica una combinazione convessa:
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   dove $I$ è l'operatore identità e $s_j \in [0, 1]$ è un parametro di fusione.

Il parametro di fusione $s_j$ è determinato da una funzione di spegnimento liscia e localmente definita $h(\eta_j)$:
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
Qui, $\eta_j$ è un indicatore normalizzato dell'entità massima dell'aggiornamento assoluto all'interno del vicinato $M_j$ (specificamente $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$).

• Regime a Aggiornamento Zero: Se tutti gli aggiornamenti nel vicinato sono alla precisione di macchina ($\Delta U_i \approx 0$), allora $\eta_j = 0$, $s_j = 0$, e il metodo si riduce esattamente all'identità ($U^{n+1} = U^*$).
• Regime Attivo: Se si verificano aggiornamenti significativi (ad esempio, vicino a discontinuità o quando la rigidità delle celle piccole è attiva), $\eta_j \approx 1$, $s_j \approx 1$, e il metodo si comporta come una SRD pesata completa.

Contributi Analitici Chiave
L'articolo fornisce un'analisi rigorosa della UM-SRD per un modello di advezione lineare unidimensionale con una singola cella tagliata piccola:

• Conservazione: L'operatore fuso è dimostrato essere conservativo per costruzione per qualsiasi $s_j \in [0, 1]$.
• Diminuzione della Variazione Totale (TVD): Sotto la stessa condizione CFL dello schema upwind del primo ordine di base, la UM-SRD è dimostrata essere TVD sulla mesh originale a celle tagliate. Ciò si basa sulla combinazione convessa dell'identità e dell'operatore SRD TVD.
• Conservazione dello Stato Stazionario: È dimostrato che la UM-SRD preserva esattamente gli stati stazionari dello schema di base. Se l'aggiornamento a volumi finiti è zero, $s_j=0$, e non avviene alcuna ridistribuzione, evitando la deriva spuria intrinseca nella SRD standard.
• Accuratezza: Nelle regioni lisce, la UM-SRD è mostrata essere una perturbazione dello schema di base. Se il parametro di fusione scala appropriatamente con la dimensione della mesh (o tramite l'indicatore normalizzato usato nella pratica), il contributo aggiuntivo all'errore di troncamento locale è di ordine superiore rispetto allo schema del primo ordine di base. Di conseguenza, la UM-SRD mantiene il tasso di convergenza del primo ordine del metodo upwind sottostante.

Risultati Numerici
Esperimenti numerici su problemi modello unidimensionali e bidimensionali confermano le previsioni teoriche:

• Convergenza: La UM-SRD esibisce convergenza del primo ordine nelle norme $L_1$ e $L_\infty$ per l'advezione liscia, con errori indistinguibili dallo schema di base e dalla SRD pesata standard (differenze < 0,01%).
• Comportamento di Spegnimento: Nei test con impulsi compatti, il parametro di fusione $s_j$ scende esattamente a zero quando l'impulso è assente, recuperando lo schema di base, e sale a $\approx 0,99$ quando l'impulso interagisce con la cella piccola.
• Stabilità: La UM-SRD stabilizza lo schema di base sotto passi temporali a mesh completa (dove il CFL locale alla cella piccola è 10 volte il limite di stabilità), prevenendo la divergenza esponenziale osservata nello schema upwind di base.
• Proprietà Non Oscillatorie: Test con dati discontinui mostrano che la UM-SRD rimane non oscillatoria, corrispondendo al comportamento della SRD standard vicino alla cella piccola.
• Conservazione dello Stato Stazionario: Esperimenti confermano che mentre la SRD standard introduce errori $O(1)$ nei profili statici, la UM-SRD mantiene il profilo iniziale alla precisione di macchina.
• Estensione 2D: L'algoritmo è stato applicato con successo a una griglia cartesiana 2D con un confine immerso, dimostrando convergenza del primo ordine senza modifiche alla logica di fusione.

Significato
L'articolo afferma che la UM-SRD offre una modifica pratica e minima all'infrastruttura esistente della SRD pesata. Riutilizzando le strutture di vicinato e pesi consolidate e aggiungendo solo un meccanismo di fusione locale, essa dota la SRD di una proprietà di aggiornamento che svanisce. Ciò combina i benefici di stabilità della ridistribuzione dello stato per celle tagliate arbitrariamente piccole con le proprietà di conservazione dell'accuratezza dei metodi di ridistribuzione del flusso, abilitando specificamente la conservazione esatta degli stati stazionari e riducendo la dissipazione non necessaria nelle regioni lisce. Gli autori notano che, sebbene l'analisi corrente sia limitata all'advezione lineare scalare 1D, la struttura algoritmica è indipendente dalla dimensione e pronta per l'integrazione in codici di produzione a celle tagliate.