🔢 mathematics

Blow-up of solutions to the Euler-Poisson-Darbox equation with critical power nonlinearity

This paper establishes the finite-time blow-up and provides an upper bound for the lifespan of solutions to the singular Cauchy problem of the semilinear Euler-Poisson-Darboux equation with critical power nonlinearity by employing an improved test function to derive enhanced lower bounds, thereby partially resolving an open problem posed by D'Abbicco.

Original authors: Mengting Fan, Ning-An Lai, Hiroyuki Takamura

Published 2026-03-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mengting Fan, Ning-An Lai, Hiroyuki Takamura

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The "Ticking Bomb" Problem

Imagine you have a giant, invisible balloon floating in space. This balloon represents a physical system (like a sound wave or a gravitational field) that changes over time.

The Euler-Poisson-Darboux (EPD) equation is the rulebook for how this balloon behaves. It has two main forces acting on it:

The Push (Nonlinearity): The balloon wants to expand and get bigger on its own. The faster it grows, the harder it pushes itself to grow even faster. This is the $|u|^p$ part of the equation.
The Brake (Damping): There is a "friction" or "air resistance" that tries to slow the balloon down. In this specific problem, the brake is a bit weird: it gets weaker as time goes on (represented by $\frac{\mu}{t}$ ). At the very beginning ( $t=0$ ), the brake is infinitely strong, but as time passes, it fades away.

The Goal: Mathematicians want to know: Will this balloon pop (blow up) in finite time, or will it settle down and exist forever?

The "Critical" Moment

In the world of these equations, there is a "Goldilocks" zone called the Critical Power.

Too Weak (Subcritical): If the balloon's self-pushing force is weak, the brake wins. The balloon survives forever.
Too Strong (Supercritical): If the self-pushing force is too strong, the balloon pops almost instantly.
Just Right (Critical): This is the dangerous middle ground. The forces are perfectly balanced. It's like balancing a pencil on its tip. In this state, the balloon might survive, or it might pop, but it's incredibly hard to tell.

What this paper does: The authors prove that for this specific "Critical" case, the balloon will eventually pop. They also calculate exactly how long it takes to pop (the "lifespan").

The Analogy: The Leaky Bucket and the Logarithmic Trick

To prove the balloon pops, the authors had to use a clever trick because standard methods failed at this "Critical" balance point.

1. The Old Way (Subcritical):
Previously, the authors (and others) showed that if the balloon pushes too hard, it pops quickly. They used a "test function" (a mathematical probe) to measure the energy. It was like shining a flashlight on the balloon; if the light got too bright, they knew it was about to explode.

2. The New Problem (Critical):
At the critical point, the flashlight isn't bright enough. The balloon is so stable that the standard probe can't detect the explosion coming. The math says the lifespan is "exponential," meaning it could take a very long time to pop, but it will pop.

3. The Authors' Solution (The Improved Probe):
The authors invented a new, super-sensitive probe.

The "Rescaled" Lens: Instead of just looking at the balloon once, they looked at it through a lens that zoomed in and out at different speeds (using "rescaled functions"). This allowed them to see tiny details of the growth that were previously invisible.
The "Logarithmic" Growth: This is the most important part. In the old subcritical cases, the energy grew like a straight line. In this critical case, the energy grows very slowly, like a logarithm (think of how a password strength meter slowly ticks up).
- The Metaphor: Imagine the balloon is filling with water. In the "easy" cases, the water level rises fast. In this "critical" case, the water rises so slowly you might think it will never overflow. But the authors proved that even though it rises slowly, it does rise. They found an "extra" bit of growth (the logarithmic term) that guarantees the water will eventually reach the top and spill over.

The "Singular" Hurdle

There was one more snag. The "brake" in the equation ( $\frac{\mu}{t}$ ) is broken at the very start ( $t=0$ ). It's like trying to drive a car where the brakes are welded shut at the moment you start the engine.

The Fix: The authors used a special mathematical tool called a Modified Bessel Function. Think of this as a "magic shield" that wraps around the start of the timeline. It smooths out the broken brake so the math doesn't crash, allowing them to calculate the future behavior without getting stuck at the beginning.

The Result: How Long Until It Pops?

The paper concludes with a formula for the Lifespan ( $T$ ).

If you start with a very small amount of energy (a tiny initial push, represented by $\epsilon$ ), the time it takes to blow up is roughly:
$T \approx e^{1/\epsilon^{\text{something}}}$

In plain English:
If you make the initial push tiny, the balloon will survive for an incredibly long time (exponentially long). However, no matter how small you make that push, it will eventually explode. It cannot exist forever.

Why Does This Matter?

This paper solves a specific "Open Problem" (Problem 1) left by a mathematician named D'Abbicco.

Before: We knew the balloon pops if pushed hard, and survives if pushed soft. We didn't know for sure what happened in the "perfect balance" zone.
Now: We know that even in the perfect balance, the "push" eventually wins. The system is inherently unstable in the long run.

This helps physicists and mathematicians understand the limits of stability in systems involving waves, gravity, and sound, ensuring we know when a system is destined to collapse, even if it looks stable for a long time.

1. Problem Statement

The paper investigates the singular Cauchy problem for the semilinear Euler-Poisson-Darboux (EPD) equation in $\mathbb{R}^n$ ( $n \ge 1$ ):

$\begin{cases} u_{tt} - \Delta u + \frac{\mu}{t} u_t = t^\alpha |u|^p, & t > 0, \, x \in \mathbb{R}^n, \\ u(0, x) = u_0(x), \\ u_t(0, x) = 0, & x \in \mathbb{R}^n, \end{cases}$

Key Parameters:

$\mu > 0$ : Damping coefficient (singular at $t=0$ ).
$\alpha \ge 0$ : Power of the time-weight in the nonlinear source term.
$p > 1$ : Nonlinear exponent.
Critical Case: The authors focus on the critical exponent $p = p_S(n + \mu, \alpha)$ , defined as the positive root of the quadratic equation:
$\gamma(n, \mu, \alpha, p) := 2 + (n + \mu + 1 + 2\alpha)p - (n + \mu - 1)p^2 = 0.$

Context:
Previous work by the authors established finite-time blow-up for the subcritical case ( $p < p_S$ ). This paper addresses the critical case, which is significantly more delicate because it represents the borderline between finite-time blow-up and global existence. In the critical case, the lifespan of solutions is expected to be of exponential type (much longer than the polynomial lifespan in the subcritical case), making the proof of blow-up technically challenging.

2. Methodology

The authors employ the test function method combined with an iteration argument, but with specific innovations tailored to the critical case and the singular damping term.

A. Construction of a Special Test Function

To handle the critical exponent, the authors construct a refined test function $b_q(t, x)$ by integrating scaled solutions of a related ODE with respect to a scale parameter $\lambda$ :

Base Functions: They utilize the modified Bessel function of the second kind, $K_{\frac{\mu-1}{2}}(t)$ , to solve the ODE $h'' - (\frac{\mu}{t}h)' - h = 0$ . This choice is crucial for handling the singularity at $t=0$ .
Scaling: They define $h_\lambda(t) = h(\lambda t)$ and a spatial function $\phi_\lambda(x)$ related to the fundamental solution of the wave operator.
Integration: The test function is defined as:
$b_q(t, x) = \int_0^1 h_\lambda(t) \phi_\lambda(x) \lambda^{q-1} d\lambda,$
where $q$ is chosen specifically to match the critical exponent $p$ .
Properties: This construction yields a function satisfying the homogeneous adjoint equation and possessing specific asymptotic behaviors ( $b_q \sim t^{-q}$ and $\partial_t b_q \sim t^{\mu/2 - (n+1)/2}(t+2-r)^{1/p - 1}$ ) that are superior to standard test functions for critical problems.

B. Handling the Singularity at $t=0$

A major difficulty in singular EPD equations is the term $\frac{\mu}{t}u_t$ when integrating by parts near $t=0$ .

The authors use the asymptotic expansion and recurrence relations of the modified Bessel function $K_{\frac{\mu-1}{2}}(t)$ to show that the singular terms cancel out or remain bounded as $t \to 0$ .
Specifically, they prove that $\lim_{t \to 0^+} (-\partial_t b_q + \frac{\mu}{t}b_q) \ge 0$ , ensuring the boundary term at $t=0$ is non-negative and contributes to the blow-up mechanism.

C. Enhanced Lower Bound via Logarithmic Growth

The core innovation for the critical case is obtaining an additional logarithmic growth in the lower bound of the functional.

They define a functional $Y(M)$ involving the spacetime integral of the nonlinear term weighted by the test function.
By carefully estimating the derivative $Y'(M)$ and utilizing the specific asymptotic properties of $b_q$ , they derive:
$Y'(M) \ge C \varepsilon^p M^{-1}.$
Integrating this yields a logarithmic lower bound: $Y(M) \ge C \varepsilon^p \ln M$ .
This logarithmic factor is essential for closing the iteration argument in the critical case, distinguishing it from the polynomial bounds used in subcritical cases.

3. Key Contributions

Resolution of an Open Problem: The paper provides a partial answer to Open Problem 1 posed by D'Abbicco (2021) regarding the blow-up behavior of the semilinear EPD equation at the critical exponent.
Critical Blow-up Result: It establishes the finite-time blow-up for the critical exponent $p = p_S(n + \mu, \alpha)$ under specific conditions on $\mu$ and $\alpha$ .
Lifespan Estimate: The authors derive an upper bound for the lifespan $T(\varepsilon)$ of the solution:
$T(\varepsilon) \le \exp\left( C \varepsilon^{-p(p-1)} \right).$
This confirms the exponential nature of the lifespan in the critical regime.
Technical Innovation: The introduction of the integrated test function $b_q$ and the rigorous handling of the $t=0$ singularity using modified Bessel functions provide a robust framework for analyzing critical singular wave equations.

4. Main Results (Theorem 1.2)

Under the assumptions:

$n \ge 3$ (or $n=2$ with $\alpha > 0$ ),
$0 < \mu \le \mu^*(n, \alpha)$ (where $\mu^*$ is a specific threshold ensuring $p_S \ge p_F$ ),
$p = p_S(n + \mu, \alpha)$ ,
Initial data $u_0 \in H^1(\mathbb{R}^n)$ is positive with compact support.

Conclusion: There exists a small $\varepsilon_0 > 0$ such that for any $0 < \varepsilon \le \varepsilon_0$ , the energy solution to the Cauchy problem blows up in finite time, and the lifespan satisfies the exponential upper bound mentioned above.

5. Significance

Theoretical Advancement: This work bridges the gap between subcritical and supercritical regimes for singular damped wave equations. It confirms that the "critical" exponent defined by the Strauss-type equation indeed acts as a threshold for blow-up, even in the presence of singular damping and time-weighted nonlinearity.
Methodological Impact: The technique of using scale-integrated test functions based on Bessel functions offers a new tool for researchers dealing with critical exponents in hyperbolic PDEs with singular coefficients.
Physical Relevance: The Euler-Poisson-Darboux equation models phenomena in acoustics, surface theory, and gravitational fields. Understanding the blow-up behavior at critical thresholds is vital for predicting the stability and lifespan of solutions in these physical contexts.

In summary, this paper successfully extends the blow-up theory to the critical case for the singular EPD equation by overcoming the technical hurdles of the singularity at $t=0$ and the delicate nature of the critical exponent through a novel test function construction.