## A note on Newton's problem of minimal resistance for convex bodies

Alexander Plakhov

Center for R&amp;D in Mathematics and Applications, Department of Mathematics, University of Aveiro, 3810-193, Portugal and Institute for Information Transmission Problems, Moscow, 127051, Russia

## Abstract

We consider the following problem: minimize the functional ∫ Ω f ( ∇ u x ( )) dx in the class of concave functions u : Ω → [0 , M ], where Ω ⊂ R 2 is a convex body and M &gt; 0. If f ( x ) = 1 (1 + / | x | 2 ) and Ω is a circle, the problem is called Newton's problem of least resistance. It is known that the problem admits at least one solution. We prove that if all points of ∂ Ω are regular and (1 + | x | ) f ( x / ) ( | y f | ( y )) → + ∞ as (1 + | x | ) / y | | → 0 then a solution u to the problem satisfies u /floorright ∂ Ω = 0. This result proves the conjecture stated in 1993 in the paper by Buttazzo and Kawohl [6] for Newton's problem.

Mathematics subject classifications: 52A15, 49Q10

Key words and phrases: Convex bodies, Newton's problem of minimal resistance

## 1 Introduction

Consider the following simple mechanical model. A convex body moves forward in a homogeneous medium composed of point particles. The medium is extremely rare, so as mutual interaction of particles is neglected. There is no thermal motion of particles, that is, the particles are initially at rest. When colliding with the body, each particle is reflected elastically. As a result of collisions, there appears the drag force that acts on the body and slows down its motion.

Take a coordinate system in R 3 with the coordinates x = ( x , x 1 2 ) , z connected with the body such that the z -axis is parallel and co-directional to the velocity of the body. Let the upper part of the body's surface be the graph of a concave function u : Ω → R , where Ω is the projection of the body on the x -plane. Then the z -component of the drag

force equals -2 ρv F 2 ( u ), where v is the scalar velocity of the body, ρ is the density of the medium, and

$$F ( u ) = \int _ { \Omega } f ( \nabla u ( x ) ) \, d x,$$

where f ( x ) = 1 / (1 + | x | 2 ). F u ( ) is called resistance of the body. See Fig. 1.

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-30893bf7-b0d4-43ea-b3c6-36fa9c6d8c61.jpg)

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

Figure 1: A convex body moving in a rarefied medium.

Remark 1. In more realistic models, where the particle-body interaction is not perfectly elastic, the resistance also has the form (1) with a different function f . Assume, for example, that the velocity of the reflected particle lies in the tangent plane to the body's surface at the point of impact and is co-directional with the projection of the velocity of incidence on this plane. Assume also that the moduli of the velocities of incidence and reflection coincide. Then we have f ( x ) = 1 -| | x / √ 1 + | x | 2 .

Another generalization of the model concerns a medium with thermal motion of particles. The resistance in such a model is the sum of two integrals ∫ Ω f 1 ( ∇ u 1 ( x )) dx + ∫ Ω f 2 ( ∇ u 2 ( x )) dx , where the functions u 1 and u 2 represent the front and rear parts of the body's surface and the functions f 1 and f 2 are determined by the temperature and composition of the medium and by the velocity of the body. A more detailed description of this model (including the study of some minimal resistance problems) can be found in [22].

Isaac Newton in his Principia [17] studied the problem of minimizing the resistance in the class of rotationally symmetric bodies with fixed length and width. In modern terms, one should minimize the functional (1) in the class of radially symmetric functions u x ( ) = φ ( | x | ), where the unknown function φ : [0 , 1] → R is concave and monotone non-increasing, and satisfies the inequalities 0 ≤ φ ξ ( ) ≤ M . Here M &gt; 0 is the parameter of the problem.

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

The paper by Buttazzo and Kawohl [6] gave rise to studying the problem of minimal resistance in various functional classes: classes of concave (not necessarily radially symmetric) functions [4, 7, 5, 8, 11, 12, 25, 14, 23], classes of functions satisfying the condition of single reflection [4, 9, 10, 21, 20, 16], and some other classes [2, 3, 13, 11, 22, 1, 19, 18, 15].

In this paper we deal with the following minimization problem. Let Ω ⊂ R 2 be a convex body, that is, a compact bounded convex set with nonempty interior, and let f : R 2 → R be a positive continuous function. Fix a value M &gt; 0. The problem is:

$$F ( u ) = \int _ { \Omega } f ( \nabla u ( x ) ) \, d x \rightarrow \inf$$

in the class of concave functions u : Ω → R satisfying 0 ≤ u x ( ) ≤ M . 1

The problem has been intensively studied, especially in the case when Ω is a circle and f ( x ) = 1 / (1 + | x | 2 ). However, despite the apparent simplicity of the statement, this problem is still open. It is well known that the problem has a solution [4, 5].

Here we answer positively the question stated in 1993 in [6]: is it true that if u is a solution to Problem (2) then u /floorright ∂ Ω = 0?

Suppose that the function f satisfies the following condition:

- (A) (1+ | x | ) f ( x ) | y f | ( y ) → + ∞ as (1 + | x | ) / y | | → 0.

Here x = ( x , x 1 2 ) and y = ( y , y 1 2 ) are points of R 2 . By ( x, z ) = ( x , x 1 2 , z ) we denote points of R 3 .

This condition should be understood in the usual sense: for any N &gt; 0 there exists δ &gt; 0 such that for all x, y satisfying (1 + | x | ) / y | | &lt; δ we have (1+ | x | ) f ( x ) | y f | ( y ) &gt; N .

Let x 0 ∈ ∂ Ω. We consider the following condition.

One easily derives that a function f of the form f ( x ) = c x | | -α (1 + o (1)) , x →∞ with c &gt; 0 satisfies Condition A, if α &gt; 1, and does not satisfy, if 0 &lt; α ≤ 1. In particular, f ( x ) = 1 / (1 + | x | 2 ) satisfies this condition.

(B) x 0 is a regular point of the boundary, that is, the support line l at x 0 is unique. Additionally, the endpoints of the segment l ∩ ∂ Ω are also regular.

The segment l ∩ ∂ Ω may degenerate to the point x 0 ; in this case the second part of Condition (B) is not needed.

Theorem 1. Let u be a solution to Problem (2) , where f satisfies Condition A, and let the point x 0 ∈ ∂ Ω satisfy Condition B. Then u x ( 0 ) = 0 .

The following Theorem 2 is a direct consequence of Theorem 1.

Theorem 2. Let u be a solution to Problem (2) , where f satisfies Condition A and all points of ∂ Ω are regular. Then u /floorright ∂ Ω = 0 .

1 Notice that since u is concave, the gradient ∇ u x ( ) exist for almost all x , and therefore, the integral in (2) is well defined.

## 2 Proof of Theorem 1

Assume the contrary: u x ( 0 ) &gt; 0. We are going to come to a contradiction.

Without loss of generality we assume that u is upper semicontinuous. Otherwise we substitute u with the function cl u defined on Ω whose subgraph is the closure of the subgraph of u . Such a function exists and is of course unique, it is concave and upper semicontinuous, and satisfies the inequalities cl u x ( ) ≥ u x ( ) and cl u x ( ) = u x ( ) in the interior of Ω. As a consequence, one has 0 ≤ cl u x ( ) ≤ M and cl u x ( 0 ) ≥ u x ( 0 ) &gt; 0. Proofs of analogous properties for convex , rather than concave , function can be found in the book of Rockafellar [24] (Part II, Section 7).

Since u is upper semicontinuous and l ∩ ∂ Ω is compact, the maximum value of u /floorright l ∩ ∂ Ω is attained at a certain point; without loss of generality we assume that

$$u ( x _ { 0 } ) = \max _ { x \in l \cap \partial \Omega } u ( x ) = \colon z _ { 0 } > 0.$$

Denote by n the outward normal to ∂ Ω at x 0 . The line l is then determined by the equation ( x 0 -x, n ) = 0, and Ω lies in the half-plane ( x 0 -x, n ) ≥ 0; here and in what follows, ( · , · ) means the scalar product.

Take k &gt; 0 and consider the plane of the equation z = k x ( 0 -x, n ). This plane contains the line l × { } 0 , has slope k , and separates the domain of u in the horizontal plane Ω × { } 0 and the vertical segment { x 0 } × [0 , z 0 ]. Consider the auxiliary function u ( k ) ( x ) = min { u x , k x ( ) ( 0 -x, n ) } ; the subgraph of u ( k ) is the part of the subgraph of u located below the plane.

We are going to prove that for k sufficiently large, F u ( ( k ) ) &lt; F u ( ), in contradiction with optimality of u .

Consider the planar convex body

$$\Omega _ { k } & = \{ x \, \colon u ( x ) \geq k ( x _ { 0 } - x, n ) \}. \\. &. &. &. &. &. &. &. &. &.$$

We have ∇ u ( k ) ( x ) = -kn for x in the interior of Ω k . Outside Ω k the function u coincides with u ( k ) , therefore

$$F ( u ) - F ( u ^ { ( k ) } ( x ) ) = \int _ { \Omega _ { k } } f ( \nabla u ( x ) ) \, d x - & \int _ { \Omega _ { k } } f ( u ^ { ( k ) } ( x ) ) \, d x = \int _ { \Omega _ { k } } f ( \nabla u ( x ) ) \, d x - f ( - k n ) \, | \Omega _ { k } | ; \\ \iota _ { \dots \, \dots \, \dots \, \Delta _ { k } } \colon = \int _ { \Omega _ { k } } \alpha _ { \dots \, \Delta _ { k } } + \dots \, \Delta _ { k } \colon = \int _ { \Omega _ { k } } \alpha _ { \dots \, \Delta _ { k } } + \dots \, \Delta _ { k } \colon = \int _ { \Omega _ { k } } \alpha _ { \dots \, \Delta _ { k } } + \dots \, \Delta _ { k } \colon = \int _ { \Omega _ { k } } \alpha _ { \dots \, \Delta _ { k } } + \dots \, \Delta _ { k } \colon = \int _ { \Omega _ { k } } \alpha _ { \dots \, \Delta _ { k } } + \dots \, \Delta _ { k }$$

here and in what follows, | · | means the area of a planar figure. It remains to show that the right hand side of this expression is positive for k sufficiently large.

$$\frac { 1 } { k } \sup _ { x \in \tilde { \Omega } _ { k } } | \nabla u ( x ) | \to 0 \quad \text{ as } \ k \to + \infty ;$$

We are going to find a family of convex bodies Ω ˜ k ⊂ Ω k satisfying the asymptotic relations

$$\inf _ { x \in \tilde { \Omega } _ { k } } \left ( \frac { | \tilde { \Omega } _ { k } | / | \Omega _ { k } | } { ( 1 + | \nabla u ( x ) | ) / k } \right ) \geq \text{const} > 0 \quad \text{for $k$ sufficiently large.}$$

Loosely speaking, we require that first, the maximum of |∇ u | in Ω ˜ k is asymptotically much smaller than the value of |∇ u ( k ) | in Ω k (which is equal to k ) and second, the relative area of Ω ˜ k in Ω k decreases not too rapidly. Let us show that (4) and (5) lead to a contradiction.

The right hand side of (3) can be transformed as follows,

$$x$$

According to (5), the former infimum in (6)is greater than or equal to a positive constant. Taking for any k &gt; 0 a point x k ∈ ˜ Ω , by (4) one has (1 + k |∇ u x ( k ) ) | /k → 0 as k → + ∞ , and hence, by Condition A,

$$I \text{ the right hand side of } ( 3 ) \, can be transformed as tolows, \\ \frac { 1 } { f ( - k n ) \, | \Omega _ { k } | } \Big ( \int _ { \Omega _ { k } } f ( \nabla u ( x ) ) \, d x - f ( - k n ) \, | \Omega _ { k } | \Big ) \geq \\ \frac { 1 } { | \tilde { \Omega } _ { k } | } \int _ { \tilde { \Omega } _ { k } } \frac { | \tilde { \Omega } _ { k } | / | \Omega _ { k } | } { ( 1 + | \nabla u ( x ) | ) / k } \, \frac { ( 1 + | \nabla u ( x ) | ) f ( \nabla u ( x ) ) } { k f ( - k n ) } \, d x - 1 \geq \\ \inf _ { x \in \tilde { \Omega } _ { k } } \Big ( \frac { | \tilde { \Omega } _ { k } | / | \Omega _ { k } | } { ( 1 + | \nabla u ( x ) | ) / k } \Big ) \, \inf _ { x \in \tilde { \Omega } _ { k } } \Big ( \frac { ( 1 + | \nabla u ( x ) | ) f ( \nabla u ( x ) ) } { k f ( - k n ) } \Big ) - 1. \\ \text{xcording to } ( 5 ), the former infimum in } ( 6 ) i s \, \text{greater than or equal to a positive constant.}$$

$$\frac { ( 1 + | \nabla u ( x _ { k } ) | ) f ( \nabla u ( x _ { k } ) ) } { k f ( - k n ) } \to + \infty \quad \text{ as } \ k \to + \infty.$$

Taking the infima over all x k ∈ ˜ Ω , one obtains that the latter infimum in (6) goes to k infinity. It follows that the expression in (6) tends to + ∞ as k → + ∞ and therefore, the right hand side of (3) is positive for k sufficiently large.

It remains to choose Ω ˜ k in such a way that (4) and (5) are satisfied. This will finish the proof of the theorem.

Our construction is illustrated in Figs. 2 and 3 corresponding to the cases when l ∩ ∂ Ω is a point and a line segment, respectively. In both figures the point x 0 is indicated by the letter O and the line OC is orthogonal to l (and its director vector is n ).

Note that Ω k lies in the intersection of Ω with the band { x : 0 ≤ ( x 0 -x, n ) ≤ M/k } . It follows that Ω k lies in an ε k -neighborhood of l ∩ ∂ Ω, with lim k → ∞ + ε k = 0.

Denote

$$\alpha _ { k } \coloneqq \sup _ { x \in \Omega _ { k } } u ( x ) - z _ { 0 }.$$

↦

The sets Ω k , k &gt; 0 form a nested family, Ω k 1 ⊂ Ω k 2 for k 1 ≥ k 2 , therefore the function k → α k is monotone decreasing. Since z 0 = u x ( 0 ) and x 0 ∈ Ω , this function is nonk negative. Further, since z 0 is the maximal value of u /floorright l ∩ ∂ Ω and u is upper semicontinuous, one has α k → 0 as k → + ∞ . It may happen than α k 0 = 0 for a certain value k 0 ; then α k equals zero for all k ≥ k 0 .

Consider the set

$$\omega _ { k } \coloneqq \{ ( x _ { 0 } - x, n ) \, \colon x \in \Omega _ { k } \}.$$

It is a closed segment contained in the positive semiaxis [0 , + ∞ ). Note that the orthogonal projection of Ω k on the line x 0 + ξn, ξ ∈ R (the line OC in Figures 2 and 3) is the line segment x 0 -ω n k .

Taking x = x 0 ∈ Ω , one concludes that the lower endpoint of k ω k is 0, and therefore, ω k has the form

$$\omega _ { k } = [ 0, \, ( z _ { 0 } + \beta _ { k } ) / k ].$$

Since for all x ∈ Ω , ( k x 0 -x, n ) ≤ u x /k ( ) ≤ ( z 0 + α k ) /k , one sees that the upper endpoint of ω k does not exceed ( z 0 + α k ) /k , and so, β k ≤ α k .

Fix ε &gt; 0. For x = x 0 -z 0 -ε k n one has

$$u ( x ) - k ( x _ { 0 } - x, n ) & = \left [ u ( x _ { 0 } - \frac { z _ { 0 } - \varepsilon } { k } n ) - z _ { 0 } \right ] + \varepsilon. \\ \dots \dots \vdots \vdots \dots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots.$$

For k sufficiently large the expression in the right hand side of this formula is positive. This means that x ∈ Ω , and therefore, ( k z 0 -ε /k ) lies in ω k . Hence for k sufficiently large, β k ≥ -ε . It follows that β k → 0 as k → + ∞ .

/Bullet

/Bullet

/Bullet

/Bullet

Figure 2: In this figure, the circular arc through A B , , and O is a part of ∂ Ω, the ellipse is Ω , the shadowed domain is Ω , the points k ˜ k x 0 and x are marked by O and X , respectively. The figure corresponds to the case when l ∩ ∂ Ω is a singleton and the point X lies below the line OC .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-14b01be0-5e04-450d-9e89-b3eff0ff8161.jpg)

/Bullet

/Bullet

Figure 3: The case when l ∩ ∂ Ω is a non-degenerate segment ( PQ in the figure) and X lies below the line OC .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-0eb2b236-56bb-45b8-9b7c-54c445a79597.jpg)

Take t &gt; 0 and consider the segment Ω ∩{ x : ( x 0 -x, n ) = t } . It is orthogonal to the line x 0 + ξn, ξ ∈ R (the line OC in Figs. 2 and 3). If t is sufficiently small, the segment is divided by this line into two non-degenerated segments; say the lower and the upper ones. Let the lengths of the lower and upper segments be a t ( ) and b t ( ). Both functions are concave, non-negative, and monotone increasing for t sufficiently small.

Let a point x is marked by X in Fig. 2 and let t = | C O x | = ( x 0 -x, n ); then a t ( ) = | A C x x | and b t ( ) = | B C x x | .

The following formula will be needed later on. Due to concavity of a , for z 0 2 k ≤ ξ ≤ z 0 + β k k we have a ( z 0 2 k ) ≥ ( 1 -z 0 2 kξ ) a (0) + z 0 2 kξ a ξ ( ), hence

The length of the interval l ∩ ∂ Ω equals a (0) + b (0). If the interval degenerates to the point x 0 , one has a (0) = b (0) = 0. By condition B the endpoints of this interval are regular points of ∂ Ω, hence a ′ (0) = b ′ (0) = + ∞ . It follows that lim t → 0 ( a t /t ( ) ) = + ∞ and lim t → 0 ( b t /t ( ) ) = + ∞ . We assume that k is sufficiently large, so that both functions a t ( ) and b t ( ) are monotone increasing for t ≤ ( z 0 + β k ) /k .

$$a ( \xi ) \leq \frac { 2 k \xi } { z _ { 0 } } \, a \left ( \frac { z _ { 0 } } { 2 k } \right ) - \frac { 2 k \xi - z _ { 0 } } { z _ { 0 } } \, a ( 0 ) \leq 2 \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } \, a \left ( \frac { z _ { 0 } } { 2 k } \right ).$$

A similar formula holds for the function b .

For 0 &lt; θ &lt; 1 consider the linear map T θ : R 2 → R 2 that leaves points of the line OC unchanged and moves other points in the direction orthogonal to OC , so as for all x ∈ R 2 the distance from T x θ to OC is θ times the distance from x to OC . It is defined by the equation T x θ = θx +(1 -θ )(( x -x , n n 0 ) + x 0 ). This map is a compression with the ratio

/Bullet

θ in the direction orthogonal to n .

We take a positive function k → θ k satisfying the conditions

↦

and (in the case when α k is always positive)

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \left ( i \right ) \theta _ { k } \rightarrow 0 \text{ as } k \rightarrow + \infty, \text{ (ii)} \left | \text{$k$} \left ( \frac { Z _ { 0 } } { 2 k } \right ) \theta _ { k } \rightarrow + \infty \text{ and } k b \left ( \frac { Z _ { 0 } } { 2 k } \right ) \theta _ { k } \rightarrow + \infty \text{ as } k \rightarrow + \infty, \text{ (8)} \\ \text{and (in the case when $\alpha$ is always positive)}$$

$$( \text{ii} ) \ \frac { \theta _ { k } } { \alpha _ { k } } \to + \infty \ \text{ as } k \to + \infty.$$

One can take, for example, θ k = 1 / √ ka ( z 0 2 k ) +1 / √ kb ( z 0 2 k ) + √ α k .

Let Ω + k = Ω k ∩ { x : ( x 0 -x, n ) ≥ z 0 2 k } , and let Ω ˜ k = T θ k (Ω + k ); that is, Ω ˜ k is the image of Ω + k under the compression with the ratio θ k . In Figs. 2 and 3, Ω ˜ k is shown shadowed. The set Ω + k is contained between the lines ( x 0 -x, n ) = z 0 + β k k and ( x 0 -x, n ) = z 0 2 k , and the set Ω k \ Ω + k is contained between the line ( x 0 -x, n ) = z 0 2 k and the line l of the equation ( x 0 -x, n ) = 0; see Fig. 4. Let the intersection of Ω k with the line ( x 0 -x, n ) = z 0 2 k be the segment MN . The intersection of Ω k with the line ( x 0 -x, n ) = z 0 + β k k is nonempty; take a point P in this intersection. Let M ′ and N ′ be the points of intersection of l with the lines PM and PN , respectively.

Figure 4: The sets Ω k and Ω k \ Ω . + k

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-abfc7be8-8cc2-4e0e-a577-96908716fe85.jpg)

The area of the triangle PMN equals

$$| \triangle P M N | = \frac { 1 } { 2 } \, \frac { z _ { 0 } + 2 \beta _ { k } } { 2 k } \, | M N |.$$

One easily sees that | M N ′ ′ | = 2( z 0 + β k ) z 0 +2 β k | MN | , and therefore, the area of the trapezoid MNN M ′ ′ is

$$| \Box M N N ^ { \prime } M ^ { \prime } | = \frac { z _ { 0 } } { 2 k } \left ( \frac { 1 } { 2 } + \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } + 2 \beta _ { k } } \right ) | M N |.$$

The set Ω + k contains the triangle PMN , and the set Ω k \ Ω + k is contained in the trapezoid MNN M ′ ′ , therefore

$$& \frac { | \Omega _ { k } \, \sum \Omega _ { k } ^ { + } | } { | \Omega _ { k } ^ { + } | } = \frac { | \Omega _ { k } | } { | \Omega _ { k } ^ { + } | } - 1 \leq \frac { z _ { 0 } ( 3 z _ { 0 } + 4 \beta _ { k } ) } { ( z _ { 0 } + 2 \beta _ { k } ) ^ { 2 } } \quad \Rightarrow \quad \frac { | \Omega _ { k } | } { | \Omega _ { k } ^ { + } | } \leq \frac { 4 ( z _ { 0 } + \beta _ { k } ) ^ { 2 } } { ( z _ { 0 } + 2 \beta _ { k } ) ^ { 2 } }. \\ & \sigma, \quad \text{$\quad $},$$

Thus, one has

$$& \frac { | \tilde { \Omega } _ { k } | } { | \Omega _ { k } | } = \theta _ { k } \frac { | \Omega _ { k } ^ { + } | } { | \Omega _ { k } | } \geq \frac { \theta _ { k } } { 4 } \, \frac { ( z _ { 0 } + 2 \beta _ { k } ) ^ { 2 } } { ( z _ { 0 } + \beta _ { k } ) ^ { 2 } } = \frac { \theta _ { k } } { 4 } \left ( 1 + o ( 1 ) \right ) \quad \text{as } k \to + \infty. \quad \text{(10)} \\ & \quad \text{$\sim}$$

Let x ∈ ˜ Ω k be a regular point of u ; it is indicated by letter X in Figs. 2 and 3. If the tangent plane to the graph of u at x is not horizontal, then its intersection with the horizontal plane z = 0 is a straight line, say l x , and

$$| \nabla u ( x ) | = \frac { u ( x ) } { \text{dist} ( x, l _ { x } ) }.$$

Of course l x does not intersect the interior of Ω.

The intersection of this tangent plane with the vertical plane through the line XO (and therefore, through the points ( x, 0) and ( x , 0 0)) is a straight line, say σ x , that contains the point ( x, u x ( )) and lies above the point ( x , u x 0 ( 0 )); see Fig. 5. The slope of this line in the direction - -→ XO is greater than or equal to ( u x ( 0 ) -u x ( )) / x | 0 -x | .

It may happen that (i) the line σ x intersects the ray with the vertex at x and with the director vector x 0 -x (the ray XO in Figs. 2, 3, and 5), or (ii) it does not intersect this ray. If the case (i) is realized, one has u x ( ) &gt; u x ( 0 ) (and therefore α k &gt; 0), and the

/Bullet

/Bullet

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-7589f345-e213-4b73-ade8-dbacf04f3f51.jpg)

/Bullet

/Bullet

/Bullet

/Bullet

Figure 5: The section of graph( u ) and the tangent plane to graph( u ) at x by the vertical plane through X and O . (i) The line σ x intersects the ray XO ; (ii) σ x does not intersect this ray.

slope of the line σ x is greater than or equal to -α / x k | 0 -x | . It follows that the distance between x 0 and the point of intersection of σ x with the ray XO is greater than or equal

/Bullet

to z 0 α k | x 0 - | x . As a consequence, the point of intersection lies on the ray behind the point x 0 + z 0 α k ( x 0 -x ) (the point F in Figs. 2, 3 and 5). For the point F we have the proportion

$$| O F | = \frac { z _ { 0 } } { \alpha _ { k } } \, | X O |.$$

Hence l x does not intersect the interior of the convex hull of Ω and the segment [ XF ],

$$l _ { x } \cap & \int \left ( \text{Conv} ( \Omega \cup [ X F ] ) \right ) = \emptyset. & ( 1 3 ) \\ \text{is realized then $1$ does not intersect the interior of the convev hull of}$$

If the case (ii) is realized then l x does not intersect the interior of the convex hull of the union of Ω and the ray XO . In particular, formula (13) remains true, where in the case α k = 0 the point F can be imagined as the infinitely remote point on the ray and [ XF ] should be understood as the ray XO .

Let A x and B x be the endpoints of the segment { x ′ : ( x , n ′ ) = ( x, n ) } ∩ Ω (see Fig. 2), and let ∂ x Ω be the part of ∂ Ω on the left of the line A B x x . In other words, ∂ x Ω := ∂ Ω ∩ { x ′ : ( x , n ′ ) ≤ ( x, n ) } . Let [ A F x ] and [ B F x ] denote the corresponding segments, A F x and B F x be the lines containing these segments, and [ A FB x x ] be the union of the segments. (If α k = 0, [ A F x ] and [ B F x ] should be understood as the rays co-directional with the ray XO with the vertices at A x and B x , respectively, and [ A FB x x ] as the union of these rays.)

The open domain bounded by the union of the curve ∂ x Ω and the broken line (or the union of rays) [ A FB x x ] (see Fig. 6) is contained in Conv(Ω ∪ [ XF ]), and therefore, does not intersect the line l x . Therefore we have

$$\text{dist} ( x, l _ { x } ) \geq \text{dist} ( x, \partial _ { x } \Omega \cup [ A _ { x } F B _ { x } ] ) = \min \{ \text{dist} ( x, \partial _ { x } \Omega ), \text{dist} ( x, [ A _ { x } F ] ), \text{dist} ( x, [ B _ { x } F ] ) \}.$$

Take t 0 sufficiently small, so as the functions a t ( ) and b t ( ) are monotone increasing for 0 ≤ t ≤ t 0 . Consider the rectangle located on the left of the line A B x x such that the segment [ A B x x ] is one of its vertical sides (see Figs. 6 and 3) and the length of its horizontal sides equals t 0 / 2. One easily sees that for k sufficiently large (namely, for k satisfying ( z 0 + β k ) /k ≤ t 0 / 2) and for x ∈ ˜ Ω k this rectangle is contained in Ω.

The distance from x to the union of the three sides of the rectangle (except A B x x ) equals min { t 0 / 2 , | XA , XB x | | x |} . It follows that dist( x, ∂ x Ω) ≥ min { t 0 / 2 , | xA x | , | xB x |} ≥ min { t 0 / 2 , dist( x, [ A F x ]) , dist( x, [ B F x ]) } .

Hence

$$\text{dist} ( x, l _ { x } ) \geq \min \{ t _ { 0 } / 2, \text{dist} ( x, [ A _ { x } F ] ), \text{dist} ( x, [ B _ { x } F ] ) \}. \quad \text{(14)} \\ \text{$\mu$ $\nu$} \quad \tau k \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{matrix} \quad \text{$z_{0}$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\nu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$} \quad \text{$\mu$}$$

Let C = C k indicate the point x 0 -z 0 2 k n . Denote by C x the projection of X on the line CO , and by E x and E the projections of F on the lines A B x x and l , respectively 2 (see Fig. 2). Note that | CO | = z 0 2 k .

2 If α k = 0, the points E x and E are not defined.

/Bullet

/Bullet

Figure 6: The domain bounded by the curve ∂ x Ω and the broken line [ A FB x x ]. In this figure, ∂ x Ω is the arc of circumference with the endpoints A x and B x located to the left of the line A B x x .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-4ac0c854-a125-4d81-8f69-0683106280e8.jpg)

Assuming that α k &gt; 0, using (12) and taking into account that the triangles XC O x and OEF are similar, one obtains

$$| E _ { x } F | & \geq | E F | = | C _ { x } O | \, \frac { | O F | } { | X O | } \geq \frac { z _ { 0 } } { 2 k } \, \frac { z _ { 0 } } { \alpha _ { k } } = \frac { z _ { 0 } ^ { 2 } } { 2 \alpha _ { k } k }. \\. & \quad \sim \, \left ( \alpha \cap \right ) \quad, \quad \, z _ { \alpha } \, \gtrsim \, \gtrsim z _ { \alpha } + \beta _ { k } \quad,$$

Denote ξ = | C O x | ; we have z 0 2 k ≤ ξ ≤ z 0 + β k k and

$$| A _ { x } C _ { x } | = a ( \xi ) \geq a ( \frac { z _ { 0 } } { 2 k } ).$$

Using formula (7), one obtains

$$| A _ { x } C _ { x } | \leq 2 \, \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } \, a ( \frac { z _ { 0 } } { 2 k } ).$$

Consider two cases.

(a) X lies below the line OC , that is, in the same half-plane bounded by OC as A x . First assume that (i) the line σ x intersects the ray XO . Using (17) we obtain

$$| X C _ { x } | \leq \theta _ { k } | A _ { x } C _ { x } | \leq 2 \theta _ { k } \, \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } \, a ( \frac { z _ { 0 } } { 2 k } ).$$

It follows that

Further, using (16) one has

$$& \frac { | X E _ { x } | } { | E _ { x } F | } = \frac { | X C _ { x } | } { | C _ { x } O | } \leq \frac { | X C _ { x } | } { | C O | } \leq 4 \theta _ { k } \frac { k a \left ( \frac { z _ { 0 } } { 2 k } \right ) } { z _ { 0 } } \ \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } }. \\ & \text{her. using (16) one has}$$

$$| A _ { x } X | = | A _ { x } C _ { x } | - | C _ { x } X | \geq ( 1 - \theta _ { k } ) | A _ { x } C _ { x } | \geq ( 1 - \theta _ { k } ) a ( \frac { z _ { 0 } } { 2 k } ).$$

Let us evaluate the distance h between X and the line A F x . Consider the triangle A XF x . The area S of this triangle can be calculated in two ways,

$$S = \frac { 1 } { 2 } \left | A _ { x } X \right | \left | E _ { x } F \right | \text{ \ and \ } S = \frac { 1 } { 2 } \, h \left | A _ { x } F \right | \leq \frac { 1 } { 2 } \, h \left ( \left | A _ { x } X \right | + \left | X E _ { x } \right | + \left | E _ { x } F \right | \right ),$$

and using inequalities (15), (18), and (19) one obtains

$$h \geq \frac { | A _ { x } X | \, | E _ { x } F | } { | A _ { x } X | + | X E _ { x } | + | E _ { x } F | } & = \frac { 1 } { \frac { 1 } { | E _ { x } F | } + \frac { 1 } { | A _ { x } X | } \left ( 1 + \frac { | X E _ { x } | } { | E _ { x } F | } \right ) } \\ & 1 & 1$$

$$& \geq \frac { 1 } { \frac { 2 \alpha _ { k } k } { z _ { 0 } ^ { 2 } } + \frac { 1 } { 1 - \theta _ { k } } \ \frac { 1 } { a \left ( \frac { z _ { 0 } } { 2 k } \right ) } \left ( 1 + 4 \theta _ { k } \, \frac { k a \left ( \frac { z _ { 0 } } { 2 k } \right ) } { z _ { 0 } } \, \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } \right ) } } = \frac { 1 } { k \theta _ { k } } \ \frac { 1 } { \frac { 2 \alpha _ { k } } { z _ { 0 } ^ { 2 } \theta _ { k } } + \frac { 1 } { 1 - \theta _ { k } } \ \frac { 1 } { k a \left ( \frac { z _ { 0 } } { 2 k } \right ) \theta _ { k } } + \frac { 4 } { 1 - \theta _ { k } } \ \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } ^ { 2 } } }. \\ & \text{It follows that} & \quad \ 1 \quad.$$

$$\mathfrak { i t } ( x, [ A _ { x } F ] ) \geq \text{dist} ( x, A _ { x } F ) \geq \frac { 1 } { k \theta _ { k } } \, H _ { k } ^ { 1 },$$

It follows that where

Note that H 1 k is chosen to be invariant with respect to exchanging a and b . Due to (8) and (9), the denominator in this expression tends to 4 /z 0 , hence H 1 k → z / 0 4 as k → + ∞ .

$$H _ { k } ^ { 1 } & = \frac { 1 } { \frac { 2 \alpha _ { k } } { z _ { 0 } ^ { 2 } \theta _ { k } } + \frac { 1 } { 1 - \theta _ { k } } \left ( \frac { 1 } { k a \left ( \frac { z _ { 0 } } { 2 k } \right ) \theta _ { k } } + \frac { 1 } { 1 - \theta _ { k } } \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } ^ { 2 } } } }. \\ \intertext { t a t h } H _ { k } ^ { 1 } & \text{ is chosen to be invariant with respect to exchanging a and b. } \ D u e \, to \, ( 8 )$$

Now assume that (ii) σ x does not intersect the ray XO . If α k &gt; 0 then the argument above in the case (a) remains valid, and formulae (20) and (21) hold true. If α k = 0, the argument is valid when F is replaced with an arbitrary point x 0 + z 0 α ( x 0 -x ), α &gt; 0 on the ray XO . It follows that inequality (20) is true when α k is replaced with arbitrary α &gt; 0 in (21). Taking the limit α → 0, one concludes that (20) is also true when we have α = α k = 0 in (21).

- (b) X lies above the line OC . Here again we consider two cases.

(b ) 1 E x lies below A x or (in the case α k = 0) does not exist; see Fig. 7 (b ). 1 In this case, using (16) and (21), we have

$$\text{dist} ( x, [ A _ { x } F ] ) = | X A _ { x } | \geq | A _ { x } C _ { x } | = a ( \xi ) \geq a ( \frac { z _ { 0 } } { 2 k } ) > \frac { 1 } { k \theta _ { k } } \, H _ { k } ^ { 1 }.$$

(b ) 2 E x lies between the points A x and C x ; see Fig. 7 (b 2 ). Due to (17) and (16), one has

$$| A _ { x } E _ { x } | \leq | A _ { x } C _ { x } | \leq 2 \, \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } \, a ( \frac { z _ { 0 } } { 2 k } ) \quad \text{ and } \quad | A _ { x } X | \geq | A _ { x } C _ { x } | \geq a ( \frac { z _ { 0 } } { 2 k } ).$$

/Bullet

/Bullet

/Bullet

x

/Bullet

Figure 7: X lies above the line OC . The two figures correspond to the cases (b ) when 1 E x lies below A x and (b ) when 2 E x lies between C x and A x .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01042v2-figure-b222438d-aa8c-4d4d-b31b-ba6e07d4fe88.jpg)

Let us consider again the area S of the triangle A XF x . One has

$$S = \frac { 1 } { 2 } \left | A _ { x } X \right | \left | E _ { x } F \right | \text{ \ and \ } S = \frac { 1 } { 2 } \, h \left | A _ { x } F \right | \leq \frac { 1 } { 2 } \, h \left ( \left | A _ { x } E _ { x } \right | + \left | E _ { x } F \right | \right ),$$

and using (15) and (22) one obtains

$$h \geq \frac { | A _ { x } X | \, | E _ { x } F | } { | A _ { x } E _ { x } | + | E _ { x } F | } = \frac { | A _ { x } X | } { \frac { | A _ { x } E _ { x } | } { | E _ { x } F | } + 1 } \geq \frac { a \left ( \frac { z _ { 0 } } { 2 k } \right ) } { 2 \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } } } a \left ( \frac { z _ { 0 } } { 2 k } \right ) \frac { 2 \alpha _ { k } k } { z _ { 0 } ^ { 2 }$$

where

H 2 k is also invariant with respect to exchanging a and b .

$$H _ { k } ^ { 2 } & = \frac { 1 } { \frac { 4 \alpha _ { k } } { \theta _ { k } } \frac { z _ { 0 } + \beta _ { k } } { z _ { 0 } ^ { 3 } } + \frac { 1 } { k a \left ( \frac { z _ { 0 } } { 2 k } \right ) \theta _ { k } } + \frac { 1 } { k b \left ( \frac { z _ { 0 } } { 2 k } \right ) \theta _ { k } } } \rightarrow + \infty \text{ as } k \rightarrow + \infty. \\ \text{s also invariant with respect to exchanging $a$ and $b$}. \quad \quad \quad$$

Thus, in all cases one has dist( x, [ A F x ]) ≥ 1 kθ k H , k where H k = min { H , H 1 k 2 k } → z / 0 4 as k → + ∞ .

$$\frac { 1 } { k \theta _ { k } } = \frac { a \left ( \frac { z _ { 0 } } { 2 k } \right ) } { k \, a \left ( \frac { z _ { 0 } } { 2 k } \right ) \theta _ { k } } \rightarrow 0 \quad \text{as } k \rightarrow + \infty, \\ \intertext { large enough we have }$$

In a completely similar way one derives the estimate dist( x, [ B F x ]) ≥ 1 kθ k H k . Using that the function a is bounded and that by (8) k a ( z 0 2 k ) θ k goes to infinity, we have hence by (14) for k large enough we have

$$\text{dist} ( x, l _ { x } ) \geq \min \left \{ \frac { t _ { 0 } } { 2 }, \, \frac { 1 } { k \theta _ { k } } H _ { k } \right \} = \frac { 1 } { k \theta _ { k } } H _ { k }.$$

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

/Bullet

By (11), for a regular x ∈ ˜ Ω k one has

$$\frac { 1 } { k } \left | \nabla u ( x ) \right | = \frac { 1 } { k } \ \frac { u ( x ) } { \text{dist} ( x, l _ { x } ) } \leq \frac { M } { H _ { k } } \, \theta _ { k } \rightarrow 0 \ \text{ as } k \rightarrow + \infty.$$

Thus, equation (4) is true. Further, using (10) and (23) one obtains

$$\frac { | \tilde { \Omega } _ { k } | / | \Omega _ { k } | } { ( 1 + | \nabla u ( x ) | ) / k } \geq \frac { 1 } { 4 } \, \frac { ( z _ { 0 } + 2 \beta _ { k } ) ^ { 2 } } { ( z _ { 0 } + \beta _ { k } ) ^ { 2 } } \, \frac { 1 } { \frac { 1 } { k \theta _ { k } } + \frac { M } { H _ { k } } } = \frac { z _ { 0 } } { 1 6 M } \, ( 1 + o ( 1 ) ) \quad \text{as } k \to + \infty. \\ -. \quad \, \cdot \quad \, \cdot \quad \, \cdot \quad \, \cdot \quad \, \cdot \quad \, \cdot$$

Thus, equation (5) is also true. This completes the proof of the theorem.

## Acknowledgements

This work was supported by Foundation for Science and Technology (FCT), within project UID/MAT/04106/2019 (CIDMA). I am very grateful to G. Wachsmuth for a discussion, which helped to fix an error in the text.

## References

- [1] A Aleksenko and A Plakhov. Bodies of zero resistance and bodies invisible in one direction . Nonlinearity 22 , 1247-1258 (2009).
- [2] M. Belloni and B. Kawohl. A paper of Legendre revisited . Forum Math. 9 , 655-668 (1997).
- [3] M. Belloni and A. Wagner. Newtons problem of minimal resistance in the class of bodies with prescribed volume . J. Convex Anal. 10 , 491500 (2003).
- [4] G. Buttazzo, V. Ferone, B. Kawohl. Minimum problems over sets of concave functions and related questions . Math. Nachr. 173 , 71-89 (1995).
- [5] G. Buttazzo and P. Guasoni. Shape optimization problems over classes of convex domains , J. Convex Anal. 4 , No.2, 343-351 (1997).
- [6] G. Buttazzo, B. Kawohl. On Newton's problem of minimal resistance . Math. Intell. 15 , 7-12 (1993).
- [7] F. Brock, V. Ferone and B. Kawohl. A symmetry problem in the calculus of variations . Calc. Var. Partial Differ. Equ. 4 , 593-599 (1996).
- [8] G. Carlier and T. Lachand-Robert. Convex bodies of optimal shape . J. Convex Anal. 10 , No. 1, 265-273 (2003).

| [9]   | M. Comte, T. Lachand-Robert. Newton's problem of the body of minimal resistance under a single-impact assumption . Calc. Var. Partial Differ. Equ. 12 , 173-211 (2001).          |
|-------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [10]  | M. Comte, T. Lachand-Robert. Existence of minimizers for Newton's problem of the body of minimal resistance under a single-impact assumption . J. Anal. Math. 83 313-335 (2001). |
| [11]  | B. Kawohl. Some nonconvex shape optimization problems . Optimal shape design (Tr´ia, o 1998), 7-46, Lecture Notes in Math., 1740 , Springer, Berlin, 2000.                       |
| [12]  | T. Lachand-Robert and E. Oudet. Minimizing within convex bodies using a convex hull method . SIAM J. Optim. 16 , 368-379 (2006).                                                 |
| [13]  | T. Lachand-Robert, M. A. Peletier. Newton's problem of the body of minimal resis- tance in the class of convex developable functions . Math. Nachr. 226 , 153-176 (2001).        |
| [14]  | L. V. Lokutsievskiy and M. I. Zelikin. Hessian measures in the aerodynamic Newton problem . J. Dynam. Control Syst. 24 , 475-495 (2018).                                         |
| [15]  | L. V. Lokutsievskiy and M. I. Zelikin. The analytical solution of newton's aero- dynamic problem without the assumption of axial symmetry . arXiv:1905.02028 [math.OC]           |
| [16]  | E. Mainini, M. Monteverde, E. Oudet, and D. Percivale. The minimal resistance problem in a class of non convex bodies . ESAIM: COCV bf 25, pp. 1-36 (2019).                      |
| [17]  | I. Newton. Philosophiae naturalis principia mathematica . 1687.                                                                                                                  |
| [18]  | A. Plakhov. Billiards and two-dimensional problems of optimal resistance . Arch. Ra- tion. Mech. Anal. 194 , 349-382 (2009).                                                     |
| [19]  | A. Plakhov. Optimal roughening of convex bodies . Canad. J. Math. 64 , 1058-1074 (2012).                                                                                         |
| [20]  | A. Plakhov. Problems of minimal resistance and the Kakeya problem. SIAM Review 57 , 421-434 (2015).                                                                              |
| [21]  | A. Plakhov. Newtons problem of minimal resistance under the single-impact assump- tion . Nonlinearity 29 , 465-488 (2016).                                                       |
| [22]  | A. Plakhov and D. Torres. Newton's aerodynamic problem in media of chaotically moving particles . Sbornik: Math. 196 , 885-933 (2005).                                           |
| [23]  | A. Plakhov. Local properties of the surface measure of convex bodies . J. Convex Anal. (in press)                                                                                |

- [24] R. T. Rockafellar. Convex Analysis. Vol. 28 of Princeton Math. Series, Princeton Univ. Press, 1970 (470 pp).
- [25] G. Wachsmuth. The numerical solution of Newtons problem of least resistance . Math. Program. A 147 , 331-350 (2014).
Alexander Plakhov
A note on Newton's problem of minimal resistance for convex bodies
We consider the following problem: minimize the functional $\int_\Omega
f(\nabla u(x))\, dx$ in the class of concave functions $u: \Omega \to [0,M]$,
where $\Omega \subset \mathbb{R}^2$ is a convex body and $M > 0$. If $f(x) =
1/(1 + |x|^2)$ and $\Omega$ is a circle, the problem is called Newton's problem
of least resistance. It is known that the problem admits at least one solution.
We prove that if all points of $\partial\Omega$ are regular and
${|x|f(x)}/(|y|f(y)) \to +\infty$ as $|x|/|y| \to 0$ then a solution $u$ to the
problem satisfies $u\rfloor_{\partial\Omega} = 0$. This result proves the
conjecture stated in 1993 for Newton's problem.