[论文解读] Two-dimensional Newton's Problem of Minimal Resistance

Cristiana J. Silva, Delfim F. M. Torres|ArXiv.org|Jul 7, 2006

Spacecraft Dynamics and Control参考文献 13被引用 6

一句话总结

本文研究了牛顿经典三维最小阻力问题的二维版本，表明与三维情况不同，在高度 H 超过底面半径 r 的 √3/3 倍时，无约束问题在二维情况下是适定的，且最小化器为三角形形状。此外，研究进一步证明，在受限情况下，当 H ≤ r 时，存在无穷多个不同的最小化器，此时最小阻力值为 r − H/2。

ABSTRACT

Newton's problem of minimal resistance is one of the first problems of optimal control: it was proposed, and its solution given, by Isaac Newton in his masterful Principia Mathematica, in 1686. The problem consists of determining, in dimension three, the shape of an axis-symmetric body, with assigned radius and height, which offers minimum resistance when it is moving in a resistant medium. The problem has a very rich history and is well documented in the literature. Of course, at first glance, one suspects that the two dimensional case should be well known. Nevertheless, we have looked into numerous references and ask at least as many experts on the problem, and we have not been able to identify a single source. Solution was always plausible to everyone who thought about the problem, and writing it down was always thought not to be worthwhile. Here we show that this is not the case: the two-dimensional problem is more rich than the classical one, being, in some sense, more interesting. Novelties include: (i) while in the classical three-dimensional problem only the restricted case makes sense (without restriction on the monotonicity of admissible functions the problem doesn't admit a local minimum), we prove that in dimension two the unrestricted problem is also well-posed when the ratio height versus radius of base is greater than a given quantity; (ii) while in three dimensions the (restricted) problem has a unique solution, we show that in the restricted two-dimensional problem the minimizer is not always unique - when the height of the body is less or equal than its base radius, there exists infinitely many minimizing functions.

研究动机与目标

• 研究牛顿经典三维最小阻力问题在二维情形的类比。
• 确定在二维情况下，无单调性约束（即导数无单调性限制）的无约束问题在何种条件下是适定的。
• 分析受限二维情况下最小化器的结构与唯一性。
• 将二维问题与经典的三维版本进行比较，突出其在适定性与解唯一性方面的关键差异。
• 在各种几何约束（H 与 r 的关系）下，对最小化形状与阻力值进行全面表征。

提出的方法

• 将二维阻力泛函表述为 R[y] = ∫₀ʳ x / (1 + y′(x)²) dx，模拟物体在稀薄介质中运动的情形。
• 应用变分法与最优控制理论分析最小化器的必要条件，包括庞特里亚金最大值原理。
• 在无约束与受限（y′(x) ≥ 0）两种设定下，推导阻力泛函的临界点。
• 构造一类分段线性控制 uₙ(x)（形式如式 (21)），在交替区间上具有零斜率与单位斜率，以生成候选最小化器。
• 通过解析优化证明，当 H ≤ r 时，阻力 Rₙ = r − H/2 对所有 n ∈ ℕ 恒定，从而证明其在细化过程中的不变性。
• 通过二阶条件建立局部极小值的存在性，并验证泛函有下界且能取到下确界。

实验结果

研究问题

• RQ1 在何种几何条件下，二维无约束牛顿最小阻力问题是适定的？
• RQ2 在二维情况下，若导数无单调性约束，是否能获得定义良好的解？这与三维情况形成对比。
• RQ3 当 H ≤ r 时，受限二维问题的解结构与 H > r 时有何不同？
• RQ4 在二维受限问题中，是否可能存在多个不同的最小化器？若存在，其数量是多少？
• RQ5 二维情况下能达到的最小阻力值是多少？其如何依赖于 H/r 的比值？

主要发现

• 当 H > √3/3 r 时，二维无约束牛顿最小阻力问题适定，且存在唯一的局部最小化器，其形式为线性函数 y(x) = (H/r)x。
• 当 H > r 时，最小阻力值为 r³ / (r² + H²)，由三角形形状的物体实现。
• 当 H ≤ r 时，受限问题存在无穷多个不同的最小化器，所有最小化器均给出相同的最小阻力值 r − H/2。
• 所有 H ≤ r 情况下的最小化器均为形式 (21) 的分段线性函数，在交替区间上具有零斜率与单位斜率。
• 阻力值 r − H/2 在所有此类最小化器中保持不变，无论区间数 n 如何。
• 经典三维问题要求单调性（y′(x) ≥ 0）以保证适定性，但在二维情况下，对于某些 H/r 比值，该限制并非必要。