🔢 mathematics

Weak and dissipative solutions for the Hasegawa-Mima equation

This paper establishes the existence of dissipative solutions for the Hasegawa-Mima equation in its Euler-like velocity form for any $L^2$ divergence-free initial condition on both the 2D torus and bounded $\mathcal{C}^1$ domains, by adapting the framework originally developed by Lions for the Euler equations.

Original authors: Michele Gorini

Published 2026-07-02✓ Author reviewed ⓘ

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Michele Gorini

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 by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Chaotic Dance in a Magnetic Field

Imagine a giant, invisible dance floor inside a nuclear fusion reactor (like a tokamak). On this floor, charged particles (plasma) are swirling around. Because there is a massive magnetic field holding them in place, they don't just move randomly; they dance in a very specific, complex pattern.

The Hasegawa-Mima equation is the mathematical rulebook that describes how this dance evolves. It tells us how the "potential" (the energy or pressure of the dance) changes over time as the particles swirl and bump into each other.

However, there's a problem. In the real world, these dances can get incredibly chaotic. The particles might swirl so violently that the math breaks down. The smooth, predictable curves we expect turn into jagged, unpredictable spikes. In mathematical terms, the "perfect" solution stops existing, or we can't prove it exists using standard tools.

The Problem: When the Math Gets "Sticky"

The author, Michele Gorini, is tackling a specific headache: What happens when the dance gets too messy for the standard rules?

In physics, we usually look for "weak solutions." Think of a weak solution as a "good enough" answer. It's like looking at a blurry photo of a fast-moving car. You can't see the license plate (the fine details), but you can tell the car is moving and in what direction.

The trouble with the Hasegawa-Mima equation is that when you try to create these "blurry" solutions by smoothing out the chaos, the math sometimes fails to settle down. The nonlinear part of the equation (the part where the particles interact with each other) creates "oscillation defects."

The Analogy: Imagine trying to blend a smoothie. If you have a few hard chunks of ice (the chaos), the blender might just spin them around without crushing them, or the mixture might splatter everywhere. Standard math says, "If you can't crush the ice perfectly, the smoothie doesn't exist." Gorini wants to say, "Actually, the smoothie exists, it's just a bit lumpy, and we can still describe its overall texture."

The Solution: "Dissipative" Solutions

To solve this, Gorini adapts a concept from a famous mathematician named Pierre-Louis Lions, who worked on the Euler equations (which describe how air or water flows). He introduces the idea of "Dissipative Solutions."

Here is the metaphor for a dissipative solution:

Imagine you are watching a chaotic crowd of people running through a hallway.

The Perfect Solution: You track every single person's exact path. This is impossible if the crowd is too dense.
The Weak Solution: You try to guess the average flow. But sometimes, the crowd creates a "traffic jam" that the average flow can't predict.
The Dissipative Solution: Instead of trying to track the crowd, you compare the chaotic crowd to a perfectly organized, smooth-flowing guide (a "test flow").

Gorini's method asks: "How much does the chaotic crowd deviate from this perfect guide?"

He proves that even if the crowd is messy, the energy of the difference between the real crowd and the perfect guide stays under control. It's like saying, "Even if the dancers are tripping over each other, the total amount of 'tripping' energy won't explode out of control."

If a perfect, smooth solution does exist, the dissipative solution will match it exactly. But if the smooth solution breaks down, the dissipative solution is the "best possible" backup that still obeys the laws of physics (specifically, energy conservation).

How the Paper Works (The Steps)

Changing the View: The paper starts by rewriting the equation. Instead of looking at the "potential" (the height of the energy waves), it looks at the "velocity" (how fast the particles are moving). This is like switching from looking at the waves on the ocean to looking at the speed of the water current. It makes the math look more like the famous Navier-Stokes equations used for weather and water.
The "Viscous" Trick: To prove these solutions exist, the author first adds a little bit of "syrup" (viscosity) to the equation. This makes the fluid thicker and easier to handle mathematically, smoothing out the chaos. He proves that with syrup, the solution exists and behaves well.
Removing the Syrup: Then, he slowly removes the syrup (letting it go to zero). He shows that even as the fluid becomes thin and chaotic again, the "dissipative" behavior remains stable. The solution doesn't vanish; it just becomes a "dissipative" one.
The Result: He proves that for any starting condition (any initial arrangement of the particles) that has a finite amount of energy, a dissipative solution exists forever.

Why "Dissipative"?

The term "dissipative" might sound like energy is being lost, but in this context, it's about stability.

Analogy: Think of a spinning top. If it's perfect, it spins forever. If it's slightly wobbly, it eventually falls. A "dissipative" solution is like a top that is allowed to wobble and lose a tiny bit of energy to friction, but we have a mathematical guarantee that it won't suddenly explode or fly off the table. It stays within the rules of the game.

Summary of Claims

The Goal: To prove that the Hasegawa-Mima equation (which models plasma in fusion reactors) always has a valid solution, even when the plasma gets extremely chaotic.
The Method: Using a "relative energy" comparison against smooth, ideal flows to define a new type of solution called a "dissipative solution."
The Result: The paper proves that these solutions exist for all time, starting from any reasonable initial state.
The Safety Net: If a perfect, smooth solution exists, this new "dissipative" solution is exactly the same as it. If the smooth one breaks, the dissipative one takes over without breaking the laws of physics.

In short, Gorini has built a safety net for the math describing plasma turbulence. Even when the particles go wild, the math still holds up, ensuring we can always predict the "big picture" of the dance, even if we can't track every single step.

Technical Summary: Weak and Dissipative Solutions for the Hasegawa-Mima Equation

Problem Statement
The Hasegawa-Mima (HM) equation is a fundamental model for drift-wave turbulence in strongly magnetized plasmas, describing the nonlinear evolution of electrostatic potential fluctuations. While the equation is closely related to the two-dimensional Euler equations and quasi-geostrophic models, its nonlinear structure presents significant analytical challenges, particularly regarding the compactness of approximate solutions. In the context of the incompressible Euler equations, weak limits of approximate solutions (e.g., via vanishing viscosity or spectral truncation) may fail to satisfy the equation in the classical weak sense due to oscillation defects in quadratic terms. This paper addresses the lack of a robust existence theory for the HM equation that accommodates these defects, specifically for arbitrary $L^2$ divergence-free initial data in both periodic ( $T^2$ ) and bounded ( $D \subset \mathbb{R}^2$ ) domains.

Methodology
The paper adapts the framework of "dissipative solutions," originally introduced by P.-L. Lions for the Euler equations, to the HM equation. The methodology proceeds through the following steps:

Formulation: The author derives the "velocity form" of the HM equation, $\partial_t(u - \Delta^{-1}u) + (u \cdot \nabla)u - u^\perp \log n_0 = 0$ , where $u = \nabla^\perp \phi$ is the drift velocity and $n_0$ is a time-independent particle density profile. This form is shown to be equivalent to the standard potential form under incompressibility constraints.
Analysis of Weak Solutions: The paper demonstrates that standard definitions of weak solutions (requiring $L^2$ or $H^1$ regularity) are unstable under weak convergence due to the nonlinearity. Specifically, the quadratic term $(u \cdot \nabla)u$ prevents the passage to the limit in standard weak formulations without additional regularity assumptions.
Dissipative Solution Framework: To overcome this, the author introduces the notion of dissipative solutions. These are defined not by satisfying the equation in a distributional sense, but by satisfying a relative energy inequality with respect to arbitrary smooth test flows.
- For the potential form, the relative energy is defined via the $H^1$ norm of the difference between the solution $\phi$ and a smooth test function $\psi$ .
- For the velocity form, the relative energy is defined via a norm involving the $L^2$ norm and a fractional Laplacian term, $\|u - v\|_{H_{\Delta L}}^2 = \|u - v\|_{L^2}^2 + \|(-\Delta)^{-1/2}(u - v)\|_{L^2}^2$ . Note that the Laplacian term in this norm is equivalent to the $H^{-1}$ norm.
Viscous Regularization and Compactness: The existence proof relies on the viscous regularization of the velocity-form HM equation (analogous to the Navier-Stokes equations). The author proves a Leray-like result for this viscous system, establishing the existence of global weak solutions satisfying an energy inequality.
Limit Passage: Using uniform energy estimates and compactness arguments (specifically the Ascoli-Arzelà theorem applied to mollified sequences), the author extracts a subsequence of viscous solutions that converges weakly. The limit is shown to satisfy the relative energy inequality, thereby constituting a dissipative solution.

Key Contributions and Results

Definition of Dissipative Solutions: The paper provides rigorous definitions for dissipative solutions for both the potential form (Definition 4) and the velocity form (Definition 5) of the HM equation.
Equivalence and Uniqueness: It is proven that dissipative solutions of the velocity form correspond to those of the potential form (Proposition 1). Furthermore, the framework satisfies a weak-strong uniqueness property (Proposition 2): if a classical (smooth) solution exists, it is unique within the class of dissipative solutions.
Global Existence Theorem: The main result (Theorem 2) establishes the global-in-time existence of dissipative solutions for the velocity-form HM equation for any divergence-free initial data $u_0 \in L^2(D)$ , assuming $\log n_0 \in L^{2+\epsilon}$ .
Leray-like Result for HM: A specific contribution is the proof (Theorem 1) that the viscous form of the HM equation admits global weak solutions satisfying an energy inequality, analogous to Leray's result for the Navier-Stokes equations. This serves as the foundational step for the existence proof of dissipative solutions.
Domain Generality: The analysis is extended to bounded domains $D \subset \mathbb{R}^2$ with $C^1$ boundaries, addressing technicalities regarding Sobolev embeddings, the spectral properties of the Laplacian with Dirichlet boundary conditions, and the definition of fractional Laplacians.

Significance
The paper claims that the dissipative formulation provides a natural and stable setting for studying the Hasegawa-Mima equation, particularly in regimes where small-scale oscillations may develop and classical weak convergence methods fail to identify a limiting equation without defect measures. By adapting Lions' framework, the work bridges the gap between the analytical theory of fluid equations (Euler/Navier-Stokes) and plasma physics models, offering a robust existence theory that does not rely on the conservation of energy in the classical sense but rather on a relative energy bound. This approach allows for the study of singular limits and approximation schemes where standard weak solutions are insufficient.