🔢 mathematics

Existence and partial regularity of suitable weak solutions to the 3D Navier-Stokes-Vlasov-Fokker-Planck equations

This paper establishes the existence of a new class of suitable weak solutions to the three-dimensional incompressible Navier-Stokes-Vlasov-Fokker-Planck system by proving strong convergence of the density function via novel a priori quantities and compactness methods, and subsequently characterizes the Hausdorff dimension of the singularity set and the Hölder continuity of the density at regular points.

Original authors: Renjun Duan, Fengqiang Shi, Wendong Wang, Jianbo Yu

Published 2026-02-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Renjun Duan, Fengqiang Shi, Wendong Wang, Jianbo Yu

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 a giant, invisible dance floor where two very different groups of performers are trying to move together without tripping over each other. One group is a thick, sticky fluid (like honey or oil), and the other group is a swarm of tiny particles or bubbles floating inside it.

This paper is about understanding the rules of this dance, specifically when things get messy and chaotic. The authors, Renjun Duan and his team, have written a mathematical "rulebook" for a system called the Navier-Stokes-Vlasov-Fokker-Planck (NSVFP) equations.

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

1. The Dance Floor: The Two-Phase Mixture

Think of the fluid as a crowd of people moving through a hallway (the fluid velocity, $u$ ). Now, imagine thousands of tiny, invisible dust motes floating in that crowd (the particles, described by $f$ ).

The Fluid: It wants to flow smoothly, but it's sticky (viscous).
The Particles: They drift with the wind (the fluid) but also bounce around randomly (like particles in a gas) and get dragged by the fluid.
The Friction: When the fluid moves faster than a particle, or vice versa, they rub against each other. This friction is the "glue" that tries to make them move at the same speed.

2. The Problem: When the Dance Goes Wrong

In the real world, fluids don't always move perfectly smoothly. Sometimes, they get turbulent, swirl into crazy patterns, or even "break" (mathematically speaking, they develop singularities or points of infinite chaos).

For decades, mathematicians have struggled to prove that a solution exists for these complex interactions in 3D space without the math breaking down. The authors of this paper asked: "Can we prove that a 'good enough' solution exists, even if it's not perfectly smooth everywhere?"

3. The Solution: "Suitable Weak Solutions"

Instead of demanding a perfect, smooth dance for every single particle at every single moment (which might be impossible), the authors constructed a class of "Suitable Weak Solutions."

Think of this like a safety net.

A "Strong Solution" is a perfect dancer who never stumbles.
A "Weak Solution" is a dancer who might stumble but keeps the rhythm.
A "Suitable Weak Solution" is a special kind of weak solution that obeys specific energy rules. It's like a dancer who is allowed to stumble, but only in a way that doesn't violate the laws of physics (conservation of energy).

The authors proved that for this specific fluid-particle mix, such a "Suitable Weak Solution" always exists in 3D space, no matter how you start the dance.

4. The Big Challenge: The "Infinite" Space

The paper is set in "whole space," meaning the dance floor is infinite. This makes the math incredibly hard because you have to account for particles that might be infinitely far away or moving infinitely fast.

To solve this, the authors used a clever trick:

The Regularization (The Training Wheels): First, they pretended the dance floor was slightly smaller and the particles were slightly smoother (adding "regularization"). They proved the dance works perfectly here.
The Removal (Taking off the Training Wheels): Then, they slowly removed the training wheels (letting the size go to infinity and the smoothing go away). The hard part was proving that as they removed the training wheels, the dance didn't collapse into chaos. They used advanced mathematical tools (like Tao's decomposition and DiPerna-Lions compactness) to show that the particles and fluid stayed "glued" together in a way that allowed the math to hold up.

5. The Result: Mapping the "Cracks"

The most exciting part of the paper is about Partial Regularity.

The authors didn't just prove the solution exists; they figured out where the solution might be "broken" (singular).

The Analogy: Imagine a cracked windshield. The glass is mostly clear (regular), but there are tiny cracks (singularities).
The Discovery: They proved that while the fluid velocity might have these "cracks," the set of points where the cracks exist is incredibly small. In mathematical terms, the Hausdorff dimension of these bad spots is zero.
- Translation: If you were to shine a light on the "bad spots," they would be so thin and sparse that they would effectively disappear. The fluid is smooth almost everywhere.

Furthermore, they showed that at the points where the fluid is smooth, the particle density is also smooth and behaves nicely (it's Hölder continuous, meaning it changes gradually without sudden jumps).

Summary

In plain English, this paper says:

"We have mathematically proven that a mixture of fluid and particles can always be described by a set of rules, even in a chaotic 3D world. While there might be tiny, isolated spots where the math gets weird, the vast majority of the system behaves smoothly and predictably. We built a new type of 'safety net' solution that captures the energy of the system correctly, ensuring the physics makes sense even when the flow gets turbulent."

This is a foundational result in fluid dynamics, ensuring that our mathematical models for things like pollution dispersion, blood flow with cells, or industrial sprays are built on solid ground.

Technical Summary: Existence and Partial Regularity of Suitable Weak Solutions to the 3D Navier-Stokes-Vlasov-Fokker-Planck Equations

Problem Statement
The paper investigates the three-dimensional incompressible Navier-Stokes-Vlasov-Fokker-Planck (NSVFP) system, which models the dynamics of a two-phase mixture consisting of a viscous incompressible fluid and dispersed particles or bubbles. The system couples the Navier-Stokes equations for the fluid velocity $u$ and pressure $P$ with the Vlasov-Fokker-Planck equation for the particle density distribution function $f$ . The coupling occurs through a frictional force term proportional to the relative velocity $(v-u)$ and a feedback term in the kinetic equation. The primary mathematical challenge addressed is the global existence of "suitable weak solutions" in the whole space $\mathbb{R}^3$ and the subsequent characterization of the partial regularity of these solutions, specifically the Hausdorff dimension of the singular set of the fluid velocity.

Methodology
The authors employ a regularization approach combined with compactness arguments to construct solutions and analyze their properties:

Regularized System: The authors introduce a regularized system (1.3) involving mollifiers ( $\theta_\varepsilon$ ) for the fluid velocity and a cutoff function ( $\gamma_\delta$ ) for the velocity dependence in the friction term. This system admits global strong solutions $(u_{\varepsilon,\delta}, f_{\varepsilon,\delta})$ via a fixed-point argument using the Banach contraction principle.
Uniform Estimates: A critical component of the methodology is establishing uniform estimates independent of the regularization parameters $\varepsilon$ $ε$ and $\delta$ $δ$ .
- Energy Estimates: Global and local energy inequalities are derived for the regularized system.
- Moment Estimates: The authors prove bounds on high-order moments of the distribution function, specifically $|v|^\kappa f$ for $\kappa \ge 3$ .
- Novel Velocity Estimates: To handle the whole-space domain and the non-local frictional force, the authors utilize Tao's Littlewood-Paley decomposition to derive $L^1_t L^p_x$ a priori estimates for the velocity $u$ . This allows them to control the convergence of the density function $f$ and the frictional force term.
Compactness and Convergence: The DiPerna-Lions compactness method is applied to establish the strong convergence of the distribution function $f_{\varepsilon,\delta}$ and its moments in $L^p$ spaces. The Aubin-Lions lemma is used for the strong convergence of the velocity field. Crucially, the authors prove the strong convergence of the entropy term $f \log f$ and weighted entropy terms, which is necessary to pass to the limit in the local energy inequalities.
Verification of Suitable Weak Solution Conditions: The authors rigorously verify that the limit pair $(u, f)$ satisfies all conditions of a "suitable weak solution," including the global energy inequality and two distinct types of local energy inequalities.

Key Contributions and Results

Existence of Suitable Weak Solutions (Theorem 1.2): The paper establishes the global existence of suitable weak solutions for the 3D NSVFP system under initial data conditions requiring finite mass, energy, entropy, and moments up to order $\kappa \ge 3$ . These solutions satisfy the standard distributional formulation and specific energy inequalities.
New Local Energy Inequalities: The authors construct and verify two types of local energy inequalities (Definition 1.1, items (v) and (vi)). These inequalities are essential for partial regularity theory and involve specific test functions and the handling of the non-local coupling terms. The derivation of these inequalities relies heavily on the strong convergence of $f$ , $|v|^\kappa f$ , and $f \log f$ .
Partial Regularity of Fluid Velocity (Theorem 1.5): Under the additional assumption that the initial data possesses moments up to order $15+\kappa$ (where $\kappa > 0$ ), the authors prove that the set of singular points $S_{u, \text{sing}}$ for the fluid velocity $u$ has a parabolic Hausdorff dimension of zero. This result extends the celebrated Caffarelli-Kohn-Nirenberg theory for the pure Navier-Stokes equations to the coupled kinetic-fluid system.
Hölder Continuity of Distribution Function (Theorem 1.6): At any regular point $z_0$ of the fluid velocity $u$ , the authors establish the $\alpha$ -Hölder continuity of the distribution function $f$ . This demonstrates that the regularity of the fluid velocity propagates to the regularity of the particle distribution.

Significance and Claims
The paper claims to present the first result establishing the existence of suitable weak solutions and quantifying the Hausdorff dimension of the singular set for the full 3D NSVFP system. The authors emphasize that the derivation of the local energy inequalities is non-trivial due to the difficulties in establishing the convergence of the density function $f$ and the non-local frictional force term in unbounded domains.

The requirement for high-order moments (up to $15+\kappa$ ) in the partial regularity theorem is described as technical but necessary. It ensures that the frictional force term belongs to $L^q_{t,x}$ with $q > 5/2$ , allowing the coupling term to be treated as a perturbation that does not degrade the parabolic Hausdorff dimension of the singular set. The work bridges the gap between the well-developed partial regularity theory for the Navier-Stokes equations and the more complex kinetic-fluid models, providing a rigorous framework for analyzing singularities in two-phase flows.