Radial Viscous Shallow Water Equations

• Radially symmetric viscous shallow water equations are compressible fluid models with radial symmetry, incorporating density-dependent viscosity and BD entropy structures.
• The system is derived from the compressible isentropic Navier–Stokes equations with variable viscosity, leading to effective simplifications in both analysis and numerical implementation.
• Recent results confirm global existence and instantaneous shock regularization, demonstrating the robustness of the model even with large, discontinuous initial data.

The radially symmetric viscous shallow water equations describe the evolution of compressible, viscous fluid flows with surface layers, subject to the closure constraint of radial (spherical or cylindrical) symmetry. These systems are derived from the compressible isentropic Navier–Stokes equations with density-dependent viscosities, representing the endpoint case of the Bresch–Desjardins (BD) entropy structure for the specific scaling $\mu(\rho)=\rho$, $\lambda(\rho)=0$. Recent research has established the global existence and remarkable regularization properties of solutions to these equations in both two- and three-dimensional radially symmetric settings, including the persistence and smoothing of initial discontinuities and the global regularity for arbitrarily large initial data (Haspot, 2019, Huanga et al., 17 Dec 2025).

1. Mathematical Formulation and Symmetry Reduction

The global dynamics are governed by the isentropic compressible Navier–Stokes system with variable viscosity: $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$ with pressure law $p(\rho) = \rho^\gamma$ ($\gamma > 1$), viscosity coefficients $\mu(\rho)=\rho$, $\lambda(\rho)=0$.

Imposing radial symmetry,

$\rho(x,t) = H(r,t), \quad \bu(x,t) = u(r,t)\frac{x}{r},\quad r=|x|,$

transforms the system to one spatial dimension (in $r$). For $N=2$ (planar/radial) and $\lambda(\rho)=0$0 (spherical) symmetries, the equations respectively become: $\lambda(\rho)=0$1

$\lambda(\rho)=0$2

This reduction is fundamental: the vorticity $\lambda(\rho)=0$3 vanishes, and the system is closed under scalar or radial vector fields, simplifying both the analysis and numerical implementation.

2. Entropy Structures and Effective Velocity

Key analytical control mechanisms arise from the total energy and BD (Bresch–Desjardins) entropy. The effective velocity is introduced as

$\lambda(\rho)=0$4

For radially symmetric flows with $\lambda(\rho)=0$5, the BD entropy becomes

$\lambda(\rho)=0$6

and dissipates with time,

$\lambda(\rho)=0$7

Both classical energy and BD entropy facilitate global estimates and regularization effects unique to the variable viscosity case (Haspot, 2019, Huanga et al., 17 Dec 2025). The case $\lambda(\rho)=0$8 in the viscous law $\lambda(\rho)=0$9 is a critical endpoint: it allows for strong entropy dissipation without further structure on the bulk viscosity.

3. Global Existence and Regularity Results

Extensive recent advances have established the global existence of weak and classical solutions of the radially symmetric viscous shallow water system, even for arbitrarily large, possibly discontinuous initial data:

• For initial data $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$0, $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$1, $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$2, there exists a unique global classical solution $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$3 with $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$4 $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$5 (Huanga et al., 17 Dec 2025).
• Global weak solutions exist when the initial momentum $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$6 and the initial density is merely bounded, potentially allowing for initial shocks (Haspot, 2019).
• For well-prepared, sufficiently regular radially symmetric initial data, strong solutions persist globally in time without smallness conditions (Haspot, 2019, Huanga et al., 17 Dec 2025).

A significant feature is the instantaneous Lipschitz regularization of the density: for any $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$7, $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$8. This implies that any initial shock in the density is smoothed out immediately, an effect not present in the constant-viscosity case.

4. Analytical Techniques and Functional Framework

The proofs invoke several advanced analytical tools:

• Parabolic $\begin{aligned} &\partial_t\rho + \div(\rho \bu) = 0\ &\partial_t(\rho\bu) + \div(\rho\bu\otimes\bu) + \nabla p(\rho) = \div(\mu(\rho)\nabla\bu) + \nabla(\lambda(\rho)\div\bu) \end{aligned}$9 smoothing: The continuity equation with an effective drift and the parabolic regularization structure yield instantaneous gain of spatial regularity in density.
• Entropy inequalities: Energy and BD entropy yield uniform-in-time bounds for relevant Sobolev and Orlicz norms.
• Koch–Tataru framework: For weak solutions, the initial momentum is placed in $p(\rho) = \rho^\gamma$0 and controlled in the corresponding $p(\rho) = \rho^\gamma$1 space using heat-kernel energy estimates and fixed-point iteration.
• Weighted Sobolev and interpolation inequalities: In 2D, crucial radial-embedding theorems of the form $p(\rho) = \rho^\gamma$2 for all $p(\rho) = \rho^\gamma$3 enable $p(\rho) = \rho^\gamma$4 control.
• Lagrangian coordinates: Transformation to mass coordinates provides sharp lower bounds on the fluid depth $p(\rho) = \rho^\gamma$5 and enables control near vacuum.

The radial symmetry plays a critical role, simplifying drift terms and making the maximum principle and certain compactness arguments effective for large data.

5. Central Estimates and Their Implications

The analysis is underpinned by a set of robust a priori inequalities:

• Uniform $p(\rho) = \rho^\gamma$6 bounds for density via maximum principles and control of effective velocity divergence.
• Upper bounds for energy and entropy with time-integrated dissipation.
• Interpolation estimates (e.g., for $p(\rho) = \rho^\gamma$7 norms with $p(\rho) = \rho^\gamma$8-weighted $p(\rho) = \rho^\gamma$9 derivatives in radial coordinates).
• Lower bounds for positivity of density via Lagrangian arguments.

These allow continuation of solutions to arbitrary times, precluding finite-time breakdown from vacuum, blow-up of $\gamma > 1$0, or loss of integrability of the momentum.

Estimate	Main Ingredient	Role in Analysis
$\gamma > 1$1 density	Parabolic max principle	Prevents blow-up, ensures positivity
$\gamma > 1$2 regularity	De Giorgi/Schauder, entropy	Smoothing of initial shocks
$\gamma > 1$3 momentum bounds	Koch–Tataru iteration	Controls infinite-energy initial data
BD entropy	Variable viscosity + structure	Ensures regularity, compactness

6. Comparison with Prior Work

Early theories of compressible Navier–Stokes equations, e.g., the Lions–Feireisl framework, require stronger integrability and regularity (such as $\gamma > 1$4) and constant viscosity, thereby excluding both large discontinuities in density and the BD entropy mechanism (Haspot, 2019). Hoff seminally treated $\gamma > 1$5 initial data but without Lipschitz regularization or BD entropy. Danchin, in the context of critical-Besov spaces, addresses small data but not global well-posedness for large or discontinuous data.

The introduction of variable viscosity and effective velocity, coupled with radially symmetric formulations, enables the treatment of both large and rough data, with strong regularization of shocks—a substantial generalized extension.

7. Open Problems and Future Prospects

Current directions and open problems include:

• Extension to general (non-radial) flows, which requires advanced vorticity control and handling the coupling between vorticity and the effective velocity.
• Analysis of shock front dynamics and discontinuity propagation in higher dimensions.
• Degenerate viscosity laws beyond the $\gamma > 1$6 endpoint and the inclusion and effect of vacuum.
• Asymptotic stability, convergence rates to equilibrium, and the large-time structure of global solutions.

The radially symmetric viscous shallow water system provides a canonical and technically tractable platform for investigating these deep questions. The intersection of BD entropy, variable viscosities, and radial symmetry currently represents the most advanced understanding of regularization, global existence, and shock dynamics in compressible viscid flows (Haspot, 2019, Huanga et al., 17 Dec 2025).