Gelfand-type problem for turbulent jets: sharp $L^\infty$ bound on extremal solutions

技术摘要：湍流射流中 Gelfand 型问题极值解的锐利 $L^\infty$ 界限

问题陈述
本文研究了由 [GHH] 推导出的、用于模拟反应性湍流射流中热爆炸的边界值问题。该问题在单位圆盘 $B \subset \mathbb{R}^2$ 上构建，属于 Gelfand 型问题范畴。其控制方程为：
$$\begin{cases} -\Delta u - \alpha r \phi(r) \frac{\partial u}{\partial r} = \lambda \psi(r) f(u) & \text{在 } B \text{ 内}, \\ u > 0 & \text{在 } B \text{ 内}, \\ u = 0 & \text{在 } \partial B \text{ 上}, \end{cases}$$
其中 $u$ 代表归一化温度，$\lambda$ 是 Frank-Kamenetskii 参数，$\alpha$ 是注入速度。函数 $\phi$ 和 $\psi$ 分别描述流速和反应组分的分布轮廓，而 $f$ 是正的、递增且凸的反应速率，满足积分条件 $\int_0^\infty ds/f(s) < \infty$。

前人工作 [GMN] 确定了临界参数 $\lambda^*(\alpha)$ 以及极值解 $u^*_\alpha$（即当 $\lambda \nearrow \lambda^*$ 时最小解的极限）。虽然 [GMN] 证明了当 $\alpha \to \infty$ 时 $u^*_\alpha(0) \to \infty$，但他们在对 $f$ 施加额外假设的情况下，仅能建立 $u^*_\alpha(0)$ 的上界。本文的核心目标是通过推导在强流极限（$\alpha \to \infty$）下 $u^*_\alpha(0)$ 的匹配上下界，从而填补这一空白，且不对 $f$ 施加额外的限制条件。

方法论
作者利用极值解的径向对称性，将偏微分方程简化为常微分方程。通过引入权重函数 $\mu(r) = \exp(\alpha \int_0^r s\phi(s)ds)$ 并进行变量代换 $z = \alpha r^2$，将问题转化为关于缩放解 $\bar{u}(z)$ 的积分方程：
$$\bar{u}'(z) = -\frac{1}{z\bar{\mu}(z)} \int_0^z \bar{\mu}(s)\bar{\psi}(s)F(\bar{u}(s))ds,$$
其中 $F = \frac{\lambda^*}{4\alpha}f$。

分析的核心依赖于极值解所满足的半稳定性条件（semi-stability condition）。在缩放坐标下，该条件表现为一种加权 Poincaré 型不等式：
$$\int_0^\alpha z (\bar{\eta}')^2 \bar{\mu} dz \geq \int_0^\alpha \bar{\eta}^2 \bar{\psi} \bar{\mu} F'(\bar{u}) dz,$$
对于所有 $\eta \in H^1_0(B)$ 成立。

证明策略包含两个组成部分：

1. 上界（定理 1.1）： 作者构造了一个特定的测试函数 $\bar{\eta}_z$（分段常数/对数型）来利用半稳定性不等式。这使得他们能够通过一个对数项来限制反应项导数的积分，即 $\int F'(\bar{u})$。结合 ODE 的积分形式以及 $f$ 的凸性，他们推导出了一个关于 $f'(u(0))$ 与 $\log \alpha$ 关系的不等式。
2. 下界（定理 1.2）： 作者采用了反证法。假设下界失效，这意味着存在一个扰动 $v > 0$ 使得 $\bar{u} + v$ 仍是该方程的解。通过分析线性算子并利用基于 $\mu$ 和 $\psi$ 特性的核估计（引理 4.2），他们表明对于足够大的 $\alpha$，这种正扰动在不违反 $\lambda^*$ 的极值性前提下是不可能存在的。具体而言，如果该界限失效，则可以构造出一个对应于参数 $\lambda > \lambda^*$ 的解，这与临界参数的定义相矛盾。

主要贡献与结果
本文确立了在 $\alpha \to \infty$ 时，极值解在原点处的尖锐渐近界限 $u^*_\alpha(0)$：

• 定理 1.1（上界）： 存在常数 $\alpha_0, c, C > 0$，使得对于所有 $\alpha \geq \alpha_0$，
  $$u^*_\alpha(0) \leq c A,$$
  其中 $A$ 是满足 $f'(A) = C \log \alpha$ 的唯一解。

• 定理 1.2（下界）： 存在常数 $\tilde{\alpha}_0, \tilde{c}, \tilde{C} > 0$，使得对于所有 $\alpha \geq \tilde{\alpha}_0$，
  $$u^*_\alpha(0) \geq \tilde{c} \tilde{A},$$
  其中 $\tilde{A}$ 是满足 $f'(\tilde{A}) = \tilde{C} \log \alpha$ 的唯一解。

这些结果证实了极值解峰值温度的增长受控于反应速率在流速对数尺度下的导数的反函数。

意义
本文声称解决了 [GMN] 留下的关于湍流射流模型中极值解 $L^\infty$ 行为的部分开放性问题。其主要意义在于移除了此前为了获得上界而必须对非线性反应率施加的“额外假设”。通过证明对于一般的凸反应速率类（满足标准积分条件），匹配的上下界均成立，作者完整地刻画了强流极限下极值解的爆发行为。这一结果巩固了 Gelfand 框架下反应性湍流射流中热爆炸极限的理论理解。