A $C^0$ weak Galerkin method with preconditioning for constrained optimal control problems with general tracking

Sintesi Tecnica: Metodo Weak Galerkin $C^0$ con Precondizionamento per Problemi di Controllo Ottimo Vincolati

1. Il Problema

Il lavoro si concentra sulla risoluzione numerica di Problemi di Controllo Ottimo (OCP) governati da equazioni alle derivate parziali (PDE) ellittiche. Nello specifico, il problema affrontato presenta due caratteristiche critiche che ne aumentano la complessità analitica e numerica:

• Funzionale di costo di "tracking" generale: L'obiettivo è far coincidere lo stato del sistema con un target desiderato ($y_d$) non solo in tutto il dominio, ma anche su punti specifici, curve o sottodomini. Questo riduce la regolarità delle soluzioni e dei moltiplicatori di Lagrange associati.
• Vincoli puntuali sullo stato: Lo stato $y(x)$ deve rimanere all'interno di limiti inferiori e superiori ($\psi^-(x) \le y(x) \le \psi^+(x)$) quasi ovunque nel dominio.

La presenza di questi vincoli trasforma il problema in una disuguaglianza variazionale del quarto ordine per la variabile di stato. La regolarità ridotta della soluzione ottima e la natura ill-conditioned (mal condizionata) dei sistemi lineari risultanti rendono difficile ottenere ordini di convergenza ottimali e garantire la stabilità numerica.

2. Metodologia

Gli autori propongono un approccio ibrido che combina una nuova discretizzazione spaziale e una strategia di precondizionamento:

• Metodo Weak Galerkin $C^0$ (C0-WG):

  • Viene utilizzato un elemento finito di Lagrange quadratico globale continuo ($P_2$).
  • A differenza dei metodi tradizionali che richiedono elementi $C^1$ (complessi da implementare) o metodi di penalità interna simmetrica ($C^0$-IP) che necessitano di parametri di penalità e termini di traccia complessi, il metodo C0-WG definisce un Laplaciano debole discreto.
  • Il Laplaciano debole è calcolato localmente come una funzione costante a tratti all'interno di ogni elemento, utilizzando i valori delle funzioni e dei gradienti sui barycentri delle interfacce condivise.
  • La forma bilineare risultante è composta dal prodotto scalare $L^2$ dei Laplaciani deboli e da un termine di stabilizzazione privo di parametri (parameter-free), eliminando la necessità di calibrare parametri di penalità.

• Precondizionatore di Schwarz Additivo:

  • Per affrontare l'ill-conditioning dei sistemi lineari derivanti dalle PDE del quarto ordine (specialmente con termini di tracking generale), viene sviluppato un precondizionatore di Schwarz additivo a un livello.
  • Questa tecnica di decomposizione del dominio suddivide il problema globale in sottoproblemi sovrapposti risolti indipendentemente, combinando poi le soluzioni additivamente. Questo approccio è altamente parallelo e migliora drasticamente il numero di condizione del sistema.

3. Contributi Chiave

1. Sviluppo del Metodo C0-WG:

   • Introduzione di una formulazione basata su elementi $P_2$ continui che semplifica l'assemblaggio della matrice di rigidezza eliminando i termini di traccia presenti nei metodi $C^0$-IP.
   • Dimostrazione che il metodo è indipendente dai parametri (nessun parametro di penalità da regolare), garantendo stabilità e facilità di implementazione.
   • Fornitura di un'analisi di errore rigorosa che stabilisce la ben posta e la precisione del metodo per OCP vincolati, collegandolo ai metodi esistenti.

2. Progettazione del Precondizionatore:

   • Sviluppo di un precondizionatore di Schwarz additivo specifico per la formulazione C0-WG.
   • Analisi teorica che dimostra come il precondizionatore riduca il numero di condizione del sistema da $O(h^{-4})$ (senza precondizionamento) a un ordine significativamente inferiore (dipendente dalla sovrapposizione $\delta$), rendendo possibili solutori iterativi scalabili.

3. Analisi di Convergenza:

   • Stime di errore che confermano l'ordine di convergenza ottimale nella norma dipendente dalla mesh, tenendo conto della ridotta regolarità della soluzione dovuta ai vincoli e al tracking generale.

4. Risultati Numerici

Gli esperimenti numerici, condotti sia su problemi biarmonici che su OCP vincolati in 2D, confermano l'efficacia del metodo:

• Accuratezza: Il metodo C0-WG mostra ordini di convergenza superiori o uguali al metodo $C^0$-IP in diverse norme ($L^\infty, L^2, H^1$), producendo approssimazioni più accurate dello stato ottimo e del controllo.
• Efficienza del Precondizionatore:
  • Senza precondizionamento, il numero di condizione cresce come $O(h^{-4})$.
  • Con il precondizionatore di Schwarz additivo, il numero di condizione viene ridotto di oltre due ordini di grandezza.
  • L'uso di una sovrapposizione maggiore ($\delta = 2h$ rispetto a $\delta = h$) migliora ulteriormente la stabilità, mantenendo la crescita del numero di condizione sotto $O(h^{-3})$.
• Robustezza: Il metodo dimostra robustezza sia per problemi biarmonici standard che per OCP complessi con vincoli puntuali.

5. Significato e Impatto

Questo lavoro presenta un quadro computazionale robusto ed efficiente per una classe di problemi di controllo ottimo di alto ordine che sono tradizionalmente difficili da trattare numericamente.

• Semplificazione Implementativa: Rimuovendo la necessità di parametri di penalità e termini di traccia complessi, il metodo C0-WG riduce l'overhead computazionale e il rischio di errori di implementazione rispetto ai metodi $C^0$-IP.
• Scalabilità: L'integrazione con il precondizionatore di Schwarz additivo rende il metodo adatto per calcoli su larga scala e paralleli, risolvendo il collo di bottiglia principale rappresentato dai sistemi mal condizionati.
• Generalità: La metodologia proposta, pur essendo dimostrata in 2D, è estendibile a 3D e fornisce una base solida per futuri sviluppi in contesti dipendenti dal tempo o non lineari.

In sintesi, la combinazione del metodo Weak Galerkin $C^0$ con un precondizionatore di dominio offre una soluzione superiore per la simulazione di problemi di controllo ottimo vincolati, bilanciando accuratezza teorica, efficienza numerica e facilità d'uso.