Consistent control of energy dissipation in non-spherical particle contact via a structure-preserving formulation

This paper resolves the challenge of controlling energy dissipation in non-spherical particle contact by deriving a structure-preserving damping formulation from first principles that ensures consistent contact-point restitution through a unique damping law aligned with the system's intrinsic translational-rotational coupling.

The Gist

The Big Problem: Why Round Balls Are Easy, But Weird Shapes Are Hard

Imagine you are playing with a bag of marbles. When two round marbles hit each other, they bounce in a very predictable way. You can easily calculate how much energy they lose (dissipation) because they are perfectly symmetrical. It's like a simple spring: you push it, it squishes, and it pops back.

Now, imagine replacing those marbles with irregular rocks, eggs, or jagged pebbles. When these hit a wall or each other, things get messy.

• They don't just bounce up and down; they spin.
• The point where they touch moves as they roll.
• The "weight" they feel at the point of impact changes constantly.

For decades, computer simulations (called the Discrete Element Method, or DEM) have struggled with this. They tried to use the same simple "spring and shock absorber" math used for marbles on these weird shapes. The result? The simulations were inconsistent. Sometimes the rocks bounced too high, sometimes too low, and the energy loss seemed to depend on the angle of the hit rather than the material itself.

The Solution: Changing the Perspective

The author, Y.T. Feng, argues that the old math is structurally broken for non-spherical particles. It's like trying to measure the speed of a spinning top by only looking at its shadow; you miss the rotation.

To fix this, the paper proposes three main ideas:

1. The "Breathing Mass" (The Inflating Balloon)

In a simple marble collision, the "effective mass" (how hard it is to stop the object) is constant. But for a weird shape, the effective mass at the contact point breathes.

• The Analogy: Imagine a gymnast doing a handstand. If they are perfectly straight, their weight is distributed evenly. But if they start to twist or lean, the force on their hands changes instantly.
• In the paper: As an egg hits a wall, the point of contact slides. The lever arm (the distance from the center of the egg to the contact point) changes. This makes the "effective mass" at the contact point shrink and grow during the split-second of impact. The old math assumed this mass was a fixed rock; the new math realizes it's a breathing balloon that changes size every millisecond.

2. The "Coupled Dance" (Translational vs. Rotational)

When a round ball hits a wall, it just bounces back. When a jagged rock hits a wall, the impact doesn't just stop its forward motion; it often spins it up.

• The Analogy: Think of a cue ball hitting a rack of pool balls. If you hit the cue ball dead center, it stops. If you hit it slightly off-center, it spins and the balls scatter.
• In the paper: The energy doesn't just disappear as heat (dissipation). It gets transferred from "moving forward" to "spinning." The old models confused this energy transfer with energy loss. The new model separates them:
  • Material Dissipation: Energy lost to heat/deformation (the "bounciness" of the material).
  • Geometric Transfer: Energy moved from sliding to spinning (the "spin" caused by the shape).

3. The New Rule: Measure the Contact Point, Not the Whole Object

The biggest breakthrough is redefining what "restitution" (bounciness) actually means.

• Old Way: We measured how much total energy the whole particle kept after the crash.
  • Problem: If a rock hits a wall and starts spinning wildly, it might keep 80% of its total energy, even though the material itself was very "dead" (low bounciness). The spin masked the energy loss.
• New Way: We must measure the bounciness only at the specific point of contact.
  • The Metaphor: Imagine a dancer spinning on a stage. If you want to know how slippery the floor is, you shouldn't look at how fast the dancer is spinning (which depends on their pose). You should look at how fast their feet are sliding against the floor.
  • The paper defines a new variable, $e_{cn}$, which is the "bounciness of the contact point." This is the true material property. The total energy bounciness ($e_E$) is just a side effect that changes depending on the angle and shape.

The "Structure-Preserving" Fix

The author didn't just invent a new number; they built a new mathematical framework that respects the physics.

• The Transformation: They used a clever mathematical trick (an "energy-phase transformation") to turn the messy, changing, non-linear collision of a weird shape into a simple, clean, linear spring system.
• The Result: Once transformed, they could apply a damping rule that is perfectly consistent. It ensures that the "contact point" always loses the exact amount of energy you intended, regardless of how the object is spinning or what shape it is.

What This Means for the Real World

1. Better Simulations: Engineers simulating grain silos, pharmaceutical powders, or asteroid impacts can now trust their results. They won't get weird artifacts where rocks bounce too high just because of the angle of impact.
2. Rethinking Experiments: If you drop a weird-shaped object and measure how high it bounces, the height depends on how much it spins. The paper suggests that scientists should stop trying to find a single "bounciness number" for a rock. Instead, they should measure the contact point behavior and let the spin be a natural result of the shape.
3. The "Coupling Offset": The paper explains that the difference between "how bouncy the material is" and "how high the object bounces" is called the coupling offset. It's not a mistake; it's physics. A jagged rock hitting a wall at an angle will naturally convert some of its forward speed into spin, making it bounce higher than a smooth ball of the same material would.

Summary in One Sentence

This paper fixes computer simulations of weird-shaped particles by realizing that their "weight" changes as they hit, they trade forward speed for spin, and to measure true bounciness, we must look at the specific point of contact, not the whole spinning object.

Technical Summary

1. Problem Statement

The control of energy dissipation in Discrete Element Method (DEM) simulations for non-spherical particles remains a fundamental unresolved issue.

• The Limitation of Classical Models: For spherical particles, contact reduces to a one-dimensional normal oscillator with constant effective mass and stiffness. Dissipation is easily calibrated using a constant damping ratio linked to a target coefficient of restitution ($e$).
• The Non-Spherical Challenge: For non-spherical particles, the contact geometry evolves during impact. This leads to:
  1. Configuration-Dependent Inertia: The effective mass projected onto the contact normal is not constant but varies as the contact point migrates (a phenomenon termed "breathing mass").
  2. Intrinsic Coupling: Translational and rotational degrees of freedom are coupled via the lever arm at the contact point. A normal force induces rotation, and rotational motion affects normal dynamics.
• The Consequence: Classical damping formulations (e.g., linear spring-dashpot with constant parameters) are structurally incompatible with these dynamics. They fail to deliver consistent energy loss or prescribed restitution across different impact geometries. Furthermore, the standard definition of the coefficient of restitution (based on total kinetic energy) becomes ambiguous because it conflates material dissipation with geometric energy transfer between translation and rotation.

2. Methodology

The author addresses the problem from first principles using a contact-centric formulation and structure-preserving mechanics.

A. Contact-Centric Dynamics

Instead of solving equations for individual bodies, the dynamics are projected onto the contact point.

• Effective Mass Tensor: The interaction is governed by an effective mass tensor ($M_{eff}$) that maps contact forces to accelerations. For non-spherical particles, this tensor is configuration-dependent and contains off-diagonal terms representing translational-rotational coupling.
• Breathing Mass: The projected normal mass ($m^*_n$) evolves during the collision. For oblique impacts, it starts at a minimum (when the lever arm is largest) and increases as the particle reorients.

B. Structure-Preserving Dissipation

To ensure consistent dissipation, the author employs an energy-phase transformation:

• Mapping to a Harmonic Oscillator: Any monotone conservative contact potential $W(\delta)$ is mapped onto an equivalent linear harmonic oscillator via a coordinate transformation ($x$) and time reparametrization ($\tau$).
• Deriving the Damping Law: In the transformed space, precise control of energy loss (velocity-independent restitution) requires linear damping. Transforming this back to the physical space yields a unique, non-empirical damping law:
  $$C(\delta) \propto \frac{W'(\delta)}{\sqrt{W(\delta)}}$$
  This law ensures that the damping force adapts to the evolving stiffness and geometry.
• Adaptive Implementation: The damping coefficient $\eta(t)$ is calculated instantaneously using the current projected mass $m^*_n(t)$ and effective stiffness $k_{eff}(t)$:
  $$\eta(t) = 2\xi \sqrt{m^*_n(t) k_{eff}(t)}$$

C. Redefinition of Restitution

The paper argues that for coupled systems, the total energy restitution ($e_E$) is not a material property but a geometric outcome. The correct input parameter is the contact-point normal restitution ($e_{cn}$), defined as:
$$e_{cn} = -\frac{v_{exit}^{cn}}{v_{entry}^{cn}}$$
where $v^{cn}$ is the relative velocity of the material points at the contact interface projected onto the normal.

3. Key Contributions

1. Identification of Structural Incompatibility: Demonstrated that classical damping fails for non-spherical particles because it ignores the "breathing mass" and translational-rotational coupling inherent in the effective mass tensor.
2. Structure-Preserving Formulation: Developed a universal damping law derived from energy-phase transformation that preserves the underlying harmonic structure of the contact dynamics, ensuring consistent dissipation control regardless of contact geometry.
3. Theoretical Redefinition of Restitution: Established that $e_{cn}$ (contact-point restitution) is the fundamental material parameter, while $e_E$ (total energy restitution) is a derived quantity dependent on impact angle and particle shape due to energy redistribution.
4. Coupling Offset Concept: Quantified the difference between $e_E$ and $e_{cn}$ as a "coupling offset," explaining why measured restitution often varies with impact angle even for homogeneous materials.

4. Numerical Results

The theory was validated using smooth, single-contact impacts of an ellipsoid against a rigid wall.

• Energy Conservation: In the undamped case, the formulation preserved total energy with errors $< 5 \times 10^{-6}$, confirming the symplectic nature of the integration.
• Control of $e_{cn}$: The adaptive scheme controlled the contact-point restitution to within 1% of the target value across a wide range of aspect ratios (1.0 to 2.5) and impact angles (0° to 60°).
• Prediction of $e_E$: The total energy restitution $e_E$ varied significantly with geometry (e.g., rising from 0.5 to 0.84 as aspect ratio increased), matching the predictions of an impulsive closed-form solution. This confirmed that the "loss" of energy dissipation control in classical models was actually a geometric energy transfer to rotation.
• Approximation Accuracy: The instantaneous application of the damping law (ignoring the time-derivative of $m^*_n$) was found to be highly accurate for typical DEM stiffnesses where relative penetration $\delta_{max}/\rho < 2\%$, resulting in $e_{cn}$ errors $< 1\%$.

5. Significance and Implications

• Resolution of a Long-Standing Problem: Provides a mathematically consistent framework for modeling energy dissipation in non-spherical DEM, moving beyond empirical calibration.
• Guidance for Simulation Practice:
  • DEM users should specify $e_{cn}$ as the input parameter, not $e_E$.
  • Variation in measured restitution with impact angle is a physical consequence of coupling, not necessarily a change in material properties or friction.
• Experimental Interpretation: Offers a mechanistic explanation for experimental observations where restitution appears to depend on impact angle. It suggests that a single fitted value of $e_{cn}$ can explain the angular dependence of measured $e_E$ for smooth, single-contact impacts.
• Stability: The analysis highlights that the "breathing mass" increases the maximum contact frequency, suggesting that standard time-steps based on decoupled estimates may be too large for non-spherical simulations, potentially causing instability.

In conclusion, the paper establishes that dissipation control and restitution definition are inseparable in non-spherical contact. By adopting a structure-preserving formulation and defining restitution at the contact point, one can achieve consistent energy dissipation control that accurately reflects the complex physics of coupled translational-rotational impacts.