Unfitted Multi-Level hp Refinement for Localized and Moving Solution Features

1. Enunciato del Problema

L'analisi agli elementi finiti accurata di problemi con caratteristiche localizzate (singolarità, gradienti ripidi, discontinuità) e sorgenti in movimento (ad esempio, nella manifattura additiva o nella propagazione di cricche) è computazionalmente impegnativa.

• Limitazioni dei Metodi Convenzionali: Le strategie adattive standard (rifinimento $h$ conforme) richiedono la generazione ripetuta di mesh, la modifica della topologia e la proiezione di variabili di stato. Ciò introduce un significativo sovraccarico algoritmico e limita la robustezza nelle simulazioni transitorie.
• Limitazioni dei Metodi di Sovrapposizione Esistenti: Sebbene il "rifinimento per sovrapposizione" (sovrapponendo mesh) eviti la rimeshatura, i metodi $hp$ multilivello allineati (fitted) esistenti impongono vincoli topologici rigorosi. Richiedono che i confini delle sovrapposizioni si allineino perfettamente con i confini degli elementi sottostanti e spesso necessitano di complesse regole di disattivazione per garantire l'indipendenza lineare e la continuità. Ciò limita la flessibilità nel posizionamento delle zone di rifinimento e aumenta l'accoppiamento tra i livelli.

2. Metodologia: Rifinimento $hp$ Multilivello Non Allineato (Unfitted)

Gli autori propongono un nuovo framework che estende il rifinimento per sovrapposizione a mesh non allineate e non corrispondenti.

Concetto Chiave

• Strategia di Sovrapposizione: Lo spazio di approssimazione globale ($V_h$) è costruito mediante la sovrapposizione di una discretizzazione di base grossolana e multiple mesh di sovrapposizione posizionate indipendentemente.
  $$u_h(\mathbf{x}) = \sum_{k=0}^{L} u_h^{(k)}(\mathbf{x})$$
• Natura Non Allineata: A differenza dei metodi allineati, i confini delle mesh di sovrapposizione non devono allinearsi con i confini degli elementi della base o di altri livelli di sovrapposizione. Le mesh di sovrapposizione possono coprire parzialmente gli elementi di base.
• Imposizione della Continuità: La continuità globale $C^0$ è imposta applicando vincoli di Dirichlet omogenei sui confini artificiali dei domini di sovrapposizione ($\Gamma_o^{(k)}$). Ciò garantisce che il contributo della sovrapposizione si annulli al bordo della regione di rifinimento, prevenendo discontinuità nella soluzione globale.
• Indipendenza Lineare: Il metodo rilassa le rigide regole di disattivazione richieste negli schemi gerarchici allineati. Sebbene dipendenze patologiche siano teoricamente possibili, gli esperimenti numerici non hanno mostrato problemi, suggerendo che i vincoli per la continuità siano sufficienti per le applicazioni pratiche.

Aspetti Implementativi

• Regioni di Integrazione: Per gestire mesh non corrispondenti, l'intersezione delle partizioni degli elementi è suddivisa in regioni di integrazione $C^\infty$ (sotto-celle). Un algoritmo indipendente dalla dimensione raccoglie le coordinate uniche dei vertici da tutte le mesh sovrapposte per definire queste sotto-regioni.
• Elementi Ausiliari: Viene introdotta una "sottile strato di sovrapposizione". Per ogni regione di integrazione, un elemento ausiliario aggrega le funzioni di forma da tutti gli elementi genitori contribuenti attraverso diversi livelli di mesh. Ciò consente il riutilizzo delle routine standard agli elementi finiti (interpolazione, valutazione del gradiente, assemblaggio della matrice) senza modificare il solver principale.
• Caratteristiche in Movimento: Per problemi transitori, le mesh di sovrapposizione possono essere inserite, traslate o rimosse dinamicamente. Il trasferimento di stato tra i passi temporali viene eseguito tramite proiezione variazionale (ad esempio, proiezione $L^2$) piuttosto che interpolazione puntuale per mantenere la coerenza.

3. Contributi Chiave

1. Disaccoppiamento della Topologia: Il metodo rimuove il requisito che i confini delle sovrapposizioni corrispondano ai confini degli elementi di base, permettendo un posizionamento e una dimensionamento arbitrari delle zone di rifinimento.
2. Ridotta Complessità: Evitando il tracciamento topologico complesso e le regole di disattivazione intrinseche ai metodi $hp$ multilivello allineati, l'implementazione è più leggera e flessibile.
3. Controllo dei Parametri: L'introduzione di un parametro di dimensione della mesh ($\alpha$) per i domini di sovrapposizione fornisce un ulteriore grado di libertà per ottimizzare la qualità dell'approssimazione e i tassi di convergenza.
4. Adattabilità Dinamica: Il framework supporta caratteristiche in movimento senza rimeshatura globale, rendendolo altamente adatto per simulazioni transitorie come la manifattura additiva.

4. Risultati Numerici

Il metodo è stato validato su problemi di riferimento con soluzioni in forma chiusa e su un problema pratico di sorgente termica in movimento.

• Barra Elastica 1D con Deformazioni Discontinue:
  • Risultato: Il metodo non allineato ha raggiunto una convergenza esponenziale.
  • Efficienza: Per un errore target di $\approx 4.34 \times 10^{-5}\%$, il metodo non allineato ha richiesto 211 incognite, mentre il metodo allineato ha richiesto 436 incognite (errore $2.17 \times 10^{-3}\%$). L'approccio non allineato ha prodotto errori significativamente inferiori per grado di libertà.
• Problema dell'Angolo Singolare 2D:
  • Risultato: Il metodo non allineato ha evitato la creazione di elementi non necessari nelle regioni lontane dalla singolarità (a differenza del metodo allineato che biseca interi elementi).
  • Efficienza: L'approccio non allineato ha raggiunto un errore di $1.54783\%$ con 1.177 incognite, rispetto al $1.54778\%$ con 3.013 incognite per il metodo allineato. Ulteriore ottimizzazione della dimensione della sovrapposizione ($\alpha=0.2$) ha ridotto l'errore allo $0.127\%$ con sole 1.160 incognite.
• Condizionamento della Piccola Sovrapposizione:
  • Lo studio ha analizzato l'effetto di sovrapposizioni molto piccole tra elementi sovrapposti. Sebbene i numeri di condizione aumentino al diminuire della sovrapposizione (specialmente per alti gradi polinomiali $p \geq 6$), le simulazioni con sovrapposizione in movimento non hanno incontrato configurazioni in cui ciò dominasse il costo computazionale.
• Sorgente Termica in Movimento (Simulazione di Manifattura Additiva):
  • Configurazione: Una sorgente termica si è spostata lungo un percorso circolare.
  • Confronto: Una mesh non allineata localmente rifinita (1.841 incognite) ha corrisposto a una soluzione di riferimento ad alta risoluzione (250.000 incognite) con alta accuratezza.
  • Accuratezza: Il modello rifinito ha catturato temperature di picco e gradienti con deviazioni relative <1%, mentre una mesh di ordine superiore non rifinita ha sottostimato i picchi del 32-71%. Il modello rifinito ha utilizzato <1% dei gradi di libertà della mesh di riferimento.

5. Significato e Conclusione

Questo lavoro presenta un'alternativa robusta ed efficiente alla meshatura adattiva tradizionale per problemi con caratteristiche localizzate o in movimento.

• Prestazioni: Mantiene i tassi di convergenza esponenziale dei metodi $hp$ riducendo drasticamente il numero di gradi di libertà richiesti per alta accuratezza.
• Flessibilità: La natura "non allineata" consente di spostare e ridimensionare dinamicamente le zone di rifinimento senza rimeshatura o aggiornamenti topologici complessi.
• Applicabilità: Il metodo è particolarmente significativo per le simulazioni transitorie (ad esempio, manifattura additiva, propagazione di cricche) dove il costo della rimeshatura ripetuta e del trasferimento di stato è proibitivo.
• Lavori Futuri: Gli autori suggeriscono ricerche future su indicatori di rifinimento automatizzati, analisi sistematica della sensibilità per piccole sovrapposizioni e tecniche di precondizionamento avanzate per sistemi mal condizionati derivanti da intersezioni minime.

In sintesi, il metodo $hp$ multilivello non allineato offre un rapporto errore-costi superiore rispetto agli approcci multilivello allineati, abilitando simulazioni ad alta fedeltà di fenomeni fisici complessi ed evolutivi con un sovraccarico computazionale minimo.