2025-11-16T20:16:12.924612

On linear waves with memory in a Bessel-like medium

Giusti, Colombaro, Mentrelli
We discuss the propagation of harmonic and transient waves for systems governed by a wave equation with memory whose integral kernel involves ratios of modified Bessel functions of the first kind in the Laplace domain. In particular, the investigation of transient waves is carried out by means of a fully numerical approach based on the Talbot method for the numerical inversion of Laplace transforms.
academic

On linear waves with memory in a Bessel-like medium

Basic Information

  • Paper ID: 2510.11493
  • Title: On linear waves with memory in a Bessel-like medium
  • Authors: A. Giusti, I. Colombaro, A. Mentrelli
  • Classification: math-ph math.MP
  • Publication Date: October 14, 2025
  • Paper Link: https://arxiv.org/abs/2510.11493

Abstract

This paper discusses the propagation of harmonic and transient waves in systems governed by wave equations with memory, where the integral kernel in the Laplace domain involves the ratio of modified Bessel functions of the first kind. In particular, the study of transient waves is conducted through a fully numerical method based on the Talbot method, which is employed for numerical inversion of the Laplace transform.

Research Background and Motivation

Problem Description

The core problem addressed in this research concerns wave propagation behavior in Bessel-like media with memory effects. Specifically, the research object is a wave equation with memory of the following form:

∂2Y∂t2−c2[1−Φν(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi_\nu(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

where the memory kernel in the Laplace domain is expressed as: Φ~ν(s)=2(ν+1)sτIν+1(sτ)Iν(sτ)\tilde{\Phi}_\nu(s) = \frac{2(\nu+1)}{\sqrt{s\tau}}\frac{I_{\nu+1}(\sqrt{s\tau})}{I_\nu(\sqrt{s\tau})}

Significance and Applications

This class of equations has important applications in multiple physical domains:

  1. Hemodynamics: Describing wave propagation in fluid-filled elastic tubes
  2. Viscoelastic mechanics: Modeling memory effects in materials
  3. Electrical systems: Analysis of circuits with memory characteristics

Limitations of Existing Methods

Previous research primarily employed the Buchen-Mainardi semi-analytical algorithm, which is reliable only near the wavefront and cannot provide a complete spatiotemporal evolution picture.

Research Motivation

This paper aims to overcome the limitations of existing methods through numerical approaches, providing a comprehensive analysis of wave propagation behavior throughout the entire spatiotemporal domain.

Core Contributions

  1. Established a complete dispersion analysis framework: For the ν=0 case, detailed analysis of the dispersion relation for spatially attenuated temporally periodic (SATP) harmonic waves
  2. Derived exact dispersion laws and attenuation factors: Provided explicit expressions through Kelvin functions
  3. Implemented fully numerical transient wave analysis: Employed the Talbot method for numerical computation of Laplace inverse transforms
  4. Revealed anomalous dispersion characteristics: Demonstrated that the system exhibits anomalous dispersion where phase velocity is less than group velocity

Methodology Details

Task Definition

Investigation of wave propagation characteristics in Bessel-like memory media, including:

  • Dispersion and attenuation analysis of harmonic waves
  • Step response calculation for transient waves
  • Determination of phase velocity and group velocity

Model Architecture

1. Fundamental Wave Equation

Considering the special case of ν=0: ∂2Y∂t2−c2[1−Φ(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

where: Φ~(s)=2sτI1(sτ)I0(sτ)\tilde{\Phi}(s) = \frac{2}{\sqrt{s\tau}}\frac{I_1(\sqrt{s\tau})}{I_0(\sqrt{s\tau})}

2. Dispersion Analysis Method

Employing the standard harmonic wave solution form: Y(t,x)=Aei(ωt−kx)Y(t,x) = A e^{i(\omega t - kx)}

yielding the dispersion relation: (iω)2+k2c2[1−Φ^(ω)]=0(i\omega)^2 + k^2 c^2[1 - \hat{\Phi}(\omega)] = 0

3. SATP Wave Analysis

For spatially attenuated temporally periodic waves, setting:

  • ω∈R\omega \in \mathbb{R} (real frequency)
  • k=κ−iδattk = \kappa - i\delta_{att} (complex wavenumber)

Solving the nonlinear system:

\kappa^2 - \delta_{att}^2 = A(\omega) \\ -2\kappa\delta_{att} = B(\omega) \end{cases}$$ where A(ω) and B(ω) are expressed through Kelvin functions. ### Technical Innovations #### 1. Kelvin Function Representation Converting the ratio of complex Bessel functions into combinations of real Kelvin functions, avoiding numerical instability from complex arithmetic. #### 2. Talbot Method Application Employing the Talbot method for Laplace inverse transform: $$Y(t,x) = \frac{1}{2\pi i}\int_C e^{st}\tilde{Y}(s,x)ds$$ Converting the Bromwich integral into an integral along the real axis through contour deformation, enhancing numerical stability. #### 3. Full-Domain Numerical Solution Overcoming the limitation of traditional methods that are only valid near the wavefront, achieving accurate computation throughout the entire spatiotemporal domain. ## Experimental Setup ### Initial-Boundary Value Problem Considering the step response problem on a semi-infinite domain: $$\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}$$ ### Numerical Implementation - Symbolic computation using Wolfram Mathematica® - Finite difference method for numerical differentiation of group velocity - High-precision numerical integration via the Talbot method ### Parameter Settings Using dimensionless parameters: - Frequency parameter: $\omega\tau$ - Spatial parameter: $x/(c\tau)$ - Time parameter: $(ct-x)/(c\tau)$ ## Experimental Results ### Main Results #### 1. Dispersion Law and Attenuation Factor Figure 1 displays the dimensionless dispersion law $c\tau\kappa$ and attenuation factor $c\tau\delta_{att}$ as functions of $\omega\tau$: - Dispersion law exhibits nonlinear growth - Attenuation factor peaks at intermediate frequencies then decreases #### 2. Phase Velocity and Group Velocity Figure 2 presents a comparison of phase and group velocities: - Phase velocity: $v_p(\omega) = \frac{\omega}{\kappa(\omega)}$ - Group velocity: $v_g(\omega) = \left[\frac{d\kappa(\omega)}{d\omega}\right]^{-1}$ - Numerical results show: $v_p(\omega) < v_g(\omega) < c$, confirming anomalous dispersion characteristics #### 3. Transient Wave Propagation Figure 3 displays step responses at different locations: - Wavefront propagates at velocity c - Wave shape undergoes significant dispersion during propagation - Regions distant from the wavefront exhibit complex oscillatory decay behavior ### Key Findings #### 1. Confirmation of Anomalous Dispersion The system exhibits typical anomalous dispersion characteristics: - $v_p < v_g < c$ - Both velocities approach c in the high-frequency limit #### 2. Effects of Memory The memory term induces: - Wave attenuation and dispersion - Complex transient behavior - Non-exponential wave evolution ## Related Work ### Historical Development 1. **Mainardi et al. (1982-1983)**: First introduction of this class of models in hemodynamics 2. **Giusti & Mainardi (2016)**: Established dynamic viscoelasticity analogy 3. **Colombaro et al. (2017)**: Bessel function-based linear viscoelasticity model ### Method Comparison - **Buchen-Mainardi algorithm**: Semi-analytical, valid only near wavefront - **Present method**: Fully numerical, valid throughout domain - **Asymptotic methods**: High-frequency approximation, strong limitations ## Conclusions and Discussion ### Main Conclusions 1. **Dispersion characteristics**: The system exhibits anomalous dispersion, with phase velocity consistently less than group velocity, both not exceeding wavefront velocity c 2. **Numerical method validity**: The Talbot method successfully overcomes limitations of traditional methods 3. **Memory effects**: Bessel-type memory kernels induce complex wave propagation behavior ### Limitations 1. **Computational complexity**: Transient phase requires numerous nodes, high computational cost 2. **Parameter restrictions**: Only the special case ν=0 is considered 3. **Insufficient theoretical analysis**: Rigorous properties of group velocity are difficult to prove analytically ### Future Directions 1. **Hybrid methods**: Combining semi-analytical methods for wavefront treatment with numerical methods for far-field 2. **Generalization**: Extension to arbitrary ν values 3. **Application expansion**: Verification through applications in concrete physical systems ## In-Depth Evaluation ### Strengths 1. **Methodological innovation**: First application of the Talbot method to wave propagation in Bessel memory media 2. **Theoretical completeness**: Provides a complete framework from dispersion analysis to transient computation 3. **Numerical accuracy**: Overcomes limitations of existing methods, achieving full-domain high-precision computation 4. **Physical insight**: Reveals important physical phenomena such as anomalous dispersion ### Weaknesses 1. **Computational efficiency**: High complexity of transient computation limits practical applications 2. **Parameter range**: Limited to specific memory kernel forms 3. **Experimental validation**: Lacks comparison with actual physical experiments ### Impact 1. **Theoretical contribution**: Provides new numerical analysis tools for wave propagation in memory media 2. **Application value**: Potential applications in biomedical and materials science fields 3. **Methodological significance**: Application of the Talbot method provides reference for similar problems ### Applicable Scenarios 1. **Academic research**: Wave propagation theory, memory material modeling 2. **Engineering applications**: Pipeline flow analysis, vibration control 3. **Numerical computation**: High-precision computation of Laplace inverse transforms ## References The paper cites 23 important references, covering: - Bessel function theory (Abramowitz & Stegun, 1972) - Wave propagation theory (Mainardi, 2022; Whitham, 2011) - Numerical methods (Talbot, 1979; Garrappa et al., 2021) - Application background (hemodynamics, viscoelasticity-related literature) --- This paper makes significant contributions in both theoretical analysis and numerical methods, providing new tools and insights for the study of wave propagation in memory media. Despite certain limitations in computational efficiency and parameter range, its methodological value and physical insights make it an important advance in this field.