- The paper presents a fast numerical scheme that reformulates the viscoelastic constitutive law using an SOE approximation to reduce computational costs.
- It integrates the SOE approximation with Newmark-beta and mixed finite element methods to achieve significant gains in memory and computational efficiency.
- Numerical experiments validate the scheme’s accuracy and superior performance over traditional L1-based methods in long-time wave propagation simulations.
An Efficient Numerical Scheme for Fractional Viscoelastic Wave Propagation Models
Introduction
The paper of viscoelastic materials has gained prominence due to their unique properties, which exhibit a combination of both elastic and viscous behaviors. This paper proposes a fast numerical scheme for fractional viscoelastic models of wave propagation by leveraging an efficient Sum-of-Exponentials (SOE) approximation for the Mittag-Leffler function. The approach significantly improves computational and memory efficiency compared to traditional methods.
Background
Viscoelastic models are characterized by their ability to describe materials that exhibit both elastic (spring-like) and viscous (dashpot-like) behaviors. Classical models like Kelvin-Voigt, Maxwell, and Zener employ integer-order derivatives to capture these behaviors. However, fractional differential operators, which generalize integer-order operators, provide a more precise characterization of viscoelastic materials, particularly those with complex rheological properties.
The fractional viscoelastic model studied in this paper is described by a second-order temporal derivative for wave propagation, coupled with a constitutive equation involving a fractional derivative. The numerical challenge arises from the nonlocal nature of fractional derivatives, which demand high computational and memory resources when using traditional numerical methods such as the L1 scheme.
Contributions and Methodology
The authors' principal contributions include:
- Constitutive Law Reformulation: Conversion of the constitutive equation into an integro-differential form with a Mittag-Leffler function as the convolution kernel.
- SOE Approximation of Mittag-Leffler Function: Usage of Gaussian quadrature to approximate the Mittag-Leffler function efficiently, significantly reducing computational complexity.
- Error Estimation: Provision of a theoretical error estimate for the SOE approximation, which was previously lacking in the literature.
- Fast Numerical Scheme Development: Integration of the SOE approximation with the Newmark-beta scheme and mixed finite element methods to develop a fast scheme for solving the fractional viscoelastic wave propagation model.
By transforming the constitutive equation into a form amenable to SOE approximation, the authors reduce both the memory complexity from O(NsN) to O(NsNexp) and the computational complexity from O(NsN2) to O(NsNexpN). Here, N denotes the number of temporal grid points, Nexp represents the number of terms in the SOE approximation, and Ns captures the complexity related to spatial discretization.
Numerical Results
Numerical experiments are presented to verify the theoretical results and demonstrate the efficiency and accuracy of the proposed scheme. Two key tests include varying the parameters l and q, while keeping the fractional order α and tolerance error ϵ fixed, and varying ϵ while holding α, l, and q constant. The results confirm that the SOE approximation achieves the desired accuracy (O(ϵ)) with significantly fewer terms compared to traditional methods.
Moreover, when applied to a model problem representing wave propagation in a viscoelastic medium, the fast scheme provided computational and memory advantages without compromising accuracy. The results indicate that the fast scheme consistently outperforms the L1-Newmark scheme, particularly for long-time simulations, making it highly suitable for practical applications involving complex viscoelastic materials.
Implications and Future Work
The proposed scheme's implications are notable for both theoretical and practical applications in computational mechanics and material science. The significant reduction in computational and memory resources opens up possibilities for high-fidelity simulations of complex viscoelastic systems over extended periods, which are otherwise infeasible with traditional methods.
Future research may explore extending the fast scheme to three-dimensional problems and other types of fractional differential equations arising in different scientific domains. Additionally, adaptive strategies for optimizing the choice of parameters in the SOE approximation could further enhance computational efficiency.
Conclusion
This paper provides a significant advancement in the numerical simulation of fractional viscoelastic wave propagation models by proposing a fast and efficient scheme based on an SOE approximation of the Mittag-Leffler function. The scheme's theoretical foundations are rigorously validated through numerical experiments, demonstrating its superior performance over traditional methods. This work promises to facilitate more efficient and accurate simulations of complex viscoelastic materials, contributing to advancements in various fields of science and engineering.