生成时间: 2025年10月15日 07:28:23

重点关注论文数量: 4

了解领域论文数量: 16

不相关论文数量: 15

重点关注论文(完整分析)

1. Nonlinear instability for the 3D MHD equations around the Taylor-Couette flow

三维磁流体动力学方程在泰勒-库埃特流附近的非线性不稳定性

作者: Víctor Navarro-Fernández, David Villringer

类别: math.AP

发布日期: 2025-10-14

arXiv ID: 2510.12707v1

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

详细分析

中文标题: 三维磁流体动力学方程在泰勒-库埃特流附近的非线性不稳定性**


1. 研究对象和背景

本文研究的是三维不可压缩磁流体动力学(MHD)方程组在环形圆柱域 $\Omega \subset \mathbb{R}^3$ 上围绕一个显式稳态解 $(\mathbf{u}_{TC}, \mathbf{0})$ 的非线性不稳定性问题。该方程组描述了导电流体中速度场 $\mathbf{u}(t,\mathbf{x})$ 与磁场 $\mathbf{B}(t,\mathbf{x})$ 的耦合演化,其数学形式如下:

$$ \begin{aligned} \partial_t \mathbf{u} + (\mathbf{u} \cdot \nabla) \mathbf{u} + \nabla p &= \nu \Delta \mathbf{u} + (\mathbf{B} \cdot \nabla) \mathbf{B}, \\ \partial_t \mathbf{B} + (\mathbf{u} \cdot \nabla) \mathbf{B} - (\mathbf{B} \cdot \nabla) \mathbf{u} &= \varepsilon \Delta \mathbf{B}, \\ \nabla \cdot \mathbf{u} &= 0, \\ \nabla \cdot \mathbf{B} &= 0, \end{aligned} $$

其中 $(t,\mathbf{x}) \in (0,\infty) \times \Omega$,$p(t,\mathbf{x})$ 为流体压力,$\nu > 0$ 为运动粘性系数,$\varepsilon > 0$ 为磁扩散率(电阻率)。该系统在物理上描述了导电流体(如等离子体或液态金属)在磁场作用下的动力学行为,广泛应用于天体物理、核聚变装置和地球发电机理论中。

研究区域 $\Omega$ 是一个环形圆柱体,用柱坐标 $(r,\theta,z)$ 表示为:

$$ \Omega = \{(r,\theta,z) \in (0,\infty) \times \mathbb{T} \times \mathbb{T} : R_1 \leq r \leq R_2\}, $$

其中 $0 < R_1 < R_2$ 分别为内、外半径,$\mathbb{T} = \mathbb{R}/2\pi\mathbb{Z}$ 表示周期性角向与轴向方向。该几何结构在实验上对应于经典的泰勒-库埃特(Taylor-Couette)装置,其边界条件为物理刚性壁:流体无法穿透边界,且内外圆柱分别以恒定角速度旋转。具体地,速度场满足狄利克雷边界条件:

$$ \mathbf{u} \cdot \hat{\mathbf{n}}|_{\partial\Omega} = 0, \quad u_\theta|_{r=R_1} = 0, \quad u_z|_{r=R_1} = 0, \quad u_\theta|_{r=R_2} = \beta_1, \quad u_z|_{r=R_2} = \beta_2, $$

其中 $\hat{\mathbf{n}}$ 为边界外法向单位向量,$\beta_1, \beta_2 \in \mathbb{R}\setminus\{0\}$ 为常数。磁场则满足完美导体边界条件:

$$ \mathbf{B} \cdot \hat{\mathbf{n}}|_{\partial\Omega} = 0, \quad \hat{\mathbf{n}} \times (\nabla \times \mathbf{B})|_{\partial\Omega} = 0, $$

这等价于磁场切向分量连续、法向分量为零,是等离子体物理和磁流体实验中的标准条件。

稳态解 $(\mathbf{u}_{TC}, \mathbf{0})$ 是由泰勒-库埃特速度场与零磁场构成的显式解。其速度场在柱坐标下为:

$$ \mathbf{u}_{TC}(r) = \left( \frac{a_1}{r} + a_2 r \right) \hat{\boldsymbol{\theta}} + (a_3 \log r + a_4) \hat{\mathbf{z}}, $$

其中系数 $a_1, a_2, a_3, a_4 \in \mathbb{R}$ 由边界条件 (1.3) 确定,且 $a_1, a_2 \neq 0$。该流场是无外力下纳维-斯托克斯方程的精确稳态解,具有径向剪切与轴向流动,是流体动力学中著名的旋转流结构。其稳定性研究已有大量工作,但将其与磁流体耦合后,特别是当磁场初始为零时,其非线性稳定性尚无严格数学结果。

本文关注的核心物理现象是“发电机效应”(dynamo effect):即流体运动通过感应机制将动能转化为磁能,导致磁场指数增长。在经典物理文献中,该现象被广泛猜想为磁星、地球磁场等的起源机制,但数学上缺乏在无外力、三维、无强制项情形下的严格证明。本文首次在无外力、三维、完全无强制的MHD系统中,证明了当速度场为泰勒-库埃特流、磁场初始为零时,即使初始扰动任意小,磁场也会因线性不稳定模的放大而指数增长,从而导致整个系统在 $L^p$ 空间中非线性不稳定。该研究填补了MHD稳定性理论中的关键空白,为发电机机制提供了首个严格的数学验证。


2. 主要定理或主要结果

本文的核心结果是定理 1.1,其内容如下:

定理 1.1(非线性不稳定性):设 $\mathbf{u}_{TC}$ 为泰勒-库埃特速度场(公式 (1.5))。对任意 $s \geq 0$、$\nu > 0$ 和足够小的磁扩散率 $\varepsilon > 0$,存在常数 $\chi > 0$、$c > 0$,使得对任意小的扰动尺度 $\delta > 0$,均存在初始扰动 $(\mathbf{u}_0, \mathbf{B}_0) \in W^{s,p}(\Omega)$,满足:

$$ \left\| \begin{pmatrix} \mathbf{u}_0 \\ \mathbf{B}_0 \end{pmatrix} - \begin{pmatrix} \mathbf{u}_{TC} \\ \mathbf{0} \end{pmatrix} \right\|_{W^{s,p}} \leq \delta, $$

且存在时间 $t^\star = c |\log \delta|$,使得对应的MHD方程解 $(\mathbf{u}(t^\star), \mathbf{B}(t^\star))$ 满足:

$$ \left\| \begin{pmatrix} \mathbf{u}(t^\star) \\ \mathbf{B}(t^\star) \end{pmatrix} - \begin{pmatrix} \mathbf{u}_{TC} \\ \mathbf{0} \end{pmatrix} \right\|_{L^p} \geq \chi. $$

此外,若粘性 $\nu$ 足够大,满足 $\nu \gg \|\mathbf{u}_{TC}\|_{W^{1,\infty}}$,则有 $\|\mathbf{B}(t^\star)\|_{L^p} \geq \chi$。

该定理的假设条件非常自然:仅要求初始扰动在任意高阶 Sobolev 空间 $W^{s,p}$ 中足够小($s \geq 0$,$p \in (1,\infty)$),磁扩散率 $\varepsilon$ 小于某个阈值(但非趋于零),粘性 $\nu$ 为任意正数。结论则表明:即使初始扰动在任意强的函数空间中任意小,系统仍会在对数时间尺度 $t^\star \sim |\log \delta|$ 内偏离稳态达到 $\mathcal{O}(1)$ 的距离,这正是“非线性不稳定”的标准定义。

特别地,定理的第二部分明确指出:当粘性 $\nu$ 足够大时,速度场的线性稳定性被保留(即 $\mathbf{u}$ 部分不增长),而磁场 $\mathbf{B}$ 的增长主导了整个系统的不稳定性。这意味着不稳定性来源于磁场的“自主”指数增长,而非速度场的失稳,从而首次在无外力三维MHD系统中证实了“能量从速度场向磁场转移”的物理猜想。

该结果的突破性在于其“非平凡性”:此前的MHD不稳定性结果要么依赖外力驱动(如 Gérard-Varet & Rousset [11,12]),要么研究的是稳定流(如Couette流在二维情形),要么仅限于线性分析。而本文首次在无外力、三维、无强制、且稳态为非均匀剪切流的情形下,证明了非线性不稳定,并给出了明确的时间尺度与增长幅度。

此外,定理所依赖的线性不稳定模的存在性由命题 2.1 保证:对任意足够小的 $\varepsilon > 0$,存在一个 $C^2$ 拓扑无散磁场 $\mathbf{B}_0$ 和正实数 $\lambda \sim \varepsilon^{1/3}$,使得 $\mathbf{B}(t) = e^{\lambda t} \mathbf{B}_0$ 是线性化磁流体方程(即“运动学发电机方程”)的解:

$$ \partial_t \mathbf{B} + (\mathbf{u}_{TC} \cdot \nabla) \mathbf{B} - (\mathbf{B} \cdot \nabla) \mathbf{u}_{TC} = \varepsilon \Delta \mathbf{B}. $$

该模是线性算子 $A_2$ 的特征函数,其特征值 $\lambda$ 的实部为正,且随 $\varepsilon \to 0$ 以 $\varepsilon^{1/3}$ 速度增长,这与实验和数值模拟中观察到的“弱扩散下更强不稳定性”现象一致。


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

本文的证明框架建立在 Friedlander–Pavlović–Shvydkoy [10] 的非线性不稳定方法之上,并将其推广至包含磁场的完整MHD系统。其核心思想是:若线性化算子存在正实部的特征值(即线性不稳定),则非线性系统在小扰动下会在有限时间内“放大”该不稳定模,从而导致解偏离稳态。

a) 核心分析工具

  1. 线性算子的谱分析与生成解析半群:作者将MHD系统线性化后,得到算子 $A = A_1 \oplus A_2$,其中 $A_1$ 为线性化Navier-Stokes算子,$A_2$ 为运动学发电机算子。关键在于证明 $A_2$ 在 $L^p$ 空间中生成解析半群,且具有正谱半径。这依赖于对Hodge-Laplacian $\Delta = -\nabla \times \nabla \times$ 的深入分析,结合完美导体边界条件下的椭圆正则性估计(引理3.2)。

  2. 分数阶算子与插值理论:为控制非线性项在 $L^p$ 中的估计,作者引入了分数阶算子 $A_\eta^\alpha$ 和 $A_\eta^{-\alpha}$,其中 $A_\eta = A - \lambda - \eta$。通过引理3.1,他们建立了关键的估计:

    $$ \|A_\eta^\alpha e^{A_\eta t}\|_{L^p \to L^p} \lesssim t^{-\alpha}, \quad \text{对 } \alpha \in (0,1), $$

    $$ \|A_\eta^{-\alpha} \nabla \mathbf{H}\|_{L^q} \lesssim \|\mathbf{H}\|_{L^p}, \quad \text{对 } \alpha \in (1/2,1), \ q \leq \frac{3p}{3-(2\alpha-1)p}. $$

    这些估计依赖于Hodge-Laplacian的插值空间理论(引理3.2)和Sobolev嵌入,使得非线性项 $N(\mathbf{v},\mathbf{v})$, $N(\mathbf{B},\mathbf{B})$, $M(\mathbf{v},\mathbf{B})$ 能被有效控制。

  3. Mild形式与积分估计:作者采用“弱解”或“mild解”形式:

    $$ \mathbf{w}(t) = e^{At} \mathbf{w}_0 + \int_0^t e^{A(t-\tau)} \mathbf{F}(\mathbf{w}(\tau)) \, d\tau, $$

    其中 $\mathbf{F}$ 为非线性项。通过将 $e^{A(t-\tau)} = e^{(\lambda+\eta)(t-\tau)} A_\eta^\alpha e^{A_\eta(t-\tau)} A_\eta^{-\alpha}$ 分解,结合上述分数阶估计,将非线性项的积分控制为 $\lesssim \delta^2 e^{2\lambda t}$,远小于线性项的 $\delta e^{\lambda t}$。

  4. 反证法与时间阈值构造:定义时间 $T = \inf\{t \geq 0 : \|\mathbf{w}(t)\|_{L^q} \geq e^{\lambda t} \delta Q\}$,利用连续性论证,证明若 $T = \infty$,则会导致矛盾。从而证明存在有限时间 $t^\star < T$,使得解脱离稳态,完成非线性不稳定证明。

b) 关键技术与难题突破

  1. 完美导体边界下Hodge Laplacian的处理:与标准Dirichlet Laplacian不同,$A_2 = -(\mathbf{u}_{TC} \cdot \nabla)\mathbf{B} + (\mathbf{B} \cdot \nabla)\mathbf{u}_{TC} + \varepsilon \Delta$ 是一个非自伴、非椭圆的扰动算子。作者通过引理3.3证明其生成解析半群,关键技巧是将 $\nabla \times (\mathbf{u}_{TC} \times \cdot)$ 视为 $\Delta$ 的低阶扰动,并利用 $\|\nabla \times \mathbf{B}\|_{L^2} \sim \|\mathbf{B}\|_{\dot{H}^1}$(由无散条件保证)建立先验估计。

  2. 分数阶 Sobolev 空间与嵌入的精确控制:为处理非线性项,作者必须在 $L^q$ 与 $L^p$ 之间建立 $\nabla$ 与 $A_\eta^{-\alpha}$ 的复合估计。他们利用引理3.4(算子扰动下的插值空间等价性)将 $A_2$ 的分数幂空间与Hodge Laplacian的分数幂空间等价化,从而能直接引用经典椭圆理论中的嵌入结果。

  3. 从 $L^q$ 到 $L^p$ 的推广:由于引理3.1的估计仅对 $q > 3$ 成立,而目标为任意 $p > 1$,作者巧妙利用域 $\Omega$ 有界,通过 $L^q \hookrightarrow L^p$($q > p$)进行嵌入,从而将结论推广至任意 $p \in (1,\infty)$,这是对 [10] 方法的重要推广。


4. 与之前工作的比较

本文在多个维度上实现了对现有文献的实质性突破:

a) 首个无外力三维MHD非线性不稳定性结果:此前所有三维MHD非线性不稳定结果均依赖外力驱动。例如,Gérard-Varet 和 Rousset [11,12] 虽研究了类似发电机机制,但其速度场由外部强迫维持,非自然稳态。本文首次在无外力、无强制、稳态为自然解(即泰勒-库埃特流)的情形下,证明了非线性不稳定,这更贴近物理现实。

b) 放宽初始数据正则性要求:文献中二维Couette流的稳定性结果(如 [7,17,20,35])要求初始扰动在 $H^s$ 中($s>0$),而本文允许任意 $s \geq 0$,甚至 $W^{s,p}$ 中的任意低正则性扰动,只要足够小即可触发不稳定性,说明该不稳定性具有“鲁棒性”。

c) 维度与几何结构的推广:此前MHD不稳定性研究多集中于二维平面流(如 $T \times \mathbb{R}$)或周期性矩形域。本文首次在三维环形域中处理具有非均匀剪切($u_\theta \sim 1/r$)和轴向流动的复杂流场,克服了柱坐标下非均匀系数带来的分析困难。

d) 从线性到非线性的严格升级:虽然 [25] 已证明线性不稳定模的存在,但从未将其升级为非线性不稳定性。本文通过系统整合 [10] 的框架与 [25] 的谱结果,完成了这一关键跃迁,使数学物理中长期猜想的“能量转移”机制获得首个严格证明。

e) 明确的时间尺度与增长速率:本文不仅证明不稳定,还给出了精确的时间尺度 $t^\star \sim |\log \delta|$ 和增长速率 $\lambda \sim \varepsilon^{1/3}$,与数值模拟中观察到的“弱电阻率下更强不稳定性”完全一致,为实验设计提供了理论依据。

综上,本文不仅在技术上完成了一次从线性到非线性、从二维到三维、从强制到无强制的全面升级,更在物理意义上首次为MHD发电机机制提供了坚实的数学基础,是该领域里程碑式的工作。


2. Two-temperature fluid models for a polyatomic gas based on kinetic theory for nearly resonant collisions

基于近共振碰撞动力学理论的多原子气体双温度流体模型

作者: Kazuo Aoki, Niclas Bernhoff

类别: math.AP, math-ph, math.MP, 76P05 (Primary) 35Q20, 82C40 (Secondary)

发布日期: 2025-10-14

arXiv ID: 2510.12464v1

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

详细分析

**)δ/2−1 (II)*δ/2−1

中文标题: 基于近共振碰撞动力学理论的多原子气体双温度流体模型


1. 研究对象和背景

本文研究的是具有弱平动-内能模式耦合的多原子理想气体的非平衡流体动力学行为。这类气体在高温、高速流动(如高超声速飞行、稀薄气体动力学、等离子体物理)中普遍存在,其分子内部存在振动、转动等自由度,导致能量在平动动能与内能之间缓慢交换。传统单温度流体模型(如欧拉或纳维-斯托克斯方程)假设系统处于局部热力学平衡,即平动温度 $ T_{\text{tr}} $ 与内能温度 $ T_{\text{int}} $ 相等,这在强耦合系统中近似成立,但在弱耦合条件下失效。此时,两种温度可长期保持差异,系统呈现“双温度”特性,需采用更精细的多温度流体模型。

为系统性地推导此类模型,作者基于玻尔兹曼方程构建了一个新的碰撞算子模型。该模型的核心是将总碰撞过程分解为两类:非弹性碰撞(inelastic collisions)与共振碰撞(resonant or elastic collisions)。其中,非弹性碰撞允许平动与内能模式间能量交换,而共振碰撞则严格禁止能量转移,即满足:

$$ \frac{m}{2}|\xi|^2 + \frac{m}{2}|\xi^*|^2 = \frac{m}{2}|\xi'|^2 + \frac{m}{2}|\xi^{*'}|^2, \quad I + I^* = I' + I^{*'}, $$

其中 $ \xi, \xi^* $ 为碰撞前分子速度,$ \xi', \xi^{*'} $ 为碰撞后速度,$ I, I^* $ 为内能,$ m $ 为分子质量。该条件等价于相对速度模长守恒:$ |g| = |g'| $,其中 $ g = \xi - \xi^* $,$ g' = \xi' - \xi^{*'} $,且内能变化 $ \Delta I = I' + I^{*'} - I - I^* = 0 $。

作者提出的玻尔兹曼方程形式为:

$$ \frac{\partial f}{\partial t} + \xi \cdot \nabla_x f = Q_\theta(f, f), $$

其中 $ f(t, x, \xi, I) $ 是七维相空间中的速度-内能分布函数,$ x \in \mathbb{R}^3 $ 为空间位置,$ \xi \in \mathbb{R}^3 $ 为分子速度,$ I \in \mathbb{R}_+ $ 为单个分子的内能。碰撞算子 $ Q_\theta $ 是一个双线性算子,其结构为:

$$ Q_\theta(f, f) = \theta Q_s(f, f) + (1 - \theta) Q_r(f, f), $$

其中 $ \theta \in [0,1] $ 是控制能量交换强度的参数:当 $ \theta = 1 $ 时,系统为标准非弹性碰撞($ Q_s $),能量完全交换;当 $ \theta = 0 $ 时,系统为纯共振碰撞($ Q_r $),能量完全不交换。该模型继承了Borgnakke-Larsen随机能量分配机制,但通过参数 $ \theta $ 实现了对碰撞机制的连续调控,从而为弱耦合极限提供了数学框架。

宏观量由分布函数的矩定义,包括数密度 $ n $、质量密度 $ \rho = mn $、流速 $ \mathbf{u} $、平动温度 $ T_{\text{tr}} $、内能温度 $ T_{\text{int}} $、应力张量 $ \mathbf{p} $ 和热流矢量 $ \mathbf{q} $。特别地,内能 $ e $ 被分解为平动与内能贡献:

$$ e = e_{\text{tr}} + e_{\text{int}}, \quad e_{\text{tr}} = \frac{1}{n} \left( \frac{m}{2} |\xi - \mathbf{u}|^2, f \right), \quad e_{\text{int}} = \frac{1}{n} (I, f), $$

相应的温度定义为:

$$ T_{\text{tr}} = \frac{2}{3k_B} e_{\text{tr}}, \quad T_{\text{int}} = \frac{1}{\delta k_B} e_{\text{int}}, $$

其中 $ \delta \geq 2 $ 为内自由度(可非整数),$ k_B $ 为玻尔兹曼常数。这一双温度结构是本文模型的核心特征。

研究区域为整个空间 $ \mathbb{R}^3 $,无边界条件限制(因聚焦于宏观方程推导,而非边界层问题),函数空间为 $ L^2(\mathbb{R}^3 \times \mathbb{R}_+) $,内积为 $ (f,g) = \iint f g \, d\xi dI $。作者假设气体为稀薄理想气体,忽略多体碰撞,仅考虑二元碰撞。系统在无外力作用下演化,初始分布为接近局部平衡态的扰动,即 $ f = M(1 + h) $,其中 $ M $ 为对应于 $ \theta $ 的局部平衡分布。

该模型的物理背景是:在高温多原子气体(如氮气、氧气、二氧化碳)中,分子内能激发(如振动模式)的弛豫时间远大于分子碰撞时间,导致 $ T_{\text{tr}} \neq T_{\text{int}} $。例如,在高超声速再入飞行器周围,气体被剧烈压缩加热,但内能模式尚未充分激发,形成“温度非平衡区”。本文旨在从动力学层面系统推导描述此类现象的宏观流体方程,填补了从玻尔兹曼方程直接推导双温度NS模型的理论空白。


2. 主要定理或主要结果

本文的核心结果是通过Chapman-Enskog展开,从所提出的双组分碰撞玻尔兹曼方程中,系统性地推导出两类双温度流体方程,分别对应于 $ \theta $ 与Knudsen数 $ \text{Kn} $ 的不同量级关系。

主要定理 1($ \theta = O(\text{Kn}^2) $ 情形)
假设:Knudsen数 $ \text{Kn} \ll 1 $,且 $ \theta = \varepsilon^2 $,其中 $ \varepsilon = \text{Kn} $。分布函数 $ f $ 具有渐近展开:

$$ f = f_0 + \varepsilon f_1 + \varepsilon^2 f_2 + \cdots, $$

其中 $ f_0 = M $ 为局部平衡态(形式见式(28),即共振情形的分布),且满足 $ \theta Q_s(f_0, f_0) = 0 $(因 $ f_0 $ 是 $ Q_r $ 的平衡态,而非 $ Q_s $ 的)。

结论:在 $ \varepsilon \to 0 $ 极限下,得到如下双温度欧拉型方程组:

$$ \begin{aligned} &\frac{\partial n}{\partial t} + \nabla \cdot (n \mathbf{u}) = 0, \\ &\frac{\partial (n \mathbf{u})}{\partial t} + \nabla \cdot (n \mathbf{u} \otimes \mathbf{u} + p \mathbf{I}) = 0, \\ &\frac{\partial}{\partial t} \left( \frac{3}{2} n k_B T_{\text{tr}} + \frac{\delta}{2} n k_B T_{\text{int}} \right) + \nabla \cdot \left( \frac{5}{2} n k_B T_{\text{tr}} \mathbf{u} + \frac{\delta + 2}{2} n k_B T_{\text{int}} \mathbf{u} \right) = 0, \\ &\frac{\partial}{\partial t} \left( \frac{3}{2} n k_B T_{\text{tr}} \right) + \nabla \cdot \left( \frac{5}{2} n k_B T_{\text{tr}} \mathbf{u} \right) = - \frac{1}{\tau_{\text{tr}}} (T_{\text{tr}} - T_{\text{int}}), \\ &\frac{\partial}{\partial t} \left( \frac{\delta}{2} n k_B T_{\text{int}} \right) + \nabla \cdot \left( \frac{\delta + 2}{2} n k_B T_{\text{int}} \mathbf{u} \right) = + \frac{1}{\tau_{\text{tr}}} (T_{\text{tr}} - T_{\text{int}}), \end{aligned} $$

其中 $ p = n k_B (T_{\text{tr}} + T_{\text{int}} \cdot \frac{\delta}{3}) $,$ \tau_{\text{tr}} $ 为平动-内能能量弛豫时间,由碰撞核 $ Q_s $ 的谱性质决定,且 $ \tau_{\text{tr}} \sim \theta^{-1} \sim \text{Kn}^{-2} $。该方程组为无粘性、无热传导的双温度欧拉方程,但包含能量交换项(即最后两式中的松弛项),表明即使在零阶近似下,温度差异仍通过碰撞驱动演化。

主要定理 2($ \theta = O(\text{Kn}) $ 情形)
假设:$ \theta = \varepsilon $,$ \varepsilon = \text{Kn} \ll 1 $,且 $ f = f_0 + \varepsilon f_1 + \varepsilon^2 f_2 + \cdots $,$ f_0 = M_r $(共振平衡态)。

结论:在 $ \varepsilon \to 0 $ 极限下,得到双温度Navier-Stokes型方程组,包含粘性应力、热传导与能量弛豫:

$$ \begin{aligned} &\frac{\partial n}{\partial t} + \nabla \cdot (n \mathbf{u}) = 0, \\ &\frac{\partial (n \mathbf{u})}{\partial t} + \nabla \cdot (n \mathbf{u} \otimes \mathbf{u} + \mathbf{p}) = 0, \\ &\mathbf{p} = -p \mathbf{I} + \mathbf{\tau}, \quad \mathbf{\tau} = \mu \left( \nabla \mathbf{u} + (\nabla \mathbf{u})^\top - \frac{2}{3} (\nabla \cdot \mathbf{u}) \mathbf{I} \right), \\ &\frac{\partial}{\partial t} \left( \frac{3}{2} n k_B T_{\text{tr}} + \frac{\delta}{2} n k_B T_{\text{int}} \right) + \nabla \cdot \left( \frac{5}{2} n k_B T_{\text{tr}} \mathbf{u} + \frac{\delta + 2}{2} n k_B T_{\text{int}} \mathbf{u} + \mathbf{q}_{\text{tr}} + \mathbf{q}_{\text{int}} \right) = 0, \\ &\frac{\partial}{\partial t} \left( \frac{3}{2} n k_B T_{\text{tr}} \right) + \nabla \cdot \left( \frac{5}{2} n k_B T_{\text{tr}} \mathbf{u} + \mathbf{q}_{\text{tr}} \right) = - \frac{1}{\tau_{\text{tr}}} (T_{\text{tr}} - T_{\text{int}}), \\ &\frac{\partial}{\partial t} \left( \frac{\delta}{2} n k_B T_{\text{int}} \right) + \nabla \cdot \left( \frac{\delta + 2}{2} n k_B T_{\text{int}} \mathbf{u} + \mathbf{q}_{\text{int}} \right) = + \frac{1}{\tau_{\text{tr}}} (T_{\text{tr}} - T_{\text{int}}), \end{aligned} $$

其中热流由傅里叶定律给出:

$$ \mathbf{q}_{\text{tr}} = - \kappa_{\text{tr}} \nabla T_{\text{tr}}, \quad \mathbf{q}_{\text{int}} = - \kappa_{\text{int}} \nabla T_{\text{int}}, $$

粘性系数 $ \mu \sim \text{Kn}^{-1} \theta^{-1} \sim \text{Kn}^{-2} $,热导率 $ \kappa_{\text{tr}}, \kappa_{\text{int}} $ 也依赖于 $ \theta $ 与 $ \text{Kn} $。松弛时间 $ \tau_{\text{tr}} \sim \theta^{-1} \sim \text{Kn}^{-1} $,表明温度弛豫与粘性效应同阶,是主导的非平衡效应。

这两个定理表明:当碰撞中共振成分占优时,系统的宏观行为由温度弛豫主导,且其强度直接由 $ \theta $ 与 $ \text{Kn} $ 的相对大小决定。当 $ \theta \ll \text{Kn} $ 时,弛豫效应在NS方程中体现为一阶项;当 $ \theta \ll \text{Kn}^2 $ 时,弛豫仅在欧拉方程中作为源项出现,而粘性与热传导仍为零阶。这一区分是本文的理论突破。

此外,作者证明了碰撞算子 $ Q_\theta $ 的平衡态唯一性(Proposition 4)、H定理(Remark 1)及线性化算子 $ L_\theta $ 的Fredholm性质(Proposition 3, Lemma 2),为Chapman-Enskog展开的数学严谨性提供了基础。


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

本文的核心方法是Chapman-Enskog展开(Chapman–Enskog expansion),这是一种从玻尔兹曼方程推导宏观流体方程的标准渐近方法,通过将分布函数展开为Knudsen数的幂级数,并逐阶求解积分方程。

a) 核心分析学工具:

  1. 线性化碰撞算子与Fredholm理论:作者将玻尔兹曼方程线性化为 $ f = M(1 + h) $,得到线性算子 $ L_\theta h = \nu_\theta h - K_\theta(h) $。为求解高阶项(如 $ f_1, f_2 $),需求解形如 $ L_\theta h = \phi $ 的方程。为此,作者证明了 $ L_\theta $ 的Fredholm性质(Proposition 3, Lemma 2),即其核空间由碰撞不变量张成(Proposition 2),且值域为核的正交补。这保证了在相容性条件满足时,方程有唯一解(模去核方向),这是Chapman-Enskog展开可行性的数学基石。

  2. 不变测度与对称性分析:作者引入了碰撞不变测度 $ dA_\theta $(式(23))与线性化后的测度 $ \widetilde{dA}_\theta $(Lemma 2),并证明其在碰撞变量交换下不变(式(24))。这一对称性被用于将双线性形式 $ (Q_\theta(f,F), g) $ 转化为对称形式(Proposition 1),从而简化了能量矩计算与不变量识别。

  3. 广义Dirac函数与积分变换:为处理共振碰撞中 $ \delta_1(\Delta I) $ 的奇异性,作者借助Dirac函数的积分性质,将碰撞核 $ W_r $ 从原始变量 $ (\xi, \xi^*, I, I^*) $ 转换为Borgnakke-Larsen变量 $ (R, r, \sigma) $(式(14b), (21)),从而将积分转化为可计算的显式形式。这一技巧避免了直接处理高维奇异积分的困难。

  4. 矩方法与平衡态投影:在Chapman-Enskog展开中,每阶方程需投影到碰撞不变量空间(质量、动量、能量)上。作者利用平衡态分布 $ M_s $ 与 $ M_r $ 的显式形式(式(27)-(28))精确计算了所有矩(如 $ (|\xi|^2, M_r) $, $ (I, M_r) $),从而直接导出温度定义与能量守恒方程。

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

  1. 双参数耦合的尺度分离:本文的关键挑战是同时处理两个小参数:Knudsen数 $ \text{Kn} $ 与能量交换参数 $ \theta $。作者的突破在于区分两种极限情形:$ \theta = O(\text{Kn}^2) $ 与 $ \theta = O(\text{Kn}) $。在第一种情形下,$ \theta \ll \text{Kn} $,碰撞算子中非弹性部分为高阶小量,仅在能量方程中作为源项出现;而在第二种情形下,$ \theta \sim \text{Kn} $,非弹性碰撞与粘性效应同阶,需在NS方程中同时保留粘性与弛豫项。这种分层处理避免了传统模型中“一揽子”近似导致的物理失真。

  2. 非平衡态平衡分布的构造:传统模型(如BGK)使用单温度Maxwell分布,但本文需构造双温度平衡态 $ M_r $(式(28)),其形式为:

$$ M_r = \frac{n I^{\delta/2 - 1}}{(2\pi k_B T_{\text{tr}}/m)^{3/2} (k_B T_{\text{int}})^{\delta/2} \Gamma(\delta/2)} \exp\left( -\frac{m|\xi - \mathbf{u}|^2}{2k_B T_{\text{tr}}} - \frac{I}{k_B T_{\text{int}}} \right), $$

该分布是两个独立Maxwell-Boltzmann分布的乘积,反映了平动与内能模式的解耦。作者证明其为 $ Q_r $ 的唯一平衡态(Proposition 4),并验证其满足H定理,为展开提供稳定基底。

  1. 碰撞核的显式分解:通过将 $ \sigma_\theta = \theta \sigma_s + (1-\theta) \sigma_r $(式(12))代入,作者将总碰撞算子分解为标准与共振部分(式(13)),并分别推导其Borgnakke-Larsen形式(式(18),(21))。这一分解使Chapman-Enskog展开中每项的物理来源清晰可辨,避免了传统模型中“黑箱式”碰撞项的模糊性。

4. 与之前工作的比较

本文在多个维度上显著超越了现有文献。

a) 理论突破:从松弛模型到玻尔兹曼型模型

此前的双温度模型主要基于弛豫型模型(如ES模型、BGK模型),如文献[2]基于Ellipsoidal Statistical模型推导出双温度NS方程。这类模型的优势是显式、易数值实现,但其本质是人为构造的简化算子,缺乏从真实玻尔兹曼方程推导的物理根基。本文首次在真实玻尔兹曼方程框架下,通过构造具有明确物理含义的碰撞核 $ Q_\theta $,系统推导出双温度方程。这使得模型参数(如弛豫时间 $ \tau_{\text{tr}} $、热导率)具有明确的分子碰撞起源,而非经验拟合。

b) 放宽假设条件:从强耦合到弱耦合极限

现有研究(如[57,27,28])通常假设碰撞为完全非弹性($ \theta \sim 1 $),或仅处理温度趋于平衡的快速弛豫情形。本文首次系统研究了弱耦合极限($ \theta \to 0 $),即内能弛豫缓慢、$ T_{\text{tr}} \neq T_{\text{int}} $ 可长期维持的情形。这在物理上对应于高能分子(如N₂的振动激发)与平动模式间耦合极弱的场景(如高层大气、电离层)。此前模型无法描述此类现象,因它们隐含假设 $ T_{\text{tr}} \approx T_{\text{int}} $。

c) 结果强度:首次获得双温度NS方程的完整形式

本文不仅导出了双温度欧拉方程,更重要的是在 $ \theta = O(\text{Kn}) $ 情形下,首次从玻尔兹曼方程推导出完整的双温度Navier-Stokes方程组,包括:

  • 粘性应力张量 $ \mathbf{\tau} $,
  • 两个独立热流 $ \mathbf{q}_{\text{tr}}, \mathbf{q}_{\text{int}} $,
  • 两个耦合的松弛源项。

这些项均具有明确的解析表达式(依赖于 $ \delta, \theta, \text{Kn} $),而此前基于ES模型的成果(如[2])虽有显式参数依赖,但未考虑热传导与粘性在双温度下的耦合结构。本文的方程组在数学上是自洽、完备、物理可解释的,是该领域里程碑式的结果。

d) 与高阶矩模型的比较

文献[5,51,52,56]等基于矩方法或不可逆热力学构造了更高阶的双温度方程(如13或21矩模型),但它们依赖于人工闭合假设(如假设高阶矩与低阶矩的代数关系),缺乏动力学基础。本文的模型则直接从玻尔兹曼方程逐阶展开,闭合由碰撞算子的谱结构自然决定,因此具有更强的物理一致性与数学严谨性。

综上,本文通过引入“近共振”碰撞模型与双参数渐近分析,实现了从微观碰撞动力学到宏观流体方程的首次完整、系统、物理清晰的桥梁构建,为高温多原子气体非平衡流动的精确模拟奠定了坚实的理论基础。


3. Random dynamics and invariant measures for a class of non-Newtonian fluids of differential type on 2D and 3D Poincaré domains

二维与三维Poincaré域上一类微分型非牛顿流体的随机动力学与不变测度

作者: Kush Kinra, Fernanda Cipriano

类别: math.PR, math.AP, 76A05, 35R60, 35Q35, 37L30, 35B40

发布日期: 2025-06-05

arXiv ID: 2506.04801v2

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

详细分析

**
中文标题: 二维与三维Poincaré域上一类微分型非牛顿流体的随机动力学与不变测度


1. 研究对象和背景

本文研究的是一类具有随机扰动的三阶非牛顿流体方程,该方程是经典Navier-Stokes方程的非线性推广,用于描述具有复杂粘弹性行为的流体,如血液、纳米流体、聚合物溶液等。与牛顿流体(剪切应力与剪切速率呈线性关系)不同,非牛顿流体的本构关系包含高阶应变率张量,从而引入更强的非线性项。论文所考虑的系统为:

$$ \begin{cases} dv = \left[ f - \nabla P + \nu \Delta v - (v \cdot \nabla)v + \alpha \operatorname{div}((A(v))^2) + \beta \operatorname{div}(|A(v)|^2 A(v)) \right] dt + dW, & \text{in } \mathcal{O} \times (0,\infty), \\ \operatorname{div} v = 0, & \text{in } \mathcal{O} \times [0,\infty), \\ v(x,0) = v_0(x), & \text{in } \mathcal{O}, \\ v = 0, & \text{on } \partial\mathcal{O} \times [0,\infty), \end{cases} \tag{1.5} $$

其中,$\mathcal{O} \subset \mathbb{R}^d$($d=2$或$3$)为Poincaré域(有界或无界),$v(x,t)$为流体速度场,$P(x,t)$为压力,$f$为外力,$W(t)$为取值于Hilbert空间$H$的双边柱状Wiener过程,代表无穷维随机噪声。$A(v) = \nabla v + (\nabla v)^T$为第一Rivlin-Ericksen张量,表征流体的应变率。参数$\nu > 0$为动力粘度,$\alpha$和$\beta$为材料模量,满足热力学相容性条件:

$$ \nu > 0, \quad \alpha = 0, \quad |\alpha| < \sqrt{2\nu\beta}, \quad \beta > 0. \tag{1.4} $$

该条件排除了二阶流体($\alpha_1 \neq 0$)的不稳定性问题,同时保留了三阶流体中关键的非线性项$\operatorname{div}(|A(v)|^2 A(v))$,该非线性项具有重要的正则化作用,保证了三维情形下解的唯一性。

论文研究的区域$\mathcal{O}$为Poincaré域,即满足Poincaré不等式:

$$ \lambda_0 \int_{\mathcal{O}} |\psi(x)|^2 dx \leq \int_{\mathcal{O}} |\nabla \psi(x)|^2 dx, \quad \forall \psi \in H^1(\mathcal{O}), \tag{1.8} $$

其中$\lambda_0 > 0$为常数。这一条件在有界域上自动成立(由Poincaré不等式标准理论),但在无界域(如$\mathcal{O} = \mathbb{R} \times (-L,L)$或$\mathbb{R}^2 \times (-L,L)$)上需额外假设。Poincaré不等式是控制低阶范数、建立能量估计的关键工具。

函数空间的设定也具有特殊性。作者定义了如下空间:

  • $H = \overline{V}^{L^2(\mathcal{O})} = \{ v \in L^2(\mathcal{O}) : \operatorname{div} v = 0, \, v \cdot n = 0 \text{ on } \partial\mathcal{O} \}$,
  • $V = \overline{V}^{H^1(\mathcal{O})} = \{ v \in H^1(\mathcal{O}) : \operatorname{div} v = 0 \}$,
  • $X = W^{1,4}(\mathcal{O}) \cap V$,配备范数$\| \cdot \|_X = \| \cdot \|_{W^{1,4}} + \| \cdot \|_V$。

特别地,噪声的Cameron-Martin空间$K$需满足Hypothesis 1.1:$K \subset H \cap W^{1,4}(\mathcal{O})$,且存在$\delta \in (0,1/2)$使得$A^{-\delta}: K \to H \cap W^{1,4}(\mathcal{O})$为$\gamma$-radonifying算子。这一条件确保了Ornstein-Uhlenbeck过程(用于将随机方程转化为确定性系统)能取值于$W^{1,4}$空间,从而支持后续的非线性项估计。在有界域中,该条件等价于噪声协方差算子的谱衰减足够快(即特征值$\lambda_j \sim j^{2/d}$,要求$s > d/4$),而在无界域中则需通过核函数的精细构造实现。

边界条件为Dirichlet齐次边界条件$v=0$,这在数学分析中可保证速度场在边界处被“固定”,便于使用Poincaré不等式和紧嵌入。外力项$f$在有界域中属于$H^{-1}(\mathcal{O}) + W^{-1,4/3}(\mathcal{O})$,而在无界域中被限制为$f \in L^2(\mathcal{O})$,这是由于在无界域中缺乏$H^1 \hookrightarrow L^2$的紧性,迫使作者采用更强的正则性假设以维持能量控制。

该研究的物理背景源于纳米流体技术的发展。纳米流体是将纳米颗粒悬浮于基础流体(如水或乙二醇)中形成的胶体,其热导率显著高于基础流体,广泛应用于微电子冷却、生物医学等领域。三阶流体模型能更精确地描述这类复杂流体的非线性粘弹性响应,因此对其随机动力学行为的研究具有重要工程意义。


2. 主要定理或主要结果

本文的核心成果体现在以下四个主要定理中,其证明分别在第3、4、5、6节完成。

定理1.4(主定理):在假设(1.4)和Hypotheses 1.1–1.2成立的条件下,系统(1.5)满足:

(i) 存在唯一弱解:对于任意初始数据$v_0 \in H$和外力$f \in H^{-1}(\mathcal{O})$,系统(1.5)存在唯一的弱解(解析意义下),即存在一个适应于$\mathcal{F}_t$的$H$-值随机过程$v(t)$,满足对任意$\phi \in V$,几乎必然有:

$$ (v(t), \phi) + \int_0^t \left[ \nu (\nabla v(s), \nabla \phi) - ((v(s) \cdot \nabla)v(s), \phi) + \alpha (A(v(s))^2, \nabla \phi) + \beta (|A(v(s))|^2 A(v(s)), \nabla \phi) \right] ds = (v_0, \phi) + (f, \phi) t + (W(t), \phi), $$

且$v \in L^2_{\text{loc}}([0,\infty); V) \cap C([0,\infty); H)$,并满足能量不等式。

(ii) 有界域上的唯一随机吸引子:当$\mathcal{O}$为有界域时,系统(1.5)生成一个随机动力系统$\Psi$,且存在唯一的随机吸引子$\mathcal{A} \subset H$,满足:

  • $\Psi(t,\omega)\mathcal{A}(\omega) = \mathcal{A}(\theta_t \omega)$(不变性);
  • 对任意有界集$B \subset H$,有$\lim_{t \to \infty} \operatorname{dist}_H(\Psi(t,\theta_{-t}\omega)B, \mathcal{A}(\omega)) = 0$(吸引性);
  • $\mathcal{A}(\omega)$为紧集(在$H$中)。

(iii) 无界域上的唯一随机吸引子:当$\mathcal{O}$为无界Poincaré域时,若外力$f \in L^2(\mathcal{O})$,系统(1.5)仍存在唯一的随机吸引子$\mathcal{A}$,尽管$H^1(\mathcal{O}) \not\hookrightarrow L^2(\mathcal{O})$不紧。

(iv) 不变测度的存在性:在有界域和无界域上,系统(1.5)分别存在不变测度$\mu$,满足对任意有界连续函数$\varphi: H \to \mathbb{R}$,有:

$$ \int_H \mathbb{E}[\varphi(v(t; v_0))] d\mu(v_0) = \int_H \varphi(v_0) d\mu(v_0), \quad \forall t \geq 0. $$

这些结果的证明依赖于一系列关键引理。例如,为证明吸引子存在,作者首先构造了随机吸引子存在的抽象框架(基于[10]),需验证两点:(1) 存在随机吸收集;(2) 随机动力系统$\Psi$是渐近紧的(asymptotically compact)。在有界域中,(2)通过紧嵌入$H^1(\mathcal{O}) \hookrightarrow L^2(\mathcal{O})$直接获得;在无界域中,则必须采用一致尾估计(uniform-tail estimates)方法,通过控制解在远离原点区域的能量衰减来恢复紧性。

此外,定理中特别指出:在无界域中,外力必须属于$L^2(\mathcal{O})$,而非$H^{-1}(\mathcal{O})$,这是由于无法证明解关于初值的弱连续性,从而无法使用能量等式方法。这一限制是当前方法的必然代价。

最终,由于随机吸引子$\mathcal{A}$是紧的、不变的随机集,根据Crauel等人的理论([28]),其存在直接推出不变测度的存在性,无需额外构造。


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

本文的分析框架高度依赖于随机动力系统(RDS)理论与非线性偏微分方程的现代分析工具,其技术路线可归纳为以下三部分:

a) 核心分析工具

  1. Ornstein-Uhlenbeck过程与Doss-Sussmann变换:为处理乘性噪声,作者引入一个取值于$X = W^{1,4} \cap V$的Ornstein-Uhlenbeck过程$z(t)$,其满足:
$$ dz + Az dt = dW, \quad z(0) = 0. $$

随后,通过Doss-Sussmann变换:

$$ u(t) = v(t) - z(t), $$

将随机方程(1.5)转化为确定性随机偏微分方程:

$$ \frac{du}{dt} + Au + G(u) = f, \quad u(0) = v_0 - z(0), \tag{3.18} $$

其中$G(u) = \nu A u + B(u+z) + \alpha J(u+z) + \beta K(u+z)$。此变换将随机项“吸收”进新变量$u$的非线性项中,使后续分析可在确定性框架下进行,是处理非线性噪声的关键步骤。

  1. 单调性方法与Minty-Browder技巧:为证明系统(3.18)的全局弱解存在唯一性,作者利用算子$G$的单调性(引理2.3):
$$ \langle G(y_1) - G(y_2), y_1 - y_2 \rangle + \frac{C_{S,3} C_K}{\nu \varepsilon_0} \|A(y_2 + z)\|_4^2 \|y_1 - y_2\|_2^2 \geq \frac{\nu \varepsilon_0}{4} \|A(y_1 - y_2)\|_2^2 + \frac{\beta \varepsilon_0}{4} \|A(y_1 - y_2)\|_4^4 \geq 0. $$

该不等式表明$G$是强单调且次线性增长的,满足Minty-Browder定理的条件,从而通过逼近序列的弱收敛与强收敛(利用紧嵌入和能量估计)证明解的存在唯一性。

  1. 能量估计与一致尾估计:在有界域中,通过标准Gagliardo-Nirenberg不等式和Poincaré不等式,可获得解在$H$和$V$中的一致有界性,从而构造吸收集。在无界域中,由于$H^1 \not\hookrightarrow L^2$紧,作者引入一致尾估计(uniform-tail estimates):对任意$\varepsilon > 0$,存在$R > 0$,使得对任意$t$足够大,有:
$$ \int_{\mathcal{O} \setminus B_R(0)} |u(t,\omega)|^2 dx < \varepsilon, $$

其中$B_R(0)$为半径$R$的球。该估计通过精细控制非线性项(如$\operatorname{div}(|A(v)|^2 A(v))$)和压力项$P$在远场的贡献实现,最终结合紧性引理(Lemma 5.3)证明$\Psi$是渐近紧的。

  1. $\gamma$-radonifying算子与Helmholtz投影:为处理噪声在$W^{1,4}$空间中的正则性,作者使用$\gamma$-radonifying算子理论,确保$A^{-\delta}$将噪声空间$K$嵌入到$W^{1,4}$中。Helmholtz投影$P$用于分解压力项,使其在无界域中仍保持有界性(Lemma 2.2)。

b) 技术性难题的突破技巧

  1. 处理非线性项$\operatorname{div}(|A(v)|^2 A(v))$:该非线性项具有四次增长性,且含二阶导数,常规方法难以控制。作者利用其结构:
$$ \langle K(v), v \rangle = \frac{1}{2} \|A(v)\|_4^4, $$

从而在能量估计中获得$A(v)$的$L^4$控制,结合Korn不等式((2.2))控制$\nabla v$,实现对高阶非线性项的有效控制。

  1. 克服无界域中缺乏紧性:传统方法(如能量等式)因无法证明解关于初值的弱连续性而失效。作者放弃该路径,转而采用Wang提出的一致尾估计,通过在截断区域上构造能量衰减,绕过紧嵌入的缺失,成为处理无界域SPDE的突破性策略。

  2. 压力项的处理:在尾估计中,压力$P$不能通过散度定理直接消去(因边界在无穷远处)。作者通过将压力项写为$P = \Delta^{-1} \operatorname{div}(\text{非线性项})$,并利用Helmholtz投影的有界性,结合Sobolev嵌入,控制其在$L^{4/3}$或$L^2$中的范数。


4. 与之前工作的比较

本文在多个维度上实现了对现有文献的实质性突破:

a) 首次在无界域中研究三阶非牛顿流体的随机吸引子:此前研究(如[38, 59])仅限于有界域或确定性情形。文献[13, 42]虽研究了随机Navier-Stokes和Brinkman-Forchheimer方程在无界域中的吸引子,但其非线性项为二次型(如$(v\cdot\nabla)v$),而本文的$\operatorname{div}(|A(v)|^2 A(v))$为四次非线性,结构更复杂,且缺乏能量等式(因解不满足弱连续性),使得分析难度剧增。本文首次将统一尾估计方法成功应用于此类高阶非线性流体模型,填补了空白。

b) 放宽对噪声的限制:多数RDS研究(如[49, 76])仅处理加性噪声或线性乘性噪声(如$dW = \sigma dW_t$),而本文的噪声$W$为一般柱状Wiener过程,其Cameron-Martin空间$K$只需满足Hypothesis 1.1,允许更广泛的噪声类型(如分数阶噪声、非对称噪声),突破了“仅限于特定结构噪声”的传统限制。

c) 维度与正则性推广:文献[69]研究了有界域中三阶流体的乘性Lipschitz噪声,但未涉及随机吸引子或不变测度。本文不仅将结果推广至三维,还首次在无界域中建立了不变测度的存在性。此外,初始数据仅需$v_0 \in H$($L^2$正则性),而此前工作(如[2, 16])要求$H^2$或$H^3$初值,本文通过单调性方法将正则性要求降至最优水平。

d) 方法论创新:作者放弃了依赖能量等式的传统方法(在无界域中失效),转而采用“一致尾估计+渐近紧性”这一更普适的框架,为研究其他高阶非线性SPDE(如Oldroyd-B模型、微极流体)提供了新范式。此外,通过Doss-Sussmann变换将噪声“确定化”,避免了对噪声系数的对称性假设,是处理一般非线性噪声的里程碑式尝试。

综上,本文不仅是对三阶非牛顿流体随机动力学的首次系统研究,更在方法论层面推动了无界域SPDE吸引子理论的发展,为复杂流体的长期行为建模提供了坚实的数学基础。


4. High-capillarity limit and smoothing effect of large solutions for a multi-dimensional generic non-conservative compressible two-fluid model

高毛细管极限下多维非守恒可压缩两相流大解的平滑效应

作者: Ling-Yun Shou, Jiayan Wu, Lei Yao, Yinghui Zhang

类别: math.AP, 76T10, 76N10

发布日期: 2025-09-20

arXiv ID: 2509.16607v2

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

详细分析

中文标题: 高毛细管极限下多维非守恒可压缩两相流大解的平滑效应


1. 研究对象和背景

本文研究的是一个多维非守恒可压缩两相流模型,该模型描述了两种不同流体(标记为“+”和“−”相)在空间 $\mathbb{R}^d$($d \geq 2$)中的动态演化,其核心是考虑毛细管效应(capillarity)对界面稳定性和流体行为的调控作用。该系统由以下非守恒形式的偏微分方程组描述:

$$ \begin{cases} \partial_t (\alpha^\pm \rho^\pm) + \operatorname{div}(\alpha^\pm \rho^\pm \mathbf{u}^\pm) = 0, \\ \partial_t (\alpha^\pm \rho^\pm \mathbf{u}^\pm) + \operatorname{div}(\alpha^\pm \rho^\pm \mathbf{u}^\pm \otimes \mathbf{u}^\pm) + \alpha^\pm \nabla P^\pm(\rho^\pm) = \operatorname{div}(\alpha^\pm \boldsymbol{\tau}^\pm + \mathbf{K}^\pm), \end{cases} $$

其中 $\alpha^\pm(x,t) \in [0,1]$ 是两相的体积分数,满足 $\alpha^+ + \alpha^- = 1$;$\rho^\pm(x,t) \geq 0$ 是密度,$\mathbf{u}^\pm(x,t)$ 是速度场,$P^\pm(\rho^\pm)$ 是各自相的压强函数。粘性应力张量 $\boldsymbol{\tau}^\pm$ 定义为:

$$ \boldsymbol{\tau}^\pm := 2\mu^\pm(\rho^\pm) D\mathbf{u}^\pm + \lambda^\pm(\rho^\pm) (\operatorname{div} \mathbf{u}^\pm) I_d, $$

其中 $D\mathbf{u}^\pm = \nabla \mathbf{u}^\pm + (\nabla \mathbf{u}^\pm)^\top$ 是变形率张量,且 $\mu^\pm, \lambda^\pm \in C^\infty(\mathbb{R}^+)$,满足 $\mu^\pm > 0$ 和 $2\mu^\pm + \lambda^\pm > 0$,确保粘性项的椭圆性。关键的毛细管项 $\mathbf{K}^\pm$ 是 Korteweg 型张量,由下式给出:

$$ \mathbf{K}^\pm := \kappa \alpha^\pm \rho^\pm \operatorname{div}\left( m(\alpha^\pm \rho^\pm) \nabla (\alpha^\pm \rho^\pm) \right) I_d + \kappa^\pm \left( m(\alpha^\pm \rho^\pm) - \alpha^\pm \rho^\pm m'(\alpha^\pm \rho^\pm) \right) |\nabla (\alpha^\pm \rho^\pm)|^2 I_d - \kappa^\pm m(\alpha^\pm \rho^\pm) \nabla (\alpha^\pm \rho^\pm) \otimes \nabla (\alpha^\pm \rho^\pm), $$

其中 $\kappa > 0$ 是统一的毛细管系数,$m(s)$ 是一个光滑函数,控制毛细管力的空间分布。为简化分析,作者假设 $m(s) = s$,从而 $\mathbf{K}^\pm$ 可简化为与密度梯度相关的二阶微分算子。

系统中最重要的物理约束是毛细管压强关系:

$$ P^+(\rho^+) - P^-(\rho^-) = f(\alpha^- \rho^-), $$

其中 $f \in C^3([0,\infty))$ 是严格递减函数,刻画了两相间界面的毛细压差。该关系是系统非守恒结构的核心来源,因为它将两相的压强通过界面变量 $\alpha^- \rho^-$ 耦合,导致方程组不能写成标准守恒形式。

初始条件设定为:

$$ (\alpha^\pm, \rho^\pm, \mathbf{u}^\pm)(x,0) \to (\bar{\alpha}^\pm, \bar{\rho}^\pm, 0), \quad |x| \to \infty, $$

其中 $\bar{\alpha}^\pm, \bar{\rho}^\pm > 0$ 是常数平衡态。为便于分析,作者归一化处理:$\bar{\alpha}^\pm \bar{\rho}^\pm = 1$,并设 $\nu = m(1) = 1$。

该模型的数学挑战主要体现在三个方面:
(1)部分耗散性:质量守恒方程无耗散项,仅动量方程含粘性耗散;
(2)退化性:当 $\alpha^\pm$ 或 $\rho^\pm$ 趋近于零时,系统可能退化为单相流,导致系数奇异;
(3)非守恒结构:动量方程中的 $\alpha^\pm \nabla P^\pm$ 项无法通过标准守恒系统理论处理,需引入“弱解”或“重构变量”方法。

研究区域为整个 $\mathbb{R}^d$($d \geq 2$),无边界条件,属于Cauchy 问题。函数空间的选择至关重要:作者使用临界 Besov 空间 $B^s_{p,q}$,特别是 $B^0_{2,1}$ 和 $B^{d/2-1}_{2,1}$,这些空间在尺度变换下不变,是研究可压缩流体全局适定性的自然框架。该空间选择确保了分析的最优性,避免了对初始数据小性的依赖。

该研究的物理背景源于工业与自然现象中的多相流,如核反应堆冷却、油气开采、喷雾燃烧等,其中界面张力(毛细管效应)显著影响混合与分离过程。本文的核心物理现象是:当毛细管系数 $\kappa \to \infty$ 时,两相界面变得极其平滑,系统行为从可压缩的声波主导转向不可压缩的Navier-Stokes扩散主导,即“高毛细管极限”(high-capillarity limit)。这一现象在单相流中已被广泛研究(如 [1,2,8,9]),但在两相流中尚属空白。


2. 主要定理或主要结果

本文的核心成果由四个主要定理构成,分别对应全局存在性、极限收敛、平滑效应与衰减率。

定理 1.1(二维全局大解存在性与收敛):
设 $d = 2$,假设毛细管压强函数满足尖锐稳定性条件:

$$ -\frac{s_-^2(1,1)}{\alpha^-(1,1)} < f'(1) < 0, $$

其中 $s_-^2 = \partial_\rho P^-(\bar{\rho}^-)$ 是负相声速平方。若初始扰动 $(n_{0,\kappa}^\pm, \mathbf{u}_{0,\kappa}^\pm)$ 满足:

$$ (\kappa^{-1/2} n_{0,\kappa}^\pm, \nabla n_{0,\kappa}^\pm, \mathbf{u}_{0,\kappa}^\pm) \in \dot{B}^{0}_{2,1}(\mathbb{R}^2) \quad \text{一致于 } \kappa, $$

则存在 $\kappa_1 > 0$,使得当 $\kappa \geq \kappa_1$ 时,系统(1.20)存在唯一的全局解 $(n^\pm_\kappa, \mathbf{u}^\pm_\kappa)$,满足:

$$ (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, \mathbf{u}^\pm_\kappa) \in C(\mathbb{R}^+; \dot{B}^{0}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{2}_{2,1}). $$

此外,对任意 $2 < p < \infty$,有收敛速率:

$$ \| (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, Q\mathbf{u}^\pm_\kappa) \|_{L^2(\mathbb{R}^+; \dot{B}^{p}_{2,1})} + \| P\mathbf{u}^\pm_\kappa - \mathbf{v}^\pm \|_{L^\infty(\mathbb{R}^+; \dot{B}^{0}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{2}_{2,1})} \leq C \kappa^{-1/4(1 - 2/p)}. $$

定理 1.2(高维全局存在性):
设 $d \geq 3$,在相同稳定性条件(1.22)下,若初始数据满足:

$$ P\mathbf{u}_{0,\kappa}^\pm = \mathbf{v}_0^\pm \in \dot{B}^{d/2 - 1}_{2,1}, \quad (\kappa^{-1/2} n_{0,\kappa}^\pm, \nabla n_{0,\kappa}^\pm, \mathbf{u}_{0,\kappa}^\pm) \in \dot{B}^{d/2 - 1}_{2,1}, \quad \|\mathbf{v}_0^\pm\|_{\dot{B}^{d/2 - 1}_{2,1}} \leq \alpha_0, $$

其中 $\alpha_0$ 是与不可压缩NS方程全局解存在性相关的常数(见(1.31)),则存在 $\kappa_2 > 0$,使得当 $\kappa \geq \kappa_2$ 时,系统有唯一全局解,并满足:

$$ (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, \mathbf{u}^\pm_\kappa) \in C(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1}), $$

且收敛速率:

$$ \| (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, Q\mathbf{u}^\pm_\kappa) \|_{L^2(\mathbb{R}^+; \dot{B}^{p}_{d,p})} + \| P\mathbf{u}^\pm_\kappa - \mathbf{v}^\pm \|_{L^\infty(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1})} \leq C \kappa^{-1/4(d/2 - d/p)}, $$

对 $2 < p \leq 2d/(d-2)$ 成立。

定理 1.3(统一框架:依赖于NS全局解):
这是本文的核心理论框架。设 $d \geq 2$,稳定性条件(1.22)成立。若初始数据满足:

$$ \| (\kappa^{-1/2} n_{0,\kappa}^\pm, \nabla n_{0,\kappa}^\pm, \mathbf{u}_{0,\kappa}^\pm) \|_{\dot{B}^{d/2 - 1}_{2,1}} \leq M_0, $$

且对应的不可压缩NS方程:

$$ \partial_t \mathbf{v}^\pm + \mathbf{v}^\pm \cdot \nabla \mathbf{v}^\pm - \Delta \mathbf{v}^\pm + \nabla \pi^\pm = 0, \quad \operatorname{div} \mathbf{v}^\pm = 0, $$

有全局解 $\mathbf{v}^\pm \in C(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1})$,则存在 $\kappa_0 = \kappa_0(M_0, \nu, f, P^\pm) > 0$,使得当 $\kappa \geq \kappa_0$ 且 $\|P\mathbf{u}_{0,\kappa}^\pm - \mathbf{v}_0^\pm\|_{\dot{B}^{d/2 - 1}_{2,1}} \leq \kappa^{-\delta}$ 时,两相系统有唯一全局解,满足:

$$ (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, \mathbf{u}^\pm_\kappa) \in C(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1}), $$

且误差估计:

$$ \| (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, Q\mathbf{u}^\pm_\kappa) \|_{L^2(\mathbb{R}^+; \dot{B}^{p}_{d,p})} + \| P\mathbf{u}^\pm_\kappa - \mathbf{v}^\pm \|_{L^\infty(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1})} \leq C \kappa^{-\delta}, $$

其中 $(\delta, p)$ 由(1.29)确定:二维时 $\delta = \frac{1}{4}(1 - \frac{2}{p})$,高维时 $\delta = \frac{1}{4}(\frac{d}{2} - \frac{d}{p})$。

定理 1.4(高阶导数的最优衰减与平滑效应):
设 $d \geq 2$,$-d \leq \sigma < d/2 - 1$,若初始扰动还满足:

$$ (\kappa^{-1/2} n_{0,\kappa}^\pm, \nabla n_{0,\kappa}^\pm, \mathbf{u}_{0,\kappa}^\pm, \mathbf{v}_0^\pm) \in \dot{B}^{\sigma}_{2,\infty}, $$

则对任意多阶导数 $\alpha$ 满足 $|\alpha| > \sigma + 1$,有:

$$ \| D^\alpha (\kappa^{-1/2} n^\pm_\kappa, \nabla n^\pm_\kappa, \mathbf{u}^\pm_\kappa)(t) \|_{L^2} \leq C (1 + t)^{-|\alpha|/2 + \sigma/2}, \quad \forall t \geq 1. $$

这是首次证明:即使初始数据仅具临界正则性,毛细管与粘性联合作用仍导致任意阶导数随时间趋于光滑,且衰减速率与热方程最优衰减一致,揭示了长期平滑效应。


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

a) 核心分析工具

本文综合运用了现代调和分析与非线性偏微分方程的前沿工具:

  1. 临界Besov空间框架:作者选择 $\dot{B}^{d/2 - 1}_{2,1}$ 作为初始空间,因其在尺度变换下不变,是研究可压缩流全局适定性的“黄金标准”。该空间允许大初值,且能精确控制非线性项的乘积结构。

  2. 线性化与谱分析:通过对系统(1.20)进行线性化,作者引入新的组合变量:

$$ N_+ := \beta_1 n_+ + (r_+ - \beta_3) n_-, \quad M_+ := \beta_1 Q\mathbf{u}_+ + (r_+ - \beta_3) Q\mathbf{u}_-, $$

其中 $r_\pm$ 是方程 $r^2 - (\beta_1 + \beta_4)r + \beta_1 \beta_4 - \beta_2 \beta_3 = 0$ 的正根(由稳定性条件保证)。这一技巧将原系统解耦为两个Navier-Stokes-Korteweg (NSK) 子系统(见(1.39)-(1.40)),每个子系统包含Schrödinger型色散项 $\sqrt{\kappa} \nabla \Delta N$ 和强耗散项 $\Delta M$。

  1. Strichartz估计与色散估计:由于毛细管项诱导了类似Gross-Pitaevskii方程的色散结构(见摘要),作者在附录B中严格推导了线性系统:
$$ \partial_t N + \sqrt{\kappa} \operatorname{div} M = 0, \quad \partial_t M + \sqrt{r} \sqrt{\kappa} \nabla N - 2\Delta M - \sqrt{\kappa} \nabla \Delta N = 0 $$

的 Strichartz 型估计:

$$ \| N \|_{L^2_t L^p_x} \lesssim \kappa^{-\delta} \| N_0 \|_{\dot{B}^{d/2 - 1}_{2,1}}, $$

这使得在非线性估计中能有效控制“高频振荡”项。

  1. 能量方法与先验估计:作者构造了加权能量泛函(见(2.20)-(2.21)),其正定性依赖于稳定性条件 $f'(1) < 0$。通过结合高频/低频分解与小波分析,得到:
$$ \frac{d}{dt} \mathcal{E}(t) + \mathcal{D}(t) \leq \text{非线性误差项}, $$

其中 $\mathcal{D}(t)$ 包含 $\| \nabla \mathbf{u} \|_{\dot{B}^{d/2}_{2,1}}^2$ 和 $\| \nabla n \|_{\dot{B}^{d/2}_{2,1}}^2$,确保耗散主导。

  1. 不动点定理:在构造解时,作者使用Banach不动点定理,在空间 $X = C(\mathbb{R}^+; \dot{B}^{d/2 - 1}_{2,1}) \cap L^1(\mathbb{R}^+; \dot{B}^{d/2 + 1}_{2,1})$ 上定义映射,证明其压缩性。

b) 关键技术与难题突破

  • 非线性耦合处理:系统中 $\nabla n \cdot \nabla \mathbf{u}$、$\nabla \Delta n \cdot \nabla \mathbf{u}$ 等项难以直接估计。作者采用Bony分解与对偶估计,将乘积项分解为“低频×高频”与“高频×低频”,利用Besov空间的嵌入性质(如 $\dot{B}^{0}_{2,1} \cdot \dot{B}^{d/2}_{2,1} \hookrightarrow \dot{B}^{d/2 - 1}_{2,1}$)控制。

  • 奇异性控制:在压强关系 $P^+ - P^- = f(\alpha^- \rho^-)$ 中,当 $\alpha^- \rho^- \to 0$ 时可能出现奇点。作者通过归一化 $R^\pm = \alpha^\pm \rho^\pm \to 1$,并定义 $\kappa^{-1/2} n^\pm = L(R^\pm)$,将变量替换为“扰动”形式,避免了奇点。

  • 高维与二维的分离处理:在 $d=2$ 时,不可压缩NS方程无需小初值即可全局存在(Danchin [18]),故作者可处理任意大的初始扰动;在 $d \geq 3$ 时,则依赖NS的小初值全局解,这是标准结果(Fujita-Kato)。

  • “ill-prepared”数据处理:作者不要求 $\operatorname{div} \mathbf{u}_{0,\kappa} \to 0$,允许初始速度有强压缩分量,这是对先前工作(如 [20])的重大突破。


4. 与之前工作的比较

本文在多个维度上实现了根本性突破,彻底改变了对两相流全局理论的认知:

  1. 首次突破“小初值”依赖:此前所有工作(如 Evje–Wang–Wen [20], Lai–Wen–Yao [27], Cui–Wang–Yao–Zhu [15])均要求初始扰动足够小,或稳定性条件中包含一个“小参数” $\eta > 0$(见(1.8)):
$$ -\frac{s_-^2}{\alpha^-} < f'(1) < -\frac{\eta s_-^2}{\alpha^-} < 0. $$

本文完全移除该小参数 $\eta$,仅保留尖锐稳定性条件 $-s_-^2 / \alpha^- < f'(1) < 0$,这是最优的,因为当 $f'(1) > 0$ 时系统已知不稳定 [37]。

  1. 维度与正则性突破:早期结果仅限于 $d=3$ 且初始数据在 $L^1 \cap H^3$ 中(Wu–Yao–Zhang [38]),或要求 $L^1$ 衰减性([34])。本文首次在任意 $d \geq 2$ 下,仅要求临界正则性 $\dot{B}^{d/2 - 1}_{2,1}$,且无需 $L^1$ 条件。

  2. 首次建立高毛细管极限的全局收敛:此前研究仅关注 $\kappa \to 0$(vanishing capillarity)或固定 $\kappa$。本文首次证明:当 $\kappa \to \infty$ 时,两相系统全局一致收敛到不可压缩Navier-Stokes方程,且给出了显式收敛率 $\kappa^{-\delta}$,适用于任意大的初始扰动。

  3. 开创“长期平滑效应”理论:定理1.4是历史上首次证明:即使初始数据仅具临界正则性,经过长时间演化,所有阶导数均趋于光滑,且衰减速率与热方程一致。这揭示了毛细管与粘性联合作用下的非平凡正则化机制,是两相流理论的重大进展。

  4. 统一框架的构建:定理1.3将两相系统与不可压缩NS方程的全局存在性直接关联,表明:只要NS方程全局存在,两相系统在高毛细管下也全局存在——这为未来研究提供了普适范式。

综上,本文不仅改进了所有已有结果,更开创了一个全新的研究范式:在大初值、高毛细管、临界正则性框架下,系统展现出色散-耗散协同主导的全新动力学行为,是多相流数学理论的里程碑式工作。


了解领域论文(摘要翻译)

1. Function spaces and trace theorems for maximally subelliptic boundary value problems

关于极大次椭圆边值问题的函数空间与迹定理

作者: Brian Street

类别: math.AP, math.CA, math.FA, 2020: 46E35 (Primary), 46E36, 35G15, and 35H20 (Secondary)

发布日期: 2025-10-14

arXiv ID: 2510.12775v1

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

摘要翻译

中文标题: 关于极大次椭圆边值问题的函数空间与迹定理

摘要翻译: 我们在具有边界的流形上引入了适应于Hörmander向量场的Besov与Triebel–Lizorkin空间,重点关注边界上的所谓非特征点邻域。我们在这些空间中建立了相应的限制算子与迹算子的最优结果,证明了这些算子为重缩映射(retractions),并获得了若干相关结论。本文是即将推出的系列论文中的第二篇,旨在构建极大次椭圆边值问题的统一理论框架,并为此理论奠定了函数空间的分析基础。


2. Rate of convergence to equilibrium for the heated string

加热弦系统趋向平衡的收敛速率

作者: Tomasz Cieślak, Jacek Jendrej, Christian Stinner

类别: math.AP, 35M13

发布日期: 2025-10-14

arXiv ID: 2510.12759v1

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

摘要翻译

中文标题: 加热弦系统趋向平衡的收敛速率

摘要翻译: 在本文中,我们计算了加热弦系统(一类混合型双曲-抛物型偏微分方程组)趋向平衡态的指数收敛速率,且该速率与初始数据无关。作为分析的副产品,我们得到了解的增强时间衰减估计。我们推理的核心工具是对线性化问题在傅里叶空间中的渐近分析,其方法沿袭了Kato的思路;该方法与此前对该系统所获得的估计高度契合。为建立线性问题与非线性原系统之间的关联,我们还需对傅里叶空间中的非线性项进行精细估计,并辅以细致的谱分析。


3. Conditional existence of maximizers for the Tomas-Stein inequality for the sphere

球面上Tomas-Stein不等式极大元的条件存在性

作者: Shuanglin Shao, Ming Wang

类别: math.CA, math.AP

发布日期: 2025-09-12

arXiv ID: 2509.10754v2

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

摘要翻译

中文标题: 球面上Tomas-Stein不等式极大元的条件存在性

摘要翻译: 对于球面 $S^d$ 的一个紧致子集 $\Gamma$,Tomas-Stein 不等式表明映射 $f \mapsto \widehat{f\sigma}$ 从 $L^2(\Gamma, \sigma)$ 到 $L^{2+4/d}(\mathbb{R}^{d+1})$ 是有界的。在假设球面情形的最佳常数严格小于薛定谔方程Strichartz不等式最佳常数的前提下,我们证明了该不等式存在极大元(极值函数),且任意极值序列均包含一个收敛于某极值函数的子序列。我们的方法基于球面上的精细Tomas-Stein不等式与剖面分解技术;其中建立剖面分解中正交性的关键在于,我们运用了Tao关于抛物面的尖锐双线性限制定理,该定理超出了经典Tomas-Stein范围。此前,Frank、Lieb与Sabin \cite{Frank-Lieb-Sabin:2007:maxi-sphere-2d} 通过“缺失质量”方法已获得类似结果。


4. Two-Sided Free Boundary Problems Arising From Branching-Selection Particle Systems

由分支-选择粒子系统引发的双边自由边界问题

作者: Jacob Mercer

类别: math.PR, math.AP, 60J80 (Primary) 35R35 (Secondary)

发布日期: 2025-10-14

arXiv ID: 2510.12701v1

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

摘要翻译

中文标题: 由分支-选择粒子系统引发的双边自由边界问题

摘要翻译: 我们提出并分析了一个源于著名$N$粒子分支布朗运动($N$-BBM)模型推广的双边自由边界问题,该推广模型记为$(N,p)$-BBM,其中选择机制可依据参数$p\in(0,1)$同时作用于最左端与最右端的粒子。本文的主要结果表明,当$N\to\infty$时,$(N,p)$-BBM的经验分布收敛于一个确定性的流体动力学极限,该极限由定义在有限区间上、具有两个移动边界的自由边界问题所描述,其边界条件为带有参数$p$的显式Neumann和Dirichlet条件。我们进一步证明,$(N,p)$-BBM的渐近速度$v_{N,p}$在$N\to\infty$时收敛至极限自由边界问题所对应的唯一行波速度$v_p$。这些结果推广了以往单边模型的结论,并与演化动力学及火焰传播中广泛出现的各类自由边界问题建立了联系。


5. How smooth are restrictions of Besov functions?

Besov 函数的限制有多光滑?

作者: Julien Brasseur

类别: math.FA, math.AP, math.CA

发布日期: 2025-09-09

arXiv ID: 2509.07420v2

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

摘要翻译

中文标题: Besov 函数的限制有多光滑?

摘要翻译: 在先前的研究中,我们证明了 Besov 空间仅在 $q \leq p$ 时才具有限制性质。具体而言,若 $p < q$,则总可以构造一个函数 $f \in B_{p,q}^s(\mathbb{R}^N)$,使得对几乎所有的 $y \in \mathbb{R}^{N-d}$,均有 $f(\cdot,y) \notin B_{p,q}^s(\mathbb{R}^d)$;而当 $q \leq p$ 时,这种“病态”现象则不会出现。我们进一步证明,当函数 $\Psi$ 满足一个仅依赖于 $p$ 和 $q$ 的简单可和性条件时,这些部分映射实际上属于广义光滑度的 Besov 空间 $B_{p,q}^{(s,\Psi)}(\mathbb{R}^d)$。本简短笔记通过证明该刻画的最优性,完善了这一结论。


6. Convolution estimates for the Boltzmann gain operator with hard spheres

硬球模型下玻尔兹曼增益算子的卷积估计

作者: Ioakeim Ampatzoglou, Tristan Léger

类别: math.AP

发布日期: 2025-05-14

arXiv ID: 2505.09554v4

链接: http://arxiv.org/abs/2505.09554v4

摘要翻译

中文标题: 硬球模型下玻尔兹曼增益算子的卷积估计

摘要翻译: 我们为具有硬势能(包括硬球临界情形)的玻尔兹曼方程增益算子建立了新的保持矩的多项式加权卷积估计。我们的方法关键依赖于一种新颖的抵消机制,用于处理能量吸收碰撞这一病态情形(即碰撞中能量仅累积至其中一个出射粒子)。这一困难为硬势能所特有,在麦克斯韦分子模型中并不存在。我们的方法从定量角度证实了如下直觉:尽管能量吸收碰撞以非平凡的概率发生,但其在统计上极为罕见,因而不会影响增益算子的整体平均行为。在技术层面,我们的证明完全基于动力学理论中的工具,如几何恒等式与角平均方法。


7. Asymptotics of the solution of the Cauchy problem for a singularly perturbed system of hyperbolic equations. Part 2. Initial conditions

奇异摄动双曲方程组柯西问题解的渐近行为 第二部分:初始条件

作者: Andrey Nesterov

类别: math.AP, 35C20 (Primary) 35E15, 35L15 (Secondary)

发布日期: 2025-10-14

arXiv ID: 2510.12550v1

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

摘要翻译

中文标题: 奇异摄动双曲方程组柯西问题解的渐近行为 第二部分:初始条件

摘要翻译: 针对描述两根刚性连接弦振动的奇异摄动双曲方程组,构建了单个柯西问题的渐近小参数展开。确定了渐近展开各项所满足的方程(如广义Korteweg-de Vries方程)及其初始条件,并在若干假设条件下,由残差估计了余项的上界。


8. Harnack inequality for Bessel operators

Bessel算子的Harnack不等式

作者: Giorgio Metafune, Luigi Negro, Chiara Spina

类别: math.AP

发布日期: 2025-10-14

arXiv ID: 2510.12529v1

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

摘要翻译

中文标题: Bessel算子的Harnack不等式

摘要翻译: 我们证明了Bessel算子
\begin{align*}
D_t-\Delta_{x} -2a\cdot\nabla_xD_y- D_{yy}- \frac{c}{y} D_y
\[1ex]&=y^{\alpha}\sum_{i,j=1}^{N+1}a_{ij}D_{ij}+y^{\alpha-1}\left(v,\nabla\right)-by^{\alpha-2}
\end{align*}
在带状区域 $[0,T]\times \mathbb{R}^{N+1}_+=\{0 \leq t \leq T, x \in \mathbb{R}^N, y>0\}$ 上满足Neumann边界条件(在 $y=0$ 处)时的唯一性结果与Harnack不等式。


9. Classification and qualitative properties of positive solutions to double-power nonlinear stationary Schrödinger equations

双幂非线性稳态薛定谔方程正解的分类与定性性质

作者: Takafumi Akahori, Slim Ibrahim, Hiroaki Kikuchi, Masataka Shibata, Juncheng Wei

类别: math.AP, 35J20, 35J61, 35B09, 34A34

发布日期: 2025-10-14

arXiv ID: 2510.12484v1

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

摘要翻译

中文标题: 双幂非线性稳态薛定谔方程正解的分类与定性性质

摘要翻译: 本文研究三维空间中双幂非线性稳态薛定谔方程的正径向解。目前已知,当频率足够小时,H¹正解可能不唯一。在适当条件下,除地面态解(其L^∞范数随频率趋于零而趋于零)外,还存在另一正解,该解极小化一个不同的约束变分问题,其L^∞范数随频率趋于零而发散(见定理1.4)。我们将所有小频率下的正解分为两类:地面态解与Aubin-Talenti型解。由此,我们确立了正解的多重性。此外,我们还进一步分析了每一正解的非退化性及其Morse指标。


10. First-order elliptic boundary value problems on manifolds with non-compact boundary

非紧边界流形上的一阶椭圆边值问题

作者: Christian Baer, Lashi Bandara

类别: math.DG, math.AP, Primary 35J56, 47A53, 58J05, 58J32, Secondary 46E35, 53C21

发布日期: 2024-01-31

arXiv ID: 2401.17784v3

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

摘要翻译

中文标题: 非紧边界流形上的一阶椭圆边值问题

摘要翻译: 我们研究作用于具有光滑边界(允许非紧)的光滑流形上向量丛的一阶椭圆微分算子。在非常温和的假设条件下,我们建立了定义在极大定义域上的截面的正则性理论。在附加的几何假设以及适配边界算子的假设下,我们得到了关于极大定义域的迹定理。这使得我们能够系统地研究局部与非局部边值条件,其中Atiyah-Patodi-Singer边值条件作为特例自然出现。此外,我们探讨了诱导半Fredholm算子与Fredholm扩张的相关情形。


11. Harnack inequalities for equations of type prescribed scalar curvature

维数 n ≥ 4 下规定标量曲率型方程的 Harnack 不等式

作者: Samy Skander Bahoura

类别: math.AP

发布日期: 2024-01-03

arXiv ID: 2401.04123v5

链接: http://arxiv.org/abs/2401.04123v5

摘要翻译

中文标题: 维数 n ≥ 4 下规定标量曲率型方程的 Harnack 不等式

摘要翻译: 我们建立了维数 n ≥ 4 下规定标量曲率型方程解的 Harnack 不等式。


12. The Master Equation in a Bounded Domain with Absorption

具有吸收边界条件的有界域中的主方程

作者: Luca Di Persio, Matteo Garbelli, Michele Ricciardi

类别: math.AP, math.PR

发布日期: 2022-03-29

arXiv ID: 2203.15583v6

链接: http://arxiv.org/abs/2203.15583v6

摘要翻译

中文标题: 具有吸收边界条件的有界域中的主方程

摘要翻译: 本文在平均场博弈(MFG)理论框架下,研究具有齐次狄利克雷边界条件的有界域中的主方程。针对N人微分博弈,当玩家触及边界时,其动态过程终止。我们对一类合适的抛物型方程,分析了主方程的适定性及其解的正则性。


13. Semi-Global Existence and Trapped Surface Formation for the Einstein-Vlasov System

爱因斯坦-弗拉索夫系统的大数据半整体存在性与捕获面形成

作者: Nikolaos Athanasiou, Puskar Mondal, Shing-Tung Yau

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

发布日期: 2025-10-14

arXiv ID: 2510.12429v1

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

摘要翻译

中文标题: 爱因斯坦-弗拉索夫系统的大数据半整体存在性与捕获面形成

摘要翻译: 我们证明了在3+1维时空下,针对爱因斯坦-无质量弗拉索夫系统,在无任何对称性假设的前提下,建立了大数据意义下的半整体存在性定理,并揭示了捕获面的动态形成机制。该分析的核心依赖于一套精细校准的估计层级,涵盖韦尔曲率、里奇系数以及分布函数的速度矩,旨在精确刻画耦合系统中几何结构与动力学结构之间的微妙相互作用。弗拉索夫场的存在带来了显著的分析挑战,尤其在可积性与正则性方面。为克服这些困难,我们引入了对里奇系数的精细重正化方法,策略性地将椭圆估计的需求限制于最少数量的里奇系数,并构建了一套与动力学时空切丛质量壳几何结构相适配的精确交换子演算。本成果首次在动能物质背景下,实现了无对称性假设的大数据动态黑洞形成构造。除直接适用于强引力相互作用且无对称性的爱因斯坦-弗拉索夫系统外,我们预期该分析框架将在更广泛的场景中发挥重要作用,例如为无质量弗拉索夫物质下闵可夫斯基时空的非线性稳定性证明提供简化的途径。


14. Nonlocal degenerate parabolic hyperbolic equations on bounded domains

有界域上非局部退化抛物-双曲方程

作者: Nathaël Alibaud, Jørgen Endal, Espen Jakobsen, Ola Mæhlen

类别: math.AP

发布日期: 2023-10-10

arXiv ID: 2310.06453v2

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

摘要翻译

中文标题: 有界域上非局部退化抛物-双曲方程

摘要翻译: 我们研究了定义在有界域上、具有广义狄利克雷边界/外条件的退化混合型抛物-双曲方程

$$ \partial_tu+\text{div}\big(f(u)\big)=\mathcal{L}[b(u)] $$


的适定性问题。该非局部扩散算子 $\mathcal{L}$ 可为任意对称勒维算子(例如分数阶拉普拉斯算子),而函数 $b$ 为单调不减函数,且允许存在退化区域(即 $b'=0$)。我们为该问题提出了一种熵解框架,并在一般假设下证明了有界熵解的唯一性。解的存在性已在另一篇论文中给出。唯一性证明基于 Kružkov 变量加倍技术,并融合了由我们的熵解框架导出的若干先验结果:$L^\infty$ 估计、能量估计、强初始迹、弱边界迹以及一个非局部边界条件。本工作既将非局部理论从全空间推广至有界域,又将局部理论从有界域拓展至非局部情形。与局部理论不同,我们的框架无需预先假定能量估计;相反,能量估计成为该框架的自然推论。然而,与以往的非局部理论不同,这些能量估计在我们的论证中扮演着关键角色。


15. Convergence to a non-explicit steady state in non-factorized kinetic Fokker-Planck equations with (very) weak velocity confinements

具有(极)弱速度约束的非因子化动力学Fokker-Planck方程收敛至非显式稳态

作者: Emeric Bouin, Luca Ziviani

类别: math.AP

发布日期: 2025-10-14

arXiv ID: 2510.12331v1

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

摘要翻译

中文标题: 具有(极)弱速度约束的非因子化动力学Fokker-Planck方程收敛至非显式稳态

摘要翻译: 在本文中,我们证明了具有强空间约束、重尾局部平衡态以及非显式全局稳态的动力学Fokker-Planck方程的若干收敛结果。我们将文献\cite{C21}的结论推广至更广泛的重尾局部平衡态类,并在一大类加权$\sfL^1$空间中获得了收敛速率。为进一步探究非显式稳态的形态,我们辅以数值模拟对结果进行了补充验证。


16. Homogenization of diffusion processes with singular drifts and potentials via unfolding method

基于展开法的具有奇异漂移与势能的扩散过程的均匀化

作者: Toshihiro Uemura, Adisak Seesanea

类别: math.AP, math.PR, 31C25 (Primary) 60J46, 35B27 (Secondary)

发布日期: 2023-06-25

arXiv ID: 2306.14307v2

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

摘要翻译

中文标题: 基于展开法的具有奇异漂移与势能的扩散过程的均匀化

摘要翻译: 本文研究如下形式的椭圆型方程的均匀化问题:

$$ \begin{cases} \mathfrak{L}_{\delta} u_{\delta} + \lambda u_{\delta} = f_{\delta} \qquad \text{在} \;\; D \text{内}, \\ \qquad \quad \;\, u = 0 \qquad \, \text{在} \;\; \partial D \text{上}, \end{cases} $$


其中 $\delta > 0$,$\lambda \in \mathbb{R}$,$D$ 为 $\mathbb{R}^{d}$ 中的有界开集,且 $f_{\delta} \in H^{-1}(D)$。算子 $\mathfrak{L}_{\delta} u = -{\rm div} \left( A^\delta \nabla u + C^\delta u \right) + B^\delta \nabla u + k^\delta u$ 中的扩散系数 $A^\delta$ 一致有界,而漂移项 $B^\delta$、$C^\delta$ 及势能项 $k^\delta$ 可能无界。本文还讨论了该方程所对应扩散过程的均匀化应用。


不相关论文(基本信息)

1. The Linearized Einstein Equations with Sources on Compact Riemannian Manifolds With Boundary

具有边界的紧致黎曼流形上带源项的线性化爱因斯坦方程

作者: Roee Leder

类别: math.DG, math.AP, Primary 35Q76, 53C21, 58A14, Secondary 35S15, 35N25, 35B60

发布日期: 2025-10-14

arXiv ID: 2510.12797v1

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

摘要: We consider compact Riemannian manifolds with arbitrary interior geometry, assuming only that the Ricci tensor vanishes on the boundary. In this setting, we prove that, under the linear Bianchi gauge and extended Cauchy data, the linearized Einstein equations with sources are globally solvable if and only if the source lies in the kernel of a lifted divergence operator and vanishes on the boundary; the solution is unique. This is the first result of its kind without rigid structural restrictions on topology, geometry, or source. The proof combines generalized Hodge theory with unique continuation results.


2. Stable bi-frequency spinor modes as Dark Matter candidates

稳定的双频旋量模式作为暗物质候选者

作者: Andrew Comech, Niranjana Kulkarni, Nabile Boussaïd, Jesús Cuevas-Maraver

类别: math.AP, hep-th, math-ph, math.MP, quant-ph, 35B32, 35B35, 35C08, 35Q41, 37K40, 37N20, 65L07, 81Q05

发布日期: 2024-12-30

arXiv ID: 2501.04027v3

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

摘要: We show that spinor systems with scalar self-interaction, such as the Dirac–Klein–Gordon system with Yukawa coupling or the Soler model, generically have bi-frequency solitary wave solutions. We develop the approach to stability properties of such waves and use the radial reduction to show that indeed the (linear) stability is available for a wide range of parameters. We show that only bi-frequency modes can be dynamically stable and suggest that stable bi-frequency modes can serve as storages of the Dark Matter. The approach is based on linear stability results of one-frequency solitary waves in (3+1)D Soler model, which we obtain as a by-product.


3. Functional a posteriori estimates for the fractional Laplacian problem

分数阶拉普拉斯问题的函数型后验估计

作者: Alexander Nazarov, Sergey Repin

类别: math.AP, 35R11, 65N99

发布日期: 2025-10-14

arXiv ID: 2510.12664v1

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

摘要: The paper is concerned with a posteriori estimates for approximations of boundary value problems generated by the spectral fractional operator. The derivation is based upon the Stinga–Torrea extension that transfers the corresponding nonlocal problem to a local problem of higher dimensionality. The estimates are fully computable and contain no conditions and constants depending on a method or mesh used to compute an approximation. They are valid for any energy admissible approximation of the extended problem.


4. Adapting Noise to Data: Generative Flows from 1D Processes

从一维过程出发的生成流:噪声自适应于数据

作者: Jannis Chemseddine, Gregor Kornhardt, Richard Duong, Gabriele Steidl

类别: stat.ML, cs.LG, math.AP

发布日期: 2025-10-14

arXiv ID: 2510.12636v1

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

摘要: We introduce a general framework for constructing generative models using one-dimensional noising processes. Beyond diffusion processes, we outline examples that demonstrate the flexibility of our approach. Motivated by this, we propose a novel framework in which the 1D processes themselves are learnable, achieved by parameterizing the noise distribution through quantile functions that adapt to the data. Our construction integrates seamlessly with standard objectives, including Flow Matching and consistency models. Learning quantile-based noise naturally captures heavy tails and compact supports when present. Numerical experiments highlight both the flexibility and the effectiveness of our method.


5. Isoperimetric Bounds for Weighted Steklov Eigenvalues with Radial Weights

具径向权重的加权Steklov特征值的等周界估计

作者: Friedemann Brock, Francesco Chiacchio

类别: math.AP, 35P15, 49R05, 49Q10

发布日期: 2025-10-14

arXiv ID: 2510.12631v1

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

摘要: We study the following class of Steklov eigenvalue problems:

$$ \nabla \cdot \bigl( w \nabla u \bigr) = 0 \quad \text{in } \Omega, \qquad \frac{\partial u}{\partial \nu} = \gamma v u \quad \text{on } \partial \Omega, $$

where $w$ and $v$ are prescribed positive radial functions, $\Omega$ is a Lipschitz domain in $\mathbb{R}^N$ with $N \geq 2$ and $\nu$ denotes its outward unit normal. Extending classical results in the unweighted case due to Weinstock, the first author, and others, we establish isoperimetric inequalities for low-order eigenvalues under suitable symmetry assumptions on the domain. In the first part, we consider the case $w(x) = |x|^{\alpha}$ and $v(x) = |x|^{\beta-\alpha}$, where the parameters $\alpha, \beta \in \mathbb{R}$ satisfy appropriate constraints. Our analysis relies on an explicit computation of the spectrum in the radial case, variational principles, and a family of weighted isoperimetric inequalities with ``double density’’. In the second part, we address the case $v \equiv 1$ and $w(x) = W(|x|)$, where $W$ is a non-decreasing, log-convex function. In this setting, the proof relies, among other tools, on a new weighted isoperimetric inequality, which may be of independent interest.


6. Metaplectic operators with quasi-diagonal kernels

具有准对角核的元辛算子

作者: Gianluca Giacchi, Luigi Rodino

类别: math.AP, math.FA, math.SG, quant-ph, 22E46, 47G30, 47F05, 81S30

发布日期: 2025-02-07

arXiv ID: 2502.05089v2

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

摘要: Metaplectic operators form a relevant class of operators appearing in different applications, in the present work we study their Schwartz kernels. Namely, diagonality of a kernel is defined by imposing rapid off-diagonal decay conditions, and quasi-diagonality by imposing the same conditions on the smoothing of the kernel through convolution with the Gaussian. Kernels of metaplectic operators are not diagonal. Nevertheless, as we shall prove, they are quasi-diagonal under suitable conditions. Motivation for our study comes from problems in time-frequency analysis, that we discuss in the last section.


7. Frequency conditions for the global stability of nonlinear delay equations with several equilibria

多平衡点非线性时滞方程全局稳定性的频率条件

作者: Mikhail Anikushin, Andrey Romanov

类别: math.DS, math.AP, 37L15, 34K20, 37L30, 37L45

发布日期: 2023-06-07

arXiv ID: 2306.04716v5

链接: http://arxiv.org/abs/2306.04716v5

摘要: In our adjacent work, we developed a spectral comparison principle for compound cocycles generated by delay equations. It allows to derive frequency inequalities for the uniform exponential stability of such cocycles by means of their comparison with stationary problems. Such inequalities are hard to verify purely analytically, and in this work we develop approximation schemes to verify some of the arising frequency inequalities. Besides some general theoretical results, in applications we stick to the case of scalar equations. By means of the Suarez–Schopf delayed oscillator and the Mackey–Glass equations, we demonstrate applications of the theory to reveal regions in the space of parameters where the absence of closed invariant contours can be guaranteed. Since the frequency inequalities are robust, so close systems also satisfy them, we expect the method to actually imply the global stability, as in known finite-dimensional results utilizing variants of the closing lemma, which is still awaiting developments in infinite dimensions.


8. Spreading properties of the Fisher–KPP equation when the intrinsic growth rate is maximal in a moving patch of bounded size

当内在增长率在有限大小的移动区域内达到最大时,Fisher-KPP方程的传播特性

作者: Thomas Giletti, Léo Girardin, Hiroshi Matano

类别: math.AP

发布日期: 2024-07-31

arXiv ID: 2407.21549v2

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

摘要: This paper is concerned with spreading properties of space-time heterogeneous Fisher–KPP equations in one space dimension. We focus on the case of everywhere favorable environment with three different zones, a left half-line with slow or intermediate growth, a central patch with fast growth and a right half-line with slow or intermediate growth. The central patch moves at various speeds. The behavior of the front changes drastically depending on the speed of the central patch. Among other things, intriguing phenomena such as nonlocal pulling and locking may occur, which would make the behavior of the front further complicated. The problem we discuss here is closely related to questions in biomathematical modelling. By considering several special cases, we illustrate the remarkable diversity of possible behaviors. In particular, when the central patch has constant size and constant speed, we provide a complete set of explicit formulas for the spreading speed.


9. Liouville results for $(p,q)$-Laplacian elliptic equations with source terms involving gradient nonlinearities

涉及梯度非线性源项的 $(p,q)$-拉普拉斯椭圆方程的刘维尔型结果

作者: Mousomi Bhakta, Anup Biswas, Roberta Filippucci

类别: math.AP, 35J60, 35J92, 35J70, 35B08

发布日期: 2025-10-14

arXiv ID: 2510.12486v1

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

摘要: In this paper, we present a series of Liouville-type theorems for a class of nonhomogeneous quasilinear elliptic equations featuring reactions that depend on the solution and its gradient. Specifically, we investigate equations of the form $-\Delta_p u - \Delta_q u = f(u,\nabla u)$ with $p > q > 1$, where the nonlinearity $f$ takes forms such as $u^s|\nabla u|^m$ or $u^s + M|\nabla u|^m$ ($s, m\geq 0$). Our approach is twofold. For cases where the reaction term satisfies $|f(u,\nabla u)|\leq g(u)|\nabla u|^m$ with $m>q$ and $g$ is continuous, we prove that every bounded solution (without sign restriction) in $\mathbb{R}^N$ is constant by means of an Ishii-Lions type technique. In the remaining scenarios, we turn to the Bernstein method. The application of this method to the nonhomogeneous operator requires a nontrivial adaptation, as, roughly speaking, constant coefficients are replaced by functions that may not be bounded from above, which enables us to establish a crucial a priori estimate for the gradient of solutions in any domain $\Omega$. This estimate, in turn, implies the desired Liouville properties on the entire space $\mathbb{R}^N$. As a consequence, we have fully extended Lions Liouville-type result for the Hamilton-Jacobi equation to the $(p,q)$-Laplacian setting, while for the $(p,q)$ generalized Lane-Emden equation, we provide an initial contribution in the direction of the classical result by Gidas and Spruck for $p=q=2$, as well as that of Serrin and Zou for $p=q$. To the best of our knowledge, this is the first paper which studies Liouville properties for equations with nonhomogeneous operator involving source gradient terms.


10. On the minimisation of the Peak-to-average ratio

关于峰均比的最小化

作者: Nikos Katzourakis

类别: math.AP, Primary 35J47, 35J60, Secondary 35D30, 35A15

发布日期: 2025-04-04

arXiv ID: 2504.03972v2

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

摘要: Let $\Omega \Subset \mathbb R^n$ and a continuous function $\mathrm H$ be given, where $n,k,N \in \mathbb N$. For $p\in [1,\infty]$, we consider the functional

$$ \mathrm E_p(u) := \big\| \mathrm H \big(\cdot,u,\mathrm D u, \ldots, \mathrm D^ku \big) \big\|_{\mathrm L^p(\Omega)},\ \ \ u\in \mathrm W^{k,p}(\Omega;\mathbb R^N). $$

We are interested in the $L^\infty$ variational problem

$$ \mathrm C_{\infty,p}(u_\infty)\, =\, \inf \Big\{\mathrm C_{\infty,p}(u) \ : \ u\in \mathrm W^{k,\infty}_\varphi(\Omega;\mathbb R^N), \ \mathrm E_1(u)\neq 0 \Big\}, $$

where $\varphi\in \mathrm W^{k,\infty}(\Omega;\mathbb R^N)$, $p$ is fixed, and

$$ \mathrm C_{\infty,p}(u)\, := \, \frac{\mathrm E_\infty(u)}{\mathrm E_p(u)} . $$

The variational problem is ill-posed. $\mathrm C_{\infty,2}$ is known as the Crest factor" and arises as the peak–to–average ratio" problem in various applications, including eg. nuclear reactors and signal processing in sound engineering. We solve it by characterising the set of minimisers as the set of strong solutions to the eigenvalue Dirichlet problem for the fully nonlinear PDE

$$ \left\{ \ \ \begin{array}{ll} \big| \mathrm H \big(\cdot,u,\mathrm D u, \ldots, \mathrm D^ku \big) \big|= \Lambda, & \text{ a.e.\ in }\Omega, \\ u = \varphi, & \text{ on }\partial \Omega,\\ \mathrm D u = \mathrm D \varphi, & \text{ on }\partial\Omega, \vdots & \vdots \\ \mathrm D^{k-1}u = \mathrm D^{k-1}\varphi, & \text{ on }\partial\Omega. \end{array} \right. $$

Under appropriate assumptions for $\mathrm H$, we show existence of infinitely-many solutions $(u,\Lambda) \in \mathrm W^{k,\infty}_\varphi(\Omega;\mathbb R^N) \times [\Lambda_*,\infty)$ for $\Lambda_*\geq0$, by utilising the Baire Category method for implicit PDEs. In the case of $k=1$ and $n=N$, these assumptions do not require quasiconvexity.


11. Metric Sobolev spaces II: dual energies and divergence measures

度量Sobolev空间II:对偶能量与散度测度

作者: Luigi Ambrosio, Toni Ikonen, Danka Lučić, Enrico Pasqualetto

类别: math.FA, math.AP, math.MG, 49J52, 46E35, 53C23, 46N10, 46E15, 31C15, 28A12

发布日期: 2025-10-14

arXiv ID: 2510.12424v1

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

摘要: This is the second of two works concerning the Sobolev calculus on metric measure spaces and its applications. In this work, we focus on several approaches to vector calculus in the non-smooth setting of complete and separable metric spaces equipped with a boundedly-finite Borel measure. More precisely, we study different notions of (co)vector fields and derivations appearing in the literature, as well as their mutual relation. We also carry forward a thorough investigation of gradients, divergence measures, and Laplacian measures, together with their applications in potential analysis (for example, regarding the condenser capacity) and in the study of duality properties of Sobolev spaces. Most of the results are obtained for the full range of exponents $p\in[1,\infty)$ and without finiteness assumption on the measure.


12. Semiclassical Resolvent Estimates for the Magnetic Schr{Ö}dinger Operator

磁薛定谔算子的半经典预解式估计

作者: Georgi Vodev

类别: math.AP

发布日期: 2025-01-13

arXiv ID: 2501.07271v4

链接: http://arxiv.org/abs/2501.07271v4

摘要: We obtain semiclassical resolvent estimates for the Schr{"o}dinger operator (ih$\nabla$ + b)^2 + V in R^d , d $\ge$ 3, where h is a semiclassical parameter, V and b are real-valued electric and magnetic potentials independent of h. Under quite general assumptions, we prove that the norm of the weighted resolvent is bounded by exp(Ch^{-2} log(h^{ -1} )) . We get better resolvent bounds for electric potentials which are H{"o}lder with respect to the radial variable and magnetic potentials which are H{"o}lder with respect to the space variable. For long-range electric potentials which are Lipschitz with respect to the radial variable and long-range magnetic potentials which are Lipschitz with respect to the space variable we obtain a resolvent bound of the form exp(Ch^{-1}) .


13. Regularity lost: the fundamental limitations and constraint qualifications in the problems of elastoplasticity

规律的丧失:弹塑性问题中的根本局限与约束条件

作者: Ivan Gudoshnikov

类别: math.AP, math-ph, math.DS, math.FA, math.MP, math.OC, 74C05, 49J40, 47J22, 47B02, 47B93

发布日期: 2024-12-17

arXiv ID: 2412.13068v2

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

摘要: We investigate the existence and non-existence of a function-valued strain solution in various models of elastoplasticity from the perspective of the constraint-based ``dual’’ formulations. We describe abstract frameworks for linear elasticity, elasticity-perfect plasticity and elasticity-hardening plasticity in terms of adjoint linear operators and convert them to equivalent formulations in terms of differential inclusions (the sweeping process in particular). Within such frameworks we consider several manually solvable examples of discrete and continuous models. Despite their simplicity, the examples show how for discrete models with perfect plasticity it is possible to find the evolution of stress and strain (elongation), yet continuum models within the same framework may not possess a function-valued strain. Although some examples with such phenomenon are already known, we demonstrate that it may appear due to displacement loading. The central idea of the paper is to explain the loss of strain regularity in the dual formulation by the lack of additivity of the normal cones and the failure of Slater’s constraint qualification. In contrast to perfect plasticity, models with hardening are known to be well-solvable for strains. We show that more advanced constraint qualifications can help to distinguish between those cases and, in the case of hardening, ensure the additivity of the normal cones, which means the existence of a function-valued strain rate.


14. Elliptic Harnack inequalities for mixed local and nonlocal $p$-energy form on metric measure spaces

度量测度空间上混合局部与非局部 $p$-能量的椭圆Harnack不等式

作者: Aobo Chen, Zhenyu Yu

类别: math.AP, math.FA, 30L15, 31C45, 31E05, 39B62, 46E36

发布日期: 2025-10-14

arXiv ID: 2510.12404v1

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

摘要: We derive the weak elliptic Harnack inequality and the elliptic Harnack inequality for mixed local and nonlocal $p$-energy form on metric measure spaces. The proof relies on the Poincar'e inequality, the cutoff Sobolev inequality, and some mild conditions on the jump measure. Our approach is based on the De Giorgi-Nash-Moser method, extending the corresponding result for Dirichlet forms without a killing part (the case $p=2$).


15. The epsilon-regularity theorem for Brakke flows near triple junctions

布拉克流在三重交点附近的ε-正则性定理

作者: Salvatore Stuvard, Yoshihiro Tonegawa

类别: math.AP, math.DG, 53E10, 35B65, 49Q15

发布日期: 2025-10-03

arXiv ID: 2510.02969v2

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

摘要: We establish the $\varepsilon$-regularity theorem for $k$-dimensional, possibly forced, Brakke flows near a static, multiplicity-one triple junction. This result provides the parabolic analogue to L. Simon’s foundational work on the singular set of stationary varifolds and confirms that the regular structure of triple junctions persists under weak mean curvature flow. The regularity holds provided the flow satisfies a mild structural assumption on its 1-dimensional slices taken orthogonal to the junction’s $(k-1)$-dimensional spine, which prohibits certain topological degeneracies. We prove that this assumption is automatically satisfied by two fundamental classes of flows where such singularities are expected: codimension-one multi-phase flows, such as the canonical $\mathrm{BV}$-Brakke flows constructed by the authors, and flows of arbitrary codimension with the structure of a mod 3 integral current, which arise from Ilmanen’s elliptic regularization. For such flows, therefore, the Simon type regularity holds unconditionally.