🔢 mathematics

A free boundary approach to the quasistatic evolution of debonding models

This paper establishes the existence of energetic solutions for the quasistatic evolution of adhesive membrane debonding by reformulating the problem as a one-phase Bernoulli free boundary problem, which simplifies the analysis by enabling a Minimizing Movements scheme in function spaces while maintaining a coupled algorithm to enforce irreversibility.

Original authors: Eleonora Maggiorelli, Filippo Riva, Edoardo Giovanni Tolotti

Published 2026-06-01

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Eleonora Maggiorelli, Filippo Riva, Edoardo Giovanni Tolotti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a piece of sticky tape (an adhesive membrane) stuck firmly to a table (a substrate). Now, imagine someone slowly pulls up one edge of that tape. As they pull, the tape peels away, creating a growing gap between the tape and the table. This process is called debonding or peeling.

This paper is a mathematical investigation into exactly how that peeling happens when it happens very slowly (so slowly that the system is always in a state of balance, or "quasistatic"). The authors, Eleonora Maggiorelli, Filippo Riva, and Edoardo Giovanni Tolotti, wanted to create a perfect mathematical recipe to predict how the peeled area grows over time.

Here is the breakdown of their work using simple analogies:

The Problem: A Messy Puzzle

In the real world, as you peel tape, the boundary between the "stuck" part and the "peeled" part moves. Mathematically, tracking this moving boundary is like trying to solve a puzzle where the pieces themselves are constantly changing shape and size. Previous attempts to solve this were very complicated, often requiring the mathematicians to look at "fuzzy" shapes or use abstract theories that were hard to apply to real, clean shapes.

The New Approach: A "Free Boundary" Trick

The authors found a clever shortcut. Instead of trying to track the shape of the peeled area directly, they decided to track the displacement (the height) of the tape itself.

The Old Way: Trying to draw the exact outline of the peeled region at every second.
The New Way: Looking at the height of the tape. If the tape is touching the table, its height is zero. If it's peeled, it lifts up. The "peeled area" is simply the region where the tape is above zero.

This turns the problem into a famous mathematical challenge known as the Bernoulli Free Boundary Problem. Think of it like this: instead of chasing a moving shadow, you just watch the object casting the shadow. If you know where the object is, you know where the shadow is.

The Rules of the Game

The authors set up two main rules that the peeling process must follow:

The "Lazy" Rule (Global Stability): At any given moment, the tape wants to settle into the position that uses the least amount of energy. It wants to peel as much as possible without spending too much energy to break the glue. It's like a hiker who wants to reach the summit but will only take a path that doesn't require climbing a mountain if a flat road is available.
The "One-Way Street" Rule (Irreversibility): Once the tape is peeled, it cannot stick back down. The peeled area can only grow, never shrink. This is crucial because real glue doesn't usually re-attach itself perfectly once broken.

The Solution: A Step-by-Step Algorithm

To prove that a solution exists (meaning, that there is a valid way the tape can peel), the authors used a method called Minimizing Movements.

Imagine the peeling process not as a smooth movie, but as a series of still photos taken very quickly.

Step 1: Take a photo of the tape.
Step 2: Ask, "If I pull the tape just a tiny bit more, what is the most energy-efficient way to do it, given that I can't un-peel anything?"
Step 3: Move to the next photo.
Repeat: Do this thousands of times.

The authors showed that if you take these steps small enough, the sequence of photos converges to a smooth, realistic peeling motion. They proved that this motion exists and satisfies all the physical laws of energy and irreversibility.

Why This Matters (According to the Paper)

Simplicity: Their method works with simple, clean "open sets" (regular shapes) rather than complicated, abstract mathematical objects. This makes the math easier to understand and verify.
Connection: They successfully linked the peeling problem to the well-studied "Bernoulli problem," allowing them to use existing mathematical tools to solve a new, difficult problem.
One-Dimensional Proof: They even tested their theory on a simple 1D line (like peeling a strip of tape) and showed that while a solution exists, it might not be unique (there could be multiple ways the tape peels) or perfectly smooth (it might jump suddenly). This highlights the complexity of the real-world phenomenon.

In Summary

The paper provides a new, simpler, and more robust mathematical framework for predicting how adhesive materials peel off surfaces. By changing the perspective from "tracking the shape of the hole" to "tracking the height of the material," they proved that a valid, energy-conserving, one-way peeling process always exists, even in complex scenarios. They didn't claim this fixes tape in real life or predicts medical outcomes; they simply proved that the mathematical rules governing this physical process are sound and solvable.

Technical Summary: A Free Boundary Approach to the Quasistatic Evolution of Debonding Models

Problem Formulation
The paper investigates the mechanical phenomenon of progressive debonding (or peeling) of an adhesive elastic membrane from a substrate. This process is modeled as a quasistatic variational evolution of sets. The system is defined on a fixed open set $\Omega \subset \mathbb{R}^d$ representing the reference configuration of the membrane, with a portion $\Gamma \subset \partial \Omega$ subject to a prescribed time-dependent vertical displacement $w(t)$ . As $w(t)$ evolves, it forces the membrane to detach, creating a debonded region $A(t) \subseteq \Omega$ .

The evolution is governed by two primary components:

Driving Energy ( $E$ ): The linearized elastic energy of the membrane, defined as the minimum Dirichlet energy among functions $v$ that satisfy the boundary condition $w(t)$ on $\Gamma$ and vanish on the adhered region $\Omega \setminus A(t)$ .
Dissipation Distance ( $D$ ): A functional representing the energy required to grow the debonded region from $A$ to $B$ . It is defined as the integral of a toughness function $\kappa$ over the new debonded area $B \setminus A$ , plus an infinite penalty if the growth is not monotonic ( $A \not\subseteq B$ ), thereby enforcing the irreversibility of the debonding process.

The goal is to find an evolution of sets $t \mapsto A(t)$ that satisfies a global stability condition (minimizing total energy at each step) and an energy balance equation (conservation of free energy, dissipated energy, and work done by external forces).

Methodology and Approach
The authors propose a novel approach that reformulates the problem in terms of displacements (scalar functions) rather than directly tracking the evolution of sets. This strategy leverages the connection between the debonding model and the one-phase Bernoulli free boundary problem (also known as the Alt–Caffarelli problem).

Key methodological steps include:

Displacement Reformulation: Instead of minimizing over sets, the authors seek an evolution of functions $t \mapsto u(t)$ that minimizes a functional combining the Dirichlet energy and a penalty term for the positivity set $\{u > 0\}$ . The debonded region is then identified as the union of the positivity sets of the solution over time: $A_u(t) = A_0 \cup \bigcup_{s \in [0,t]} \{u(s) > 0\}$ . This identification ensures the irreversibility constraint is satisfied.
Minimizing Movements Scheme: The existence of solutions is established using the classical Minimizing Movements (or time-discretization) scheme. The authors construct a sequence of discrete-time solutions by iteratively minimizing the Alt–Caffarelli functional.
Compactness and Convergence: To pass to the limit as the time step vanishes, the authors utilize a version of Helly's Theorem for functions of bounded variation in Banach spaces. They extract convergent subsequences for both the characteristic functions of the debonded regions and the displacement fields.
Handling Irreversibility: A crucial technical aspect is the explicit tracking of the debonded region at each discrete time step within the algorithm. This is necessary because, unlike reversible free boundary problems (e.g., droplet evolution), the debonded region cannot shrink, and the positivity set of the minimizer at a single time step does not fully capture the history of the evolution.

Key Contributions and Results

Existence of Energetic Solutions: The paper proves the existence of displacement energetic solutions (DES) and, consequently, shape energetic solutions (SES) for the quasistatic debonding problem.
- Theorem 2.7: Establishes existence under a specific assumption (2.10) regarding the initial configuration and the boundary data, which ensures the existence of a suitable competitor function. In this case, the energy balance holds with a simplified expression for the work of the displacement.
- Theorem 2.9: Removes the restrictive assumption (2.10), proving the existence of a globally stable evolution even for cases where the initial debonded region is empty ( $A_0 = \emptyset$ ). However, in this general case, the energy balance is initially obtained as an inequality.
- Proposition 2.10: Demonstrates that under stronger regularity assumptions on the boundary data ( $w \in W^{2,p}$ ) and a uniform global Lipschitz bound on the approximating sequence, the energy balance inequality becomes an equality, yielding a full DES.
Simplification of the Framework: By shifting the focus from sets to functions, the authors avoid the technical complexities associated with the theory of capacities and quasi-open sets used in previous literature (e.g., [7]). The new approach allows for a more accessible proof of existence within the natural class of open sets.
Equivalence of Formulations: The paper rigorously establishes the equivalence between the set-based formulation (SES) and the function-based formulation (DES) under appropriate conditions, showing that the displacement formulation is a robust tool for analyzing the problem.
One-Dimensional Analysis: The authors provide an explicit one-dimensional framework (Section 2.1) to illustrate the model. This section includes examples demonstrating that uniqueness and time-regularity of solutions are not guaranteed in general, highlighting the complexity of the problem.

Significance and Claims
The paper claims to improve and simplify previous results in the literature regarding quasistatic debonding. While prior works (such as [7]) provided existence results for capacitary measures, they often lacked complete information on the energy balance or required strong geometric assumptions to recover shape information.

The authors' approach is significant because:

It provides a complete existence result for energetic solutions within the class of open sets, avoiding the need for the more abstract theory of capacitary measures.
It offers a simplified and more accessible mathematical framework by utilizing the well-developed theory of the one-phase Bernoulli free boundary problem.
It explicitly addresses the irreversibility of the phenomenon through a coupled algorithm that tracks the history of the debonded region, a feature that distinguishes it from reversible free boundary models like droplet evolution.

The authors remain modest regarding regularity and uniqueness, noting that these properties are not expected to hold generally, as evidenced by the one-dimensional counterexamples provided. The work focuses primarily on establishing the existence of solutions that satisfy the fundamental energetic criteria of stability and energy balance.