Exercice

La distillation sous vide du résidu atmosphérique permet d’obtenir des distillats visqueux qui
peuvent être utilisés comme matière première dans la fabrication des huiles lubrifiantes (Figure
1). La pression de vapeur au sein de la tour est donnée par l’expression suivante :
0,75 1
𝑓 (𝑃) = 0,65 − 1+𝑃2 − 0,65𝑃 𝑎𝑟𝑐𝑡𝑎𝑛𝑔 (𝑃)

En considérant une précision de 10-4, minimiser la fonction objectif f par

- La méthode de la section dorée La méthode de la section dorée en considérant ]0, 3[

puis [0, 3];
- La méthode de Newton avec Po= 0,1.

Comparez entre les deux méthodes.

Figue 1 : Colonne de distillation sous vide.

Corrigé

La distillation sous vide permet de récupérer en tête de la colonne du gasoil lourd puis
latéralement on peut soutirer les huiles de différentes viscosités et par conséquent, les
lubrifiants. En pratique, la pression doit être positive, en cas de dysfonctionnement la
pression affichée peut être négative.
Fonction à minimiser : càd on cherche le minimum

Fonction objectif

f(P) = 0.65-(0.75/(1+P2))-(0.65 P arctan(1/P)) dans l’intervalle ]0, 3[

1. Méthode de la section Dorée

𝒙𝒍 = 𝟎 et 𝒙𝒖 = 𝟑

𝒇 (𝒙𝒍 ) = 𝒇(𝟎) = -0.1

𝒇 (𝒙𝒖 ) = 𝒇(𝟑) = -0.0524

Calculons le ratio d’or :

√𝟓−𝟏
𝒅= (𝟑 − 𝟎) = 1.8541
𝟐

𝒙𝟏 = 𝒙𝒍 +d = 0+ 1,8541 = 1.8541

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 𝟑 − 𝟏, 𝟖𝟓𝟒𝟏 = 𝟏. 𝟏𝟒𝟓𝟗

Calculons les fonctions f(x1) et f(x2) :

𝒙𝒍 = 𝟎 ; 𝒇 (𝒙𝒍 ) = -0.1

𝒙𝒖 = 𝟑 ; 𝒇 (𝒙𝒖 ) = -0.0524

𝒙𝟏 = 𝟏. 𝟖𝟓𝟒𝟏 ; f(𝒙𝟏 ) = -0.1151

𝒙𝟐 = 𝟏. 𝟏𝟒𝟓𝟗 ; f(𝒙𝟐 ) = -0.2087

Comme 𝒇 (𝒙𝟐 )  𝒇 (𝒙𝟏 )

𝑥2 = 1,1459
𝑥1 = 1,8541
𝒙𝒖 = 𝟑
𝒙𝒍
-0,1

-0,1151 -0.0524

-0.2087
Le nouvel intervalle : 𝒙𝒍 = 0 et 𝒙𝟏 = 𝒙𝒖 = 𝟏. 𝟖𝟓𝟒𝟏

𝒇 (𝒙𝒍 ) = 𝒇(𝟎) = 𝟎

𝒇 (𝒙𝒖 ) = 𝒇(𝟏. 𝟖𝟓𝟒𝟏) = −𝟎. 𝟏𝟏𝟓𝟏

D’où :

√𝟓−𝟏
𝒅= (𝟏. 𝟖𝟓𝟒𝟏 − 𝟎) = 1.1459
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 1.8541 − 1.1459 = 𝟎. 𝟕𝟎𝟖𝟐

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 𝟎 + 1.1459 = 𝟏, 𝟏𝟒𝟓𝟗

𝒇 (𝒙𝟐 ) = -0.2889

√𝟓−𝟏
𝒅= (𝟏. 𝟖𝟓𝟒𝟏 − 𝟎. 𝟕𝟎𝟖𝟐) = 0.7082
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 1.8541 − 0.7082 = 𝟏. 𝟏𝟒𝟓𝟗

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 + 0.7082 = 𝟏. 𝟒𝟏𝟔𝟒

𝒇 (𝒙𝟏 ) = -0.1655

4ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 = 𝟏. 𝟒𝟏𝟔𝟒

√𝟓−𝟏
𝒅= (𝟏. 𝟒𝟏𝟔𝟒 − 𝟎. 𝟕𝟎𝟖𝟐) = 0.4377
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 1.4164 − 0.4377 =0.9787

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 + 0.4377 = 𝟏. 𝟏𝟒𝟓𝟗

5ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 = 1.1459

√𝟓−𝟏
𝒅= (𝟏. 𝟏𝟒𝟓𝟗 − 𝟎. 𝟕𝟎𝟖𝟐) = 0.2705
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 1.1459 − 0.2705 = 𝟎. 𝟖𝟕𝟓𝟒

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 + 0.2705 = 𝟎. 𝟗𝟕𝟖𝟕

6ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.9787
 √𝟓−𝟏
𝒅= (𝟎. 𝟗𝟕𝟖𝟕 − 𝟎. 𝟕𝟎𝟖𝟐) = 0.1672
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 𝟎. 𝟗𝟕𝟖𝟕 −0.1672= 0.8115

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.1672 = 0.8754

7ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.8754

√𝟓−𝟏
𝒅= (0.8754 − 𝟎. 𝟕𝟎𝟖𝟐) =0.1033
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.8754 −0.1033= 0.7721

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.1033= 0.8115

8ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.8115

√𝟓−𝟏
𝒅= (0.8115 − 𝟎. 𝟕𝟎𝟖𝟐) = 0.0638
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.8115 −0.0638=0.7477

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0638= 0.7720

9ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.7721

√𝟓−𝟏
𝒅= (0.7721 − 0.7082) = 0.0395
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.7721 −0.0395= 0.7326

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0395= 0.7477

11ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.7326

√𝟓−𝟏
𝒅= (0.7326 − 0.7082) = 0.0151
𝟐
𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.7326 −0.0151= 0.7175

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0151 = 0.7233

12ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.7233

√𝟓−𝟏
𝒅= (0.7233 − 0.7082) = 0.0093
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.7233 −0.0093= 0.7140

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0093= 0.7175

13ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.7175

√𝟓−𝟏
𝒅= (0.7175 − 0.7082) = 0.0057
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.7175 −0.0057= 0.7118

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0057= 0.7139

14ème itération

𝒙𝒍 = 0.7082 et 𝒙𝒖 =0.7118

√𝟓−𝟏
𝒅= (0.7118 − 0.7082) = 0.0022
𝟐

𝒙𝟐 = 𝒙𝒖 − 𝒅 = 0.7118 −0.0022= 0.7096

𝒙𝟏 = 𝒙𝒍 + 𝒅 = 0.7082 +0.0022= 0.7104

Les résultats des différentes itérations sont donnés dans le tableau suivant :

i 𝒙𝒍 𝒇(𝒙𝒍 ) 𝒙𝟐 𝒇(𝒙𝟐 ) 𝒙𝟏 𝒇(𝒙𝟏 ) 𝒙𝒖 𝒇(𝒙𝒖 ) d

1 0 -0.1 1.1459 -0.2087 1.8541 -0.1151 3 -0.0524 1.8541

2 0 -0.1 0.7082 -0.2889 1.1459 -0.2087 1.8541 -0.1151 1.1459

3 0.7082 -0.2889 1.1459 -0.2087 1.4164 -0.1655 1.8541 -0.1151 0.7082

4 0.7082 -0.2889 0.9787 -0.2396 1.1459 -0.2087 1.4164 -0.1655 0.4377

5 0.7082 -0.2889 0.8754 -0.2593 0.9787 -0.2396 1.1459 -0.2087 0,2705

6 0.7082 -0.2889 0.8115 -0.2712 0.8754 -0.2593 0.9787 -0.2396 0.1672

7 0.7082 -0.2889 0.7721 -0.2782 0.8115 -0.2712 0.8754 -0.2593 0.1033

8 0.7082 -0.2889 0.7477 -0.2824 0.7721 -0.2782 0.8115 -0.2712 0.0638

9 0.7082 -0.2889 0.7326 -0.2850 0.7477 -0.2824 0.7721 -0.2782 0.0395

10 0.7082 -0.2889 0.7233 -0.2865 0.7326 -0.2850 0.7477 -0.2824 0.0244

11 0.7082 -0.2889 0.7175 -0.2874 0.7233 -0.2865 0.7326 -0.2850 0.0151

12 0,7082 -0.2889 0.7140 -0.2880 0.7175 -0.2874 0.7233 -0.2865 0.0093

13 0.7082 -0.2889 0.7118 -0.2883 0.7140 -0.2880 0.7175 -0.2874 0.0057

14 0.7082 -0.2889 0.7104 -0.2886 0.7118 -0.2883 0.7140 -0.2880 0.0036

15 0.7082 -0.2889 0.7096 -0.2887 0.7104 -0.2886 0.7118 -0.2883 0.0022

16 0.7082 -0.2889 0.7090 -0.2888 0.7096 -0.2887 0.7104 -0.2886 0.0014

17 0.7082 -0.2889 0.7088 -0.2888 0.7090 -0.2888 0.7096 -0.2887 0.0008

La convergence est atteinte après la 17ème itération et le min est atteint à - 0,2889
dans l’intervalle de p [0.7082 , 0.7088 ]

 (P) = 0.65-(0.75/(1+P2))-(0.65 P arctan(1/P)) dans l’intervalle [0, 3]

 Les résultats des différentes itérations sont donnés dans le tableau suivant :

i 𝒙𝒍 𝒇(𝒙𝒍 ) 𝒙𝟐 𝒇(𝒙𝟐 ) 𝒙𝟏 𝒇(𝒙𝟏 ) 𝒙𝒖 𝒇(𝒙𝒖 ) d

1 0 -0.1 1.1459 -0.2087 1.8541 -0.1151 3 -0.0524 1.8541

2 0 -0.1 0.7082 -0.2889 1.1459 -0.2087 1.8541 -0.1151 1.1459

3 0.7082 -0.2889 1.1459 -0.2087 1.4164 -0.1655 1.8541 -0.1151 0.7082

 4 0.7082 -0.2889 0.9787 -0.2396 1.1459 -0.2087 1.4164 -0.1655 0.4377

5 0.7082 -0.2889 0.8754 -0.2593 0.9787 -0.2396 1.1459 -0.2087 0,2705

La convergence est atteinte après la 5ème itération et le min est atteint à - 0,2889 dans
l’intervalle de p [0.7082 , 1,1459 ]

2) Méthode de Newton
𝒇′ (𝒙𝒊)
Xi+1 = Xi -
𝒇′′(𝒙𝒊)

f(X) = 0.65-(0.75/(1+X^2))-(0.65 X arctan(1/X))

𝟏,𝟓 𝑿 𝟎,𝟔𝟓 𝑿
f’(X) = -0.65 X arctan(1/X) +
(𝟏+𝑿𝟐 )𝟐 𝟏+𝑿𝟐

𝟏,𝟑 𝟏,𝟓−𝟏,𝟑 𝑿𝟐 𝟔𝑿
f’’(X) = 𝟐
+ -
𝟏+𝑿 (𝟏+𝑿𝟐 )𝟐 (𝟏+𝑿𝟐 )𝟑
La convergence est atteinte à la 18ème itération la plus petite valeur atteinte est de -0,31 par la
suite les valeurs augmentent pas nécessaire d’aller au-delà.

Comparaison : Précision à 10-4

Méthodes Minimum Nombre d’itérations

Section Dorée - 0,2889 17

Newton -0,31 18

Conclusion

Théoriquement, la méthode de la section dorée est plus rapide. Cependant,

pour le cas étudié et en utilisant la méthode de la section Dorée avec ]0, 3[, la
pression minimale est atteinte après 17 itérations à la valeur de -0.2889 alors
que la plus basse valeur a été obtenue après 18 itérations via la méthode de
Newton, pour le cas étudié les deux méthodes ne sont pas rapides et se valent.

Par contre, en considérant l’intervalle [0, 3] la méthode de la section est plus

rapide on atteint la convergence après 5 itérations.

Exercice
Un fluide non-Newtonien dont la circulation du débit dans une conduite cylindrique est donnée
par l’expression suivante :
Q (r, x) = r3-3r + r x2 (x est le rapport de la contrainte de l’écoulement sur la
contrainte de cisaillement)
- Déterminer les extremums de la fonction objectif (débit).

Corrigé
Fonction objectif : Q (r, x) = r3-3r + r x2

𝝏𝑸
(𝒓, 𝒙) 𝟑 𝒓𝟐 − 𝟑 + 𝒙𝟐
𝝏𝒓
𝜵 𝑸 (𝒓, 𝒙) = (𝝏𝑸 ) = ( )
(𝒓, 𝒙) 𝟐𝒓𝒙
𝝏𝒙

𝝏𝑸
(𝒓, 𝒙) = 𝟎 𝟑 𝒓𝟐 − 𝟑 + 𝒙𝟐 = 𝟎 (𝟏)
𝝏𝒓
{𝝏𝑸 ce qui donne {
(𝒓, 𝒙) = 𝟎 𝟐𝒓𝒙 = 𝟎 (𝟐)
𝝏𝒙

(2) donne r=0 ou x= 0

(1) donne 𝒙𝟐 − 𝟑 = 𝟎 ou 𝟑 𝒓𝟐 − 𝟑 = 𝟎
Les systèmes de solutions :

𝑟=0
{ ou { 𝑥 = 0
𝑥 = ∓√3 𝑟 = ∓1

On aura 4 points stationnaires

𝑃1 (0, −√3) ; 𝑃2 (0, √3) ; 𝑃3 (−1, 0) ; 𝑃4 (1, 0)

Nature des points stationnaires

Calculons la Hessienne

𝝏𝟐 𝑸 𝝏𝟐 𝑸
𝝏𝒓𝟐 𝝏𝒓𝝏𝒙 𝟔𝒓 𝟐𝒙
H =( ) =( )
𝝏𝟐 𝑸 𝝏𝟐 𝑸 𝟐𝒙 𝟐𝒓
𝝏𝒙𝝏𝒓 𝝏𝒙𝟐

La matrice hessienne est donnée par :

𝟔𝒓 𝟐𝒙
∆=( )
𝟐𝒙 𝟐𝒓
∆1 = 6𝑟

𝟔𝒓 𝟐𝒙
∆2 = ( ) = 12r2- 4x2
𝟐𝒙 𝟐𝒓
D’où on aura les possibilités suivantes :

Points 1 2 Nature de  Nature de x Q(r, x)

(0, −√3) 0 -12 - Pas 0
d’extremum
(0, √3) 0 -12 - Pas 0
d’extremum
(−1, 0) -6 12 0 Maximum 2
(1, 0)) 6 12 0 Minimum -2