🔢 mathematics

Wave-front tracking for a quasi-linear scalar conservation law with hysteresis II: the case of Preisach

This paper extends the wave-front tracking method to solve the Cauchy problem for a quasi-linear scalar conservation law with hysteresis, specifically addressing the more complex and versatile Preisach operator by analyzing the associated Riemann problem and establishing solution uniqueness via an entropy-like condition.

Original authors: Fabio Bagagiolo, Stefan Moreti

Published 2026-01-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Fabio Bagagiolo, Stefan Moreti

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 are trying to predict how a crowd of people moves through a hallway. In a normal hallway, people just walk forward, and if they bump into each other, they might slow down or speed up, but their behavior is usually straightforward: they react to what's happening right now.

But what if these people have memory? What if their decision to speed up or slow down depends not just on the current situation, but on whether they were previously pushed forward or held back? This is the world of hysteresis. It's like a thermostat that doesn't just turn on when the room gets cold, but remembers how long it's been cold and how hot it got before.

This paper is about solving a complex mathematical puzzle involving a "hallway" (a conservation law) where the people (the variables) have this specific type of memory called the Preisach model.

Here is a breakdown of what the authors did, using simple analogies:

1. The Problem: A Hallway with a "Memory" Wall

The authors are studying an equation that describes how a quantity (let's call it $u$ , like traffic density) moves. But there's a twist: the movement is influenced by a "memory term" ( $w$ ).

The Old Way: In a previous paper, the authors studied a simple memory model called the "Play" operator. Imagine a door that only opens if you push it past a certain point, and stays open until you pull it back past another point. It's a simple "on/off" switch with a delay.
The New Challenge: In this paper, they upgraded the model to the Preisach operator. Think of this not as one door, but as a giant wall made of thousands of tiny, independent switches, each with its own unique sensitivity. Some switches flip at a low push, others at a high push. The total effect is the sum of all these tiny switches flipping on and off.
Why it's hard: Because there are so many switches, the "memory" isn't just a simple on/off state; it's a complex, shifting landscape. The math becomes much messier because the "flow" of the crowd isn't just smooth or jagged; it's a mix of smooth waves and sudden jumps.

2. The Solution: "Wave-Front Tracking"

To solve this, the authors used a method called Wave-Front Tracking.

The Analogy: Imagine you are trying to predict a traffic jam. Instead of trying to calculate the exact speed of every single car (which is impossible), you break the traffic into distinct "waves."
- Some waves are Rarefaction Waves: Think of a traffic jam suddenly clearing up. The cars spread out smoothly, like a wave of water expanding.
- Some waves are Shocks: Think of a sudden stop where cars pile up. This is a sharp, jagged line where density changes instantly.
The Innovation: In the old "Play" model, the math only had to deal with these sharp "Shocks." But with the new "Preisach" model (the wall of thousands of switches), the math is forced to deal with Rarefaction Waves (smooth expansions) as well.
The authors had to invent a new way to track how these smooth waves and sharp shocks interact. It's like learning to juggle not just balls, but also slippery, expanding balloons. They had to prove that even with these complex interactions, the waves don't create an infinite mess; they eventually settle into a predictable pattern.

3. The "Internal Variables" (The Switches in the Wall)

The paper highlights that the Preisach model tracks "internal variables."

The Metaphor: Imagine the wall of switches again. As the input (the crowd) moves back and forth, different switches flip. The "state" of the wall is defined by the shape of the line separating the "ON" switches from the "OFF" switches.
The authors had to treat this "shape" as a moving object in a special mathematical space. They proved that even though this shape is complex, it behaves nicely enough that we can track it without losing control of the math.

4. The Result: Existence and Uniqueness

After all this hard work, the authors proved two main things:

Existence: A solution exists. If you set up this complex system with a specific starting crowd and a specific memory wall, there is definitely a valid way the system will evolve over time. It doesn't break or disappear.
Uniqueness: The solution is unique. There is only one correct way for the system to evolve. Even though the math is complicated, there are no "multiple realities" or ambiguous outcomes. If you know the starting point, you know exactly what happens next.

Summary

The authors took a difficult math problem about systems with memory. They moved from a simple "switch" model to a complex "thousands-of-switches" model. They showed that even though this new model creates a chaotic mix of smooth waves and sharp jumps, we can still track them, predict them, and be certain that there is only one true answer to how the system behaves. They did this by building a sophisticated "tracking system" for these waves, ensuring that the math holds up even when the memory gets complicated.

Technical Summary: Wave-front tracking for a quasi-linear scalar conservation law with hysteresis II: the case of Preisach

Problem Statement
The paper addresses the Cauchy problem for a quasi-linear scalar conservation law incorporating a hysteresis operator, specifically the Preisach operator. The governing equation is given by:
$u_t + w_t + u_x = 0, \quad w(x, t) = F[u(x, \cdot), w_0(x)](t)$
where $u$ is the state variable, $w$ is the output of the hysteresis operator $F$ , and the initial data $(u_0, w_0)$ are of bounded variation (BV). The presence of the hysteresis term $w$ introduces a non-local, rate-independent memory effect, rendering the equation non-standard. While the authors previously analyzed the case where $F$ is the Play operator, this work extends the analysis to the Preisach operator, which is widely regarded as a versatile model for hysteresis in applications involving internal variables.

Methodology
The primary methodology employed is the wave-front tracking technique, adapted to handle the specific complexities introduced by the Preisach hysteresis. The approach involves the following steps:

Riemann Problem Analysis: The authors first solve the Riemann problem for piecewise constant initial data. Unlike the Play operator case (where the flux function is piecewise constant), the Preisach operator yields a flux function $f(u)$ that is piecewise continuous and non-constant. Consequently, the solution to the Riemann problem involves not only shocks but also combinations of rarefaction waves. The analysis utilizes the geometric representation of the Preisach operator on the "Preisach plane," tracking the evolution of the switching function (an antimonotone graph) as the input $u$ varies.
Approximation Scheme: To handle general BV initial data, the authors construct a sequence of approximate solutions $(u_n, w_n)$ . The initial data are approximated by piecewise constant functions taking values in a discrete grid $A_n$ . Crucially, the initial hysteresis configurations are approximated by elements of a set $\tilde{S}^{(n)}_L$ , consisting of configurations generated by Lipschitz inputs with graphs composed of segments having vertices on a discrete grid.
Wave Interaction Analysis: The paper analyzes the interactions between waves in the approximate solutions. A key difficulty identified is the interaction between rarefaction waves, which does not occur in the classical scalar case or the Play operator case. The authors prove that the number of discontinuities remains finite and non-increasing, and that interactions between the state variable $u$ and the vertical discontinuities of the hysteresis configuration $z$ are finite.
Metric Space Framework: Since the Preisach operator's state is a function $z: P \to \{-1, 1\}$ (where $P$ is the Preisach plane), the solution $z$ takes values in a metric space of antimonotone graphs. The authors utilize properties of BV functions with values in metric spaces to establish compactness and convergence.
Limit Procedure: By applying Helly's compactness theorem adapted to metric-space-valued functions, the authors extract a convergent subsequence of the approximate solutions. They rigorously pass to the limit to show that the limit pair $(u, w)$ satisfies the weak formulation of the PDE and the hysteresis relationship.
Uniqueness: Uniqueness is established within the class of entropy solutions. An entropy condition is formulated involving a test function and a comparison state, and the authors prove that the constructed solution satisfies this condition and that any two entropy solutions satisfy an $L^1$ -contraction property.

Key Contributions and Results

Existence of Weak Solutions: The paper proves the existence of a weak solution to the Cauchy problem for the conservation law with Preisach hysteresis for initial data in $BV(\mathbb{R}) \times BV(\mathbb{R}; \tilde{S}_L)$ .
Handling Non-Standard Flux: The authors demonstrate that the flux function associated with the Preisach operator is piecewise continuous and non-constant, necessitating a more complex wave-front tracking procedure that accounts for the interaction of rarefaction waves, a feature absent in the Play operator case.
Metric Space Compactness: A significant technical contribution is the rigorous treatment of the hysteresis state as a function with values in a metric space of antimonotone graphs. The authors establish the necessary compactness results to pass to the limit in the hysteresis relationship, addressing the non-local nature of the Preisach operator.
Uniqueness via Entropy: The paper defines an appropriate entropy condition for this system and proves that the wave-front tracking solution is the unique entropy solution. This is achieved by adapting the doubling of variables technique to the hysteresis context.
Relaxation of Hysteresis: The authors employ a relaxation of the hysteresis relationship, rewriting it as a measure-dependent variational integral inequality, which is essential for handling the BV regularity of the solution.

Significance and Claims
The authors claim that this work represents a novelty in the analysis of hyperbolic conservation laws with hysteresis, particularly regarding the passage to the limit in the hysteresis relationship for the Preisach model. While previous works (e.g., Visintin [26]) studied similar equations with different hysteresis operators (like the delayed Relay) or used different approaches, this paper adopts a constructive approach focused on characteristics and their interactions.

The paper explicitly notes that while wave-front tracking has been discussed in other contexts (e.g., Fan [13], Amadori et al. [1, 2]), the specific challenges posed by the Preisach operator—namely the presence of rarefaction wave interactions and the need to work in a metric space of graphs—require a new analysis. The authors do not claim to have discovered new physical applications but rather provide a rigorous mathematical framework for solving the Cauchy problem for this specific class of non-local conservation laws, extending the scope from the Play operator to the more general and versatile Preisach operator. The results are presented as a constructive existence and uniqueness proof for the defined class of initial data.

1. The Problem: A Hallway with a "Memory" Wall

2. The Solution: "Wave-Front Tracking"

3. The "Internal Variables" (The Switches in the Wall)

4. The Result: Existence and Uniqueness

Summary

More like this