2025-11-12T02:46:09.412369

Critical curve for weakly coupled system of semilinear Euler-Poisson-Darboux-Tricomi equations

Li, Guo
This paper investigates a weakly coupled system of semilinear Euler-Poisson-Darboux-Tricomi equations (EPDTS) with power-type nonlinear terms. More precisely, in the case where the damping terms dominate over the mass terms, the critical curve in the $p-q$ plane that delineates the threshold between global existence and blow-up for the EPDTS is given by \begin{equation*} Γ_m(n,p,q,β_1,β_2)=0, \end{equation*} where $Γ_m$ is defined by (\ref{gammam}). Through the construction of new test functions, the blow-up problem is addressed when $Γ_m(n,p,q,β_1,β_2)\geq0$. Based on the $(L^1\cap L^2)-L^2$ estimates of the solution to the corresponding linear equation established in our previous work \cite{LiGuo2025}, we derive the global existence of solutions with small initial data when $Γ_m(n,p,q,β_1,β_2)<0$, provided that the damping terms prevail over the mass terms.
academic

Critical curve for weakly coupled system of semilinear Euler-Poisson-Darboux-Tricomi equations

基本信息

  • 论文ID: 2511.08084
  • 标题: Critical curve for weakly coupled system of semilinear Euler-Poisson-Darboux-Tricomi equations
  • 作者: Yuequn Li, Fei Guo
  • 分类: math.AP (Analysis of PDEs)
  • 发表时间: 2025年11月11日
  • 论文链接: https://arxiv.org/abs/2511.08084

摘要

本文研究带幂型非线性项的半线性Euler-Poisson-Darboux-Tricomi方程(EPDTS)的弱耦合系统。在阻尼项主导质量项的情况下，本文给出了在p−qp-q平面上区分全局存在性与爆破现象的临界曲线Γm(n,p,q,β1,β2)=0\Gamma_m(n,p,q,\beta_1,\beta_2)=0。通过构造新的测试函数，本文处理了Γm(n,p,q,β1,β2)≥0\Gamma_m(n,p,q,\beta_1,\beta_2)\geq0时的爆破问题。基于作者之前工作中建立的相应线性方程解的(L1∩L2)−L2(L^1\cap L^2)-L^2估计，本文导出了当Γm(n,p,q,β1,β2)<0\Gamma_m(n,p,q,\beta_1,\beta_2)<0且阻尼项占优时小初值解的全局存在性。

研究背景与动机

研究问题

本文研究以下弱耦合半线性Euler-Poisson-Darboux-Tricomi方程组：

\partial_t^2 u - t^{2m}\Delta u + \frac{\mu_1}{t}\partial_t u + \frac{\nu_1^2}{t^2}u = |v|^p, & t>1, x\in\mathbb{R}^n,\\ \partial_t^2 v - t^{2m}\Delta v + \frac{\mu_2}{t}\partial_t v + \frac{\nu_2^2}{t^2}v = |u|^q, & t>1, x\in\mathbb{R}^n,\\ (u,\partial_t u, v, \partial_t v)(1,x) = (u_0,u_1,v_0,v_1)(x), & x\in\mathbb{R}^n \end{cases}$$ 其中$m>-1$，$n\geq 1$为空间维数，$\mu_i, \nu_i^2\geq 0$为非负参数，$p,q>1$。 ### 研究重要性 1. **理论意义**：该系统推广了经典的波方程理论，涉及Strauss指数和Fujita指数的相互作用 2. **数学结构**：Gellerstedt算子$\partial_t^2 - t^{2m}\Delta$的存在使得阻尼项$\frac{\mu}{t}\partial_t u$和质量项$\frac{\nu^2}{t^2}u$具有相同的尺度不变性 3. **临界现象**：确定全局存在与爆破的临界曲线是偏微分方程理论中的核心问题 ### 现有方法的局限性 1. **单方程结果不能直接推广**：弱耦合系统的临界曲线不是简单的矩形区域 2. **测试函数构造**：Gellerstedt算子与标准波算子的差异要求构造依赖于参数$m$的新测试函数 3. **参数$\delta$的作用**：参数$\delta_i = (\mu_i-1)^2 - 4\nu_i^2$刻画阻尼项与质量项的相互作用，其大小决定了系统的抛物型或双曲型特征 ### 研究动机 1. 推广作者之前关于单方程的工作[17]到耦合系统 2. 揭示参数$\mu_i, \nu_i, m$之间的复杂相互作用 3. 建立完整的全局存在性与爆破理论 ## 核心贡献 1. **建立临界曲线**：证明了当$\delta_i$足够大时，临界曲线为 $$\Gamma_m(n,p,q,\beta_1,\beta_2) := \max\left\{\frac{p+1}{pq-1}-\frac{\beta_1-1}{2}, \frac{q+1}{pq-1}-\frac{\beta_2-1}{2}\right\} - \frac{(m+1)n}{2} = 0$$ 其中$\beta_i = \frac{\mu_i+1-\sqrt{\delta_i}}{2}$ 2. **构造新测试函数**：针对Gellerstedt算子设计了依赖于参数$m$的新测试函数，证明了$\Gamma_m\geq 0$时的爆破结果 3. **全局存在性理论**：基于$(L^1\cap L^2)-L^2$估计，在$\Gamma_m<0$时建立了小初值解的全局存在性，分为高正则性($\sigma\geq 1$)和低正则性($0<\sigma<1$)两种情况，每种情况又细分为三个子情形 4. **非平凡推广**：证明了耦合系统的结果不是单方程的简单推广，临界区域存在"偏移"现象，体现了两个阻尼项和两个质量项之间的复杂平衡 ## 方法详解 ### 任务定义 研究系统(1.1)的Cauchy问题，确定参数$(p,q)$在何种条件下解全局存在或有限时间爆破。 **能量解定义**：$(u,v)$是系统(1.1)在$[1,T)$上的能量解，如果： - 正则性：$(u,v)\in (C([1,T);H^1)\cap C^1([1,T);L^2))^2$ - 支集性质：$\text{supp}(u,v)(t,\cdot)\subset B_{\varphi_m(t)-\varphi_m(1)+M}$，其中$\varphi_m(t) = \frac{t^{m+1}}{m+1}$ - 满足积分等式(2.6)和(2.7) ### 爆破理论的方法架构 #### 测试函数构造 对于爆破问题，关键创新是构造新的测试函数： $$\psi_{d,R}(t,x) = \lambda\left(\frac{t-1}{d}\right)\varphi\left(\frac{x}{R^{m+1}}\right)$$ 其中$\lambda(t), \varphi(x)\in C_0^\infty$满足特殊性质(3.3)。 **关键技术点**： 1. 空间截断依赖于$R^{m+1}$而非$R$，这是适应Gellerstedt算子的核心 2. 利用$h_i(t) = t^{\beta_i}$满足的ODE性质(3.1) 3. 将方程改写为散度形式(3.2) #### 爆破证明策略 1. 定义积分泛函$I_{d,R}, J_{d,R}$ 2. 通过Hölder不等式和支集性质估计各项 3. 利用条件(2.9)保证积分指数的可积性 4. 导出矛盾：若$T=\infty$，则在$\Gamma_m\geq 0$时必有$u\equiv 0$或$v\equiv 0$，与初值条件(2.8)矛盾 **关键不等式**： $$I_{R,R} \lesssim (R+1)^{-2-\frac{2}{q}+(1-\frac{1}{pq})((m+1)n+\beta_1+1)}I_{R,R}^{\frac{1}{pq}}$$ $$J_{R,R} \lesssim (R+1)^{-2-\frac{2}{p}+(1-\frac{1}{pq})((m+1)n+\beta_2+1)}J_{R,R}^{\frac{1}{pq}}$$ ### 全局存在性的方法架构 #### 线性理论基础 利用作者之前工作[17]中建立的线性方程 $$\partial_t^2 u - t^{2m}\Delta u + \frac{\mu}{t}\partial_t u + \frac{\nu^2}{t^2}u = 0$$ 的$(L^1\cap L^2)-L^2$估计(Propositions 4.1-4.2)。 **核心估计**：对于$\kappa\in[0,\sigma]$， $$\|u(t,\cdot)\|_{\dot{H}^\kappa} \lesssim \|(f,g)\|_{D_\kappa} \times \begin{cases} t^{-\frac{\mu+m}{2}}, & \text{if } \kappa > \frac{\sqrt{\delta}}{2(m+1)} + \frac{1}{2} - \frac{n}{2}\\ t^{-\frac{\mu+m}{2}}(1+\log t)^{1/2}, & \text{if } \kappa = \frac{\sqrt{\delta}}{2(m+1)} + \frac{1}{2} - \frac{n}{2}\\ t^{-(m+1)(\kappa+\frac{n}{2})+\frac{\sqrt{\delta}-\mu+1}{2}}, & \text{if } \kappa < \frac{\sqrt{\delta}}{2(m+1)} + \frac{1}{2} - \frac{n}{2} \end{cases}$$ #### 非线性理论：不动点方法 通过Duhamel原理将解表示为： $$u(t,x) = E_0^{\mu_1,\nu_1}(t,1,x)*u_0 + E_1^{\mu_1,\nu_1}(t,1,x)*u_1 + \int_1^t E_1^{\mu_1,\nu_1}(t,\tau,x)*|v(\tau,x)|^p d\tau$$ **函数空间设计**（高正则性$\sigma\geq 1$）： $$X(T) = \{(u,v)\in (C([1,T];H^\sigma)\cap C^1([1,T];H^{\sigma-1}))^2: \text{supp}(u,v)\subset B_{\varphi_m(t)-\varphi_m(1)+M}\}$$ 配备范数 $$\|(u,v)\|_{X(T)} = \sup_{t\in[1,T]}(t^{-\alpha_1}M_1(t,u) + t^{-\alpha_2}M_2(t,v))$$ 其中$M_i$包含$L^2, \dot{H}^\sigma$范数及其时间导数，权重$\alpha_i$由条件$p\lessgtr \tilde{p}, q\lessgtr \tilde{q}$决定。 **关键参数**： $$\tilde{p} = \frac{(m+1)n+\beta_1+1}{(m+1)n+\beta_2-1}, \quad \tilde{q} = \frac{(m+1)n+\beta_2+1}{(m+1)n+\beta_1-1}$$ #### 三种情形的处理 1. **$p>\tilde{p}, q>\tilde{q}$**（Theorem 2.2）：$\alpha_1=\alpha_2=0$，标准压缩映射 2. **$p\leq\tilde{p}, q>\tilde{q}$**（Theorem 2.3）：$\alpha_1>0, \alpha_2=0$，需要额外技术条件(2.27) 3. **$p>\tilde{p}, q\leq\tilde{q}$**（Theorem 2.4）：对偶于情形2 ### 技术创新点 1. **测试函数的$m$-依赖性**：空间截断$R^{m+1}$精确捕捉了Gellerstedt算子的传播速度 2. **加权范数的精细设计**：权重$t^{-\alpha_i}$的引入使得在$p\leq\tilde{p}$时仍能建立全局存在性 3. **对数因子的处理**：当$\delta_i = (m+1)^2(n+2\sigma-1)^2$时，估计中出现$(1+\log t)^{1/2}$因子，需要特殊的可积性分析 4. **初值正则性的分层**： - $\sigma\geq 1$：需要控制$\partial_t u$的范数 - $0<\sigma<1$：不需要时间导数的估计，简化了函数空间 ## 实验设置 **注**：本文为纯数学理论研究，不涉及数值实验。所有结果均通过严格的数学证明获得。 ### 理论验证 论文通过构造具体例子验证理论条件的非空性： **Example 2.1**（爆破）：$m=0.1, n=2, \mu_1=3, \nu_1=0.25, \mu_2=0.5, \nu_2=0.125$ - 计算得$\delta_1=\frac{15}{4}, \delta_2=\frac{3}{16}, \beta_1=\frac{2-\sqrt{15}}{4}, \beta_2=\frac{3}{4}-\frac{\sqrt{3}}{8}$ - 选择$p=2.1>1.172, q=2.2>0.687$ - 得到$\Gamma_m\approx 0.017>0$，解爆破 **Example 2.3**（全局存在）：$\sigma=1.1, m=0.4, n=3, \mu_1=15, \nu_1=6, \mu_2=8.7, \nu_2=0.25$ - 验证$\delta_1=52, \delta_2=59.04 > 34.5744 = (m+1)^2(n+2\sigma-1)^2$ - 对任意$2<p\leq 3.5, 2<q\leq 3.5$，系统有唯一全局解 ## 实验结果 ### 主要理论结果 #### 爆破结果（Theorem 2.1） **条件**： - $\mu_i>1$（当$-1<m<0$）或$\mu_i>0$（当$m\geq 0$） - $\delta_1,\delta_2>0$ - 初值满足(2.8) - $p,q$满足(2.9)和$\Gamma_m(n,p,q,\beta_1,\beta_2)\geq 0$ **结论**：能量解在有限时间爆破 **意义**： 1. 爆破区域不是矩形$\{1<p\leq p_F^1, 1<q\leq p_F^2\}$ 2. 存在$p$或$q$大于相应Fujita指数但仍爆破的情形 3. 当$\beta_1\neq\beta_2$时，爆破区域发生"偏移" #### 全局存在性结果 **高正则性（$\sigma\geq 1$）**： 1. **Theorem 2.2**（$p>\tilde{p}, q>\tilde{q}$）： - 要求$\delta_1,\delta_2 \geq (m+1)^2(n+2\sigma-1)^2$ - 解满足衰减估计： $$\|u(t,\cdot)\|_{L^2} \lesssim t^{-\frac{(m+1)n}{2}+\frac{\sqrt{\delta_1}-\mu_1+1}{2}}\|(u_0,u_1,v_0,v_1)\|$$ $$\|u(t,\cdot)\|_{\dot{H}^\sigma} \lesssim t^{-\frac{(m+1)(\sigma+n/2)-\sqrt{\delta_1}-\mu_1+1}{2}}\ell_1(t)\|(u_0,u_1,v_0,v_1)\|$$ 2. **Theorem 2.3**（$p\leq\tilde{p}, q>\tilde{q}$）： - 额外要求$\frac{q+1}{pq-1} < \frac{(m+1)n+\beta_2-1}{2}$ - $u$的估计带额外增长因子$t^{\alpha_1}$，$v$保持标准衰减 3. **Theorem 2.4**：对偶于Theorem 2.3 **低正则性（$0<\sigma<1$）**：Theorems 2.5-2.7对应于Theorems 2.2-2.4，但不需要控制时间导数 ### 关键发现 1. **临界曲线的非平凡性**（Remark 2.9）： - $\Gamma_m\geq 0 \Leftrightarrow p\leq\tilde{p} \text{ and } q\leq\tilde{q}$ - $\Gamma_m<0$可由三种情形实现：$(p>\tilde{p}, q>\tilde{q})$或$(p\leq\tilde{p}, q>\tilde{q})$或$(p>\tilde{p}, q\leq\tilde{q})$ 2. **参数$\delta$的决定性作用**： - $\delta$刻画阻尼项与质量项的平衡 - 当$\delta$足够大时，系统呈现抛物型特征 - 临界值$\delta_i = (m+1)^2(n+2\sigma-1)^2$决定对数因子的出现 3. **Strauss指数与Fujita指数的竞争**： - 当$m=0$（标准波方程），系统退化到已知结果 - 参数$m$的存在使得临界指数在Strauss型和Fujita型之间插值 ## 相关工作 ### 经典单方程理论 1. **Strauss指数**（波方程$\partial_t^2 u - \Delta u = |u|^p$）： - 临界指数$p_S(n)$是$(n-1)p^2-(n+1)p-2=0$的正根 - Glassey, John, Strauss等人的奠基性工作 2. **Fujita指数**（热方程$\partial_t u - \Delta u = |u|^p$）： - 临界指数$p_F(n) = 1+\frac{2}{n}$ - Fujita (1969)的经典结果 3. **带阻尼波方程**（$\partial_t^2 u - \Delta u + \partial_t u = |u|^p$）： - 临界指数为Fujita指数$p_F(n)$ - Todorova-Yordanov, Zhang等人的工作 ### 耦合系统理论 1. **标准波方程耦合系统**： - Del Santo, Georgiev, Mitidieri等人研究了临界曲线$\Gamma_W(n,p,q)=0$ - 爆破区域比单方程情形扩大 2. **带阻尼的耦合系统**： - Narazaki, Nishihara-Wakasugi等人的工作 - 临界曲线$\Gamma_{DW}(n,p,q)=0$ 3. **带尺度不变阻尼和质量的系统**： - Chen-Palmieri (2019)研究了$m=0$的情形 - 本文推广到一般的$m>-1$ ### 单方程Euler-Poisson-Darboux-Tricomi方程 1. **Nascimento-Palmieri-Reissig**的工作： - 建立了$\delta\geq 0$时的爆破结果 - 证明了$\delta$足够大时的全局存在性 2. **作者之前的工作[17]**： - 建立了$(L^1\cap L^2)-L^2$估计 - 确定了单方程的临界指数$p_F((m+1)n+\frac{\mu-1-\sqrt{\delta}}{2})$ ### 本文的相对优势 1. **首次处理耦合系统**：推广单方程结果到弱耦合系统 2. **新测试函数**：构造适应Gellerstedt算子的测试函数 3. **完整理论**：涵盖爆破和全局存在性的完整刻画 4. **非平凡推广**：证明了耦合效应导致的"偏移"现象 ## 结论与讨论 ### 主要结论 1. **临界曲线的完整刻画**：在阻尼项占优（$\delta_i$足够大）时，临界曲线为 $$\Gamma_m(n,p,q,\beta_1,\beta_2) = 0$$ 2. **爆破与全局存在的二分法**： - $\Gamma_m\geq 0$：解有限时间爆破 - $\Gamma_m<0$：小初值解全局存在 3. **耦合效应的非平凡性**：临界区域不是单方程结果的简单乘积，体现了两个阻尼项和两个质量项的复杂相互作用 ### 局限性 1. **$\delta$的限制**： - 全局存在性结果要求$\delta_i \geq (m+1)^2(n+2\sigma-1)^2$ - 对于较小的$\delta$，问题仍然开放 2. **排除的情形**： - Theorem 2.3不包括$\delta_2 = (m+1)^2(n+2\sigma-1)^2$的情形（Remark 4.4） - 原因是对数因子$(1+\log\tau)^{pθ_1(2p)/2}$在$p\leq\tilde{p}$时导致可积性问题 3. **初值条件**： - 爆破结果需要初值满足特殊的符号条件(2.8) - 全局存在性要求初值小且紧支 4. **技术条件**： - Theorems 2.3-2.4需要额外的技术条件(2.27)和(2.31) - 这些条件的几何意义不够清晰 ### 未来方向 1. **扩展$\delta$的范围**：研究$0<\delta_i < (m+1)^2(n+2\sigma-1)^2$时的情形 2. **临界情形**：研究$\Gamma_m=0$时的精细行为（修正项、生命跨度估计） 3. **一般非线性项**：推广到$F(u,v)$形式的更一般非线性项 4. **数值验证**：虽然理论完备，但数值模拟可以提供直观理解 5. **物理应用**：探索该系统在流体力学、等离子体物理中的应用 ## 深度评价 ### 优点 1. **数学严谨性**： - 证明完整且技术细节清晰 - 通过具体例子验证理论条件的非空性 - 对各种情形进行了细致的分类讨论 2. **方法创新性**： - 测试函数的构造巧妙地适应了Gellerstedt算子 - 加权范数的设计使得能够处理$p\leq\tilde{p}$的情形 - 对对数因子的精细分析 3. **理论完整性**： - 涵盖了爆破和全局存在性两个方面 - 处理了高正则性和低正则性两种情形 - 针对$p,q$与$\tilde{p},\tilde{q}$的不同大小关系给出了完整分类 4. **写作质量**： - 结构清晰，逻辑严密 - 注记(Remarks)提供了丰富的直观解释 - 例子的选择具有代表性 ### 不足 1. **技术复杂性**： - 证明涉及大量技术细节，初学者可能难以理解核心思想 - 八种情形的分类讨论(Cases 1-8)显得繁琐 2. **条件的限制性**： - $\delta_i$的下界限制较强 - 技术条件(2.27)和(2.31)缺乏几何直观 3. **结果的局部性**： - 排除了$\delta_2 = (m+1)^2(n+2\sigma-1)^2$的情形 - 对于$\beta_1=\beta_2$的特殊情形，结果退化为已知结果 4. **应用背景**： - 论文未讨论该系统的物理或工程背景 - 缺少对参数物理意义的解释 ### 影响力 1. **理论贡献**： - 推进了Euler-Poisson-Darboux-Tricomi方程的理论 - 为一般的耦合双曲-抛物混合系统提供了研究范式 - 揭示了Strauss指数和Fujita指数在耦合系统中的竞争机制 2. **方法论价值**： - 测试函数的构造技术可应用于其他变系数方程 - 加权范数的设计思想具有普适性 - 线性估计与非线性迭代的结合方式值得借鉴 3. **可复现性**： - 证明细节充分，可以验证 - 例子具体，易于检验 - 依赖的线性估计已在作者之前工作中建立 4. **后续研究**： - 为临界情形$\Gamma_m=0$的研究铺平道路 - 启发对一般$m$-依赖算子的研究 - 可能推广到更高阶方程或其他几何背景 ### 适用场景 1. **理论研究**： - 偏微分方程理论，特别是双曲-抛物混合型方程 - 临界指数理论 - 爆破理论和全局适定性理论 2. **相关问题**： - 其他带变系数的波动方程 - 带记忆项或分数阶算子的方程 - 高维空间中的非线性波动问题 3. **技术应用**： - 测试函数方法可用于其他爆破问题 - 加权Sobolev空间技术可用于其他临界问题 - 线性估计技术可用于相关算子 ## 参考文献（关键文献） [2] W. Chen, A. Palmieri, *Weakly coupled system of semilinear wave equations with distinct scale-invariant terms in the linear part*, Z. Angew. Math. Phys., 70 (2019) [17] Y. Li, F. Guo, *Global existence for the semilinear regular Euler-Poisson-Darboux-Tricomi equation*, J. Differential Equations, 431 (2025) [22] W.N. Nascimento, A. Palmieri, M. Reissig, *Semi-linear wave models with power non-linearity and scale-invariant time-dependent mass and dissipation*, Math. Nachr., 290 (2017) [28] A. Palmieri, M. Reissig, *A competition between Fujita and Strauss type exponents for blow-up of semi-linear wave equations with scale-invariant damping and mass*, J. Differential Equations, 266 (2019) --- **总体评价**：这是一篇高质量的数学理论论文，在Euler-Poisson-Darboux-Tricomi方程的弱耦合系统理论方面做出了重要贡献。论文技术严谨，方法创新，结果完整。虽然存在一些技术限制，但整体上推进了该领域的理论发展，为后续研究奠定了坚实基础。论文特别值得称道的是对耦合效应非平凡性的深刻揭示，这超越了单方程结果的简单推广。