生成时间: 2026年02月04日 04:23:11

重点关注论文数量: 4

了解领域论文数量: 3

不相关论文数量: 16

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

1. 带阻尼的可压欧拉系统在混合型Besov空间中的全局适定性与松弛极限

The compressible Euler system with damping in hybrid Besov spaces: global well-posedness and relaxation limit

作者: Timothée Crin-Barat, Zihao Song

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03575v1

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

详细分析

中文标题: 带阻尼的可压欧拉系统在混合型Besov空间中的全局适定性与松弛极限


1. 研究对象和背景

本文研究的是 多维带阻尼的可压欧拉系统(compressible Euler system with damping),其原始形式为定义在 $\mathbb{R}^d$($d \geq 1$)上的如下 Cauchy 问题:

$$ \begin{cases} \partial_t \rho + \operatorname{div}(\rho u) = 0, \\ \partial_t (\rho u) + \operatorname{div}(\rho u \otimes u) + \nabla P(\rho) + \dfrac{1}{\varepsilon} \rho u = 0, \\ (\rho, u)\big|_{t=0} = (\rho_0, u_0), \end{cases} \tag{1.1} $$

其中 $\rho = \rho(x,t) > 0$ 表示密度,$u = u(x,t) \in \mathbb{R}^d$ 表示速度场,$\varepsilon > 0$ 是时间松弛参数,压力 $P$ 满足 $\gamma$-律状态方程:

$$ P(\rho) = A \rho^\gamma, \quad \gamma > 1,\; A > 0. $$

该系统描述了具有线性摩擦阻尼项 $\frac{1}{\varepsilon}\rho u$ 的可压缩流体动力学行为。物理上,它建模了气体在强耗散环境(如多孔介质、粘性主导但忽略剪切应力的近似)中的演化;数学上,它是典型的一阶拟线性双曲系统叠加线性耗散源项,属于 部分耗散双曲系统(partially dissipative hyperbolic system)的范畴。

为便于分析,作者将系统围绕常值平衡态 $(\bar{\rho}, 0)$ 进行线性化,并设 $\bar{\rho} = 1$(不失一般性)。引入声速变量 $c$ 和重缩放速度 $v$:

$$ c_\varepsilon(\varepsilon t, x) \triangleq \sqrt{\frac{1}{\gamma} \frac{\partial P}{\partial \rho}(\rho(t,x))} = \sqrt{A \gamma} \, \rho^{\frac{\gamma-1}{2}}, \quad v_\varepsilon(\varepsilon t, x) \triangleq \frac{1}{\varepsilon} u(t,x), $$

并记 $\check{\gamma} = \frac{\gamma - 1}{2}$,$\bar{c} = \sqrt{A \gamma}$。定义扰动变量 $c = c_\varepsilon - \bar{c}$、$v = v_\varepsilon$,则系统 (1.1) 等价于如下对称形式:

$$ \begin{cases} \partial_t c + v \cdot \nabla c + \check{\gamma}(c + \bar{c}) \operatorname{div} v = 0, \\ \varepsilon^2 (\partial_t v + v \cdot \nabla v) + \check{\gamma}(c + \bar{c}) \nabla c + v = 0. \end{cases} \tag{2.8} $$

该形式更利于能量估计与谱分析,且显式呈现了系统的对称双曲–耗散耦合结构。

本文的核心目标之一是研究其在 $\varepsilon \to 0$ 下的松弛极限(relaxation limit),即在所谓扩散尺度下:

$$ (\rho^\varepsilon, u^\varepsilon)(t,x) \triangleq \big( \rho(\varepsilon^{-1}t, x),\; \varepsilon^{-1} u(\varepsilon^{-1}t, x) \big), $$

此时系统变为:

$$ \begin{cases} \partial_t \rho^\varepsilon + \operatorname{div}(\rho^\varepsilon u^\varepsilon) = 0, \\ \varepsilon^2 \partial_t (\rho^\varepsilon u^\varepsilon) + \varepsilon^2 \operatorname{div}(\rho^\varepsilon u^\varepsilon \otimes u^\varepsilon) + \nabla P(\rho^\varepsilon) + \rho^\varepsilon u^\varepsilon = 0. \end{cases} \tag{1.3} $$

当 $\varepsilon \to 0$ 时,期望 $(\rho^\varepsilon, u^\varepsilon) \to (N, V)$,其中极限函数满足多孔介质方程(porous medium equation):

$$ \begin{cases} \partial_t N - \Delta P(N) = 0, \\ V = - \dfrac{\nabla P(N)}{N}. \end{cases} \tag{1.4} $$

第二式即著名的 Darcy 定律,表明极限速度由压力梯度驱动,体现强耗散下的准静态平衡。

值得注意的是:该极限系统 (1.4) 在 $L^p$ 类空间中对任意 $p \in [2,\infty)$ 都具有良好适定性(见 Proposition 5.1),而此前工作(如 [12])仅能在受限范围 $p \in \big[2, \max\{4, \frac{2d}{d-2}\}\big]$ 内建立原系统的适定性——这构成了本文方法论革新的根本动机。

在函数空间方面,本文采用混合型齐次 Besov 空间(hybrid homogeneous Besov spaces),其核心思想是:

  • 对低频部分($|\xi| \lesssim \varepsilon^{-1}$),使用 $L^p$-型范数($p \in [2,\infty)$),以匹配极限方程的自然正则性;
  • 对高频部分($|\xi| \gtrsim \varepsilon^{-1}$),使用 $L^2$-型范数,以保障双曲–耗散耦合的能量估计可行性(避免导数损失);
  • 引入频率阈值 $J_\varepsilon = \lfloor \log_2(\varepsilon^{-1}) \rfloor + k_0$,并定义分解:
    $$ f = f^{\ell,\varepsilon} + f^{h,\varepsilon}, \quad f^{\ell,\varepsilon} \triangleq \sum_{j < J_\varepsilon} \dot{\Delta}_j f, \quad f^{h,\varepsilon} \triangleq \sum_{j \geq J_\varepsilon} \dot{\Delta}_j f. $$

进一步,为精细刻画非线性相互作用,作者在低频与高频之间插入中间频率区(medium-frequency regime),构成三段式频域划分,这是突破 $p$ 上界的关键技术设计。

此外,初始数据被要求处于临界正则性(critical regularity)框架:低频部分属 $\dot{B}^{d/p}_{p,1}$,高频部分属 $\dot{B}^{d/2+1}_{2,1}$,整体属于“混合初值空间” $E^\varepsilon_0$,其范数为:

$$ X^\varepsilon_0 \triangleq \| (c_0, \varepsilon v_0) \|_{\dot{B}^{d/p}_{p,1}}^{\ell,\varepsilon} + \| (\varepsilon c_0, \varepsilon^2 v_0) \|_{\dot{B}^{d/2+1}_{2,1}}^{h,\varepsilon}. $$

该设定既保证了小初值条件下的全局存在唯一性,又为松弛极限提供足够强的收敛性控制。


2. 主要定理或主要结果

本文有两个核心定理,分别确立全局适定性与松弛极限的强收敛性及定量误差估计。

定理 2.1(全局适定性)

设 $d \geq 1$,$p \in [2,\infty)$。存在常数 $\eta > 0$ 与整数 $k_0 \in \mathbb{Z}$,使得对任意 $\varepsilon > 0$,若初始扰动满足

$$ \| (c_0, \varepsilon v_0) \|_{\dot{B}^{d/p}_{p,1}}^{\ell,\varepsilon} + \| (\varepsilon c_0, \varepsilon^2 v_0) \|_{\dot{B}^{d/2+1}_{2,1}}^{h,\varepsilon} \leq \eta, $$

则系统 (2.8) 存在唯一的全局经典解 $(c,v)$,且满足一致先验估计:

$$ X^\varepsilon(t) \leq C X^\varepsilon_0, \quad \forall t > 0, $$

其中 $C > 0$ 为与 $\varepsilon$ 无关的普适常数,而 $X^\varepsilon(t)$ 是如下混合型能量泛函:

$$ X^\varepsilon(t) \triangleq X^{\ell,\varepsilon}(t) + X^{h,\varepsilon}(t), $$

$$ X^{\ell,\varepsilon}(t) \triangleq |c|{L^\infty_T(\dot{B}^{d/p}{p,1}) \cap L^1_T(\dot{B}^{d/p+2}_{p,1})}^{\ell}

  • \varepsilon |v|{L^\infty_T(\dot{B}^{d/p}{p,1})}^{\ell}
  • |v|{L^2_T(\dot{B}^{d/p}{p,1}) \cap L^1_T(\dot{B}^{d/p+1}_{p,1})}^{\ell}
  • \frac{1}{\varepsilon} |W|{L^1_T(\dot{B}^{d/p}{p,1})}^{\ell},
    $$ 其中 $W \triangleq v + \check{\gamma}(c+\bar{c}) \nabla c$ 为**有效速度**(effective velocity);而 $$
    X^{h,\varepsilon}(t) \triangleq \varepsilon |c|{L^\infty_T(\dot{B}^{d/2+1}{2,1})}^{h}
  • \frac{1}{\varepsilon} |c|{L^1_T(\dot{B}^{d/2+1}{2,1})}^{h}
  • \varepsilon^2 |v|{L^\infty_T(\dot{B}^{d/2+1}{2,1})}^{h}
  • |v|{L^1_T(\dot{B}^{d/2+1}{2,1})}^{h}. $$

该定理首次在任意 $p \in [2,\infty)$ 范围内建立了带阻尼可压欧拉系统的小初值全局适定性,显著拓宽了先前工作的适用范围。

定理 2.2(松弛极限与误差估计)

设 $d \geq 1$,$p \in [2,\infty)$,$(c,v)$ 为定理 2.1 所给解,对应密度 $\rho^\varepsilon$ 由 (1.2) 和 (2.7) 定义;设 $N$ 是多孔介质方程 (1.4) 的全局解,初值 $N_0 \in \dot{B}^{d/p}_{p,1}$。若初始误差满足

$$ \|\rho^\varepsilon_0 - N_0\|_{\dot{B}^{d/p - \delta}_{p,1}} \leq C \varepsilon^\delta, \quad \delta \in (0,1], $$

则当 $\varepsilon \to 0$ 时,有如下强收敛性:

$$ \rho^\varepsilon - N \to 0 \quad \text{strongly in } L^\infty(\mathbb{R}_+; \dot{B}^{d/p - \delta}_{p,1}) \cap L^{1 + 2\delta}(\mathbb{R}_+; \dot{B}^{d/p + 1}_{p,1}), $$

$$ v^\varepsilon + \frac{\nabla P(\rho^\varepsilon)}{\rho^\varepsilon} \to 0 \quad \text{strongly in } L^1(\mathbb{R}_+; \dot{B}^{d/p}_{p,1}). $$

并有定量误差估计:对任意 $\delta \in (0,1]$ 与 $r \in [1,2)$,

$$ \|\rho^\varepsilon - N\|_{L^\infty_T(\dot{B}^{d/p - \delta}_{p,1}) \cap L^{1+2\delta}_T(\dot{B}^{d/p + 1}_{p,1})} \leq C \varepsilon^\delta, $$

$$ \left\| v^\varepsilon + \frac{\nabla P(\rho^\varepsilon)}{\rho^\varepsilon} \right\|_{L^r_T(\dot{B}^{d/p}_{p,1})} \leq C \varepsilon. $$

该定理不仅给出强收敛性,还提供了最优阶误差估计($\mathcal{O}(\varepsilon^\delta)$),且适用于全范围 $p \in [2,\infty)$,是对 [12] 中 $\delta = 1$、$p$ 受限情形的根本性推广。


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

a) 核心分析学工具

本文的证明基于Littlewood–Paley 理论与频域分解方法,结合以下几类关键分析工具:

  • 混合型 Besov 空间理论:利用齐次 Besov 空间 $\dot{B}^s_{p,r}$ 的定义与基本嵌入性质,特别是:

    $$ \dot{B}^{d/2}_{2,1} \hookrightarrow L^\infty, \quad \dot{B}^{d/p}_{p,1} \hookrightarrow L^\infty \quad (p < \infty), $$

    以及 Bernstein 不等式:

    $$ \|\dot{\Delta}_j f\|_{L^q} \lesssim 2^{j(d/p - d/q)} \|\dot{\Delta}_j f\|_{L^p}, \quad q \ge p. $$
  • 谱分析与模型线性化:对线性化系统

    $$ \begin{cases} \partial_t c + \operatorname{div} v = 0, \\ \varepsilon^2 \partial_t v + \nabla c + v = 0, \end{cases} $$

    进行 Fourier 变换,得矩阵 $H(\xi) = \begin{pmatrix} 0 & i\xi \\ i\xi & \varepsilon^{-1} \end{pmatrix}$,其特征值渐近行为为:

    • 当 $|\xi| \ll \varepsilon^{-1}$:$\lambda_+ \sim \varepsilon^{-1}$(快衰减模),$\lambda_- \sim \varepsilon |\xi|^2$(慢扩散模);
    • 当 $|\xi| \gg \varepsilon^{-1}$:$\Re(\lambda_\pm) \sim \varepsilon^{-1}$(全频指数衰减)。

    此揭示了低频热型行为与高频阻尼主导行为的本质差异,为混合空间设计提供物理依据。

  • Lyapunov 泛函构造(hypocoercivity 方法):在高频区,构造组合能量泛函(如 $\|c\|_{\dot{B}^{d/2+1}_{2,1}}^2 + \|v\|_{\dot{B}^{d/2+1}_{2,1}}^2 + \varepsilon \langle \nabla c, v \rangle$)以捕获交叉耗散;在低频区,则借助有效速度 $w = v + \nabla c$ 将系统改写为:

    $$ \begin{cases} \partial_t c - \Delta c = - \operatorname{div} w, \\ \varepsilon \partial_t w + w = -\varepsilon \nabla \Delta c - \varepsilon \nabla \operatorname{div} w, \end{cases} $$

    从而实现模态解耦,使两方程可分别在 $L^p$ 框架下独立估计。

  • 不动点论证(contraction mapping principle):在构造好的混合 Banach 空间 $\mathcal{X}^\varepsilon_T$(由 $X^\varepsilon(t)$ 控制)中,将原系统视为映射 $\Phi: (c,v) \mapsto (\tilde{c},\tilde{v})$,通过精细的双线性估计证明 $\Phi$ 是压缩映射,从而获得局部解;再结合先验估计得到全局解。

b) 关键技术:非线性项处理与频率交互控制

最大技术难点在于估计形如 $(f^\ell g^\ell)^h$、$(f^h g^\ell)^h$ 等跨频非线性项。传统混合空间方法受限于嵌入失效:$\dot{B}^{d/p}_{p,1} \not\hookrightarrow \dot{B}^{d/2}_{2,1}$ 当 $p > 2$,导致 $(f^\ell g^\ell)^h$ 无法直接控制。

本文提出三项突破性技巧:

  1. 高频截断消失引理(Key vanishing lemma):
    利用支撑集性质,证明存在常数 $a_0 \in (0,1)$,使得对任意 $j \ge J_\varepsilon$,只要 $j_1, j_2 < J_\varepsilon$ 且 $\max\{2^{j_1}, 2^{j_2}\} < a_0 2^{J_\varepsilon}$,则

    $$ \operatorname{supp} \mathcal{F}[\dot{\Delta}_j(\dot{\Delta}_{j_1} f \, \dot{\Delta}_{j_2} g)] = \emptyset. $$

    因此,

    $$ (f^\ell g^\ell)^h = \sum_{\substack{j \ge J_\varepsilon \\ j_1,j_2 < J_\varepsilon \\ \max\{j_1,j_2\} \ge \log_2(a_0) + J_\varepsilon}} \dot{\Delta}_j(\dot{\Delta}_{j_1} f \, \dot{\Delta}_{j_2} g). $$

    即:低频×低频 的高频投影,必含至少一个“中频因子”(频率 $\sim a_0 2^{J_\varepsilon}$)。

  2. 多段频率插值框架(Iterated hybrid decomposition):
    引入 $m$ 个中间频段 $\{J^{(k)}\}_{k=1}^m$,对应 $L^{p_k}$ 空间,满足 $2 = p_0 < p_1 < \cdots < p_m = p$。对每一段,应用 Hölder 型乘子估计:

    $$ \|fg\|_{\dot{B}^{s}_{p_k,1}} \lesssim \|f\|_{\dot{B}^{s_1}_{p_i,1}} \|g\|_{\dot{B}^{s_2}_{p_j,1}}, \quad \frac{1}{p_k} = \frac{1}{p_i} + \frac{1}{p_j}. $$

    通过迭代选择 $p_i$,可将 $p$ 的上界逐步提升至任意有限值。

  3. 精炼的乘子与交换子估计:
    建立如下新型估计(Proposition 3.1):

    $$ \|[ \dot{\Delta}_j, f ] g\|_{L^p} \lesssim 2^{-j} \|\nabla f\|_{L^{p_1}} \|g\|_{L^{p_2}}, \quad \frac{1}{p} = \frac{1}{p_1} + \frac{1}{p_2}, $$

    以及针对混合频段的乘积估计:

    $$ \|(f^\ell g^\ell)^{m}\|_{\dot{B}^{d/p_1}_{p_1,1}} \lesssim \|f^\ell\|_{\dot{B}^{d/p}_{p,1}} \|g^\ell\|_{\dot{B}^{d/p}_{p,1}}, $$

    其中 $m$ 表示中频投影。此类估计是实现全 $p$ 范围适定性的基石。


4. 与之前工作的比较

a) 突破与改进

本文相较先前工作(尤其是 [12], [10], [40])实现了三重根本性突破:

  • 正则性假设的彻底放宽:
    [12] 要求 $p \in \big[2, \max\{4, \frac{2d}{d-2}\}\big]$,其根源在于:

    • 使用 $L^p \times L^{p^*} \to L^2$ Hölder 不等式估计 $(f^\ell g^\ell)^h$,需 $1/p + 1/p^* = 1/2$,得 $p^* = 2p/(p-2)$,进而由嵌入 $\dot{B}^{d/p - d/p^*}_{p,1} \hookrightarrow L^{p^*}$ 要求 $d/p - d/p^* > 0$,即 $p \le 4$;
    • 同时为控制 $\nabla f^\ell$,需 $\dot{B}^{d/p - d/p^* + 1}_{p,1} \hookrightarrow L^{p^*}$,导出 $p \le 2d/(d-2)$(Sobolev 嵌入限制)。
      本文通过中频中介与支撑集消失性,完全绕开上述嵌入瓶颈,将 $p$ 推广至整个 $[2,\infty)$,首次实现与极限方程正则性完全匹配。
  • 维度无关性:
    [12] 中 $\frac{2d}{d-2}$ 在 $d=2$ 时发散,需单独处理;而本文 $p < \infty$ 对所有 $d \ge 1$ 统一成立,包括临界维 $d = 2$。

  • 收敛速率的统一性与最优性:
    [12] 仅给出 $\delta = 1$ 的 $\mathcal{O}(\varepsilon)$ 误差;本文给出连续依赖于初值正则性的误差阶 $\varepsilon^\delta$($\delta \in (0,1]$),且对任意 $p$ 成立,体现了更强的稳定性与鲁棒性。

综上,本文不仅拓展了适用初值类,更构建了一套可扩展的混合频域分析范式,为一般部分耗散双曲系统(满足 SK 条件)的临界正则性理论提供了新模板。


2. 关于一维可压Navier-Stokes方程弱解的研究:弱范数下解对初边值数据的Lipschitz连续依赖性及其均匀化近似的误差分析

On weak solutions to the 1d compressible Navier-Stokes equations: a Lipschitz continuous dependence on data in weaker norms and an error of their homogenization

作者: Alexander Zlotnik

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03481v1

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

详细分析

中文标题: 关于一维可压Navier-Stokes方程弱解的研究:弱范数下解对初边值数据的Lipschitz连续依赖性及其均匀化近似的误差分析


研究对象和背景

本文研究的是一维可压缩Navier-Stokes方程组(1D compressible Navier-Stokes system)在拉格朗日质量坐标 $(x,t)\in Q = \Omega \times (0,T)$ 下的全局弱解理论,其中 $\Omega = (0,X)$ 是有界区间,$T>0$ 为任意给定终时。该系统描述粘性热传导气体的一维流动,由三个守恒律耦合构成:质量守恒、动量守恒与内能守恒。其核心变量为:

  • $\eta > 0$: 比容(specific volume),即密度 $\rho = 1/\eta$ 的倒数;
  • $u$: 速度(velocity);
  • $\theta > 0$: 绝对温度(absolute temperature);
  • $x_e$: 欧拉坐标(Eulerian coordinate),满足 $D_t x_e = u$,且初始值 $x_e|_{t=0} = x_0 = I\eta^0 := \int_0^x \eta^0(\zeta)\,d\zeta$。

系统方程如下(见原文第6–7页,公式(20)–(23)):

$$ D_t \eta = D_x u, $$
$$ D_t u = D_x \sigma + g[x_e], $$
$$ c_V D_t \theta = D_x \pi + \sigma D_x u + f[x_e], $$
$$ D_t x_e = u,\quad \sigma = \nu \rho D_x u - p,\quad \rho = \frac{1}{\eta},\quad p = k \rho \theta,\quad \pi = \lambda \rho D_x \theta. $$

其中 $c_V > 0$, $\nu > 0$, $k > 0$, $\lambda > 0$ 为正常物理常数;$g(\chi,x,t)$ 和 $f(\chi,x,t) \geq 0$ 分别为体力密度与热源强度,通过复合 $h[x_e](x,t) = h(x_e(x,t),x,t)$ 引入欧拉依赖;$D_t$ 与 $D_x$ 表示偏导 $\partial_t$ 与 $\partial_x$。

该系统配备非齐次边界条件(三类标准类型,对应不同物理设定):

  • 类型 $m=1$: $u|_{x=0,X} = u_b = (u_0,u_X)$(速度固定);
  • 类型 $m=2$: $\sigma|_{x=0} = -p_0$, $u|_{x=X} = u_X$(左端压力、右端速度控制);
  • 类型 $m=3$: $\sigma|_{x=0,X} = -p_b = (-p_0,-p_X)$(两端压力控制);

以及热通量边界条件:$\pi|_{x=0,X} = \pi_b = (\pi_0,\pi_X)$。初始条件为:

$$ (\eta,u,\theta,x_e)|_{t=0} = (\eta^0,u^0,\theta^0,x_0). $$

值得注意的是,该问题的函数空间框架具有显著特殊性:

  • 解属于 $L^\infty(Q)$-有界但未必连续的函数类:$\eta,1/\eta \in L^\infty(Q)$,$u,\theta \in V_2(Q)$,其中
    $$ V_2(Q) := \left\{ \varphi \in L^{2,\infty}(Q) : D_x \varphi \in L^2(Q) \right\},\quad \|\varphi\|_{V_2(Q)} = \|\varphi\|_{L^{2,\infty}(Q)} + \|D_x \varphi\|_{L^2(Q)}; $$
  • 初始数据允许大间断(large discontinuous data):$\eta^0 \in L^\infty(\Omega)$,$u^0 \in L^\infty(\Omega)$,$\theta^0 \in L^2(\Omega)$ 或仅 $L^1(\Omega)$,远弱于经典强解所需的 $H^1(\Omega)$ 正则性;
  • 边界数据属 $W^{1,1}(0,T)$ 或 $L^{4/3}(0,T)$,不强制连续或光滑;
  • 所有估计均在弱拓扑下进行:例如解差 $(\hat\eta-\eta,\hat u-u,\hat\theta-\theta)$ 的控制范数为
    $$ \|\cdot\|_{C(0,T;L^2(\Omega))} \times \|\cdot\|_{L^2(Q)\cap L^\infty(0,T;H^{-1;m})} \times \|\cdot\|_{L^2(Q)}, $$
    其中 $H^{-1;m}$ 是 $H^1_m(\Omega)$ 的对偶空间,依边界类型 $m$ 定义(如 $m=1$ 对应零边界值子空间)。

这一设定直接面向真实物理场景中常见的不规则初边值(如激波初态、突变边界驱动),突破了传统理论对高正则性的依赖,体现了对适定性本质——即解如何稳定响应粗糙输入——的深层追问。


主要定理或主要结果

本文的核心成果是以下Lipschitz型连续依赖定理(Theorem 1,原文第9–10页):

Theorem 1. 设 $P_m$ 与 $\widehat{P}_m$ 分别为原问题与扰动问题(含额外扰动项 $\beta,\gamma,g_1,g_2,f-f'$ 等)的 IBVP,其数据满足假设 $(C_1)$–$(C_3)$ 与 $(\widehat{C}_1)$–$(\widehat{C}_3)$(见原文第7–9页)。令 $q_e \in [2,\infty]$ 满足 $(C_2)$ 中 Lipschitz 条件,$M_1 \subset M_0$ 为特定可积性指数集。则存在仅依赖于 $N$(数据上界参数)与 $T$ 的常数 $K(N)$,使得弱解差满足:

$$ \begin{aligned} &|\hat\eta - \eta|_{C(0,T;L^2(\Omega))}

  • |\hat u - u|_{L^2(Q)\cap L^\infty(0,T;H^{-1;m})}
  • |\hat\theta - \theta|_{L^2(Q)}
  • |\hat x_e - x_e|{L^{q_e,\infty}(Q)} \ &\quad + |I_t(\hat\sigma - \sigma)|{C(0,T;L^2(\Omega))} \ &\leq K(N)\Big( |\hat\eta^0 - \eta^0|_{L^2(\Omega)}
  • |\hat u^0 - u^0|_{H^{-1;m}}
  • |\widehat{E}^0 - E^0|{H^{-1;3}} \ &\quad + |\hat\beta_e|{L^{q_e}(\Omega)}
  • |\hat u_b - u_b|_{W^{1,1}(0,T)}
  • |\hat p_b - p_b|_{W^{1,1}(0,T)}
  • |\hat\pi_b - \pi_b|{L^1(0,T)} \ &\quad + \inf{(q,r)\in M_0}|\beta_1|_{L^{q,r}(Q)}
  • \inf_{(q,r)\in M_1}|I^{\langle 3\rangle}\beta_2|_{L^{q,r}(Q)}
  • |I_t I^{\langle 1\rangle}\beta_2|{L^{q_e,\infty}(Q)} \ &\quad + |\langle \beta_2\rangle\Omega|_{L^1(0,T)}
  • |\gamma|{[V{2;3^}(Q)]^}
  • |g_1|_{L^1(Q)}
  • |I^{\langle m\rangle}g_2|{L^2(Q)} \ &\quad + \delta{m3}|\langle g_2\rangle_\Omega|_{L^1(0,T)}
  • | \hat f[x_e] - f[x_e] |_{[H^{2,1;\hat\rho,m}(Q)]^*} \Big). \end{aligned} $$

该定理首次在统一框架下建立了对三类边界条件均适用的 Lipschitz 连续性,并将控制范数降至 $L^2(\Omega)\times H^{-1;m}(\Omega)\times H^{-1;3}(\Omega)$ —— 远弱于以往要求的 $L^\infty\times L^2\times L^q$。其推论包括:

  • 均匀化误差估计:当初始数据与源项含周期 $\varepsilon$ 的快速振荡时,原解与 Bakhvalov–Eglit 双尺度均匀化解之差为 $O(\varepsilon)$(在 $L^2(Q)$ 等弱范数下);
  • Hölder 型稳定性:在附加数据正则性下,可提升至更强范数(如 $C([0,T];L^2)$)下的 Hölder 连续性。

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

a) 核心分析学工具

论文构建了一套分层线性化与对偶估计技术,其支柱为三类关键工具:

  1. 加权对偶空间理论:引入空间 $[V_{2;\kappa,m^*}(Q)]^*$ 与 $[H^{2,1;\kappa,m}(Q)]^*$,其范数定义为

    $$ \|w\|_{[V_{2;\kappa,m^*}(Q)]^*} = \sup_{\substack{\varphi \in V_{2;\kappa,m^*}(Q) \\ \|\varphi\|_{V_2(Q)}=1}} \int_Q w \varphi \,dxdt, $$

    并建立与 $L^{q,r}(Q)$ 的嵌入不等式(见原文(4))及积分恒等式(如(7)):

    $$ \|w\|_{[V_{2;\kappa,m^*}(Q)]^*} \leq \|I^{\langle m\rangle}{}^* w\|_{L^2(Q)} \quad (m=2,3), $$

    使对 $L^1$ 级扰动项(如 $g,f$)的控制成为可能。

  2. 原始算子与平均算子体系:定义 $Iy(x)=\int_0^x y$, $I^{\langle 1\rangle}y = Iy - \langle Iy\rangle_\Omega$, $I^{\langle 3\rangle}y = I(y-\langle y\rangle_\Omega)$ 等,利用其对偶关系(如 $(I^{\langle 3\rangle})^* = -I^{\langle 1\rangle}$)将边界条件与平均效应显式分离,为处理 $m=3$(双压力边界)情形提供代数基础。

  3. Gronwall–Bellman 型迭代引理(Lemma 1):针对形如 $\|y\|_{L^r(0,t)} \leq N\|y\|_{L^{r_1}(0,t)} + C$ 的积分不等式,导出 $\|y\|_{L^r(0,T)} \leq K(N,T)C$,成为连接局部能量估计与全局 Lipschitz 常数的关键桥梁。

b) 技术性难题的处理技巧

  • 非线性应力项 $\sigma = \nu\rho D_x u - p$ 的线性化:通过构造扰动差 $\hat\sigma - \sigma$ 并分解为 $\nu(\hat\rho - \rho)D_x u + \nu\hat\rho(D_x\hat u - D_x u) - (p - \hat p)$,再结合 $\eta$ 的 $L^\infty$ 有界性与对数变换 $\nu\ln\eta = \nu\ln\eta^0 + I_t\sigma + I_tp$(见(43)),将非线性耦合转化为可逐项估计的线性结构。

  • 内能方程弱解定义的降阶处理:放弃传统 $V_2(Q)$-框架(需 $\sigma D_x u \in L^1$ 导致 $\varphi \in L^\infty$),转而采用 $L^2(Q)$-解定义(见(32)),测试函数取 $H^{2,1;\hat\rho,m}(Q)$,从而容纳更广的源项正则性,并启用 $H^{2,1}$ 对偶范数控制 $f$。

  • 三类边界条件的统一处理:通过引入通用算子 $I^{\langle m\rangle}$ 与边界修正项 $\sigma_\Gamma$(见(48)),将所有 $m$ 型解表示为 $I^{\langle m\rangle}u + \text{boundary terms}$,使估计公式(42)对 $m=1,2,3$ 具有形式一致性。

  • 快速振荡数据的均匀化耦合:将 $\varepsilon$-振荡初值 $\eta^0_\varepsilon(x) = \eta^0(x,x/\varepsilon)$ 视为扰动 $\beta,\gamma$,利用 Theorem 1 的 Lipschitz 常数 $K(N)$ 与 $\varepsilon$-尺度下扰动范数 $\|\beta\|_{L^{q,r}} = O(\varepsilon)$ 的乘积关系,直接导出 $O(\varepsilon)$ 误差界,避免传统双尺度展开的复杂渐近匹配。


与之前工作的比较

a) 突破与改进

本文相较 Kazhikhov(1980s)、Jiang–Zhang(2000s)、Zlotnik 自身前期工作(如 [3,44]),实现了三重实质性跃升:

  1. 数据正则性要求的系统性降低:

    • 文献 [3, Theorem 3(a)] 要求 $u^0 \in L^\infty(\Omega), \theta^0 \in L^2(\Omega)$,而本文在 Theorem 1 中仅需 $u^0 \in H^{-1;m}(\Omega), \theta^0 \in L^1(\Omega)$(见 $(\widehat{C}_1)$),且 $H^{-1;m}$ 范数可由 $L^2$ 初值 $\eta^0$ 与总能量 $E^0$ 共同控制;
    • 边界数据从 $W^{1,\infty}(0,T)$ 削弱至 $W^{1,1}(0,T)$,源项 $g,f$ 从 $L^2(Q)$ 放宽至 $[V_{2;\kappa,m^*}]^*$ 或 $[H^{2,1}]^*$,覆盖 $L^1$ 甚至部分测度级扰动。
  2. 结果强度的质的提升:

    • 以往 Lipschitz 结果(如 [44])仅对 $m=3$ 成立,且依赖 $q_e=2$;本文 Theorem 1 首次对 $m=1,2,3$ 统一成立,并允许 $q_e \in [2,\infty]$,为 Hölder 推论(Remark 1.2)提供自由度;
    • 误差估计从文献 [7] 的 $O(\varepsilon)$(仅 barotropic 情形)推广至完整热力学模型,且明确给出非平凡温度平均 $\langle \theta^0\rangle_{\text{fast}}$ 在 Bakhvalov–Eglit 方程中的作用(见摘要)。
  3. 方法论的范式拓展:

    • 将“连续依赖性”从解存在性证明的副产品(如 [14] 中用于唯一性)提升为独立主定理,并发展出可直接服务于数值分析(如多尺度算法验证)与反问题(如参数识别稳定性)的定量工具;
    • 通过对 $I^{\langle m\rangle}$ 算子系的深度挖掘,将边界类型差异转化为同一估计框架内的指标切换,为高维推广(如柱对称情形)预留了清晰路径。

综上,本文不仅填补了弱范数下 Lipschitz 连续性的理论空白,更构建了一套面向实际工程粗糙输入的鲁棒分析范式,标志着一维可压流体力学适定性理论从“理想光滑世界”向“真实不规则世界”的关键迈进。


3. 平面定常Navier–Stokes方程自相似解的局部渐近轮廓

Local profiles of self-similar solutions of the planar stationary Navier–Stokes equations

作者: Ming Li, Linyu Peng, Ping Zhang, Xin Zhang

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03341v1

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

详细分析

中文标题: 平面定常Navier–Stokes方程自相似解的局部渐近轮廓


研究对象和背景

本文研究的是二维定常不可压缩Navier–Stokes方程(Stationary Incompressible Navier–Stokes equations, SNS)在以原点为顶点的锥形区域 $\Omega \subset \mathbb{R}^2$ 上的自相似解(self-similar solutions),特别关注其在 $\Omega$ 内部的局部渐近轮廓(local profiles),即解在原点附近、不依赖于边界条件时的结构特征。该问题属于偏微分方程(PDE)中关于奇点结构与缩放不变性的经典课题,与流体力学中著名的 Jeffery–Hamel 流密切相关。

论文考虑的控制方程为无外力情形下的 SNS 方程:

$$ \begin{cases} \operatorname{div}(U \otimes U) - \Delta U + \nabla p = 0, \\ \operatorname{div} U = 0, \end{cases} \quad \text{in } \Omega, $$

其中未知量 $U = (u(x,y), v(x,y))^\top$ 表示速度场,$p = p(x,y)$ 是压力函数。符号 $a \otimes b$ 表示 $2\times2$ 矩阵 $[a\;b]$;对任一 $2\times2$ 矩阵 $A = [A_{jk}]$,向量 $\operatorname{div} A$ 的第 $j$ 分量定义为 $\sum_{k=1}^{2}\partial_k A_{jk}$。该系统在 $\mathbb{R}^2\setminus\{(0,0)\}$ 上定义,且要求解在开锥域 $\Omega$(含原点为顶点)上光滑:$(u,v,p)\in C^\infty(\Omega)^3$。

关键几何设定是:$\Omega$ 是一个以原点为顶点的开锥,即存在角度区间 $I \subset (-\pi,\pi]$,使得

$$ \Omega = \left\{ (r\cos\theta, r\sin\theta) \mid r > 0,\; \theta \in I \right\}. $$


论文未对 $\partial\Omega$ 施加任何边界条件(如无滑移或自由滑移),亦不要求 $\Omega$ 关于坐标轴对称——这显著区别于传统 Jeffery–Hamel 问题中两刚性平板夹角流动的经典设定(如 $\Omega = \{ (x,y) \mid |\theta| < \alpha \}$ 且 $u=v=0$ 在 $\theta = \pm \alpha$)。这种“无边界约束的局部分析”使得研究焦点完全转向由缩放对称性所决定的内在结构。

该系统具有基本的缩放不变性:若 $(U,p)$ 是解,则对任意 $\lambda>0$,函数

$$ U_\lambda(x,y) := \lambda U(\lambda x, \lambda y), \quad p_\lambda(x,y) := \lambda^2 p(\lambda x, \lambda y) $$


也是解。由此引出自相似解的定义:满足 $(U,p) = (U_\lambda,p_\lambda)$ 对所有 $\lambda>0$ 成立。这一条件等价于解在极坐标 $(r,\theta)$ 下具有如下形式:

$$ U = \frac{f(\theta)}{r} e_r + \frac{g(\theta)}{r} e_\theta, \quad p = \frac{p(\theta)}{r^2}, $$

其中 $e_r = (\cos\theta,\sin\theta)^\top$, $e_\theta = (-\sin\theta,\cos\theta)^\top$。代入原方程后可导出关于 $f,g,p$ 的常微分方程组(ODE system):

$$ \begin{cases} (p - 2f)' = 0, \\ -f'' + g f' - f^2 - g^2 - 2p = 0, \\ g' = 0. \end{cases} $$

由第三式知 $g \equiv \text{const}$,进而第一式推出 $p = 2f + C$;代入第二式得单个二阶 ODE:

$$ f'' - C f' + f^2 + C^2 + 4f = 0, $$

此即经典的 Jeffery–Hamel ODE。当 $C = 0$(即 $g \equiv 0$)时,对应径向流(radial flow),此时 $U \parallel e_r$,满足 $xv = yu$;当 $C \neq 0$ 时,则为非径向流,满足 $xv = yu + C_0$。

论文强调,其研究区域 $\Omega$ 可为任意开锥,开口角 $\alpha(\Omega)$ 可小于或大于 $\pi$(即“小锥”或“大锥”),且不预设解在 $\partial\Omega$ 上的行为。尤其重要的是,作者明确排除了对解整体延拓至 $\mathbb{R}^2$ 的要求,仅关注其在锥内局部的解析结构——这使问题本质上成为奇点附近的渐近分析,而非全局边值问题。

函数空间设定为 $C^\infty(\Omega)$,即要求解无限次可微,从而保证所有推导中的微分运算合法,并支持后续通过 ODE 积分、椭圆函数反演等精细构造。值得注意的是,由于锥域可能包含正半轴或负半轴(如 $\Omega \cap \{x=0\} = \emptyset$),作者采用直角坐标系下的变量替换 $z = y/x$,将 PDE 约化为关于 $z$ 的常微分方程,避免了极坐标在 $\theta = \pm\pi/2$ 处的奇异性。该约化依赖于Lie 点对称性理论,特别是由无穷小生成元 $v_3 = x\partial_x + y\partial_y - u\partial_u - v\partial_v - 2p\partial_p$ 所生成的缩放对称。

综上,本文的研究对象是:在任意顶点在原点的开锥域 $\Omega$ 上,无边界条件约束、无限可微的定常不可压缩 Navier–Stokes 方程的自相似解,核心目标是刻画其在原点附近的局部渐近形态,按径向与非径向两类分别给出显式表达式。


主要定理或主要结果

论文建立了三类核心结果,分别对应径向小锥、径向大锥及非径向情形:

Theorem 1.3(小锥径向解的局部轮廓):
设 $\Omega \subset \mathbb{R}^2$ 为以原点为顶点的开锥,且满足 $\Omega \cap \{x = 0\} = \emptyset$;设 $(u,v,p) \in C^\infty(\Omega)^3$ 是 (1.1) 的非平凡局部自相似解,并满足径向条件 $xv = yu$。则 $(u,v,p)$ 具有形式

$$ u(x,y) = \frac{x}{x^2+y^2}\kappa(x,y), \quad v(x,y) = \frac{y}{x^2+y^2}\kappa(x,y), \quad p(x,y) = \frac{2}{x^2+y^2}\kappa(x,y) + \frac{C_1}{x^2+y^2}, $$


其中 $\kappa = \kappa(x,y)$ 由表格 1 给出(对应不同参数组合下 $h(z)=(z^2+1)u(x,y)$ 的显式积分表达式),$C_1\in\mathbb{R}$ 为常数。

Theorem 1.5(大锥径向解的局部轮廓):
设 $\Omega$ 为顶点在原点、开口角 $\alpha(\Omega) \in (\pi,2\pi]$ 的开锥;$(u,v,p) \in C^\infty(\Omega)^3$ 满足径向条件。则其形式同上,但 $\kappa(x,y)$ 是 Theorem 1.3 中某一支 $\kappa$ 的实解析延拓;进一步,若该解可光滑延拓至整个 $\mathbb{R}^2$,则必为经典径向 Jeffery–Hamel 解(Definition 1.1)。

Theorem 1.6(非径向解的局部结构):
设 $\Omega$ 满足 $\Omega \cap \{x = 0\} = \emptyset$,$(u,v,p) \in C^\infty(\Omega)^3$ 是满足非径向条件 $xv = yu + C_0$($C_0 \neq 0$)的非平凡解。则存在两类解:
(i) 经典线性 Jeffery–Hamel 解 $(u_L,v_L,p_L)$;
(ii) 形如

$$ u = u_L + \frac{x}{x^2+y^2}H(\theta),\quad v = v_L + \frac{y}{x^2+y^2}H(\theta),\quad p = p_L + \frac{2}{x^2+y^2}H(\theta), $$


其中 $\theta = \arctan(y/x)$,且 $H = H(\theta)$ 满足 Lié nard 方程:

$$ H'' - C_0 H' + H^2 + (2\tilde{C}_1 + 4)H = 0. $$


在可积条件下,该方程的解可用 Weierstrass 椭圆函数 $\wp$ 显式表示(见 (4.16));退化情形下则回归初等函数解((4.17))。


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

a) 核心分析学工具

论文的核心方法是对称性约化与常微分方程定性分析。首先,利用 Lie 群理论识别系统 (2.1) 的全部 Lie 点对称,提取缩放对称生成元 $v_3$,构造其不变量:$z = y/x$,并定义缩放不变变量

$$ \tilde{u}(z) = x u(x,y),\quad \tilde{v}(z) = x v(x,y),\quad \tilde{p}(z) = x^2 p(x,y). $$


此变换将二维 PDE 完全降维为关于 $z$ 的三阶耦合 ODE 系统 (2.5),再结合连续性方程消去 $v$,最终导出单个四阶非线性 ODE (2.10)。这是典型的相似约化(similarity reduction)技术,其数学基础是对称群作用下的商流形理论。

对于径向情形($C_0 = 0$),方程 (2.11) 被进一步简化为关于 $h(z) = (z^2+1)\tilde{u}(z)$ 的二阶 ODE (2.13):

$$ (z^2+1)^2 h_{zz} + 2z(z^2+1) h_z + h^2 + 4h + 2C_1 = 0. $$


该方程经乘以 $2h_z$ 后可首次积分,得到一阶守恒律:

$$ (z^2+1)^2 h_z^2 + \frac{2}{3}h^3 + 4h^2 + 4C_1 h + C_2 = 0, $$


即

$$ h_z = \pm \frac{\sqrt{-h^3 - 6h^2 - 6C_1 h - C_2}}{z^2+1}. $$


此即标准的椭圆型一阶 ODE,其解由椭圆积分表出。作者系统运用Jacobi 椭圆函数理论:通过变量代换(如 $\theta = \arccot\sqrt{(\alpha-h)/\beta}$ 或 $\theta = \arcsin\sqrt{(c-h)/(c-b)}$)将积分化为第一类不完全椭圆积分 $F(\phi,k)$,再利用其反函数 Jacobi 幅角函数 $\operatorname{am}(m,k)$ 表达原函数 $h(z)$,最终获得 (2.21) 等显式公式。该过程严格依赖于三次多项式判别式分析($\Delta_{\text{cubic}}$)对根分布的分类(三实根、一实两共轭复根、重根),每种情形对应不同的椭圆模 $k$ 与幅角变换。

b) 技术性难题处理技巧

首要挑战是非线性与奇异性耦合。方程 (2.11) 含有 $(z^2+1)^2 h_{zz}$ 项,在 $z \to \pm\infty$(即 $\theta \to \pm\pi/2$)时系数退化,但作者通过引入 $h(z)$ 巧妙规避了直接处理 $u_{zz}$ 的奇性。其次,非径向情形下导出的 Liénard 方程 (1.13) 一般不可积,作者引入Painlevé–Gambier 可积性判据:当系数满足特定关系((4.10))时,方程可化为 Painlevé II 型,从而允许用 Weierstrass $\wp$ 函数求解。具体地,令 $H = a \wp(z-z_0) + b$,代入并匹配系数,确定 $a,b$ 及不变量 $g_2,g_3$,实现椭圆函数反演构造。此外,为处理大锥情形($\alpha > \pi$),作者采用解析延拓原理:小锥解 $\kappa$ 定义在 $z$ 的某区间上,只要其幂级数收敛半径足够大,即可唯一延拓至更大区域,确保 (1.10) 的合理性——这依赖于复分析中的唯一性定理与圆盘上解析函数延拓理论。


与之前工作的比较

a) 突破与改进

本文相较经典文献(如 Sverák [19], Hamel [8], Jeffery [10])实现了三重实质性推进:

第一,区域设定的彻底解放。 以往 Jeffery–Hamel 理论严格限定于全平面 $\mathbb{R}^2$ 或对称双壁锥域(如 $\theta \in (-\alpha,\alpha)$ 且 $u=v=0$ 在边界),而本文首次系统研究任意开锥 $\Omega$(含 $\alpha > \pi$)上无边界条件的局部解。Theorem 1.3/1.5 表明,局部轮廓完全由缩放对称性与 ODE 相图决定,与边界无关——这揭示了自相似奇点的内在普适性,为奇点分类提供新范式。

第二,非径向解的显式构造突破。 经典工作仅给出 $g \equiv \text{const}$ 的 ODE,但未解出非零 $g$ 时的显式解。本文 Theorem 1.6 首次将非径向问题归结为 Liénard 方程,并在可积条件下导出 Weierstrass $\wp$ 函数表达式((4.16)):

$$ H(\theta) = -2\wp\!\left(\theta - \theta_0; \; g_2 = \frac{4}{3}(2\tilde{C}_1+4)^2, \; g_3 = \frac{8}{27}(2\tilde{C}_1+4)^3 - \frac{C_0^2}{4}\right) - \frac{2\tilde{C}_1+4}{3}. $$


此为首个将 Jeffery–Hamel 非径向解与现代椭圆函数理论深度关联的工作,远超早期仅数值或定性讨论的水平。

第三,技术工具的精细化升级。 相比 [2,7,19] 中对 ODE 的定性分析或特殊参数求解,本文运用完整的三次方程根分布—椭圆积分类型—Jacobi 函数反演三级映射,对 (2.16) 进行穷尽式分类(Case 1–3),覆盖所有实参数组合,并给出统一框架下的显式公式。这种基于判别式的全参数显式求解法,在 Navier–Stokes 自相似解研究中属首创,显著增强了结果的完备性与应用潜力。


4. 三维Navier-Stokes–Vlasov–Fokker–Planck方程组适定弱解的存在性与部分正则性

Existence and partial regularity of suitable weak solutions to the 3D Navier-Stokes-Vlasov-Fokker-Planck equations

作者: Renjun Duan, Fengqiang Shi, Wendong Wang, Jianbo Yu

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03131v1

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

详细分析

中文标题: 三维Navier-Stokes–Vlasov–Fokker–Planck方程组适定弱解的存在性与部分正则性


研究对象和背景

本文研究的是定义在全空间 $\mathbb{R}^+ \times \mathbb{R}^3 \times \mathbb{R}^3$ 上的 三维不可压Navier-Stokes–Vlasov–Fokker–Planck耦合系统(NSVFP), 描述一种典型的两相流动力学模型:粘性不可压缩流体与悬浮于其中的大量微观粒子(或气泡)通过摩擦力项相互作用。该模型兼具宏观连续介质力学与微观粒子统计动力学特征,在多相流、气溶胶传输、生物活性流体(如趋化性-流体耦合)、工业悬浮液模拟等领域具有重要物理意义。

系统的基本方程组在引言第1节中明确给出,即如下初值问题:

$$ \begin{cases} \partial_t u + (u \cdot \nabla_x) u - \Delta_x u + \nabla_x P - \displaystyle\int_{\mathbb{R}^3_v} (v - u) f \, dv = 0, \\ \partial_t f + v \cdot \nabla_x f + \nabla_v \cdot \big[ (u - v) f - \nabla_v f \big] = 0, \\ \nabla_x \cdot u = 0, \end{cases} \tag{1.1} $$

其中:

  • $u = u(t,x) \in \mathbb{R}^3$ 是流体速度场,定义在 $(t,x) \in \mathbb{R}^+ \times \mathbb{R}^3$;
  • $f = f(t,x,v) \geq 0$ 是粒子在相空间 $\mathbb{R}^3_x \times \mathbb{R}^3_v$ 中的概率密度分布函数;
  • $P = P(t,x)$ 是标量压力;
  • 摩擦力项 $\int_{\mathbb{R}^3_v} (v - u) f \, dv$ 表征流体与粒子间的动量交换,是系统非局部耦合的核心来源;
  • Vlasov-Fokker-Planck方程第二式包含对流项 $v \cdot \nabla_x f$、漂移项 $\nabla_v \cdot[(u - v)f]$ 和扩散项 $\Delta_v f$,体现粒子受流体拖曳与热噪声共同作用下的演化。

研究区域为三维全空间 $\mathbb{R}^3_x$,无边界条件;初始数据满足:

$$ \begin{cases} u_{\text{in}} \in L^2_\sigma(\mathbb{R}^3), \\ f_{\text{in}} \in L^\infty \cap L^1(\mathbb{R}^6),\quad f_{\text{in}} \ge 0, \\ \big( |x|^2 + |v|^\kappa + |\log f_{\text{in}}| \big) f_{\text{in}} \in L^1(\mathbb{R}^6), \quad \kappa \ge 3. \end{cases} \tag{1.2} $$

此处 $L^2_\sigma(\mathbb{R}^3)$ 表示散度为零的 $L^2$ 向量场空间;$L^1(\mathbb{R}^6)$ 和 $L^\infty(\mathbb{R}^6)$ 分别对应相空间上的可积性与有界性;加权对数矩 $\int f_{\text{in}} \log f_{\text{in}} \, dxdv$ 控制熵结构,保障自由能守恒;而高阶速度矩 $|v|^\kappa f_{\text{in}} \in L^1$ 对后续控制摩擦力项的可积性至关重要。

特别地,论文强调其解空间需满足强加权可积性要求:除基本能量空间外,还需 $|v|^\kappa f \in L^\infty_{\text{loc}}([0,\infty); L^1(\mathbb{R}^6))$($\kappa \ge 3$),以及更严格的 $|v|^{15+\kappa} f_{\text{in}} \in L^1(\mathbb{R}^6)$(见Theorem 1.5假设(1.6))。这一技术性条件旨在确保摩擦力项 $\int (v-u)f \, dv$ 属于 $L^q_{t,x}$,且 $q > 5/2$,从而兼容Caffarelli-Kohn-Nirenberg型部分正则性分析框架。

此外,解需满足三个不同形式的能量不等式:全局能量不等式、第一类局部能量不等式(含测试函数 $\psi(x), \phi(x,v)$)、第二类局部能量不等式(仅含 $\psi(x)$)。这些不等式并非标准弱解所自动满足,而是定义“suitable weak solution”的关键判据(Definition 1.1),其构造与验证构成全文核心难点。


主要定理或主要结果

本文建立了以下三项主成果:

Theorem 1.2(适定弱解的存在性)
设初始数据 $(u_{\text{in}}, f_{\text{in}})$ 满足条件 (1.2),即

$$ u_{\text{in}} \in L^2_\sigma(\mathbb{R}^3),\quad f_{\text{in}} \in L^\infty \cap L^1(\mathbb{R}^6),\; f_{\text{in}} \ge 0,\quad \big(|x|^2 + |v|^\kappa + |\log f_{\text{in}}|\big) f_{\text{in}} \in L^1(\mathbb{R}^6),\; \kappa \ge 3, $$


则Cauchy问题 (1.1) 存在全局 suitable weak solution $(u,f)$,满足 Definition 1.1 中全部六条性质(包括能量估计、三类能量不等式等)。

Theorem 1.5(速度场的部分正则性)
在Theorem 1.2基础上,若进一步假设

$$ |v|^{15+\kappa} f_{\text{in}} \in L^1(\mathbb{R}^6),\quad \kappa > 0, $$


则存在 $\varepsilon_0 > 0$,使得对任意 $z_0 = (t_0,x_0) \in (0,T] \times \mathbb{R}^3$,若

$$ \limsup_{r \to 0^+} \frac{1}{r} \iint_{Q(z_0,r)} |\nabla_x u|^2 \, dxdt < \varepsilon_0, $$


其中 $Q(z_0,r) = (t_0 - r^2, t_0) \times B(x_0,r)$ 为抛物柱,则 $u$ 在 $z_0$ 处正则(即存在邻域内 $u \in L^\infty$);进而,奇点集 $S_{u,\text{sing}} = \{z : u \text{ 不正则于 } z\}$ 满足

$$ \mathcal{P}^1(S_{u,\text{sing}}) = 0, $$


即其抛物 Hausdorff 维数为零(等价于一维抛物测度为零)。

Theorem 1.6(分布函数的Hölder连续性)
在Theorem 1.2与(1.6)条件下,令 $S_{u,\text{reg}} = \{z \in (0,T] \times \Omega_x : u \text{ 正则于 } z\}$,则对任意 $z_0 = (t_0,x_0,v_0) \in S_{u,\text{reg}} \times \mathbb{R}^3_v$,解 $f(t,x,v)$ 在内部区域 $Q_{\text{int}} = (t_0 - r^2, t_0) \times B(x_0,r) \times B(v_0,r)$ 上满足

$$ \|f\|_{C^\alpha(Q_{\text{int}})} \le C \big( \|f\|_{L^2(Q_{\text{ext}})} + \|f\|_{L^\infty(Q_{\text{ext}})} \big), $$


其中 $Q_{\text{ext}} = (t_0 - r_1^2, t_0) \times B(x_0,r_1) \times B(v_0,r_1)$,$r_1 > r > 0$,且 Hölder指数 $\alpha = \alpha(\Lambda) \in (0,1)$ 依赖于某参数 $\Lambda$。


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

a) 核心分析学工具

  1. 正则化与紧性方法:
    全文采用经典的双参数正则化策略:引入 mollifier $\theta_\varepsilon$ 正则化速度场与耦合力项,并用截断函数 $\gamma_\delta(v)$ 限制粒子速度范围,构造近似系统 (1.3)。该系统具备光滑性,可借助Banach不动点定理(Section 2.1)及Stokes/Fokker-Planck基本解理论建立局部强解,再通过先验估计延拓为全局强解(Theorem 1.4)。

  2. Tao 的 $L^p$ 分解与 DiPerna-Lions 紧性框架:
    为克服 $f_{\varepsilon,\delta}$ 在全空间上缺乏紧性的困难,作者创新性地结合 Tao 的 Littlewood-Paley 分解(见 Remark 1.3 引用 [62])控制 $\|u\|_{L^1_t L^p_x}$ 范数,从而导出 $|v|^\kappa f$ 的一致 $L^\infty_t L^1_{x,v}$ 界(Lemma 3.6–3.7)。继而运用 DiPerna-Lions 的输运方程紧性引理([25]),在加权空间中获得 $f_{\varepsilon,\delta} \to f$、$|v|^\kappa f_{\varepsilon,\delta} \to |v|^\kappa f$ 及 $f_{\varepsilon,\delta} \log f_{\varepsilon,\delta} \to f \log f$ 的强收敛性(Sections 5.2–5.5)。这是验证局部能量不等式中非线性项收敛性的基石。

  3. Aubin-Lions 引理与 Lions 弱收敛法:
    对速度场 $u_{\varepsilon,\delta}$,利用 Aubin-Lions 引理(Section 5.4)从 $L^2_t H^1_x \cap L^\infty_t L^2_x$ 嵌入得 $u_{\varepsilon,\delta} \to u$ 在 $L^2_{t,x}$ 中强收敛;对非局部摩擦力项 $\int (v-u)f \, dv$,则采用 Lions 的弱收敛技巧([47])处理其在分布意义下的极限(Section 5.1)。

  4. 三类能量不等式的构造技术:

    • 全局能量不等式(Section 4.1):直接对 (1.1) 做 $L^2$ 内积并积分,结合熵结构 $\int f \log f$ 得到总能量耗散;
    • 第一类局部能量不等式(Section 4.2):以时空截断函数 $\psi(x), \phi(x,v)$ 为权重,对 $|u|^2$ 和 $\left(\frac{|v|^2}{2} + \log f\right)f$ 分别加权测试,需精细处理交叉项如 $\int fv \cdot u \phi$、压力项 $(P - \bar{P}) u \cdot \nabla_x \psi$;
    • 第二类局部能量不等式(Section 4.3):仅用 $\psi(x)$ 加权,将 $f$-相关项整体视为源项,适用于更高阶矩控制。

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

  • 非局部摩擦力项的收敛性:
    关键难点在于证明 $\int (v - u_{\varepsilon,\delta}) f_{\varepsilon,\delta} \, dv \rightharpoonup \int (v - u) f \, dv$。作者通过加权强收敛 $|v|^\kappa f_{\varepsilon,\delta} \to |v|^\kappa f$($\kappa > 3$)与 $u_{\varepsilon,\delta} \to u$ 强收敛的组合,结合 Hölder 不等式与 Lemma 2.2 的矩估计(式 (2.3)(2.4)),将该非局部项嵌入 $L^q_{t,x}$($q > 5/2$),从而规避弱收敛下乘积极限失效问题。

  • 对数非线性项 $f \log f$ 的处理:
    利用凸性与 Jensen 不等式保证熵项符号;通过加权 $L^1$ 强收敛(Section 5.5)及 Vitali 收敛定理,验证 $\iint f_{\varepsilon,\delta} \log f_{\varepsilon,\delta} \phi \to \iint f \log f \phi$。

  • 全空间高阶矩控制:
    针对 $\mathbb{R}^3_x$ 无界性导致的矩估计失效,作者借助 Tao 分解获得 $\|u\|_{L^1_t L^p_x}$ 控制($p > 3$),再由 Lemma 2.2 的插值不等式 (2.4) 将 $|v|^\kappa f$ 的 $L^1_{x,v}$ 范数与 $L^\infty_{x,v}$ 范数耦合,实现一致有界性。


与之前工作的比较

a) 突破与改进

本文首次在三维全空间中为 NSVFP 系统建立了具备Caffarelli-Kohn-Nirenberg 型部分正则性的 suitable weak solution 理论,相较既有工作有三重实质性推进:

  1. 解的概念升级:
    此前文献(如 [13,14,33,50])多关注 Leray-Hopf 类弱解或经典解的存在性,未系统构建满足多类局部能量不等式的 suitable 解。本文明确定义并验证了三类能量不等式,为奇点分析提供严格基础。

  2. 初始数据假设的优化:
    相比 Chae-Kang-Lee [14] 要求 $f_{\text{in}} \in H^1_{x,v} \cap L^1_{x,v}(|v|^2)$,本文仅需 $f_{\text{in}} \in L^\infty \cap L^1$ 及加权对数矩,显著降低正则性要求;且对速度矩仅需 $|v|^\kappa f_{\text{in}} \in L^1$($\kappa \ge 3$),而非早期工作中常见的 $|v|^2 f_{\text{in}} \in L^\infty$ 或更高阶 Sobolev 假设。

  3. 正则性结论的强化:
    Theorem 1.5 给出奇点集的抛物 Hausdorff 维数为零,优于 [17] 中 chemotaxis-NS 系统的 $\mathcal{P}^1$-measure 零结果;Theorem 1.6 首次建立 $f$ 在 $u$ 正则点处的 Hölder 连续性,揭示粒子分布对流场正则性的敏感依赖,此为 kinetic-fluid 模型中全新现象。

综上,本文不仅填补了三维 NSVFP 系统 suitable 解理论的空白,更通过融合现代调和分析、输运方程紧性与抛物正则性理论,为复杂多尺度耦合系统的数学分析树立了新范式。


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

1. 冷等离子体方程的内自由边界问题

Internal free boundary problem for cold plasma equations

作者: Lidia Gargyants, Anna Konovalova, Olga Rozanova

类别: math.AP, math-ph

发布日期: 2026-02-03

arXiv ID: 2602.03463v1

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

摘要翻译

中文标题: 冷等离子体方程的内自由边界问题

摘要翻译: 针对描述电子在静止离子场中运动的冷等离子体方程组,本文研究了在两种介质之间不可穿透界面上提出的黎曼问题。这两种介质的区别在于其恒定离子场强度的不同。介质间的界面被假定为自由边界,其位置由广义Rankine–Hugoniot条件及稳定性条件共同确定;该稳定性条件即拉格朗日质点轨迹在界面上发生交汇。


2. 通过半经典孤子系综研究中长波方程的小色散极限

The Small Dispersion Limit of the Intermediate Long Wave Equation via Semiclassical Soliton Ensembles

作者: Matthew Dominique Mitchell

类别: math.AP, math-ph

发布日期: 2026-02-03

arXiv ID: 2602.03063v1

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

摘要翻译

中文标题: 通过半经典孤子系综研究中长波方程的小色散极限

摘要翻译: 本文研究中长波(ILW)方程的小色散极限,特别针对一类性质良好的初值条件 $u_0$,其对应解中孤子个数随色散参数趋于零而无界增长。首先,我们对ILW方程的直接散射问题开展形式化的WKB型渐近分析,由此构造出特征值与归一化常数的渐近近似表达式;继而基于该近似,定义一组修正的散射数据,并严格分析其对应的逆散射问题。主要结果包括:证明在 $t = 0$ 时刻,解在 $L^2$ 范数意义下收敛于原始初值 $u_0$;而在 $0 < t < t_\mathrm{c}$ 区间内,解则在 $L^2$ 范数意义下收敛于无粘性Burgers方程的相应解,其中 $t_\mathrm{c}$ 为梯度破裂(gradient catastrophe)发生的时间。


3. 稳定旋转星–行星系统的存在性

Existence for Stable Rotating Star-Planet Systems

作者: Hangsheng Chen

类别: math.AP, math-ph

发布日期: 2026-02-02

arXiv ID: 2602.02761v1

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

摘要翻译

中文标题: 稳定旋转星–行星系统的存在性

摘要翻译: 本文研究稳定、均匀旋转的星–行星系统(即质量比足够小)的存在性与性质,其物理模型由欧拉–泊松方程组描述。沿用McCann关于双星系统所建立的变分框架[参见文献《McC06》],我们采用变分方法,在状态方程 $P(\rho) = K\rho^\gamma$ 及质量比 $m$ 足够小(对应于星–行星构型)的假设下,证明了关于Wasserstein $L^\infty$ 度量的局部能量极小元的存在性。此类极小元对应于欧拉–泊松方程组的解。我们分两种情形进行讨论:当 $\gamma > 2$ 时,我们不仅证明极小元的存在性,还通过尺度分析进一步表明——其支撑集的半径(更准确地说,支撑集的直径上界)趋于零;当 $\frac{3}{2} < \gamma \leq 2$ 时,我们估计了半径(即支撑集尺度)可能的膨胀速率之上界,并证实该情形下前述存在性结论依然成立。最后,我们给出了极小元支撑集各连通分支之间距离的下界估计,并就支撑集连通分支的个数提出一个猜想。


不相关论文(基本信息)

1. 平面上临界质量下全抛物型Keller–Segel系统的整体解存在性

Global existence for the fully parabolic Keller–Segel system with critical mass on the plane

作者: Tatsuya Hosono

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03768v1

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

摘要: We study the global existence of solutions to the Cauchy problem for the two-dimensional fully parabolic Keller–Segel system at the critical mass. It is known that global-in-time existence holds for initial data with critical mass under radial symmetry or suitable moment conditions, whereas the behavior of general solutions in the critical regime remains delicate. In this paper, we establish global-in-time existence for general initial data with critical mass, without imposing any symmetry or moment assumptions. The proof relies on the construction of a reconstructed Lyapunov functional, combined with refined regularity estimates for the associated dissipative terms, which enable us to control the solution dynamics in the critical regime.


2. 类时曲线:同伦性与决定域

Timelike curves: homotopies and domain of determinacy

作者: Jérôme Le Rousseau, Jeffrey B. Rauch

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03715v1

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

摘要: This paper studies domains of determination of linear strictly hyperbolic second order operators $P$. For an open set $\mathcal O$, a set $Z$ is a domain of determination when the values of solutions of the differential equation $Pu=0$ are determined on $Z$ by their values in $\mathcal O$. Fritz John’s global Hölmgren theorem implies that points that can be reached by deformations of noncharacteristic hypersufaces with initial surface and boundaries in $\mathcal O$ belong to a domain of determination provided that local uniqueness holds at noncharacteristic surfaces. Using spacelike hypersurfaces yields sharp finite speed results whose domains of determination are described in terms of influence curves that never exceed the local speed of propagation. This paper studies deformations of noncharacteristic nonspacelike hypersurfaces. We prove that points reachable by (repeated) deformations by noncharacteristic nonspacelike hypersurfaces coincide exactly with the set of points reachable by (repeated) homotopies of timelike arcs whose initial curves and endpoints belong to $\mathcal O$. When the set $\mathcal O$ is a small neighborhood of a forward timelike arc connecting $a$ to $b$, a natural candidate for $Z$ is the intesection of the future of $a$ with the past of $b$. This candidate is exact for D’Alembert’s equation. We prove that it is also exact when $a,b$ are points close together on a fixed timelike arc. The timelike homotopy criterion fuels the construction of surprising examples for which the domain of determination is strictly larger (resp. strictly smaller) than the future-intersect-past candidate.


3. 关于几乎周期椭圆算子的临界性与主特征值

On the criticality and the principal eigenvalue of almost periodic elliptic operators

作者: Luca Rossi

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03644v1

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

摘要: We review the notion and the properties of the generalised \pe\ for elliptic operators in unbounded domains, and we relate it with the criticality theory. We focus on operators with almost periodic coefficients. We present a Liouville-type result in dimension $N\leq2$. Next, we show with a counter-example that criticality is not equivalent to the existence of an almost periodic principal eigenvalue, even for self-adjoint operators. Finally, we exhibit an almost periodic operator which is subcritical but which admits a critical limit operator. This is a manifestation of the instability character of the criticality property in the almost periodic setting.


4. 具有孤立奇点的平面Hartree方程解的渐近行为

Asymptotic behavior of solutions to a planar Hartree equation with isolated singularities

作者: Tao Feng, Minbo Yang, Xianmei Zhou

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03559v1

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

摘要: In this paper we investigate the isolated singularities of the Hartree type equation \begin{equation*} -Δu (x)= \left(\frac{1}{|x|^α}e^u\right)e^{u(x)}\quad \text{in } B_{1}\setminus{0} , \end{equation} where $α>0$, $\displaystyle \frac{1}{|x|^α}*e^u\triangleq\int_{B_{1} \setminus \{0\}}\frac{e^u(y)}{|x-y|^α}dy$, and the punctured ball $B_{1}\setminus\{0\}\subset \mathbb{R}^2$. Under the finite total curvature condition, by establishing a representation formula for singular solutions, we obtain the asymptotic behavior of the solutions near the origin. We also extend this asymptotic behavior results to the case with a general non-negative coefficient $K(x)$, and to the higher-order Hartree-type equations in any dimension $n \geq 3$.


5. 关于高频散射中三类模型的奇异伽辽金离散化方法

On singular Galerkin discretizations for three models in high-frequency scattering

作者: T. Chaumont-Frelet, S. Sauter

类别: math.NA, math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03428v1

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

摘要: We consider three common mathematical models for time-harmonic high frequency scattering: the Helmholtz equation in two and three spatial dimensions, a transverse magnetic problem in two dimensions, and Maxwell’s equation in three dimensions with dissipative boundary conditions such that the continuous problem is well posed. In this paper, we construct meshes for popular (low order) Galerkin finite element discretizations such that the discrete system matrix becomes singular and the discrete problem is not well posed. This implies that a condition “the finite element space has to be sufficiently rich” in the form of a resolution condition - typically imposed for discrete well-posedness - is not an artifact from the proof by a compact perturbation argument but necessary for discrete stability of the Galerkin discretization.


6. 关于双曲空间上涉及分数阶GJMS算子的Poincaré–Sobolev型不等式

On Poincaré-Sobolev level involving fractional GJMS operators on hyperbolic space

作者: Huyuan Chen, Rui Chen

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03299v1

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

摘要: This paper is devoted to a qualitative analysis of the Poincaré–Sobolev level associated with the fractional GJMS operators $\mathcal{P}_s$ $\bigl(s\in(0,\tfrac n2)\setminus\mathbb N\bigr)$ on the hyperbolic space $\mathbb H^n$. In contrast to the integer-order case, when $s\notin\mathbb N$ the operator $\mathcal{P}_s$ does not enjoy the conformal covariance that allows one, in the upper half-space or ball model, to relate it to the Euclidean fractional Laplacian $(-Δ)^s$; this link is crucial for importing Euclidean theory. We therefore introduce $\widetilde{\mathcal{P}}_s$ ($s>0$), which is conformally related to $(-Δ)^s$. Our purpose in the paper is to analyze the monotonicity, attainability, and strict-gap regions of the Poincaré–Sobolev levels associated with $\mathcal{P}_s$ and with $\widetilde{\mathcal{P}}_s$ with $s\in(0,\tfrac n2)\setminus\mathbb N$. First, we reinterpret the Brezis–Nirenberg problem through the lens of Poincaré–Sobolev levels, connecting earlier results for the Euclidean Laplacian and for operators $\mathcal{P}_k$ on $\mathbb H^n$ with integer $k\in(0,\tfrac n2)$. We then establish new, explicit lower bounds for the Hardy term in fractional Hardy–Sobolev–Maz’ya inequalities involving both $\mathcal{P}_s$ and $\widetilde{\mathcal{P}}_s$. By applying the concentration–compactness principle together with a detailed analysis of the strict-gap regions for the Poincaré–Sobolev levels, we prove the existence of solutions to the Brezis–Nirenberg problem on $\mathbb H^n$ for both operators. Finally, combining the Hardy lower bounds with criteria for attainability, we obtain a complete characterization of the Poincaré–Sobolev levels $H_{n,s}$ and $\widetilde H_{n,s}$.


7. 具有常Tanaka–Webster标量曲率的非齐次完备周期接触形式

Non-homothetic complete periodic contact forms with constant Tanaka–Webster scalar curvature

作者: Jeffrey S. Case, Yuya Takeuchi

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

发布日期: 2026-02-03

arXiv ID: 2602.03241v1

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

摘要: We study the existence problem for complete contact forms with constant Tanaka–Webster scalar curvature on non-compact strictly pseudoconvex CR manifolds. We prove that, under mild assumptions, the universal cover of a compact strictly pseudoconvex CR manifold admits infinitely many non-homothetic such contact forms whenever its fundamental group has infinite profinite completion. As applications, we treat complements of real or complex spheres in the standard CR sphere, as well as circle bundles over compact Kähler manifolds and the boundary of a Reinhardt domain.


8. 双变量Hardy-Sobolev不等式及其精确稳定性

Bivariate Hardy-Sobolev Inequality and Its Sharp Stability

作者: Yingfang Zhang, Xuexiu Zhong, Wenming Zou

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03191v1

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

摘要: This paper establishes a bivariate Hardy-Sobolev inequality. Let $Ω\subset \mathbb{R}^N$ ($N \geq 3$) be an open domain, $s \in (0,2)$, $α> 1$, $β> 1$ with $α+ β= 2^*(s)$, and $κ\in \mathbb{R}$. For any functions $u, v \in D_0^{1,2}(Ω)$, we prove the inequality: \begin{multline*} \int_Ω |\nabla u|^2 , \mathrm{d}x + \int_Ω |\nabla v|^2 , \mathrm{d}x \ge S_{α,β,λ,μ}(Ω) \left( \int_Ω \Big( λ\frac{|u|^{2^(s)}}{|x|^s} + μ\frac{|v|^{2^(s)}}{|x|^s} + 2^(s) κ\frac{|u|^α|v|^β}{|x|^s} \Big), \mathrm{d}x \right)^{\frac{2}{2^(s)}}. \end{multline*} We derive the best constant $S_{α,β,λ,μ}(Ω)$ and characterize the set of minimizers. Moreover, for $Ω= \mathbb{R}^N$ and $κ> 0$, we obtain sharp stability results for nonnegative functions.


9. 关于Fokker-Planck方程相对于Wasserstein距离的定量敏感性分析

Quantitative sensitivity analysis for Fokker-Planck equation with respect to the Wasserstein distance

作者: Martin Morange

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03174v1

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

摘要: We analyze the sensitivity of solutions to the Fokker-Planck equation with respect to some unknown parameter. Our main result is to provide quantitative upper bounds for the $p$-Wasserstein distance $\mathcal{W}_p$ between two solutions with different parameters, for every $p \geq 2$. We are able to give two proofs of this result, the first relying on synchronous coupling between two solutions of an SDE, and another one that relies on the differentiation of Kantorovitch dual formulation of optimal transport. We also provide more specific bounds in the case of the overdamped Langevin process, for which we are able to compare convergence to the invariant measure and sensitivity to the parameter.


10. 具有有界指标的相变:与德乔治猜想的类比

Phase transitions with bounded index: Parallels to De Giorgi’s conjecture

作者: Enric Florit-Simon

类别: math.AP, math.DG

发布日期: 2026-02-03

arXiv ID: 2602.03136v1

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

摘要: A well-known conjecture of De Giorgi – motivated by analogy with the Bernstein problem for minimal surfaces – asserts the rigidity of monotone solutions to the Allen–Cahn equation in $\mathbb{R}^{d+1}$, with $d\leq 7$. We establish close parallels to De Giorgi’s conjecture for general solutions of bounded Morse index, far stronger than the minimal surface analogy would suggest: Namely, any finite index solution to the Allen–Cahn equation with bounded energy density in $\mathbb{R}^4$ is one-dimensional, and – conditionally on the classification of stable solutions – the same holds for all $4\leq n \leq 7$. As a geometric application, phase transitions with bounded energy and index in closed four-manifolds have smooth transition layers which behave like minimal hypersurfaces. Consequently, phase transitions exhibit a remarkably rigid behaviour in higher dimensions. This is in stark contrast with the 3D case, in which a wealth of nontrivial entire solutions with finite index (and energy density) is conversely known to exist, by work pioneered by Del Pino–Kowalczyk–Wei. The authors conjectured that any such solution must have parallel ends which are either planar or catenoidal, suggesting it as a parallel to De Giorgi’s conjecture in this framework. We confirm this picture under the bounded energy density assumption.


11. 关于具有双重相增长的某些椭圆系统之极弱解

On very weak solutions of certain elliptic systems with double phase growth

作者: Yoshiki Kaiho

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.03044v1

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

摘要: In this paper, we prove a higher integrability result for very weak solutions of higher-order elliptic systems involving a double phase operator as the principal part. As a model case, we consider \begin{equation} \int_Ω \left( |D^m u|^{p-2}D^m u + a(x)|D^m u|^{q-2}D^m u \right) \cdot D^m \varphi = 0 \quad \text{for any } \varphi \in C_c^{\infty}(Ω), \end{equation} where $n,m \in \mathbb{N},\ n\ge 2,\,1 < p \le q < \infty,\,Ω\subset \mathbb{R}^n$ is an open set and $a:Ω\rightarrow [0,\infty)$ is a measurable function. The proof is based on a construction of an appropriate test function by the Lipschitz truncation technique, a deduction of a reverse Hölder inequality and an application of Gehring’s lemma. Our contributions include estimates for weighted mean value polynomials and sharp Sobolev–Poincaré-type inequalities for the double phase operator. Our result can be viewed as a generalization with respect to the derivative order, the coefficient function and the growth conditions of the recent paper by Baasandorj, Byun and Kim (Trans. Amer. Math. Soc. 376:8733-8768,2023).


12. 麦克斯韦–诺伊曼–庞加莱算子的对称化、$\mathbf{H}(\mathrm{curl},D)$迹空间中的谱分解及表面等离激元共振模的边界局域化

Symmetrization of the Maxwell–Neumann–Poincar’e operator, spectral decomposition in $\mathbf{H}(\mathrm{curl},D)$ traces, and boundary localisation of SPRs

作者: Bochao Chen, Yixian Gao, Hongyu Liu

类别: math.AP

发布日期: 2026-02-03

arXiv ID: 2602.02998v1

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

摘要: The Neumann–Poincaré (NP) operator, a fundamental operator in potential theory, has attracted renewed attention for its central role in the analysis of surface plasmon resonances (SPRs). SPRs, characterized by non-radiative electromagnetic waves at material interfaces with opposing permittivities, underpin advanced technologies such as bio-sensing and cloaking devices. While spectral properties of the scalar NP operator and SPR dynamics for scalar waves are well-established, their vectorial counterparts in Maxwell’s framework remain poorly understood. This work bridges this gap by introducing a novel symmetrization principle for the matrix-valued Maxwell Neumann–Poincaré (MNP) operator, enabling a spectral decomposition of traces in the $\mathbf{H}(\mathrm{curl},D)$ space–a foundational advance for electromagnetic theory. Building on this framework, we rigorously characterize the quantum-ergodic localization of weak surface plasmon resonances at material boundaries in the full Maxwell system, thereby settling a long-standing question concerning their quantitative description.


13. ℝ⁴ 中的整体极小曲面必为代数曲面

Entire area-minimizing surfaces in R^4 are algebraic

作者: Nick Edelen, Luis Atzin Franco Reyna, Paul Minter

类别: math.DG, math.AP

发布日期: 2026-02-03

arXiv ID: 2602.02997v1

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

摘要: We classify entire 2-dimensional area-minimizing or stable surfaces in R^4 with quadratic area growth as algebraic, cut out by a finite union of holomorphic polynomials whose collective degrees are controlled by the density at infinity. As a consequence, we obtain bounds on the singular set size and genus in terms of the density at infinity.


14. 环面上一类具非线性非局部漂移项的椭圆方程解的不唯一性

Lack of uniqueness for an elliptic equation with nonlinear and nonlocal drift posed on a torus

作者: Adrian Muntean, Giulia Rui

类别: math.AP

发布日期: 2026-02-02

arXiv ID: 2602.02818v1

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

摘要: We study a nonlinear and nonlocal elliptic equation posed on the flat torus. While constant solutions always exist, we show that uniqueness fails in general. Using spectral analysis and the Crandall–Rabinowitz bifurcation theorem, we prove the existence of branches of non-constant periodic solutions bifurcating from constant states. This result is qualitative and non-constructive. Using a conceptually different argument, we construct explicit multiple solutions for a specific one–dimensional formulation of our target problem.


15. 集中型非线性波动方程中常微分方程型爆破解曲面的散射性质与稳定性

Scattering and stability for ODE-type blow-up surfaces for focusing nonlinear wave equations

作者: Istvan Kadar, Warren Li

类别: math.AP

发布日期: 2026-02-02

arXiv ID: 2602.02715v1

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

摘要: We study the focusing power nonlinear wave equation with any power, in Minkowski space of any spacetime dimension. We present a complete understanding of the local stability and scattering theory (both in high regularity spaces) for solutions exhibiting ODE type blow-up on spacelike hypersurfaces, with the blow-up at each point modelled by the explicit solution $φ_{\mathrm{model}} = c_p t^{-α_p}$. Given a sufficiently regular spacelike hypersurface $Σ_f$, together with auxiliary scattering data $ψ$, we construct the unique corresponding solution to the nonlinear wave equation that (locally) forms an ODE type singularity on $Σ_f$ attaining $ψ$ as scattering data. Conversely, we show that such ODE type singularities are (locally) stable to suitably regular perturbations away from the singularity, and that the blow-up surface and scattering data remain regular, in a continuously dependent manner, following such perturbations.


16. 极大次椭圆二次型的先验估计

A Priori Estimates for Maximally Subelliptic Quadratic Forms

作者: Brian Street

类别: math.AP

发布日期: 2026-02-02

arXiv ID: 2602.02441v1

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

摘要: We prove a priori subelliptic estimates, near a non-characteristic boundary point, for the heat operators associated to a wide class of maximally subelliptic quadratic forms. This is the third paper in a series devoted to studying general maximally subelliptic boundary value problems.