[论文解读] Finite-time Singularity Formation for Strong Solutions to the Boussinesq System

Tarek M. Elgindi, In‐Jee Jeong|arXiv (Cornell University)|Aug 9, 2017

Navier-Stokes equation solutions参考文献 41被引用 26

一句话总结

本文在扇形区域 $\Omega_\gamma$ 上构造了二维Boussinesq方程的有限能量、$C^\infty$-光滑解，这些解尽管二维欧拉方程始终保持全局正则性，却在有限时间内形成奇点。奇点源于由密度梯度引起的涡量放大，该现象通过尺度不变解与截断论证实现，以确保有限能量，且奇点在 $t \to 1^-$ 时出现在速度和密度的梯度中。

ABSTRACT

The global regularity problem for the Boussinesq system is a well known open problem in mathematical fluid dynamics. As a follow up to our work \cite{EJSI}, we give examples of finite-energy and Lipschitz continuous velocity field and density $(u_0,ρ_0)$ which are $C^\infty$-smooth away from the origin and belong to a natural local well-posedness class for the Boussinesq equation whose corresponding local solution becomes singular in finite time. That is, while the sup norm of the gradient of the velocity field and the density remain finite on the time interval $t\in [0,1)$, both quantities become infinite as $t ightarrow 1$. The key is to use scale-invariant solutions similar to those introduced in \cite{EJSI}. The proof consists of three parts: local well-posedness for the Boussinesq equation in critical spaces, the analysis of certain special infinite-energy solutions belonging to those critical spaces, and finally a cut-off argument to ensure finiteness of energy. All of this is done on spatial domains $\{(x_1,x_2): x_1 \ge γ|x_2|\}$ for any $γ> 0$ so that we can get arbitrarily close to the half-space case. We remark that the $2D$ Euler equation is globally well-posed in all of the situations we look at, so that the singularity is not coming from the domain or the lack of smoothness on the data but from the vorticity amplification due to the presence of a density gradient. It is conceivable that our methods can be adapted to produce finite-energy $C^\infty$ solutions on $\mathbb{R}^2_+$ which become singular in finite time.

研究动机与目标

• 证明在具有有限能量的关键正则性类中，二维Boussinesq系统强解的有限时间奇点形成。
• 分离由密度梯度引起的涡量放大导致奇点形成的机制，该机制与二维欧拉方程不同。
• 在具有锐角（$\gamma > 0$）的区域上构造 $C^\infty$-光滑、有限能量的初始数据，使其在有限时间内导致梯度爆破。
• 将尺度不变解方法扩展至临界空间，并使用截断论证确保有限能量，同时保持奇点形成特性。
• 表明奇点并非由边界不规则性或初始数据光滑性不足引起，而是由涡量与密度梯度之间的非线性耦合所致。

提出的方法

• 在扇形区域 $\Omega_\gamma = \{(x_1,x_2) : x_1 \geq \gamma |x_2|\}$（$\gamma > 0$）上构造二维Boussinesq系统的尺度不变解，其为一次齐次函数。
• 利用在原点处从完整Boussinesq方程约化得到的1D系统，捕捉涡量 $\omega(0)$ 和密度梯度 $\partial_{x_1}\rho(0), \partial_{x_2}\rho(0)$ 的演化。
• 推导关于 $A = \omega(0)$, $B = \partial_{x_1}\rho(0)$, $C = \partial_{x_2}\rho(0)$ 的常微分方程组：$A' = C$, $B' = -\frac{1}{1-\beta^2}AC$, $C' = -\frac{\beta^2}{1-\beta^2}AB$，其中 $\beta = \gamma^{-1}$。
• 证明当 $\beta < 1/2$（即 $\gamma > 2$）时，若 $C_0 > 0$，则系统在有限时间内发生爆破，这是由于 $C$ 和 $A$ 呈指数增长。
• 应用截断过程以局部化尺度不变解，确保有限能量的同时保留奇点形成机制。
• 在临界 Hölder 空间 $C^{k,\alpha}$（$\alpha < 1/\beta - k - 2$）中建立局部适定性，确保解在爆破前保持在所需正则性类中。

实验结果

研究问题

• RQ1 在二维欧拉方程全局正则的背景下，二维Boussinesq系统是否存在有限能量、$C^\infty$-光滑解在有限时间内形成奇点？
• RQ2 在无粘性条件下，密度梯度在促进涡量放大与奇点形成中起什么作用？
• RQ3 尺度不变解能否用于构造有限能量解，使其在有限时间内爆破，即使原始解为无限能量？
• RQ4 奇点形成机制在区域截断与截断操作下是否稳健？能否将其推广至上半平面 $\mathbb{R}^2_+$？
• RQ5 区域的几何形状（特别是扇形的夹角）如何影响Boussinesq系统中有限时间爆破的可能性？

主要发现

• 在扇形区域 $\Omega_\gamma$（$\gamma > 0$）上存在有限能量、$C^\infty$-光滑的初始数据，使得对应的二维Boussinesq系统强解在有限时间内形成奇点。
• 当 $t \to 1^-$ 时，速度场的梯度和密度均发生爆破，而初始数据始终保持在 $W^{1,\infty}$ 中且为Lipschitz连续。
• 奇点形成源于由密度梯度驱动的涡量非线性放大，而非边界效应或光滑性不足所致。
• 当 $\beta < 1/2$（即 $\gamma > 2$）时，若 $\partial_{x_2}\rho_0(0) > 0$，则原点处约化的ODE系统在有限时间内发生爆破，且 $A(t)$ 与 $C(t)$ 无界增长。
• 在相同条件下，二维欧拉方程仍保持全局适定，确认奇点是Boussinesq系统特有现象，源于密度-涡量耦合。
• 该方法具有足够鲁棒性，可推测在 $\mathbb{R}^2_+$ 上的有限能量 $C^\infty$ 解也可能在有限时间内形成奇点，但需进一步分析。