Rapport:

Optimisation fonctionnelle

TENBIH Wissal, ELFILALI Othmane

Mécanique des Fluides, Energétique et Environnement - MF2E 1A

Mai 2024

1
Table des matières
1 Introduction 3

2 Miroir liquide : 3
2.1 Modélisation du problème . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Expression de l’énergie potentielle . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.3 Expression de l’énergie cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.4 Expression de la fonctionnelle Λ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.5 Équation d’Euler-Lagrange pour f . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.6 Commentaire sur la forme de la surface libre . . . . . . . . . . . . . . . . . . . . . . 4
2.7 Calcul du volume de liquide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.8 Méthode du multiplicateur de Lagrange . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.9 Équation d’Euler-Lagrange modifiée . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.10 Développement du calcul . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3 Surface d’aire minimale : 7

3.1 Positionnement du problème : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.2 Expression de la fonctionnelle A[f ] : . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.3 Équation différentielle pour minimiser A[f ] : . . . . . . . . . . . . . . . . . . . . . . 7
x− L

3.4 Vérifions si f (x) = c0 cosh c0 2 est une solution lorsque r0 = rL : . . . . . . . . . 8
3.5 Vérifions que dans le cas de la surface considérée ici la courbure moyenne est nulle : 9
3.6 Proposition d un algorithme, permettant de trouver (numériquement) la forme de
l’interface minimisant la surface : . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 2

1 Introduction
Ce rapport explore deux problèmes d’optimisation fonctionnelle : le miroir liquide et la surface
d’aire minimale. Dans le premier cas, nous examinons la forme optimale d’un miroir liquide dans
un télescope, tandis que dans le second, nous étudions la minimisation de l’aire d’une surface de
révolution. Ces exemples illustrent l’application des méthodes d’optimisation dans des domaines
aussi variés que l’astronomie et la physique des interfaces. En combinant théorie et application,
nous découvrirons comment ces méthodes permettent de résoudre des problèmes concrets.

2 Miroir liquide :
2.1 Modélisation du problème
Dans cet exercice, nous vérifions si la rotation de la cuve conduit la surface libre du liquide
à adopter la forme d’un paraboloı̈de en utilisant le principe de moindre action. L’objectif est de
déterminer la forme de la configuration qui minimise l’action, donnée par l’intégrale de la différence
entre l’énergie cinétique (Ec ) et potentielle (Ep ) du système :

Figure 1 – Schéma de la cuve en rotation

Les paramètres physiques intervenant dans ce problème sont : - La densité du liquide ρ - La

vitesse angulaire de rotation ω - La gravité g - Le rayon de la cuve R - Le volume de liquide V0

2.2 Expression de l’énergie potentielle

L’énergie potentielle Ep d’un élément de volume en coordonnées cylindriques (dv = r dr dθ dz)
est donnée par :
Z Z 2π Z R Z f (r)
Ep = ρgz dv = ρg z r dz dr dθ
V 0 0 0
En intégrant par rapport à z, on obtient :
Z 2π Z R  2 f (r) Z 2π Z R
z f (r)2
Ep = ρg r dr dθ = ρg r dr dθ
0 0 2 0 0 0 2
Puis en intégrant par rapport à θ, on obtient :
Z R
Ep = πρg f (r)2 r dr
0

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 3

2.3 Expression de l’énergie cinétique
L’énergie cinétique Ec d’un élément de volume en rotation avec une vitesse angulaire ω est
donnée par :
ρω 2 2π R f (r) 2
Z Z Z Z
1 2
Ec = ρ(ωr) dv = r r dz dr dθ
2 V 2 0 0 0
En intégrant par rapport à z, on obtient :
ρω 2 2π R 3
Z Z
Ec = r f (r) dr dθ
2 0 0
Puis en intégrant par rapport à θ, on obtient :
Z R
2
Ec = πρω r3 f (r) dr
0

2.4 Expression de la fonctionnelle Λ

La fonctionnelle Λ à minimiser est donnée par la différence entre l’énergie cinétique et l’énergie
potentielle :  Z R Z R 
2 3 2
Λ = Ec − Ep = πρ ω r f (r) dr − g f (r) r dr
0 0
La forme de la surface libre f (r) sera influencée par les paramètres ρ, ω, et g selon cette fonction-
nelle.

2.5 Équation d’Euler-Lagrange pour f

L’équation d’Euler-Lagrange pour minimiser Λ est donnée par :
 
∂Λ d ∂Λ
− =0
∂f dr ∂f ′
df
où f ′ = dr
. En appliquant cette formule, on obtient l’équation différentielle gouvernant f (r).

2.6 Commentaire sur la forme de la surface libre

La forme de la surface libre f (r) est déterminée par l’équilibre entre les forces centrifuges
(proportionnelles à ω 2 r2 ) et la force de gravité (proportionnelle à gf (r)). Cela conduit typiquement
à une forme parabolique.

2.7 Calcul du volume de liquide

Le volume de liquide V0 est donné par l’intégrale du volume sous la surface f (r) :
Z 2π Z R Z f (r) Z R
V0 = r dz dr dθ = 2π rf (r) dr
0 0 0 0

2.8 Méthode du multiplicateur de Lagrange

Pour minimiser Λ tout en satisfaisant la contrainte de conservation du volume, on utilise un
multiplicateur de Lagrange λ. La fonctionnelle à minimiser devient :
Z R 
′ V0
Λ =Λ+λ rf (r) dr −
0 2π

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 4

2.9 Équation d’Euler-Lagrange modifiée
L’équation d’Euler-Lagrange correspondante est :
∂Λ′ d ∂Λ′
 
− =0
∂f dr ∂f ′
En résolvant cette équation, nous déterminons la forme de f (r) et la valeur du multiplicateur de
Lagrange λ qui satisfait la contrainte de volume.

2.10 Développement du calcul

Nous allons maintenant développer les calculs pour obtenir les expressions explicites de f (r) et
de λ.
Pour Λ, nous avons :
 Z R Z R 
2 3 2
Λ = πρ ω r f (r) dr − g f (r) r dr
0 0

et pour la contrainte de volume : Z R

V0 = 2π rf (r) dr
0
La fonctionnelle Λ′ devient alors :
 Z R Z R   Z R 
′ 2 3 2
Λ = πρ ω r f (r) dr − g f (r) r dr + λ 2π rf (r) dr − V0
0 0 0

L’équation d’Euler-Lagrange pour f (r) est :

∂Λ′ ∂Λ′
 
d
− =0
∂f dr ∂f ′
∂Λ′
Calculons ∂f
:
∂Λ′
= πρ ω 2 r3 − 2gf (r)r + 2πλr


∂f
∂Λ′
Étant donné que ∂f ′
= 0, l’équation d’Euler-Lagrange devient :

πρ ω 2 r3 − 2gf (r)r + 2πλr = 0

Simplifions cette expression :

πρω 2 r3 − 2πρgf (r)r + 2πλr = 0

Divisons par πr pour obtenir :

ρω 2 r2 − 2ρgf (r) + 2λ = 0
Réarrangeons pour obtenir f (r) :
ρω 2 r2 + 2λ
f (r) =
2ρg

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 5

 Pour déterminer λ, nous utilisons la contrainte de volume :
Z R
V0 = 2π rf (r) dr
0

R
ρω 2 r2 + 2λ
Z  
V0 = 2π r dr
0 2ρg
Z R
2π
V0 = r(ρω 2 r2 + 2λ) dr
2ρg 0
 Z R Z R 
π 2 3
V0 = ρω r dr + 2λ r dr
ρg 0 0
Calculons les intégrales :
R
R4
Z
r3 dr =
0 4
R
R2
Z
r dr =
0 2
En substituant, nous obtenons :
4
R2
 
π 2R
V0 = ρω + 2λ
ρg 4 2
π ρω 2 R4
 
2
V0 = + λR
ρg 4
πR2 ρω 2 R2
 
V0 = +λ
ρg 4
Isolons λ :
ρg ρω 2 R2
V0 = +λ
πR2 4
V0 ρg ρω 2 R2
λ= −
πR2 4
La forme finale de f (r) est donc :
 
V0 ρg ρω 2 R2
ρω 2 r2 + 2 πR2
− 4
f (r) =
2ρg
2V0 g ω 2 R2
ω 2 r2 + πR2
− 2
f (r) =
2g

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 6

3 Surface d’aire minimale :
3.1 Positionnement du problème :
Dans ce deuxième exercice, nous étudions une surface de révolution générée par une courbe
définie par la fonction r = f (x). Cette surface est caractérisée par les conditions aux limites
r0 = f (0) et rL = f (L), où r0 et rL sont des valeurs connues.
L’objectif est de déterminer la forme de la surface qui minimise son aire tout en respectant
ces conditions aux limites. Ce problème trouve une analogie physique dans la forme d’un film de
savon tendu entre deux anneaux. En équilibre, ce film adopte la forme qui minimise son énergie
de surface, proportionnelle à l’aire de la surface. Ainsi, la résolution de ce problème mathématique
nous permet de comprendre et de prédire la configuration optimale d’un tel film de savon.

Figure 2 – Illustration de la surface de révolution générée par la courbe r = f (x).

3.2 Expression de la fonctionnelle A[f ] :

La surface d’aire d’une surface de révolution générée par une courbe r = f (x) autour de l’axe x
peut être calculée en intégrant les éléments de surface circulaires sur toute la longueur L. L’élément
de surface d’une bande circulaire à une distance x de largeur infinitésimale dx et de rayon f (x) est
donné par p
2πf (x) 1 + (f ′ (x))2 dx.
L’aire totale de la surface est donc :
Z L p
A[f ] = 2π f (x) 1 + (f ′ (x))2 dx.
0

3.3 Équation différentielle pour minimiser A[f ] :

Pour minimiser cette fonctionnelle, nous devons utiliser l’équation d’Euler-Lagrange :
 
d ∂L ∂L
′
− = 0,
dx ∂f ∂f
où p
L = f (x) 1 + (f ′ (x))2 .

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 7

 Calculons les dérivées nécessaires :
∂L p
= 1 + (f ′ (x))2 ,
∂f
et
∂L f (x)f ′ (x)
= .
∂f ′
p
1 + (f ′ (x))2
Ensuite, nous devons dériver par rapport à x :
!
′
 
d ∂L d f (x)f (x)
′
= p .
dx ∂f dx 1 + (f ′ (x))2
En utilisant l’équation d’Euler-Lagrange, nous obtenons :
!
d f (x)f ′ (x) p
p − 1 + (f ′ (x))2 = 0.
dx 1 + (f ′ (x))2

Cette équation différentielle doit être résolue pour trouver la fonction f (x) qui minimise l’aire
A[f ].

x− L
 
3.4 Vérifions si f (x) = c0 cosh c0
2
est une solution lorsque r0 = rL :

Calculons f ′ (x) :
    
′ d x − L/2 x − L/2
f (x) = c0 cosh = sinh .
dx c0 c0
Calculons f (x)f ′ (x) :
   
′ x − L/2 x − L/2
f (x)f (x) = c0 cosh sinh .
c0 c0
Nous avons aussi :  
′ 2 2 x − L/2
(f (x)) = sinh .
c0
∂L
Calculons ∂f ′
:
   
x−L/2 x−L/2
∂L f (x)f (x) c ′
0 cosh c0
sinh c0
′
=p = r .
∂f 1 + (f ′ (x))2

2 x−L/2

1 + sinh c0

Utilisant l’identité hyperbolique cosh2 (y) − sinh2 (y) = 1, nous avons :

   
2 x − L/2 2 x − L/2
1 + sinh = cosh .
c0 c0
Donc,  
 
x−L/2
∂L c0 cosh c0
sinh x−L/2
c0

x − L/2

=   = c0 sinh .
∂f ′ cosh x−L/2 c0
c0

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 8

 Calculons maintenant la dérivée par rapport à x :
    
d x − L/2 x − L/2
c0 sinh = cosh .
dx c0 c0
Substituons dans l’équation différentielle :
  s      
x − L/2 2 x − L/2 x − L/2 x − L/2
cosh − 1 + sinh = cosh − cosh = 0.
c0 c0 c0 c0
 
Cela montre que f (x) = c0 cosh x−L/2 c0
satisfait l’équation différentielle.

Pour trouver l’expression de la constante c0 , utilisons les conditions aux limites f (0) = r0 et
f (L) = rL . Si r0 = rL , alors :
 
−L/2
r0 = c0 cosh
c0
et  
L/2
rL = c0 cosh .
c0
Comme r0 = rL , nous avons :
 
L/2
r0 = rL = c0 cosh .
c0
Cela donne l’expression pour c0 :
r0
c0 =  .
L/2
cosh c0

Ainsi, c0 peut être déterminé implicitement à partir de cette équation.

3.5 Vérifions que dans le cas de la surface considérée ici la courbure moyenne est
nulle :
La courbure moyenne Hm est définie comme :
 
1 1 1
Hm = +
2 R1 R2
Pour une surface de révolution r = f (x), les courbures principales R1 et R2 sont données par :

3/2 ′′
R1 = 1 + (f ′ )2 f
p
R2 = f 1 + (f ′ )2
Ainsi, la courbure moyenne Hm est :
p !
1 f ′′ (1 + (f ′ )2 )3/2 f 1 + (f ′ )2
Hm = +
2 f′ 1
 
x−L/2
Pour la solution f (x) = c0 cosh c0
, nous avons :

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 9

  
′ x − L/2
f (x) = sinh
c0
 
′′ 1 x − L/2
f (x) = cosh
c0 c0
Substituons ces expressions dans la formule de Hm :
    
2 x−L/2 x−L/2
1 c 0 cosh c0
+ c 0 cosh c0
Hm =    
2 c cosh2 x−L/2
0 c0

En simplifiant, nous obtenons :

     
x−L/2 x−L/2
1 c0 cosh c0
cosh c0
+1
Hm =    
2 c cosh 2 x−L/2
0 c0

Puisque cosh(x) est toujours positif et non nul, la moyenne des courbures est effectivement nulle
pour cette forme de f (x), prouvant que la courbure moyenne est nulle.

3.6 Proposition d un algorithme, permettant de trouver (numériquement) la forme

de l’interface minimisant la surface :
L’algorithme de minimisation de la surface pour une surface de révolution vise à trouver la
forme de l’interface qui minimise l’aire de la surface. Il peut être décrit en plusieurs étapes clés :
1. Initialisation :
— Choisir une forme initiale pour l’interface de la surface de révolution.
— Discrétiser l’intervalle de la forme de l’interface en un certain nombre de points.
— Initialiser les paramètres de l’algorithme tels que le pas de descente de gradient, le critère
d’arrêt et le nombre maximal d’itérations.
2. Calcul de l’Aire de la Surface :
— Utiliser la forme de l’interface actuelle pour calculer la fonctionnelle de l’aire de la surface.
3. Calcul du Gradient :
— Calculer le gradient de la fonctionnelle de l’aire par rapport à la forme de l’interface.
4. Mise à Jour de la Forme de l’Interface :
— Utiliser le gradient calculé pour mettre à jour la forme de l’interface en utilisant la méthode
de descente de gradient.
— Appliquer une limite pour éviter des valeurs extrêmes ou non physiques de la forme de
l’interface.
5. Vérification de Convergence :
— Vérifier si la différence entre les valeurs successives de la fonctionnelle de l’aire est
inférieure à un certain critère d’arrêt. Si oui, terminer l’algorithme.
6. Répéter :
— Répéter les étapes précédentes jusqu’à ce que le critère d’arrêt soit satisfait ou que le
nombre maximal d’itérations soit atteint.
7. Résultat :
— La forme finale de l’interface après convergence représente la forme de l’interface minimi-
sant la surface.

ENSEEIHT MF2E 1A Introduction aux méthodes d’optimisation 10