[论文解读] Near identity transformations for the Navier-Stokes equations

Peter Constantin|ArXiv.org|Dec 12, 2001

Lattice Boltzmann Simulation Studies参考文献 46被引用 24

一句话总结

本文提出了一套近恒等扩散变换的框架，用于分析纳维-斯托克斯方程，其中速度和涡量通过虚拟场在粘性修正传输下演化，利用韦伯公式和柯西公式进行重构。关键贡献在于建立了关于解正则性的时变界，表明光滑行为在由涡度能和粘性决定的函数下限所界定的时间区间内持续存在，为全局正则性问题提供了新方法。

ABSTRACT

The Navier-Stokes equations and their various approximations can be described in terms of near identity maps, that are diffusive particle path transformations of physical space. The active velocity is obtained from the diffusive path transformation and a virtual velocity using the Weber formula. The active vorticity is obtained from the diffusive path transformation and a virtual vorticity using a Cauchy formula. The virtual velocity and the virtual vorticity obey diffusive equations, which reduce to passive advection formally, if the viscosity is zero. Apart from being proportional to the viscosity, the coefficients of these diffusion equations involve second derivatives of the near identity transformation and are related to the Christoffel coefficients. If and when the near-identity transformation departs excessively from the identity, one resets the calculation. Lower bounds on the minimum time between two successive resettings are given in terms of the maximum enstrophy.

研究动机与目标

• 开发一种概念与分析框架，通过物理空间的扩散型近恒等变换来理解纳维-斯托克斯方程。
• 通过近似方法分析长时间行为，解决三维纳维-斯托克斯流中全局正则性的开放问题。
• 基于涡度能和粘性，建立解保持光滑的时间区间的定量下界。
• 在基于变换的统一形式下，整合伽辽金法、滤波型方法、涡方法及滤波公式等多种数值近似。

提出的方法

• 该方法采用从恒等映射出发并随时间演化的扩散型粒子路径变换，只要雅可比行列式保持接近1，即可维持可逆性。
• 通过变换和一个虚拟速度，利用韦伯公式重构主动速度；通过变换和一个虚拟涡量，利用柯西公式获得主动涡量。
• 虚拟速度和涡量场服从带粘性系数与变换二阶导数成比例的扩散方程，在无粘极限下退化为被动输运。
• 当变换显著偏离恒等（即雅可比行列式变化超过阈值）时，使用当前主动场作为新的初始虚拟场重置计算。
• 该框架适用于滤波型方程和涡方法，后者通过变量变换与纳维-斯托克斯方程的滤波公式相关联。
• 利用索博列夫型范数和插值不等式，推导出速度和梯度范数的先验界，其依赖于涡度能和粘性。

实验结果

研究问题

• RQ1 纳维-斯托克斯方程及其近似能否通过物理空间的近恒等扩散变换系统地描述？
• RQ2 在涡度能和粘性参数下，纳维-斯托克斯方程光滑解可保证存在的最小时间区间是多少？
• RQ3 在该变换框架下，那些精确保持能量耗散的近似与那些精确保持涡量方程的近似之间有何关系？
• RQ4 通过涡度能和变换雅可比行列式的结构，解的正则性在多大程度上可被控制？
• RQ5 该框架能否在解并非全局光滑的情况下，仍对多个时间区间内的解梯度提供一致的界？

主要发现

• 建立了光滑解保持时间区间的下界，其与 $ \nu \mathcal{E}^{-2} $ 成正比，其中 $ \mathcal{E} $ 为涡度能上界，$ \nu $ 为粘性系数。
• 对任意 $ t_1 \geq t_0 $，方程 (49) 的零初值解 $ \ell $ 满足 $ \|\ell(\cdot,t)\|_{\{A,r,1\}} \leq (t-t_1) U_r e^{c (t-t_1) U_r^2 / \nu} $，其中 $ U_r $ 依赖于涡度能和尺度。
• 梯度满足 $ \|\nabla\ell(\cdot,t)\|_{\{A;r_1,1\}} \leq g $，对 $ t \in [t_1, t_1 + \tau T] $ 成立，其中 $ \tau \sim g G^{-7} $，$ G $ 为无量纲涡度能参数。
• 对 $ \tilde{r} = (1-\gamma)\tilde{\lambda} $，速度范数满足 $ \tilde{U}_{\tilde{r}} \leq c_3 G^3 $，表明在无量纲单位下对涡度能有强依赖性。
• 二阶梯度满足 $ \|\nabla\nabla\ell(\cdot,t)\|_{\{A;r_2,1\}} \leq c_5 (\nu T)^{-1/2} G^4 g $，表明高阶正则性得到控制。
• 结果仅通过涡度能上界 $ \mathcal{E} $ 依赖于滤波截断，意味着对涡度能有界的所有近似均具有统一控制。