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

This paper presents a novel, computationally efficient integral solution for the mechanical response of a viscoelastic Kelvin-Voigt medium to pulse excitations, offering a simpler alternative to complex inverse Laplace transform methods and enabling straightforward asymptotic analysis.

The Gist

Imagine you are standing in a very long, thick hallway made of a strange, stretchy material—like a giant block of Jell-O mixed with honey. This material is what scientists call a Kelvin-Voigt medium. It's not just a solid (like a steel beam) and not just a liquid (like water); it's a bit of both. It has elasticity (it wants to snap back to shape) and viscosity (it resists moving, like thick syrup).

Now, imagine you give one end of this hallway a sharp push (a "pulse") or a steady hold (a "step"). The question is: How does that movement travel down the hallway? Does it zip through instantly? Does it fade away? Does it wiggle?

This paper is about solving that exact puzzle, but with a major upgrade to the math tools we use to do it.

The Old Way: The "Black Box" Problem

For decades, scientists have tried to predict how this "stretchy-honey" hallway reacts. They used a powerful mathematical tool called the Laplace Transform. Think of this tool as a translator that turns a messy, moving problem into a static, easy-to-solve algebra problem.

However, once you solve the algebra, you have to translate it back to the real world. This is called the Inverse Laplace Transform.

• The Problem: Translating it back is like trying to reverse-engineer a complex cake just by tasting a single crumb. It requires doing difficult calculations in a "complex plane" (a weird, imaginary mathematical space).
• The Result: To get the answer, you usually have to use a computer to do a massive, slow, and error-prone numerical integration. It's like trying to find a needle in a haystack by looking at every single piece of hay one by one.

The New Way: The "Direct Map"

The authors of this paper (González-Santander, Mainardi, and Mentrelli) said, "Let's find a shortcut." They derived a new, direct formula (an integral solution) that skips the difficult "translation" step entirely.

Instead of using the complex, imaginary math, they found a way to write the answer as a single, clean integral.

• The Analogy: Imagine you need to drive from City A to City B. The old method was to take a detour through a foggy, winding mountain pass (the complex plane) just to get there. The new method is a straight, paved highway. You can see the destination clearly, and the car (the computer) gets there much faster and with fewer mistakes.

Two Types of Pushes

The paper looks at two specific scenarios, which are like two different ways to test the hallway:

1. The "Delta Pulse" (The Sharp Tap):

   • Imagine flicking the end of the hallway with your finger. It's a quick, sharp tap that happens in a split second.
   • The authors found a new formula that tells you exactly how that sharp tap ripples out, fades, and spreads over time. They also figured out what happens if you look at the hallway immediately after the tap (very small time) or long after (very large time).

2. The "Step Pulse" (The Steady Push):

   • Imagine pushing the end of the hallway and holding it there. It's a steady, constant force.
   • They found a new formula for this too. It shows how the hallway slowly stretches out and settles into a new shape. Again, they figured out the "beginning" and "end" behaviors of this stretch.

Why Does This Matter?

You might ask, "Who cares about a stretchy hallway?"

• Earthquakes (Seismology): The Earth's crust isn't perfectly rigid; it acts a bit like this Kelvin-Voigt material. When an earthquake happens, the energy travels through the ground. Understanding exactly how these waves move helps scientists predict how strong the shaking will be in different places.
• Medical Imaging: Some tissues in the human body are viscoelastic. Understanding these waves helps improve ultrasound and other imaging techniques.
• Engineering: When designing materials for cars or planes that need to absorb shock without breaking, engineers need to know exactly how vibrations travel through them.

The Bottom Line

This paper is a mathematical toolkit upgrade.

• Before: Scientists had to use a slow, clunky, and sometimes inaccurate method to predict how waves move through stretchy materials.
• Now: They have a faster, simpler, and more accurate formula. It's like upgrading from a flip phone to a smartphone. The result is the same (you can still make calls), but the new way is so much more efficient that you can do complex things in seconds that used to take hours.

The authors also double-checked their work using powerful computer software (Mathematica) and proved that their new "highway" leads to the exact same destination as the old "mountain pass," just with a much smoother ride.

Technical Summary

Here is a detailed technical summary of the paper "Pulse waves in the viscoelastic Kelvin-Voigt model: a revisited approach" by González Santander, Mainardi, and Mentrelli.

1. Problem Statement

The paper addresses the fundamental problem of determining the mechanical response $r(x, t)$ of an initially quiescent, semi-infinite, homogeneous medium subjected to a pulse applied at the origin ($x=0$). The medium is modeled using the Kelvin-Voigt constitutive law, a specific type of linear viscoelasticity where the material exhibits both instantaneous elasticity and viscous damping.

• Context: This problem is critical in seismology and the study of transient waves in dissipative media.
• The Challenge: While the problem is well-known, standard solutions typically rely on the inverse Laplace transform, which often requires numerical integration in the complex plane. This approach can be computationally expensive and numerically unstable for certain parameter ranges.
• Goal: The authors aim to derive novel integral-form solutions for both step-pulse (Heaviside) and delta-pulse excitations that avoid complex-plane integration, offering simpler, more computationally efficient expressions and facilitating asymptotic analysis.

2. Methodology

The authors employ a rigorous analytical approach combining Laplace transform theory, special functions, and integral representations.

• Governing Equations:
  • The Kelvin-Voigt constitutive relation is $\sigma = G_e [\epsilon + t_\epsilon \frac{\partial \epsilon}{\partial t}]$, where $G_e$ is the equilibrium modulus and $t_\epsilon$ is the retardation time.
  • The wave equation in the Laplace domain is derived as:
    $$\tilde{r}(x, s) = \tilde{r}_0(s) \exp\left( -\frac{s}{c'} \sqrt{1 + t_\epsilon s} \, x \right)$$
    where $c' = \sqrt{G_e/\rho}$ is a characteristic velocity.
• Dimensionless Variables:
  To simplify the analysis, the authors introduce dimensionless variables:
  • $\xi = \frac{x}{c' t_\epsilon}$ (dimensionless distance)
  • $\tau = \frac{t}{t_\epsilon}$ (dimensionless time)
• Derivation Strategy:
  1. Inverse Laplace Transform: Instead of direct numerical inversion, the authors utilize the convolution theorem and specific properties of the Laplace transform to express the solution as a real integral.
  2. Special Functions: The derivation involves expanding exponential terms and utilizing Hermite functions ($H_\nu$) and the confluent hypergeometric function (${}_0F_1$).
  3. Kernel Construction: A fundamental kernel function $g(\xi, \tau)$ is derived, representing the impulse response. The general solution for any pulse $r_0(t)$ is then expressed as a convolution of this kernel with the input pulse.
  4. Asymptotic Analysis: The authors perform asymptotic expansions for the limits $\tau \to 0, \infty$ and $\xi \to 0, \infty$ to derive simple approximate formulas.

3. Key Contributions

The paper makes several significant theoretical and practical contributions:

• Novel Integral Representations:
  • Delta-Pulse: An explicit integral formula involving the hypergeometric function ${}_0F_1$ is derived (Eq. 2.48).
  • Step-Pulse: A corresponding integral formula is derived (Eq. 2.28) by integrating the delta-pulse solution or via integration by parts.
  • These formulas are defined over the real line, avoiding the need for complex contour integration.
• Computational Efficiency:
  • The authors compare their new integral forms against existing solutions in the literature (e.g., Hanin [1957], Dozio [1990], Morrison [1956]).
  • They demonstrate that their formulation is computationally superior, particularly for large times ($\tau \gg 1$), where previous formulations (like Morrison's) become numerically difficult to evaluate.
  • They note that Dozio's solution appears incorrect based on their numerical experiments.
• Comprehensive Asymptotic Formulas:
  The paper provides simple, closed-form asymptotic approximations for the response in four distinct regimes:
  1. Short time / Large distance ($\tau \to 0$ or $\xi \to \infty$): The solution behaves like a diffusion wave (error function form).
  2. Long time ($\tau \to \infty$): The solution approaches equilibrium with corrections expressed via generalized hypergeometric functions (${}_0F_2$).
  3. Small distance ($\xi \to 0$): Taylor series expansions are provided.
  4. Delta-pulse specific asymptotics: Explicit formulas for the impulse response in all limits.

4. Results

• Numerical Validation:
  • The authors used MATHEMATICA to verify the equivalence between their new integral solutions and the exact inverse Laplace transform solutions.
  • Step-Pulse: The maximum discrepancy between the integral formula and the inverse Laplace transform was found to be $\approx 8.3 \times 10^{-10}$.
  • Delta-Pulse: The maximum discrepancy was $\approx 5.1 \times 10^{-9}$.
  • These errors are within the expected numerical precision, confirming the mathematical correctness of the new integral forms.
• Asymptotic Accuracy:
  • The derived asymptotic formulas (e.g., Eq. 2.30, 2.40, 2.49) were numerically checked and shown to accurately approximate the exact solutions in their respective limits.
• Comparison with Literature:
  • The new step-pulse solution is shown to be more robust for numerical evaluation at large $\tau$ compared to Morrison's integral.
  • The new delta-pulse solution is more efficient than Hanin's formulation.

5. Significance

• Seismology and Geophysics: The Kelvin-Voigt model is widely used to describe seismic wave attenuation. The new solutions provide seismologists with more efficient tools to model transient wave propagation without the computational overhead of complex-plane integration.
• Mathematical Physics: The paper resolves a long-standing issue in the literature regarding the computational efficiency of solving viscoelastic wave equations. It provides a unified framework for handling both step and delta excitations using real integrals.
• Practical Application: By providing explicit asymptotic formulas, the paper allows researchers to quickly estimate wave behavior in limiting cases (very short/long times or distances) without running full numerical simulations.
• Reproducibility: The authors provide a link to a MATHEMATICA notebook containing all derivations and numerical checks, ensuring the results are transparent and reproducible.

In summary, this paper revisits a classic problem in viscoelasticity, successfully replacing complex numerical inversion methods with elegant, computationally efficient real-integral solutions and providing a complete set of asymptotic approximations for the Kelvin-Voigt model.