## Shape Optimization Course – Day 2

**Speaker – Giuseppe Buttazzo**

The problem of finding a minimal resistance body (due to Newton) consists of finding the shape of a body which travels in a straight line through a fluid when we are given a certain fixed section of it, orthogonal to the flow of the fluid. The classical problem presented here makes some assumptions about the fluid and about the movement, which are not really accurate, taking into account the physics of the fluids, but which turns to be a good approximation in the case the liquid is rare, such as the movement of the airplanes. The assumptions made are:

- the single shock property: every particle which hits the body is reflected and it doesn’t influence the behavior of other particles in the fluid, moreover, if a particle hits the body, it never touches the body after that moment.
- the part of the body below the fixed orthogonal section is neglected, which means that it is considered that its resistance is zero.

Considering a point ( denotes the fixed orthogonal section), the force with which a particle hitting the body in the point holds the body back can be calculated decomposing the normal force as in the figure to be proportional to , where is the angle made by the normal to the surface of the body with the direction of the flow. Using the equality and the fact that the slope of the tangent line to the body we are motivated to chose the resistance functional by taking the mean of all such local resistances over , resulting in the formula .

Our problem is minimizing on some class of admissible functions . Choosing the right class of admissible functions is important, because of the following:

- if we do not impose a boundedness condition on , then we can take to be the function which gives us a very long cone. For this cone, the slope of its generator is very large, and that slope is counted in . Since this is in the denominator, the resistance gets very small, as we take a higher and higher cone. Then the infimum of the resistance is zero, but is never attained, yielding no solution for the problem.
- taking as an admissible class the functions bounded above by some constant , we still don’t have a solution. This is because we can chose a shape with many zig-zag’s, i.e. imagine the picture rotated about a vertical axis. This produces a shape which comes again from a function with gradient of great size (the steeper the slopes of the zig-zag’s, the greater the size). Therefore, in this case the infimum is zero and is never achieved. At this point, one may argue that in reality, the shapes considered are obviously not optimal, because the more zig-zag’s we put, the larger the resistance will be. This strange fact is a consequence of the property of ‘single-shock’, considered in our assumptions, property which is does not hold in reality.

Moreover, there may be some problems due to the functional , which is small as . This means that our functional is not coercive, and we cannot apply the direct methods of the calculus of variation. The admissible class which assures us of the existence of the solutions is .

**Compactness lemma: **For every and every the class is compact with respect to the strong topology of .

*Proof: *Let be a sequence of elements of ; since all are concave, they are locally Lipschitz continuous on , that is , where is a suitable constant. From the fact that , the constants can be chosen independent of ; we can take . Therefore is equi-Lipschitz continuous and equi-bounded on every subset which is relatively compact in . By theorem, is compact relative to the uniform convergence in . By a diagonal argument, we can construct a subsequence of , denoted for simplicity such that uniformly on all compact subsets of for a suitable .

Since the gradients are equibounded on every , by the Lebesgue dominated convergence theorem, in order to conclude the proof it is enough to show that for a.e. .

To do this, fix an integer and a point where all are differentiable (almost all points of are of this kind). The functions are concave, so that we get for every ,

, where is the -th vector in the standard basis in . Passing to the limit for in the last inequality we obtain

.

Taking finalizes the proof:

, which is the desired result.

The result above proves the existence of a solution for the Newton problem in the class . Here are a few facts about the solution of the Newton problems. If is one dimensional then the solution is a triangle for large enough (such that the slope of the side is greater than ) or a trapezoid with lateral sides of slope . For a great period of time, the solutions of the Newton problem in the real case and is a disk, were thought to be radially symmetric. Recently, a few results and approaches prove that this is not the case. Here are some of them:

- P. Guasoni calculated the value of the functional on a shape of the form of a screwdriver (convex hull of the disk and a segment parallel to the plane of the disk) and it turns out that this is smaller than the resistance of the computed rotationally symmetric solution;
- One necessary condition of optimality for a solution of the Newton problem states that if in an open set the function is of class and does not touch the upper bound , then is . This relation is not satisfied by the rotationally symmetric solution; (for a proof see Bucur-Buttazzo,
*Variational Methods in Shape Optimization Problems,*Theorem 2.2.6) - See Theorem 2.23 from the same book as above for another proof;

Again, one necessary condition of optimality is that the set is non-void, for if then we can use a dilation , which decreases the resistance functional. There are some interesting things which can be seen in numerical computations of the optimal solution about the set . All the sets seem to be regular polygons with sides, and the number increases as increases. If we consider all sets of the form and denote with this subclass of , it can be shown that is also compact for the same topologies, and therefore, the Newton problem has a solution in this class. All such solutions are proved to be more optimal than any radially symmetric solution. Even if we know many things about the optimal solution, there a full characterization is not available even if is a disk.

Another well known shape optimization problem is the **Optimal Mixing of Two Materials**. The problem is formulated like this: We have a region and we must fill the region with two materials with conductivities such that contains the first material and contains the second material. We search for an optimal configuration (i.e. an optimal shape ) which is the most performant, with respect to a given cost functional. The volume of each material can also be prescribed. Denoting , the combined conductivity of the two materials, we get the state equation

where is the given source density, and we denote by the unique solution of this equation. It is well known that if we consider an arbitrary cost functional of the form then a general optimal configuration does not exist. However, if we add an additional perimeter penalization like , where , then the optimization problem has a solution.