Pulse waves in the viscoelastic Kelvin-Voigt model: a revisited approach

Questo articolo presenta un nuovo approccio risolutivo per le onde d'impulso nel modello viscoelastico di Kelvin-Voigt, fornendo espressioni integrali più semplici ed efficienti dal punto di vista computazionale rispetto alle precedenti, che evitano il calcolo numerico della trasformata inversa di Laplace e permettono di derivare formule asintotiche per eccitazioni a impulso e a gradino.

L'essenziale

Ecco una spiegazione semplice e creativa del paper, pensata per chiunque, anche senza un background scientifico.

🌊 Il Mistero dell'Impulso nel "Miele" Viscoelastico

Immagina di avere un materiale speciale, un po' come un miele molto denso mescolato con una molla. Se lo colpisci, non si comporta come l'acqua (che si muove subito e basta) e nemmeno come un blocco di metallo rigido (che vibra e basta). Si chiama modello di Kelvin-Voigt ed è usato per descrivere cose come la terra durante un terremoto o certi tessuti biologici.

Il problema che gli scienziati (González Santander, Mainardi e Mentrelli) hanno affrontato è questo:

"Se colpisco questo materiale in un punto (l'origine) con un impulso improvviso, come si muove l'onda attraverso il materiale nel tempo?"

🚧 Il Problema Vecchio: La "Scatola Nera" Matematica

Fino ad ora, per calcolare come si muove questa onda, gli scienziati usavano un metodo matematico molto complicato chiamato Trasformata di Laplace.
Immagina di dover risolvere un puzzle guardando solo il retro dei pezzi. Devi fare un calcolo inverso (una "trasformata inversa") che è come cercare di indovinare la forma originale di un oggetto solo guardando la sua ombra proiettata su un muro.

• Il difetto: Per ottenere il risultato, dovevi usare computer potenti per fare calcoli complessi nel "piano complesso" (un mondo matematico astratto). Era lento, difficile da programmare e prone a errori numerici.

💡 La Nuova Soluzione: La "Ricetta Diretta"

Questi tre ricercatori hanno trovato un modo migliore. Hanno derivato una nuova formula integrale.

• L'analogia: Invece di dover indovinare l'ombra per capire l'oggetto, hanno trovato una ricetta diretta. Ora puoi calcolare esattamente come si muove l'onda semplicemente integrando una funzione (un'operazione matematica più semplice e diretta) senza dover saltare nel mondo astratto delle trasformate inverse.

È come se prima dovessi cucinare una torta usando un forno a microonde che devi programmare con codici binari, e ora hai scoperto che puoi semplicemente mescolare gli ingredienti in una ciotola e infornarla. È più veloce, più preciso e molto più facile da usare.

🎯 Due Casi Pratici: Il "Colpo Secco" e la "Spinta Continua"

Gli autori hanno applicato questa nuova ricetta a due scenari comuni:

1. L'Impulso Delta (Il Colpo Secco): Immagina di dare un pugno secco al materiale. La loro nuova formula ti dice esattamente come quell'onda di pugno si allontana e si smorza nel tempo.
2. L'Impulso Gradino (La Spinta Continua): Immagina di spingere il materiale e mantenerlo premuto. Anche qui, la nuova formula descrive perfettamente come il materiale si deforma e si stabilizza.

🔍 Cosa succede all'infinito? (Le Approssimazioni)

Oltre alla ricetta precisa, hanno anche creato delle stime veloci per due situazioni estreme:

• Appena dopo il colpo (Tempo breve): Come si comporta l'onda nei primi istanti? La loro formula dice che si comporta in un modo specifico, simile a come si diffonde il calore.
• Dopo molto tempo (Tempo lungo): Cosa succede quando l'onda si è stabilizzata? La formula mostra come il materiale torna a riposo o raggiunge un equilibrio.

Hanno verificato che le loro nuove formule "veloci" danno esattamente lo stesso risultato delle vecchie formule "complesse", ma con un errore numerico quasi nullo (come dimostrato dai grafici nel paper).

🏁 In Sintesi

Questo articolo è come se gli scienziati avessero trovato una mappa più chiara per navigare in un territorio matematico accidentato.

• Prima: "Dobbiamo calcolare l'inverso di una trasformata complessa per vedere l'onda." (Difficile, lento).
• Ora: "Ecco una formula integrale diretta che ci dà il risultato immediatamente." (Semplice, veloce, efficiente).

Questo è fondamentale per campi come la sismologia (studio dei terremoti), perché permette di simulare come le onde si muovono attraverso la Terra in modo molto più rapido e affidabile, aiutandoci a capire meglio i rischi e il comportamento dei materiali sotto stress.

Il messaggio finale: A volte, per risolvere un problema vecchio di decenni, non serve inventare una nuova fisica, ma solo trovare un modo più intelligente e semplice per scrivere la matematica che la descrive.

Sommario tecnico

Ecco una sintesi tecnica dettagliata del paper "Pulse waves in the viscoelastic Kelvin-Voigt model: a revisited approach", redatta in italiano.

Titolo: Onde d'impulso nel modello viscoelastico di Kelvin-Voigt: un approccio rivisitato

Autori: Juan Luis González-Santander, Francesco Mainardi, Andrea Mentrelli.
Pubblicazione: Mathematics (MDPI), Vol. 14 (2025).

---

1. Il Problema

Il lavoro affronta il problema fondamentale della meccanica dei solidi viscoelastici: determinare la risposta meccanica $r(x, t)$ di un mezzo omogeneo semi-infinito, inizialmente a riposo, soggetto a un impulso applicato all'origine ($x=0$).
Il modello specifico utilizzato è il modello di Kelvin-Voigt, una delle configurazioni più semplici di viscoelasticità lineare, rilevante in sismologia e nella propagazione di onde transitorie.
Il problema classico richiede l'inversione della trasformata di Laplace di una funzione complessa che coinvolge la radice quadrata di un polinomio in $s$. Le soluzioni esistenti in letteratura spesso richiedono:

• Calcoli numerici complessi dell'inversione di Laplace nel piano complesso.
• Integrali che risultano computazionalmente inefficienti o difficili da gestire per l'analisi asintotica (specialmente per tempi lunghi o distanze grandi).

2. Metodologia

Gli autori adottano un approccio analitico rigoroso basato sulle seguenti fasi:

• Formulazione nel dominio di Laplace: Partendo dalle equazioni costitutive del modello di Kelvin-Voigt ($\sigma = G_e [\varepsilon + t_\varepsilon \dot{\varepsilon}]$), viene derivata l'equazione differenziale nel dominio di Laplace per la risposta $r(x,s)$. La soluzione formale è data da:
  $$\tilde{r}(x, s) = \tilde{r}_0(s) \exp\left(-\frac{s}{c'} \sqrt{1 + t_\varepsilon s} \, x\right)$$
  dove $c'$ è una velocità caratteristica e $t_\varepsilon$ il tempo di ritardo.

• Introduzione di variabili adimensionali: Per semplificare l'analisi, vengono introdotte le variabili adimensionali:
  $$\xi = \frac{x}{c' t_\varepsilon}, \quad \tau = \frac{t}{t_\varepsilon}$$

• Derivazione della soluzione integrale: Invece di calcolare direttamente l'inversione di Laplace numerica, gli autori utilizzano il teorema di convoluzione e sviluppi in serie di funzioni speciali (funzioni di Hermite e funzioni ipergeometriche confluenti $0F_1$).
  Vengono definiti un nucleo di Green $g(\xi, \tau)$ e una funzione ausiliaria $f(\xi, \tau)$, permettendo di esprimere la risposta generale come una convoluzione temporale:
  $$r(\xi, \tau) = \int_0^\tau r_0(t_\varepsilon(\tau - \tau')) e^{-\tau'} g(\xi, \tau') d\tau'$$

• Casi specifici: La metodologia viene applicata a due casi di eccitazione fondamentali:

  1. Impulso Delta ($r_0(t) = \delta(t)$): Risoluzione esplicita tramite la derivata temporale del nucleo.
  2. Impulso Gradino ($r_0(t) = \theta(t)$): Risoluzione tramite integrazione per parti e utilizzo di funzioni complementari dell'errore (erfc).

• Analisi Asintotica: Vengono derivati sviluppi asintotici rigorosi per i limiti $\tau \to 0, \infty$ e $\xi \to 0, \infty$, utilizzando espansioni in serie delle funzioni ipergeometriche generalizzate ($0F_2$).

3. Risultati Chiave

A. Soluzioni Integrali Nuove ed Efficienti

Gli autori hanno derivato nuove forme integrali per la risposta $r(\xi, \tau)$ che evitano l'integrazione numerica nel piano complesso.

• Per l'impulso gradino:
  $$r(\xi, \tau) = \int_\xi^\infty {}_0F_1\left(-; 1; \xi(v-\xi)\right) \left[ \frac{e^{-\frac{v^2}{4\tau}-\tau}}{\sqrt{\pi\tau}} + \frac{e^{-v}}{2}\text{erfc}\left(\frac{v-2\tau}{2\sqrt{\tau}}\right) - \frac{e^{v}}{2}\text{erfc}\left(\frac{v+2\tau}{2\sqrt{\tau}}\right) \right] dv$$
• Per l'impulso delta:
  $$r(\xi, \tau) = \frac{e^{-\tau}}{4\sqrt{\pi}\tau^{5/2}} \int_\xi^\infty (v^2 - 2\tau) e^{-\frac{v^2}{4\tau}} {}_0F_1\left(-; 1; \xi(v-\xi)\right) dv$$

B. Confronto Numerico ed Efficienza

Le nuove formule integrali sono state validate numericamente contro:

1. La soluzione esatta tramite inversione numerica di Laplace.
2. Soluzioni storiche di Hanin (1957) e Dozio (1990).

• Risultato: Le nuove formule mostrano una discrepanza numerica trascurabile ($< 10^{-9}$) rispetto all'inversione di Laplace, ma sono computazionalmente più efficienti e stabili, specialmente per grandi valori di $\tau$ (tempi lunghi), dove le formule precedenti (es. Morrison, 1956) diventano instabili o difficili da calcolare.
• È stato notato che la soluzione proposta da Dozio (1990) presenta errori numerici significativi rispetto ai risultati ottenuti in questo studio.

C. Formule Asintotiche

Il lavoro fornisce formule asintotiche semplici e chiuse per i limiti estremi, utili per l'analisi fisica rapida:

• Per $\tau \to 0$ o $\xi \to \infty$ (comportamento iniziale/diffusivo):
  $$r(\xi, \tau) \approx \text{erfc}\left(\frac{\xi}{2\sqrt{\tau}}\right) \quad (\text{gradino})$$
  $$r(\xi, \tau) \approx \frac{\xi}{2\sqrt{\pi}\tau^{3/2}} e^{-\frac{\xi^2}{4\tau} - \tau} \quad (\text{delta})$$
• Per $\tau \to \infty$ (comportamento a regime):
  Sono state ottenute espressioni in termini di funzioni ipergeometriche generalizzate ${}_0F_2$ che descrivono accuratamente la coda della risposta.

4. Significato e Contributi

Questo lavoro rivisita un problema classico della viscoelasticità con un approccio moderno e rigoroso, offrendo contributi significativi:

1. Efficienza Computazionale: Le nuove forme integrali superano le limitazioni delle soluzioni storiche, rendendo il calcolo della risposta transitoria più veloce e stabile, specialmente per simulazioni numeriche su lunghi intervalli temporali.
2. Chiarezza Analitica: La derivazione di formule asintotiche esplicite per tutti i limiti ($\xi, \tau \to 0, \infty$) facilita l'interpretazione fisica del comportamento dell'onda (diffusione istantanea vs. propagazione con fronte d'onda) senza ricorrere a simulazioni numeriche pesanti.
3. Validazione Critica: Lo studio mette in discussione e corregge implicitamente alcune formulazioni precedenti (come quella di Dozio) attraverso un confronto numerico rigoroso.
4. Strumento per la Sismologia: Poiché il modello di Kelvin-Voigt è ampiamente usato in sismologia per modellare l'attenuazione delle onde, queste nuove soluzioni forniscono strumenti più precisi per l'analisi della propagazione delle onde in mezzi dissipativi.

In sintesi, il paper fornisce un "nuovo standard" per il calcolo e l'analisi delle onde transitorie nel modello di Kelvin-Voigt, combinando eleganza matematica (uso di funzioni speciali) con robustezza numerica.