Numerical stabilization for a mixture system with kind damping

Titolo: Stabilizzazione Numerica per un Sistema di Miscela con Smorzamento di Tipo Frazionario

Autori: K. Ammari, V. Komornik, M. Sepúlveda-Cortés, O. Vera-Villagrán.

---

1. Problema e Contesto

Il lavoro si concentra sull'analisi numerica e teorica della stabilità asintotica di un sistema di equazioni alle derivate parziali che modella la dinamica di una miscela di due solidi rigidi con porosità.

• Modello Fisico: Il sistema descrive l'interazione tra un campo di porosità ($u$) e un campo di temperatura ($v$) all'interno di un corpo elastico isotropo omogeneo.
• Equazioni Governative: Il modello originale include termini dissipativi classici (derivate prime rispetto al tempo moltiplicate per il Laplaciano). Tuttavia, gli autori sostituiscono questi termini con derivate frazionarie di tipo Caputo (con un peso esponenziale), indicate come $\partial^{\alpha, \eta}_t u$ e $\partial^{\alpha, \eta}_t v$.
• Obiettivo: Stabilire le condizioni per la stabilità forte e quantificare il tasso di decadimento polinomiale delle soluzioni, validando teoricamente e numericamente l'efficacia di questo tipo di smorzamento "di tipo frazionario" (spesso associato a materiali viscoelastici come le miscele bituminose).

2. Metodologia

A. Approccio Teorico (Analisi Funzionale)

1. Riformulazione del Modello (Modello Augmentato):

   • La presenza di derivate frazionarie rende difficile l'analisi diretta nello spazio di Hilbert standard. Gli autori utilizzano una tecnica di "augmentazione" basata su un teorema di rappresentazione (Teorema 2.1, riferito a [26]).
   • Le derivate frazionarie vengono sostituite da un sistema accoppiato che include equazioni di diffusione in una variabile ausiliaria $\xi \in \mathbb{R}$. Il sistema originale viene trasformato in un sistema di evoluzione del primo ordine in un spazio di Hilbert esteso $\mathcal{H}$.
   • Il nuovo sistema include variabili aggiuntive $\phi_1$ e $\phi_2$ che soddisfano equazioni di tipo parabolico: $\phi_t + (\xi^2 + \eta)\phi = \mu(\xi)u_t$.

2. Ben-posedness e Semigruppo:

   • Viene definito un operatore lineare non limitato $A$ sullo spazio di stato $\mathcal{H}$.
   • Si dimostra che $A$ è il generatore di un semigruppo di contrazioni $C_0$, garantendo l'esistenza e l'unicità delle soluzioni per dati iniziali appropriati (Teorema 3.2).

3. Analisi di Stabilità:

   • Stabilità Forte: Viene applicato il teorema di Arendt-Batty. Si dimostra che lo spettro dell'operatore $A$ non contiene autovalori sull'asse immaginario puro (eccetto eventualmente l'origine) e che l'intersezione con l'asse immaginario è al più numerabile. Questo garantisce che le soluzioni convergano asintoticamente a zero.
   • Stabilità Polinomiale: Si analizza il tasso di decadimento dell'energia.
     • Per $\eta > 0$: Il decadimento è polinomiale con tasso $O(t^{-1/(1-\alpha)})$.
     • Per $\eta = 0$: Si utilizza un risultato di [8] basato sulla stima della risolvente $(i\lambda I - A)^{-1}$ per dimostrare un decadimento polinomiale, confermando che il sistema non è esponenzialmente stabile (tipico dei sistemi con smorzamento frazionario).

B. Approccio Numerico

1. Discretizzazione Spaziale:

   • Utilizzo del Metodo dei Volumi Finiti (FVM) per discretizzare le variabili spaziali su un dominio bidimensionale.
   • Le derivate spaziali (Laplaciano) sono approssimate tramite differenze finite.

2. Discretizzazione Temporale:

   • Adozione dello schema Newmark-$\beta$ per le equazioni del moto lineari, scelto per preservare l'energia e garantire stabilità numerica di secondo ordine.
   • Le equazioni ausiliarie per le variabili frazionarie ($\phi$) sono discretizzate con il metodo Crank-Nicolson per mantenere la conservazione dell'energia o la sua non-crescita.

3. Approssimazione delle Derivate Frazionarie:

   • L'integrale che definisce l'operatore frazionario è approssimato mediante una quadratura numerica su una griglia discreta della variabile $\xi$.

3. Risultati Chiave

Risultati Teorici

• Stabilità Asintotica: È stato provato rigorosamente che il sistema è fortemente stabile; l'energia del sistema decade a zero per $t \to \infty$.
• Tasso di Decadimento: Il sistema non è esponenzialmente stabile. Il tasso di decadimento è polinomiale.
  • La presenza del parametro $\eta$ influenza il tasso: per $\eta > 0$ il decadimento è più rapido rispetto al caso $\eta = 0$, ma entrambi seguono una legge di potenza $t^{-\gamma}$.
• Spettro dell'Operatore: L'analisi spettrale conferma che gli autovalori si accumulano vicino all'asse immaginario, il che è coerente con la mancanza di stabilità esponenziale e la presenza di decadimento polinomiale.

Risultati Numerici

• Conferma del Decadimento: Le simulazioni numeriche mostrano chiaramente che l'energia del sistema decade nel tempo quando sono presenti i termini dissipativi frazionari, mentre rimane costante (sistema conservativo) in loro assenza.
• Analisi Spettrale: La visualizzazione numerica dello spettro dell'operatore $A$ (Figura 3) mostra che tutti gli autovalori giacciono nel semipiano sinistro del piano complesso, confermando la stabilità. Gli autovalori con modulo più grande sono quelli più vicini all'asse immaginario.
• Decadimento Polinomiale: L'analisi log-log dell'energia nel tempo (Figura 5) mostra una relazione lineare, confermando il comportamento polinomiale.
  • Per $\alpha = 0.5$, i risultati numerici si allineano bene con le curve teoriche $C/t$ (per $\eta=0$) e $C/t^2$ (per $\eta>0$), sebbene le simulazioni mostrino una tendenza generale verso il comportamento $1/t$ indipendentemente da $\eta$ in certi regimi, confermando la natura sub-esponenziale del decadimento.
• Interazione delle Onde: Le simulazioni mostrano l'interazione di onde gaussiane iniziali. La differenza qualitativa tra il caso con e senza smorzamento è sottile nella forma d'onda istantanea, ma evidente nel comportamento energetico a lungo termine.

4. Contributi Principali

1. Modellizzazione Ibrida: Integrazione efficace di derivate frazionarie in un sistema di miscele di solidi porosi, offrendo un modello più realistico per materiali viscoelastici complessi.
2. Dimostrazione Rigorosa: Fornisce una prova matematica completa della stabilità forte e del tasso di decadimento polinomiale per questo specifico sistema accoppiato, utilizzando la teoria dei semigruppi e l'analisi spettrale.
3. Validazione Numerica: Sviluppo di uno schema numerico conservativo (FVM + Newmark + Crank-Nicolson) che riesce a catturare correttamente la dinamica a lungo termine e il decadimento polinomiale, validando le previsioni teoriche.
4. Analisi Parametrica: Studio dell'influenza del parametro $\eta$ (peso esponenziale) e dell'ordine frazionario $\alpha$ sulla dinamica del sistema.

5. Significato e Implicazioni

Questo lavoro è significativo perché:

• Teoria dei Materiali: Fornisce una base matematica solida per l'uso di derivate frazionarie nella modellazione di materiali porosi e viscoelastici, spiegando perché tali materiali non mostrano un decadimento esponenziale rapido delle vibrazioni ma un decadimento più lento (polinomiale).
• Ingegneria Strutturale: I risultati sono rilevanti per la progettazione di strutture in materiali compositi o bituminosi (es. pavimentazioni stradali), dove la dissipazione dell'energia è cruciale per la stabilità a lungo termine.
• Metodologia Numerica: Dimostra che è possibile simulare accuratamente sistemi con derivate frazionarie accoppiate a equazioni d'onda senza perdere la proprietà di stabilità energetica, offrendo uno strumento affidabile per la previsione del comportamento dinamico di tali sistemi.

In sintesi, il paper collega teoria analitica avanzata e simulazione numerica per dimostrare che lo smorzamento di tipo frazionario garantisce la stabilità di sistemi di miscele complesse, sebbene con un tasso di decadimento più lento rispetto alla dissipazione classica.