生成时间: 2025年10月27日 14:53:53

重点关注论文数量: 3

了解领域论文数量: 3

不相关论文数量: 22

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

1. Global Weak Solutions of a Thermodynamically Consistent Diffuse Interface Model for Nonhomogeneous Incompressible Two-phase Flows with a Soluble Surfactant

一类具有可溶性表面活性剂的非均质不可压缩两相流热力学一致扩散界面模型的整体弱解

作者: Bohan Ouyang, Maurizio Grasselli, Hao Wu

类别: math.AP, 35K52, 35Q35, 76D45, 76T06

发布日期: 2025-05-24

arXiv ID: 2505.18662v2

链接: http://arxiv.org/abs/2505.18662v2

详细分析

中文标题: 一类具有可溶性表面活性剂的非均质不可压缩两相流热力学一致扩散界面模型的整体弱解

1. 研究对象和背景

本文研究的是一类描述含有可溶性表面活性剂（soluble surfactant）的非均质不可压缩两相流（nonhomogeneous incompressible two-phase flows）的动力学行为的热力学一致扩散界面模型。该模型适用于两个黏性、不可压缩、牛顿流体在有界区域 $\Omega \subset \mathbb{R}^d$（$d=2$ 或 $3$）中的运动，且两相流体密度不匹配（unmatched densities），并包含一种可在两相中溶解并在界面处吸附的表面活性剂。

这类系统广泛存在于自然界与工业应用中，如石油开采、食品加工、微流控技术等。表面活性剂因其两亲性分子结构（hydrophilic head 和 hydrophobic tail）倾向于富集于油水等两相流体之间的自由界面，从而显著降低界面张力（interfacial tension），促进乳化或泡沫稳定；同时，其浓度的空间不均匀性还会引发马兰戈尼力（Marangoni forces），影响界面稳定性与流动动力学。

为数学建模此类复杂多尺度问题，作者采用扩散界面方法（diffuse interface model），将原本尖锐的界面视为一个有限厚度的过渡层，并引入相场变量（phase-field variables）来区分不同流体相及表面活性剂浓度。相比经典锐界面模型（sharp interface model）需要显式追踪移动边界，扩散界面法能自然处理拓扑变化（如液滴破裂、融合等），更适用于数值模拟与理论分析。

所研究的完整偏微分方程组如下，在 $\Omega \times (0, \infty)$ 上定义：

$$ \partial_t (\rho(\phi)\mathbf{u}) + \mathrm{div}(\mathbf{u} \otimes (\rho(\phi)\mathbf{u} + \mathbf{J})) - \mathrm{div}(\nu(\phi) D\mathbf{u}) + \nabla \pi = 0, $$
$$ \mathrm{div}\, \mathbf{u} = 0, $$
$$ \partial_t \phi + \mathbf{u} \cdot \nabla \phi + \sigma_1(\phi)(\phi - c) = \mathrm{div}(m_\phi(\phi) \nabla \mu_\phi), $$
$$ \mu_\phi = -\Delta \phi + F'_\phi(\phi) + \sigma_2 N(\phi - \overline{\phi}) + \partial_\phi G(\phi, \psi), $$
$$ \partial_t \psi + \mathbf{u} \cdot \nabla \psi = \mathrm{div}(m_\psi(\psi) \nabla \mu_\psi), $$
$$ \mu_\psi = -\beta \Delta \psi + F'_\psi(\psi) + \partial_\psi G(\phi, \psi). $$

其中：

  • $\mathbf{u}: \Omega \times [0,\infty) \to \mathbb{R}^d$ 是体积平均速度场（满足 $\mathrm{div}\, \mathbf{u} = 0$）；
  • $\pi: \Omega \times [0,\infty) \to \mathbb{R}$ 是压力；
  • $\phi: \Omega \times [0,\infty) \to [-1,1]$ 表示两种流体体积分数之差（相场变量）；
  • $\psi: \Omega \times [0,\infty) \to [0,1]$ 表示表面活性剂的体积分数；
  • $\rho(\phi) = \frac{\rho_1 - \rho_2}{2} \phi + \frac{\rho_1 + \rho_2}{2}$ 为平均密度（依赖于 $\phi$）；
  • $\nu(\phi)$ 为平均黏度函数；
  • $D\mathbf{u} = \frac{1}{2}(\nabla \mathbf{u} + \nabla \mathbf{u}^T)$ 为对称梯度；
  • $\mathbf{J} = -\frac{\rho_1 - \rho_2}{2} m_\phi(\phi) \nabla \mu_\phi$ 为由于密度不匹配引起的相对通量，是保证热力学一致性的重要项；
  • $m_\phi(\phi), m_\psi(\psi)$ 分别为二元流体与表面活性剂的迁移率函数（mobility coefficients），允许退化（degenerate）情形；
  • $F_\phi, F_\psi$ 为自由能密度函数，取为物理相关的奇异势（singular potentials），如 Flory–Huggins 型：
    $$ F_\phi(s) = \frac{\theta_{c1}}{2}[(1+s)\ln(1+s)+(1-s)\ln(1-s)] + \frac{\hat{\theta}_1}{2}(1 - s^2), \quad s \in (-1,1) $$
    $$ F_\psi(s) = \frac{\theta_{c2}}{2}[s\ln s+(1-s)\ln(1-s)] + \frac{\hat{\theta}_2}{2}s(1-s), \quad s \in (0,1) $$
    这些奇异势确保了解 $\phi(t,x) \in [-1,1]$, $\psi(t,x) \in [0,1]$ 几乎处处成立，符合物理意义。
  • $G(\phi,\psi)$ 为流体混合物与表面活性剂之间的相互作用能量密度，典型形式为：
    $$ G(\phi,\psi) = \frac{\gamma_1}{2} \psi \phi^2 - \frac{\gamma_2}{4} \psi (1 - \phi^2)^2, $$
    第一项表示体相中游离表面活性剂的能量代价，第二项促使表面活性剂聚集于界面（$\phi \approx 0$ 区域）。
  • $N$ 为带齐次诺伊曼边界的负拉普拉斯算子的逆，即 $-\Delta N f = f$ 且 $\partial_n N f = 0$，用于刻画非局部相互作用（Ohta–Kawasaki 型泛函）；
  • $\sigma_1(\phi)$ 为化学反应系数，对应 Oono 类型的质量源项，允许组分间发生可逆化学反应；
  • $c \in (-1,1)$ 为目标平衡值。

系统配备如下边界条件：

$$ \mathbf{u} = 0, \quad \partial_n \phi = m_\phi(\phi) \partial_n \mu_\phi = \partial_n \psi = m_\psi(\psi) \partial_n \mu_\psi = 0, \quad \text{on } \partial\Omega \times (0,\infty), $$

即无滑移边界条件（no-slip）用于速度场，齐次诺伊曼边界条件用于所有相场变量及其化学势，体现封闭系统。

初始条件为：

$$ \mathbf{u}|_{t=0} = \mathbf{u}_0(x), \quad \phi|_{t=0} = \phi_0(x), \quad \psi|_{t=0} = \psi_0(x), \quad \text{in } \Omega. $$

文中特别强调了系统的热力学一致性：通过第一定律、Onsager 变分原理推导得出，系统满足质量守恒与能量耗散律。总能量定义为：

$$ E_{\mathrm{tot}}(\mathbf{u}, \phi, \psi) = \int_\Omega \frac{\rho(\phi)}{2} |\mathbf{u}|^2 dx + E_{\mathrm{free}}(\phi, \psi), $$

其中自由能 $E_{\mathrm{free}}$ 由下式给出：

$$ E_{\mathrm{free}}(\phi, \psi) = \int_\Omega \left[ \frac{1}{2}|\nabla \phi|^2 + F_\phi(\phi) + \frac{\sigma_2}{2} |\nabla N(\phi - \overline{\phi})|^2 + \frac{\beta}{2} |\nabla \psi|^2 + F_\psi(\psi) + G(\phi, \psi) \right] dx. $$

对于光滑解，有以下能量恒等式：

$$ \frac{d}{dt} E_{\mathrm{tot}} + \int_\Omega \left( \nu(\phi)|D\mathbf{u}|^2 + m_\phi(\phi)|\nabla \mu_\phi|^2 + m_\psi(\psi)|\nabla \mu_\psi|^2 \right) dx + \int_\Omega (\phi - c)\sigma_1(\phi)\left( \mu_\phi - \frac{\rho_1 - \rho_2}{4} |\mathbf{u}|^2 \right) dx = 0. $$

这表明系统能量随时间单调递减（若最后一项非正），体现了耗散机制。

函数空间方面，使用标准 Sobolev 空间 $L^p(\Omega), H^k(\Omega)$ 等，并引入向量空间如 $L^p_\sigma(\Omega)$ 表示零散度的 $L^p$ 向量场闭包。关键算子包括 Neumann 拉普拉斯 $A_N$ 及其逆 $N$，以及 Helmholtz 投影 $P_\sigma$，用于处理不可压 Navier–Stokes 部分。

综上所述，该模型是一个高度耦合的非线性偏微分方程组，涉及非均质 Navier–Stokes 方程与耦合 Cahn–Hilliard 系统，具有变密度、变黏度、退化迁移率、奇异势、非局部项、质量源项等多个数学挑战点，研究其整体弱解的存在性具有重要理论价值与实际意义。

2. 主要定理或主要结果

本文证明了两类情形下的整体弱解存在性：

定理 2.1（非退化迁移率情形的整体弱解）

假设条件：

  • (H0)–(H2) 成立；
  • 迁移率函数 $m_\phi, m_\psi$ 非退化，即存在常数 $m_*, m^* > 0$ 使得
    $$ m_* \leq m_\phi(s) \leq m^*, \quad m_* \leq m_\psi(s) \leq m^*, \quad \forall s \in [-1,1] \text{（或 } [0,1]). $$
  • 初始数据满足：
    $$ \mathbf{u}_0 \in \mathbf{H}, \quad \phi_0 \in H^1(\Omega), \quad \psi_0 \in H^1(\Omega), \quad F_\phi(\phi_0) \in L^1(\Omega), \quad F_\psi(\psi_0) \in L^1(\Omega). $$

结论：存在一对 $(\mathbf{u}, \phi, \psi)$ 满足适当的正则性（见原文定义 2.1），使得对任意测试函数，弱形式成立，且满足能量不等式：

$$ E_{\mathrm{tot}}(\mathbf{u}(t), \phi(t), \psi(t)) + \int_0^t \int_\Omega \left( \nu(\phi)|D\mathbf{u}|^2 + m_\phi|\nabla \mu_\phi|^2 + m_\psi|\nabla \mu_\psi|^2 \right) dx d\tau \leq E_{\mathrm{tot}}(\mathbf{u}_0, \phi_0, \psi_0). $$

定理 2.2（退化迁移率情形的整体弱解）

假设条件：

  • (H0)–(H2) 成立；
  • 迁移率函数退化于纯相态：
    $$ m_\phi(s) = (1 - s^2)^k, \quad k \geq 1; \quad m_\psi(s) = s^k (1 - s)^k, \quad k \geq 1; $$
    并假设 $\sigma_1(s)(s - c)$ 在某种积分意义下可控（见原文 (2.8)）；
  • 初始数据同上。

结论：存在一个整体弱解 $(\mathbf{u}, \phi, \psi)$，满足修正后的弱形式与能量不等式。

这两个定理共同构成了论文的核心成果：首次在一般退化迁移率与双重奇异势条件下，建立了该复杂两相流系统的整体弱解存在性。

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

a) 核心分析学工具

  1. 半隐式时间离散化（Semi-implicit time discretization）
    对于非退化情形（Theorem 2.1），作者采用时间步进法构造近似解序列。每一步求解一个耦合的椭圆系统，其中非线性项（如 $\rho(\phi)$, $\nu(\phi)$, $G(\phi,\psi)$）在前一时刻取值，保持系统线性化。此方法继承自 Abels–Garcke–Grün 的工作，但需调整以适应表面活性剂与质量源项。

  2. Leray–Schauder 不动点定理
    在每一时间步，通过构造一个映射并验证其紧性与先验估计，利用该不动点定理证明离散问题解的存在性。这是处理非线性椭圆系统的经典工具。

  3. 能量估计与熵估计（Energy and entropy estimates）
    能量估计基于前述能量恒等式，控制速度、相场梯度与化学势。特别地，对于退化迁移率情形，仅靠能量估计不足以获得足够紧性，故引入熵估计（Lemma 4.1）：

    $$ \int_\Omega \left( F_\phi(\phi) + F_\psi(\psi) \right) dx + \int_0^t \int_\Omega \left( m_\phi(\phi) |\nabla \mu_\phi|^2 + m_\psi(\psi) |\nabla \mu_\psi|^2 \right) dx d\tau \leq C, $$

    它利用了奇异势在边界趋于无穷的性质，提供对相场变量接近纯相时行为的控制。

  4. 紧性原理与弱收敛
    利用 Aubin–Lions 引理、Banach–Alaoglu 定理等，从近似解序列中提取弱/弱*收敛子列，并结合紧嵌入（如 $H^1 \hookrightarrow L^2$ 紧）证明强收敛性，以处理非线性项的极限。

  5. 逼近策略（Approximation scheme）
    对于退化迁移率情形（Theorem 2.2），作者设计了一套新的逼近方案：

    • 将退化迁移率 $m_\phi^\delta(s) = m_\phi(s) + \delta$ 正则化为非退化；
    • 将奇异势 $F_\phi$ 逼近为光滑势 $F_\phi^\delta = F_\phi^{reg} + \delta F_\phi^{sing}$，保留奇异性但使其可微；
    • 对每个 $\delta > 0$，应用 Theorem 2.1 得到近似解 $(\mathbf{u}^\delta, \phi^\delta, \psi^\delta)$；
    • 利用统一的能量/熵估计，令 $\delta \to 0$ 取极限。

  6. Lebesgue 控制收敛定理
    在极限过程中，某些非线性项（如 $m_\phi(\phi^\delta)\nabla \mu_\phi^\delta$）的收敛性依赖于该定理，需构造合适的控制函数。

  7. Korn 不等式与 Poincaré–Wirtinger 不等式
    用于控制速度场的 $H^1$ 范数，特别是在处理能量恒等式中交叉项 $(\phi - c)\sigma_1(\phi)|\mathbf{u}|^2$ 时，需分时间段估计。

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

  1. 处理质量源项的影响
    质量源项 $(\phi - c)\sigma_1(\phi)$ 可能破坏能量有界性。作者通过假设 $\sigma_1(\phi)(\phi - c)$ 具有适当衰减性（见 (2.8)），并在估计时将时间区间分为 $[0,T]$ 与 $[T,\infty)$ 两部分：前者用 Gronwall 引理，后者利用其衰减性与能量耗散控制。

  2. 处理退化迁移率的奇异性
    传统方法要求迁移率满足特定单调性（如 [12]），本文通过精心设计的逼近格式绕开此限制。关键是让逼近势 $F_\phi^\delta$ 仍保持“结构性”，使得熵估计可一致成立。

  3. 处理双奇异势的耦合效应
    同时存在 $F_\phi$ 与 $F_\psi$ 的奇异性，需确保两者在逼近过程中协同处理。作者对两者采用类似逼近策略，保证极限中仍维持物理约束 $\phi \in [-1,1], \psi \in [0,1]$。

  4. 恢复压力项
    在弱解框架下，压力 $\pi$ 通常作为 Lagrange 乘子出现。作者在附录中通过 Nečas 定理或 Bogovskiǐ 算子恢复其存在性，属于技术性但必要的步骤。

4. 与之前工作的比较

a) 突破与改进

本文相较于已有工作有多项显著推进：

  1. 推广至可溶性表面活性剂系统
    此前关于非均质两相流的工作（如 Abels–Garcke–Grün 模型 [4]）未考虑表面活性剂。本文首次将 AGG 模型拓展至含可溶性表面活性剂情形，增加了物理真实性与数学复杂性。

  2. 允许一般退化迁移率
    文献 [2] 仅处理线性退化迁移率 $m_\phi(s) = 1 - s^2$ 且使用正则势。本文允许更高阶退化（如 $(1-s^2)^k, k \geq 1$）并结合奇异势，更贴近物理实际（纯相中扩散趋于停止）。

  3. 采用物理相关奇异势而非正则近似
    多数前期工作为简化分析使用四次多项式势，无法保证 $\phi \in [-1,1]$。本文坚持使用 Flory–Huggins 型对数奇异势，确保解自动落在物理区间内，结果更具物理意义。

  4. 无需高阶正则化项
    相比于 [20] 中需添加 $\gamma |\Delta \phi|^2$ 来闭合估计，本文通过精细的能量与熵估计，在无额外正则化的前提下完成证明，说明模型本身已具备足够耗散性。

  5. 处理非恒定质量源项
    允许 $\sigma_1(\phi)$ 为非常数函数，增强了模型灵活性，可用于模拟非均匀化学反应过程。

  6. 统一处理二维与三维情形
    解的存在性在 $d=2,3$ 下均成立，而强解通常只在二维唯一。

综上，本文在模型完整性、数学假设宽松性、物理保真度三方面均有重要突破，为后续研究此类复杂多相流系统的适定性、渐近行为与数值方法奠定了坚实基础。

2. On the continuity in time and the exponential decay of solutions to a generalized Navier–Stokes–Fourier system in 2D

关于二维广义Navier-Stokes-Fourier系统解在时间上的连续性及指数衰减性

作者: Miroslav Bulíček, Petr Kaplický, Lucie Wintrová

类别: math.AP, 35Q30, 35K61, 35K92, 37L15, 76D03, 76E30

发布日期: 2025-10-24

arXiv ID: 2510.21218v1

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

详细分析

中文标题: 关于二维广义Navier-Stokes-Fourier系统解在时间上的连续性及指数衰减性

1. 研究对象和背景

本文研究的是一个描述二维不可压缩、非牛顿、热传导流体的广义 Navier–Stokes–Fourier (NSF) 系统。该系统由动量守恒、质量守恒、内能平衡以及熵方程构成，适用于具有温度依赖材料参数的复杂流体动力学建模。

所考虑的数学模型定义在一个有界 Lipschitz 域 $\Omega \subset \mathbb{R}^2$ 上，并在时间区间 $(0,T)$ 内演化。系统的未知函数包括：速度场 $uuu: Q \to \mathbb{R}^2$（其中 $Q := (0,T)\times\Omega$），压力 $\pi: Q \to \mathbb{R}$，温度 $\vartheta: Q \to \mathbb{R}$ 和熵 $\eta = \log\vartheta$。此外，记号 $Duuu := (\nabla uuu + (\nabla uuu)^T)/2$ 表示速度梯度的对称部分。

系统的基本控制方程如下：

$$ \partial_t uuu + \mathrm{div}(uuu \otimes uuu) - \mathrm{div}\, S + \nabla \pi = fff, $$
$$ \mathrm{div}\, uuu = 0, $$
$$ \partial_t \vartheta + \mathrm{div}(\vartheta uuu) + \mathrm{div}\, qqq = S : Duuu, $$
$$ \partial_t \eta + \mathrm{div}(\eta uuu) + \mathrm{div}\left( \frac{qqq}{\vartheta} \right) = - \frac{S : Duuu}{\vartheta} - \frac{qqq \cdot \nabla \vartheta}{\vartheta^2}. $$

这些方程在区域 $Q = (0,T)\times\Omega$ 中几乎处处成立。边界条件为：

  • 对于速度：齐次 Dirichlet 边界条件，即 $uuu = 0$ 在 $\partial\Omega\times(0,T)$；
  • 对于温度：非齐次 Dirichlet 边界条件，即 $\vartheta = \vartheta_b$ 在 $\partial\Omega\times(0,T)$。

初始条件设为：

$$ uuu(0) = uuu_0,\quad \vartheta(0) = \vartheta_0\quad \text{在 } \Omega \text{ 中}. $$

值得注意的是，尽管系统包含四个偏微分方程，但由于熵 $\eta = \log\vartheta$ 是温度的函数，且假设比热容恒为1，因此第(3)与第(4)个方程在光滑解情形下是形式等价的。然而，在弱解框架中，它们不再自动等价，必须分别处理。

本研究的关键物理特性体现在以下两个方面：

非牛顿流变行为

粘性 Cauchy 应力张量 $S$ 满足一个依赖于温度和应变率的幂律关系，即 $S = S_*(\vartheta, Duuu)$，满足以下结构条件（见原文(1.8)）：

  • 单调性：$(S_*(\vartheta,D_1)-S_*(\vartheta,D_2)):(D_1-D_2) \geq 0$，
  • 强增长性：$S_*(\vartheta,D):D \geq \nu|D|^p - \nu$，
  • 增长阶估计：$|S_*(\vartheta,D)| \leq \nu(1+|D|)^{p-1}$，

其中 $p \geq 2$，表明这是一个剪切增稠流体（shear-thickening fluid）。这类流体在高剪切速率下粘度增加，常见于泥浆、聚合物溶液等工业材料。

温度相关的热传导机制

热通量 $qqq$ 服从Fourier 定律：

$$ qqq = qqq_*(\vartheta) = -\kappa(\vartheta)\nabla \vartheta, $$

其中热导率 $\kappa(\vartheta)$ 是正的连续函数，并满足上下界：

$$ 0 < \underline{\kappa} \leq \kappa(\vartheta) \leq \overline{\kappa} < \infty. $$

这种设定允许热导率随温度变化，但避免了奇异性或退化过强的情形。

函数空间与边值条件的特殊性

论文采用的标准函数空间包括：

  • $V := W^{1,p}_{0,\mathrm{div}}(\Omega;\mathbb{R}^2)$：零迹、散度自由、Sobolev 类型的速度空间；
  • $L^2_{0,\mathrm{div}}$：其在 $L^2$ 中的闭包；
  • Bochner 空间如 $L^p(0,T;V)$、$C([0,T];L^2)$ 等用于描述时空正则性。

特别地，由于采用了齐次 Dirichlet 速度边界条件，结合 Korn 不等式，可有效控制速度的 Sobolev 范数。而温度边界条件虽为非齐次，但通过引入参考函数 $\hat\vartheta$（满足 $\mathrm{div}(\kappa(\hat\vartheta)\nabla\hat\vartheta)=0$ 且 $\hat\vartheta|_{\partial\Omega}=\vartheta_b$），将问题转化为关于 $\vartheta - \hat\vartheta$ 的齐次问题，从而便于分析。

此外，作者强调初始数据只需满足有限能量条件（如 $uuu_0 \in L^2_{0,\mathrm{div}}, \vartheta_0 \in L^1(\Omega)$），并未要求高阶正则性，这使得结果更具普适性。

综上所述，本文研究的是一个高度耦合、非线性、跨尺度的偏微分方程组，其难点在于：

  • 动量方程中的非线性应力项；
  • 能量方程中的耗散项 $S:Duuu$；
  • 温度与速度之间的强耦合；
  • 弱解框架下的初值实现问题；
  • 特别是对温度在时间方向上的连续性这一长期未被严格证明的性质。

这些问题在以往关于非牛顿热传导流体的研究中往往被忽略或仅以不等式形式处理（如熵不等式），而本文首次在二维情形下建立了满足熵等式的全局弱解，并进一步证明了其时间连续性。

2. 主要定理或主要结果

论文的核心成果总结于以下定理（原文 Theorem 1.1）：

定理 1.1（满足熵等式的解的存在性）
设 $\Omega \subset \mathbb{R}^2$ 为有界 Lipschitz 域，$T>0$。假设 $S_* = S_*(\vartheta, Duuu)$ 和 $\kappa(\vartheta)$ 满足条件 (1.7)–(1.8)，外力 $fff$、初值 $uuu_0, \vartheta_0$ 及边界温度诱导的参考函数 $\hat\vartheta$ 满足 (1.13)–(1.15)。

则存在四元组 $(uuu, S, \vartheta, \eta)$，满足以下正则性：

$$ uuu \in C([0,T]; L^2_{0,\mathrm{div}}) \cap L^p(0,T; V), $$
$$ \partial_t uuu \in L^{p'}(0,T; V^*),\quad S \in L^{p'}(Q; \mathbb{R}^{2\times2}_{\mathrm{sym}}), $$
$$ \vartheta \in C([0,T]; L^1(\Omega)),\quad (\vartheta)^\alpha \in L^2(0,T; W^{1,2}(\Omega)) \quad \forall \alpha \in (0,1/2), $$
$$ \vartheta \in L^r(Q) \quad \forall r \in [1,2), $$
$$ \vartheta - \hat\vartheta \in L^s(0,T; W^{1,s}(\Omega)) \quad \forall s \in [1,4/3), $$
$$ \eta = \log\vartheta \in L^2(0,T; W^{1,2}(\Omega)) \cap L^q(Q) \quad \forall q \in [1,\infty), $$

并且该四元组在如下意义下满足系统 (1.1)–(1.5)：

  • 动量方程：$S = S_*(\vartheta, Duuu)$ a.e. in $Q$，初值 $uuu(0)=uuu_0$ 成立于 $L^2$，且对所有测试函数 $www \in L^p(0,T; V)$ 有：

    $$ \int_0^T \langle \partial_t uuu, www \rangle \,dt + \iint_Q S : Dwww \,dxdt = \iint_Q (uuu\otimes uuu):Dwww\,dxdt + \int_0^T \langle fff, www \rangle \,dt; $$

  • 内能平衡：$\vartheta \geq \mu > 0$ a.e. in $Q$，初值 $\vartheta(0)=\vartheta_0$ 成立于 $L^1$，且对所有 $\varphi \in C^\infty_0((-\infty,T)\times\Omega)$ 有：

    $$ -\iint_Q \vartheta \partial_t \varphi \,dxdt - \iint_Q \vartheta uuu \cdot \nabla \varphi \,dxdt + \iint_Q \kappa(\vartheta)\nabla\vartheta \cdot \nabla\varphi \,dxdt = \iint_Q S:Duuu\,\varphi\,dxdt + \int_\Omega \vartheta_0 \varphi(0)\,dx; $$

  • 熵等式：$\eta = \log\vartheta$ a.e. in $Q$，且对所有 $\varphi \in C^\infty_0((-\infty,T)\times\Omega)$ 有：

    $$ -\iint_Q \eta \partial_t \varphi \,dxdt - \iint_Q \eta uuu \cdot \nabla \varphi \,dxdt + \iint_Q \kappa(\vartheta)\nabla\eta \cdot \nabla\varphi \,dxdt = \iint_Q \frac{S:Duuu}{\vartheta} \varphi \,dxdt + \iint_Q \frac{\kappa(\vartheta)|\nabla\vartheta|^2}{\vartheta^2} \varphi \,dxdt + \int_\Omega \eta_0 \varphi(0)\,dx. $$

核心创新点在于：

  1. 首次严格证明了温度 $\vartheta$ 在时间上属于 $C([0,T]; L^1(\Omega))$，即关于时间连续取值于 $L^1$；
  2. 解满足精确的熵等式而非通常的熵不等式；
  3. 利用上述两点构造了一个合适的 Lyapunov 泛函，进而推出稳态解的非线性稳定性，并最终得到所有合适弱解向稳态指数衰减的结果。

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

a) 核心分析学工具

本文综合运用了多种现代偏微分方程分析技术，主要包括：

  • Galerkin 近似法：用于构造近似解序列。选取 $W^{2,2}_{0,\mathrm{div}}$ 和 $L^2$ 的正交基，构建有限维 ODE 系统，保证局部解存在。
  • Aubin–Lions 引理：这是紧性论证的核心。利用速度和温度的时间导数在对偶空间中的有界性，结合空间嵌入，获得强收敛性。例如从 $\partial_t uuu^N \in L^{p'}((W^{2,2})^*)$ 和 $uuu^N \in L^p(V)$ 推出 $uuu^N \to uuu$ 强于 $L^r(Q)$。
  • 单调算子理论与 Minty 技巧：用于识别极限中的非线性项 $S = S_*(\vartheta, Duuu)$。通过不等式
    $$ \iint_Q (S_N - S_*(\vartheta, D)) : (Duuu^N - D) \,dxdt \geq 0 $$
    并令 $N\to\infty$，结合弱收敛和连续性，最终得出 $S = S_*(\vartheta, Duuu)$。
  • 能量估计与 Gronwall 不等式：通过对近似系统进行逐项测试（如用 $uuu^N$ 测试动量方程，用 $\vartheta^N - \hat\vartheta$ 测试能量方程），建立一致先验估计，确保解在整个时间区间上延拓。
  • 截断函数与凹变换技巧：为了处理低正则性的温度方程，作者引入形如 $\psi = 1 - \left( \frac{K}{K + \vartheta^N - \hat\vartheta} \right)^\alpha$ 的测试函数（其中 $K = \|\hat\vartheta\|_\infty$），从而导出关于 $(\vartheta^N)^{1-\alpha/2}$ 的梯度估计，这是获得 $L^2(W^{1,2})$ 正则性的关键。

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

  • 克服初值实现问题：传统方法常无法保证初值在强拓扑中实现。本文通过证明 $\vartheta \in C([0,T]; L^1)$，直接得出 $\lim_{t\to0^+} \|\vartheta(t) - \vartheta_0\|_{L^1} = 0$，解决了长期以来弱解中“初值是否可达”的难题。
  • 绕开经典重正规化解框架：以往工作（如 DiPerna-Lions）多依赖重正规化技巧处理测度初值或低正则解。本文因证明了熵等式和时间连续性，无需使用重正规化程序，极大简化了理论结构。
  • 构造 Lyapunov 泛函的基础——熵等式：由于解满足精确熵等式而非不等式，作者得以定义基于相对熵的 Lyapunov 泛函：
    $$ \mathcal{L}(t) := \int_\Omega \left( \frac{|uuu|^2}{2} + \vartheta - \bar\vartheta - \bar\vartheta \log\frac{\vartheta}{\bar\vartheta} \right) dx, $$
    其中 $\bar\vartheta$ 为稳态温度。利用熵等式可证 $\mathcal{L}'(t) \leq -c \mathcal{L}(t)$，从而推出指数衰减。
  • 精细的 Hölder 插值与梯度估计：为估计 $\|\nabla \vartheta^N\|_{L^t(Q)}$，作者使用 Hölder 不等式分解：
    $$ \int |\nabla \vartheta^N|^t = \int \frac{|\nabla \vartheta^N|^2}{\vartheta^{\alpha+1}} \cdot (\vartheta^N)^{t(\alpha+1)/2}, $$
    第一项由熵相关估计控制，第二项由 $L^r(Q)$ 估计插值得到，最终得 $t < 4/3$ 时的一致界。

4. 与之前工作的比较

a) 突破与改进

相比已有文献，本文做出了以下几个关键突破：

  • 首次证明二维 NSF 系统弱解的时间连续性：此前即使在 Newtonian 情形，也极少能证明 $\vartheta \in C([0,T]; L^1)$。本文通过熵等式与精细估计实现了这一点。
  • 从熵不等式到熵等式：多数现有工作（如 Consiglieri [10,11], Bulíček et al. [6,9]）只能证明熵不等式，丢失了信息。本文构造的解满足等式，为稳定性分析提供了基础。
  • 无需重正规化框架：Abbatiello 等人 [1,2] 虽在三维中建立熵等式解，但仍依赖复杂的逼近过程。本文在二维中因正则性更高，完全避开重正规化，路径更简洁。
  • 获得指数稳定性：这是本文最深刻的推论。以往结果多止步于存在性或渐近趋稳，本文利用 Lyapunov 方法证明了解以指数速率趋于稳态，这是对长期行为的强有力刻画。

b) 放宽假设与增强结论

  • 初始数据要求极低：仅需 $uuu_0 \in L^2$, $\vartheta_0 \in L^1$，远低于 $H^1$ 或 $W^{1,r}$ 要求。
  • 允许温度依赖的材料系数：$\kappa(\vartheta)$ 和 $S_*(\vartheta, D)$ 均可依赖温度，增强了模型现实意义。
  • 边界温度可非均匀：$\vartheta_b$ 不必为常数，只需 $\hat\vartheta \in W^{1,2} \cap L^\infty$，适用范围广。
  • 结论更强：不仅存在解，而且该解是时间连续的、满足熵等式、吸引其他解并指数衰减。

综上，本文在二维情形下完成了从“存在性”到“高质量解”再到“动态稳定性”的完整链条，标志着非牛顿热流体理论的一个重要进展。

3. Axisymmetric self-similar solutions to the MHD equations

磁流体动力学方程的轴对称自相似解

作者: Shaoheng Zhang

类别: math.AP

发布日期: 2025-10-24

arXiv ID: 2510.21194v1

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

详细分析

中文标题: 磁流体动力学方程的轴对称自相似解

1. 研究对象和背景

本文研究的是三维空间中 磁流体动力学（Magnetohydrodynamics, MHD）方程 的 轴对称自相似解。这类方程描述了导电流体（如等离子体）在磁场作用下的运动行为，是将 Navier-Stokes 方程 与 Maxwell 方程组 耦合而成的一类重要偏微分方程系统。

论文考虑的是 稳态不可压缩 MHD 方程，其形式如下：

$$ \begin{cases} -\Delta \mathbf{u} + (\mathbf{u} \cdot \nabla)\mathbf{u} - (\mathbf{B} \cdot \nabla)\mathbf{B} + \nabla p = 0, \\ -\Delta \mathbf{B} + (\mathbf{u} \cdot \nabla)\mathbf{B} - (\mathbf{B} \cdot \nabla)\mathbf{u} = 0, \\ \mathrm{div}\,\mathbf{u} = \mathrm{div}\,\mathbf{B} = 0, \end{cases} $$

其中 $\mathbf{u}$ 是速度场，$\mathbf{B}$ 是磁场，$p$ 是压力标量函数。该方程组定义在区域 $\Omega \subseteq \mathbb{R}^3$ 上。

特别地，作者关注两个几何情形：

  • 全空间去掉原点：$\mathbb{R}^3 \setminus \{0\}$；
  • 半空间：$\mathbb{R}^3_+ = \{(x_1,x_2,z) \in \mathbb{R}^3 \mid z > 0\}$。

在这两类区域上，分别附加不同的边界条件：

  • 在 $\mathbb{R}^3 \setminus \{0\}$ 中无边界，但要求解在原点具有奇异性；
  • 在 $\mathbb{R}^3_+$ 上则施加 无滑移边界条件 或 Navier 滑移边界条件，具体将在后文详述。

一个关键的数学结构是 尺度不变性（scaling invariance）。MHD 方程在如下变换下保持形式不变：

$$ \mathbf{u}(x) \mapsto \lambda \mathbf{u}(\lambda x),\quad \mathbf{B}(x) \mapsto \lambda \mathbf{B}(\lambda x),\quad p(x) \mapsto \lambda^2 p(\lambda x), \quad \lambda > 0. $$

若某个解 $(\mathbf{u}, \mathbf{B}, p)$ 满足上述变换对其自身成立，即对于所有 $\lambda > 0$ 都有 $\mathbf{u}(x) = \lambda \mathbf{u}(\lambda x)$，则称其为 自相似解（self-similar solution）。进一步，若解还满足关于 $z$-轴的旋转对称性，则称为 轴对称自相似解（axisymmetric self-similar solution）。

这类解在流体力学中具有重要意义，因为它们可能对应于某些奇点附近的渐近行为或全局解的渐近形态。例如，在 Navier-Stokes 方程中，著名的 Landau 解 就是一类非平凡的、在 $\mathbb{R}^3 \setminus \{0\}$ 上定义的轴对称自相似解，它由 L. Landau 在 1944 年发现，并被广泛研究。

本文的研究对象正是在上述背景下，探讨当磁场 $\mathbf{B}$ 只保留 旋量分量（swirl component） ——即仅含角向分量 $B^\theta(r,z)\mathbf{e}_\theta$——时，是否存在非平凡的轴对称自相似解。

更精确地说，设在柱坐标系 $(r,\theta,z)$ 下：

$$ \mathbf{u} = u^r(r,z)\mathbf{e}_r + u^\theta(r,z)\mathbf{e}_\theta + u^z(r,z)\mathbf{e}_z,\quad \mathbf{B} = B^\theta(r,z)\mathbf{e}_\theta, $$

且 $\mathbf{u}, \mathbf{B}$ 均为光滑的、在 $\mathbb{R}^3 \setminus \{0\}$ 或 $\mathbb{R}^3_+$ 上定义的 (−1)-齐次函数，即满足：

$$ \mathbf{u}(\lambda x) = \lambda^{-1} \mathbf{u}(x),\quad \mathbf{B}(\lambda x) = \lambda^{-1} \mathbf{B}(x). $$

这一齐次性保证了解的自相似性质。同时，由于只保留角向磁场分量，系统的复杂度有所降低，但仍保留了非线性耦合的本质特征。

在 $\mathbb{R}^3 \setminus \{0\}$ 的情形下，前人已有部分结果。例如 Zhang, Wang, and Wang [ZWW25] 曾证明：如果速度大小满足 $|\mathbf{u}(x)| \leq \varepsilon |x|^{-1}$ 且 $\varepsilon$ 足够小，则必有 $\mathbf{B} \equiv 0$，且 $\mathbf{u}$ 为 Landau 解。而 Zhang [Zha25] 则去除了小性假设，但在没有磁扩散项的前提下得出相同结论。

本文的主要突破在于：通过引入一个新的上界条件——径向速度 $u^r$ 在单位球面 $\partial B_1$ 上满足

$$ u^r(r,z) < \frac{1}{3r} + \frac{2r}{3}, $$

即可推出 $\mathbf{B} \equiv 0$ 且 $\mathbf{u}$ 为 Landau 解。这个条件比之前的小范数假设更弱，更具几何意义。

此外，在半空间 $\mathbb{R}^3_+$ 上，作者进一步研究了带有边界的自相似解。此时需附加适当的边界条件：

  • 无滑移边界条件：
    $$ \mathbf{u} = 0,\quad \mathbf{B} \cdot \mathbf{n} = 0,\quad \mathrm{curl}\,\mathbf{B} \times \mathbf{n} = 0,\quad \text{on } \partial\mathbb{R}^3_+ \setminus \{0\}; $$
  • Navier 滑移边界条件：
    $$ \mathbf{u} \cdot \mathbf{n} = 0,\quad (D\mathbf{u} \cdot \mathbf{n}) \cdot \boldsymbol{\tau} = 0,\quad \mathbf{B} \cdot \mathbf{n} = 0,\quad \mathrm{curl}\,\mathbf{B} \times \mathbf{n} = 0, $$
    其中 $\mathbf{n} = -\mathbf{e}_z$ 是外法向，$\boldsymbol{\tau}$ 是切向向量，$D\mathbf{u} = \frac{1}{2}(\nabla \mathbf{u} + \nabla \mathbf{u}^T)$ 是应变张量。

这些边界条件体现了物理上的合理性：磁场满足无源通量和旋度约束，速度满足不同层次的滑移机制。

综上所述，本文的研究对象是一个具有尺度不变性和轴对称性的稳态 MHD 系统，其解允许在原点或边界处出现奇点，目标是刻画在此类结构限制下唯一可能存在的自相似解类型。

2. 主要定理或主要结果

定理 1.1

设 $(\mathbf{u}, \mathbf{B})$ 是定义在 $\mathbb{R}^3 \setminus \{0\}$ 上的光滑轴对称自相似解，满足：

  • $\mathbf{u} = u^r(r,z)\mathbf{e}_r + u^\theta(r,z)\mathbf{e}_\theta + u^z(r,z)\mathbf{e}_z$，
  • $\mathbf{B} = B^\theta(r,z)\mathbf{e}_\theta$，
  • 径向速度在单位球面上满足严格不等式：
    $$ u^r(r,z) < \frac{1}{3r} + \frac{2r}{3},\quad \forall (r,z) \in \partial B_1. $$

则结论为：

$$ \mathbf{u} \text{ 是一个 Landau 解},\quad \mathbf{B} \equiv 0. $$

此定理表明，在给定对称性和磁场结构下，只要径向速度在单位球上被适当控制，磁场就必须消失，速度退化为经典的 Landau 解。

推论 1.1

若在 $\partial B_1$ 上有 $u^r(x) \leq \alpha$，且常数 $\alpha < \frac{2\sqrt{2}}{3}$，则同样可得 $\mathbf{u}$ 为 Landau 解，$\mathbf{B} \equiv 0$。

注意：函数 $\frac{1}{3r} + \frac{2r}{3}$ 在 $\partial B_1$（即 $r^2 + z^2 = 1$）上的最小值恰好为 $\frac{2\sqrt{2}}{3}$，因此该推论是定理 1.1 的直接应用。

定理 1.2

设 $(\mathbf{u}, \mathbf{B})$ 是定义在半空间 $\mathbb{R}^3_+$ 上的光滑轴对称自相似解，满足无滑移或 Navier 滑移边界条件，且 $\mathbf{u}, \mathbf{B}$ 形式同上。若在 $\partial B_1 \cap \mathbb{R}^3_+$ 上满足

$$ u^r(r,z) < \frac{1}{3r} + \frac{2r}{3}, $$

则必有：

$$ \mathbf{u} = 0,\quad \mathbf{B} = 0. $$

这说明在半空间中，即使允许自相似奇点存在，也无法构造非平凡的此类解，即 唯一可能是零解。

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

a) 核心分析工具

本文的核心技术依赖于 能量型积分估计 和 齐次函数的几何恒等式，结合 分部积分 与 coarea 公式 来处理自相似结构带来的缩放对称性。

最关键的一步是对角向磁场分量 $B^\theta$ 所满足的方程进行加权积分。从原始 MHD 方程出发，利用轴对称性和仅有 $B^\theta$ 分量的设定，可推导出 $B^\theta$ 满足的简化方程：

$$ -\Delta B^\theta + (\mathbf{u} \cdot \nabla) B^\theta + \frac{1}{r^2} B^\theta - \frac{u^r B^\theta}{r} = 0. $$

这是通过对 MHD 第二个方程投影到 $\mathbf{e}_\theta$ 方向并利用柱坐标下 Laplace 算子表达式得到的。

接下来，选取测试函数 $r B^\theta$，并在环形区域 $\Omega = B_{\rho_2} \setminus \overline{B}_{\rho_1}$ 上乘以上述方程并积分，得到恒等式：

$$ \int_\Omega \left[ (-\Delta B^\theta) r B^\theta + r (\mathbf{u} \cdot \nabla)\left( \frac{1}{2} |B^\theta|^2 \right) + \frac{1}{r} |B^\theta|^2 - u^r |B^\theta|^2 \right] dx = 0. $$

随后使用标准的 分部积分技巧 处理每一项。特别是对拉普拉斯项：

$$ \int_\Omega (-\Delta B^\theta)(r B^\theta)\,dx = \int_\Omega r |\nabla B^\theta|^2 dx + \int_{\partial \Omega} \left( -\frac{1}{2} \partial_\nu |B^\theta|^2 + r B^\theta \partial_\nu B^\theta \right) dS. $$

但由于 $B^\theta$ 是 $(-1)$-齐次的，满足：

$$ x \cdot \nabla B^\theta(x) = -B^\theta(x), $$

从而可以计算边界积分为零（利用 coarea 公式和齐次性），进而消去边界项。

另一个关键恒等式是：

$$ \int_\Omega \partial_r |B^\theta|^2 dx = -\int_\Omega \frac{1}{r} |B^\theta|^2 dx, $$

这是通过散度定理和 $\mathrm{div}( |B^\theta|^2 \mathbf{e}_r )$ 的展开得到的，依赖于轴对称性。

最终整合所有项后，得到主恒等式：

$$ \int_\Omega \left( r |\nabla B^\theta|^2 + \frac{1}{2r} |B^\theta|^2 - \frac{3}{2} u^r |B^\theta|^2 \right) dx = 0. $$

再利用梯度估计：

$$ |\nabla B^\theta(x)| \geq \frac{|B^\theta(x)|}{|x|}, $$

代入得：

$$ \int_\Omega \left( \frac{r}{|x|^2} + \frac{1}{2r} - \frac{3}{2} u^r \right) |B^\theta|^2 dx \leq 0. $$

最后应用 coarea 公式 将体积积分转化为球面积分：

$$ \int_{\rho_1}^{\rho_2} \frac{1}{\rho} \left( \int_{\partial B_\rho} \left( \frac{r}{\rho^2} + \frac{1}{2r} - \frac{3}{2} u^r \right) |B^\theta|^2 dS \right) d\rho \leq 0. $$

由于被积函数在每个球面上连续，且若在 $\partial B_1$ 上有

$$ \frac{r}{1} + \frac{1}{2r} - \frac{3}{2} u^r > 0, $$

即等价于 $u^r < \frac{1}{3r} + \frac{2r}{3}$，则整个积分非负，只有当 $B^\theta \equiv 0$ 时才能为零。

因此，正性条件导致矛盾，除非 $B^\theta = 0$。

一旦 $B^\theta \equiv 0$，则 $\mathbf{B} \equiv 0$，MHD 系统退化为 Navier-Stokes 方程。此时已知结果（如 Tian-Xin、Cannone-Karch、Sverák）表明：在 $\mathbb{R}^3 \setminus \{0\}$ 上的轴对称自相似解必为 Landau 解。

在半空间情形，边界项虽不再为零，但借助边界条件仍可证明其贡献为零。例如，由 $\mathrm{curl}\,\mathbf{B} \times \mathbf{n} = 0$ 推出 $\partial_z B^\theta = 0$；由无滑移或滑移条件推出 $u \cdot \mathbf{n} = 0$ 或切向应力为零，从而使得边界积分为零。

因此，整体框架可概括为：

通过精心选择加权乘子 $r B^\theta$，结合齐次性、轴对称性与边界条件，建立一个包含 $B^\theta$ 的非负能量泛函，其符号由 $u^r$ 控制；一旦 $u^r$ 被适当上界限制，则强制 $B^\theta \equiv 0$，进而归约至 Navier-Stokes 情形。

b) 技术难点与解决技巧

  • 非线性耦合项的处理：$(\mathbf{u} \cdot \nabla)\mathbf{B}$ 与 $(\mathbf{B} \cdot \nabla)\mathbf{u}$ 的交叉项通常难以控制。但本文利用了磁场仅有 $B^\theta$ 分量的特殊结构，使得许多高阶项自动消失，大大简化了方程。

  • 奇点处理：解在原点不光滑，甚至发散（如 $O(|x|^{-1})$）。然而，由于齐次性，可通过 coarea 公式将局部信息“平均”到单位球上，避免直接分析奇点邻域。

  • 边界效应：在半空间中，边界积分不可忽略。作者巧妙利用边界条件（尤其是 $\partial_z B^\theta = 0$ 和 $u_z = 0$）来抵消边界项，体现了对物理边界条件的深刻理解。

  • 比较函数的构造：虽然未显式构造 Lyapunov 泛函，但本质上使用了一个类似于“加权 Dirichlet 能量”的量，并通过其正定性得出唯一性。

4. 与之前工作的比较

a) 改进与突破

本文相较于先前工作（特别是 [ZWW25] 和 [Zha25]）做出了显著改进：

  1. 放宽了对速度的整体小性假设：

    • 文献 [ZWW25] 要求在整个 $\partial B_1$ 上 $|\mathbf{u}| \leq \varepsilon$ 且 $\varepsilon$ 足够小，这是一种很强的范数控制；
    • 本文仅要求 径向分量 $u^r$ 满足一个显式的上界 $u^r < \frac{1}{3r} + \frac{2r}{3}$，这是一个逐点条件，且允许速度其他分量较大。

    这意味着新的条件 更弱、更具几何意义，因为它只限制特定方向的速度增长速率。

  2. 无需小参数论证：

    • 以往结果往往依赖于摄动方法或隐含常数的存在性；
    • 本文给出的是 显式可验证的判据，并且临界值 $\frac{2\sqrt{2}}{3} \approx 0.9428$ 是明确的，增强了结果的应用价值。

  3. 推广到边界情形：

    • 此前关于 MHD 自相似解的研究多集中于全空间；
    • 本文首次在 半空间 中讨论此类问题，并证明在自然边界条件下，不存在非平凡轴对称自相似解。

    这一结果呼应了 Navier-Stokes 情形下 Kang-Miura-Tsai 的工作，将其推广到了 MHD 情形。

  4. 去除了对磁扩散项的依赖：

    • 虽然 [Zha25] 已在无磁扩散情形下得到类似结论，但本文在有扩散项的情况下仍能得出相同结果，说明磁场衰减是由对流主导而非耗散机制。

b) 结果强度对比

工作区域对 $\mathbf{u}$ 的假设对 $\mathbf{B}$ 的结构主要结论
[ZWW25]$\mathbb{R}^3\setminus\{0\}$$\mathbf{u}\leq \varepsilon \ll 1$
[Zha25]$\mathbb{R}^3\setminus\{0\}$无小性假设（无磁扩散）同上同上
本文$\mathbb{R}^3\setminus\{0\}$$u^r < \frac{1}{3r}+\frac{2r}{3}$ on $\partial B_1$同上同上，但条件更弱
本文$\mathbb{R}^3_+$同上 + 边界条件同上$\mathbf{u}=\mathbf{B}=0$

由此可见，本文不仅 统一并强化了现有理论，而且将研究范围扩展到了带边区域，填补了 MHD 自相似解理论的一个空白。

更重要的是，它揭示了一个深层现象：即使允许磁场存在，只要其结构受限（仅旋量），且速度径向分量不过大，则磁场必然湮灭。这暗示了 MHD 系统中磁场与速度之间的强耦合抑制机制。

综上，本文在技术上采用了简洁而深刻的能量估计方法，在理论上推进了对自相似解分类的理解，是一篇具有基础性意义的工作。

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

1. A general Frobenius’ Theorem via the Transport of Currents

通过流的输运推广弗罗贝尼乌斯定理

作者: Paolo Bonicatto

类别: math.AP, 49Q15, 35Q49

发布日期: 2025-10-24

arXiv ID: 2510.21478v1

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

摘要翻译

中文标题: 通过流的输运推广弗罗贝尼乌斯定理

摘要翻译: 微分几何中的一个经典结果指出，两个光滑向量场的流相容当且仅当它们的李括号为零。本文将这一结果推广到更一般的情形：其中一个向量场是有界且满足Lipschitz条件的，而另一个可以是奇异的向量值测度，即一个正规的1-流。该结果是通过对两种不同的演化型偏微分方程——描述向量量输运的向量对流方程与几何输运方程——的研究而获得的。此外，我们指出磁流体力学中阿尔文（Alfvén）的一个著名定理可被解释为弗罗贝尼乌斯定理的一个适当的时间依赖版本。本研究的方法建立在近年来关于流的几何输运方程研究进展的基础之上[5, 6]。

2. Semiclassical limit of cubic nonlinear Schrödinger equations for mixed states

混合态下三次非线性薛定谔方程的半经典极限

作者: Daniel Han-Kwan, Frédéric Rousset

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

发布日期: 2025-10-24

arXiv ID: 2510.21313v1

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

摘要翻译

中文标题: 混合态下三次非线性薛定谔方程的半经典极限

摘要翻译: 本文研究了混合态下三次非线性薛定谔方程的半经典极限。我们证明了该极限可趋近于一个奇异的Vlasov方程（其中力场与密度梯度成正比），适用于具有有限Sobolev正则性的初值数据，且其速度分布满足量子Penrose稳定性条件。这一稳定性条件在聚焦和散焦情形下对于小初值数据均成立（小性条件不依赖于半经典参数）；此外，在散焦情形下，对于一大类物理上重要的例子（如局部Maxwellian型分布）的小扰动也满足该条件。

3. On the Propulsion of a Rigid Body in a Viscous Liquid by Time-Periodic Force with a Zero Average

在零平均周期力作用下刚体于粘性液体中的推进机制研究

作者: Joris Edelmann, Giovanni P. Galdi, Mher M. Karakouzian, Thomas Richter

类别: math.AP, physics.flu-dyn

发布日期: 2025-10-23

arXiv ID: 2510.20982v1

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

摘要翻译

中文标题: 在零平均周期力作用下刚体于粘性液体中的推进机制研究

摘要翻译: 本文对在周期性外力（其在一个周期内的平均值为零）作用下，刚体在粘性流体中推进的现象进行了理论分析与数值模拟。该通用理论框架特别适用于一类重要情形，即推进力由位于刚体内部空腔中的质量振荡所产生。我们严格证明了在力的量级二阶近似下实现推进的必要且充分条件。对于不具有前后对称性的物体，这些条件可通过数值实验加以实施并得到验证；而对于具有前后对称性的物体（如旋转椭球），这些条件不提供任何推进效应，意味着此类情形下的推进仅可能出现在高于二阶的近似阶次中。针对后一情况，本文通过数值求解完整的控制方程组进行研究，结果表明推进现象确实发生，从而为后续的理论分析开辟了新的研究方向。

不相关论文（基本信息）

1. Recovering Initial States in Certain Quasilinear Parabolic Problems from Time Averages

从时间平均值中恢复某些拟线性抛物型问题的初始状态

作者: Bogdan-Vasile Matioc, Christoph Walker

类别: math.AP

发布日期: 2025-10-24

arXiv ID: 2510.21687v1

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

摘要: The inverse problem of reconstructing the initial state in quasilinear parabolic equations from time averages is investigated. Under suitable regularity assumptions on the quasilinear structure and a superlinear growth condition near zero for the semilinear part, it is shown that the initial state can be uniquely recovered from small time averages taken over an arbitrary time period. The applicability of the result is demonstrated for certain chemotaxis models and reaction-diffusion systems.

2. Local boundedness of weak solutions to elliptic equations under unbalanced Orlicz growth conditions

在非平衡Orlicz增长条件下椭圆型方程弱解的局部有界性

作者: Gabriele Giannone

类别: math.AP, 35J60 (Primary) 46E30 (Secondary)

发布日期: 2025-10-21

arXiv ID: 2510.18979v2

链接: http://arxiv.org/abs/2510.18979v2

摘要: We establish sufficient conditions for the local boundedness of weak solutions to a broad class of nonlinear elliptic equations in divergence form, under unbalanced growth conditions on the stress field. Our analysis is carried out in a non-variational setting, with no symmetry or structural assumptions on the operator. The ellipticity and growth are prescribed via distinct Young functions, leading to a Orlicz-type setting that captures a wide class of nonstandard behaviors. As a special case, the theory encompasses and extends the known results on equations with $p,q$-growth.

3. Describing smooth small-data solutions to a quasilinear hyperbolic-parabolic system by $W^{1,p}$ energy analysis

通过 $W^{1,p}$ 能量分析描述拟线性双曲-抛物系统的光滑小初值解

作者: Leander Claes, Michael Winkler

类别: math.AP, 35B40

发布日期: 2025-10-24

arXiv ID: 2510.21660v1

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

摘要: In bounded $n$-dimensonal domains with $n\ge 1$, this manuscript examines an initial-boundary value problem for the system

$$ \left\{ \begin{array}{l} u_{tt} = \nabla \cdot (\gamma(\Theta) \nabla u_t) + a \nabla \cdot (\gamma(\Theta) \nabla u) + \nabla\cdot f(\Theta), \Theta_t = D\Delta\Theta + \Gamma(\Theta) |\nabla u_t|^2 + F(\Theta)\cdot \nabla u_t, \end{array} \right. $$

which in the case $n=1$ and with $\gamma\equiv \Gamma$ as well as $f\equiv F$ reduces to the classical model for the evolution of strains and temperatures in thermoviscoelasticity. Unlike in previous related studies, the focus here is on situations in which besides $f$ and $F$, also the core ingredients $\gamma$ and $\Gamma$ may depend on the temperature variable $\Theta$. Firstly, a statement on local existence of classical solutions is derived for arbitrary $a>0, D>0$ as well as $0<\gamma\in C^2([0,\infty))$ and $0\le\Gamma\in C^1([0,\infty))$, for functions $f\in C^2([0,\infty);{\mathbb{R}}^n)$ and $F\in C^1([0,\infty);{\mathbb{R}}^n)$ with $F(0)=0$, and for suitably regular initial data of arbitrary size. Secondly, it is seen that for each $p\ge 2$ such that $p>n$ there exists $\delta(p)>0$ with the property that whenever in addition to the above we have

$$ \frac{a}{\gamma(0)} \le \delta(p) \qquad \mbox{and} \qquad \frac{|f'(\Theta_\star)| \cdot |F(\Theta_\star)|}{D \cdot \gamma(\Theta_\star)} \le \delta(p), $$

for initial data suitably close to the constant level given by $u=0$ and $\Theta=\Theta_\star$, with any fixed $\Theta_\star\ge 0$, these solutions are actually global in time and have the property that $\nabla u_t, \nabla u$ and $\nabla\Theta$ decay exponentially fast in $L^p$. This is achieved by detecting suitable dissipative properties of functionals involving norms of these gradients in $L^p$ spaces.

4. Nash’s $G$ bound for the Kolmogorov equation

纳什关于柯尔莫哥洛夫方程的 $G$ 界

作者: Helge Dietert, Lukas Niebel

类别: math.AP

发布日期: 2025-10-24

arXiv ID: 2510.21621v1

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

摘要: We prove Nash’s $G$ bound for the Kolmogorov equation with rough coefficients. Our proof is inspired by the treatment of the parabolic problem by Nash (1958) and Fabes and Stroock (1986). To transfer their ideas to the kinetic setting, we employ critical kinetic trajectories. From Nash’s $G$ bound, we recover the sharp lower bound on the fundamental solution and thus provide an alternative proof of the Harnack inequality for the Kolmogorov equation.

5. Multilevel Picard scheme for solving high-dimensional drift control problems with state constraints

求解带状态约束的高维漂移控制问题的多层Picard格式

作者: Yuan Zhong

类别: math.OC, cs.NA, math.AP, math.NA, math.PR

发布日期: 2025-10-24

arXiv ID: 2510.21607v1

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

摘要: Motivated by applications to the dynamic control of queueing networks, we develop a simulation-based scheme, the so-called multilevel Picard (MLP) approximation, for solving high-dimensional drift control problems whose states are constrained to stay within the nonnegative orthant, over a finite time horizon. We prove that under suitable conditions, the MLP approximation overcomes the curse of dimensionality in the following sense: To approximate the value function and its gradient evaluated at a given time and state to within a prescribed accuracy $\varepsilon$, the computational complexity grows at most polynomially in the problem dimension $d$ and $1/\varepsilon$. To illustrate the effectiveness of the scheme, we carry out numerical experiments for a class of test problems that are related to the dynamic scheduling problem of parallel server systems in heavy traffic, and demonstrate that the scheme is computationally feasible up to dimension at least $20$.

6. The Gross-Pitaewskii equation with time and space dependent coefficients

具有时空依赖系数的Gross-Pitaevskii方程

作者: Federico Lai

类别: math.AP

发布日期: 2025-09-07

arXiv ID: 2509.06001v2

链接: http://arxiv.org/abs/2509.06001v2

摘要: We study the existence of weak solutions of a generalized Gross-Pitaewskii equation, with time and space dependent coefficients that could blow up or vanish asymptotically in time, with initial data not necessarily segregated. We also study the asymptotic behavior in time of these solutions, including cases in which the segregation phenomenon appears.

7. On the uniqueness of even $L^p$ Minkowski problem

关于偶 $L^p$ Minkowski 问题解的唯一性

作者: Weiyong He, Junbang Liu

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

发布日期: 2025-10-24

arXiv ID: 2510.21530v1

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

摘要: We prove that there is a unique $p_0\in [0,1)$, which can be characterized by the eigenvalue of Hilbert operator related to a convex body, that the even $L^p$ Minkowski problem has a unique solution for $p\geq p_0$, and the uniqueness fails for infinitely many convex bodies if $p<p_0$. The previous results by many experts in the field assert that the uniqueness holds for $p>p_0$.

8. The asymptotic behavior of solutions to a doubly degenerate chemotaxis-consumption system in the two-dimensional setting

二维情形下一类双重退化趋化-消耗系统解的渐近行为

作者: Duan Wu

类别: math.AP

发布日期: 2024-09-18

arXiv ID: 2409.12083v2

链接: http://arxiv.org/abs/2409.12083v2

摘要: The present work proceeds to consider the convergence of the solutions to the following doubly degenerate chemotaxis-consumption system \begin{align*} \left{ \begin{array}{r@{,}l@{\quad}l@{,}c} &u_{t}=\nabla\cdot\big(u^{m-1}v\nabla v\big)-\nabla\cdot\big(f(u)v\nabla v\big)+\ell uv,\ &v_{t}=\Delta v-uv, \end{array}\right.%} \end{align*} under no-flux boundary conditions in a smoothly bounded convex domain $\Omega\subset \R^2$, where the nonnegative function $f\in C^1([0,\infty))$ is asked to satisfy $f(s)\le C_fs^{\al}$ with $\al, C_f>0$ for all $s\ge 1$. The global existence of weak solutions or classical solutions to the above system has been established in both one- and two-dimensional bounded convex domains in previous works. However, the results concerning the large time behavior are still constrained to one dimension due to the lack of a Harnack-type inequality in the two-dimensional case. In this note, we complement this result by using the Moser iteration technique and building a new Harnack-type inequality.

9. Constant sign and nodal solutions for singular quasilinear elliptic systems

奇异拟线性椭圆系统的一类恒号解与节点解

作者: Nouredine Medjoudj, Abdelkrim Moussaoui

类别: math.AP

发布日期: 2025-10-24

arXiv ID: 2510.21489v1

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

摘要: We consider singular quasilinear elliptic systems with homogeneous Dirichlet boundary condition. Using Leray-Schauder topological degree, combined with the sub-supersolutions method and suitable truncation arguments, we establish the existence of at least three nontrivial solutions, two of which are of opposite constant sign. The third solution is nodal and exhibits components of at least opposite constant sign. In the case of a sign-coupled system, these components are of changing and synchronized sign.

10. Higher Hölder regularity for fractional $(p,q)$-Laplace equations

分数阶$(p,q)$-Laplace方程的高阶Hölder正则性

作者: Prashanta Garain, Erik Lindgren

类别: math.AP

发布日期: 2025-09-19

arXiv ID: 2509.15988v2

链接: http://arxiv.org/abs/2509.15988v2

摘要: We study the fractional $(p,q)$-Laplace equation

$$ (-\Delta_p)^s u +(-\Delta_q)^t u= 0 $$

for $s,t\in(0,1)$ and $p,q\in(1,\infty)$. We establish H"older estimates with an explicit exponent. As a consequence, we derive a Liouville-type theorem. Our approach builds on techniques previously developed for the fractional $p$-Laplace equation, relying on a Moser-type iteration for difference quotients.

11. Interior Hessian estimates for sum Hessian quotient equation

关于和曲率商方程的内部分部Hessian估计

作者: Changyu Ren, Ziyi Wang

类别: math.AP, 35B45

发布日期: 2025-10-24

arXiv ID: 2510.21301v1

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

摘要: This paper is devoted to the interior $C^2$ estimates for a class of sum Hessian quotient equations. For $0\leq l<k<n$, we establish the interior estimates and the Pogorelov type estimates. In the case $k=n$, we obtain a weaker Pogorelov type estimate for $0\leq l<n-1$.

12. Infinitely many saturated travelling waves for a degenerate Fisher-KPP equation not in divergence form

非散度形式退化Fisher-KPP方程的无穷多个饱和行波解

作者: Matthieu Alfaro, Maxime Herda, Andrea Natale

类别: math.AP, 35K65, 35C07, 92D30

发布日期: 2025-02-17

arXiv ID: 2502.11589v2

链接: http://arxiv.org/abs/2502.11589v2

摘要: We consider an epidemic model with distributed-contacts. When the contact kernel concentrates, one formally reaches a very degenerate Fisher-KPP equation with a diffusion term that is not in divergence form. We make an exhaustive study of its travelling waves. For every admissible speed, there exist not only a unique non-saturated (smooth) wave but also infinitely many saturated (sharp) ones. Furthermore their tails may differ from what is usually expected. These results are thus in sharp contrast with their counterparts on related models.

13. Existence and stability of curved fronts for spatially periodic combustion reaction-diffusion equations in $\mathbb{R}^N$

$\mathbb{R}^N$ 中具空间周期性的燃烧反应扩散方程的曲面前端的存在性与稳定性

作者: Wei-Jie Sheng, Xin-Tian Zhang

类别: math.AP, 35B10, 35K10, 35K15, 35K57

发布日期: 2025-10-24

arXiv ID: 2510.21163v1

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

摘要: This paper is concerned with curved fronts of combustion reaction-diffusion equations in spatially periodic media in $\mathbb{R}^N$ $(N\geq2)$. Under the assumption that there are moving pulsating fronts for any given propagation direction $e \in \mathbb{S}^{N-1}$, and by constructing suitable super- and sub-solutions, we prove the existence of a curved front with polytope-like shape in $\mathbb{R}^N$. Then we show that the curved front is unique and asymptotically stable.

14. Propagation of Singularities for the Damped Stochastic Klein-Gordon Equation

阻尼随机Klein-Gordon方程奇性的传播

作者: Hongyi Chen, Cheuk Yin Lee

类别: math.PR, math.AP

发布日期: 2025-08-19

arXiv ID: 2508.13671v2

链接: http://arxiv.org/abs/2508.13671v2

摘要: For the $1+1$ dimensional damped stochastic Klein-Gordon equation, we show that random singularities associated with the law of the iterated logarithm exist and propogate in the same way as the stochastic wave equation. This provides evidence for possible connections to microlocal analysis, ie. the exact regularity and singularities described in this paper should admit wavefront set type descriptions whose propagation is determined by the highest order terms of the linear operator. Despite the results being exactly the same as those of the wave equation, our proofs are significantly different than the proofs for the wave equation. Miraculously, proving our results for the critically damped equation implies them for the general equation, which significantly simplifies the problem. Even after this simplification, many important parts of the proof are significantly different than (and we think are more intuitive from the PDE viewpoint compared to) existing proofs for the wave equation.

15. $L_p$-estimates of the conormal derivative problem for parabolic equations with time measurable coefficients and $A_p$-weights

具有时变可测系数和 $A_p$ 权的抛物型方程共法向导数问题的 $L_p$ 估计

作者: Hongjie Dong, Pilgyu Jung, Doyoon Kim

类别: math.AP

发布日期: 2025-10-24

arXiv ID: 2510.21139v1

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

摘要: This paper investigates weighted mixed-norm estimates for divergence-type parabolic equations on Reifenberg-flat domains with the conormal derivative boundary condition. The leading coefficients are assumed to be merely measurable in the time variable and to have small mean oscillations in the spatial variables. In deriving the boundary estimates, we overcome a regularity issue by employing half-time derivative estimates.

16. Regularity and structure of non-planar $p$-elasticae

非平面 $p$-弹性曲线的正则性与结构

作者: Florian Gruen, Tatsuya Miura

类别: math.AP, math.DG, 49Q10, 53A04, and 53E40

发布日期: 2025-01-14

arXiv ID: 2501.07987v3

链接: http://arxiv.org/abs/2501.07987v3

摘要: We prove regularity and structure results for $p$-elasticae in $\mathbb{R}^n$, with arbitrary $p\in (1,\infty)$ and $n\geq2$. Planar $p$-elasticae are already classified and known to lose regularity. In this paper, we show that every non-planar $p$-elastica is analytic and three-dimensional, with the only exception of flat-core solutions of arbitrary dimensions. Subsequently, we classify pinned $p$-elasticae in $\mathbb{R}^n$ and, as an application, establish a Li-Yau type inequality for the $p$-bending energy of closed curves in $\mathbb{R}^n$. This extends previous works for $p=2$ and $n\geq2$ as well as for $p\in (1,\infty)$ and $n=2$.

17. Asymptotic behavior at infinity and existence of solutions to the Lagrangian mean curvature flow in $\mathbb R^{n+1}_-$

关于 $\mathbb{R}^{n+1}_-$ 中拉格朗日平均曲率流解的无穷远处渐近行为及存在性研究

作者: Jiguang Bao, Zixiao Liu

类别: math.AP, 35B40, 35A01, 35K55, 53E10

发布日期: 2025-07-22

arXiv ID: 2507.16129v2

链接: http://arxiv.org/abs/2507.16129v2

摘要: This paper investigates the asymptotic behavior at infinity of ancient solutions to the Lagrangian mean curvature flow. Under conditions that admit Liouville type rigidity theorems, we prove that every classical solution converges at infinity to the sum of a quadratic polynomial in $x$ and a linear function in $t$, with an explicitly derived exponential rate of convergence. As a critical part of the proof framework of this paper, we establish the existence of a global viscosity solution with prescribed asymptotic behavior at infinity, featuring two key innovations: (i) applicability to all dimensions $n\geq 2$, and (ii) no requirement that the Hessian matrix of the prescribed quadratic term be positive definite or close to a scalar multiple of the identity matrix. These results establish the relationship between Liouville type rigidity, asymptotic analysis at infinity, and the existence of viscosity solutions.

18. Low energy resolvent asymptotics of the multipole Aharonov–Bohm Hamiltonian

多极Aharonov-Bohm哈密顿量的低能预解核渐近性

作者: T. J. Christiansen, K. Datchev, M. Yang

类别: math.AP, math-ph, math.MP, math.SP

发布日期: 2024-08-06

arXiv ID: 2408.03233v2

链接: http://arxiv.org/abs/2408.03233v2

摘要: We compute low energy asymptotics for the resolvent of the Aharonov–Bohm Hamiltonian with multiple poles for both integer and non-integer total fluxes. For integral total flux we reduce to prior results in black-box scattering while for non-integral total flux we build on the corresponding techniques using an appropriately chosen model resolvent. The resolvent expansion can be used to obtain long-time wave asymptotics for the Aharonov–Bohm Hamiltonian with multiple poles. An interesting phenomenon is that if the total flux is an integer then the scattering resembles even-dimensional Euclidean scattering, while if it is half an odd integer then it resembles odd-dimensional Euclidean scattering. The behavior for other values of total flux thus provides an `interpolation’ between these.

19. A uniform-in-time nonlocal approximation of the standard Fokker-Planck equation

标准福克-普朗克方程的一致非局部逼近

作者: José A. Cañizo, Niccolò Tassi

类别: math.AP, math.PR, 35B40, 35C10, 35R09, 45K05, 60F15, 60J76

发布日期: 2024-07-04

arXiv ID: 2407.03870v3

链接: http://arxiv.org/abs/2407.03870v3

摘要: We study a nonlocal approximation of the Fokker-Planck equation in which we can estimate the speed of convergence to equilibrium in a way which does not degenerate as we approach the local limit of the equation. This uniform estimate cannot be easily obtained with standard inequalities or entropy methods, but can be obtained through the use of Harris’s theorem, finding interesting links to quantitative versions of the central limit theorem in probability. As a consequence one can prove that solutions of this nonlocal approximation converge to solutions of the usual Fokker-Planck equation uniformly in time-hence we show the approximation is asymptotic-preserving in this sense. The associated equilibrium has some interesting tail and regularity properties, which we also study.

20. On the surjectivity of the Cauchy-Riemann and Laplace operators on weighted spaces of smooth functions

关于Cauchy-Riemann算子和Laplace算子在光滑函数加权空间上的满射性

作者: Andreas Debrouwere, Quinten Van Boxstael, Jasson Vindas

类别: math.AP, math.CV, math.FA, Primary. 35A01, 46E10. Secondary. 35A05, 35A35, 35J05

发布日期: 2025-01-07

arXiv ID: 2501.03751v2

链接: http://arxiv.org/abs/2501.03751v2

摘要: We study the surjectivity of the Cauchy-Riemann and Laplace operators on certain weighted spaces of smooth functions of rapid decay on strip-like domains in the complex plane that are defined via weight function systems. We fully characterize when these operators are surjective on such function spaces in terms of a growth condition on the defining weight function systems.

21. Robust Synchronization of Time-Fractional Memristive Hopfield Neural Networks

时滞分数阶忆阻Hopfield神经网络的鲁棒同步

作者: Yuncheng You

类别: math.AP, 26A33, 34A08, 45J05, 68T07, 92B20

发布日期: 2025-10-23

arXiv ID: 2510.20949v1

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

摘要: In this paper we study robust synchronization of time-fractional Hopfield neural networks with memristive couplings and Hebbian learning rules. This new model of artificial neural networks exhibits strong memory and long-range path-dependence in learning processes. Through scaled group estimates it is proved that under rather general assumptions the solution dynamics is globally dissipative. The main result established a threshold condition for achieving robust synchronization of the neural networks if it is satisfied by the interneuron coupling strength coefficient. The synchronizing threshold is explicitly computable in terms of the original parameters and strictly decreasing for the fractional order $\alpha \in (0, 1)$.

22. The linearized Floer equation in a chart

图表中的线性化Floer方程

作者: Urs Frauenfelder, Joa Weber

类别: math.SG, math.AP, math.FA, 57R58 47A53 46B70

发布日期: 2025-10-23

arXiv ID: 2510.20945v1

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

摘要: In this article we are considering the Hessian of the area functional in a non-Darboux chart. This does not seem to have been considered before and leads to an interesting new mathematical structure which we introduce in this article and refer to as \emph{almost extendable weak Hessian~field}. Our main result is a Fredholm theorem for Robbin-Salamon operators associated to \emph{non-continuous} Hessians which we prove by taking advantage of this new structure.