On the Profile of Singularity Formation for the Incompressible Hydrostatic Boussinesq System

Basic Information

• Paper ID: 2510.10090
• Title: On the Profile of Singularity Formation for the Incompressible Hydrostatic Boussinesq System
• Authors: S. Ibrahim, Q. Lin, L. Qian, E. S. Titi
• Classification: math.AP (Analysis of PDEs)
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2510.10090v1

Abstract

The primitive equations (PEs) are employed to model large-scale ocean and atmospheric dynamics on planetary scales. While it has been established that smooth solutions of the inviscid primitive equations with constant temperature (isothermal case, also known as the hydrostatic Euler equations) develop stable singularities in finite time, the impact of non-constant temperature on singularity formation remains undetermined. This paper investigates the stability of singularity formation under non-constant temperature, considering two scenarios: temperature without diffusion and temperature dynamics with added vertical diffusion. For both cases, the results demonstrate that temperature variations neither affect the formation of singularities nor the stability of singularities in the velocity field.

Research Background and Motivation

Problem Background

1. Importance of Primitive Equations: The primitive equations serve as fundamental models in ocean and atmospheric science. Through the hydrostatic balance approximation, they neglect vertical acceleration of fluid motion, providing a useful approximation to the Navier-Stokes equations.
2. Limitations of Previous Results: Prior research has primarily focused on the isothermal case (constant temperature), where it has been proven that smooth solutions of the inviscid primitive equations develop stable singularities in finite time.
3. Demands of Practical Applications: In real atmospheric science, temperature gradients play a crucial role in the formation of flow patterns, inducing instabilities and leading to chaotic and turbulent behavior.

Research Motivation

1. Theoretical Completeness: Understanding the impact of non-constant temperature on singularity formation mechanisms is necessary to fill theoretical gaps.
2. Physical Significance: Temperature variations are ubiquitous in geophysical fluids, and understanding their effect on singularity stability has important physical implications.
3. Mathematical Challenges: The introduction of non-constant temperature increases system complexity, requiring novel mathematical techniques for analysis.

Core Contributions

1. Theoretical Breakthrough: First rigorous proof of the impact of non-constant temperature on singularity formation and stability in inviscid primitive equations.
2. Two Temperature Models: Systematic investigation of two distinct temperature dynamics models:
   • Non-diffusive temperature (σ = 0): Temperature governed by transport equation
   • Diffusive temperature (σ = 1): Temperature governed by convection-diffusion equation
3. Stability Results: Proof that temperature variations neither affect the formation of singularities in the velocity field nor their stability.
4. Explicit Singularity Profile: Provision of explicit asymptotic behavior and blowup rates for singularity formation.

Detailed Methodology

Problem Formulation

Investigation of singularity formation in three-dimensional incompressible inviscid primitive equations under non-constant temperature:

$\begin{cases} u_t + uu_X + vu_Y + wu_Z + p_X - \Omega v = 0 \\ v_t + uv_X + vv_Y + wv_Z + p_Y + \Omega u = 0 \\ p_Z + \theta = 0 \\ \theta_t + u\theta_X + v\theta_Y + w\theta_Z - \sigma\theta_{ZZ} = 0 \\ u_X + v_Y + w_Z = 0 \end{cases}$

where $\sigma \in \{0,1\}$ corresponds to non-diffusive and diffusive temperature cases, respectively.

Model Architecture

1. Symmetry Reduction

Through imposition of specific symmetry conditions:

• Solutions are independent of the Y coordinate
• Velocity field possesses specific parity properties
• Coriolis force is set to zero ($\Omega = 0$)

The 3D problem is reduced to a 1D evolution equation: $\begin{cases} a_t - a^2 + (\partial_Z^{-1} a)a_Z + \partial_Z^{-1} c + \int_0^1 (2a^2 - \partial_Z^{-1} c)dZ = 0 \\ c_t - 2ac + (\partial_Z^{-1} a)c_Z - \sigma c_{ZZ} = 0 \end{cases}$

where $a(Z,t) = -u_X(X=0,Z,t)$ and $c(Z,t) = \theta_{XX}(X=0,Z,t)$.

2. Self-Similar Transformation

Introduction of spacetime rescaling: $z = \frac{Z}{\nu(t)}, \quad \frac{ds}{dt} = \frac{1}{\lambda(t)}$

Adoption of self-similar decomposition: $a(t,Z) = \frac{1}{\lambda(s)}(\phi(z) + \tilde{a}(s,z))$$c(t,Z) = \frac{1}{\lambda(s)^{1+\sigma}}\tilde{c}(s,z)$

where $\phi(z) = e^{-z}$ is the approximate profile function.

3. Perturbation Equations

Through self-similar transformation, evolution equations for perturbation terms $(\tilde{a}, \tilde{c})$ and modulation parameters $(\lambda, \nu)$ are derived.

Technical Innovations

1. Profile Selection: Selection of smooth profile $\phi(z) = e^{-z}$ rather than conventional self-similar profiles, avoiding complex computations.
2. Differentiated Scaling Strategies:
   • Non-diffusive case: Temperature and velocity employ identical scaling
   • Diffusive case: Temperature employs heat equation scaling $\lambda^{-2}$
3. Vanishing Rate Conditions: Establishment of time-invariant vanishing rate conditions at the origin, ensuring well-posedness of weighted energy.
4. Bootstrap Argument: Combination strategy of interior energy methods and exterior maximum principle.

Experimental Setup

Mathematical Framework

This paper employs pure mathematical analysis methods without numerical experiments. Primary techniques include:

1. Function Spaces: Hölder spaces $C^{k,\alpha}$, Sobolev spaces, etc.
2. Weighted Hardy Inequalities: For handling singular weights
3. Maximum Principle: For exterior estimates
4. Gronwall Inequality: For solving differential inequalities

Parameter Settings

For the non-diffusive case, critical parameters satisfy:

• $\alpha_0 < \alpha < 3$, $1 < \gamma < 3$ (weight parameters)
• $1 < h_a < 2$, $0 < h_c < \gamma - 1$ (decay rates)

For the diffusive case:

• $\alpha_0 < \alpha < 5/2$
• $\eta_0 \geq 2$, $\eta_0 \in \mathbb{N}$ (regularity parameters)

Experimental Results

Main Results

Theorem 1 (Non-Diffusive Temperature Case)

There exists $\lambda_0^* > 0$ such that for initial data satisfying the conditions, the solution exhibits blowup at finite time $T$:

$a(t,Z) = (T-t)^{-1+Z} + O((T-t)^{-1}|\log(T-t)|^{-1/3})$

$c(t,Z) = \begin{cases} 0 & \text{if } Z = 0 \\ O((T-t)^{-3/4}) & \text{if } Z \in (0,1] \end{cases}$

Theorem 2 (Diffusive Temperature Case)

Similarly, the solution exhibits blowup:

$a(t,Z) = (T-t)^{-1+Z} + O((T-t)^{-1}|\log(T-t)|^{-1/2})$

$c(t,Z) = \begin{cases} 0 & \text{if } Z = 0 \\ O((T-t)^{-29/16}|\log(T-t)|^{3/4}) & \text{if } Z \in (0,1] \end{cases}$

Key Findings

1. Singularity Location: Singularities form at $Z = 0$ with blowup rate $(T-t)^{-1}$.
2. Temperature Effect: Temperature variations do not affect the singularity formation mechanism or stability in the velocity field.
3. Diffusion Effect: Diffusive terms alter the blowup rate of the temperature field but do not affect the dominant velocity singularity.
4. Stability: Under appropriate perturbations, singularity formation is stable.

Related Work

Historical Development

1. Lions-Temam-Wang (1990s): Established mathematical foundations of primitive equations
2. Cao-Titi (2005): Proved well-posedness of fully viscous primitive equations
3. Wong (2014): First proof of blowup in 2D hydrostatic Euler equations
4. Cao et al. (2015): Proved existence of singularities in 3D isothermal case
5. Collot et al. (2023): Established stable singularity theory for isothermal case

Contribution of This Paper

This paper is the first systematic investigation of singularity stability in the non-isothermal case, filling an important gap in the field.

Conclusions and Discussion

Main Conclusions

1. Universality: Temperature variations do not affect fundamental properties of velocity field singularities in primitive equations
2. Stability: Singularities remain stable under small perturbations
3. Precise Description: Exact asymptotic behavior of singularities is provided

Limitations

1. Special Geometry: Results are restricted to specific symmetry settings
2. Initial Value Conditions: Initial data must be close to specific profiles
3. Inviscid Assumption: Effects of viscosity are not considered

Future Directions

1. Investigation of singularity behavior in more general geometric settings
2. Consideration of viscous effects on singularity stability
3. Exploration of effects of other types of temperature boundary conditions
4. Numerical verification of theoretical predictions

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Mathematical proofs are complete and rigorous with refined technical treatment
2. Methodological Innovation: Different scaling strategies are employed for different temperature models
3. Complete Results: Exact blowup profiles and rates are provided
4. Physical Significance: Results have important implications for understanding singularities in geophysical fluids

Weaknesses

1. Limited Scope: Results are restricted to specific symmetry assumptions
2. Computational Complexity: Proof procedures are extremely technical and difficult to generalize
3. Lack of Numerical Verification: Absence of numerical experiments supporting theoretical results

Impact

1. Theoretical Contribution: Important supplement to singularity theory of primitive equations
2. Methodological Value: Developed techniques can be applied to other geophysical fluid models
3. Practical Significance: Provides theoretical foundation for understanding extreme events in atmosphere and ocean

Applicable Scenarios

This research is primarily applicable to:

1. Theoretical fluid mechanics research
2. Geophysical fluid dynamics
3. Singularity theory of partial differential equations
4. Mathematical modeling in atmospheric and oceanic science

References

The paper includes 22 important references covering classical and recent developments in primitive equation theory, singularity analysis, Hardy inequalities, and related fields. Key references include pioneering work by Lions et al., well-posedness theory by Cao-Titi, and stable singularity research by Collot et al. for the isothermal case.