A domain decomposition strategy for natural imposition of mixed boundary conditions in port-Hamiltonian systems

Sintesi Tecnica: Una strategia di decomposizione del dominio per l'imposizione naturale di condizioni al contorno miste in sistemi port-Hamiltoniani

Enunciato del Problema
L'articolo affronta la sfida di imporre condizioni al contorno miste (condizioni simultanee di Dirichlet e Neumann) in sistemi iperbolici port-Hamiltoniani (pH) utilizzando il metodo degli elementi finiti. Nelle formulazioni standard degli elementi finiti, le condizioni al contorno vengono tipicamente imposte o fortemente (limitando gli spazi di funzioni discrete) o debolmente (tramite termini di contorno naturali nella formulazione variazionale). L'imposizione forte porta spesso a equazioni differenziali-algebriche (DAE), che sono computazionalmente più difficili da risolvere rispetto alle equazioni differenziali ordinarie (ODE). Al contrario, l'imposizione debole richiede solitamente l'introduzione di moltiplicatori di Lagrange per imporre le condizioni al contorno essenziali, aumentando la dimensione e la complessità del sistema. Gli autori osservano che, sebbene il framework port-Hamiltoniano gestisca naturalmente i flussi di confine tramite la struttura Stokes-Dirac, gli schemi di discretizzazione esistenti spesso faticano a incorporare esplicitamente le condizioni miste senza ricorrere ai moltiplicatori di Lagrange, pur preservando la sottostante struttura geometrica.

Metodologia
La strategia proposta impiega un approccio di decomposizione del dominio combinato con la Finite Element Exterior Calculus (FEEC) e una formulazione primale-duale. La metodologia centrale prevede i seguenti passaggi:

1. Decomposizione del Dominio: Il dominio spaziale $\Omega$ viene suddiviso in due sottodomini non sovrapposti, $\Omega_1$ e $\Omega_2$, separati da un'interfaccia interna arbitraria $\Gamma_{int}$.
2. Rappresentazione del Campo Duale: Vengono applicate simultaneamente due distinte formulazioni di elementi finiti misti:
   • Sistema 1 (su $\Omega_1$): Utilizza una formulazione in cui la condizione al contorno su $\Gamma_1$ è imposta naturalmente (debolmente), e la condizione sull'interfaccia è trattata come un input. Questa formulazione utilizza spazi di elementi finiti derivati da un sottocomplesso di Hilbert in cui lo spazio della variabile duale è un sottoinsieme dell'immagine dello spazio primale sotto l'operatore aggiunto.
   • Sistema 2 (su $\Omega_2$): Utilizza la formulazione duale in cui la condizione al contorno su $\Gamma_2$ è imposta naturalmente, e la condizione sull'interfaccia è trattata come un input. Gli spazi degli elementi finiti sono selezionati da un sottocomplesso di Hilbert aggiunto.
3. Interconnessione di Feedback: I due sottodomini sono accoppiati all'interfaccia $\Gamma_{int}$ tramite un'interconnessione a giratore (un tipo specifico di feedback nei sistemi port-Hamiltoniani). Questa interconnessione impone la continuità degli spazi degli elementi finiti in modo debole. Nello specifico, l'input di un sottodominio è impostato uguale all'output (variabile coniugata) dell'altro, e viceversa.
4. Discretizzazione: Il metodo utilizza sottocomplessi di elementi finiti (ad esempio, complessi de Rham) per garantire stabilità e preservazione della struttura. Ad esempio, nell'equazione delle onde, $\Omega_1$ potrebbe utilizzare elementi Discontinuous Galerkin (DG) e Raviart-Thomas (RT), mentre $\Omega_2$ utilizza elementi Continuous Galerkin (CG) e Nédélec.
5. Integrazione Temporale: L'articolo valuta due schemi di integrazione temporale:
   • Midpoint Implicito: Preserva l'energia complessiva e la struttura di Poisson del sistema, ma richiede la risoluzione di un sistema monolitico.
   • Störmer-Verlet: Permette il disaccoppiamento dei due sottodomini (risolvendoli separatamente), ma non preserva strettamente l'energia totale, sebbene mantenga il bilancio di potenza all'interno di ciascun sottodominio.

Contributi Chiave
L'articolo formalizza e dimostra i seguenti contributi:

• Condizioni al Contorno Miste Senza Lagrange: La strategia impone con successo condizioni al contorno miste esplicitamente senza introdurre moltiplicatori di Lagrange, convertendo il problema in un sistema di ODE anziché DAE.
• Decomposizione del Dominio Primale-Duale: Diversamente dai classici metodi Schwarz alternati Dirichlet-Neumann che si affidano a un accoppiamento iterativo, questo metodo utilizza una formulazione primale-duale simultanea per accoppiare i sottodomini tramite un'interconnessione di feedback diretta.
• Preservazione della Struttura: L'uso di sottocomplessi di Hilbert assicura che la discretizzazione preservi la struttura geometrica sottostante (struttura Stokes-Dirac) e le proprietà fisiche (ad esempio, il bilancio di potenza) del sistema port-Hamiltoniano.
• Robustezza al Shear Locking: Si dimostra che la formulazione mista evita i fenomeni di shear locking nei problemi di meccanica strutturale, un problema comune nelle discretizzazioni standard degli elementi finiti per piastre e travi.
• Estendibilità: Viene dimostrato che la metodologia è applicabile a sistemi semilineari contenenti non linearità algebriche (dove la non linearità non interferisce con gli operatori differenziali) ed è teoricamente estendibile a problemi ellittici statici e fenomeni elettromagnetici.

Risultati Numerici
Gli autori validano il metodo attraverso quattro esempi numerici:

1. Trave Intrinseca Geometricamente Esatta (1D): Un modello semilineare che dimostra la capacità del metodo di gestire non linearità algebriche ed evitare lo shear locking in uno scenario di roll-up quasi-statico.
2. Equazione delle Onde 2D: Uno studio di convergenza conferma i tassi di convergenza teorici ($h^k$). L'analisi spettrale mostra che il metodo proposto produce frequenze proprie con accuratezza comparabile o superiore alle classiche discretizzazioni di elementi Lagrange lineari. Il bilancio di potenza è preservato fino alla precisione di macchina all'interno di ciascun sottodominio.
3. Elastodinamica Lineare: L'analisi modale su un quadrato unitario dimostra che la nostra discretizzazione mista raggiunge un'accuratezza comparabile agli elementi Lagrange quadratici classici.
4. Piastra di Mindlin: Un problema di benchmark utilizzato per testare lo shear locking. I risultati mostrano che la proposta discretizzazione mista supera significativamente le discretizzazioni primali classiche in termini di accuratezza per le frequenze proprie, eliminando efficacemente lo shear locking anche con mesh grossolane.

Significato e Rivendicazioni
L'articolo sostiene che questa strategia offre un'alternativa robusta per la modellazione di sistemi multifisici complessi dove la descrizione modulare e la preservazione della struttura sono critiche. Evitando i moltiplicatori di Lagrange, il metodo semplifica la struttura algebrica delle equazioni risultanti, il che è benefico non solo per la simulazione ma anche per l'ottimizzazione numerica dove i vincoli algebrici possono essere problematici. Gli autori pongono questo lavoro come una "prova di concetto" dimostrando che le formulazioni miste primale-duali possono essere utilizzate simultaneamente per gestire naturalmente le condizioni al contorno miste. Suggeriscono che, sebbene l'attuale implementazione richieda la definizione di un'interfaccia, l'approccio potrebbe essere esteso a metodi Discontinuous Galerkin completamente sviluppati dove l'interconnessione avviene tra ogni elemento finito, riducendo potenzialmente il costo computazionale della risoluzione di grandi sistemi monolitici. Il lavoro è progettato principalmente per sistemi port-Hamiltoniani iperbolici, ma nota esplicitamente il suo potenziale di estensione a problemi statici e semilineari.