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

Basic Information

  • Paper ID: 2511.08084
  • Title: Critical curve for weakly coupled system of semilinear Euler-Poisson-Darboux-Tricomi equations
  • Authors: Yuequn Li, Fei Guo
  • Classification: math.AP (Analysis of PDEs)
  • Publication Date: November 11, 2025
  • Paper Link: https://arxiv.org/abs/2511.08084

Abstract

This paper investigates weakly coupled systems of semilinear Euler-Poisson-Darboux-Tricomi equations (EPDTS) with power-type nonlinearities. In the regime where damping terms dominate mass terms, the authors establish a critical curve Γm(n,p,q,β1,β2)=0\Gamma_m(n,p,q,\beta_1,\beta_2)=0 on the p−qp-q plane that distinguishes between global existence and finite-time blow-up. By constructing novel test functions, the paper addresses blow-up phenomena when Γm(n,p,q,β1,β2)≥0\Gamma_m(n,p,q,\beta_1,\beta_2)\geq 0. Based on (L1∩L2)−L2(L^1\cap L^2)-L^2 estimates for the corresponding linear equations established in the authors' prior work, global existence of small-data solutions is derived when Γm(n,p,q,β1,β2)<0\Gamma_m(n,p,q,\beta_1,\beta_2)<0 and damping dominates.

Research Background and Motivation

Research Problem

This paper studies the following weakly coupled semilinear Euler-Poisson-Darboux-Tricomi system:

\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}$$ where $m>-1$, $n\geq 1$ is the spatial dimension, $\mu_i, \nu_i^2\geq 0$ are non-negative parameters, and $p,q>1$. ### Significance of the Research 1. **Theoretical Significance**: This system generalizes classical wave equation theory, involving the interplay between Strauss and Fujita exponents. 2. **Mathematical Structure**: The presence of the Gellerstedt operator $\partial_t^2 - t^{2m}\Delta$ ensures that the damping term $\frac{\mu}{t}\partial_t u$ and mass term $\frac{\nu^2}{t^2}u$ share the same scale-invariance property. 3. **Critical Phenomena**: Determining the critical curve that distinguishes global existence from blow-up is a central problem in PDE theory. ### Limitations of Existing Methods 1. **Non-trivial Extension**: Results for single equations cannot be directly extended to weakly coupled systems; the critical curve is not simply a rectangular region. 2. **Test Function Construction**: The difference between the Gellerstedt operator and the standard wave operator requires constructing new test functions that depend on the parameter $m$. 3. **Role of Parameter $\delta$**: The parameter $\delta_i = (\mu_i-1)^2 - 4\nu_i^2$ characterizes the interplay between damping and mass terms, with its magnitude determining whether the system exhibits parabolic or hyperbolic characteristics. ### Research Motivation 1. Extend the authors' previous work on single equations [17] to coupled systems. 2. Reveal the complex interplay between parameters $\mu_i, \nu_i, m$. 3. Establish a complete theory of global existence and blow-up. ## Core Contributions 1. **Establishment of the Critical Curve**: When $\delta_i$ is sufficiently large, the critical curve is given by $$\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$$ where $\beta_i = \frac{\mu_i+1-\sqrt{\delta_i}}{2}$. 2. **Construction of Novel Test Functions**: New test functions dependent on parameter $m$ are designed specifically for the Gellerstedt operator, establishing blow-up results when $\Gamma_m\geq 0$. 3. **Global Existence Theory**: Based on $(L^1\cap L^2)-L^2$ estimates, global existence of small-data solutions is established when $\Gamma_m<0$, with two cases of regularity (high regularity $\sigma\geq 1$ and low regularity $0<\sigma<1$), each subdivided into three subcases. 4. **Non-trivial Extension**: The results for coupled systems are shown to be non-trivial extensions of single-equation results, with "shift" phenomena in the critical region reflecting the complex balance between two damping and two mass terms. ## Detailed Methodology ### Problem Formulation The Cauchy problem for system (1.1) is studied to determine under what conditions on parameters $(p,q)$ solutions exist globally or blow up in finite time. **Energy Solution Definition**: $(u,v)$ is an energy solution to system (1.1) on $[1,T)$ if: - Regularity: $(u,v)\in (C([1,T);H^1)\cap C^1([1,T);L^2))^2$ - Support property: $\text{supp}(u,v)(t,\cdot)\subset B_{\varphi_m(t)-\varphi_m(1)+M}$, where $\varphi_m(t) = \frac{t^{m+1}}{m+1}$ - Satisfies integral identities (2.6) and (2.7) ### Blow-up Theory Methodology #### Test Function Construction For blow-up problems, the key innovation is constructing new test functions: $$\psi_{d,R}(t,x) = \lambda\left(\frac{t-1}{d}\right)\varphi\left(\frac{x}{R^{m+1}}\right)$$ where $\lambda(t), \varphi(x)\in C_0^\infty$ satisfy special properties (3.3). **Key Technical Points**: 1. Spatial truncation depends on $R^{m+1}$ rather than $R$, which is crucial for adapting to the Gellerstedt operator. 2. Utilization of ODE properties satisfied by $h_i(t) = t^{\beta_i}$ (equation 3.1). 3. Reformulation of the equation in divergence form (3.2). #### Blow-up Proof Strategy 1. Define integral functionals $I_{d,R}, J_{d,R}$. 2. Estimate each term using Hölder's inequality and support properties. 3. Employ condition (2.9) to ensure integrability of integral exponents. 4. Derive contradiction: if $T=\infty$, then under condition $\Gamma_m\geq 0$, we must have $u\equiv 0$ or $v\equiv 0$, contradicting initial conditions (2.8). **Key Inequalities**: $$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}}$$ ### Global Existence Methodology #### Linear Theory Foundation The authors utilize $(L^1\cap L^2)-L^2$ estimates from their prior work [17] for the linear equation $$\partial_t^2 u - t^{2m}\Delta u + \frac{\mu}{t}\partial_t u + \frac{\nu^2}{t^2}u = 0$$ (Propositions 4.1-4.2). **Core Estimates**: For $\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}$$ #### Nonlinear Theory: Fixed Point Method Via the Duhamel principle, the solution is represented as: $$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$$ **Function Space Design** (High Regularity $\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}\}$$ equipped with norm $$\|(u,v)\|_{X(T)} = \sup_{t\in[1,T]}(t^{-\alpha_1}M_1(t,u) + t^{-\alpha_2}M_2(t,v))$$ where $M_i$ contains $L^2, \dot{H}^\sigma$ norms and their time derivatives, with weights $\alpha_i$ determined by conditions $p\lessgtr \tilde{p}, q\lessgtr \tilde{q}$. **Critical Parameters**: $$\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}$$ #### Three Cases Treatment 1. **$p>\tilde{p}, q>\tilde{q}$** (Theorem 2.2): $\alpha_1=\alpha_2=0$, standard contraction mapping. 2. **$p\leq\tilde{p}, q>\tilde{q}$** (Theorem 2.3): $\alpha_1>0, \alpha_2=0$, requires additional technical condition (2.27). 3. **$p>\tilde{p}, q\leq\tilde{q}$** (Theorem 2.4): Dual to Case 2. ### Technical Innovations 1. **$m$-Dependent Test Functions**: Spatial truncation $R^{m+1}$ precisely captures the propagation speed of the Gellerstedt operator. 2. **Refined Weighted Norm Design**: Introduction of weights $t^{-\alpha_i}$ enables establishment of global existence even when $p\leq\tilde{p}$. 3. **Logarithmic Factor Handling**: When $\delta_i = (m+1)^2(n+2\sigma-1)^2$, estimates contain $(1+\log t)^{1/2}$ factors, requiring special integrability analysis. 4. **Stratification of Initial Regularity**: - $\sigma\geq 1$: Control of $\partial_t u$ norms required. - $0<\sigma<1$: Time derivative estimates unnecessary, simplifying the function space. ## Experimental Setup **Note**: This is pure mathematical theoretical research with no numerical experiments. All results are obtained through rigorous mathematical proofs. ### Theoretical Verification The paper validates the non-emptiness of theoretical conditions through concrete examples: **Example 2.1** (Blow-up): $m=0.1, n=2, \mu_1=3, \nu_1=0.25, \mu_2=0.5, \nu_2=0.125$ - Computed: $\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}$ - Choice: $p=2.1>1.172, q=2.2>0.687$ - Result: $\Gamma_m\approx 0.017>0$, solution blows up. **Example 2.3** (Global Existence): $\sigma=1.1, m=0.4, n=3, \mu_1=15, \nu_1=6, \mu_2=8.7, \nu_2=0.25$ - Verification: $\delta_1=52, \delta_2=59.04 > 34.5744 = (m+1)^2(n+2\sigma-1)^2$ - For any $2<p\leq 3.5, 2<q\leq 3.5$, the system has a unique global solution. ## Experimental Results ### Main Theoretical Results #### Blow-up Result (Theorem 2.1) **Conditions**: - $\mu_i>1$ (when $-1<m<0$) or $\mu_i>0$ (when $m\geq 0$) - $\delta_1,\delta_2>0$ - Initial data satisfying (2.8) - $p,q$ satisfying (2.9) and $\Gamma_m(n,p,q,\beta_1,\beta_2)\geq 0$ **Conclusion**: Energy solutions blow up in finite time. **Significance**: 1. The blow-up region is not a rectangle $\{1<p\leq p_F^1, 1<q\leq p_F^2\}$. 2. There exist cases where $p$ or $q$ exceeds the corresponding Fujita exponent yet blow-up still occurs. 3. When $\beta_1\neq\beta_2$, the blow-up region exhibits a "shift" phenomenon. #### Global Existence Results **High Regularity ($\sigma\geq 1$)**: 1. **Theorem 2.2** ($p>\tilde{p}, q>\tilde{q}$): - Requires $\delta_1,\delta_2 \geq (m+1)^2(n+2\sigma-1)^2$ - Solutions satisfy decay estimates: $$\|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}$): - Additional requirement: $\frac{q+1}{pq-1} < \frac{(m+1)n+\beta_2-1}{2}$ - Estimates for $u$ contain additional growth factor $t^{\alpha_1}$; $v$ maintains standard decay. 3. **Theorem 2.4**: Dual to Theorem 2.3. **Low Regularity ($0<\sigma<1$)**: Theorems 2.5-2.7 correspond to Theorems 2.2-2.4 but do not require control of time derivatives. ### Key Findings 1. **Non-triviality of the Critical Curve** (Remark 2.9): - $\Gamma_m\geq 0 \Leftrightarrow p\leq\tilde{p} \text{ and } q\leq\tilde{q}$ - $\Gamma_m<0$ can be realized by three cases: $(p>\tilde{p}, q>\tilde{q})$ or $(p\leq\tilde{p}, q>\tilde{q})$ or $(p>\tilde{p}, q\leq\tilde{q})$ 2. **Decisive Role of Parameter $\delta$**: - $\delta$ characterizes the balance between damping and mass terms. - When $\delta$ is sufficiently large, the system exhibits parabolic characteristics. - The critical value $\delta_i = (m+1)^2(n+2\sigma-1)^2$ determines the appearance of logarithmic factors. 3. **Competition Between Strauss and Fujita Exponents**: - When $m=0$ (standard wave equation), the system reduces to known results. - The presence of parameter $m$ causes critical exponents to interpolate between Strauss and Fujita types. ## Related Work ### Classical Single-Equation Theory 1. **Strauss Exponent** (Wave equation $\partial_t^2 u - \Delta u = |u|^p$): - Critical exponent $p_S(n)$ is the positive root of $(n-1)p^2-(n+1)p-2=0$ - Foundational work by Glassey, John, Strauss, and others. 2. **Fujita Exponent** (Heat equation $\partial_t u - \Delta u = |u|^p$): - Critical exponent $p_F(n) = 1+\frac{2}{n}$ - Classical result by Fujita (1969). 3. **Damped Wave Equations** ($\partial_t^2 u - \Delta u + \partial_t u = |u|^p$): - Critical exponent is the Fujita exponent $p_F(n)$ - Work by Todorova-Yordanov, Zhang, and others. ### Coupled Systems Theory 1. **Standard Wave Equation Coupled Systems**: - Del Santo, Georgiev, Mitidieri, and others studied the critical curve $\Gamma_W(n,p,q)=0$ - The blow-up region is larger than for single equations. 2. **Damped Coupled Systems**: - Work by Narazaki, Nishihara-Wakasugi, and others. - Critical curve $\Gamma_{DW}(n,p,q)=0$. 3. **Systems with Scale-Invariant Damping and Mass**: - Chen-Palmieri (2019) studied the case $m=0$ - This paper extends to general $m>-1$. ### Single-Equation Euler-Poisson-Darboux-Tricomi Equations 1. **Work by Nascimento-Palmieri-Reissig**: - Established blow-up results for $\delta\geq 0$ - Proved global existence when $\delta$ is sufficiently large. 2. **Authors' Prior Work [17]**: - Established $(L^1\cap L^2)-L^2$ estimates - Determined the critical exponent $p_F((m+1)n+\frac{\mu-1-\sqrt{\delta}}{2})$ for single equations. ### Relative Advantages of This Paper 1. **First Treatment of Coupled Systems**: Extends single-equation results to weakly coupled systems. 2. **Novel Test Functions**: Constructs test functions adapted to the Gellerstedt operator. 3. **Complete Theory**: Provides complete characterization of both blow-up and global existence. 4. **Non-trivial Extension**: Demonstrates "shift" phenomena in the critical region due to coupling effects. ## Conclusions and Discussion ### Main Conclusions 1. **Complete Characterization of the Critical Curve**: When damping dominates (sufficiently large $\delta_i$), the critical curve is $$\Gamma_m(n,p,q,\beta_1,\beta_2) = 0$$ 2. **Dichotomy Between Blow-up and Global Existence**: - $\Gamma_m\geq 0$: Solutions blow up in finite time. - $\Gamma_m<0$: Small-data solutions exist globally. 3. **Non-triviality of Coupling Effects**: The critical region is not a simple product of single-equation results, reflecting the complex interplay between two damping and two mass terms. ### Limitations 1. **Restriction on $\delta$**: - Global existence results require $\delta_i \geq (m+1)^2(n+2\sigma-1)^2$ - The problem remains open for smaller $\delta$. 2. **Excluded Cases**: - Theorem 2.3 excludes the case $\delta_2 = (m+1)^2(n+2\sigma-1)^2$ (Remark 4.4) - Reason: logarithmic factor $(1+\log\tau)^{pθ_1(2p)/2}$ causes integrability issues when $p\leq\tilde{p}$. 3. **Initial Data Conditions**: - Blow-up results require initial data satisfying special sign conditions (2.8) - Global existence requires small and compactly supported initial data. 4. **Technical Conditions**: - Theorems 2.3-2.4 require additional technical conditions (2.27) and (2.31) - The geometric meaning of these conditions is not entirely clear. ### Future Directions 1. **Extension of $\delta$ Range**: Study the case $0<\delta_i < (m+1)^2(n+2\sigma-1)^2$. 2. **Critical Case**: Investigate fine behavior when $\Gamma_m=0$ (correction terms, lifespan estimates). 3. **General Nonlinearities**: Extend to more general nonlinear terms of the form $F(u,v)$. 4. **Numerical Verification**: Although the theory is complete, numerical simulations could provide intuitive understanding. 5. **Physical Applications**: Explore applications in fluid mechanics and plasma physics. ## In-Depth Evaluation ### Strengths 1. **Mathematical Rigor**: - Complete proofs with clear technical details - Concrete examples verify non-emptiness of theoretical conditions - Careful case-by-case analysis for various scenarios. 2. **Methodological Innovation**: - Test function construction cleverly adapts to the Gellerstedt operator - Weighted norm design enables handling of $p\leq\tilde{p}$ cases - Refined analysis of logarithmic factors. 3. **Theoretical Completeness**: - Covers both blow-up and global existence - Addresses both high and low regularity cases - Complete classification based on relationships between $p,q$ and $\tilde{p},\tilde{q}$. 4. **Writing Quality**: - Clear structure with rigorous logic - Remarks provide rich intuitive explanations - Representative choice of examples. ### Weaknesses 1. **Technical Complexity**: - Proofs involve extensive technical details, potentially difficult for newcomers - Eight-case classification (Cases 1-8) appears cumbersome. 2. **Restrictive Conditions**: - Strong lower bound on $\delta_i$ - Technical conditions (2.27) and (2.31) lack geometric intuition. 3. **Local Nature of Results**: - Excludes the case $\delta_2 = (m+1)^2(n+2\sigma-1)^2$ - Results degenerate to known results when $\beta_1=\beta_2$. 4. **Application Background**: - Paper does not discuss physical or engineering background - Lacks explanation of physical meaning of parameters. ### Impact 1. **Theoretical Contribution**: - Advances theory of Euler-Poisson-Darboux-Tricomi equations - Provides research paradigm for general coupled hyperbolic-parabolic mixed systems - Reveals competition mechanism between Strauss and Fujita exponents in coupled systems. 2. **Methodological Value**: - Test function construction techniques applicable to other variable-coefficient equations - Weighted norm design ideas have universal applicability - Combination of linear estimates and nonlinear iteration worthy of emulation. 3. **Reproducibility**: - Sufficient proof details for verification - Concrete examples easy to check - Dependent linear estimates established in authors' prior work. 4. **Subsequent Research**: - Paves the way for studying critical case $\Gamma_m=0$ - Inspires research on general $m$-dependent operators - May extend to higher-order equations or other geometric settings. ### Applicable Scenarios 1. **Theoretical Research**: - PDE theory, especially hyperbolic-parabolic mixed-type equations - Critical exponent theory - Blow-up theory and global well-posedness theory. 2. **Related Problems**: - Other variable-coefficient wave equations - Equations with memory terms or fractional operators - Nonlinear wave problems in higher dimensions. 3. **Technical Applications**: - Test function method applicable to other blow-up problems - Weighted Sobolev space techniques for other critical problems - Linear estimate techniques for related operators. ## Key References [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) --- **Overall Assessment**: This is a high-quality mathematical theory paper making important contributions to the theory of weakly coupled Euler-Poisson-Darboux-Tricomi systems. The paper is technically rigorous, methodologically innovative, and theoretically complete. While certain technical limitations exist, the work overall advances the field and provides a solid foundation for subsequent research. Particularly noteworthy is the deep revelation of the non-triviality of coupling effects, which transcends simple extensions of single-equation results.