Sulle onde lineari con memoria in un mezzo di tipo Bessel

Informazioni Fondamentali

• ID Articolo: 2510.11493
• Titolo: On linear waves with memory in a Bessel-like medium
• Autori: A. Giusti, I. Colombaro, A. Mentrelli
• Classificazione: math-ph math.MP
• Data di Pubblicazione: 14 ottobre 2025
• Link Articolo: https://arxiv.org/abs/2510.11493

Riassunto

Il presente articolo affronta la propagazione di onde armoniche e transienti in sistemi governati da equazioni d'onda con memoria, dove il nucleo integrale nel dominio di Laplace coinvolge il rapporto tra funzioni di Bessel modificate di primo tipo. In particolare, lo studio delle onde transienti viene condotto mediante un metodo completamente numerico basato sul metodo di Talbot, impiegato per l'inversione numerica della trasformata di Laplace.

Contesto di Ricerca e Motivazione

Descrizione del Problema

Il problema centrale affrontato in questa ricerca riguarda il comportamento della propagazione ondosa in mezzi di tipo Bessel con effetti di memoria. Specificamente, l'oggetto di studio è l'equazione d'onda con memoria della forma:

$\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi_\nu(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0$

dove il nucleo di memoria nel dominio di Laplace è rappresentato come: $\tilde{\Phi}_\nu(s) = \frac{2(\nu+1)}{\sqrt{s\tau}}\frac{I_{\nu+1}(\sqrt{s\tau})}{I_\nu(\sqrt{s\tau})}$

Importanza e Applicazioni

Questa classe di equazioni presenta importanti applicazioni in molteplici campi fisici:

1. Emodinamica: descrizione della propagazione ondosa in tubi elastici pieni di liquido
2. Meccanica Viscoelastica: modellazione degli effetti di memoria nei materiali
3. Sistemi Elettrici: analisi di circuiti con caratteristiche di memoria

Limitazioni dei Metodi Esistenti

Le ricerche precedenti si basavano principalmente sull'algoritmo semi-analitico di Buchen-Mainardi, il quale risulta affidabile solo in prossimità del fronte d'onda e non fornisce un'immagine completa dell'evoluzione spazio-temporale.

Motivazione della Ricerca

Il presente lavoro mira a superare le limitazioni dei metodi esistenti attraverso approcci numerici, fornendo un'analisi completa del comportamento della propagazione ondosa nell'intero dominio spazio-temporale.

Contributi Fondamentali

1. Sviluppo di un quadro completo di analisi dispersiva: per il caso ν=0, analisi dettagliata della relazione di dispersione per onde armoniche con attenuazione spaziale e periodicità temporale (SATP)
2. Derivazione di leggi di dispersione esatte e fattori di attenuazione: mediante espressioni esplicite fornite dalle funzioni di Kelvin
3. Implementazione di analisi numerica completa delle onde transienti: utilizzo del metodo di Talbot per il calcolo numerico dell'antitrasformata di Laplace
4. Rivelazione di caratteristiche di dispersione anomala: dimostrazione del fenomeno di dispersione anomala dove la velocità di fase è inferiore alla velocità di gruppo

Descrizione Dettagliata dei Metodi

Definizione del Compito

Lo studio delle caratteristiche di propagazione ondosa in mezzi di memoria di tipo Bessel, includendo:

• Analisi delle proprietà dispersive e di attenuazione delle onde armoniche
• Calcolo della risposta al gradino per onde transienti
• Determinazione della velocità di fase e della velocità di gruppo

Architettura del Modello

1. Equazione d'Onda Fondamentale

Considerando il caso particolare ν=0: $\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0$

dove: $\tilde{\Phi}(s) = \frac{2}{\sqrt{s\tau}}\frac{I_1(\sqrt{s\tau})}{I_0(\sqrt{s\tau})}$

2. Metodo di Analisi Dispersiva

Utilizzo della forma standard di soluzione armonica: $Y(t,x) = A e^{i(\omega t - kx)}$

ottenendo la relazione di dispersione: $(i\omega)^2 + k^2 c^2[1 - \hat{\Phi}(\omega)] = 0$

3. Analisi di Onde SATP

Per onde con attenuazione spaziale e periodicità temporale, definendo:

• $\omega \in \mathbb{R}$ (frequenza reale)
• $k = \kappa - i\delta_{att}$ (numero d'onda complesso)

Risolvendo il sistema di equazioni non lineari: $\begin{cases} \kappa^2 - \delta_{att}^2 = A(\omega) \\ -2\kappa\delta_{att} = B(\omega) \end{cases}$

dove A(ω) e B(ω) sono espressi mediante funzioni di Kelvin.

Punti di Innovazione Tecnica

1. Rappresentazione mediante Funzioni di Kelvin

Conversione del rapporto di funzioni di Bessel complesse in combinazioni di funzioni di Kelvin reali, evitando l'instabilità numerica dei calcoli complessi.

2. Applicazione del Metodo di Talbot

Utilizzo del metodo di Talbot per l'antitrasformata di Laplace: $Y(t,x) = \frac{1}{2\pi i}\int_C e^{st}\tilde{Y}(s,x)ds$

Mediante deformazione del percorso di integrazione, convertendo l'integrale di Bromwich in un integrale sull'asse reale, migliorando la stabilità numerica.

3. Soluzione Numerica Globale

Superamento della limitazione dei metodi tradizionali validi solo in prossimità del fronte d'onda, realizzando calcoli precisi nell'intero dominio spazio-temporale.

Configurazione Sperimentale

Problema ai Valori Iniziali e al Contorno

Considerando il problema di risposta al gradino su un dominio semi-infinito: $\begin{cases} Y_{tt} - c^2[1 - \Phi(t)*]Y_{xx} = 0, & x > 0, t > 0 \\ Y(t,0) = H(t), & t > 0 \\ Y(0,x) = 0, Y_t(0,x) = 0, & x > 0 \end{cases}$

Implementazione Numerica

• Utilizzo di Wolfram Mathematica® per calcoli simbolici
• Metodo delle differenze finite per il calcolo numerico della derivata della velocità di gruppo
• Integrazione numerica ad alta precisione del metodo di Talbot

Impostazione dei Parametri

Utilizzo di parametri adimensionali:

• Parametro di frequenza: $\omega\tau$
• Parametro spaziale: $x/(c\tau)$
• Parametro temporale: $(ct-x)/(c\tau)$

Risultati Sperimentali

Risultati Principali

1. Legge di Dispersione e Fattore di Attenuazione

La Figura 1 mostra la legge di dispersione adimensionale $c\tau\kappa$ e il fattore di attenuazione $c\tau\delta_{att}$ in funzione di $\omega\tau$:

• La legge di dispersione presenta una crescita non lineare
• Il fattore di attenuazione raggiunge un picco a frequenze intermedie, quindi decresce

2. Velocità di Fase e Velocità di Gruppo

La Figura 2 presenta il confronto tra velocità di fase e velocità di gruppo:

• Velocità di fase: $v_p(\omega) = \frac{\omega}{\kappa(\omega)}$
• Velocità di gruppo: $v_g(\omega) = \left[\frac{d\kappa(\omega)}{d\omega}\right]^{-1}$
• I risultati numerici mostrano: $v_p(\omega) < v_g(\omega) < c$, confermando le caratteristiche di dispersione anomala

3. Propagazione di Onde Transienti

La Figura 3 mostra la risposta al gradino a diverse posizioni:

• Il fronte d'onda si propaga alla velocità c
• La forma d'onda subisce dispersione significativa durante la propagazione
• Le regioni lontane dal fronte d'onda presentano comportamenti oscillatori complessi con attenuazione

Scoperte Importanti

1. Conferma della Dispersione Anomala

Il sistema presenta caratteristiche tipiche di dispersione anomala:

• $v_p < v_g < c$
• Entrambe le velocità tendono a c nel limite di alta frequenza

2. Effetti della Memoria

Il termine di memoria produce:

• Attenuazione e dispersione dell'onda
• Comportamento transienti complesso
• Evoluzione della forma d'onda di tipo non esponenziale

Lavori Correlati

Sviluppo Storico

1. Mainardi et al. (1982-1983): introduzione iniziale di questa classe di modelli in emodinamica
2. Giusti & Mainardi (2016): stabilimento dell'analogia con viscoelasticità dinamica
3. Colombaro et al. (2017): modelli di viscoelasticità lineare basati su funzioni di Bessel

Confronto dei Metodi

• Algoritmo di Buchen-Mainardi: semi-analitico, valido solo presso il fronte d'onda
• Metodo del presente articolo: completamente numerico, valido nell'intero dominio
• Metodi asintotici: approssimazione ad alta frequenza, limitazioni significative

Conclusioni e Discussione

Conclusioni Principali

1. Caratteristiche Dispersive: il sistema presenta dispersione anomala, con velocità di fase sempre inferiore alla velocità di gruppo e entrambe non superiori alla velocità del fronte d'onda c
2. Validità del Metodo Numerico: il metodo di Talbot supera con successo le limitazioni dei metodi tradizionali
3. Effetti di Memoria: il nucleo di memoria di tipo Bessel produce comportamenti complessi di propagazione ondosa

Limitazioni

1. Complessità Computazionale: la fase transienti richiede un numero elevato di nodi, con costi computazionali significativi
2. Restrizioni Parametriche: considerazione limitata al caso particolare ν=0
3. Insufficienza dell'Analisi Teorica: difficoltà nella dimostrazione analitica rigorosa delle proprietà della velocità di gruppo

Direzioni Future

1. Metodi Ibridi: combinazione di metodi semi-analitici per il fronte d'onda e metodi numerici per il campo lontano
2. Generalizzazione: estensione al caso di valori arbitrari di ν
3. Ampliamento Applicativo: verifica dell'applicazione in sistemi fisici concreti

Valutazione Approfondita

Punti di Forza

1. Innovazione Metodologica: prima applicazione del metodo di Talbot al problema della propagazione ondosa in mezzi di memoria di tipo Bessel
2. Completezza Teorica: fornisce un quadro completo dall'analisi dispersiva al calcolo transienti
3. Precisione Numerica: supera le limitazioni dei metodi esistenti, realizzando calcoli ad alta precisione nell'intero dominio
4. Intuizioni Fisiche: rivela fenomeni fisici importanti come la dispersione anomala

Carenze

1. Efficienza Computazionale: l'elevata complessità del calcolo transienti limita le applicazioni pratiche
2. Intervallo Parametrico: limitazione a forme specifiche di nuclei di memoria
3. Verifica Sperimentale: mancanza di confronto con esperimenti fisici reali

Impatto

1. Contributo Teorico: fornisce nuovi strumenti di analisi numerica per la propagazione ondosa in mezzi con memoria
2. Valore Applicativo: potenziali applicazioni in biomedicina, scienza dei materiali e altri campi
3. Significato Metodologico: l'applicazione del metodo di Talbot fornisce un riferimento per problemi analoghi

Scenari Applicabili

1. Ricerca Accademica: teoria della propagazione ondosa, modellazione di materiali con memoria
2. Applicazioni Ingegneristiche: analisi del flusso in condotte, controllo delle vibrazioni
3. Calcolo Numerico: calcolo ad alta precisione dell'antitrasformata di Laplace

Bibliografia

L'articolo cita 23 importanti riferimenti bibliografici, che coprono:

• Teoria delle funzioni di Bessel (Abramowitz & Stegun, 1972)
• Teoria della propagazione ondosa (Mainardi, 2022; Whitham, 2011)
• Metodi numerici (Talbot, 1979; Garrappa et al., 2021)
• Contesti applicativi (letteratura relativa a emodinamica, viscoelasticità, etc.)

Il presente articolo fornisce contributi significativi sia nell'analisi teorica che nei metodi numerici, offrendo nuovi strumenti e intuizioni per la ricerca sulla propagazione ondosa in mezzi con memoria. Sebbene presenti alcune limitazioni in termini di efficienza computazionale e intervallo parametrico, il suo valore metodologico e le sue intuizioni fisiche lo rendono un progresso importante in questo campo di ricerca.