arxiv: 2605.11641 · v1 · submitted 2026-05-12 · 🧮 math.AP · math.DG

Recognition: no theorem link

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

Rafael L\'opez

Pith reviewed 2026-05-13 01:12 UTC · model grok-4.3

classification 🧮 math.AP math.DG

keywords Newton's minimal resistanceexponential densityradial solutionslocal existencefixed-point theoremcalculus of variationsoptimal shape

0 comments

The pith

Radial solutions to Newton's minimal resistance problem with exponential fluid density exist locally but terminate at slope 1/sqrt(3).

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends Newton's classic problem of the body shape minimizing resistance while traveling through a fluid, now allowing the fluid density to decrease exponentially with altitude. It establishes that radial solutions u(r) with flat initial conditions at the origin exist and remain regular on a neighborhood of zero by means of a fixed-point argument. These solutions cannot be continued indefinitely: their maximal interval is finite and ends exactly when the slope reaches 1/sqrt(3). A reader cares because the result supplies an explicit mathematical description of how non-uniform density limits the size of the optimal shape in an axisymmetric setting.

Core claim

We prove the local existence and regularity of radial solutions u(r) satisfying the initial conditions u(0)=u'(0)=0 using a fixed-point theorem. We show that the maximal domain of the solution is finite, [0, r_M), terminating at a critical slope u'(r_M) = 1/√3.

What carries the argument

Fixed-point theorem applied to the integral equation for the radial profile u(r) derived from the minimal-resistance variational problem with exponential density.

If this is right

The optimal radial profile is regular up to its maximal finite radius.
The domain of definition ends precisely when the slope condition u' = 1/sqrt(3) is attained.
Existence holds in a neighborhood of the axis for the given exponential density law.
The critical slope value marks the boundary beyond which the radial solution cannot be continued.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The finite termination suggests that, under exponential stratification, no arbitrarily wide minimal-resistance bodies exist in radial symmetry.
The same fixed-point technique might adapt to other radially symmetric density profiles, such as power-law decays.
Physically, the critical slope may correspond to the point where the body surface becomes too steep for the resistance model to remain valid.
Non-radial perturbations could potentially allow solutions to extend farther than the radial case permits.

Load-bearing premise

The fluid density decreases exactly exponentially with altitude and the search is restricted to radial solutions whose governing nonlinearity permits a fixed-point argument.

What would settle it

A numerical integration or shooting method for the radial ODE that either produces a C^1 solution extending past the radius where the slope would equal 1/sqrt(3) or fails to locate any solution satisfying the initial conditions u(0)=u'(0)=0 near the origin.

Figures

Figures reproduced from arXiv: 2605.11641 by Rafael L\'opez.

read the original abstract

This paper investigates the Newton's problem of minimal resistance for a body moving through a fluid whose density decreases exponentially with altitude. We prove the local existence and regularity of radial solutions $u(r)$ satisfying the initial conditions $u(0)=u'(0)=0$ using a fixed-point theorem. We show that the maximal domain of the solution is finite, $[0, r_M)$, terminating at a critical slope $u'(r_M) = \frac{1}{\sqrt{3}}$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper proves local existence of radial solutions for Newton's minimal resistance problem with exponential density decay and shows termination at the classical critical slope.

read the letter

The punchline is that this paper gives a local existence proof for the radial solutions in Newton's minimal resistance problem with exponentially decaying density, and the solution terminates at the same critical slope 1/sqrt(3) as in the constant density version. They derive the Euler-Lagrange ODE from the resistance functional with the variable density, rewrite it as an integral equation to apply a fixed-point theorem, and get local C^2 solutions starting from the axis. The maximal domain is finite because the slope is increasing and cannot cross the point where the leading coefficient vanishes. That point is independent of the density as long as it is positive and smooth, which the exponential is. The work is honest and direct. The fixed-point application is standard but fits the setting, and the argument for finite termination is clean. The stress-test note is right that no major obstruction appears. The limitations are that it stops at local existence. There is no global solution constructed, no numerical plots, and no discussion of the actual resistance achieved or how it compares to constant density. The radial assumption means it is only for bodies of revolution. The exponential density is a modeling choice that works mathematically but may not capture all atmospheric variations. This is for researchers who know the constant-density Newton's problem and want to see the variable-density extension. It would be worth a reading group if your group does variational problems in mechanics. I would not cite it in my own work unless I was building on this exact model. It deserves peer review. The math is coherent and the result is new for this setting, even if the scope is limited.

Referee Report

2 major / 2 minor

Summary. The paper investigates Newton's minimal resistance problem for a body moving through a fluid with density decreasing exponentially with altitude. It proves the local existence and C^2 regularity of radial solutions u(r) to the Euler-Lagrange equation with initial conditions u(0) = u'(0) = 0 using a fixed-point theorem applied to an integral formulation of the ODE. Additionally, it shows that the maximal interval of existence [0, r_M) is finite, with the solution terminating as the slope u'(r_M) approaches 1/√3.

Significance. If the fixed-point argument is complete, this extends the classical constant-density Newton's problem to a variable-density setting relevant for atmospheric or stratified fluids. The persistence of the critical slope termination at 1/√3 indicates robustness of the qualitative behavior. The approach using fixed-point for local existence in the presence of the r=0 singularity is appropriate, and the finite maximal domain provides a concrete prediction for the shape.

major comments (2)

[Proof of local existence (likely §3 or equivalent)] The application of the fixed-point theorem requires explicit verification that the integral operator is a contraction mapping on a suitable complete metric space (e.g., a ball in C[0,R] or appropriate function space). The abstract states the result but the manuscript must provide the estimates showing the Lipschitz constant is less than 1, accounting for the exponential density factor; without these, the hypotheses are not confirmed.
[Maximal domain and termination (likely §4)] The argument for finite r_M relies on showing that u' is strictly increasing and cannot reach or exceed 1/√3. This should be supported by a specific a-priori estimate or comparison with the constant-density case, deriving from the ODE form; the post-hoc observation that the coefficient vanishes at that point needs to be turned into a rigorous bound preventing crossing.

minor comments (2)

[Abstract] The abstract claims 'regularity' but the body specifies C^2; ensure consistency and state the precise Sobolev or classical regularity obtained.
[Notation] Define the density function ρ(z) = exp(-α z) explicitly at the beginning, including the value of α if normalized to 1.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will revise the paper to strengthen the presentation of the proofs.

read point-by-point responses

Referee: [Proof of local existence (likely §3 or equivalent)] The application of the fixed-point theorem requires explicit verification that the integral operator is a contraction mapping on a suitable complete metric space (e.g., a ball in C[0,R] or appropriate function space). The abstract states the result but the manuscript must provide the estimates showing the Lipschitz constant is less than 1, accounting for the exponential density factor; without these, the hypotheses are not confirmed.

Authors: We agree that the contraction-mapping estimates require more explicit verification, including the role of the exponential density. In the revised manuscript we will add a dedicated paragraph (or subsection) that specifies the complete metric space (a closed ball in C[0,R] equipped with the sup norm), derives the Lipschitz constant of the integral operator while keeping the factor e^{-α u(r)} explicit, and shows that this constant is strictly less than 1 for sufficiently small R>0. This will confirm that the Banach fixed-point theorem applies directly and that the local C^2 solution exists on [0,R]. revision: yes
Referee: [Maximal domain and termination (likely §4)] The argument for finite r_M relies on showing that u' is strictly increasing and cannot reach or exceed 1/√3. This should be supported by a specific a-priori estimate or comparison with the constant-density case, deriving from the ODE form; the post-hoc observation that the coefficient vanishes at that point needs to be turned into a rigorous bound preventing crossing.

Authors: We accept that the termination argument needs to be upgraded from an observation to a rigorous a-priori bound. In the revision we will insert a comparison argument: after rewriting the Euler-Lagrange equation in the form (r u' / sqrt(1+(u')^2))' = -r ρ(u) f(u'), we multiply by a positive test function and integrate to obtain an energy-type inequality that prevents u' from attaining or crossing 1/√3 on any finite interval. The same inequality also yields an explicit upper bound on the length of the existence interval, proving that r_M is finite. This approach mirrors the constant-density case while accounting for the strictly positive density factor. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper starts from the variational formulation of the resistance functional with exponentially decaying density, derives the corresponding Euler-Lagrange ODE, rewrites it in integral form to remove the r=0 singularity, and invokes a standard fixed-point theorem on a suitable operator to obtain local existence and C^2 regularity for the radial solution satisfying the initial conditions at r=0. The finite maximal interval is then established by an a priori bound showing u' is strictly increasing and cannot be continued past the value 1/√3 where the coefficient of u'' vanishes; this bound follows directly from the structure of the ODE coefficients (which remain bounded for positive smooth density) and does not rely on any fitted parameters, self-referential definitions, or load-bearing self-citations. All steps are independent applications of classical ODE theory to the derived equation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The proof rests on the standard Banach or Schauder fixed-point theorem applied to an integral operator derived from the Euler-Lagrange equation of the resistance functional; no free parameters or new entities are introduced.

axioms (1)

standard math Banach or Schauder fixed-point theorem
Invoked to obtain local existence of the radial solution to the ODE.

pith-pipeline@v0.9.0 · 5360 in / 1100 out tokens · 50931 ms · 2026-05-13T01:12:14.360718+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

9 extracted references · 9 canonical work pages

[1]

J. D. Anderson Jr., Fundamentals of Aerodynamics. McGraw-Hill Education. New York, 2017

work page 2017
[2]

Buttazzo, A survey on the Newton problem of optimal profiles

G. Buttazzo, A survey on the Newton problem of optimal profiles. In Variational Analysis and Aerospace Engineering, Springer Optim. Appl., Springer, New York, 33 (2009), 33--48

work page 2009
[3]

Buttazzo, V

G. Buttazzo, V. Ferone and B. Kawohl, Minimum problems over sets of concave functions and related questions. Math. Nachr. 173 (1993), 71--89

work page 1993
[4]

Buttazzo and B

G. Buttazzo and B. Kawohl, On Newton's problem of minimal resistance. Math. Intelligencer, 15 (1993), 7--12

work page 1993
[5]

Donald Ahrens, R

C. Donald Ahrens, R. Henson, Meteorology Today: An Introduction to Weather, Climate, and the Environment. Cengage Learning. 12th ed, Boston, 2019

work page 2019
[6]

H. H. Goldstine, A History of the Calculus of Variations from the 17th through the 19th Century. Heidelberg: Springer-Verlag, 1980

work page 1980
[7]

Plakhov, Exterior Billiards: Systems with Impacts Outside Bounded Domains

A. Plakhov, Exterior Billiards: Systems with Impacts Outside Bounded Domains. Springer, Berlin 2012

work page 2012
[8]

Plakhov, D

A. Plakhov, D. Torres, Newton's aerodynamic problem in media of chaotically moving particles. Sb. Math. 196 (2005) 885--933

work page 2005
[9]

J. M. Wallace, P. V. Hobbs, Atmospheric Science: An Introductory Survey. Elsevier, Amsterdam, 2006

work page 2006