生成时间: 2025年12月30日 03:53:55

重点关注论文数量: 7

了解领域论文数量: 6

不相关论文数量: 25

重点关注论文（完整分析）

1. 带速度阻尼项的理想不可压缩磁流体动力学方程组的全局稳定性与渐近行为

Global stability and asymptotic behavior for incompressible ideal MHD equations with velocity damping term

作者: Hui Fang, Pingping Gui, Yanping Zhou

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23263v1

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

详细分析

中文标题: 带速度阻尼项的理想不可压缩磁流体动力学方程组的全局稳定性与渐近行为

研究对象和背景

本文研究的对象是带有速度阻尼项的多维理想不可压缩磁流体动力学（MHD）方程组。MHD方程组描述了导电流体（如液态金属、等离子体）在磁场中的运动，在地球物理和天体物理流体中具有基础性作用。当电阻效应极其微弱，磁扩散可以忽略不计时，系统退化为非电阻性MHD系统。本文重点关注以下方程组（即文中的方程(1.1)），其中参数 $\alpha \ge 0$ 代表速度耗散的强度，$\alpha=0$ 对应于纯阻尼效应：

其中 $v(x,t)$ 表示流体速度，$b(x,t)$ 表示磁场，$p(x,t)$ 表示压力。作者研究的特定模型对应于 $\alpha = 0$ 的情况，即带有纯速度阻尼项 $v$ 的理想MHD系统（文中的方程(1.4)）。该模型在物理上对应于等离子体受到强烈的碰撞阻尼或外在的摩擦松弛机制。

研究区域与边值条件：该研究在 $n$ 维环面 $\mathbb{T}^n$ ($n \ge 2$) 上进行，即具有周期性边界条件。此外，研究假设流体和磁场满足无散度条件 ($\nabla \cdot v = \nabla \cdot b = 0$) 且在全空间上的积分为零（零平均条件），即 $\int_{\mathbb{T}^n} v dx = 0, \int_{\mathbb{T}^n} b dx = 0$。

现象验证与背景： 该论文旨在探究一个非平凡的稳态磁场背景下的稳定性问题。背景磁场记为 $\bar{b}$，它满足 $v=0, b=\bar{b}$ 是系统的稳态解。物理直觉和数值模拟都表明磁场具有“稳定化”效应，即能够抑制流体的湍流和扰动的增长。本文数学上刻画了背景磁场的这种稳定化效应，并研究了小扰动在长时间下的渐近行为。

为了实现这一目标，论文假设背景磁场 $\bar{b}$ 满足 丢番图条件 (Diophantine condition) (方程(1.3))：

其中 $r > n-1, c > 0$。这一条件保证了背景磁场在傅里叶空间中具有非共振特性，从而避免了特征值的共振重合，使得线性化算子具有良好的谱性质和色散效应（类似于相位混合）。

主要定理或主要结果

本文主要包含两个核心定理，分别描述了线性化系统的衰减行为和非线性系统的全局稳定性及渐近行为。

定理 1.1 (线性化系统的衰减) 设 $n \ge 2$，$r > n-1$，背景磁场 $\bar{b}$ 满足丢番图条件。假设初始数据 $(V_0, H_0) \in H^m(\mathbb{T}^n)$，整数 $m \ge 1$，且满足零平均条件。则线性化系统 (1.5) 的解 $(V, H)$ 满足如下衰减估计：

定理 1.2 (非线性系统的全局稳定性与衰减) 设 $m \in \mathbb{N}$ 满足以下正则性要求：

假设 $\bar{b}$ 满足丢番图条件，初始数据 $(v_0, b_0) \in H^m(\mathbb{T}^n)$ 满足零平均条件且足够小。则非线性 MHD 系统 (1.4) 存在唯一的全局经典解 $(v,b)$，且满足能量估计：

并且解具有衰减估计：

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

本文的证明方法主要依赖于傅里叶分析和精细的能量估计。作者构建了一个通用的分析框架，适用于许多具有部分耗散性的流体模型。以下是证明中的关键技术细节：

a) 核心分析学工具

1. 谱分析与矩阵指数 (Spectral Analysis and Matrix Exponential)： 为了处理线性部分，作者将线性化系统写成矩阵形式：

   $$ \partial_t \hat{\psi} + Q \hat{\psi} + \hat{N} = 0, \quad \text{其中 } Q = \begin{pmatrix} 1 & -i\bar{b}\cdot j \\ -i\bar{b}\cdot j & 0 \end{pmatrix}, \ \hat{\psi} = (\hat{v}, \hat{b})^T. $$

   通过对矩阵 $Q$ 进行对角化（利用其特征值 $\lambda_1, \lambda_2$ 和特征向量），可以显式写出解的积分表示（Duhamel 公式）：

   $$ \hat{\psi} = e^{-Qt} \hat{\psi}_0 - \int_0^t e^{-Q(t-\tau)} \hat{N}(\tau) d\tau. $$

   这一步将微分方程转化为积分方程，是进行逐点估计的基础。

2. Leray 投影 (Helmholtz-Leray Projection)： 在处理非线性项之前，必须剔除压力项。通过 Leray 投影算子 $P$ 将速度方程投影到无散度向量场空间中，从而消去压力 $\nabla p$，确保算子的结构适合进行能量估计。

3. Sobolev 空间中的能量估计： 核心工具是在 $H^s$ 空间中构建能量泛函。由于存在阻尼项 $v$，这与理想流体或无阻尼 MHD 有本质不同。作者通过对时间微分 $\frac{1}{2}\frac{d}{dt}\|U\|_{H^s}^2$ 进行估计，并利用阻尼项提供的正定项 $\|v\|_{H^s}^2$ 来控制非线性项和耦合项。

b) 解决技术性难题的技巧

1. 频率空间的分裂与核函数估计 (Splitting of Frequency Space)： 这是处理丢番图条件和特征值奇点的关键技巧。

   • 特征值的正则性：矩阵 $Q$ 的特征值 $\lambda_1, \lambda_2$ 取决于 $\sqrt{1-4|\bar{b}\cdot j|^2}$。当 $|\bar{b}\cdot j| \approx 1/2$ 时，根号内接近 0，产生奇点。为了解决这个问题，作者引入了辅助核函数 $G_2$（如 Proposition 3.1 所示），通过分析 $(e^{-\lambda_2 t} - e^{-\lambda_1 t})/(\lambda_2 - \lambda_1)$ 的有界性来平滑奇点。
   • 低频与高频分离：为了得到最优的衰减率，作者将频率变量 $j$ 分裂为不同的区域（如 $|j| \le J_0$ 和 $|j| > J_0$，或者根据 $|\bar{b}\cdot j|$ 的大小）。在低频区利用 $L^\infty$ 估计和 $L^2$ 能量转移；在高频区利用对角化后的显式衰减率进行逐点估计。

2. 丢番图条件的应用 (Application of Diophantine Condition)： 丢番图条件保证了对于非零整数向量 $j$，$\bar{b}\cdot j$ 不会过小。这一性质被用于两个方面：

   • 线性 Poincaré 不等式 (Lemma 2.3)：$\|g\|_{H^s} \le c \|\bar{b}\cdot \nabla g\|_{H^{s+r+1}}$。这是连接速度场和磁场的核心不等式，它表明由于背景场的存在，磁场的导数可以控制速度场的高阶范数，反之亦然。这体现了背景磁场提供的“恢复力”。
   • 衰减率的优化：丢番图条件中的指数 $r$ 直接出现在最终的衰减指数 $-(m-s)/(2(1+r))$ 中。条件越强（$r$ 越小），衰减越快。

3. Besov 空间与正则性指标的精细处理： 论文在处理非线性项 $\Lambda^s (v\cdot \nabla v)$ 等项时，使用了交换子估计 (Lemma 2.4)。在证明非线性稳定性（Section 5）时，作者需要处理不同 Sobolev 指标之间的插值和转移。特别是，定理 1.2 中给出了能量估计 $\int \|b\|_{H^{m-r-1}}^2 dt$，这比 $\int \|b\|_{H^m}^2 dt$ 弱。这意味着在构造整体解的能量估计时，必须容忍磁场在低正则性空间中的耗散，或者通过精细的迭代（Picard 迭代）将低阶能量逐步提升到高阶能量。

4. 加权能量估计 (Weighted Energy Estimates)： 为了证明长时间的衰减行为，论文结合了短时间的能量有界性和长时间的代数衰减。通过将解的频率分量乘以 $(1+t)^\gamma$ 并对其微分，可以证明加权能量的递减性，从而导出 $(1+t)^{-\alpha}$ 形式的衰减。

与之前工作的比较

本文在多个维度上显著改进了前人的工作，特别是在降低正则性要求、推广维度以及细化衰减行为方面。

a) 假设条件的放宽与维度的推广

1. 从 $n=3$ 推广到 $n \ge 2$： 之前的代表性工作，如 Zhao 和 Zhai (2021) [45]，主要关注三维周期域 $\mathbb{T}^3$ 上的模型。本文将结果推广到了 $n$ 维环面 $\mathbb{T}^n$ ($n \ge 2$)，使得理论框架更加普适。

2. 初始数据正则性要求的显著降低：

   • Zhao & Zhai (2023)：要求初始数据属于 $H^m(\mathbb{T}^3)$，且 $m \ge 4r+7$。这意味着如果 $r$ 较大（即丢番图条件较弱），则需要非常高阶的导数存在。
   • 本文 (Theorem 1.2)：要求 $m > 1 + 2r + n/2$。

   对比分析：对于三维情况 ($n=3$)，本文要求 $m > 1 + 2r + 1.5 = 2.5 + 2r$。相比于 Zhao & Zhai 的 $m \ge 4r+7$，本文所需的正则性指标 $m$ 显著降低。例如，当 $r=2$ 时，前人需要 $m \ge 15$，而本文仅需 $m > 6.5$（即 $m \ge 7$）。这种降低归功于更精细的核函数分析和丢番图条件的深入利用，减少了对高阶 Sobolev 范数的依赖。

b) 结果的强化：衰减行为的完整性与速度

1. 衰减范围的扩展 (Completeness of Decay Range)：

   • 前人结果：Zhao & Zhai 证明的衰减估计 (1.13) $\|(v,b)\|_{H^s} \le C(1+t)^{-3(m-s)/2(m-r-4)}$ 仅在 $s \in [r+4, m)$ 范围内成立。这意味着对于较低阶的 Sobolev 范数（如 $H^s$ 当 $s < r+4$），前人没有给出衰减界。
   • 本文结果：定理 1.2 给出的衰减估计 (1.12) $\|(v,b)\|_{H^s} \le C(1+t)^{-\frac{m-s}{2(1+r)}}$ 对 所有 $s \in [0, m]$ 都成立。这提供了一个完整的从低正则性（弱范数）到高正则性（强范数）的衰减剖面，表明解在任意 Sobolev 指标下均以代数速率衰减。

2. 衰减速度的“饱和”与“发散”现象：

   • 前人结果：当 $m \to +\infty$ 时，Zhao & Zhai 的衰减指数 $\frac{3(m-s)}{2(m-r-4)}$ 趋向于 $\frac{3}{2}$（针对 $s$ 固定的情况，或者更准确地说，针对 $H^s$ 的衰减率有上限）。这意味随着初始数据正则性的提高，长时间衰减速度会趋于饱和，不会无限加快。
   • 本文结果：当 $m \to +\infty$ 时，本文的衰减指数 $\frac{m-s}{2(1+r)}$ 趋向于 $+\infty$。

   物理意义与数学优势：这揭示了一个深刻的数学结构——在带有丢番图背景磁场和速度阻尼的系统中，初始数据的高阶正则性直接转化为更快的时间衰减率。而在经典抛物型系统（如热方程）中，虽然高阶范数衰减也快，但通常受限于传播子的性质。本文证明了该系统比前人发现的系统具有“本质上更强的代数衰减”。

3. 与粘性耗散 ($\alpha=1$) 的对比： 论文在 Remark 1.1 中对比了速度阻尼 ($\alpha=0$) 与粘性扩散 ($\alpha=1$) 的线性衰减率。结果表明，速度阻尼产生的衰减比粘性扩散更快（对于所有 $s \in [0, m-1]$）。这是因为直接阻尼均匀作用于所有频率模式，而粘性扩散在高频端由于与磁场的耦合而效率降低。这一观察在数学上解释了为何阻尼项在稳定流体方面可能比粘性项更有效。

综上所述，本文不仅将前人的结果推广到了更一般的维度，大幅降低了所需的初始正则性，更重要的是，它揭示了该类系统衰减行为的完整图景以及正则性与衰减速度之间的深刻联系，这些是之前文献未曾发现的。

2. 无温变可压缩向列液晶流在尺度不变小性条件下的全局强解

Global strong solutions for non-isothermal compressible nematic liquid crystal flows under a scaling-invariant smallness condition

作者: Lin Xu, Xin Zhong

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23197v1

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

详细分析

中文标题: 无温变可压缩向列液晶流在尺度不变小性条件下的全局强解

1. 研究对象和背景

本文研究了三维带远场真空的**非等温可压缩向列液晶（Non-isothermal compressible nematic liquid crystal）**流动的柯西问题。液晶材料介于普通液体和固态晶体之间，其分子具有取向性。在宏观尺度上，Ericksen-Leslie 连续介质理论描述了这类流动，该理论将流体速度与分子取向（即导向场）耦合在一起。

论文关注的数学模型是一个耦合的非线性偏微分方程组，包含质量守恒方程、动量守恒方程、能量守恒方程以及导向场演化方程。在设定常数 $R, c_v, \kappa$ 均为 1 之后，该系统（即文中的方程组 (1.1)）可以写为：

这里未知函数分别为：

• $\rho \ge 0$：流体密度；
• $u \in \mathbb{R}^3$：速度场；
• $\theta \ge 0$：绝对温度；
• $d \in S^2$：表示液晶分子平均取向的单位向量场（导向场）。

状态方程为 $P = \rho \theta$（假设 $R=1$）。$\mu > 0$ 和 $\lambda$ 分别是剪切粘性系数和体粘性系数，满足物理条件 $\mu > 0, 2\mu + 3\lambda \ge 0$。$\kappa=1$ 是热传导系数，$c_v=1$ 是热容。

研究背景涉及液体晶体动力学中非等温效应的建模与分析。温度 $\theta$ 不仅通过状态方程影响压力，还通过能量方程中的粘性耗散和热传导与速度场及导向场发生复杂的耦合。导向场方程本质上是一个具有输运项的热流形式的谐波映射方程（transported heat flow of harmonic maps）。

论文研究的是柯西问题（Cauchy problem），即考虑全空间 $\mathbb{R}^3$ 上的演化。边界条件体现为远场行为：

这意味着流体在无穷远处趋于静止且处于低温状态，同时导向场趋于一个固定的单位向量 $\mathbf{1}$。

此外，问题的一个重要特征是允许初始真空（initial vacuum），即 $\rho_0$ 可能在某些区域为零。这使得动量方程和能量方程的系数退化，给正则性估计带来很大困难。作者使用的函数空间主要是标准的 Sobolev 空间，如 $L^q, H^k, W^{k,q}$ 以及 $D^{k,q}$ 等，并考虑了强解（Strong solution）的概念。

2. 主要定理或主要结果

本文的核心定理是 Theorem 1.1，它确立了全局强解的存在唯一性。

假设条件：

1. 初始正则性：初始数据 $(\rho_0, u_0, \theta_0, d_0)$ 满足：
   • $\rho_0 \in H^1 \cap W^{1,q}$ (其中 $q \in (3,6]$), $\rho_0 \ge 0$;
   • $\sqrt{\rho_0} u_0, \sqrt{\rho_0} \theta_0 \in L^2$;
   • $u_0, \theta_0 \in D^1 \cap D^2$;
   • $\nabla d_0 \in H^2$, $|d_0|=1$.
2. 相容性条件：存在 $g_1, g_2 \in L^2(\mathbb{R}^3)$ 使得：
   $$ \begin{cases} -\mu \Delta u_0 - (\lambda + \mu) \nabla \text{div} u_0 + \nabla (\rho_0 \theta_0) + \nabla d_0 \cdot \Delta d_0 = \rho_0 g_1, \\ -\Delta \theta_0 - Q(\nabla u_0) - |\Delta d_0 + |\nabla d_0|^2 d_0|^2 = \sqrt{\rho_0} g_2. \end{cases} $$
3. 尺度不变的小性条件：定义了如下的尺度不变量 $S_0$：
   $$ S_0 \triangleq \left( 1 + \bar{\rho} + \frac{1}{\bar{\rho}} \right) S' S'' $$
   其中 $\bar{\rho} = \|\rho_0\|_{L^\infty}$，
   $$ S' = \|\rho_0\|_{L^3} + (\bar{\rho}^2 + \bar{\rho}) (\|\sqrt{\rho_0} u_0\|_{L^2}^2 + \|\nabla d_0\|_{L^2}^2), $$
   $$ S'' = \|\nabla u_0\|_{L^2}^2 + (\bar{\rho} + 1) \|\sqrt{\rho_0} \theta_0\|_{L^2}^2 + \|\nabla^2 d_0\|_{L^2}^2 + \|\nabla d_0\|_{L^4}^4. $$
   要求 $S_0 \le \varepsilon_0$，其中 $\varepsilon_0$ 是一个仅依赖于物理参数（$R, \mu, \lambda, \kappa, c_v$）而不依赖于初始数据的小常数。

结论： 在上述假设下，柯西问题 (1.1)-(1.3) 存在唯一的全局强解 $(\rho, u, \theta, d)$。

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

本文的证明采用了典型的**先验估计（a priori estimates）结合连续性论证（bootstrap argument）**的方法。由于直接证明全局存在性极其困难，作者首先假设解在有限时间内存在（由局部存在性定理保证），然后推导出解的各种能量估计，证明只要初始数据的某个尺度不变量足够小，这些估计就能被初始值控制住，从而排除解在有限时间内爆破（blow-up）的可能性。

以下是证明中使用的核心分析工具和技术细节：

a) 核心分析学工具

1. 能量估计（Energy Estimates）： 这是最基础的工具。通过对动量方程点乘 $u$、温度方程点乘 $\theta$、以及利用导向场的几何约束 $|d|=1$ 进行估计，获得如 $\|\sqrt{\rho} u\|_{L^2}$, $\|\nabla \theta\|_{L^2}$, $\|\nabla d\|_{L^2}$ 等的 $L^2$ 范数控制。
2. 椭圆估计（Elliptic Estimates）： 由于动量方程和能量方程中含有高阶项，作者经常将时间导数项转化为椭圆系统求解。例如，文中引入 $F = (2\mu + \lambda)\text{div} u - P$，并推导了方程组 (2.20)：
   $$ \begin{cases} \Delta F = \text{div}(\rho \dot{u} + \nabla d \cdot \Delta d), \\ \mu \Delta \text{curl} u = \text{curl}(\rho \dot{u} + \nabla d \cdot \Delta d). \end{cases} $$
   利用 $L^p$ 理论（或极大正则性）将 $F$ 和 $\text{curl} u$ 的 $L^3$ 范数用 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 和 $\|\nabla d\|_{L^3} \|\Delta d\|_{L^2}$ 等来控制。这是在不依赖 $2\mu > \lambda$ 条件下控制压力和速度的关键。
3. Sobolev 嵌入与插值不等式： 用于在不同 $L^p$ 空间和 Sobolev 空间之间转换范数，特别是处理非线性项时必不可少。例如，处理 $|\nabla d|^4$ 项时利用 $\|\nabla d\|_{L^4}^4 \le C \|\nabla d\|_{L^2} \|\nabla^2 d\|_{L^2}^3$ 这种插值。

b) 解决技术性难题的主要技巧

1. 处理密度真空退化与粘性系数限制（核心创新点）：

   • 难题： 在动量方程中，项 $\rho \dot{u}$ 是退化的。标准方法通常需要假设 $2\mu > \lambda$ 来获得关于 $u$ 的高阶能量估计（如 $\|\nabla u\|_{L^2}$ 的耗散估计），从而利用 Gronwall 不等式。然而，物理上 $\lambda$ 可以非常大，作者希望移除这一限制。
   • 技巧： 作者没有直接利用含时导数的高阶估计，而是采用了类似于处理不可压缩流或特定可压缩流的策略。
     • 压力解耦与 $L^3$ 估计： 作者引入了 $\|\rho\|_{L^3}$ 的估计（Lemma 2.3）。通过将方程作用于 $(-\Delta)^{-1}\text{div}$ 并结合 $\partial_t \rho^3 + \text{div}(\rho^3 u) = ...$，得到了 $\|\rho(t)\|_{L^3}$ 的估计，其增长受控于 $(\bar{\rho}^2 + \bar{\rho}) \int_0^T (\|\nabla u\|_{L^2}^2 + \|\nabla^2 d\|_{L^2}^2) dt$。这是关闭 Bootstrap 循环的关键一环。
     • 温度方程的直接处理： 与之前工作不同，作者没有使用总能量方程（通常更简单，但可能丢失细节），而是直接处理温度方程 (2.18)，并利用 $L^3$ 估计中的 $\|\rho\|_{L^3}$ 来控制与速度相关的项。
     • 精细的耦合项估计： 详细分析了导向场项 $\nabla d \cdot \Delta d$ 对动量的影响。引入了辅助量 $K_1(t)$ 和 $K_2(t)$ 来捕捉这些项与温度梯度、速度梯度的相互作用。

2. 尺度不变性（Scaling Invariance）的利用：

   • 定义： 系统在变换 $\rho_k(x,t) = \rho(kx, k^2 t), u_k(x,t) = k u(kx, k^2 t), \theta_k(x,t) = k^2 \theta(kx, k^2 t), d_k(x,t) = d(kx, k^2 t)$ 下保持不变。
   • 构造尺度不变泛函： 作者精心构造了 $S_t$（随时间演化的量）和 $S_0$（初始量）。
     $$ S_t \triangleq \left(1 + \bar{\rho}_t + \frac{1}{\bar{\rho}_t}\right) S'_t S''_t $$
     其中 $\bar{\rho}_t = \|\rho(t)\|_{L^\infty}$。证明的核心在于建立 $S_t$ 与初始 $S_0$ 之间的直接不等式关系。

3. Bootstrap 论证（Bootstrap Argument）：

   • 策略： 作者定义了一个组合能量泛函 $B(t)$（见 Proposition 2.1），它包含了密度的 $L^\infty$ 范数、温度耗散积分、以及高阶速度/导向场的非线性积分项。
   • 步骤：
     1. 假设在时间区间 $[0,T]$ 上 $B(t) \le 2\epsilon$。
     2. 利用 Lemmas 2.2-2.6 中的估计，推导出 $B(t) \le \frac{3}{2}\epsilon$（即 $B(t) \le \epsilon$）。
     3. 这表明只要 $\epsilon$ 足够小，解的范数就不会增长，从而可以延拓到全局。

4. 复杂的非线性项处理：

   • 在温度方程估计中，出现了形如 $\|\nabla u\|_{L^2}^2 \|\rho\|_{L^3}$ 的项。为了控制它，作者将其与 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 联系起来，并利用了 $F$ 和 $\text{curl} u$ 的椭圆估计。
   • 导向场项极其复杂，例如 $|\Delta d + |\nabla d|^2 d|^2$。作者利用几何恒等式 $|\Delta d + |\nabla d|^2 d|^2 = |\Delta d|^2 - |\nabla d|^4$，并结合 Holder 不等式和 Sobolev 嵌入，将其转化为涉及 $\|\nabla^2 d\|_{L^2}, \|\nabla^3 d\|_{L^2}, \|\nabla d\|_{L^4}$ 的项，最终归入 $S''$ 中进行控制。

4. 与之前工作的比较

本文在 Commun. Math. Sci. 21 (2023), 1455–1486（由本文第二作者参与完成）的基础上做出了显著的改进，主要体现在移除粘性系数限制和引入新的更广泛的尺度不变量两个方面。

a) 放宽的假设条件

1. 移除粘性系数限制 $2\mu > \lambda$：
   • 先前工作： 文中引用的 [17]（即之前的那篇工作）为了获得全局解，必须要求 $2\mu > \lambda$。这是一个数学上的技术限制，用于保证动量方程中压力项的正定性以及获得速度的高阶能量估计。
   • 本文改进： 本文的 Theorem 1.1 完全没有对 $\lambda$ 的大小做任何额外假设，只要满足物理条件 $2\mu + 3\lambda \ge 0$ 即可。这使得结论更符合物理实际，因为许多流体的体积粘性系数 $\lambda$ 远大于剪切粘性 $\mu$。
   • 技术原因： 这一改进归功于作者对压力 $P$ 和函数 $F$ 的精细处理。通过将 $P$ 与 $\text{div} u$ 分离，并利用 $L^3$ 范数的估计（Lemma 2.3），作者避免了在估计 $\int \|\nabla u\|_{L^2}^2$ 时对系数矩阵正定性的依赖。

b) 结果的强度与新颖性

1. 新的尺度不变量：

   • 先前工作： [17] 定义的尺度不变量为 $N = \hat{\rho} [\dots] [\dots]$，其中 $\hat{\rho} = \|\rho_0\|_{L^\infty} + 1$。
   • 本文改进： 本文定义了一个全新的量 $S_0 = (1 + \bar{\rho} + \frac{1}{\bar{\rho}}) S' S''$。
   • 比较： 观察 $S_0$ 中的因子 $(1 + \bar{\rho} + \frac{1}{\bar{\rho}})$。在 [17] 中使用的是 $(1 + \bar{\rho})$。当 $\bar{\rho}$ 很大时，$(1 + \bar{\rho} + \frac{1}{\bar{\rho}})$ 与 $(1 + \bar{\rho})$ 表现类似；但当 $\bar{\rho}$ 趋于 0（真空区域较大）时，$\frac{1}{\bar{\rho}}$ 项变得很大。作者在文中提到，虽然这个形式看起来更复杂（Remark 1.2），但它是为了配合新的先验估计（特别是涉及 $\|\rho\|_{L^3}$ 的估计）而构造的。这个新量在同样满足尺度不变性的同时，允许作者在不假设 $2\mu > \lambda$ 的情况下关闭能量估计。

2. 包含更复杂的耦合项：

   • 在 $S''$ 的定义中，包含了 $\|\nabla d_0\|_{L^4}^4$。这直接对应于导向场方程中的高阶非线性项。之前的估计可能更多依赖于 $L^2$ 或 $L^\infty$ 范数，而本文直接将 $L^4$ 范数纳入尺度不变框架，这使得估计更加精细，能够处理更强的非线性耦合。

总结来说，这篇论文通过引入一个新的、包含 $1/\bar{\rho}$ 项的尺度不变量，配合对压力和温度方程的精细分析，成功解决了非等温可压缩向列液晶流中长期存在的粘性系数限制问题，推广了全局强解的存在唯一性理论。

3. 二维粘弹性 Giesekus 模型的收敛数值格式

Convergent numerical schemes for the viscoelastic Giesekus model in two dimensions

作者: Endre Süli, Dennis Trautwein

类别: math.NA, math.AP

发布日期: 2025-12-28

arXiv ID: 2512.22831v1

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

详细分析

中文标题: 二维粘弹性 Giesekus 模型的收敛数值格式

1. 研究对象和背景

本文主要研究二维空间中不可压缩粘弹性流体的 Giesekus 模型 (Giesekus model) 的数值逼近及其收敛性分析。该模型属于非牛顿流体力学的范畴，广泛应用于描述聚合物溶液、生物流体等复杂流体的行为。这些流体的显著特征是同时具备粘性（耗散能量）和弹性（存储能量），即所谓的粘弹性 (viscoelasticity)。与牛顿流体不同，粘弹性流体会表现出如韦森伯格效应 (Weissenberg effect, 即爬杆效应)、挤出胀大 (die swell) 和减阻 (drag reduction) 等独特的物理现象。

在宏观层面，描述这些流体的常用模型是在不可压缩 Navier-Stokes 方程的基础上引入一个额外的应力张量 (extra stress tensor) 来表征弹性效应。在 Giesekus 模型中，这个额外的应力通常与构象张量 (conformation tensor) $B$ 有关，后者描述了聚合物链的拉伸和取向。然而，本文作者选择使用弹性变形梯度 $F$ 来表述控制方程，即 $B = F F^\top$。这种表述方式能够利用更多的几何结构，有助于理论分析。

数学模型

设 $T > 0$ 为时间上限，$\Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$) 为具有 Lipschitz 边界的有界区域。我们需要求解函数 $(v, p, F): \Omega \times (0, T) \to \mathbb{R}^d \times \mathbb{R} \times \mathbb{R}^{d \times d}$，其中 $v$ 是流体速度，$p$ 是压力，$F$ 是弹性变形梯度。系统由以下方程组控制：

且要求满足约束条件：

这里 $\rho, \nu, \mu, \lambda > 0$ 分别是质量密度、溶剂粘度、弹性模量和弹性松弛时间。方程 (1.1a) 是动量方程，(1.1b) 是不可压缩条件，(1.1c) 是变形梯度的演化方程，其中包含对流项和非线性松弛项。

边值条件与函数空间

系统配备了齐次 Dirichlet 边界条件：

以及初始条件：

为了进行严格的数学分析，文章定义了弱解 (weak solution)。对于二维情况 ($d=2$)，弱解属于特定的 Sobolev 空间和 Bochner 空间：

• $v \in C([0,T]; L^2_{\text{div}}(\Omega; \mathbb{R}^2)) \cap L^2(0,T; H^1_{0,\text{div}}(\Omega; \mathbb{R}^2))$
• $F \in C([0,T]; L^2(\Omega; \mathbb{R}^{2\times 2})) \cap L^4(\Omega_T; \mathbb{R}^{2\times 2})$

其中 $\Omega_T = \Omega \times (0, T)$。弱解的定义满足变分形式（测试函数分别为 $w \in H^1_{0,\text{div}}(\Omega; \mathbb{R}^2)$ 和 $G \in H^1(\Omega; \mathbb{R}^{2\times 2})$）。

2. 主要定理或主要结果

本文的核心数学结果是证明了所提出的数值格式在二维情形下的收敛性，并由此给出了 Giesekus 模型全局弱解存在性的一个替代性证明。文章主要包含以下定理：

定理 3.5 (离散解的存在性与稳定性) 在二维或三维空间中，对于给定的离散初始数据 $v_{0,h} \in U_{h,k+1,\text{div}}$ 和 $F_{0,h} \in S_{h,m}$，存在离散解序列 $(v_h^n, p_h^n, F_h^n)$ 满足离散数值格式 (3.3)。此外，该离散解序列满足离散的能量不等式，即能量随时间衰减，且在固定时间步长 $\Delta t$ 和网格尺寸 $h$ 下保持一致有界。

定理 4.1 (空间离散极限 $h \to 0$) 固定时间步长 $\Delta t > 0$，当空间网格尺寸 $h \to 0$ 时，离散解序列 $(v_h^n, F_h^n)$ 的子序列收敛于连续时间、连续空间的函数 $(v^\Delta, F^\Delta)$，这些函数满足经过修正的、带有“人工应力扩散”项的连续时间方程。

定理 4.4 (整体收敛性，二维) 设 $\Omega \subset \mathbb{R}^2$ 为多边形 Lipschitz 区域，初始数据 $(v_0, F_0)$ 满足 $v_0 \in L^2_{\text{div}}(\Omega)$，$F_0 \in L^2(\Omega; \mathbb{R}^{2\times 2})$ 且 $\det(F_0) > 0$，$\ln \det(F_0) \in L^1(\Omega)$。在时间步长 $\Delta t$ 和空间网格尺寸 $h$ 充分小的限制下，数值方法产生的解序列 $(v_h^n, F_h^n)$ 存在子序列，当 $h \to 0$ 且随后 $\Delta t \to 0$ 时，收敛于 Giesekus 模型 (1.1) 的一个全局弱解 $(v, F)$。该弱解满足定义 2.1 的所有要求，包括正性条件 $\det(F) > 0$。

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

本文的研究方法结合了有限元方法 (Finite Element Method, FEM) 的离散化技术与现代非线性偏微分方程 (PDE) 的紧性理论。证明过程极其复杂，主要分为三个阶段：离散格式的构建与性质分析、空间极限 ($h \to 0$)、时间极限 ($\Delta t \to 0$)。

a) 核心分析学工具

1. 能量估计 (Energy Estimates)： 这是整个分析的基石。在连续层面，系统存在一个能量耗散律（公式 1.6），包含动能和弹性能。在离散层面，作者精心设计了数值格式 (3.3)，特别是动量方程中的数值对流项和变形梯度演化方程中的项 $c_h(\cdot, \cdot, \cdot)$，使得离散解满足类似于连续情形的能量不等式。这种离散能量稳定性保证了数值格式的长期稳定性，并为后续的先验估计提供了 $L^2$ 和 $L^p$ 范数控制。

2. Bochner 空间与紧性定理： 为了处理时间依赖的 PDE，证明中大量使用了 Bochner 空间理论，特别是 $L^p(0, T; X)$ 空间。核心工具是 Aubin-Lions 引理 (及其变体)。该引理用于从一致有界性（通常在能量空间中获得）推导出强收敛性。例如，在 $h \to 0$ 的极限过程中，需要处理非线性项 $\mu F F^\top$，这需要速度场 $v$ 和变形梯度 $F$ 的强收敛性。证明中利用了 $F$ 在某种正则性空间中的紧嵌入。

3. 对偶性方法与弱收敛： 在通过先验估计得到解序列的一致有界性后，作者使用弱收敛（或弱星收敛）来提取子序列的极限。通过取极限，非线性项可能会产生“缺陷项”（defect measures）。为了证明这些缺陷项为零（即极限满足原始方程），需要使用对偶性，即通过测试函数空间的密度来逼近。

b) 解决技术性难题的技巧

本文面临的最大挑战是 Giesekus 模型的高度非线性以及缺乏先验的强正则性（尤其是在没有应力扩散项的情况下）。作者采用了以下核心技巧来克服这些困难：

1. 两阶段极限过程 ($h \to 0$ 先于 $\Delta t \to 0$)： 这是一个关键的策略。如果同时取极限，很难控制非线性项。

   • 第一步 ($h \to 0$)：首先固定 $\Delta t$，证明空间离散解收敛到一组连续时间、空间的函数 $(v^\Delta, F^\Delta)$。这一步利用了有限元的离散紧性，但此时得到的方程中引入了一个稳定化项 $\Delta t \Delta F^\Delta$（应力扩散）。
   • 第二步 ($\Delta t \to 0$)：现在问题转化为研究一组连续时间、但时间离散（半离散）的解序列。此时可以利用 2022 年 Bulíček 等人的分析框架，证明该序列收敛于 Giesekus 模型的弱解。这种“先空间后时间”的顺序避免了直接在完全离散格式下处理极度复杂的非线性紧性。

2. 人工应力扩散 (Artificial Stress Diffusion) 的巧妙运用： 在离散格式 (3.3c) 中，作者引入了一个项 $\Delta t (\nabla F_h^n, \nabla G_h)_{L^2}$。这在物理上对应于应力扩散，但在数学上起到了至关重要的正则化作用。它提供了变形梯度 $F_h$ 在 $H^1$ 空间中的均匀估计。没有这一项，很难证明 $F$ 的紧性。作者在文中也讨论了这种扩散的系数可以随 $\Delta t$ 变化（例如 $\phi(\Delta t)$），只要 $\phi(0)=0$，这样就保证了当 $\Delta t \to 0$ 时，扩散效应消失，不干扰原方程的极限。

3. 数值对流项 $c_h$ 的三种构造与离散链式法则： 由于有限元空间的限制（例如，速度场是连续的，变形梯度可能是分片常数或线性的），直接离散 $(v \cdot \nabla)F$ 会破坏能量守恒律（如 $\int (v \cdot \nabla) |F|^2 = 0$），导致数值不稳定。作者为此设计了三种不同的离散对流项 $c_h(v, F, G)$：

   • 反对称形式 (3.8)：适用于高阶元。
   • 基于链式法则的形式 (3.10)：这是最精妙的部分。当 $F$ 是线性有限元时，$|F|^2$ 是二次的，不在速度的对偶空间中。作者定义了一组辅助量 $\Lambda_{i,j}(F)$，使得 $\sum \Lambda_{i,j} \partial_{x_j} G$ 能够模拟 $\partial_{x_j} G$ 且满足离散的链式法则 $\sum \Lambda_{i,j} \partial_{x_j} F = \partial_{x_i} I_h[|F|^2]$。这保证了 $c_h(v, F, F) = 0$，从而保护了能量结构。

4. Renormalization (重整化) 技术： 在证明 $F$ 强收敛并识别非线性极限（即证明 $\lim F_h \cdot F_h = (\lim F_h) \cdot (\lim F_h)$）时，使用了重整化方法。通过测试 $\det(F)$ 或 $|F|^2$ 的演化方程，分析差分 $\langle F_h^2 - (\lim F_h)^2, G \rangle$，并利用单调性或凸性证明这些项在极限下消失。这通常涉及到构造复杂的测试函数，并利用 Giesekus 模型中非线性松弛项 $FF^\top F - F$ 的特定结构。

5. 局部压力重构 (Local Pressure Reconstruction)： 为了处理不可压缩条件，文章利用了 Taylor-Hood 元素的 inf-sup 条件。在极限过程中，为了测试非散度自由的函数，作者使用了压力重构技术，这允许在不假设压力正则性极高的情况下，通过局部压力来处理动量方程中的压力梯度项。

4. 与之前工作的比较

本文在数值分析和理论推导两个方面都对现有文献做出了显著贡献，特别是在处理无应力扩散的 Giesekus 模型方面。

a) 之前的工作及其局限性

• 高 Weissenberg 数问题 (High Weissenberg Number Problem)：传统的数值方法（如纯拉格朗日法或简单的欧拉法）在 Weissenberg 数 $Wi$ 较大时往往失效，即所谓的“高 Weissenberg 数崩溃”。这通常归因于无法准确追踪聚合物构象的剧烈变化。
• 理论存在性结果的缺失：在本文所依托的理论工作（Bulíček et al., 2022）之前，对于无应力扩散的 Giesekus 模型，仅在二维情况下证明了全局弱解的存在性。三维情况及高阶正则性仍然是公开难题。
• 数值收敛性的缺乏：之前的许多数值研究依赖于人工应力扩散（即在方程中显式加入 $\Delta \sigma$ 项）来保证稳定性，但往往缺乏严格的收敛性证明，或者只能证明收敛到正则化后的方程的解，而不是原方程的解。此外，许多数值格式在能量守恒方面表现不佳。

b) 本文的改进与突破

1. 构造了无需“截断”或“正则化”的收敛格式： 本文最大的亮点是证明了数值解收敛到原模型的弱解，而不是某个带有额外人工粘性项的修正模型。许多现有的收敛性证明为了处理非线性项，不得不假设某种截断函数（cut-off）来保证解的有界性，或者必须在方程中加入 $\varepsilon \Delta F$ 项并令 $\varepsilon$ 保持为正。本文通过上述的两阶段极限过程，成功地在 $\Delta t \to 0$ 之前消除了人为引入的扩散，实现了真正的“无损”逼近。

2. 针对低正则性空间 ($L^4$ 范数) 的处理： 二维 Giesekus 模型的弱解仅能保证 $F$ 在 $L^4(\Omega_T)$ 中有界（相比 $L^2$ 是次优的）。这对数值格式的稳定性分析提出了极高要求。本文设计的格式能够兼容这种低正则性，并且证明了在这种正则性下的收敛性。这比那些需要 $H^1$ 或更高正则性的数值方法要强得多。

3. 全面的数值方法分类与分析： 文章不仅仅提出了一种方法，而是系统地分析了三种针对变形梯度 $F$ 的离散策略（反对称形式、高阶元近似、离散链式法则近似）。这种分类揭示了在有限元离散中保持数学结构（如能量守恒和几何结构）的重要性，为后续研究者提供了明确的设计指南。

4. 对偶性证明与存在性定理的补充： 本文的收敛性证明实际上重构了弱解的存在性。通过证明数值近似解的极限存在，且满足弱解的定义，文章为 Bulíček 等人的纯分析结果提供了一个构造性证明。这在计算流体力学（CFN）领域具有重要意义，因为它连接了理论分析与实际计算。

综上所述，这篇论文通过引入精细的有限元离散技术和两阶段极限策略，解决了二维 Giesekus 模型在无应力扩散假设下的数值收敛难题，填补了从严格理论分析到鲁棒数值实现之间的关键空白。

4. 随机等熵欧拉方程的全局鞅熵解

Global Martingale Entropy Solutions to the Stochastic Isentropic Euler Equations

作者: Gui-Qiang G. Chen, Feimin Huang, Danli Wang

类别: math.AP, math-ph, math.FA, math.PR

发布日期: 2025-12-27

arXiv ID: 2512.22719v1

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

详细分析

中文标题: 随机等熵欧拉方程的全局鞅熵解

研究对象和背景

本文研究的是带有随机外力项的一维等熵欧拉方程组（Stochastic Isentropic Euler Equations）。在流体力学中，欧拉方程描述了理想流体（忽略粘性）的运动。等熵假设意味着流体的熵保持不变，从而可以通过压力 $P(\rho)$ 仅作为密度 $\rho$ 的函数来描述系统。本文考虑的方程组如下：

其中 $\rho(x,t)$ 表示流体密度，$m(x,t) = \rho(x,t) u(x,t)$ 表示动量，$u$ 为速度。$W(t)$ 是柱状维纳过程（cylindrical Wiener process），代表随机外力。$\Phi$ 是噪声系数。该系统属于严格双曲型方程组，在 $\rho > 0$ 时具有 genuinely nonlinearity（本质非线性）。

背景与假设：

1. 压力律（Pressure Law）：论文处理的是一般压力律，这比通常的 $\gamma$-律（$P(\rho) = \kappa \rho^\gamma$）更广泛。压力函数 $P(\rho)$ 需满足特定的渐近行为：
   • 当 $\rho$ 较小时，表现为 $\rho^{\gamma_1}$ ($1 < \gamma_1 \le \gamma_2$)；
   • 当 $\rho$ 较大时，表现为 $\rho^{\gamma_2}$ ($\gamma_2 < 3/2$)。 这种假设允许压力函数在不同密度区间有不同的幂次增长，处理了真空（$\rho \to 0$）和高密度区域的行为。
2. 定义域与边界：空间域为全空间 $\mathbb{R}$。研究允许正的远场密度（positive far-field density），即当 $|x| \to \infty$ 时，$\rho(x,t) \to \rho_\infty > 0$，且速度 $u \to 0$。这与某些假设密度在无穷远处趋于零的工作不同，更符合某些物理场景。
3. 解的概念：由于欧拉方程激波的存在以及随机项的引入，经典的光滑解不存在。论文引入**全局鞅熵解（Global Martingale Entropy Solutions）**的概念。这种解在分布意义下满足方程，并且满足针对所有凸熵对 $(\eta, q)$ 的熵不等式（由随机 Itô 公式推导），其中的熵产生项是一个关于时间的鞅测度。
4. 函数空间与能量：解具有有限相对能量（finite relative-energy）。相对能量的定义基于远场状态 $(\rho_\infty, 0)$，形式为：
   $$ \eta^*(\rho, m) = \frac{m^2}{2\rho} + e^*(\rho, \rho_\infty) $$
   其中 $e^*$ 是相对内能。这意味着初始数据 $(\rho_0, m_0)$ 不需要具备高阶矩估计（higher moment estimates），仅需满足相对能量有限，这大大放宽了对初始数据的要求。

主要定理或主要结果

论文的核心成果包含两个主要定理，分别针对一般压力律和多激式（Polytropic）压力律。

1. 针对一般压力律（General Pressure Law）： 定理 2.1 和 2.2 建立了全局鞅熵解的存在性与紧性。

• 存在性：对于具有有限相对能量的初始数据 $(\rho_0, m_0)$，存在一个概率空间 $(\tilde{\Omega}, \tilde{\mathcal{F}}, \tilde{\mathbb{P}})$ 和定义在其上的随机过程 $(\tilde{\rho}, \tilde{m})$，它是随机抛物逼近方程（通过粘性消失法引入）的极限，且满足随机欧拉方程的分布意义下的方程及熵不等式。
• 紧性：对于上述解序列，存在某种紧性，使得我们可以从逼近解中提取收敛子列。

2. 针对多激式压力律（Polytropic Pressure Case, $P(\rho) = \kappa \rho^\gamma$）： 定理 7.1 和 7.2 给出了更精细的结果。

• 局部机械能不等式：对于任意 $\gamma > 1$，全局解满足局部机械能不等式（Local Mechanical Energy Inequality）。值得注意的是，初始数据仅需有限相对能量，而不需要像早期工作那样要求熵的高阶矩估计。
• 高阶相对能量保持：如果初始数据具有有限高阶相对能量，那么解也保持有限高阶相对能量。
• 熵不等式与紧性：利用高阶能量估计，证明了解对更广泛的凸熵对（包括机械能及其通量）满足熵不等式，并建立了随机欧拉系统解序列的紧性。

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

本文的证明过程极其复杂，融合了现代偏微分方程（PDE）理论与随机分析的深层技术。主要挑战在于随机项破坏了经典 PDE 中常用的 $L^\infty$ 估计（例如通过不变区域法获得），这使得传统的补偿紧性（Compensated Compactness）方法失效。

a) 核心分析学工具

1. 粘性消失法 (Vanishing Viscosity Method)： 论文构造了含粘性项 $\epsilon \partial_{xx} u$ 和随机项的抛物型逼近方程：

   $$ d\rho^\epsilon + \partial_x m^\epsilon dt = 0 $$
   $$ dm^\epsilon + \partial_x \left( \frac{(m^\epsilon)^2}{\rho^\epsilon} + P(\rho^\epsilon) \right) dt = \epsilon \partial_{xx} u^\epsilon dt + \Phi(\rho^\epsilon, m^\epsilon) dW(t) $$

   目标是证明当 $\epsilon \to 0$ 时，逼近解收敛于欧拉方程的解。

2. 能量估计与高阶可积性 (Energy Estimates & Higher Integrability)： 由于缺乏 $L^\infty$ 边界，必须通过能量估计获得 $L^p$ ($p < \infty$) 空间中的高阶可积性。文章利用 Itô 公式推导了密度 $\rho$ 和速度 $u$ 的 $L^p$ 范数估计。例如，对于密度，通常的守恒律推导 $L^p$ 守恒，但在随机情形下，需要处理噪声项带来的耗散或增长，通过适当的熵函数（如 $S(\rho) = \rho^\gamma$）来控制。

3. 概率空间中的紧性理论 (Compactness in Probability Spaces)： 由于解是随机的，不能直接使用确定性 PDE 中的弱收敛。论文使用了 Jakubowski-Skorokhod 表示定理（或其推广形式，因为某些空间不是度量空间）。该定理允许在新的概率空间上构造与原序列同分布的随机变量序列，使得该新序列几乎必然收敛。这是将 PDE 的紧性结果推广到随机环境的关键桥梁。

4. 相对能量 (Relative Energy)： 这是一个核心概念，用于克服一般压力律缺乏显式熵核（Entropy Kernel）的困难。相对能量泛函 $\mathcal{R}(\rho, m; \rho_\infty, 0)$ 提供了对解的稳定性的度量。通过控制相对能量，可以绕过对 $\rho$ 的 $L^\infty$ 控制，直接获得 $L^p$ 估计。

b) 技术性难题的解决技巧

1. 随机 $L^p$ 补偿紧性框架 (Stochastic $L^p$ Compensated Compactness Framework)： 这是本文最大的技术创新点。在确定性情形（Chen-Perepelitsa），$L^p$ 补偿紧性通过高阶可积性和 Murat 引理（或类似的 H-分布理论）来获得紧性。在随机情形下，作者开发了随机版本：

   • 困难：随机项破坏了方程的微分结构，使得直接应用补偿紧性变得困难。
   • 解决：首先利用抛物逼近的强解存在性（Theorem 3.1），通过截断技术（Cut-off technique）证明强解存在，从而可以对固定的空间变量 $x$ 应用有限维 Itô 公式。这保证了高阶能量估计的合法性。
   • 紧性传递：利用 Skorokhod 表示定理将弱收敛转化为几乎必然收敛，然后在新的概率空间上应用随机版本的 Div-Curl Lemma（广义旋度-散度引理）和 Murat 引理。

2. 处理缺乏显式熵核的困难 (General Pressure Law)：

   • 问题：对于多激式气体（$\gamma$-law），熵核 $\chi$ 有显式公式，允许精细计算。对于一般压力律，没有显式公式，且 $\rho \to 0$ 时方程退化。
   • 解决：作者引用并扩展了 Chen-Huang-Wang (2011) 中的技术，构造了特殊的非凸熵对。通过分析这些熵对的生成函数和奇异性质，推导出速度 $u$ 的高阶可积性（Higher integrability of velocity）。这是建立补偿紧性的基石，因为 $L^p$ 补偿紧性依赖于 $\nabla u \in L^2$ 和 $\rho$ 的高阶可积性。

3. Young 测度的约化 (Reduction of Young Measure)：

   • 问题：在极限过程中，非线性项 $\rho_\epsilon u_\epsilon$ 可能收敛到一个 Young 测度 $\nu_{x,t}$，而不是点态收敛。
   • 解决：利用随机版本的 Div-Curl Lemma 推导出 Tartar 交换关系（Commutation Relation）。由于缺乏 $L^\infty$ 边界，传统的“紧支集 Young 测度约化”方法（基于支撑有界）失效。作者利用建立的高阶可积性证明了 $\nu_{x,t}$ 是 $L^r$ 随机 Dirac 测度（即几乎处处点态收敛）。这一步依赖于对随机熵产生测度的**紧性（Tightness）**分析，特别是在负阶 Sobolev 空间中的紧性。

4. 随机近似解的一致估计 (Uniform Estimates for Stochastic Approximate Solutions)： 由于不能使用不变区域法（Invariant Region Method）获得 $L^\infty$ 估计，作者通过细致的能量估计来获得 $L^p$ 估计。

   • 技巧：为了处理非线性项中的 Itô 积分，作者证明了逼近方程具有全局强解（Theorem 3.1）。这使得可以对 $\rho^\epsilon$ 和 $u^\epsilon$ 进行逐点或高阶导数的 Itô 公式操作，进而通过 Grönwall 不等式类型的论证得到一致的 $L^p$ 上界。

5. 高阶相对能量与熵不等式： 文章证明了如果初始数据具有有限高阶相对能量，解也保持该性质。这使得可以对生成函数具有亚三次增长（sub-cubic growth）的凸熵对建立熵不等式。这对于证明解的唯一性（在某种意义下）和进一步的紧性至关重要。

与之前工作的比较

本文在多个维度上显著改进了现有文献，特别是对比 Berthelin-Vovelle (2021) 和 Chen-Perepelitsa (2009) 的工作。

a) 初始数据正则性要求的显著降低

• 先前工作 (Berthelin-Vovelle)：为了处理随机项导致的 $L^\infty$ 边界缺失，Berthelin-Vovelle 依赖于对熵生成函数进行高阶矩估计（higher moment estimates）。这要求初始数据必须满足对应熵的高阶矩（例如 $\mathbb{E}[\int \rho^p dx] < \infty$，其中 $p$ 较大）。这在物理上是苛刻的，因为它要求初始密度在无穷远处的衰减必须非常快。
• 本文改进：本文仅要求初始数据具有有限相对能量（Finite Relative Energy）。相对能量只涉及 $\rho$ 的某些幂次（与压力律有关）和动能，不涉及高阶矩。这意味着初始数据可以更粗糙，允许在无穷远处有较慢的衰减（只要密度趋于 $\rho_\infty > 0$）。作者通过建立逼近解的高阶相对能量估计（而不直接要求初始数据的高阶矩）实现了这一突破。

b) 压力律适用范围的推广

• 先前工作：大多数随机欧拉方程的研究（如 Berthelin-Vovelle）集中在多激式压力律 $P(\rho) = \kappa \rho^\gamma$，且主要针对 $\gamma$ 接近 1 或某些特定范围。确定性领域中，Chen-Le Floch 等人研究了广义压力律，但随机领域尚未覆盖。
• 本文改进：本文处理了一般压力律，允许压力函数在低密度和高密度区域有不同的幂次行为（$\gamma_1$ 和 $\gamma_2$）。这是对物理模型的极大贴近，因为真实气体在不同密度下表现不同。由于没有显式熵核，作者必须开发新的数学工具（如基于奇异性的分数阶导数分析和非凸熵对构造）来处理这一广义情形。

c) 解的性质与物理意义的增强

• 机械能不等式：在 Berthelin-Vovelle 的工作中，为了证明存在性，需要对熵进行高阶控制，这掩盖了物理能量的直接性质。本文证明了，即便在有限相对能量假设下，解依然满足局部机械能不等式。这是流体力学中的基本物理性质，意味着能量的耗散是由粘性逼近产生的，而在激波处满足 Rankine-Hugoniot 条件和熵条件。
• 紧性结果：本文不仅证明了单个解的存在性，还证明了解序列的紧性（Compactness of solutions）。这意味着从不同初始数据或不同逼近参数出发的解族具有某种连续性，这是建立适定性理论（Well-posedness）的基础。

d) 随机分析与 PDE 理论的深度融合

• 比较：之前的随机守恒律研究（如 Feng-Nualart 对标量方程）主要依赖于图-函数（Graph-transform）或随机紧性方法。对于方程组，这些方法失效。
• 本文贡献：本文将确定性 PDE 中强大的 $L^p$ 补偿紧性框架 成功移植到随机世界。这需要解决随机噪声带来的额外困难：
  1. 随机项作为扰动而非噪声：通过将随机项视为源项并在能量估计中控制它，利用 Itô 等式来获得耗散/增长界限。
  2. 随机 Young 测度的识别：利用 Skorokhod 表示定理后，识别新概率空间上的极限测度为 Dirac 测度。这需要结合 PDE 的 Tartar 交换关系和概率论的分布性质（Law of the variable）。

综上所述，这篇论文在随机流体力学领域迈出了重要一步，它将处理确定性非线性双曲方程组的复杂工具（如 $L^p$ 补偿紧性和广义相对能量）成功应用于随机环境，并极大地放宽了物理上的限制条件。

5. 具有趋化性和质量输运的不可压缩两相流的热力学一致的弥散界面模型研究

On a Thermodynamically Consistent Diffuse-Interface Model for Incompressible Two-Phase Flows with Chemotaxis and Mass Transport

作者: Andrea Giorgini, Jingning He, Hao Wu

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22585v1

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

详细分析

中文标题: 具有趋化性和质量输运的不可压缩两相流的热力学一致的弥散界面模型研究

研究对象和背景

该论文研究的是一个描述两相流体运动的流体动力学系统，具体为 Navier-Stokes/Cahn-Hilliard 类型的耦合方程组。该系统描述了两种密度不匹配（unmatched densities）的不可压缩粘性牛顿流体的混合流动，同时该混合流体中还含有一种可溶的化学物质。

论文的研究对象是定义在二维有界光滑区域 $\Omega$ 上的初边值问题。系统涉及四个主要的状态变量：

• $v(x,t)$：体积平均速度场；
• $P(x,t)$：流体混合物的压力；
• $\phi(x,t)$：两流体组分的体积分数差（相场变量），取值在 $[-1,1]$ 之间；
• $\sigma(x,t)$：化学物质的密度（或浓度），非负。

该模型基于 Onsager 变分原理 推导而来，具有热力学一致性。它包含了由化学物质引起的趋化效应（chemotaxis）以及混合物内的质量输运过程。界面是“弥散”的，即通过相场变量 $\phi$ 平滑过渡，参数 $\varepsilon$（在论文后续分析中常设为1，但在模型推导中代表界面厚度）控制界面的厚度。

论文在引言部分（Section 1.1）给出的完整数学方程组如下：

在论文的分析中，经过适当的重新解释压力项，方程组被重写为标准的对流-扩散形式（方程 2.8），这对于数值分析和存在性证明更为方便。该系统的边界条件为齐次 Neumann 边界条件（即无通量）：

• $v = 0$ on $\partial \Omega$ (无滑移边界)；
• $\partial_n \phi = m(\phi) \partial_n \mu = \sigma \partial_n w = 0$ on $\partial \Omega$。

初始条件为 $v|_{t=0}=v_0, \phi|_{t=0}=\phi_0, \sigma|_{t=0}=\sigma_0$。

该系统的核心物理机制在于：

1. 动量守恒：受变密度的 Navier-Stokes 方程控制，包含由于密度不匹配引起的相对通量 $J$。
2. 相分离：由 Cahn-Hilliard 方程控制，其中化学势 $\mu$ 包含双阱势 $\Psi$，保证 $\phi$ 维持在物理区间 $[-1,1]$。
3. 趋化与质量输运：化学物质的演化遵循对流-扩散方程，其中扩散项由 $w$ 的梯度驱动，而 $w$ 与 $\phi$ 耦合（通过 $\beta(\phi)$）。这种耦合导致了交叉扩散效应，类似于经典的 Keller-Segel 模型，但此处的扩散是由四阶 Cahn-Hilliard 动力学驱动的。

主要定理或主要结果

该论文的主要结果分为三个层次：弱解的存在性、强解的存在唯一性以及正则性的提升。所有结果均针对二维有界区域 $\Omega$。

1. 全局弱解与有限能量解的存在性 (Theorem 2.1)

在假设 (H1)-(H4) 下，对于给定的初始数据 $(v_0, \phi_0, \sigma_0)$，系统 (2.8) 满足以下性质：

• 如果初始能量有限，即 $E(v_0, \phi_0, \sigma_0) < \infty$，则存在一个全局有限能量解 $(v, \phi, \mu, \sigma)$。该解满足能量不等式：
  $$ \frac{d}{dt} E(v, \phi, \sigma) + \int_{\Omega} \left( 2\nu(\phi)|Dv|^2 + m(\phi)|\nabla \mu|^2 + \sigma |\nabla w|^2 \right) dx \le 0. $$
• 如果额外假设 $\sigma_0 \in L^2(\Omega)$，则存在一个全局弱解。

2. 全局强解的存在唯一性与正则性提升 (Theorem 2.3 & Corollary 5.1)

在更强的假设下（初始数据足够光滑，$v_0 \in V_{0,\sigma}$, $\phi_0 \in H^1(\Omega)$, $\sigma_0 \in L^2(\Omega) \cap L^\infty(\Omega)$ 等），论文证明了：

• 存在性与唯一性 (Theorem 2.2)：系统存在唯一的全局强解 $(v, \phi, \mu, \sigma)$。
• 正则性提升 (Theorem 2.3)：任意由 Theorem 2.1 得到的全局弱解，如果其初始数据满足强解的正则性要求，则该解在 $t>0$ 的任意时刻自动成为全局强解（即瞬时正则化）。
• 有界性 (Corollary 5.1)：对于全局强解，如果初始密度 $\sigma_0$ 是有界的（$\sigma_0 \in L^\infty(\Omega)$），则化学物质密度 $\sigma$ 对于所有时间 $t \ge 0$ 保持有界。这一点与经典的二维 Keller-Segel 系统有本质区别，后者可能出现质量集中奇点（blow-up）。

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

该论文的分析结合了现代偏微分方程理论中的多种高级工具，特别是在处理变密度流体、奇异势和交叉扩散时展示了复杂的技巧。

a) 核心分析学工具

1. 能量方法 (Energy Methods)： 这是分析的基石。作者利用系统固有的热力学结构（Onsager 变分原理导出的结构），推导出总能量 $E(v, \phi, \sigma)$ 的耗散律。这提供了全局 $L^2$ 类型的先验估计。关键的困难在于处理总能量中的相互作用项 $\int \beta(\phi)\sigma dx$。由于 $\beta$ 和 $\sigma$ 的符号不确定，作者利用了 $\phi \in [-1,1]$ 的物理约束（由奇异势 $\Psi$ 保证）以及熵项 $\sigma \ln \sigma$ 的凸性，通过 Young 不等式等手段控制了该项。

2. Galerkin 逼近与紧性方法 (Semi-Galerkin Scheme & Compactness Methods)： 为了构造弱解，论文采用了半 Galerkin 逼近方案。具体而言，对于 Navier-Stokes/Cahn-Hilliard 子系统，利用有限维子空间逼近速度场和相场变量。通过 Schauder 不动点定理证明局部解的存在，然后结合能量估计延拓为全局解。在取极限 $n \to \infty$ 时，利用了经典紧嵌入定理（如 Rellich-Kondrachov 定理）和弱收敛技术。为了通过非线性项（如 $\rho(\phi)v \otimes v$），需要利用能量耗散律导出的更高阶估计来获得更强的紧性，例如利用 Aubin-Lions 引理。

3. Stokes 算子与 Helmholtz-Weyl 分解： 在处理不可压缩 Navier-Stokes 部分时，论文使用了 Helmholtz-Weyl 分解将速度场分解为无散部分和梯度部分，利用 Leray 投影 $P$ 将方程投影到无散空间 $L^2_{0,\sigma}(\Omega)$ 上。这允许将动量方程转化为关于速度场的演化方程，其中压力项被消除。

b) 解决技术性难题的技巧

1. 奇异势的处理与正规化 (Approximation of Singular Potential)： 由于 Cahn-Hilliard 方程中的双阱势 $\Psi(\phi)$ 在 $\phi = \pm 1$ 处趋于无穷大（奇异），直接处理非常困难。论文采用了正规化策略：

   • 对 $\Psi$ 进行多项式型的正规化（截断），使其在 $\mathbb{R}$ 上光滑，从而可以利用标准的抛物理论处理正规化系统。
   • 关键修正：由于正规化后的相场变量 $\phi$ 可能不再严格位于 $[-1,1]$ 内，导致物理参数（如密度 $\rho(\phi)$）可能取到负值（无物理意义）。为了解决这个问题，论文引入了一个非线性的密度扩展 $\hat{\rho}(\phi)$ 和非线性修正项 $\hat{R}\hat{v}/2$ 到 Navier-Stokes 系统中。这些修正项保证了正规化系统的能量耗散结构，从而使得在取极限时能够恢复原始系统的性质。

2. 分解策略 (Decomposition Strategy)： 由于系统 (2.8) 具有高度非线性耦合（特别是 $\sigma$ 方程中的 $\text{div}(\beta'(\phi)\sigma\nabla\phi)$ 项），直接构造强解极其困难。作者采取了分解-重构的策略：

   • 将系统分解为三个解耦的子系统：速度场子系统（带漂移的 Stokes 方程）、相场子系统（带漂移的 Cahn-Hilliard 方程）和化学物质子系统（带漂移的扩散方程）。
   • 论文在 Section 5 中单独分析了这些子系统的正则性。例如，证明了如果漂移项 $v$ 和 $\phi$ 具有足够的正则性，那么 $\sigma$ 的方程具有全局解且保持正则性。

3. Bootstrap 论证 (Bootstrap Argument)： 这是证明强解存在性和正则性提升的核心技术。

   • 从全局弱解出发，利用弱解已经具有的基本能量估计。
   • 通过分解后的子系统分析，获得更高阶的正则性反馈。例如，利用 $\sigma$ 的 $L^\infty$ 估计（若 $\sigma_0 \in L^\infty$）或者对 $\nabla \sigma$ 的 $L^q$ 估计，结合耦合项 $\beta'(\phi)\sigma\nabla\phi$ 的结构，反过来提升 $\phi$ 和 $v$ 的正则性，然后再进一步提升 $\sigma$ 的正则性。
   • 这种迭代过程最终证明了时间 $t>0$ 后解的光滑性。

4. 交叉扩散项的控制 (Controlling Cross-Diffusion)： 化学物质 $\sigma$ 方程中的交叉扩散项 $-\text{div}(\beta'(\phi)\sigma\nabla\phi)$ 是导致 Keller-Segel 类型奇点产生的根源。论文证明的关键在于，此处的扩散并非由 Laplace 算子直接驱动，而是由 Cahn-Hilliard 算子 $-\Delta \phi + ...$ 驱动。Cahn-Hilliard 方程提供了更好的空间正则性（相对于二阶抛物方程），这使得 $\nabla \phi$ 更光滑，从而有效地抑制了 $\sigma$ 的爆炸性增长。

与之前工作的比较

该论文在前人工作的基础上做出了显著的改进，特别是在模型的一般性、无需额外正则化假设以及揭示与经典模型的区别方面。

1. 模型设置的推广 (Generalization of Model Setting)

• 变密度 (Unmatched Densities)：许多先前关于流体-趋化耦合的工作（如 [40]）假设流体是恒密度的（$\rho_1 = \rho_2$）。这简化了动量方程中的对流项。本文处理了密度不匹配的情况，这引入了额外的相对通量项 $J = -\frac{\rho_1-\rho_2}{2} m(\phi) \nabla \mu$。这在数学上处理起来更困难，因为它导致动量方程中的密度和速度是非对称耦合的。
• 非线性系数与一般相互作用 (General Coefficients)：论文假设粘度 $\nu(\phi)$ 和迁移率 $m(\phi)$ 是 Lipschitz 连续的有界函数，不需要是常数或简单的线性函数。更重要的是，相互作用函数 $\beta(\phi)$ 被假设为非线性且有界的（只要求 $\beta \in C^2(\mathbb{R})$ 且在区间外延拓为0），这比之前许多工作（如 [40, 62]）中假设的线性形式 $\beta(\phi) = \chi(1-\phi)$ 或 $\chi \phi$ 更为一般。这种非线性 $\beta$ 对应于更复杂的物理相互作用，且在数学上要求更精细的估计来控制 $\sigma$ 方程中的交叉扩散项。

2. 摆脱正则化假设 (Absence of Regularization Terms)

这是本文相对于 [40, 62, 64] 等工作的一个重大突破。

• 无 logistic 源项：经典的趋化模型（Keller-Segel）及其流体耦合变体通常依赖于 logistic 源项（如 $-\sigma^2$）来控制浓度的增长，防止有限时间爆破。本文的模型没有这种源项，完全依赖于系统固有的扩散机制。
• 无退化灵敏度：一些工作（如 [64]）使用退化的趋化灵敏度函数（如 $\frac{\sigma}{1+\sigma^q}$）来正则化模型。本文使用的是标准的线性灵敏度结构（在 $w = \ln \sigma + \beta(\phi)$ 中体现为 $\sigma \nabla w$），没有引入人为的退化。
• 无截断：在处理 $\sigma$ 时，本文没有像 [62] 那样使用截断技术，而是利用强最大原理和近似解理论保持 $\sigma$ 的严格正性。

3. 结果的强度与物理洞察 (Strength of Results and Physical Insight)

• 全局有界性 (Global Boundedness)：论文证明了在二维情况下，如果初始数据有界，则解 $\sigma$ 保持有界。这与经典的二维 Keller-Segel 系统形成鲜明对比。在经典的 KS 系统中，存在一个“临界质量”现象：当质量超过某个阈值时，解会在有限时间内爆破（形成奇点）。然而，在本文研究的两相流模型中，由于四阶 Cahn-Hilliard 动力学的引入，即使没有 logistic 项，$\sigma$ 的扩散也是由化学势 $\mu$ 的梯度驱动的，而 $\mu$ 包含 $-\Delta \phi$。这种高阶扩散机制有效地抑制了聚集，防止了爆破。
• 正则性提升 (Propagation of Regularity)：论文证明了弱解的“瞬时正则化”性质，即从较低正则性的初始数据出发，解在 $t>0$ 后立即变得光滑。这一结果的证明依赖于对分解子系统的细致分析和 Bootstrap 论证，比之前仅获得弱解或局部强解的结果要强。

总结来说，该论文通过引入更贴近物理的假设（变密度、一般非线性相互作用）并结合精细的分析技巧（分解策略、Bootstrap、正规化修正），在没有人工正则化项的情况下，完整刻画了二维不可压缩两相流趋化系统的全局适定性，并揭示了该系统独特的正则化机制。

6. 光滑欧拉流中涡破碎现象的波动动力学机制

Wave dynamics governing vortex breakdown in smooth Euler flows

作者: Tsuyoshi Yoneda

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22543v1

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

详细分析

中文标题: 光滑欧拉流中涡破碎现象的波动动力学机制

研究对象和背景

本论文 Tsuyoshi Yoneda 的研究对象是三维不可压缩欧拉方程（3D incompressible Euler equations）所描述的流体动力学系统，特别是该系统中**涡破碎（vortex breakdown）**这一现象的数学机制。涡破碎是流体力学中一种重要的不稳定性现象，通常表现为涡核处速度轴的停滞、回流区的形成或者螺旋结构的出现，它在高雷诺数湍流中起着关键作用。

论文聚焦于一个特定的物理场景：一个涡旋（vortex）被一个预先给定的拉格朗日流（prescribed Lagrangian flow）所输运。这意味着背景流场已经存在，而研究的主体是叠加在该背景流上的涡旋的演化。作者试图回答的核心问题是：在没有奇异涡核模型（如经典的涡丝模型）的情况下，光滑的欧拉方程中是否自然涌现出描述涡动力学的波动方程？

作者将速度场表示为流映射 $\Phi(t, x)$ 的时间导数，即 $v = \partial_t \Phi$。在拉格朗日形式下，三维不可压缩欧拉方程可以写为（隐含了体积保持约束）：

为了精细地描述涡旋轴附近的动力学，作者引入了一套动态的费涅尔-塞雷（Frenet-Serret）标架。令 $\ell \subset \mathbb{R}^3$ 表示初始时刻的涡轴曲线，$\Phi(t, \ell)$ 为输运后的曲线。定义弧长参数 $z$，使得切向量 $\tau = \partial_z \Phi$。基于 Frenet-Serret 公式，引入曲率 $\kappa(z)$、挠率 $T(z)$ 以及相应的法向量 $n(z)$ 和副法向量 $b(z)$。

对于涡轴附近的粒子，其位置可以展开为：

其中 $Z$ 描述了沿轴向的位移，$R_1, R_2$ 描述了横向位移。作者重点研究了线性化后的横向运动，引入参数 $\alpha_1(t)$ 和 $\alpha_2(t)$ 来描述涡轴的旋转和弯曲动力学。

主要定理或主要结果

论文的主要结果定理 1 (Theorem 1) 建立了描述涡轴几何特性的时间演化方程，揭示了涡破碎现象与波动动力学之间的直接联系。

定理 1 (Theorem 1): 支配涡破碎发生的量 $\alpha_1$ 和 $\alpha_2$ 满足如下的二阶波动方程：

这里，$v(t) = |\partial_t \Phi(t, \ell(s))|$ 是沿涡轴的切向速度，$\kappa(t)$ 是涡轴的曲率。

该定理表明，在光滑欧拉流中，描述涡轴偏离行为的量 $\alpha_1, \alpha_2$ 遵循线性波动方程，其“波速”由切向速度 $v(t)$ 及其二阶导数决定，且 $\alpha_1$ 的方程中含有与背景曲率 $\kappa$ 及其导数相关的非齐次强迫项。这一结果在数学上严格证明了即使在没有奇异假设的光滑流体中，涡动力学也具有本质的波动特征，这些波动方程控制着涡轴的演化，进而决定了涡破碎是否发生。

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

本论文的证明风格偏向于几何分析与形式化的微分几何计算，而非传统的偏微分方程估计理论。其核心思想是通过构造特殊的拉格朗日坐标系来消除压力项的奇异性，从而将欧拉方程的动力学转化为几何量的演化方程。

a) 核心分析学工具

1. Frenet-Serret 标架与微分几何计算： 论文的基础工具是经典微分几何中的 Frenet-Serret 公式。作者将三维欧拉方程的动力学限制在由输运涡轴生成的曲面上，并沿用 Frenet 标架（切向量 $\tau$，法向量 $n$，副法向量 $b$）来分解向量场。这使得复杂的流体运动可以分解为沿轴向的运动和绕轴向的旋转运动。 Frenet-Serret 公式为：

   $$ \begin{pmatrix} \tau' \\ n' \\ b' \end{pmatrix} = \begin{pmatrix} 0 & \kappa & 0 \\ -\kappa & 0 & T \\ 0 & -T & 0 \end{pmatrix} \begin{pmatrix} \tau \\ n \\ b \end{pmatrix} $$

   其中 $'$ 表示对弧长 $z$ 的导数。这些公式被反复用于计算 $\partial_z \Phi$, $\partial_z^2 \Phi$ 等高阶导数在各基底上的投影。

2. 拉格朗日坐标与逆函数定理： 为了处理涡轴附近的邻域，作者使用了隐函数定理（Inverse Function Theorem）。通过将粒子的初始坐标 $(r_1, r_2)$ 映射为当前时刻的横向坐标 $(R_1, R_2)$，作者构造了一个有效的拉格朗日坐标系 $(s, R_1, R_2)$。这使得对流场的微分运算（如 $\partial_{R_1}, \partial_{R_2}$）变得定义良好，且能反映流体微元的变形历史。

3. Lie Bracket (李括号) 为零的性质： 这是证明中最精妙的分析工具。在坐标变换 $(r_1, r_2) \to (R_1, R_2)$ 下，导数算子 $\partial_{R_1}$ 和 $\partial_{R_2}$ 代表了在时刻 $t$ 固定的情况下的空间切向方向。而 $\partial_z$ 沿着涡轴。作者指出，在选定的有效拉格朗日坐标系中，算子 $\partial_z$ 与 $\partial_{R_1}, \partial_{R_2}$ 是可交换的，即李括号为零：

   $$ [\partial_z, \partial_{R_1}] = \partial_z \partial_{R_1} - \partial_{R_1} \partial_z = 0. $$

   这一性质保证了在计算高阶混合导数时，可以自由交换求导顺序，极大地简化了压力梯度项的计算。

b) 解决技术性难题的主要技巧

1. 规避压力项的奇异积分表示： 传统的涡丝理论（如 Biot-Savart 定律或局部诱导近似 LIA）通常需要计算奇异积分来获得压力梯度或诱导速度，这要求涡量高度集中（即涡核半径趋于 0）。本论文的作者避免了这种做法。他不直接求解压力的积分表达式，而是将压力梯度项 $\nabla p$ 通过欧拉方程 $\partial_t^2 \Phi = -\nabla p$ 直接替换为加速度项。 技巧在于：作者计算 $\nabla p$ 在 Frenet 标架上的投影（例如 $\nabla p \cdot \tau$），但利用坐标变换将其表达为 $R_1, R_2$ 的导数。通过利用上述的“可交换性”，作者将对压力项的求导转化为对加速度项（即 $\partial_t^2 \Phi$）的几何投影计算。这成功地在不涉及奇异积分的情况下，解耦了压力场对涡轴动力学的影响。

2. 导数的层级剥离与闭合系统： 证明的核心在于对高阶项的处理。作者计算了 $\partial_z \partial_t^2 \Phi$ 等混合导数。 通过计算两条路径，作者得到了“闭合方程”：

   • 路径一（几何展开）：直接对 $\Phi$ 的几何展开式求导，利用 Frenet 公式得到 $\partial_{R_1} \partial_z \Phi$ 和 $\partial_{R_1} \partial_z^2 \Phi$ 的表达式，进而得到 $\partial_{R_1} (\partial_t^2 \Phi \cdot \tau)$。
   • 路径二（压力投影）：利用欧拉方程 $\partial_t^2 \Phi = -\nabla p$，计算 $\partial_{R_1} (\nabla p \cdot \tau)$，然后利用坐标变换和可交换性展开。 作者证明了这两条路径的结果必须相等。在比较这两者时，$O(R_1, R_2)$ 的低阶项相互抵消，留下了关于 $\alpha_1, \alpha_2$ 的二阶常微分方程（ODE）或波动方程。特别是，压力项的贡献与几何加速度项的贡献在 $R_1, R_2 \to 0$ 的极限下精确吻合，从而推导出了定理中的方程。

3. 处理对流项与非线性的几何化： 虽然欧拉方程本质上是非线性的，但在涡轴附近的局部线性化分析中，非线性效应被转化为背景流场几何量（曲率 $\kappa$、挠率 $T$）的演化。论文成功地将非线性的对流效应剥离为背景参数 $v(t)$ 和 $\kappa(t)$ 的函数，从而使得 $\alpha_1, \alpha_2$ 满足的方程呈现为（变系数的）线性波动方程。这是一种将非线性 PDE 约化为低维动力学系统的降维技巧。

与之前工作的比较

本论文的主要贡献在于将通常存在于奇异模型中的“波动子结构”推广到了光滑的欧拉框架中，填补了理论上的一个空白。

1. 突破了“奇异模型”的限制：

   • 前人工作：在流体力学和数学物理中，关于涡丝波动性质的经典结果大多基于奇异近似。最著名的是 H. Hasimoto (1972) 的发现，他证明了在局部诱导近似（LIA）和 Biot-Savart 定律下，涡丝的包络演化满足非线性薛定谔方程（NLS）。此外，许多关于涡动力学的研究都假设涡量集中在细丝上，或者使用涡核半径趋于零的渐近展开。
   • 本文改进：Yoneda 指出，这些奇异模型虽然直观且在某些物理极限下有效，但它们是否是欧拉方程的固有性质一直不清楚。本文证明了不需要奇异假设，在完全光滑的设定下，波动动力学依然存在。通过构造精确的拉格朗日坐标系，作者证明了波动方程是从原始的光滑欧拉方程中严格导出的，而不是奇异逼近的产物。

2. 对“涡破碎”机制的数学解释：

   • 前人工作：涡破碎通常被视为一种分岔现象，往往需要复杂的数值模拟或实验观察来描述。理论分析多集中在寻找临界雷诺数或特定的流动不稳定性。
   • 本文改进：本文提供了一个明确的数学判据。涡破碎定义为涡轴 $\zeta^*(t)$ 与由波动方程决定的轴 $\zeta(t)$ 之间的角度发散。由于 $\alpha_1, \alpha_2$ 满足波动方程，这意味着涡轴的演化具有传播性和干涉性。涡破碎的发生被归结为波动方程解的动力学性质，从而从纯分析的角度解释了涡核处的几何拓扑变化。

3. 方法论的创新：

   • 前人工作：传统上处理这类问题要么依赖于流体静力学近似，要么依赖于谱方法或能量估计来处理 Navier-Stokes 或 Euler 方程。
   • 本文改进：本文采用了一种基于几何李群结构的方法。通过强制在有效坐标系中 $[\partial_z, \partial_{R_i}] = 0$，作者利用微分同胚的性质消除了坐标系带来的复杂性。这种“无旋坐标系”或“随动坐标系”的构造，是将欧拉方程的复杂性转化为几何量演化的关键。这与之前基于 Fourier 估计或能量方法的工作截然不同，展示了微分几何在流体 PDE 分析中的强大能力。

总结来说，这篇论文在一个看似经典的领域（3D Euler 方程）中，通过精细的微分几何构造，发现并严格证明了涡动力学中隐含的波动机制，将其从奇异模型的特例提升为光滑流体的普适性质。

7. 具有壁面涡流边界条件的Navier-Stokes-αβ方程解的正则性

Regularity of solutions of the Navier-Stokes-αβ equations with wall-eddy boundary conditions

作者: Nella Rotundo, Gantumur Tsogtgerel

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22436v1

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

详细分析

中文标题: 具有壁面涡流边界条件的Navier-Stokes-αβ方程解的正则性

研究对象和背景

本文研究的对象是Navier-Stokes-αβ (NS-αβ) 方程组，这是一类描述不可压缩流体运动的正则化模型。该模型由Nella Rotundo和Gantumur Tsogtgerel提出，旨在解决经典α模型在处理**近壁湍流（near-wall turbulence）**物理机制时的不足。与传统的Navier-Stokes方程或仅引入单一过滤尺度的α模型（如Leray-α或Navier-Stokes-α）不同，NS-αβ方程组引入了两个长度尺度：一个是体过滤尺度 $\alpha$，另一个是壁面涡流尺度 $\beta$。

该研究的核心创新在于采用了由Fried和Gurtin (2008, 2009) 提出的壁面涡流边界条件（wall-eddy boundary conditions）。这种边界条件不仅包含经典的无滑移条件（$u=0$），还引入了正比于壁面涡量（wall vorticity）的切向涡量牵引力，从而在连续介质力学的框架下为近壁湍流提供了一个更合理的数学模型。

该论文主要关注该系统的稳态线性问题及非线性演化问题的全局适定性和正则性。研究区域为三维有界光滑区域 $\Omega \subset \mathbb{R}^3$。

1. 基本方程组

该模型的演化方程组（Navier-Stokes-αβ方程）如下：

其中，$v$ 是过滤后的速度场，$u$ 是流体的实际速度场，$p$ 是压力，$\alpha$ 和 $\beta$ 是常数，满足 $\alpha > \beta > 0$。速度场之间的关系由下式给出：

2. 边界条件

系统耦合了Fried和Gurtin提出的壁面涡流边界条件：

以及（针对涡量 $\omega = \nabla \times u$）：

这里，$n$ 是 $\partial\Omega$ 的单位外法向量，$\gamma$ 和 $\ell$ 是参数，满足 $|\gamma| \le 1$ 且 $\ell > 0$。

3. 函数空间与变分结构

论文首先在无散度向量场空间中建立了问题的变分形式。令 $V$ 为在 $H^1(\Omega)^3$ 中 $\nabla \cdot u = 0$ 且 $u=0$ 于 $\partial\Omega$ 的向量场的闭包。论文定义了一个双线性形式 $a(u, \phi)$：

该形式同时编码了体内的耗散项和壁面涡流的边界贡献。该研究的理论基础建立在对该双线性形式的性质分析之上，包括其连续性、对称性以及关键的Gårding不等式（Gårding inequality）。

主要定理或主要结果

论文的主要贡献在于证明了该非线性系统的全局适定性（global well-posedness）和正则性（regularity）。以下是论文中确立的核心结果：

定理 1 (稳态线性问题的正则性与Gårding不等式):

1. 对称性与连续性：双线性形式 $a(\cdot, \cdot)$ 在 $H^2(\Omega)^3$ 空间上是连续且对称的。
2. Gårding不等式：存在常数 $c_0 > 0$ 和 $\gamma_0 \ge 0$，使得对于所有 $u \in V_2$（即 $H^2$ 中的无散度向量场），满足：
   $$ a(u, u) + \gamma_0 \|u\|_{L^2}^2 \ge c_0 \|u\|_{H^2}^2 $$
   这表明该形式控制了 $H^2$ 范数，是后续建立算子理论的基础。
3. 椭圆算子理论：存在唯一的自伴算子 $A$，使得 $Au = f$ 等价于稳态系统的弱形式。利用谱理论，建立了相关的分数幂空间结构。

定理 2 (线性系统的全正则性): 在 $\partial\Omega \in C^4$ 及右端项 $f \in L^2$ 的假设下，通过验证 Douglis-Nirenberg 椭圆性 和 Lopatinskii-Shapiro 覆盖条件，利用 Agmon-Douglis-Nirenberg (ADN) 理论 证明了该线性系统的全正则性。具体来说，若 $u \in V_2$ 是弱解，则 $u \in H^4(\Omega)^3$ 且压力 $p \in H^3(\Omega)$，并满足强解意义下的边界条件。

定理 3 (非线性演化方程的全局适定性与正则性): 对于非线性演化问题，论文建立了全局存在性、唯一性、稳定性和正则性。 这通过建立**能量估计的层级（hierarchy of energy estimates）**来实现。具体结论包括：

• 存在全局解 $u(t)$。
• 解具有高阶正则性，即对于任意时间 $t > 0$，解保持在高阶 Sobolev 空间中（即 $u(t) \in H^s$ for high $s$）。
• 解是唯一的，且对初始数据具有连续依赖性。

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

该论文采用了现代偏微分方程理论中的高级分析工具，结合变分方法、算子半群理论和ADN椭圆正则性理论，系统地解决了NS-αβ系统的正则性问题。

a) 核心分析学工具

1. 变分形式与算子半群理论：

   • 论文首先放弃了传统的直接处理PDE的方法，转而研究由双线性形式 $a(u, \phi)$ 定义的变分问题。这种方法的优势在于可以利用Hilbert空间中的算子理论。
   • 通过Gårding不等式，作者构造了与 $a(\cdot, \cdot)$ 相关的自伴算子 $A$。这使得可以将稳态问题转化为算子方程 $Au = f$。
   • 对于演化方程 $\partial_t u + A u = \text{Nonlinear}(u)$，论文利用了线性算子生成的解析半群理论。通过Hille-Yosida定理或Lumer-Phillips定理，确保了线性部分的适定性，进而为非线性扰动分析提供了基础。

2. Agmon-Douglis-Nirenberg (ADN) 正则性理论：

   • 这是处理高阶椭圆系统（特别是强耦合系统）的标准工具。由于NS-αβ系统包含四阶导数（来自 $\Delta^2 u$）以及压力项，且边界条件混合了涡量及其梯度，传统的Schauder估计或$L^p$理论难以直接应用。
   • 论文的关键步骤是将系统写成ADN理论适用的形式（Equation 69）。这要求仔细标记每个方程和未知量的“阶数（weight）”以确保系统的椭圆性。
   • 为了应用ADN理论，必须验证两个技术性条件：
     • Douglis-Nirenberg 椭圆性：需要对系统的主象征（principal symbol）进行详细计算，确保其满足椭圆性条件。对于变系数（由边界几何形状引起的），这涉及到复杂的矩阵分析。
     • Lopatinskii-Shapiro 条件：这是边值问题正则性的核心。它要求线性化后的系统在边界上没有非平凡的衰减解。对于壁面涡流边界条件，由于其涉及切向导数和涡量关系，该条件的验证非常困难。论文通过对边界算子的符号分析和反射法完成了这一验证。

3. 先验能量估计（A Priori Energy Estimates）：

   • 对于非线性演化方程，论文建立了一个**“先验估计的层级”**。这是一种通过数学归纳法形式展示的估计过程。
   • 证明逻辑如下：
     • 假设解已经达到了一定的正则性（例如 $H^k$）。
     • 利用双线性形式的对称性和Gårding不等式，推导出关于高阶范数的微分不等式。
     • 使用Gronwall不等式证明这些高阶范数可以被初始数据和时间区间控制。
     • 由于估计不依赖于时间区间的长度，从而证明了全局存在性和正则性。

b) 解决技术性难题的技巧

1. 处理复杂的边界条件（Wall-Eddy BCs）：

   • 难点：标准的Navier-Stokes方程通常只处理 $u=0$ 或 $u=0$ 结合应力条件。Fried-Gurtin的边界条件引入了涡量的切向牵引力 $\beta^2 (1 - n \otimes n) \nabla \omega + \gamma (\nabla \omega)^\top n = \ell \omega$。这使得边界算子极其复杂，混合了二阶导数，且耦合了法向和切向分量。
   • 技巧：
     • 迹定理与密度论证：在Theorem 2.2中，作者证明了弱解实际上满足强边界条件。这通过证明边界积分项 $\int (kn \times \omega - n \times G n) \cdot \partial_n \phi \, dS = 0$ 对所有无散度测试函数 $\phi$ 成立，从而利用Sobolev空间中迹算子的稠密性得出边界算子必须消失。这巧妙地绕过了直接处理高阶迹的困难。
     • 几何分解：在验证Lopatinskii-Shapiro条件时，必须将边界算子分解为切向和法向部分，分别分析其在半空间上的符号计算。

2. 处理非线性项 $(\nabla v) u + (\nabla u)^\top v$：

   • 难点：这是非守恒形式的对流项，标准的能量估计会产生 $\int u \cdot \nabla u \cdot u$ 这样的项，这是零（如果 $u$ 足够光滑）。但在弱解正则性提升的过程中，如何控制非线性项的增长率是关键。
   • 技巧：利用乘积估计（Product Estimates）。在高阶能量估计中，将非线性项视为扰动。由于 $v = (1-\alpha^2\Delta)u$，非线性项实际上是 $u$ 和其高阶导数的乘积。论文利用Sobolev嵌入定理（特别是 $H^2 \subset L^\infty$ 在三维中的性质）和迹定理，确保非线性项相对于耗散项（主部）是高阶小量，从而可以用Gronwall不等式吸收。

3. 压力项的重构：

   • 难点：在变分形式中，压力 $p$ 是通过消除法向力得到的，但在正则性分析中需要证明 $p$ 具有足够的光滑性。
   • 技巧：应用 de Rham 定理 和 正则性提升。首先在分布意义下恢复压力 $\nabla p = f - \Delta^2 u$。然后，利用ADN理论得到的椭圆正则性，直接推断出 $p \in H^3(\Omega)$。这一步证明了压力与速度场的正则性是兼容的。

与之前工作的比较

该论文在Navier-Stokes正则化模型的文献中占据了独特的位置，其主要突破体现在以下几个方面：

1. 首次对Fried-Gurtin壁面涡流模型的完整分析：

   • 背景：Fried和Gurtin在2008-2009年提出了该边界条件作为近壁湍流的连续介质力学模型。虽然物理动机很强，但在该论文之前，缺乏严格的PDE分析。之前的文献主要关注物理推导或数值模拟，没有数学家证明过该系统的全局适定性。
   • 突破：论文声称这是“对该模型的第一个完整解析处理（first complete analytical treatment）”。它不仅证明了弱解的存在，还证明了高阶正则性和全局唯一性，填补了理论上的空白。

2. 改进了边界条件的处理方式：

   • 常规做法：传统的高阶流体模型（如 $u, v, \omega$ 关系的分析）通常假设周期性边界条件或简单的Dirichlet/Neumann条件。这使得傅里叶分析或简单的谱方法可行。
   • 本文改进：论文处理的是有界域上的复杂非齐次（虽然齐次化了）边界条件。它展示了如何将复杂的几何边界条件嵌入到ADN理论框架中。这比之前处理简单边界条件的α模型论文（如Holst-Lunasin-Tsogtgerel的抽象框架）更具几何复杂性。特别是对于 $\gamma$ 参数（$-1 \le \gamma \le 1$）的处理，涵盖了多种物理情形（如 $\gamma=1$ 时的对称梯度形式），这要求对双线性形式进行细致的Case-by-Case分析。

3. Gårding不等式的精细构建：

   • 技术贡献：对于NS-αβ系统，直接的能量估计失效，因为边界项涉及 $\partial_n u$。论文证明了（Theorem 2.6）即使存在这种复杂的边界项，仍然可以建立控制 $H^2$ 范数的Gårding不等式。
   • 比较：相比于只在体内有耗散的算子，该工作必须处理边界上的涡量牵引力带来的正负贡献。特别是当 $\gamma=1$ 时，作者利用第二Korn不等式来处理对称梯度，这比简单的Poincaré不等式更具技巧性，展示了该方法在处理流体边界摩擦模型时的鲁棒性。

4. 从线性正则性到非线性全局适定性的跨越：

   • 背景：许多正则化模型的线性部分是好理解的，但非线性项往往会导致解在有限时间内爆破（blow-up），或者需要非常小的初始数据。
   • 本文结论：该论文证明了对于任意初始数据（只要在适当的空间），解都是全局存在的且保持正则性。这得益于双线性形式的良好性质和能量估计层级的建立。相比于经典的Navier-Stokes方程（3D情况下的全局正则性仍是未解难题），NS-αβ模型作为正则化模型，确实表现出了更好的全局性质，而本文为这种性质提供了坚实的数学依据。

总结来说，这篇论文通过将复杂的物理边界条件（Fried-Gurtin）转化为严格的数学结构，并综合利用变分法、算子半群和ADN理论，解决了该特定流体模型的全局正则性问题，是连接连续介质力学模型与现代偏微分方程理论的桥梁式工作。

了解领域论文（摘要翻译）

1. 带有一般边界条件的半空间上时周期 Cahn-Hilliard-Gurtin 系统作为混合阶系统的研究

The Time-Periodic Cahn-Hilliard-Gurtin System on the Half Space as a Mixed-Order System with General Boundary Conditions

作者: Guillaume Neuttiens, Jonas Sauer

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23582v1

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

摘要翻译

中文标题: 带有一般边界条件的半空间上时周期 Cahn-Hilliard-Gurtin 系统作为混合阶系统的研究

摘要翻译: 本文证明了半空间上时周期 Cahn-Hilliard-Gurtin 系统的适定性与最大正则性结果。为此，我们引入了一类新颖的补余边界条件，将椭圆与抛物型理论中的经典 Lopatinskiĭ-Shapiro 条件推广到了带有一般边界条件的时周期混合阶系统。此外，我们证明了经典 Lopatinskiĭ-Shapiro 条件在一般情况下对于混合阶系统的适定性是不充分的。

2. 不可压缩粘弹性流体的扩散波现象与最优时间衰减

Diffusion wave phenomena and optimal time decay for incompressible viscoelastic flows

作者: Shenghan Li, Yong Wang

类别: math.AP

发布日期: 2025-12-28

arXiv ID: 2512.22921v1

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

摘要翻译

中文标题: 不可压缩粘弹性流体的扩散波现象与最优时间衰减

摘要翻译: 受 D. Hoff 与 K. Zumbrun (Indiana Univ. Math. J. 44: 603-676, 1995) 工作的启发，我们研究了三维不可压缩粘弹性流体中的扩散波现象。通过利用波动方程的表示公式以及球面 $\mathbb{S}^{d-1}$ 上的驻相法，我们建立了该解在整个范围 $1\leq p \leq \infty$ 上的 $L^p$ 衰减估计，这揭示了不可压缩粘弹性流体的双曲性质。

3. 当壳体变为薄板时非线性浅壳模型的渐近行为

Asymptotic behavior of a nonlinear shallow shell model when the shell becomes a plate

作者: Trung Hieu Giang, Ngoc Quynh Nguyen

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22677v1

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

摘要翻译

中文标题: 当壳体变为薄板时非线性浅壳模型的渐近行为

摘要翻译: 本文研究了由 Donnell、Vlasov、Mushtari、Galimov 和 Koiter 提出的一个非线性浅壳模型。更具体地说，我们探讨了关于极小化解渐近行为的问题。我们的结果可以应用于一般的外加力。因此，它极大地推广了 \cite{oana2} 中的结果，后者假设外加力的切向分量为零。

4. 具有准周期势的一维离散Klein-Gordon方程的色散估计

Dispersive estimates for discrete Klein-Gordon equations on one-dimensional lattice with quasi-periodic potentials

作者: Zhiqiang Wan, Heng Zhang

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22613v1

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

摘要翻译

好的，请看以下翻译：

中文标题: 具有准周期势的一维离散Klein-Gordon方程的色散估计

摘要翻译: 我们证明了在$\mathbb Z$上、具有小的实解析准周期势的离散Klein-Gordon方程的$\ell^{1}\!\to\!\ell^{\infty}$色散估计，并表明其时间衰减率保持为$(\tfrac13)^{-}$。作为应用，我们推导出相应的Strichartz估计，并为相关的非线性离散Klein-Gordon方程建立了小数据全局适定性。

5. 具有双线性控制的非线性复 Ginzburg-Landau 方程的小时间近似能控性

Small-time approximate controllability for the nonlinear complex Ginzburg-Landau equation with bilinear control

作者: Xingwu Zeng, Can Zhang

类别: math.OC, math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22512v1

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

摘要翻译

中文标题: 具有双线性控制的非线性复 Ginzburg-Landau 方程的小时间近似能控性

摘要翻译: 本文研究了环面上具有任意整数次幂型非线性的复 Ginzburg-Landau (CGL) 方程的双线性近似能控性。在控制算子的饱和性假设下，我们证明了 CGL 方程的小时间全局能控性。该证明是通过发展 Agrachev 和 Sarychev 在 \cite{AS05,AS06} 中引入的几何控制方法的一种乘法版本而获得的。

6. 空间周期性水波的长时间动力学

Long time dynamics of space periodic water waves

作者: Massimiliano Berti

类别: math.AP

发布日期: 2025-12-26

arXiv ID: 2512.22115v1

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

摘要翻译

好的，请看以下翻译：

中文标题: 空间周期性水波的长时间动力学

摘要翻译: 本文回顾了关于空间周期性水波长时间动力学研究的最新进展，重点关注以下三个方面：1）定常波和行波的准周期解的分岔；2）长时间适定性结果；3）斯托克斯波的调制不稳定性。这些结果的得出，依赖于针对哈密顿半线性偏微分方程的KAM理论与Birkhoff正规形理论的非传统方法，以及针对可逆算子和哈密顿算子分离谱的辛Kato扰动理论。

不相关论文（基本信息）

1. 完全非线性外几何流的旋转对称平移孤子：分类与应用

Rotationally symmetric translating solitons of fully nonlinear extrinsic geometric flows: Classification and Applications

作者: José Torres Santaella

类别: math.DG, math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23623v1

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

摘要: We study rotationally symmetric translators for fully nonlinear extrinsic geometric flows driven by a curvature function, and we establish the fine asymptotics of bowl-type evolutions and, when admissible, the construction and classification of catenoidal-type solutions, together with their asymptotic behavior. Under natural structural and convexity assumptions, we also prove rigidity and uniqueness results within appropriate classes of graphical translators of such curvature flows.

2. 具有退化系数的混合阶非线性分数阶方程的正则性

Regularity for mixed-order nonlinear fractional equations with degenerate coefficients

作者: Ho-Sik Lee, Jihoon Ok, Kyeong Song

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23359v1

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

摘要: We consider a class of nonlinear integro-differential equations whose leading operator is obtained as a superposition of $(-Δ_{p})^{s}$ and $(-Δ_{p})^{t}$, where $0<s<t<1<p<\infty$, weighted via two possibly degenerate coefficients $a(\cdot,\cdot),b(\cdot,\cdot) \ge 0$. We prove local boundedness and Hölder regularity of its weak solutions under natural assumptions on the coefficients $a(\cdot,\cdot)$, $b(\cdot,\cdot)$ and the powers $s,t$, and $p$. Moreover, when $a(\cdot,\cdot) \equiv 1$, we also prove a Harnack inequality for weak solutions.

3. 一类双线性四阶抛物型方程的小时间全局可控性

Small-time global controllability of a class of bilinear fourth-order parabolic equations

作者: Subrata Majumdar, Debanjit Mondal

类别: math.OC, math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23339v1

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

摘要: In this work, we investigate the small-time global controllability properties of a class of fourth-order nonlinear parabolic equations driven by a bilinear control posed on the one-dimensional torus. The controls depend only on time and act through a prescribed family of spatial profiles. Our first result establishes the small-time global approximate controllability of the system using three scalar controls, between states that share the same sign. This property is obtained by adapting the geometric control approach to the fourth-order setting, using a finite family of frequency-localized controls. We then study the small-time global exact controllability to non-zero constant states for the concerned system. This second result is achieved by analyzing the null controllability of an appropriate linearized fourth-order system and by deducing the controllability of the nonlinear bilinear model through a fixed-point argument together with the small-time global approximate control property.

4. 奇异边界附近的奇异 Yamabe 问题的解

Solutions of the singular Yamabe problem near singular boundaries

作者: Weiming Shen, Zhehui Wang, Jiongduo Xie

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23331v1

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

摘要: In this paper, we investigate the asymptotic behaviors of solutions to the singular Yamabe problem with negative constant scalar curvature near singular boundaries and derive optimal estimates, where the background metrics are not assumed to be conformally flat. Specifically, we demonstrate that for a wide class of Lipschitz domains with asymptotic conical structure, the local positive solutions are well approximated by the positive solutions in the tangent cones at singular boundary points. This extends the results of [10, 12, 26].

5. 度量图上非线性磁Schrödinger方程的归一化解

Normalized solutions of nonlinear magnetic Schrödinger equations on metric graphs

作者: Pietro d’Avenia, Zhentao He, Chao Ji

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23321v1

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

摘要: In this paper we first establish the theory of a magnetic Sobolev space $H^1_A(\mathcal{G},\mathbb{C})$ on metric graphs $\mathcal{G}$ and we prove the self-adjointness of its corresponding magnetic Schrödinger operator. Then, in this setting, we investigate the existence and multiplicity of normalized solutions to nonlinear magnetic Schrödinger equations on compact metric graphs and on noncompact metric graphs with localized nonlinearities or nonlinearities acting on whole metric graphs, covering the mass-subcritical, mass-critical, and mass-supercritical cases.

6. 开书上非线性 Schrödinger 方程的基态及到度量图的维度约化

Ground States for the Nonlinear Schr{ö}dinger Equation on Open Books and Dimensional Reduction to Metric Graphs

作者: Stefan Le Coz, Boris Shakarov

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23286v1

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

摘要: In this work, we study the dimensional reduction of stationary states in the shrinking limit for a broad class of two-dimensional domains, called open books, to their counterparts on metric graphs. An open book is a two-dimensional structure formed by rectangular domains sharing common boundaries. We first develop a functional-analytic framework suited to variational problems on open books and establish the existence of solutions as constrained action minimizers. For graph-based open books (i.e., those isomorphic to the product of a graph with an interval) we prove the existence of a sharp transition in the dimensionality of ground states. Specifically, there exists a critical transverse width: below this threshold, all ground states coincide with the ground states on the underlying graph trivially extended in the transverse direction; above it, ground states become genuinely two-dimensional.

7. 具有Sobolev超临界非线性的Hénon抛物方程的爆破速率

On blow-up rate for the Hénon parabolic equation with Sobolev supercritical nonlinearity

作者: Kotaro Hisa, Yukihiro Seki

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23271v1

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

摘要: We discuss the Hénon parabolic equation $\partial_t u = Δu + |x|^σu^p$ in a finite ball in $\mathbb{R}^N$ under the Dirichlet boundary condition, where $N\ge1$, $p>1$, and $σ>0$. We assume that the exponent $p$ is supercritical in the Sobolev sense. Since the spatial potential term $|x|^σ$ vanishes at the origin, solutions seem less likely to blow up at the origin. We construct a solution that blows up at the origin and also carry out an analysis of blow-up rate of solutions. In particular, if $p$ is less than the Joseph–Lundgren exponent, all blow-ups are shown to be of Type I. The lower bound corresponding to Type I rate is also shown for some particular blow-up solutions. As by products, we present a basic result on classification to threshold solutions for every $p>1+σ/N$.

8. 具有非光滑非线性的非线性标量场方程的无穷多正解

Infinitely many positive solutions to nonlinear scalar field equation with nonsmooth nonlinearity

作者: Tianhao Liu, Juncheng Wei, Wenming Zou

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23207v1

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

摘要: This paper investigates the existence of infinitely many positive solutions for the logarithmic scalar field equation \begin{equation} \tag{$P$} \label{equ1} -Δu+ V(x) u= u\log u^2, \quad u\in H^1(\mathbb{R}^N), \end{equation} and its counterpart with prescribed $L^2$-norms \begin{align}\label{equ2} \tag{$P_N$} & -Δu+ V(x) u +λu= u\log u^2, \quad u\in H^1(\mathbb{R}^N), &\int_{\mathbb{R}^N} u^2 ~\mathrm{d}x=a^2>0, \end{align} which come from physically relevant situations. Here, $N\geq 2$, $V:\mathbb{R}^N\to \mathbb{R}$ is a non-symmetric and non-periodic potential satisfying certain decay conditions, $ a $ is prescribed constant, and $λ$ arises as an unknown Lagrange multipliers. For problem \eqref{equ1}, using purely variational methods, we establish the existence of multi-bump positive solutions with either finitely or infinitely many bumps. For normalized problem \eqref{equ2}, we prove the existence of normalized multi-bump positive solutions with a finite number of bumps. The main difficulty comes from the nonsmooth nature of logarithmic nonlinearity, which introduces some challenges to the variational framework. In particular, the corresponding energy functional is not of class $C^1$ on $H^1(\mathbb{R}^N)$, which prevents the direct application of standard critical point theory for $C^1$ functional or any reduction methods for $C^{1+σ}$ nonlinearity. The main ingredients in this paper are nonsmooth critical point theory, localized variational methods and a max-min argument. To the best of our knowledge, this paper appears to be the first successful application of the localized variational method to nonsmooth functionals.

9. 带自然边界条件的 $\mathbb{S}^2$ 上薛定谔流的局部适定性

Local well-posedness of the Schrödinger flow into $\mathbb{S}^2$ with natural boundary conditions

作者: Bo Chen, Youde Wang

类别: math.AP, math.DG

发布日期: 2025-12-29

arXiv ID: 2512.23201v1

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

摘要: In this paper, we develop a new approximation scheme to resolve the local well-posedness problem for the Landau-Lifshitz equation (i.e., the Schrödinger flow into the standard unit 2-sphere $\mathbb{S}^2\subset \mathbb{R}^3$) with natural boundary conditions.

10. 基尔霍夫-鲁斯函数临界点的定性分析

Qualitative analysis on the critical points of the Kirchhoff-Routh function

作者: Francesca Gladiali, Massimo Grossi, Peng Luo, Shusen Yan

类别: math.AP

发布日期: 2025-12-29

arXiv ID: 2512.23172v1

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

摘要: In this paper, we study the number of critical points of the Kirchhoff-Routh function \begin{equation*} \mathcal{KR}_D(x,y)=Λ_1^2\mathcal{R}_D(x)+Λ_2^2\mathcal{R}_D(y)-2Λ_1Λ_2G_D(x,y), \end{equation*} where $D$ is a bounded domain in $\mathbb{R}^2$, $x,y\in D$, $Λ_1,Λ_2>0$, $\mathcal{R}_D$ is the Robin function, and $G_D$ is the Green function of the operator $-Δ$ with $0$ Dirichlet boundary condition on $D$. This function arises from concentration phenomena in nonlinear elliptic problems and from the de-singularization problem for the steady Euler equation. For domains with a small hole, we establish not only the exact number and the location of the critical points of $\mathcal{KR}_D$, but also their nondegeneracy. We show that the location of the hole plays a crucial role. Finally in the context of elliptic problems, we establish the existence of multiple two-peak solutions.

11. 磁流与温控流的吉亚马鲁正规算子

Guillarmou’s Normal Operator for Magnetic and Thermostat Flows

作者: Sebastián Muñoz-Thon, Sean Richardson

类别: math.AP, math.DG, math.DS

发布日期: 2025-12-28

arXiv ID: 2512.23106v1

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

摘要: Guillarmou’s normal operator over a closed Anosov manifold is analogous to the classical normal operator of the geodesic X-ray transform over manifolds with boundary. In this paper, we generalize this normal operator, under some dynamical assumptions, to thermostat flows as well as to the case of the magnetic flows. In particular, we show that these generalized normal operators are elliptic pseudodifferential operators of order -1 in each case. As an application, we prove a stability estimate for the magnetic X-ray transform.

12. 基于均匀化方法的电穿孔模型非线性时变宏观电导率推导

Derivation of nonlinear time-dependent macroscopic conductivity for an electropermeabilization model via homogenization

作者: Tobias Gebäck, Ioanna Motschan-Armen, Irina Pettersson

类别: math.AP

发布日期: 2025-12-28

arXiv ID: 2512.23007v1

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

摘要: We study a phenomenological electropermeabilization model in a periodic medium representing biological tissue. Starting from a cell-level model describing the electric potential and the degree of porosity, we perform dimension analysis to identify a relevant scaling in terms of a small parameter $\ve$ - the ratio between the cell and the tissue size. The electric potential satisfies electrostatic equations in the extra- and intracellular domains, while its jump across the cell membrane evolves according to a nonlinear law coupled with an ordinary differential equation for the porosity degree. We prove the well-posedness of the microscopic problem, derive a priori estimates, obtain formal asymptotics, and rigorously justify the expansion combining two-scale convergence with monotonicity arguments. The resulting macroscopic model exhibits memory effects and a nonlinear, time-dependent effective current. It captures the nontrivial evolution of effective conductivity, including a characteristic drop reflecting the capacitive behavior of the lipid bilayer, in agreement with experimental data. Numerical computations of the effective conductivity confirm that, although microscopic conductivity is constant, tissue conductivity varies nonlinearly with electric field strength, showing a sigmoid trend. This suggests a rigorous mathematical explanation for experimentally observed conductivity dynamics.

13. Anderson-Gross-Pitaevskii方程的驻波

Standing waves of the Anderson-Gross-Pitaevskii equation

作者: Samaël Mackowiak

类别: math.AP, math-ph, math.PR

发布日期: 2025-12-28

arXiv ID: 2512.22960v1

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

摘要: In this paper, we study standing waves for the Anderson-Gross-Pitaevskii equation in dimension 1 and 2. The Anderson-Gross-Pitaevskii equation is a nonlinear Schrödinger equation with a confining potential and a multiplicative spatial white noise. Standing waves are characterized by a profile which is invariant by the dynamic and solves a nonlinear elliptic equation with spatial white noise potential. We construct such solutions via variational methods and obtain some results on their regularity, localization and stability.

14. 跨扩散-捕食者-趋化性模型中栖息地异常的确定

Determining habitat anomalies in cross-diffusion predator-prey chemotaxis models

作者: Yuhan Li, Hongyu Liu, Catharine W. K. Lo

类别: math.AP, q-bio.CB

发布日期: 2025-12-28

arXiv ID: 2512.22946v1

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

摘要: This paper addresses an open inverse problem at the interface of mathematical analysis and spatial ecology: the unique identification of unknown spatial anomalies – interpreted as zones of habitat degradation – and their associated ecological parameters in multi-species predator-prey systems with multiple chemical signals, using only boundary measurements. We formulate the problem as the simultaneous recovery of an unknown interior subdomain and discontinuous ecological interaction rules across its boundary. A unified theooretical framework is developed that unique determines both the anomaly’s geometry and discontinuous coefficients characterizing the altered interactions within the degraded region. Our results cover smooth anomalies in time-dependent systems and are extended to non-smooth polyhedral inclusions in stationary regimes. This work bridges a gap between ecological sensing and the quantitative inference of internal habitat heterogeneity, offering a mathamtical basis for detecting and characterizing habitat degradation from limited external data.

15. 大平流主特征值问题的精细极限剖面

Refined Limiting Profiles of the Principal Eigenvalue Problems with Large Advection

作者: Yujin Guo, Yuan Lou, Hongfei Zhang

类别: math.AP

发布日期: 2025-12-28

arXiv ID: 2512.22918v1

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

摘要: In this paper, we are concerned with the following eigenvalue problem with an advection term: \begin{equation}\label{0.1} \left{ \begin{split} -εΔφ-2α\nabla m(x)\cdot\nabla φ+V(x)φ&=λφ \ \text{in}\ \ Ω,\ φ&=0\ \ \hbox{on}\ \ \partialΩ, ~~~\text{(0.1)} \end{split} \right. \end{equation} where $Ω\subset\mathbb{R}^N~(N\geq1)$ satisfying $\partialΩ\in C^{2}$ is a bounded domain and contains the origin as an interior point, the constants $ε>0$ and $α>0$ are the diffusive and advection coefficients, respectively, and $m(x)\in C^{2}(\barΩ)$, $V (x)\in C^γ(\barΩ)~(0<γ<1)$ are given functions. We analyze the refined limiting profiles of the principal eigenpair $(λ, φ)$ for (0.1) as $α\rightarrow\infty$, which display the visible effect of the large advection on $(λ, φ)$. It expects that our argument is applicable to investigating the refined expansions of the general principal eigenvalue problems.

16. 空间形式中Steklov特征值的缩放不等式及扭曲积流形上的尖锐特征值估计

Scaling inequalities for Steklov eigenvalues in space forms and sharp eigenvalue estimates on warped product manifolds

作者: Zongyi Lv, Changwei Xiong, Yuxun Zou

类别: math.DG, math.AP, math.SP

发布日期: 2025-12-28

arXiv ID: 2512.22885v1

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

摘要: In the first part, we derive monotonicity of the normalized spectra for the second-order Steklov problem and two fourth-order Steklov problems on the $2$-dimensional geodesic disks with respect to the geodesic radius in the sphere and the hyperbolic space. The normalizations are made using four natural geometric factors. As corollaries, we get Escobar-type bounds for Steklov eigenvalues on $2$-dimensional geodesic disks with varying curvature in space forms. We also get two monotonicity results for higher-dimensional cases. In the second part, we obtain some sharp bounds concerning the spectra of the two fourth-order Steklov problems on warped product manifolds with non-negative Ricci curvature and a strictly convex boundary. In particular, we confirm Qiaoling Wang and Changyu Xia’s conjecture (2018) on the sharp lower bound of the first non-zero eigenvalue of a fourth-order Steklov problem in the case of $3$-dimensional warped product manifolds.

17. Blume-Emery-Griffiths模型中离散界面的结晶运动：部分润湿

Crystalline Motion of discrete interfaces in the Blume-Emery-Griffiths Model: partial wetting

作者: Marco Cicalese, Giuliana Fusco, Giovanni Savaré

类别: math.AP, cond-mat.soft

发布日期: 2025-12-28

arXiv ID: 2512.22870v1

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

摘要: We continue the variational study of the discrete-to-continuum evolution of lattice systems of Blume-Emery-Griffith type which model two immiscible phases in the presence of a surfactant. In our previous work \cite{CFS}, we analyzed the case of a completely wetted crystal and described how the interplay between surfactant evaporation and mass conservation leads to a transition between crystalline mean curvature flow and pinned evolutions. In the present paper, we extend the analysis to the regime of partial wetting, where the surfactant occupies only a portion of the interface. Within the minimizing-movements scheme, we rigorously derive the continuum evolution and show how partial wetting introduces a complex coupling between interfacial motion and redistribution of surfactant. The resulting evolution exhibits new features absent in the fully wetted case, including the coexistence of moving and pinned facets or the emergence and long-lived metastable states. This provides, to our knowledge, the first discrete-to-continuum variational description of partially wetted crystalline interfaces, bridging the gap between microscopic lattice models and experimentally observed surfactant-induced pinning phenomena in immiscible systems.

18. 具有排斥势的旋转非线性薛定谔方程的爆破率

Blowup rate for rotational NLS with a repulsive potential

作者: Yi Hu, Yongki Lee, Shijun Zheng

类别: math.AP, math-ph

发布日期: 2025-12-28

arXiv ID: 2512.22821v1

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

摘要: In this paper we give an analytical proof of the ``$\log$-$\log$’’ blowup rate for mass-critical nonlinear Schrödinger equation (NLS) with a rotation ($Ω\neq 0$) and a repulsive harmonic potential $V_γ(x) = \textrm{sgn}(γ) γ^2 |x|^2$, $γ< 0$ when the initial data has a mass slightly above that of $Q$, the ground state solution to the free NLS. The proof is based on a virial identity and an $\mathcal{R}_γ$-transform, a pseudo-conformal transform in this setting. Further, we obtain a limiting behavior description concerning the mass concentration near blowup time. A remarkable finding is that increasing the value $|γ|$ for the repulsive potential $V_γ$ can give rise to global in time solution for the focusing RNLS, which is in contrast to the case where $γ$ is positive. This kind of phenomenon was earlier observed in the non-rotational case $Ω= 0$ in Carles’ work. In addition, we provide numerical simulations to partially illustrate the blowup profile along with the blowup rate using dynamic rescaling and adaptive mesh refinement method.

19. 针对一般自由数据求解约束方程

Solving the constraint equation for general free data

作者: Xuantao Chen, Sergiu Klainerman

类别: gr-qc, math.AP, math.DG

发布日期: 2025-12-27

arXiv ID: 2512.22704v1

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

摘要: We revisit the problem of solving the Einstein constraint equations in vacuum by a new method, which allows us to prescribe four scalar quantities, representing the full dynamical degrees of freedom of the constraint system. We show that once appropriate gauge conditions have been chosen and four scalars freely specified (modulo $\ell\leq 1$ modes), we can rewrite the constraint equations as a well-posed system of coupled transport and elliptic equations on $2$-spheres, which we solve by an iteration procedure. Our method provides a large class of exterior solutions of the constraint equations that can be matched to given interior solutions, according to the existing gluing techniques. As such, it can be applied to provide a large class of initial Cauchy data sets evolving to black holes, generalizing the well-known result of the formation of trapped surfaces due to Li and Yu. Though in our main theorem, we only specify conditions consistent with $g-g_{Schw}=O(r^{-1-δ})$, $k=O(r^{-2-δ})$, the method is flexible enough to be applied in many other situations. It can, in particular, be easily adapted to construct arbitrarily fast decaying data. We expect, moreover, that our method can also be applied to construct data with slower decay, such as that used by Shen. In fact, an important motivation for developing our method is to show that the result of Shen is sharp, i.e., construct small, smooth initial data sets which violate Shen’s decay conditions, and for which the stability of the Minkowski space result is wrong.

20. 与四阶引力耦合的薛定谔方程的基态——第一部分：$K_{a, b} \leq 0$ 情形

Ground states of the Schrödinger equation coupled with fourth-order gravitation – Part 1: the case $K_{a, b} \leq 0$

作者: Gustavo de Paula Ramos

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22619v1

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

摘要: We are interested in the existence and asymptotic behavior of ground states of the following normalized nonlocal semilinear problem:

where

$0 \leq a, b \leq \infty$; $V$ denotes a singular potential that vanishes at infinity and the unknowns are $ω\in \mathbb{R}$, $u \colon \mathbb{R}^3 \to \mathbb{R}$. This problem is obtained by looking for standing waves of the Schrödinger equation coupled with the nonrelativistic gravitational potential prescribed by a family of fourth-order gravity theories. In this paper, (i) we obtain a complete picture of the existence/nonexistence of ground states of the associated autonomous problem for every possible geometry of $K_{a, b}$, (ii) we obtain conditions that ensure the existence of ground states of the nonautonomous problem when $K_{a, b} \leq 0$ and (iii) we prove that as

ground states of this problem respectively converge to a ground state of (1) the Schrödinger equation, (2) the Choquard equation and (3) a rescaling of the Choquard equation.

21. 周期性磁性赝微分算子的布洛赫-弗洛凯分解中的纤维算子

The fibre operators in the Bloch-Floquet decomposition of periodic magnetic pseudo-differential operators

作者: Horia D. Cornean, Bernard Helffer, Radu Purice

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22547v1

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

摘要: We study the structure of the fibre operators corresponding to periodic magnetic pseudo-differential operators having periodic magnetic potentials. We obtain explicit formulas for their distribution kernel, both when these fibres are seen as operators on the $d$-dimensional torus, and also when they are seen as infinite matrices acting on a discrete $\ell^2$ space via a discrete Fourier transform. Moreover, using these distribution kernels we prove that the fibre operators are toroidal pseudo-differential operators.

22. 由p-Laplacian算子驱动的反应-扩散-对流方程中的有限传播与饱和现象

Finite propagation and saturation in reaction-diffusion-advection equations governed by p-Laplacian operator

作者: Cristina Marcelli

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22493v1

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

摘要: The paper concerns front propagation for the following mono-stable reaction-diffusion-advection equation

Besides existence and non-existence results for traveling wave solutions, the main focus is their classification: we provide criteria to establish if they attain one or both the equilibria at a finite time and in this case, if they are continuable as $C^1$-solutions or if they are sharp solutions.

23. 拓扑绝缘体波导的逆散射研究

Inverse scattering for waveguides in topological insulators

作者: Guillaume Bal, Xixian Wang, Zhongjian Wang

类别: math-ph, math.AP, math.NA

发布日期: 2025-12-27

arXiv ID: 2512.22480v1

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

摘要: This paper concerns the inverse scattering problem of a topologically non-trivial waveguide separating two-dimensional topological insulators. We consider the specific model of a Dirac system. We show that a short-range perturbation can be fully reconstructed from scattering data in a linearized setting and in a finite-dimensional setting under a smallness constraint. We also provide a stability result in appropriate topologies. We then solve the problem numerically by means of a standard adjoint method and illustrate our theoretical findings with several numerical simulations.

24. 海森堡群上通过移动球法研究的CR Yamabe方程

CR Yamabe Equation on the Heisenberg Group via the method of moving spheres

作者: Congwen Liu

类别: math.AP

发布日期: 2025-12-27

arXiv ID: 2512.22458v1

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

摘要: In this paper, we classify positive solutions to the CR Yamabe equation on the Heisenberg group $\mathbb{H}^n$. We show that all such solutions are Jerison-Lee bubbles, without imposing any finite-energy or a priori symmetry assumptions. This result can be regarded as an analogue for $\mathbb{H}^n$ of the celebrated Caffarelli-Gidas-Spruck classification theorem in $\mathbb{R}^n$. To establish this, we develop a systematic approach to implement the method of moving spheres in the setting of the Heisenberg group.

25. 时间重参数化，而非分数阶微积分：对可导数导数的重新评估

Time Reparametrization, Not Fractional Calculus: A Reassessment of the Conformable Derivative

作者: Aziz El Ghazouani, Fouad Ibrahim Abdou Amir, Khoulane Mohamed, M’hamed Elomari

类别: math.AP, math-ph, math.DS, math.NA

发布日期: 2025-12-26

arXiv ID: 2512.22366v1

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

摘要: The conformable derivative has been promoted in numerous publications as a new fractional derivative operator. This article provides a critical reassessment of this claim. We demonstrate that the conformable derivative is not a fractional operator but a useful computational tool for systems with power-law time scaling, equivalent to classical differentiation under a nonlinear time reparametrization. Several results presented in the literature as novel fractional contributions can be reinterpreted within a classical framework. We show that problems formulated using the conformable derivative can be transformed into classical formulations via a change of variable. The solution is derived classically and then transformed back, this reformulation highlights the absence of genuinely nonlocal fractional effects. We provide a theoretical analysis, numerical simulations comparing conformable, classical, and truly fractional (Caputo) models, and discuss the reasons why this misconception persists. Our results suggest that classical derivatives, as well as established fractional derivatives, offer a more faithful framework for modeling memory-dependent phenomena.