生成时间: 2026年08月05日 05:34:34

运行状态: 完成

计划处理论文数量: 24

已完成论文数量: 24

未完成论文数量: 0

重点关注论文数量: 4

了解领域论文数量: 11

不相关论文数量: 9

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

1. Brinkman 项在具有给定通量率或压降的穿孔管道中定常 Navier-Stokes 方程中的出现

Brinkman’s term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

作者: Richard M. Höfer, Amina Mecherbet, Gianmarco Sperone

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03672v1

摘要翻译: 通过带混合边界条件的 Navier-Stokes 方程,对粘性不可压缩流体在含若干直径为 $\eps^3$、相互间距为 $\eps$ 的小颗粒的扭曲管道中的定常运动建模。除颗粒表面的非齐次 Dirichlet 边界条件外,这些边界条件还包含管口入口与出口处的 Bernoulli 压力和切向速度,且沿管道方向要么给定横向通量率,要么给定压降。运用均匀化理论中的能量方法,我们在不加任何数据大小限制的条件下研究了当 $\eps \to 0$ 时该系统解的渐近行为,并证明有效方程中出现一个附加的 Brinkman 项。本工作的关键特征在于所需的一致估计:在给定通量问题的情形下,该估计是通过基于定常 Euler 方程解的 Bernoulli 定律所构造的反证法实现的。

Comment: 无

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

1. 研究对象和背景

本文研究三维穿孔管道中稳态不可压缩粘性流体的运动,其数学模型为具有混合边界条件的Navier-Stokes方程。设Ω⊂R³为满足定义2.1的admissible domain(即由两个圆柱形端口Ω₁=Θ₁×(0,ℓ₁)、Ω₂=Θ₂×(0,ℓ₂)与中间连接区域Ω₀构成的有界单连通区域),其边界分解为入口Γ_I = Θ₁×{0}、出口Γ_O = Θ₂×{ℓ₂}及侧壁Γ_W。对小参数ε∈(0,ε∗),在Ω内嵌入N(ε)个固态颗粒Kⁿ_ε = ξⁿ_ε + ε³Kₙ,其中Kₙ⊂B(0,δ₀)为一致John域,且满足互不相交的分离条件:dist(Kⁿ_ε,Kᵐ_ε)≥δ₁ε,dist(Kⁿ_ε,∂Ω)≥δ₂ε(δ₂>δ₁>0)。由此定义穿孔流体域Ω_ε = Ω\∪_{n=1}^{N(ε)} Kⁿ_ε。流体速度场u_ε:Ω_ε→R³与Bernoulli压力Φ_ε = p_ε + ½|u_ε|²满足如下两类边值问题:

(i)预定通量问题(系统(2.8)):

$$ \begin{cases} -\Delta u_ε + (\nabla \times u_ε) \times u_ε + \nabla \Phi_ε = f & \text{in } \Omega_ε, \\ \nabla \cdot u_ε = 0 & \text{in } \Omega_ε, \\ u_ε = 0 & \text{on } \Gamma_W, \\ u_ε = \muⁿ_ε & \text{on } \partial Kⁿ_ε, \; n=1,\dots,N(ε), \\ u_ε \times \nu = 0, \; -\frac{\partial}{\partial \nu}(u_ε \cdot \nu) + \Phi_ε = 0 & \text{on } \Gamma_I, \\ u_ε \times \nu = 0, \; -\frac{\partial}{\partial \nu}(u_ε \cdot \nu) + \Phi_ε = p⁺_ε & \text{on } \Gamma_O, \\ \int_{\Gamma_O} u_ε \cdot \nu = F & \end{cases} $$

其中f∈L²(Ω;R³)为外力,F∈R为预定横向通量,μⁿ_ε∈R³为第n个颗粒的恒定速度,p⁺_ε∈R为待定出口Bernoulli压力。

(ii)预定压降问题(系统(2.9)):同上动量与连续性方程,但边界条件中Γ_I与Γ_O上的Bernoulli压力被指定为常数p⁻,p⁺∈R,即

$$ -\frac{\partial}{\partial \nu}(u_ε \cdot \nu) + \Phi_ε = p^∓ \quad \text{on } \Gamma_{I/O}. $$

此时通量F_ε = ∫{Γ_O} u_ε·ν成为解依赖量。该问题设定源于实际工程场景:在石油输送、血流动力学及微流控芯片中,常需在固定压差下求解流量,或在固定流量下反演压降,而非经典Dirichlet问题中给定速度边界。传统homogenization理论对Navier-Stokes方程的处理通常依赖于小数据假设(如Leray解的存在性要求‖f‖{H^{-1}}足够小),或仅适用于Stokes线性情形;而本文所考虑的混合边界条件(含Bernoulli压力与切向无滑条件)导致标准能量估计失效——因测试函数u_ε本身不满足u_ε×ν=0在Γ_{I/O}上,故无法直接获得ε-一致先验估计。这一障碍使得既往工作(如Allaire[2]、Sanchez-Palencia[57])无法直接推广。本文核心动机在于填补以下空白:在临界穿孔尺度α=3(即颗粒直径∝ε³、间距∝ε)下,对非小数据、非周期、非零颗粒速度的稳态Navier-Stokes方程建立rigorous homogenization理论,并证明有效方程中出现Brinkman摩擦项Gu−J,其中G为Stokes阻力张量密度极限,J为颗粒运动诱导的源项。此结果将[59]中亚临界情形(α>3)的结论推进至物理上更相关的临界情形,且首次在混合边界条件下实现无小数据限制的均匀估计。

2. 主要定理或主要结果

本文主要结果由定理2.4与2.5给出,在无小数据假设下确立了临界穿孔尺度下的Brinkman homogenization极限。设Ω⊂R³为admissible domain,f∈L²(Ω;R³),(μⁿ_ε){n=1}^{N(ε)}⊂R³满足动能均匀有界条件(2.11):S∗ = sup{ε∈(0,ε∗)} (1/N(ε))∑_{n=1}^{N(ε)} |μⁿ_ε|² < ∞;并假设Stokes容量密度收敛性(2.13):存在G∈L²(Ω;R³×³)与J∈L²(Ω;R³),使得对任意M∈W¹,∞(Ω;R³×³)与φ∈W¹,∞(Ω;R³),有

$$ \lim_{ε→0⁺} ε³ \sum_{n=1}^{N(ε)} G_n M(ξⁿ_ε) = \int_Ω G M, \quad \lim_{ε→0⁺} ε³ \sum_{n=1}^{N(ε)} G_n μⁿ_ε · φ(ξⁿ_ε) = \int_Ω J · φ, $$

其中Gₙ∈R³×³为颗粒Kₙ的Stokes阻力矩阵(见附录A.1)。在此假设下:

• 定理2.4(预定通量问题):设(u_ε,Φ_ε)∈W¹,²(Ω_ε;R³)×L²(Ω_ε;R)为系统(2.8)的弱解,将其延拓至全Ω为(ũ_ε,˜Φ_ε),其中ũ_ε=u_ε在Ω_ε上、ũ_ε=μⁿ_ε在Kⁿ_ε上,˜Φ_ε=Φ_ε在Ω_ε上、˜Φ_ε=0在K_ε上。若(ũ_ε,˜Φ_ε)在W¹,²(Ω;R³)×L²(Ω;R)中有弱聚点(u,Φ),则(u,Φ)是Ω上预定通量Brinkman方程的弱解:

$$ \begin{cases} -\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f, & \nabla \cdot u = 0 \text{ in } \Omega, \\ u = 0 & \text{on } \Gamma_W, \\ u × \nu = 0, \; -\frac{\partial}{\partial \nu}(u · \nu) + \Phi = 0 & \text{on } \Gamma_I, \\ v × \nu = 0, \; -\frac{\partial}{\partial \nu}(u · \nu) + \Phi = p^+ & \text{on } \Gamma_O, \\ \int_{\Gamma_O} u · \nu = F & \end{cases} $$

其中p⁺∈R为某常数。

• 定理2.5(预定压降问题):设(u_ε,Φ_ε)为系统(2.9)的弱解,作相同延拓。若(ũ_ε,˜Φ_ε)有弱聚点(u,Φ),则(u,Φ)满足预定压降Brinkman方程:

$$ \begin{cases} -\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f, & \nabla \cdot u = 0 \text{ in } \Omega, \\ v = 0 & \text{on } \Gamma_W, \\ v × \nu = 0, \; -\frac{\partial}{\partial \nu}(u · \nu) + \Phi = p^- & \text{on } \Gamma_I, \\ v × \nu = 0, \; -\frac{\partial}{\partial \nu}(u · \nu) + \Phi = p^+ & \text{on } \Gamma_O. \end{cases} $$

关键技术细节包括:(1)均匀估计成立的ε范围是I∗=(0,ε∗),且不依赖于f、F或p±的大小;(2)Brinkman项Gu−J的系数G,J由颗粒阻力矩阵Gₙ与速度μⁿ_ε经ε³尺度加权求和后取极限得到,体现集体阻力效应;(3)预定通量问题需额外假设颗粒远离Γ_I与Γ_O(条件(4.3)),以确保通量载体构造的正则性;(4)所有结论均基于弱解框架,解空间为W¹,²(Ω;R³)×L²(Ω;R),压力Φ为Bernoulli压力而非原始压力p。

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

本文的核心技术障碍在于:在临界穿孔尺度α=3且具混合边界条件下,无法通过标准能量测试(即用u_ε自身测试动量方程)获得ε-一致先验估计,因为u_ε在Γ_I/Γ_O上不满足u_ε·ν=const,导致边界项∫_{Γ_I/Γ_O} (½|u_ε|²+p_ε)(u_ε·ν)无法控制。作者的真正创新点在于将Bernoulli定律与矛盾论证相结合,构建了一种新型的ε-一致估计策略。具体而言:

(1)Bernoulli矛盾论证框架:针对预定通量问题(定理4.3),作者假设‖u_ε‖{W¹,²(Ω_ε)}无界,从而存在子列使‖u_ε‖→∞。通过缩放v_ε = u_ε/‖u_ε‖,利用Divergence Theorem与边界条件推导出∫{Γ_I} v_ε·ν = −F/‖u_ε‖→0,结合Bernoulli定律在Euler方程解上的极值性质(即对无粘不可压缩流,Bernoulli常数沿流线守恒),导出矛盾:若v_ε弱收敛到某非零极限v,则v必须满足Euler方程且Bernoulli常数在Γ_I/Γ_O上相等,但这与通量约束F≠0冲突。此论证巧妙规避了对u_ε的直接能量估计,转而利用流体物理定律的内在约束。

(2)Stokes校正器与限制算子:区别于Allaire的经典振荡测试函数法,本文构造了基于Stokes校正器的限制算子R_ε(定理3.6)。对任意φ∈C_c^∞(Ω;R³)满足φ×ν=0在Γ_I/Γ_O上,R_εφ被定义为φ减去在每个Kⁿ_ε邻域内由Stokes阻力矩阵Gₙ生成的局部校正场,使得R_εφ在∂Kⁿ_ε上精确为零。该算子满足:‖R_εφ‖{W¹,²(Ω_ε)} ≤ C‖φ‖{W¹,²(Ω)}且R_εφ→φ强收敛于W¹,²(Ω)。其关键作用在于:(i)将Ω上的测试函数φ“限制”为Ω_ε上的容许测试函数,从而在步骤(c)极限过程中避免振荡;(ii)校正场的构造显式包含Gₙ,使Brinkman项Gu自然涌现于极限方程;(iii)当φ已在∂Kⁿ_ε上为零时R_εφ=φ,保证算子的真实性。

(3)容量势与通量载体:为处理非齐次Dirichlet条件(u_ε=μⁿ_ε),作者利用相对容量势ϕ_ε(满足Δϕ_ε=0在Ω_ε中,ϕ_ε=1在Kⁿ_ε上,ϕ_ε=0在∂Ω上)构造通量载体w_ε∈W²,²(Ω_ε;R³),使得∇·w_ε=0且w_ε=μⁿ_ε在∂Kⁿ_ε上(引理3.4)。关键估计‖ϕ_ε‖{L^∞}+‖∇ϕ_ε‖{L²}≤C∗(引理3.1)依赖于穿孔几何的均匀John性质与分离条件(2.3),确保容量势在ε→0时保持有界,这是临界尺度下Brinkman项非退化的基础。

(4)压力项收敛处理:与多数homogenization工作不同,本文完整证明了Bernoulli压力Φ_ε在L²(Ω)中的弱收敛性(而非仅分布意义),这得益于限制算子R_ε对测试函数的精确控制及Stokes校正器对压力梯度的补偿作用。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开:

(1)尺度临界性与Brinkman项来源:对比Allaire[2]与Sanchez-Palencia[57]的经典结果,其Brinkman项仅出现在α=3的临界情形,且假设颗粒速度为零(μⁿ_ε=0),此时J≡0,Brinkman项仅为Gu。本文将此推广至μⁿ_ε≠0情形,证明J项源自颗粒运动通量极限ε³∑Gₙμⁿ_ε·φ(ξⁿ_ε),这与[12,18,30]中随机颗粒运动的Brinkman模型一致,但本文首次在确定性混合边界条件下实现。相较[59](亚临界α>3),其Brinkman项退化为零(因容量密度→0),故无需Stokes校正器;而本文α=3时容量密度有界(引理3.1),迫使采用Stokes校正器而非标量容量校正器,这是方法论的根本差异。

(2)边界条件与估计技术:Leray通量问题[43]虽也处理非齐次Dirichlet条件,但要求Ω为环形域且使用小数据假设;本文在一般admissible pipe上取消小数据限制,代价是引入Bernoulli矛盾论证——此论证在[54](二维)与[59](亚临界)中已有雏形,但本文首次将其与Stokes校正器耦合,以应对临界尺度下的强非线性。对比[28](随机穿孔),其限制算子未处理非零颗粒速度,且未收敛压力项;本文算子R_ε明确设计为保持φ×ν=0在Γ_I/Γ_O上,确保混合边界条件在极限中精确保留。

(3)适用范围与局限:本文允许非周期、多分散颗粒分布(Kₙ形状各异),优于周期homogenization框架[2,3];但要求颗粒严格分离(dist(Kⁿ_ε,Kᵐ_ε)≥δ₁ε),无法处理接触或聚集情形(如[27,34])。预定通量问题需颗粒远离Γ_I/Γ_O(条件(4.3)),这是为构造光滑通量载体w_ε所需,而预定压降问题无此要求,体现两类问题本质差异。此外,解的存在性依赖于Leray-Schauder原理,未给出唯一性或正则性结果,亦未处理时间演化情形(如[23,33])。

(4)物理建模精度:相较于Darcy律(α<3)或纯Navier-Stokes(α>3),Brinkman方程(α=3)更适配微米级颗粒悬浮液(如生物流体、纳米流体),本文通过混合边界条件(Bernoulli压力+切向无滑)更真实刻画入口/出口处的压力-速度耦合,而经典工作多采用纯Dirichlet或Neumann条件。然而,本文未考虑颗粒动力学(即μⁿ_ε为给定常数而非由流体-颗粒耦合方程决定),此完全耦合问题被作者明确认为‘extremely challenging’(引言末段),从提供内容中无法确认其后续进展。


2. 具有多个涡度空洞的均匀旋转欧拉构型

Uniformly Rotating Euler Configurations with Multiple Vorticity Holes

作者: Vittorio Baroncini, Juan Carlos Cantero, Claudia García, Zineb Hassainia, Joan Mateu

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03654v1

摘要翻译: 我们构造了二维不可压欧拉方程的一类新的均匀旋转涡块解族,该解由一个单连通的外部涡块与多个内部界面组成,这些内部界面在几何上可解释为空洞。更精确地说,每个解由一个单一的外部涡块包围 $\mathbf m\geq2$ 个全同且高度集中的内部组分构成,这些内部组分位于正 $\mathbf m$-边形的顶点上;整个构型以刚性方式顺时针旋转。

Comment: 38 pages, 3 figures

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

1. 研究对象和背景

本文研究二维不可压无粘欧拉方程的涡旋块(vortex-patch)解,其核心对象是如下vorticity形式的初值问题:

$$ \begin{cases} \partial_t \omega + u \cdot \nabla \omega = 0, & (t,x) \in [0,\infty) \times \mathbb{R}^2, \\ u = K * \omega, \quad K(x) = \frac{1}{2\pi} \frac{x^\perp}{|x|^2}, & x^\perp = (-x_2,x_1), \\ \omega(0,x) = \omega_{0,\varepsilon}^\sigma(x), \end{cases} \tag{1.4} $$

其中初始涡量 $\omega_{0,\varepsilon}^\sigma$ 是一个分段常值、紧支、符号变化的构型,由单连通外域 $D_1^\varepsilon$ 和 $m \geq 2$ 个内域 $D_{2,j}^{\varepsilon,\sigma}$ 构成:

$$ \omega_{0,\varepsilon}^\sigma(x) = \mathbf{1}_{D_1^\varepsilon}(x) - \frac{1}{\pi \varepsilon^2} \sum_{j=1}^m \mathbf{1}_{D_{2,j}^{\varepsilon,\sigma}}(x). \tag{1.6} $$

内域通过缩放与旋转定义:

$$ D_{2,j}^{\varepsilon,\sigma} = e^{i \pi \sigma / m} e^{2\pi i (j-1)/m} \left( \varepsilon D_2^{\varepsilon,\sigma} + \varepsilon^\alpha e_1 \right), \quad e_1 = (1,0), \tag{1.7} $$

即 $m$ 个全等区域以正 $m$-边形顶点为几何中心,每个中心距原点距离为 $\varepsilon^\alpha$,自身直径为 $O(\varepsilon)$。参数 $\alpha \in \left( \frac{1}{m+2}, \frac{1}{2} \right)$(对 $m=5$ 可延拓至 $\frac{2}{3}$)控制内域收缩速率与相对分离尺度之间的精细平衡。

该构型的物理意义在于:它代表一种“多空洞涡旋”——并非经典意义下零涡量的几何孔洞(hole),而是负涡量集中区,故称“涡量空洞”。当 $\varepsilon \to 0^+$,外边界 $\partial D_1^\varepsilon$ 在 $C^{1,\nu}$ 拓扑下收敛于单位圆 $\partial \mathbb{D}$,而所有内域收缩并坍缩至原点,整体涡量测度弱收敛为

$$ \omega_{0,\varepsilon}^\sigma \xrightarrow{\text{meas}} \omega_{0,0} := \mathbf{1}_{\mathbb{D}} - m \delta_0, \tag{1.8} $$

即一个单位圆盘上的Rankine涡(恒定正涡量)叠加一个位于原点、总环量为 $-m$ 的点涡。此极限态是径向对称的,因此作为涡量分布是静止的,但其角速度在构造中发散:$\Omega_\varepsilon \sim \frac{1-m}{4\pi \varepsilon^{2\alpha}} \to -\infty$(因 $m \geq 2$)。这揭示了非平凡旋转分支在奇异极限下的退化现象:有限 $\varepsilon$ 下,内域构成一个刚性旋转的正 $m$-边形;极限时该结构完全消失,仅剩旋转不变态。

已有理论对此类构型的研究存在明确空白。经典V-state理论(如Burbea、Hmidi–Mateu–Verdera、Castro–Córdoba–Gómez-Serrano)聚焦于单连通或双连通(annular)构型,其内边界(若存在)为单一闭曲线,且不坍缩;而de la Hoz–Hmidi–Mateu–Verdera对双连通V-states的构造中,内环与外环保持固定分离距离,不趋于重合。点涡反演(desingularization)工作(如Hmidi–Mateu、Hassainia–Hmidi、García、Hassainia–Wheeler)均处理的是多个独立外域(如共转/反转涡对、von Kármán街、Thomson多边形),各涡块彼此分离、不被同一外域包裹。本文首次实现“嵌套式多空洞”的严格构造:所有 $m$ 个内涡块被同一个外涡块包含,并在 $\varepsilon \to 0$ 时同步坍缩至同一中心点。这一设定在数学上引发三重困难:(i) 内–内相互作用项在尺度 $\varepsilon^{1-2\alpha}$ 处主导并破坏线性化算子的Fredholm性质;(ii) 内–外与外–外耦合项需在不同尺度($\varepsilon^\alpha$ vs $1$)下精确匹配;(iii) 对称性破缺(phase $\sigma = 0$ 或 $1$)导致边界方程非平凡耦合,要求精细的Fourier模态分析。因此,本文填补的是“多重嵌套坍缩型涡旋构型”这一长期悬而未决的分析空白。

2. 主要定理或主要结果

定理1.1(主结果):设整数 $m \geq 2$,Hölder指数 $\nu \in (0,1)$,且参数 $\alpha$ 满足 $\frac{1}{m+2} < \alpha < \frac{1}{2}$(当 $m=5$ 时可取 $\alpha < \frac{2}{3}$)。则对相位 $\sigma \in \{0,1\}$,存在 $\varepsilon_0 > 0$,使得对任意 $\varepsilon \in (0,\varepsilon_0)$,存在实数 $\Lambda_\varepsilon$ 及边界属于 $C^{1,\nu}$ 类的有界单连通区域 $D_1^\varepsilon$、$D_2^{\varepsilon,\sigma}$,满足: (1) $D_1^\varepsilon$ 具有 $m$ 重旋转对称性及关于实轴的反射对称性; (2) 由 (1.7) 定义的 $D_{2,j}^{\varepsilon,\sigma}$ 两两不交且严格包含于 $D_1^\varepsilon$; (3) 当 $\varepsilon \to 0^+$ 时,生成域 $D_1^\varepsilon$ 与 $D_2^{\varepsilon,\sigma}$ 在 $C^{1,\nu}$ 拓扑下收敛于单位圆盘 $\mathbb{D}$,且 $\Lambda_\varepsilon \to 0$; (4) 初始涡量 $\omega_{0,\varepsilon}^\sigma$ 生成 (1.4) 的弱解 $\omega_\varepsilon(t,x) = \omega_{0,\varepsilon}^\sigma(e^{-i\Omega_\varepsilon t} x)$,其角速度为

$$ \Omega_\varepsilon = \frac{1-m}{4\pi \varepsilon^{2\alpha}} + \frac{\Lambda_\varepsilon}{\varepsilon^\alpha}. $$

关键技术细节:(i) 角速度主项 $\Omega_0 = \frac{1-m}{4\pi}$ 精确等于 $m$ 个环量为 $1$ 的点涡构成的正 $m$-边形的Kirchhoff–Routh角速度(见 (1.5)–(1.6) 推导);(ii) 修正项 $\Lambda_\varepsilon$ 属于 $C^1$ 类且 $|\Lambda_\varepsilon| = O(\varepsilon^{\beta})$,其中 $\beta > 0$ 依赖于 $\alpha$ 和 $\nu$,由隐函数定理保证;(iii) 边界正则性 $C^{1,\nu}$ 来自对称适配Hölder空间 $C^{1,\nu}_{\text{sym}}(\mathbb{T})$ 中的解存在性;(iv) 对 $m=5$,因五边形第二阶Fourier系数的精确抵消(Remark 1.2),允许更大的 $\alpha$ 上界,此时内–内相互作用余项阶为 $\varepsilon^{2-3\alpha}$ 而非 $\varepsilon^{1-2\alpha}$;(v) 解的局部唯一性由线性化算子 $D_{(\Lambda,f)} F^\sigma(0,0,0)$ 的同构性保证,其显式形式为

$$ D_{(\Lambda,f)} F^\sigma(0,0,0)[\lambda,h] = \left( \Omega_0 h_1', \; -\lambda \sin(\cdot) + \frac{1}{2\pi} (h_2' + H[h_2]) \right), $$

其中 $H$ 为周期Hilbert变换,该算子在Fourier基下对角化,且第一分量控制径向扰动,第二分量中 $\lambda$ 提供 $\sin\theta$ 模,$h_2$ 控制所有更高阶模。

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

证明的核心障碍在于:当 $\varepsilon \to 0$ 时,内域间相互作用产生的奇异性远超经典单涡块或双涡块情形。具体而言,$m$ 个内域中心间距为 $O(\varepsilon^\alpha)$,而其自身尺度为 $O(\varepsilon)$,导致内–内速度场贡献在边界方程中以 $O(\varepsilon^{1-2\alpha})$ 阶主导(参见Section 3推导),该阶次在 $\alpha < 1/2$ 下趋于无穷,彻底破坏标准线性化算子的紧性与Fredholm性质。作者的真正创新点在于构建了一个三重尺度分离的奇异重标定框架:

第一,对称性约化与参量化:利用 $m$ 重旋转对称性,将原始 $m+1$ 个边界方程(外边界 + $m$ 个内边界)约化为仅含两个未知函数的系统:外边界 $f_1(\theta)$ 与参考内边界 $f_2(\theta)$(均表示为极坐标径向图)。相位 $\sigma$ 通过旋转因子 $e^{i\pi\sigma/m}$ 编码进核函数,使非线性算子 $F^\sigma(\varepsilon,\Lambda,f)$ 显式依赖于 $\sigma$,但其零阶极限 $F^\sigma(0,0,0)=0$ 与 $\sigma$ 无关。

第二,奇异重标定与主导项消去:将角速度设为 $\Omega_{\varepsilon,\Lambda} = \frac{\Omega_0}{\varepsilon^{2\alpha}} + \frac{\Lambda}{\varepsilon^\alpha}$,其中 $\Omega_0 = \frac{1-m}{4\pi}$ 的选取精确抵消了内边界方程中由 $m$ 个点涡正多边形产生的 $O(\varepsilon^{-2\alpha})$ 奇异项。此步是构造可行性的基石——若 $\Omega_0$ 偏离该值,余项将无法被后续Hölder估计控制。

第三,连续参数隐函数定理:由于重标定引入 $\varepsilon^{\alpha}$、$\varepsilon^{2\alpha}$ 等非整数幂,参数空间取为度量空间 $[0,\varepsilon_0)$,应用Dontchev–Rockafellar的参数依赖隐函数定理(Theorem 5F.4)。其关键前提是验证 $F^\sigma$ 在 $\varepsilon=0$ 处的连续延拓及Fréchet可微性,这依赖于附录A中建立的对称适配奇异积分估计:包括对Biot–Savart核在 $m$-对称域上卷积的精细渐近展开、Hilbert变换在 $C^{1,\nu}_{\text{sym}}$ 空间中的有界性,以及内–外、外–外、内–内三类相互作用项在不同Fourier模上的衰减率分析(例如,内–内项在第 $k$ 次谐波上衰减为 $O(\varepsilon^{|k|-2\alpha})$)。

第四,线性化算子的显式谱分析:在极限配置 $(\varepsilon,\Lambda,f)=(0,0,0)$ 处,线性化 $D_{(\Lambda,f)}F^\sigma$ 的第二分量含Hilbert变换 $H$,其在 $L^2$ 上谱为 $-i\,\text{sgn}(k)$,结合 $\frac{1}{2\pi}(h_2'+H[h_2])$ 项,确保对所有 $k \neq \pm 1$ 的高阶模可逆;而 $\lambda \sin\theta$ 项专门补偿 $k=1$ 模的零空间,使整个算子成为Banach空间间的同构。此结构依赖于 $\sigma=0$ 或 $1$ 下反射对称性对Fourier系数的约束(仅保留偶次谐波),从而排除了其他可能的零模。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度:

1. 几何构型与坍缩机制:此前所有rigorous desingularization构造(Hmidi–Mateu [26], Hassainia–Hmidi [20], García [12], Hassainia–Wheeler [22])均处理分离式多涡块(disjoint patches),即每个涡块位于不同连通分支,彼此距离 $O(1)$。而本文是首个实现嵌套式多坍缩(nested multi-collapse)的构造:所有 $m$ 个内涡块被同一外域包裹,且其中心距原点为 $\varepsilon^\alpha$,自身直径为 $\varepsilon$,故当 $\varepsilon \to 0$ 时它们同步坍缩至同一点。这与双连通V-states(de la Hoz et al. [7])有根本不同:后者内环半径 $b>0$ 固定,不坍缩;其极限 $b \to 0$ 导致分支退化(Remark 1.4),但该退化不产生新构型,而本文的 $\varepsilon \to 0$ 是主动设计的构造参数,旨在生成新解族。

2. 假设放宽与技术代价:经典V-state构造(Burbea [2], Hmidi–Mateu–Verdera [27])要求 $m$-fold对称性及小变形,但仅适用于单连通或环状;本文虽同样假设 $m$-对称,却成功处理了更复杂的嵌套拓扑。代价是:必须限制 $\alpha$ 的范围($\frac{1}{m+2}<\alpha<\frac{1}{2}$),此约束源于内–内相互作用项在极限方程中需趋于零(否则破坏线性化),而该条件在现有技术下无法进一步放宽。对比Hassainia–Wheeler [22]对一般点涡构型的构造,其非退化条件是代数的(Jacobian非零),而本文的非退化性是解析的,依赖于 $\alpha$ 的精细选择及 $\sigma$ 的离散性(仅 $\sigma=0,1$ 有效)。

3. 时间尺度与动力学含义:本文解的角速度 $\Omega_\varepsilon \sim \varepsilon^{-2\alpha}$ 发散,对应时间尺度 $|\Omega_\varepsilon|^{-1} \sim \varepsilon^{2\alpha} \to 0$,即内多边形旋转越来越快。这与Hmidi–Mateu [26]中涡对的角速度 $\Omega = \Omega_0 + O(\varepsilon^2)$ 收敛形成鲜明对比。该发散性并非缺陷,而是几何退化的必然体现:极限态 $\mathbf{1}_\mathbb{D} - m\delta_0$ 是径向对称的,其“角速度”无几何意义;发散的 $\Omega_\varepsilon$ 正是有限 $\varepsilon$ 下非平凡旋转结构的定量刻画。

4. 解的类型与物理诠释:本文构造的 $\omega_{0,\varepsilon}^\sigma$ 在内域取值 $1 - \frac{1}{\pi\varepsilon^2} < 0$(因 $\varepsilon$ 小),故为负涡量区,非零涡量空洞。Remark 1.3指出,当 $\varepsilon_h = 1/\sqrt{\pi}$ 时,内域涡量恰为零,此时解成为真正的 $m$-连通V-state。但本文仅保证局部存在性($\varepsilon < \varepsilon_0$),从提供内容中无法确认能否全局延拓至 $\varepsilon_h$ 而不发生界面碰撞或退化——这正是未来工作的关键挑战,也凸显了本文构造的前沿性与局限性。


3. 基于 $\varphi$-散度的矩闭合方法分析:适用于具有二元碰撞的稀薄动力学及其Galerkin离散化

Analysis of Moment Closures Using $\varphi$-Divergences for Rarefied Dynamics with Binary Collisions and Their Galerkin Discretizations

作者: Michael R. A. Abdelmalik, Irene M. Gamba, Torsten Kessler, Sergej Rjasanow

类别: math.NA, math.AP, physics.comp-ph, physics.flu-dyn

发布日期: 2026-08-04

arXiv ID: 2608.03640v1

摘要翻译: 本工作引入一种鲁棒的确定性框架,通过Galerkin方法对Boltzmann方程(含二元碰撞)解在时间、位置与速度变量上的依赖关系进行离散化,从而逼近其解。借助一族基于速度空间中 $\varphi$-散度的参数化Galerkin闭合模型,我们推导出严格控制流体动力学变量的矩方程层级体系。针对此类闭合模型本身不能保证对真实二元碰撞算子耗散 $\varphi$-散度熵这一局限,我们通过为每种闭合模型专门构造一个相容的近似碰撞算子,恢复该耗散性质。所构造的算子内在地保持高保真流动模拟所需的基本物理性质,包括伽利略不变性、质量、动量与能量的精确守恒,以及 $\varphi$-散度熵的严格耗散。此外,我们证明所得闭合矩系统是对称-耗散的,从而其对应的Cauchy问题在局部时间上适定。为将此数学基础转化为高效计算工具,我们采用熵稳定间断Galerkin(DG)有限元方法对位置与时间变量进行离散化。该全隐式、熵稳定的时空方法允许时间步长远超典型CFL限制步长,并支持稳态解的直接计算。该方法的鲁棒性与精度通过氩气超声速喷管流、通道内质量流及平行壁间热传导等数值模拟予以验证与确认,结果与解析基准解、实验测量及随机粒子模拟均吻合良好。

Comment: 35 pages, 8 figures, 2 tables

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

1. 研究对象和背景

本文研究对象是单原子气体在低Knudsen数过渡区的稀薄动力学行为,其数学模型为全空间速度变量上的Boltzmann方程:

$$ \partial_t f + v \cdot \nabla_x f = C(f), \quad f : (0,\infty) \times \Omega \subset \mathbb{R}^3 \times \mathbb{R}^3 \to (0,\infty), $$

其中二元碰撞算子 $C$ 具有标准Boltzmann–Grad形式:

$$ C(f)(v) = \int_{\mathbb{R}^3} \int_{S^2} B(v - v_*, \sigma) \big( f(\tilde{v}_*) f(\tilde{v}) - f(v_*) f(v) \big) \, d\sigma \, dv_*, $$

其中后碰撞速度由守恒律定义:

$$ \tilde{v} = \frac{1}{2}(v + v_*) + \frac{1}{2}|v - v_*|\sigma, \quad \tilde{v}_* = \frac{1}{2}(v + v_*) - \frac{1}{2}|v - v_*|\sigma, \quad \sigma \in S^2. $$

碰撞核取幂律形式 $B(v - v_*, \sigma) = |v - v_*|^\lambda b(\hat{u} \cdot \sigma)$,其中 $\lambda \in (-3,1]$,涵盖硬球($\lambda=1$)、Maxwell分子($\lambda=0$)及库仑极限($\lambda \to -3$)。该模型的核心物理结构包括:(i)Galilean不变性——对任意平移 $T_z f(v) = f(v - z)$ 和正交旋转 $T_O f(v) = f(O^\top v)$,满足 $\int_{\mathbb{R}^3} \psi(v) C(f)(v)\,dv = \int_{\mathbb{R}^3} \psi(v) C((T_z T_O T_{-z})^* f)(v)\,dv$;(ii)五维碰撞不变量空间 $I = \operatorname{span}\{1, v_1, v_2, v_3, |v|^2\}$,满足 $\int_{\mathbb{R}^3} \psi(v) C(f)(v)\,dv = 0$ 当且仅当 $\psi \in I$;(iii)Boltzmann熵耗散:对 $\eta(f) = f \ln f - f$,有 $\int_{\mathbb{R}^3} \ln f \, C(f)\,dv \leq 0$,等号成立当且仅当 $f = M_{\rho,u,\theta}(v) = \rho (2\pi\theta)^{-3/2} \exp\big(-|v-u|^2/(2\theta)\big)$。这些结构共同保证了从动理学到流体力学的渐近一致性(如Euler与Navier–Stokes–Fourier极限),并构成数值方法设计的刚性约束。

矩方法试图通过有限维试函数空间 $\mathcal{M} = \operatorname{span}\{m_1,\dots,m_K\} \supset I$ 构造闭包:对给定 $f$,定义其矩闭包 $f_{\mathrm{KL}}$ 为最小化 $\varphi$-散度 $D_\varphi(f_{\mathrm{KL}} \| M) = \int_{\mathbb{R}^3} M \varphi(f_{\mathrm{KL}}/M)\,dv$ 的解,满足矩约束 $\int m_j f_{\mathrm{KL}}\,dv = \int m_j f\,dv$。经典最大熵闭包取 $\varphi(z) = z \ln z - z$,对应指数型重构 $f_{\mathrm{KL}} = M e^{\lambda \cdot m}$,但已被证明在超二次多项式空间中存在实可实现性(realizability)缺陷:存在邻近平衡态的矩向量无法被任何 $e^{\lambda \cdot m}$ 表示,导致优化问题无解、初值问题不适定(ill-posed),且输运项在平衡点附近可能奇异。本文聚焦于一类由 [26] 提出的正则化 $\varphi$-散度族:

$$ \varphi_N(z) = z^N \left( \frac{z^{1/N} N}{N+1} - 1 \right) + \frac{N}{N+1}, \quad N \in \mathbb{N}\ \text{odd}, $$

其对应重构映射为 $\beta_N(g) = M (1 + g/N)^N$,该族满足 $\lim_{N\to\infty} \varphi_N(z) = z \ln z - z + 1$ 且对任意多项式空间 $\mathcal{M}$,重构方程 $\int m \, \beta_N(\lambda \cdot m)\,dv = \int m f\,dv$ 存在唯一解(Theorem 1)。然而,关键空白在于:尽管 $\beta_N$ 保证了闭包的存在性与唯一性,但 $\int \varphi'_N(f/M) C(f)\,dv$ 的符号性质完全未知——即经典碰撞算子 $C$ 并不耗散 $\varphi_N$-熵,这使得由此导出的矩系统缺乏热力学稳定性保障,无法支撑长时间稳定数值模拟。本文旨在严格填补这一理论空白:构造一族与 $\varphi_N$-熵严格兼容的近似碰撞算子 $C_N$,使其同时满足Galilean不变性、五维不变量守恒、以及 $\int \varphi'_N(f/M) C_N(f)\,dv \leq 0$,从而建立具有对称耗散结构的闭合矩系统。

2. 主要定理或主要结果

本文在以下核心假设下建立了完整的理论框架:(1)速度空间上固定背景Maxwellian $M(v) = \rho_0 (2\pi\theta_0)^{-3/2} \exp(-|v-u_0|^2/(2\theta_0))$;(2)试函数空间 $\mathcal{M} = \operatorname{span}\{m_1,\dots,m_K\}$ 包含碰撞不变量 $I$ 且由至多总次数 $K$ 的三维多项式张成;(3)$N \in \mathbb{N}$ 为奇数;(4)初始分布 $f_0 \in L^1_K(\mathbb{R}^3)$ 满足 $f_0 > 0$ a.e. 且矩向量实可实现。

主要结果如下:

Theorem 2: 对任意 $u \in \mathbb{R}^3$(即 $M$ 的平均速度),近似碰撞算子 $C_N$ 相对于移动背景具有Galilean不变性:设 $T = \tau_{-u} \circ T_O \circ \tau_u$,其中 $\tau_u$ 为平移算子、$T_O$ 为旋转算子,则 $C_N \circ T = T \circ C_N$。

Theorem 3: $C_N$ 的弱形式 $\int_{\mathbb{R}^3} \psi C_N(f)\,dv = 0$ 对所有 $f \in D(C)$ 成立当且仅当 $\psi(v) = a|v|^2 + b \cdot v + c$,即 $C_N$ 的碰撞不变量空间严格等于 $I$,不引入额外不变量。

Theorem 4: $C_N$ 严格耗散 $\varphi_N$-熵:对任意 $f \in D(C)$,有

$$ \int_{\mathbb{R}^3} \varphi'_N\left(\frac{f(v)}{M(v)}\right) C_N(f)(v)\,dv \leq 0, $$

且等号成立当且仅当 $\varphi'_N(f/M) \in I$,即 $f$ 为近似Maxwellian $f(v) = M(v) \big(1 + (a + b \cdot v + c|v|^2)/N\big)^N$(Corollary 1)。

Theorem 5: 所有 $C_N$ 在背景 $M$ 处的线性化相同,且等于经典线性化算子 $L$:

$$ Lg(v) = \int_{\mathbb{R}^3} \int_{S^2} B(v-v_*,\sigma) \big( M(\tilde{v}) g(\tilde{v}_*) + M(\tilde{v}_*) g(\tilde{v}) - M(v) g(v_*) - M(v_*) g(v) \big)\,d\sigma\,dv_*; $$

特别地,$C_1 = L$ 是线性的。

Theorem 6: 序列 $\{C_N\}_{N\in2\mathbb{N}+1}$ 弱收敛于真实碰撞算子 $C$:对任意 $f \in D(C)$ 满足 $f>0$ a.e. 及适当测试函数 $\psi \in C_c(\mathbb{R}^3)$,有

$$ \lim_{N\to\infty} \int_{\mathbb{R}^3} \psi(v) C_N(f)(v)\,dv = \int_{\mathbb{R}^3} \psi(v) C(f)(v)\,dv. $$

Theorem 8(数值离散): 对闭合矩系统采用全隐式时空间断Galerkin(DG)离散,若时间步长 $\Delta t$ 满足 $\Delta t \leq \Delta t_{\max}$(由能量范数与DG跳跃项决定),则离散熵 $\mathcal{E}_h^N(t) = \sum_{K} \int_K \int_{\mathbb{R}^3} M \varphi_N\big( f_h^N / M \big)\,dv\,dx$ 满足

$$ \mathcal{E}_h^N(t^{n+1}) \leq \mathcal{E}_h^N(t^n), \quad \forall n \geq 0, $$

即离散格式严格熵稳定。此外,该格式允许 $\Delta t$ 远超CFL限制,并可直接求解稳态解。

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

本文的核心技术障碍在于:经典矩闭包 $f_{\mathrm{KL}} = \beta_N(\lambda \cdot m)$ 虽保证了实可实现性与唯一性,但其导出的 $\varphi_N$-熵 $\int \varphi_N(f_{\mathrm{KL}}/M)\,dv$ 并非 $C$ 的李雅普诺夫泛函——即 $\int \varphi'_N(f_{\mathrm{KL}}/M) C(f_{\mathrm{KL}})\,dv$ 符号不定,导致闭合系统缺乏内在耗散结构,无法保证适定性与数值稳定性。作者的真正创新点在于提出一种“熵兼容算子构造范式”:不是修改闭包本身,而是为每个 $\varphi_N$-熵定制一个近似碰撞算子 $C_N$,使其在保持物理结构的同时严格耗散该熵。该构造的关键工具是“$N$阶近似乘积” $P_N(z_1,z_2)$,定义为

$$ P_N(z_1,z_2) = \left( z_1^{1/N} + z_2^{1/N} - 1 \right)^N, \quad z^{1/N} = -|z|^{1/N} \text{ for } z < 0. $$

此定义巧妙地将 $\varphi_N$ 的导数与乘积结构耦合:Lemma 2 证明了函数方程 $\varphi'_N(P_N(z_1,z_2)) = \varphi'_N(z_1) + \varphi'_N(z_2)$,这正是熵耗散所需的可加性。基于此,$C_N$ 被定义为(Definition 2)

$$ C_N(f)(v) = \int_{\mathbb{R}^3} \int_{S^2} B(v-v_*,\sigma) M(v) M(v_*) \left[ P_N\left( \frac{f(\tilde{v})}{M(\tilde{v})}, \frac{f(\tilde{v}_*)}{M(\tilde{v}_*)} \right) - P_N\left( \frac{f(v)}{M(v)}, \frac{f(v_*)}{M(v_*)} \right) \right] d\sigma dv_*. $$

该定义的精妙之处在于:它将 $C$ 中的非线性项 $f(\tilde{v}) f(\tilde{v}_*)$ 替换为 $P_N$,而 $P_N$ 的构造确保了 $\varphi'_N$ 的可加性,从而直接导出 Theorem 4 的耗散不等式。第二个关键技术是弱收敛性证明(Theorem 6):作者利用 $P_N(z_1,z_2)$ 在 $N\to\infty$ 时一致收敛于 $z_1 z_2$(因 $\lim_{N\to\infty} (a^{1/N} + b^{1/N} -1)^N = ab$),结合碰撞核 $B$ 的可积性与 $f$ 的正则性假设,通过控制收敛定理(Dominated Convergence Theorem)证得 $\int \psi C_N(f)\,dv \to \int \psi C(f)\,dv$。第三个关键技术是适定性证明框架:作者将闭合矩系统识别为“对称耗散系统”(symmetric-dissipative system)[33],其标准形式为 $A^0(u) u_t + \sum_i A^i(u) u_{x_i} = Q(u)$,其中 $A^0$ 正定、$A^i$ 对称、$Q$ 满足 $u^\top Q(u) \leq 0$。此处,$u = \lambda$ 为Lagrange乘子向量,$A^0_{ij} = \int m_i m_j \beta'_N(\lambda \cdot m) M\,dv$,$A^i_{ij} = \int v_i m_i m_j \beta'_N(\lambda \cdot m) M\,dv$,$Q_i = \int m_i C_N(\beta_N(\lambda \cdot m))\,dv$。Theorem 4 确保了 $\lambda^\top Q(\lambda) \leq 0$,而 $\beta'_N > 0$ 保证了 $A^0$ 正定,从而应用 [33] 的局部适定性定理得到有限时间解的存在唯一性。最后,熵稳定DG离散的核心工具是“熵修正通量”(entropy-stable flux):在单元边界上,作者采用以 $\varphi_N$-熵为凸势的Lax–Friedrichs型数值通量,并严格验证其满足离散熵不等式(Theorem 8),这依赖于 $C_N$ 的耗散性与DG投影算子的正交性。

4. 与之前工作的比较

本文与既有工作存在本质性推进,体现在四个维度:

1. 与经典矩方法(Grad, Levermore)的比较:Grad的13矩方法($\mathcal{M} = \operatorname{span}\{1,v_i,v_i v_j,|v|^2\}$)虽具物理直观性,但其指数型闭包 $f = \exp(\lambda \cdot m)$ 在高阶矩空间中存在实可实现性缺陷([30,35]),导致初值问题在平衡态附近不适定,且无法保证熵耗散。本文通过 $\beta_N(g) = M(1+g/N)^N$ 彻底解决了此问题:Theorem 1 证明其对任意奇数 $N$ 均存在唯一重构解,且 $N=1$ 时退化为Galerkin–Petrov方法,$N\to\infty$ 时收敛于最大熵闭包。代价是引入参数 $N$,需在精度(大 $N$)与计算成本(高次积分)间权衡。

2. 与BGK/Shakhov模型的比较:BGK算子 $C_{\mathrm{BGK}} = \nu (M_f - f)$ 虽满足质量/动量/能量守恒与熵耗散,但其弛豫率 $\nu$ 为标量,无法刻画二元碰撞的真实非线性动力学(如各向异性弛豫、高阶矩演化)。本文的 $C_N$ 严格继承 $C$ 的双线性结构与 $B(v-v_*,\sigma)$ 核,Theorem 5 证明其线性化与 $C$ 完全一致,Theorem 6 证明其弱收敛于 $C$,因此在物理保真度上远超BGK类模型。代价是 $C_N$ 的计算复杂度与 $C$ 相当(仍需三重积分),但本文通过Galerkin离散与高效求积规避了此问题。

3. 与Fourier-Galerkin(FGM)和离散速度法(DVM)的比较:FGM虽能快速计算 $C$,但会破坏Galilean不变性与熵耗散([17,20]);DVM则受限于速度域截断与网格对齐问题([22,23])。本文的矩方法作用于全速度空间 $\mathbb{R}^3$,避免截断误差;而 $C_N$ 的构造(Theorem 2,3,4)显式保证了Galilean不变性、五维不变量守恒与 $\varphi_N$-熵耗散,这是FGM/DVM无法提供的结构性保障。

4. 与现有熵稳定数值方法的比较:Tadmor等人的熵稳定格式([43,44])针对守恒律设计,但未解决动理学矩系统的内在耗散缺失问题;LeFloch–Tzavaras框架([45,46])要求源项满足特定符号条件。本文首次为 $\varphi_N$-矩系统构造了严格耗散的 $C_N$,并将此结构无缝嵌入全隐式时空DG框架(Theorem 8),实现了 $\Delta t$ 远超CFL限制的稳定大步长模拟,这是DSMC(受统计噪声限制)与传统显式动理学求解器(受CFL限制)均无法达到的。从提供内容中无法确认本文方法与最近基于机器学习的矩闭包(如[49])的定量比较,作者亦未显式展开此类比较。


4. 调和映射热流进入扰动球面的稳定爆破

Stable blowup for the harmonic map heat flow into perturbed spheres

作者: Alexander Wittenstein

类别: math.AP, math.FA

发布日期: 2026-08-04

arXiv ID: 2608.03332v1

摘要翻译: 我们证明了在维数 $3$ 到 $6$ 时,对目标流形为标准球面的微小扰动的情形,调和映射热流存在稳定的自相似爆破解。我们从球面上已知的爆破解出发,为此类目标流形构造出自相似爆破解,并证明这些解是渐近非线性稳定的。

Comment: 无

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

1. 研究对象和背景

本文研究 $d$ 维欧氏空间 $\mathbb{R}^d$ 到一类微小几何扰动球面 $S^d_\varepsilon$ 的调和映射热流(harmonic map heat flow),其演化方程为

$$ \partial_t U = \tau(U) \quad \text{on } [0,\infty) \times \mathbb{R}^d, $$

其中 $U : [0,\infty) \times \mathbb{R}^d \to S^d_\varepsilon \subset \mathbb{R}^{d+1}$,且 $\tau(U)$ 是目标流形 $(S^d_\varepsilon, h)$ 上的张力场:

$$ \tau(U)^\alpha = \Delta_{\delta} U^\alpha + \Gamma^\alpha_{\beta\gamma}(U) \, \partial_i U^\beta \, \partial_i U^\gamma, \quad \alpha = 1,\dots,d+1. $$

这里 $\Delta_{\delta}$ 是 $\mathbb{R}^d$ 上的标准拉普拉斯,$\Gamma^\alpha_{\beta\gamma}$ 是 $S^d_\varepsilon$ 的 Christoffel 符号。目标流形 $S^d_\varepsilon$ 定义为具有度量 $h = du^2 + w_\varepsilon(u)^2 d\Omega^2$ 的扭曲乘积 $S^d_\varepsilon = (0,\pi) \times_{w_\varepsilon} S^{d-1}$,其中扭曲函数 $w_\varepsilon(u) = \sin u \, (1 + \varepsilon \, \alpha(u))$,$\alpha \in C^\infty(\mathbb{R})$ 是偶的、$2\pi$-周期、满足 $\alpha(0)=\alpha(\pi)=0$ 的非平凡函数,且 $|\varepsilon| \leq \varepsilon_0 := (2 \|\alpha\|_{L^\infty})^{-1}$。当 $\varepsilon = 0$ 时,$w_0(u) = \sin u$,$S^d_0$ 即标准单位球面 $S^d$。在径向对称(corotational)设定下,映射形如 $U(t,x) = u(t,|x|) \, x/|x|$,其中 $u : [0,\infty) \times [0,\infty) \to \mathbb{R}$ 满足一维径向方程

$$ \left( \partial_t - \partial_r^2 - \frac{d-1}{r} \partial_r \right) u(t,r) + \frac{d-1}{r^2} \, w_\varepsilon(u(t,r)) \, w'_\varepsilon(u(t,r)) = 0. \tag{1.10} $$

该方程在 $d \geq 3$ 时属能量超临界情形(因 Dirichlet 能量 $E[U] = \frac{1}{2} \int_{\mathbb{R}^d} |\nabla U|^2 dx$ 在缩放 $U_\lambda(t,x) = U(t/\lambda^2, x/\lambda)$ 下变为 $E[U_\lambda](t) = \lambda^{d-2} E[U](t/\lambda^2)$,故 $d > 2$ 时能量不守恒)。已知在 $d = 3,4,5,6$ 时,$\varepsilon = 0$ 对应的球面情形存在无穷多个自相似爆破解 $u^T_0(t,r) = \hat f_0(r / \sqrt{T-t})$,其中 $\hat f_0$ 满足常微分方程

$$ f''(\rho) + \left( \frac{d-1}{\rho} - \frac{\rho}{2} \right) f'(\rho) - \frac{d-1}{\rho^2} \sin f(\rho) \cos f(\rho) = 0. $$

Biernat–Donninger–Schörkhuber($d=3$)及后续工作($d=4,5,6$)已证明这些解是渐近非线性稳定的:对足够接近 $u^T_0$ 的初值,解仍于有限时间 $T$ 爆破,且在相似变量 $\rho = r/\sqrt{T-t}$ 下收敛至同一 $\hat f_0$。然而,所有既有稳定性结果均严格依赖于目标流形为精确球面 $S^d$ 的对称性与显式结构;一旦几何发生扰动(即使 $C^\infty$ 小扰动),Christoffel 符号、非线性项及谱结构均改变,原有稳定性论证失效。本文旨在填补这一关键空白:在不假设目标流形具有任何显式对称性或可积结构的前提下,证明当 $S^d$ 被 $C^\infty$ 小扭曲(即 $|\varepsilon| \ll 1$)时,自相似爆破机制仍以 Lipschitz 连续方式延续,并保持其渐近非线性稳定性。这触及了 PDE 奇点理论的核心问题——爆破机制对几何背景的鲁棒性。若该机制脆弱,则其物理意义受限;若鲁棒,则表明此类奇点是能量超临界调和映射热流的内在结构性现象,而非特定对称性的偶然产物。

2. 主要定理或主要结果

定理 1.2(存在性):设 $3 \leq d \leq 6$,则存在 $\varepsilon_* > 0$,使得对任意 $|\varepsilon| \leq \varepsilon_*$ 及任意 $T > 0$,方程 (1.10) 存在自相似爆破解 $u^T_\varepsilon(t,r) = \hat f_\varepsilon(r / \sqrt{T-t})$,满足 $\lim_{t \to T^-} |\partial_r u^T_\varepsilon(t,0)| = \infty$。进一步,$\hat f_\varepsilon$ 可分解为 $\hat f_\varepsilon(\rho) = \hat f_0(\rho) + \rho \, \tilde\phi_\varepsilon(\rho)$,其中 $\hat f_0$ 是 $\varepsilon = 0$ 时的已知自相似剖面,且扰动项 $\tilde\phi_\varepsilon$ 关于 $\varepsilon$ Lipschitz 连续:$\| \tilde\phi_\varepsilon - \tilde\phi_\kappa \|_{W^{2,\infty}([0,\infty))} \lesssim |\varepsilon - \kappa|$。此外,$\hat f_\varepsilon$ 为奇函数,满足 $0 < \hat f_\varepsilon(\rho) < \pi$($\forall \rho > 0$),$\lim_{\rho \to \infty} \hat f_\varepsilon(\rho)$ 存在,且具快速衰减:对任意 $k \in \mathbb{N}$,存在 $C_{\varepsilon,k} > 0$ 使得 $|\hat f^{(k)}_\varepsilon(\rho)| \leq C_{\varepsilon,k} \langle \rho \rangle^{-2-k}$($\forall \rho \geq 0$)。定理 1.3(稳定性):设 $3 \leq d \leq 6$,取 $\varepsilon_*$ 同上。令 $U^T_\varepsilon(t,x) = \hat f_\varepsilon(|x|/\sqrt{T-t}) \, x/|x|$,并考虑初值 $U_0 = U^1_\varepsilon(0,\cdot) + \nu_0$,其中 $\nu_0(x) = x v_0(|x|)$,$v_0 \in \mathcal{S}(\mathbb{R}^d)$。设 $(s,k) \in \mathbb{R} \times \mathbb{N}$ 满足

$$ \frac{d}{2} < s \leq \frac{d}{2} + \frac{1}{2d + 2}, \quad k = d + 3. $$

则存在 $\varepsilon > 0$($\leq \varepsilon_*$),及常数 $\delta > 0$、$M_0 > 1$,使得对任意满足 $\|\nu_0\|_{\dot{H}^s \cap \dot{H}^k(\mathbb{R}^d; \mathbb{R}^d)} < \delta / M_0$ 的 $\nu_0$,存在唯一 $T_\varepsilon \in [1-\delta, 1+\delta]$ 及唯一径向对称解 $U \in C^\infty([0,T_\varepsilon) \times \mathbb{R}^d; \mathbb{R}^d)$ 满足 (1.8),其梯度在 $(T_\varepsilon, 0)$ 处爆破,且可分解为

$$ U(t,x) = U^{T_\varepsilon}_\varepsilon(t,x) + \nu\left(t, \frac{x}{\sqrt{T_\varepsilon - t}}\right), $$

其中余项 $\nu : [0,T_\varepsilon) \times \mathbb{R}^d \to \mathbb{R}^d$ 满足:对任意 $r \in [s,k]$,有 $\|\nu(t,\cdot)\|_{\dot{H}^r(\mathbb{R}^d; \mathbb{R}^d)} \to 0$(当 $t \to T_\varepsilon^-$);且 $U(t, \sqrt{T_\varepsilon - t} \, x) \to \hat f_\varepsilon(|x|) \, x/|x|$ 在 $\mathbb{R}^d$ 的任意紧集上一致成立。定理 1.5(相似坐标下的稳定性):令 $n = d + 2$(故 $5 \leq n \leq 8$),并设 $(s,k)$ 满足

$$ \frac{n}{2} - 1 < s \leq \frac{n}{2} - 1 + \frac{1}{2n - 2}, \quad k = n + 1. $$

则存在 $\omega > 0$ 及 $0 < \varepsilon \leq \varepsilon_*$,使得对任意 $|\varepsilon| \leq \varepsilon$,存在 $\delta > 0$、$M > 1$,使得对任意实值 $\varphi_0 \in \mathcal{S}_r(\mathbb{R}^n)$ 满足 $\|\varphi_0\|_{\dot{H}^s \cap \dot{H}^k(\mathbb{R}^n)} < \delta / M$,存在 $T_\varepsilon \in [1-\delta, 1+\delta]$ 及唯一径向解 $v \in C^\infty([0,T_\varepsilon) \times \mathbb{R}^n)$ 满足 (1.16),且可分解为

$$ v(t,x) = v^{T_\varepsilon}_\varepsilon(t,x) + \frac{1}{\sqrt{T_\varepsilon - t}} \, \varphi\left( \log\frac{T_\varepsilon}{T_\varepsilon - t}, \frac{x}{\sqrt{T_\varepsilon - t}} \right), $$

其中 $\varphi \in C^\infty([0,\infty) \times \mathbb{R}^n)$ 径向对称,且对任意 $r \in [s,k]$,有 $\|\varphi(-\log(T_\varepsilon - t) + \log T_\varepsilon, \cdot)\|_{\dot{H}^r(\mathbb{R}^n)} \lesssim \delta (T_\varepsilon - t)^\omega$。

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

核心障碍在于:当 $\varepsilon \neq 0$ 时,非线性项 $w_\varepsilon(u) w'_\varepsilon(u)$ 不再具有 $\sin u \cos u$ 的精确三角结构,导致线性化算子 $L_\varepsilon$ 的谱无法显式计算,且传统基于特殊函数的谱分析(如超几何函数解)完全失效。作者的真正创新点是构建一个参数化、几何自洽的扰动框架,将 $\varepsilon$ 视为小参数,而非试图求解新方程,而是将问题转化为对已知 $\varepsilon = 0$ 解的稳定性分析。具体技术分三步:(1)维度提升与相似坐标变换:通过变量替换 $\tilde v(t,r) = r^{-1} u(t,r)$,将 $d$ 维径向方程 (1.10) 提升为 $n = d+2$ 维的径向热方程 (1.16),其非线性项 $\frac{n-3}{|x|^3} (|x| v - w_\varepsilon(|x|v) w'_\varepsilon(|x|v))$ 在 $|x| \to 0$ 处光滑(因 $w_\varepsilon(z) = z - \frac{z^3}{6} + O(\varepsilon z^3)$),从而避免了原方程在 $r=0$ 处的奇异性。此提升使线性部分成为标准热算子 $\partial_t - \Delta_x$,便于应用现代半群理论。(2)线性化算子的谱扰动分析:在相似坐标 $(\tau, y)$ 下,静态解满足 $0 = \widetilde{L} \psi_\varepsilon + N_\varepsilon(\psi_\varepsilon)$。作者令 $\psi_\varepsilon = \psi_0 + \phi_\varepsilon$,其中 $\psi_0(y) = |y|^{-1} \hat f_0(|y|)$ 是 $\varepsilon = 0$ 的解。代入后得到扰动方程 (2.6): $0 = \widetilde{L}_0 \phi_\varepsilon + V_\varepsilon(\psi_0) \phi_\varepsilon + \widetilde{N}_\varepsilon(\phi_\varepsilon) + R_\varepsilon(\psi_0)$,其中 $\widetilde{L}_0 = \widetilde{L} + V_0(\psi_0)$ 是 $\varepsilon = 0$ 时的精确线性化算子,而 $V_\varepsilon(\psi_0) = \frac{n-3}{|y|^2} (\eta'_\varepsilon(|y|\psi_0) - \eta'_0(|y|\psi_0))$、$R_\varepsilon(\psi_0)$ 和 $\widetilde{N}_\varepsilon$ 均为 $O(\varepsilon)$ 小量(由 $w_\varepsilon$ 的定义保证)。关键突破在于:作者利用 $\alpha$ 的周期性与零点条件,证明 $\|V_\varepsilon(\psi_0)\|_{L^\infty} \lesssim \varepsilon$,且 $R_\varepsilon(\psi_0)$、$\widetilde{N}_\varepsilon$ 在交集空间 $X^k_s = \dot{H}^s_r(\mathbb{R}^n) \cap \dot{H}^k_r(\mathbb{R}^n)$ 上满足 Lipschitz 连续性(引理 3.2),从而将存在性问题归约为 $\widetilde{L}_0$ 的可逆性与 Banach 不动点定理。(3)不稳定模态的鲁棒控制:稳定性分析的核心是处理线性化算子 $L_\varepsilon$ 的谱。已知 $L_0$ 有单重不稳定本征值 $\lambda = 1$(对应于爆破时间平移不变性),其余谱位于 $\{ \Re \lambda \leq -\sigma_0 < 0 \}$。作者采用 Neumann 级数 + Riesz 投影连续性 的双重扰动策略:首先,对任意 $\mu$ 满足 $\Re \mu > -\sigma_0/2$,证明 $(\mu - L_\varepsilon)^{-1} = (\mu - L_0)^{-1} \sum_{j=0}^\infty \left[ (L_\varepsilon - L_0)(\mu - L_0)^{-1} \right]^j$ 收敛(因 $\|L_\varepsilon - L_0\|_{X^k_s \to X^k_s} \lesssim \varepsilon$),从而 $\sigma(L_\varepsilon)$ 在 $\{ \Re \lambda \leq -\sigma_0/2 \}$ 外至多位于以 $\lambda = 1$ 为中心的固定紧集内;其次,在该紧集内,利用 Riesz 投影 $P_\varepsilon$ 关于 $\varepsilon$ 的 Lipschitz 连续性([16, Lemma 4.10])及 $\dim \operatorname{Ran} P_0 = 1$,推出 $\dim \operatorname{Ran} P_\varepsilon = 1$,故 $\lambda = 1$ 是 $L_\varepsilon$ 唯一的不稳定本征值。此论证完全规避了对 $L_\varepsilon$ 的显式谱计算,仅依赖于 $\varepsilon$ 的小性与已知 $L_0$ 的谱隙结构。

4. 与之前工作的比较

与既有工作相比,本文的推进体现在三个根本维度:(1)几何假设的实质性放宽:此前所有稳定性结果(Biernat–Donninger–Schörkhuber [4,5], [1], Glogić–Kistner–Schörkhuber [15])均要求目标流形具有高度对称性:或是精确球面 $S^d$,或是特定旋转对称紧致流形(如 $S^1 \times S^{d-1}$)。本文首次证明,只要目标流形是 $S^d$ 的 $C^\infty$ 小扰动(即 $w_\varepsilon(u) = \sin u (1 + \varepsilon \alpha(u))$),稳定性即成立。这种扰动破坏了所有全局对称性(如 $SO(d+1)$ 不变性),且 $\alpha$ 的任意性意味着曲率张量 $\operatorname{Rm}_{S^d_\varepsilon}$ 不再为常数,而本文结论对此无额外要求。(2)技术路径的根本转换:以往工作依赖于 $\varepsilon = 0$ 时的显式结构:例如,$d=3$ 时 $\hat f_0$ 可用超几何函数表示,其线性化算子 $L_0$ 的谱可通过 Sturm–Liouville 理论精确刻画;$d=4,5,6$ 时虽无闭式解,但 [1] 仍依赖于 $S^d$ 的对称性导出精细的加权估计。本文则完全放弃显式解,转而发展一套纯泛函分析的扰动范式:将 $\varepsilon$ 视为参数,利用 $L_0$ 的谱隙与 $L_\varepsilon - L_0$ 的小性,通过算子扰动理论(而非 PDE 估计)控制谱。这使方法原则上可推广至更一般的几何扰动。(3)适用范围与代价:本文将稳定性结论从 $d=3$ 推广至 $d=3,4,5,6$,与 [1] 一致,但代价是初值正则性要求更严格:定理 1.3 中 $s$ 的上界为 $\frac{d}{2} + \frac{1}{2d+2}$,而 [1] 在 $d=3$ 时允许 $s < 2$;本文的 $k = d+3$ 也高于 [1] 的 $k = d+2$,反映了扰动引入的额外技术复杂性。此外,本文未处理 $d \geq 7$,因此时 $\widetilde{L}_0$ 的谱隙可能消失(见 [6]),且 $w_\varepsilon$ 的扰动可能导致 $\eta_\varepsilon$ 的高阶导数增长过快,破坏 $X^k_s$ 上的 Lipschitz 性。最后,本文明确指出(Remark 1.4),稳定性结论不保证扰动解始终位于单坐标卡 $u \in (0,\pi)$ 内;若剖面穿过 $u=\pi$(即南极大),需切换坐标卡,这与 $\varepsilon = 0$ 时的全局覆盖不同,体现了几何扰动带来的拓扑复杂性。从提供内容中无法确认本文是否改进了爆破速率或误差衰减率 $\omega$ 的具体数值,作者亦未显式展开与波映射(wave maps)文献 [9] 的技术细节比较。


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

1. 聚焦非线性薛定谔方程在外部势作用下的近似行波解

Approximate traveling waves for the focusing NLS with an external potential

作者: Luke Baker, Xuemei Li, Jason Murphy

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03915v1

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

摘要翻译

中文标题: 聚焦非线性薛定谔方程在外部势作用下的近似行波解

摘要翻译: 我们构造了受快速衰减外部势扰动的 $3d$ 聚焦三次非线性薛定谔方程(NLS)的近似行波解。作为应用,我们证明:对未扰动方程而言,在任意高于基态阈值的质量-能量下,存在一个正向全局、$H^1$ 有界且不散射的解。


2. 全抛物型Keller–Segel系统中有界性的显式Logistic阻尼判据

An Explicit Logistic Damping Criterion for Boundedness in a Fully Parabolic Keller–Segel System

作者: Jie Jiang

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03869v1

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

摘要翻译

中文标题: 全抛物型Keller–Segel系统中有界性的显式Logistic阻尼判据

摘要翻译: 我们研究定义在有界光滑凸区域上的全抛物型Keller–Segel系统

$$ u_t=Δu-χ\nabla\!\cdot(u\nabla v)+λu-μu^2, \qquad τv_t=Δv-v+u $$

对任意固定的 $τ>0$,我们证明了条件

$$ μ>\frac{Nχ}{4} $$

可保证整体解的存在性及关于时间的一致有界性。该系数显式的充分条件与 $τ$ 无关,且不涉及任何嵌入常数或极大正则性常数。据我们所知,这是首个在任意空间维数下对所有 $τ>0$ 均保持不变的系数显式有界性判据。

证明基于一个新颖的辅助比较函数

$$ Y_τ =u+\frac{χτ}{2}|\nabla v|^2-(τ-1)Δv, $$

该函数对每个 $τ>0$ 满足一个封闭的标量抛物型不等式。当 $τ\ge1$ 时,该不等式直接导出逐点比较及 $u$ 的显式上界;当 $0<τ<1$ 时,它给出 $v$ 的一致上界。随后对变换 $z=e^{-χv/2}$ 应用抛物型挤压论证,得到 $v$ 的一致Hölder界。再结合Hölder–Sobolev插值及加权极大 $L^p$-正则性,可获得对足够大的 $p$ 成立的 $u$ 的 $L^p$-界,标准抛物型正则化最终完成论证。


3. 矩形区域上磁性 Neumann Laplacian 的谱不等式

Spectral inequalities for magnetic Neumann Laplacian on rectangles

作者: Tuyen Vu

类别: math.SP, math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03843v1

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

摘要翻译

中文标题: 矩形区域上磁性 Neumann Laplacian 的谱不等式

摘要翻译: 我们考虑矩形区域上的磁性 Neumann Laplacian,其磁场为非均匀场。我们证明:在固定面积的所有矩形中,正方形是弱磁场下最低特征值的局部极小元。我们猜想:在面积约束或周长约束下,正方形均为该最低特征值的全局极小元;并部分验证该猜想,即建立了主特征值的上下界估计。


4. 含时依赖阻尼与质量项的半线性波动方程的爆破准则与寿命估计

Blow-up criteria and lifespan estimates for semilinear wave equations with time-dependent damping and mass

作者: Wenhui Chen, Mohamed Ali Hamza

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03713v1

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

摘要翻译

中文标题: 含时依赖阻尼与质量项的半线性波动方程的爆破准则与寿命估计

摘要翻译: 我们考虑具时依赖阻尼与质量项及导数型非线性项的半线性波动方程。通过应用 Liouville 变换并求解一个从无穷远处出发的 Volterra 积分方程,我们在不借助线性传播子显式表示的情况下构造了伴随方程的一个正精确解。该解被用于推导具有有限传播速度的能量解的爆破准则,并给出其寿命的上界估计。若阻尼系数的原函数至多以对数速率增长,则我们得到完整的平移 Glassey 范围,包括临界指数,以及相应的多项式型与指数型寿命估计。该结果独立于尺度不变判别式 $\frac{\mu^2}{4} - m^2$ 的符号,因而涵盖质量主导情形。此外,当阻尼系数的非可积振荡摄动项之贡献被相应质量项所抵消时,该结果亦适用。


5. Morrey 问题在 $\mathbb{R}^{2\times m}$ 中的一个解

A solution to Morrey’s problem in $\mathbb{R}^{2\times m}$

作者: Gabriele Cassese

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03488v1

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

摘要翻译

中文标题: Morrey 问题在 $\mathbb{R}^{2\times m}$ 中的一个解

摘要翻译: 我们构造了齐次的秩一凸被积函数 $F\colon \mathbb{R}^{2\times m} \to \mathbb R$,当 $m$ 充分大时,该函数处处不是拟凸的。


6. 从 $\mathrm{RCD}$ 空间到 $\mathrm{CAT}(0)$ 空间的调和映照热流的 Lipschitz 正则性

Lipschitz regularity of harmonic map heat flows from $\mathrm{RCD}$ spaces into $\mathrm{CAT}(0)$ spaces

作者: Bang-Xian Han, Hui-Chun Zhang, Xi-Ping Zhu

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

发布日期: 2026-08-04

arXiv ID: 2608.03442v1

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

摘要翻译

中文标题: 从 $\mathrm{RCD}$ 空间到 $\mathrm{CAT}(0)$ 空间的调和映照热流的 Lipschitz 正则性

摘要翻译: 即使目标空间具有奇异性,进入完备 $\mathrm{CAT}(0)$ 空间的调和映照热流仍允许变分构造,但当源空间亦非光滑时,其正则性理解尚不充分。我们在初始数据具有有界像的假设下,证明了从有限维 $\mathrm{RCD}(K,N)$ 空间到完备 $\mathrm{CAT}(0)$ 目标空间的 Korevaar–Schoen 能量的 $\mathrm{EVI}$ 梯度流具有正时间局部 Lipschitz 正则性。所得代表元在空间与时间变量上联合局部 Lipschitz。证明结合了度量 Hamilton–Jacobi 论证与施加于 $\mathrm{RCD}$ 源空间适当时间截面上的椭圆接触估计。


7. 混合维数边界的泊松-狄利克雷问题的可解性

Solvability of the Poisson-Dirichlet problem in domains with boundaries of mixed dimension

作者: Léo Mandô

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03373v1

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

摘要翻译

中文标题: 混合维数边界的泊松-狄利克雷问题的可解性

摘要翻译: 我们将齐次狄利克雷问题的 $L^p$-可解性的若干刻画推广至定义在 $\mathbb{R}^n$ 中一大类开集 $Ω$ 上的(可能退化的)椭圆算子 $L=-\textrm{div}(wA\nabla)$。该框架不仅涵盖 Lipschitz 区域上的具一致椭圆性的算子,亦包含 Caffarelli-Sylvestre 型算子,且边界 $\partialΩ$ 不必是 $(n-1)$-维的,例如。我们证明:对 $1<p<\infty$,如下陈述彼此等价:(i) 齐次狄利克雷问题在 $L^p$ 意义下可解;(ii) 非切向极大函数与非切向迹满足 $L^p$ 估计;(iii) 对所有 $f\in L^p(\partialΩ)$,有 $\|N(u_f)\|_{L^p(\partialΩ)}\lesssim \|f\|_{L^p(\partialΩ)}$,其中 $u_f$ 是以 $f$ 为边界数据的解;(iv) 层势算子 $\mathcal{S}:L^p(\partialΩ)\to L^p(\partialΩ)$ 在 $L^p(\partialΩ)$ 上有界;(v) 单层势 $\mathcal{S}$ 的法向导数 $\partial_\nu \mathcal{S}$ 在 $L^p(\partialΩ)$ 上可逆;(vi) 双层势 $\mathcal{D}$ 在 $L^p(\partialΩ)$ 上有界。此外,我们给出一个适用于任意 $p\in (1,\infty)$ 的充分条件,即权重 $w$ 满足 Muckenhoupt $A_2$ 条件,且系数矩阵 $A$ 满足某种 Carleson 测度条件;该条件保证了上述等价性成立。


8. 长时区间上拟周期函数的Strichartz估计:一种算术方法

Strichartz estimates for quasi-periodic functions on long time intervals: An arithmetic approach

作者: Kotaro Inami

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

发布日期: 2026-08-04

arXiv ID: 2608.03194v1

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

摘要翻译

中文标题: 长时区间上拟周期函数的Strichartz估计:一种算术方法

摘要翻译: 我们研究具拟周期初值的一维Schrödinger方程的长时间Strichartz估计。对具有代数频率比的双频率初值,我们观察到线性Schrödinger演化的行为依赖于该频率比的代数次数。利用这一观察,我们在长时区间上改进了Strichartz估计。我们还证明了一个端点 $L^4$ Strichartz估计。我们的证明使用了Roth型Diophantine不等式以及Parsell–Vinogradov系统的Vinogradov型均值估计。


9. 关于相对论性动能薛德尔估计的一则注记

A note on relativistic kinetic Schauder estimates

作者: Weinan Wang

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03126v1

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

摘要翻译

中文标题: 关于相对论性动能薛德尔估计的一则注记

摘要翻译: 在本注记中,受 \cite{DW26} 启发,我们为具有相对论性输运的线性动能方程构造了仅关于动量的薛德尔估计(momentum-only Schauder estimates)的反例。该构造利用了非线性动量-速度映射与一对固定的光滑扩散系数和漂移系数之间的结构相容性:经共轭变换后,扩散项化为拉普拉斯算子,且所有诱导的一阶项完全抵消。因此,该障碍由一个固定的算子实现——该算子在每个有界动量集上一致椭圆,而非通过一族依赖于频率的系数实现。我们进一步证明,即使在外力为零的情形下,该失效现象依然存在:我们构造了具有一致正性的稳态解,其隐含的空间振荡被编码在一个有界的零阶系数中,而该系数关于动量的 Hölder 半范数趋于零。因此,\cite{HSTT} 中所采用的完整洛伦兹型 Hölder 控制一般不能被仅关于动量的 Hölder 控制所替代。


10. Sobolev 空间上的分数阶半 Dirac 算子

The Fractional Half Dirac Operator Over Sobolev Spaces

作者: Dejenie A. Lakew

类别: math.AP, math.FA

发布日期: 2026-08-03

arXiv ID: 2608.02871v1

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

摘要翻译

中文标题: Sobolev 空间上的分数阶半 Dirac 算子

摘要翻译: 设 $\Omega$ 为 $n$ 维欧几里得空间中的有界光滑区域,并考虑半阶 Dirac 算子及其幂的边值问题:其中 $f$ 为解在区域内的半阶导数值,$g$ 为解在边界上的迹值。我们将证明:该问题的解是一个广义函数 $u$,属于半阶 Hilbert 空间,且该空间可正交分解为两部分之直和:一部分是由分数阶半 Dirac 算子在区域内的取值所诱导的区域积分解;另一部分是由解在边界(或迹)上的取值所诱导的边界积分解。我们在 Sobolev 空间中对子空间及其正交子空间中的函数分别采用两种不同的符号来表示正交直和。


11. 具并发输入与状态时滞的一阶偏积分-微分方程的稳定化

Stabilization of First-Order Partial Integro-Differential Equations with Concurrent Input and State Delays

作者: Sanguan Zhong, Jie Qi

类别: eess.SY, math.AP, math.OC

发布日期: 2026-08-03

arXiv ID: 2608.02851v1

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

摘要翻译

中文标题: 具并发输入与状态时滞的一阶偏积分-微分方程的稳定化

摘要翻译: 本文研究一类具并发输入与状态时滞的一阶双曲型偏积分-微分方程(PIDE)的边界稳定化问题。这两类时滞的同时存在使得控制设计复杂化,尤其在输入时滞较大时,需预测更多状态信息。本文基于反步法构造了一种边界控制器,以实现系统稳定化与时滞补偿。该设计依赖于两个仿射 Volterra 变换,其核函数同时包含 Fredholm 型与 Volterra 型积分项,从而导出四组 PIDE 核方程。为建立这些核方程的适定性,将核定义域沿特征线划分为有限个三角子区域,并在前一子区域中逐次构造核解至下一子区域。本文建立了所得闭环系统的有限时间稳定性。数值模拟结果验证了所提出控制器的有效性。


不相关论文(基本信息)

1. 电磁反问题的降阶建模:分层介质基准

Reduced-order modeling for electromagnetic inverse problems: a layered medium benchmark

作者: Konstantinos Alexopoulos, Josselin Garnier

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

发布日期: 2026-08-04

arXiv ID: 2608.03996v1

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

摘要: We study reduced-order modeling for inverse problems in layered media, focusing on the recovery of impedance profiles from time-domain measurements. Using the Goupillaud structure, we formulate the forward problem as a discrete dynamical system and introduce a ROM-based objective defined at the level of reduced operators. Through numerical experiments on a 5-layer medium under various structured perturbations, we compare the ROM-based objective with a classical data misfit. The results show that the ROM-based inversion provides a more favorable reconstruction in the clean setting and in certain structured perturbation regimes, while remaining consistently competitive with the classical approach across all cases considered.


2. 时空介质中的折射定律

Refraction laws in spatio-temporal media

作者: Cristian E. Gutiérrez, Eric Stachura

类别: math.AP, math-ph, physics.optics

发布日期: 2026-08-04

arXiv ID: 2608.03995v1

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

摘要: We study the time-dependent Maxwell system, formulated in the sense of distributions, for electromagnetic waves propagating through media with temporal and spatial material interfaces. Under explicit trace and regularity assumptions on the permittivity and permeability, we derive the jump conditions produced by temporal discontinuities and by subsequent spatial interfaces. These conditions are used to obtain generalized Snell laws for the wave vectors generated by temporal splitting, reflection, and transmission. We also derive the associated amplitude equations and organize them as a finite-dimensional linear system for a space-time slab geometry. Finally, we provide an explicit oblique-incidence example. Our approach does not impose a smooth-field ansatz and allows material parameters that need not be constant away from the interfaces.


3. 关于 $q$-阶对偶曲率积分的 Brunn-Minkowski 不等式(其中 $q>n$)

On the Brunn-Minkowski inequality for $q$-th dual quermassintegrals with $q>n$

作者: Haizhong Li, Yao Wan, Yijia Zhang

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

发布日期: 2026-08-04

arXiv ID: 2608.03949v1

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

摘要: In this paper, we study the Brunn-Minkowski inequality for $q$-th dual quermassintegrals with $q>n$. This problem was recently posed by Sadovsky and Zhang. First, by a second variation argument and a dimension reduction construction, we show that the inequality fails for arbitrary convex bodies when $q>n$, and fails even in the origin-symmetric class when $q>n+2$. Secondly, we prove the endpoint case $q=n+2$ for origin-symmetric convex bodies via Hadwiger’s inequality for the polar moment of inertia. Finally, for unconditional convex bodies, we establish the inequality in the full range $0<q\le n+1$ by using a singular weighted Reilly formula and a coordinate-slice Hardy inequality. As applications, we derive several uniqueness results for the corresponding dual curvature measures.


4. 具有线性增长的拟凸泛函的严格 BV 松弛能量的显式公式

Explicit formulas of strict BV relaxed energies for polyconvex functionals with linear growth

作者: Domanico Mucci, Riccardo Scala

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03821v1

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

摘要: We consider polyconvex functionals defined on vector valued functions, the model case being the graph area functional, and we analyze the relaxed energy with respect to the strict convergence in $BV$. In several cases where the relaxed area is known, by exploiting a continuity result by Reshetnyak we are able to find an explicit formula for wide classes of integrands with linear growth in the minors of the gradient. We preliminarily extend to the $BV$ setting a continuity property observed by Acerbi–Dal Maso in the Sobolev case. Finally, partial results concerning the relaxation of the vortex map with respect to the $L^1$ convergence are obtained.


5. 关于反极大值原理的注记

Remarks on the antimaximum principle

作者: Vladimir Bobkov

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03460v1

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

摘要: We present several observations on the antimaximum principle (AMP) for the model problem $-Δ_p u = λ|u|^{p-2} u + f$ in a bounded smooth domain $Ω$, subject to the zero Dirichlet boundary conditions, and where the source function $f$ is nontrivial, nonnegative, and sufficiently regular. Denote by $λ_f$ the endpoint of validity of the AMP, so that every solution of the problem is negative in $Ω$ for any $λ\in (λ_1,λ_f)$. Our discussion covers the following aspects: identification of a class of sources over which the AMP is uniform, lower semicontinuity of the map $f \mapsto λ_f$, bounds on $λ_f$, the nonexistence of negative solutions for sufficiently large $λ$ (extended AMP), the anticomparison principle, and the weakening of the source regularity from the Lebesgue to Morrey spaces. Some of the results are stated only in the linear case $p=2$. As a part of the discussion, we provide a few related open problems.


6. 临界指数双调和Choquard方程多泡解的集中性、局部唯一性与Morse指标

Concentration, Local Uniqueness, and Morse Index of Multi-bubble Solutions for a Critical Exponential Biharmonic Choquard Equation

作者: Wenjing Chen, Shengbing Deng

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03321v1

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

摘要: Let $Ω\subset\mathbb R^4$ be a bounded domain of class $C^6$, let $0<α<4$, and let $K\in C^4(\overlineΩ)$ be positive. We study the Navier problem

$$ Δ^2u=\varepsilon^{8-α}K(x)e^{u(x)} \left(\int_Ω\frac{K(y)e^{u(y)}}{|x-y|^α}\,dy\right), \qquad u=Δu=0\quad\text{on }\partialΩ. $$

Let $G$ be the Navier Green function, let $H$ be its regular part, and put $M_α=8π^2(8-α)$. The concentration points are governed by

$$ \mathcal F_m(\boldsymbolξ) =\sum_{i=1}^m \left[\log K(ξ_i)+\frac{M_α}{2}H(ξ_i,ξ_i)\right] +M_α\sum_{i<j}G(ξ_i,ξ_j). $$

Every $C^1$-stable critical point of $\mathcal F_m$ produces a positive $m$-bubble solution whose scales are of order $\varepsilon^{-1}$ and whose nonlinear source converges to $M_α\sum_iδ_{ξ_i^*}$. If the critical point is nondegenerate, the corresponding $m$-bubble solution is locally unique, modulo permutations, in a fixed scaled modulation neighborhood. The linearized operator is nondegenerate on $H^2(Ω)\cap H_0^1(Ω)$, and

$$ \operatorname{ind}(u_\varepsilon) =m+\operatorname{ind}\!\left(-D^2\mathcal F_m(\boldsymbolξ^*)\right). $$

A critical four-dimensional capacity controls the scale directions. The dilation block of the reduced Hessian is positive and equals $8π^2b_α^2|\log\varepsilon|^{-1}I_m+o(|\log\varepsilon|^{-1})$, where $b_α=(8-α)/2$.


7. 三维Heisenberg Bieberbach流形上次拉普拉斯算子与Reeb向量场的联合谱

Joint spectrum of the sub-Laplacian and the Reeb vector field on three-dimensional Heisenberg Bieberbach manifolds

作者: Yoshiaki Suzuki, Yuya Takeuchi

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

发布日期: 2026-08-04

arXiv ID: 2608.03265v1

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

摘要: A Heisenberg Bieberbach manifold is a compact quotient of the Heisenberg group by a discrete torsion-free subgroup of its pseudo-Hermitian automorphism group. In this paper, we study the joint spectrum of the sub-Laplacian and the Reeb vector field on three-dimensional Heisenberg Bieberbach manifolds. In particular, we explicitly compute the multiplicities of the joint spectrum using theta functions and the character theory of finite groups.


8. 各向异性情形下电导率成像的稳定性及其在深度Caldeón方法中的应用

Stability of Electrical Impedance Tomography with Anisotropies and its Application to the Deep Caldeón Method

作者: Tianhao Hu, Bangti Jin, Yiran Wang

类别: math.AP

发布日期: 2026-08-04

arXiv ID: 2608.03053v1

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

摘要: In this work, we establish new conditional Lipschitz stability results for electrical impedance tomography (EIT) with anisotropies, of recovering the conductivity in a conformal class of a known anisotropic conductivity in both two- and multi-dimensional cases. Then we employ the stability theory to understand the property of the deep Calderón method, one deep learning-based technique for image reconstruction in EIT that has shown promising empirical results, but still lacks theoretical underpinnings. Specifically, we relate the stability theory to the robustness of the method with the proper choice of the training data, and present numerical results in two-dimension to complement the theoretical analysis.


9. 各向异性Maxwell系统中介电常数与磁导率的通用恢复

Generic Recovery of Permittivity and Permeability in Anisotropic Maxwell Systems

作者: Antonio Cocan, Maarten V. de Hoop, Joonas Ilmavirta, Matti Lassas, Anthony Várilly-Alvarado

类别: math.AP, math.AG

发布日期: 2026-08-03

arXiv ID: 2608.02900v1

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

摘要: We study the inverse problem of recovering the constitutive tensors of a homogeneous anisotropic electromagnetic medium without magnetoelectric coupling (non-chiral) from its Fresnel surface, the characteristic variety of Maxwell’s equations governing electromagnetic wave propagation. For known isotropic permeability, normalized to $μ= I$, we prove that the Fresnel surface uniquely determines the permittivity tensor $\varepsilon$, and that the associated Fresnel polynomial is reducible precisely when $\varepsilon$ has a repeated eigenvalue. For general, positive-definite symmetric tensors $(\varepsilon,μ)$, we prove that the Fresnel polynomial is generically irreducible over $\mathbb{C}$ and we identify the natural gauge symmetry under which it is invariant. Using geometric invariant theory, a powerful tool of modern algebraic geometry, we construct an affine quotient of the parameter space by this gauge action and prove that the induced Fresnel-polynomial map is birational onto its image. We deduce that, outside a proper real algebraic exceptional set, the real Fresnel surface determines $(\varepsilon,μ)$ up to gauge. This establishes generic uniqueness for the inverse problem and, to our knowledge, provides a new application of affine geometric invariant theory to gauge freedom in a PDE inverse problem.