光滑区域上欧拉斑块方程的全局正则性与小尺度快速生成

Global Regularity and Fast Small Scale Formation for Euler Patch Equation in a Smooth Domain

作者: Alexander Kiselev, Chao Li

类别: math.AP

发布日期: 2018-06-19

arXiv ID: 1806.07744v1

链接: http://arxiv.org/abs/1806.07744v1

详细分析

中文标题: 光滑区域上欧拉斑块方程的全局正则性与小尺度快速生成

研究对象和背景

本文的核心研究对象是二维有界光滑区域上的欧拉斑块方程，这是不可压缩理想流体力学的经典模型，描述涡度为分片常数的弱解的演化规律，其基础是二维不可压缩欧拉方程。考虑紧区域$D\subset\mathbb{R}^2$上的无穿透边值欧拉方程，形式如下：

$$ \begin{cases} \partial_t u + (u\cdot\nabla)u = -\nabla p, \\ \nabla\cdot u = 0, \\ u\cdot n|_{\partial D} = 0, \end{cases} $$

其中$u=(u_1,u_2)$是流体速度场，$p$是压强，$n$是区域$\partial D$的单位外法向量。为了简化分析，通常引入涡度$\omega=\text{curl}u=\partial_1 u_2-\partial_2 u_1$，将欧拉方程转化为涡度输运形式：

$$ \partial_t \omega + (u\cdot\nabla)\omega = 0, $$

速度场由Biot-Savart律重构：

$$ u = \nabla^\perp (-\Delta_D)^{-1}\omega, $$

这里$\nabla^\perp=(\partial_2,-\partial_1)$，$\Delta_D$是$D$上的Dirichlet拉普拉斯算子，$(-\Delta_D)^{-1}$表示其逆算子。该系统的解满足涡度沿流场输运的性质：若流映射$\Phi_t(x)$满足$\frac{d}{dt}\Phi_t(x)=u(\Phi_t(x),t),\Phi_0(x)=x$，则$\omega(x,t)=\omega_0(\Phi_t^{-1}(x))$，其中$\omega_0$是初始涡度。对于$\omega_0\in L^\infty\cap L^1(D)$的情形，Yudovich已证明存在唯一的全局弱解，称为Yudovich解，其流映射是全局Holder连续的。

欧拉斑块是一类特殊的Yudovich解，其涡度为分片常数，形式为$\omega(x,t)=\sum_{k=1}^N \theta_k\chi_{\Omega_k(t)}(x)$，其中$\theta_k$为非零常数，$\Omega_k(t)$是随时间演化的有界开集，闭包两两不交，边界$\partial\Omega_k(t)$为简单闭曲线。本文关注的核心问题是斑块边界的全局正则性，即初始足够光滑的边界是否会在演化中保持正则性，以及边界曲率的增长速度——曲率的快速增长对应流体高梯度小尺度结构的涌现，即小尺度快速生成现象，这是流体湍流起源研究的核心基础问题之一。

在本文发表之前，欧拉斑块的正则性理论已经有了一定进展：全空间$\mathbb{R}^2$下，Chemin（1993）、Bertozzi和Constantin（1993）分别独立证明，只要初始边界是$C^{1,\gamma}(\gamma\in(0,1))$的，斑块边界将永远保持该正则性，即全空间下欧拉斑块具有全局正则性。但有界区域的情形由于固壁边界的存在，Biot-Savart律的结构发生变化，会引入额外的边界奇性，相关理论非常不完备：Depauw（1999）证明半平面下初始不接触边界的单斑块具有全局$C^{1,\gamma}$正则性，但初始接触边界的斑块仅能保持有限时间$C^{1,\gamma}$正则；Dutrifoy（2003）证明初始接触边界的斑块存在全局解，但正则性会损失到$C^{1,s}(s<\gamma)$；Kiselev等人（2016）证明了半平面下多斑块的全局$C^{1,\gamma}$正则性，没有正则性损失，但一般光滑有界域下的欧拉斑块正则性仍是完全开放的问题，且所有现有结果都没有给出边界曲率增长的定量上界，也未验证上界的最优性，即有界域下的小尺度生成速度问题没有得到解决。

本文的研究设定为：求解区域是$C^4$光滑有界域，边值条件为无穿透固壁边界，斑块初始边界为$C^{1,\gamma}$正则，核心目标是建立一般光滑域下欧拉斑块的全局正则性理论，给出曲率增长的定量上界，并构造例子验证上界的最优性，刻画小尺度快速生成的速度。

主要定理或主要结果

本文共证明了四个核心定理，覆盖单斑块、多斑块的全局正则性，以及对称情形下的上界最优性：

  1. 定理1.3（单斑块全局正则性）：设$D$为$C^4$有界光滑区域，$\gamma\in(0,1)$，对任意$C^{1,\gamma}$单斑块初始数据$\omega_0$，存在唯一的全局$C^{1,\gamma}$斑块解$\omega$，其斑块边界的曲率增长速度最多为双指数。
  2. 定理1.4（多斑块全局正则性）：设$D$为$C^4$有界光滑区域，$\gamma\in(0,1)$，对任意$C^{1,\gamma}$多斑块初始数据$\omega_0$，存在唯一的全局$C^{1,\gamma}$斑块解$\omega$，其斑块边界的曲率增长速度最多为三指数。
  3. 定理1.5（单位圆盘对称双斑块的上界）：设$D$为单位圆盘，$\gamma\in(0,1)$，初始数据满足奇对称性$\omega_0(x)=\chi_{\Omega_1}(x)-\chi_{\Omega_2}(x)$，其中$\Omega_2$是$\Omega_1$关于$x_2$轴的反射，则对任意$C^{1,\gamma}$的该类初始数据，存在唯一的全局$C^{1,\gamma}$斑块解，其边界曲率增长速度最多为双指数。
  4. 定理1.6（对称情形的最优性例子）：在定理1.5的设定下，存在$C^{1,\gamma}$的初始数据，使得对应的斑块边界曲率确实以双指数速度增长，说明定理1.5的上界是定性最优的。

研究方法、关键技术和核心工具

核心分析学工具

本文的证明融合了调和分析、偏微分方程正则性估计、几何分析等多类工具，核心工具包括三类： 第一类是Dirichlet格林函数的分解与椭圆估计，这是处理有界域Biot-Savart律的核心基础。$D$上的Dirichlet格林函数$G_D(x,y)$满足$\Delta_x G_D(x,y)=\delta(x-y)$，$G_D(x,y)|_{x\in\partial D}=0$，首先可以分解为全空间奇性部分和调和修正项：

$$ G_D(x,y) = \frac{1}{2\pi}\log|x-y| + h(x,y), $$

其中$h(x,y)$是关于$x$的调和函数，边界值为$-\frac{1}{2\pi}\log|x-y|$，其高阶导数在区域内部有一致上界。为了处理边界附近的奇性，作者进一步引入边界管状邻域$T(r)=\{y\in\mathbb{R}^2:d(y,\partial D)\leq r\}$（$r$为仅依赖$D$的足够小常数），在该邻域内格林函数可以分解为镜像形式：

$$ G_D(x,y) = \frac{1}{2\pi}\left(\log|x-y| - \log|x-\tilde{y}|\right) + B(x,y), $$

其中$\tilde{y}=2e(y)-y$是$y$关于边界的镜像点，$e(y)$是$y$到边界的最近点，$B(x,y)$是正则项，满足对任意$\omega\in L^\infty(D)$，$\int_{T(r)\cap D}B(x,y)\omega(y)dy\in C^{2,\alpha}(T(r))$，其二阶导数有一致上界，该分解将边界附近的奇性转化为可处理的镜像项和无奇异的正则项，是本文最关键的技术工具。

第二类是Hölder空间的演化估计框架，借鉴了Bertozzi和Constantin的全空间斑块估计思路，将斑块边界的正则性转化为切向量场的范数估计。通过水平集表示斑块$\Omega(t)=\{x:\phi(x,t)>0\}$，定义切于边界的向量场$w=\nabla^\perp\phi$，再引入三个刻画边界正则性的核心量：

$$ \begin{aligned} A_\gamma(t) &= \|w(\cdot,t)\|_{\dot{C}^\gamma(\Omega(t))}, \\ A_\infty(t) &= \|w(\cdot,t)\|_{L^\infty(\Omega(t))}, \\ A_{\text{inf}}(t) &= \inf_{x\in\partial\Omega(t)}|w(x,t)|, \end{aligned} $$

若这三个量在任意有限时间内一致有界，则斑块边界保持$C^{1,\gamma}$正则。通过输运方程可以推导三个量的时间导数满足的演化不等式：

$$ \begin{aligned} A_\infty'(t) &\leq C(\gamma)A_\infty(t)\|\nabla u(\cdot,t)\|_{L^\infty}, \\ A_{\text{inf}}'(t) &\geq -C(\gamma)A_{\text{inf}}(t)\|\nabla u(\cdot,t)\|_{L^\infty}, \\ A_\gamma'(t) &\leq C(\gamma)\|\nabla u(\cdot,t)\|_{L^\infty}A_\gamma(t) + \|(\nabla u)w\|_{\dot{C}^\gamma(\Omega(t))}, \end{aligned} $$

只需估计出$\|\nabla u\|_{L^\infty}$和$\|(\nabla u)w\|_{\dot{C}^\gamma}$的上界，即可得到三个量的增长率，进而得到曲率的增长速度。

第三类是几何分析工具，包括边界管状邻域的构造、镜像映射的正则性分析、处理Hölder半范的路径积分技巧，以及刻画斑块边界附近几何性质的几何引理。例如为了避免直接估计穿过斑块边界的Hölder差，作者构造等边三角形的积分路径，保证路径上所有点到斑块边界的距离有一致下界；再例如Lemma2.13的几何引理，用来区分边界附近的点属于斑块内部还是外部，从而精准控制格林函数奇性的积分范围，避免发散。

技术性难题的解决技巧

本文的核心技术难点是处理固壁边界引入的额外奇性，以及多斑块之间的相互作用，作者采用的主要技巧包括： 第一，边界奇性的拆分与迁移技巧：通过格林函数的镜像分解，将速度场拆分为三部分：原斑块的全空间速度贡献、镜像斑块的速度贡献、正则项的贡献。正则项的导数有一致上界，无需额外处理；镜像斑块的贡献通过变量替换转化为带光滑权重的原斑块积分，权重的$C^2$正则性足够高，不会引入新的奇性，从而可以将全空间的Hölder估计技巧迁移到有界域情形，仅需额外处理光滑权重的影响，避免了边界奇性导致的估计发散。 第二，多斑块相互作用的尺度分离技巧：多斑块情形下，不同斑块之间的最小距离$\delta(t)$随时间最多双指数衰减，因此不同斑块对目标斑块的速度梯度贡献，可以通过距离的负幂次估计得到$\|\nabla u_i\|_{L^\infty}\leq C\delta(t)^{-3}$，将其代入演化不等式后，即可得到三指数的增长上界。 第三，对称情形的奇性抵消技巧：在单位圆盘的对称双斑块情形下，利用初始数据的奇对称性，镜像斑块的贡献恰好等于另一个斑块贡献的相反数，从而抵消了多斑块相互作用的高阶衰减项，将上界从三指数降回双指数。 第四，小尺度生成的构造技巧：借鉴Kiselev和Sverák的全平面小尺度生成构造思路，利用对称情形下边界附近的点会被速度场以越来越快的速度推向对称轴的性质，构造初始斑块使得边界上的点以双指数速度趋近对称轴，导致边界的曲率双指数增长，验证了上界的最优性。

与之前工作的比较

本文是有界域欧拉斑块正则性领域的突破性工作，相比之前的研究，其改进和创新主要体现在以下几个方面：

假设条件的放宽

首先是区域假设的大幅放宽：之前的有界域欧拉斑块正则性结果仅针对半平面这类特殊的无界区域，而本文首次将结果推广到任意$C^4$光滑有界域，覆盖了所有实际应用中常见的光滑有界流动区域，解决了一般光滑域下欧拉斑块全局正则性的长期开放问题。半平面的镜像映射是线性的，估计难度较低，而一般光滑域的镜像映射是非线性的，需要额外处理镜像映射的正则性和权重的影响，本文发展的格林函数镜像分解技巧，为此类一般有界域的有源标方程正则性研究提供了通用框架。 其次是初始条件的放宽：之前半平面下初始接触边界的斑块，仅能得到有限时间$C^{1,\gamma}$正则，或者全局更低的$C^{1,s}(s<\gamma)$正则，存在明显的正则性损失。本文证明无论初始斑块是否接触边界，都能永久保持初始的$C^{1,\gamma}$正则性，没有任何正则性损失，直接改进了Depauw和Dutrifoy的半平面结果，说明固壁边界的存在不会导致斑块边界的正则性损失，仅会加快曲率的增长速度。

结果强度的提升

首先是定量结果的突破：之前的所有正则性结果仅定性证明了解的存在性和正则性，完全没有给出边界曲率的增长速度，本文首次给出了单斑块双指数、多斑块三指数的定量上界，并且在对称情形下构造了双指数增长的例子，证明了上界的定性最优性，首次在有界域下严格证明了欧拉斑块的小尺度快速生成现象，填补了该领域的空白。 其次是小尺度生成结果的创新：之前Kiselev和Sverák的小尺度生成结果是全平面光滑初值的涡度梯度双指数增长，而本文是有界域下斑块弱解的边界曲率双指数增长，是第一个有界域下欧拉方程小尺度快速生成的严格结果，而且针对的是更贴近实际流动的斑块模型，对湍流的小尺度生成理论有重要的支撑意义。 第三是多斑块结果的首创性：无论是全空间还是有界域，之前的研究都没有给出多斑块情形的曲率增长上界，本文首次给出了多斑块的三指数上界，尽管该上界的最优性尚未验证，但为后续多斑块相互作用的研究提供了第一个定量基准。

技术方法的通用性

本文提出的边界附近格林函数镜像分解+正则项估计的方法，具有很强的通用性，可以推广到其他有源标方程的有界域问题，例如表面准地转方程的斑块正则性问题、耗散欧拉方程的边界问题、有源标方程的边界层问题等，为后续有界域流体方程的正则性研究提供了新的技术路径。 当然本文也存在一定的局限性：多斑块的三指数上界是否最优尚未得到验证，目前没有构造出对应的下界例子；区域的正则性假设是$C^4$，未来可以尝试降低到更低的正则性，例如$C^{2,\alpha}$；本文的结果仅针对二维情形，三维欧拉斑块的正则性仍是完全开放的问题，本文的方法是否可以推广到三维也值得进一步研究。