The Newton's problem assuming non-constant density of the fluid

技术摘要：非恒定密度下的牛顿最小阻力问题

问题陈述
本文研究了牛顿经典最小阻力问题的一个变体，该问题旨在寻找以恒定速度在流体中运动且阻力最小的固体形状。虽然经典模型假设流体是均匀的，但本研究针对的是更为现实的场景，即流体密度随高度呈指数衰减，这是地球大气的典型特征。密度建模为 $\rho(z) = \rho_0 e^{-cz}$，其中 $z$ 代表高度。物体沿 $z$ 坐标轴运动，阻力泛函基于正弦平方压力定律推导得出，假设无摩擦且粒子在撞击后不会与表面再次碰撞。

对于由图 $u: \Omega \to \mathbb{R}$ 定义的曲面，所得的阻力泛函为：
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
与经典牛顿泛函不同，指数权重 $e^{-u}$ 的存在给相关的欧拉 - 拉格朗日方程带来了显著的结构变化。本文明确指出，其重点在于经典 $C^2$ 极值（径向解）的存在性与正则性，而非全局最小值的确定，并指出在经典设定下，无约束的最小阻力问题往往缺乏 $C^2$ 正则性。

方法论
分析聚焦于径向解 $u(r)$，其中曲面是绕 $z$ 轴旋转而成的。泛函 $E[u]$ 的欧拉 - 拉格朗日方程简化为一个二阶非线性常微分方程（ODE）：
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
本研究处理的是初始值问题，条件为 $u(0) = 0$ 和 $u'(0) = 0$。主要挑战在于方程在 $r=0$ 处的奇异性，这阻碍了经典 ODE 存在定理的直接应用。

为克服这一困难，作者采用了基于巴拿赫不动点定理的不动点论证。该 ODE 被重构为一个积分方程，涉及特定单调函数 $f(x) = x/(1+x^2)^2$ 的逆。通过在巴拿赫空间 $C^1([0, R])$ 的闭球上定义算子 $T$，作者证明了对于足够小的 $R$，$T$ 是一个压缩映射，从而确立了 $C^2$ 解的局部存在性。

为了分析解的全局行为及其最大定义域，问题被转化为相平面分析。曲线由弧长 $s$ 重新参数化，引入角度 $\theta(s)$，使得 $u'(r) = \tan \theta$。这将二阶 ODE 转化为关于变量 $x(s) = r$ 和 $\theta(s)$ 的一阶方程自治系统。通过考察相平面 $\Theta$ 中的轨迹，特别是平衡点及临界线附近的流，研究了解的行为。

主要贡献与结果

1. 局部存在性与正则性：本文证明（定理 2.1），初始值问题 $u(0)=u'(0)=0$ 在某个 $R > 0$ 下 admits 一个解 $u \in C^2([0, R])$。一个关键发现是，解与旋转轴正交相交（$u'(0)=0$）。这与经典恒定密度模型形成鲜明对比，在经典模型中，径向解除非是平凡的常数函数，否则无法正交到达旋转轴。证明进一步确立了 $u''(0) = 1/4$。

2. 单调性与凸性：分析表明，斜率 $u'(r)$ 在解的最大定义域内严格为正且递增。因此，径向解 $u(r)$ 是凸的。

3. 有限最大定义域与临界斜率：定理 2.2 确立了解的最大定义域是一个有限区间 $[0, r_M)$。解终止于 ODE 中 $u''$ 的系数消失的点，对应于一个临界斜率：
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   数值计算表明 $r_M \approx 1.230$ 且 $u(r_M) \approx 0.228$。相平面分析显示，轨迹渐近地趋近于直线 $\theta = \pi/6$（其中 $\tan \theta = 1/\sqrt{3}$），且定义域在有限半径处结束。

4. 相平面结构：该研究刻画了相关自治系统的相平面，识别了平衡点，并表明特定点周围的流由拓扑圆组成。研究表明，从原点出发的轨迹终止于水平线 $\theta = \pi/6$，且 $x$ 坐标为有限值。

意义与主张
本文声称，流体密度的指数衰减从根本上改变了最小阻力体的几何形状，与经典均匀模型相比尤为显著。具体而言，存在与对称轴正交相交的光滑（$C^2$）径向解是一个新颖的结果。在经典模型中，对于非平凡解，轴处的这种光滑性是不可能的，通常需要“平头”圆盘（导数的不连续性）来最小化阻力。

作者提出，这些光滑的径向解为分层介质中的最优形状提供了理论基础。斜率在定义域边缘被 $1/\sqrt{3}$ 限制的几何约束，被强调为设计穿越大气层物体时的重要特征。本文结论指出，虽然这些 $C^2$ 极值存在且具有正则性，但它们是否代表全局最小值而非经典非光滑最小值，仍是一个有待未来研究解决的开放问题。该工作被呈现为对这些特定极值的存在性和性质的严谨数学调查，并未提出新的实验应用或声称立即的工程实施。