[Paper Review] Two-Dimensional Problems of Minimal Resistance in a Medium of Positive Temperature

A. Yu. Plakhov, Delfim F. M. Torres|ArXiv.org|Apr 9, 2004

Thermoelastic and Magnetoelastic Phenomena11 references3 citations

TL;DR

This paper extends Newton's classical minimal resistance problem to two-dimensional bodies moving through a rarefied medium with particles in thermal motion, modeled by a Gaussian velocity distribution. It provides analytical formulas for resistance and identifies four distinct minimizer types, with numerical results showing that resistance scales as $\tilde{R}(V,h) \to 1 - h/2$ for high velocities and diverges as $V \to 0$, revealing the critical role of temperature in shaping optimal body geometries.

ABSTRACT

We study the Newton-like problem of minimal resistance for a two-dimensional body moving with constant velocity in a homogeneous rarefied medium of moving particles. The distribution of the particles over velocities is centrally symmetric. The problem is solved analytically; the minimizers are shown to be of four different types. Numerical results are obtained for the physically significant case of gaussian circular distribution of velocities, which corresponds to a homogeneous ideal gas of positive temperature.

Motivation & Objective

• To extend Newton's minimal resistance problem to a medium with particles in thermal motion, modeling a gas at positive temperature.
• To derive analytical expressions for resistance in the two-dimensional case under a general spherically symmetric velocity distribution.
• To analyze the structure of minimizers and classify them into four distinct types based on geometric and dynamic properties.
• To study the physically relevant case of a Gaussian circular velocity distribution, corresponding to an ideal gas at positive temperature.
• To perform numerical simulations and derive asymptotic resistance behaviors for both low and high-velocity regimes.

Proposed method

• Model the medium as a flux of particles with a spherically symmetric velocity distribution centered at $-Ve_d$, where $V > 0$ is the body's speed.
• Use the pressure formula $\pi(n_x) = -\int (v|n_x)_-^2 \rho(v)\, dv \cdot n_x$ to compute force on the body's boundary.
• Reduce the body to a convex, axisymmetric shape defined by functions $f_\pm(|x'|)$, with fixed height $h = -f_+(0) - f_-(0)$.
• Derive the resistance functional $R(\mathcal{B}) = \int_{\partial\mathcal{B}} \pi(n_x)\, d\mathcal{H}^{d-1}(x)$ and specialize to $d=2$.
• For the Gaussian case, express pressure components $p_\pm(u,V)$ via a double integral over velocity angles and magnitudes, then simplify using Bessel functions and error functions.
• Perform numerical simulations using Maple to compute reduced resistance $\tilde{R}(V,h) = R(V,h)/V^2$ and analyze asymptotic limits as $V \to 0$ and $V \to \infty$.

Experimental results

Research questions

• RQ1 How does thermal motion of particles (modeled by a Gaussian velocity distribution) affect the minimal resistance of a 2D body compared to the classical zero-temperature case?
• RQ2 What are the analytical forms of the resistance functional for a 2D convex, axisymmetric body in a medium with thermal particle motion?
• RQ3 How many distinct types of minimizers exist in the 2D problem under positive-temperature conditions, and what characterizes each type?
• RQ4 What is the asymptotic behavior of the minimal resistance as the body's velocity $V \to 0$ and $V \to \infty$?
• RQ5 How does the reduced resistance $\tilde{R}(V,h)$ depend on body height $h$ and velocity $V$ in the Gaussian case?

Key findings

• The reduced resistance $\tilde{R}(V,h)$ tends to $1 - h/2$ as $V \to \infty$ for $h \leq 1$, and to $1/(1 + h^2)$ for $h \geq 1$, indicating a transition in optimal shape at $h=1$.
• As $V \to 0$, $\tilde{R}(V,h)$ diverges, with $\sqrt{2/\pi} \cdot \lim_{V \to 0} (V \tilde{R}(V,h)) = 2 - h/a^5$ for $h \leq 2a$, where $a \approx 1.27$, showing strong dependence on body height at low speeds.
• For small $V$, the leading-order asymptotic behavior of the pressure is $p_\pm(u,V) = \pm \frac{1}{2} + \sqrt{2/\pi} \cdot \frac{V}{\sqrt{1+u^2}} + O(V^2)$, indicating a linear correction to the zero-temperature limit.
• The pressure $p_+(u,V)$ grows as $V^2/(1+u^2)$ for large $V$, matching the Newtonian (zero-temperature) result up to a factor of $V^2$, confirming the recovery of classical behavior at high speeds.
• Numerical simulations confirm four distinct minimizer types, separated by curves in the $V$-$h$ plane, with the lower curve approaching 1 as $V \to \infty$.
• The analytical resistance formulas, while complex, are amenable to numerical evaluation using modern computational tools like Maple, enabling full characterization of the optimal shapes across parameter regimes.