On the application of the SCD semismooth* Newton method to solving Stokes problem with stick-slip boundary conditions

Titolo: Applicazione del metodo Newton SCD semismooth* alla risoluzione del problema di Stokes con condizioni al contorno stick-slip

1. Il Problema

Il documento affronta la risoluzione numerica del problema di Stokes tridimensionale per fluidi incomprimibili, caratterizzato da condizioni al contorno di tipo stick-slip (aderenza-scivolamento) di Navier-Tresca.

• Contesto fisico: In molte applicazioni reali, la velocità del fluido al confine non è necessariamente zero (condizione di "no-slip"), ma può presentare scivolamento se lo sforzo di taglio supera una certa soglia.
• Formulazione matematica: Il problema è modellato come un'equazione di Stokes accoppiata a condizioni di non-penetrazione e a una condizione di attrito non liscia (Tresca/Navier-Tresca).
• Sfida numerica: La formulazione debole del problema porta a una disuguaglianza variazionale di secondo tipo accoppiata a un vincolo di uguaglianza (incomprimibilità). La natura non liscia delle condizioni di stick-slip rende difficile l'applicazione diretta dei metodi Newton classici, che richiedono la differenziabilità delle funzioni.

2. Metodologia

Gli autori propongono un approccio basato sul recente metodo Newton SCD semismooth* (denominato SSSN), sviluppato da Gfrerer e Outrata.

• Discretizzazione:

  • Viene utilizzata la metodologia agli elementi finiti misti con la coppia di elementi P1-bubble/P1 (noti anche come elementi MINI). Questa scelta garantisce la stabilità inf-sup (condizione di Babuška-Brezzi) e riduce i gradi di libertà rispetto ad altre coppie (es. P2-P1).
  • Il problema continuo viene trasformato in un sistema algebrico di equazioni generalizzate (GE).

• Riformulazione come Equazione Generalizzata (GE):

  • Il problema viene scritto nella forma $0 \in F(x) = H(x) + \partial q(x)$, dove $H$ è un operatore affine lineare e $\partial q$ è il sottodifferenziale di una funzione convessa propria e semicontinua inferiormente (lsc), che incorpora le condizioni di stick-slip e di impermeabilità.
  • A differenza dei metodi tradizionali che richiedono una riformulazione o una regolarizzazione (smoothing) del problema non liscio, l'approccio SSSN lavora direttamente sulla struttura non liscia.

• Il Metodo SCD Semismooth:*

  • Il metodo si basa sul concetto di derivata contenente sottospazi (SCD - Subspace Containing Derivative).
  • Invece di calcolare un Jacobiano classico, il metodo costruisce un'approssimazione lineare efficiente utilizzando sottospazi che contengono la derivata del mappaggio multivalore ($\partial q$).
  • Passaggi dell'algoritmo:
    1. Passo di approssimazione: Calcolo di un punto tramite l'applicazione dell'operatore di prossimità (proximal mapping) su una funzione separabile.
    2. Passo Newton: Risoluzione di un sistema lineare che utilizza una rappresentazione basata su sottospazi (matrici $P$ e $W$) per approssimare la derivata del sottodifferenziale.
  • Globalizzazione: Per garantire la convergenza da punti di partenza arbitrari, il metodo è integrato con una tecnica ibrida che utilizza il metodo di Douglas-Rachford e una ricerca lineare basata su un residuo specifico.

3. Contributi Chiave

• Nuovo approccio numerico: Introduzione del metodo SSSN al problema di Stokes con condizioni di Navier-Tresca, evitando la necessità di trasformare il problema in un sistema di equazioni lisce o di utilizzare moltiplicatori di Lagrange complessi come nei metodi a set attivo tradizionali.
• Analisi teorica: Dimostrazione che il sottodifferenziale associato alle condizioni di stick-slip possiede la proprietà SCD e che il metodo Newton SCD semismooth* converge superlinearmente localmente.
• Implementazione efficiente: Sviluppo di un algoritmo dettagliato che sfrutta la struttura a blocchi delle matrici derivate dalle condizioni al contorno, permettendo l'eliminazione delle variabili interne e la risoluzione efficiente dei sistemi lineari tramite GMRES con precondizionamento.
• Confronto con lo stato dell'arte: Il metodo non richiede la riformulazione del problema in forma di equazione, a differenza di molte tecniche esistenti.

4. Risultati Sperimentali

Gli autori hanno testato l'efficienza del metodo su due casi di studio:

1. Cubo Unitario (Test di convergenza):

   • Un dominio cubico con diverse configurazioni dei parametri di adesione ($\kappa$) e limite di scivolamento ($g$).
   • Sono stati testati tre scenari: puro scivolamento ($g=0$), misto stick-slip ($g=5$) e pura aderenza ($g=10$).
   • Risultati: Il numero di iterazioni Newton (SSSN) è risultato indipendente dalla dimensione della mesh, dimostrando la robustezza del metodo. La convergenza è stata raggiunta in un numero molto ridotto di iterazioni (tra 3 e 8) per diverse risoluzioni della griglia.

2. Aneurisma Cerebrale (Geometria complessa):

   • Applicazione su una geometria 3D complessa tratta da un database medico reale.
   • Sono stati monitorati i parametri di scivolamento e le iterazioni GMRES interne.
   • Risultati: Il metodo ha mostrato una rapida convergenza anche in geometrie complesse. Il confronto con il metodo Newton semismooth* presentato in letteratura ([9]) indica che l'approccio SSSN offre tassi di convergenza leggermente migliori e non necessita di castare il problema in forma di equazione.

5. Significato e Implicazioni

• Efficienza Computazionale: Il metodo SSSN si dimostra altamente efficiente per problemi di fluidodinamica con condizioni al contorno non lisce, riducendo il numero di iterazioni necessarie rispetto ai metodi tradizionali.
• Versatilità: L'algoritmo proposto è modulare e può essere esteso con modifiche relativamente semplici per risolvere problemi di Navier-Stokes (inclusi quelli non stazionari), aprendo la strada a simulazioni più realistiche di flussi complessi.
• Rilevanza Teorica: Conferma l'applicabilità della teoria delle mappe SCD e dei metodi Newton generalizzati a problemi di meccanica dei fluidi con vincoli di contatto e attrito, offrendo un'alternativa solida ai metodi di regolarizzazione che possono introdurre errori di approssimazione.

In sintesi, il paper presenta un avanzamento significativo nella risoluzione numerica dei problemi di Stokes con condizioni di stick-slip, combinando una rigorosa analisi teorica con un'implementazione pratica che dimostra superiorità in termini di velocità di convergenza e indipendenza dalla mesh.