Écoulements incompressibles et non

visqueux - Partie 1

MEC 757 - Introduction à l’Aérodynamique

Hiver 2025
Séance 3

21 Janvier 2025
Organisation de la séance

Équations de Navier-Stokes

Vorticité et circulation

Équations d’Euler

Équations de Bernoulli

Théorème de Kelvin

Écoulements potentiels

Conclusions

Annexes

MEC757 - Hiver 2025 - Séance 3 2

 Équations de Navier-Stokes

MEC757 - Hiver 2025 - Séance 3 3

Rappel : Bilan de la quantité de mouvement

ˆ #» 1 #» #»2 ÿ #»
˚ ‹
#»
(1) flV dv + flV V · n dS = F extævc
ˆt V vc A vc

I Les forces externes qui s’appliquent sur le volume de contrôle sont :

I Forces volumiques :
I Force d’accélération gravitationnelle (poids)
I Forces de surface :
I Force résultante de pression du fluide externe
I Force de cisaillement visqueux du fluide externe
I Force de réaction aux points de contact

ÿ #» #» #» #» ÿ #»
(2) F vc = F grav + F pres + F visq + F reac
A vc
ÿ #» ˚ ‹
#» ÿ #»
(3) F extævc = fl #»
g dv ≠ p n dS + F visq +
#» F reac
Vvc A vc A vc

MEC757 - Hiver 2025 - Séance 3 4

Tenseur de cisaillement visqueux
I Pour un fluide Newtonien, le cisaillement visqueux est proportionnel au
tenseur de déformation du champ de vitesse
3 4
ˆu ˆv
(4) ·xy = ·yx = µ +
ˆy ˆx
1 2
ˆu ˆw
(5) ·xz = ·zx =µ +
ˆz ˆx
3 4
ˆv ˆw
(6) ·yz = ·zy =µ +
ˆz ˆy

I Si les gradients de vitesse sont importants,

la viscosité induit une force normale en
I ⁄ coefficient de viscosité
plus de la pression
de volume
! #» #»" ˆu
(7) ·xx = ⁄ Ò · V + 2µ I Sous l’hypothèse de Stokes
ˆx
! #» #»" ˆv
(8) ·yy = ⁄ Ò · V + 2µ 2
ˆy (10) ⁄=≠ µ
! #» #»" 3
ˆw
(9) ·zz = ⁄ Ò · V + 2µ
ˆz
MEC757 - Hiver 2025 - Séance 3 5
Forces visqueuses
I Les forces visqueuses s’expriment :
(11) S T
‹
! " #» ·xx ·xy ·xz
#»
F visq = #»
· · n · dS avec : · = U·yx ·yy ·yz V .
A vc ·zx ·zy ·zz

MEC757 - Hiver 2025 - Séance 3 6

Équations de Navier-Stokes
I Ainsi, � (1), (3) et (11), en négligeant les forces volumiques (poids
de l’air négligeables) et en considérant le fluide isolé (pas de force de
réaction), on peut écrire :
ˆ #» 1 #» #»2 ! #»" #»
˚ ‹ ‹ ‹
#»
flV dv + flV V · n dS = ≠ p n dS +
#» · · n · dS
ˆt V vc A vc A vc A vc

I En appliquant le théorème de la divergence et du gradient :

ˆflu #» 1 #»2 ˆp ˆ·xx ˆ·yx ˆ·zx

(12) + Ò · fluV = ≠ + + +
ˆt ˆx ˆx ˆy ˆz
ˆflv #» 1 #»2 ˆp ˆ·xy ˆ·yy ˆ·zy
(13) + Ò · flv V = ≠ + + +
ˆt ˆy ˆx ˆy ˆz
ˆflw #» 1 #»2 ˆp ˆ·xz ˆ·yz ˆ·zz
(14) + Ò · flw V = ≠ + + +
ˆt ˆz ˆx ˆy ˆz

MEC757 - Hiver 2025 - Séance 3 7

Équations de Navier-Stokes en incompressible

ˆu ˆv ˆw
(15) + + =0
ˆx ˆy ˆz
3 2 2 2
4
ˆu ˆu ˆu ˆu 1 ˆp ˆ u ˆ u ˆ u
(16) +u +v +w =≠ +‹ + 2 + 2
ˆt ˆx ˆy ˆz fl ˆx ˆx2 ˆy ˆz
3 2 2 2
4
ˆv ˆv ˆv ˆv 1 ˆp ˆ v ˆ v ˆ v
(17) +u +v +w =≠ +‹ 2
+ 2+ 2
ˆt ˆx ˆy ˆz fl ˆy ˆx ˆy ˆz
3 2 2 2
4
ˆw ˆw ˆw ˆw 1 ˆp ˆ w ˆ w ˆ w
(18) +u +v +w =≠ +‹ + +
ˆt ˆx ˆy ˆz fl ˆz ˆx2 ˆy 2 ˆz 2

MEC757 - Hiver 2025 - Séance 3 8

 Vorticité et circulation

MEC757 - Hiver 2025 - Séance 3 9

Déformation d’un élément fluide

I Lors de son déplacement un volume fluide se déforme et tourne.

#»
I Le gradient de vitesse est noté ÒV et se décompose en une partie
symétrique S et une partie anti-symétrique :

S T
ˆu ˆu ˆu
ˆx ˆy ˆz
#» W ˆv ˆv X
(19) ÒV = U ˆx
ˆv
ˆy ˆz V
ˆw ˆw ˆw
ˆx ˆy ˆz
#»
(20) ÒV = S + .

MEC757 - Hiver 2025 - Séance 3 10

Déformation d’un élément fluide
I Le tenseur des vitesses de déformation S conduit à une modification
de la forme : allongement et cisaillement.
I Le tenseur des vitesses de rotation conduit à un changement de
son orientation.

S 1 2 ! " T
1 1 ˆu
ˆu
ˆx 2
ˆu
ˆy + ˆv
ˆx 2 1 ˆz + ˆw
ˆx 2X
W 1 2
W 1 ˆu ˆv 1 X
(21) S = W 2 ˆy + ˆx ˆv
ˆy 2
ˆv
ˆz + ˆw
ˆy X
U ! " 1 2 V
1 ˆu 1
2 ˆz + ˆx +
ˆw ˆv ˆw ˆw
2 ˆz ˆy ˆz

S 1 2 ! "T
1 1 ˆu
0 2
ˆu
ˆy ≠ ˆv
ˆx 2 1 ˆz ≠ ˆw
ˆx 2X
W 1 2
W 1 ˆv 1 X
(22) = W 2 ˆx ≠ ˆu
ˆy 0 2
ˆv
ˆz ≠ ˆw
ˆy X
U ! " 1 2 V
1 ˆw 1
ˆu
2 ˆx ≠ ˆz 2
ˆw
ˆy ≠ ˆv
ˆz 0

MEC757 - Hiver 2025 - Séance 3 11

Écoulement visqueux et vorticité
I Le tenseur des contraintes de viscosité est lié aux déformations de
vitesse :
2 1 #» #»2
(23) · = 2µ S ≠ µ Ò · V I comportement de Newton-Stokes
3
I La vorticité (rotationnel de vitesse) est créée par le cisaillement
visqueux.
# » #» # » #» #» #» #» #»
(24) rot V = curl V = Ò ◊ V = 2 Ê #» #» ◊ A
Ê = ·A
I Le cisaillement crée de la vorticité. Les effets visqueux sont des
mécanismes qui génèrent de la vorticité � rotationnel.

MEC757 - Hiver 2025 - Séance 3 12

Expressions du rotationnel
I Attention aux calculs des composantes de l’opérateur rotationnel dans
les différents systèmes de coordonnées �
I En cartésien
3 4 1 2
# » #» ˆfz ˆfy ˆfx ˆfz
(25) curl F = ≠ ex +
#» ≠ #»
ey
ˆy ˆz ˆz ˆx
3 4
ˆfy ˆfx
+ ≠ #»
ez
ˆx ˆy
I En cylindrique
1 2 1 2
# » #» 1 ˆfz ˆf◊ #» ˆfr ˆfz #»
(26) curl F = ≠ er + ≠ e◊
r ˆ◊ ˆz ˆz ˆr
3 4
1 ˆ(rf◊ ) ˆfr #»
+ ≠ ez
r ˆr ˆ◊
I En sphérique
3 4
# » #» 1 ˆ(f„ sin ◊) ˆf◊
(27) curl F = ≠ #»
er
r sin ◊ ˆ◊ ˆ„
3 4 3 4
1 ˆfr 1 ˆ(rf„ ) 1 ˆ(rf◊ ) ˆfr
+ ≠ e◊ +
#» ≠ #»
e„
r sin ◊ ˆ„ r ˆr r ˆr ˆ◊
MEC757 - Hiver 2025 - Séance 3 13
Théorème de Stokes
I Essentiel en aérodynamique : la circulation est un outil qui permet de
calculer la portance aérodynamique.
I Circulation : Opposé (par convention) de l’intégrale linéique de la
vitesse sur un contour fermé C
I Par application du théorème de Stokes, pour toute surface S ouverte
s’appuyant sur le contour C :
˛ ¨
#» #» #» #» # »
(28) =≠ V · dl = ≠ Ò ◊ V · dS
C S

� La circulation de vitesse sur C correspond

au flux de vorticité au travers d’une surface
qui s’appuie sur C.
I En rapetissant C à une taille infinitésimale,
on peut écrire d la circulation par unité de
surface selon :
#» #»
(29) d = ≠Ò ◊ V · #»
n dS .
MEC757 - Hiver 2025 - Séance 3 14
 Équations d’Euler

MEC757 - Hiver 2025 - Séance 3 15

Retour sur les équations de Navier-Stokes
I Forme algébrique, en cartésien, pour un écoulement incompressible

ˆu ˆv ˆw
(30) + + =0
ˆx ˆy ˆz
3 2 2 2
4
ˆu ˆu ˆu ˆu 1 ˆp ˆ u ˆ u ˆ u
(31) +u +v +w =≠ +‹ + 2 + 2
ˆt ˆx ˆy ˆz fl ˆx ˆx2 ˆy ˆz
3 2 2 2
4
ˆv ˆv ˆv ˆv 1 ˆp ˆ v ˆ v ˆ v
(32) +u +v +w =≠ +‹ + +
ˆt ˆx ˆy ˆz fl ˆy ˆx2 ˆy 2 ˆz 2
3 2 2 2
4
ˆw ˆw ˆw ˆw 1 ˆp ˆ w ˆ w ˆ w
(33) +u +v +w =≠ +‹ + +
ˆt ˆx ˆy ˆz fl ˆz ˆx2 ˆy 2 ˆz 2

I Forme vectorielle, en cartésien, pour un écoulement incompressible

#» #»
(34) Ò·V =0
#»
DV 1 #» #»
(35) = ≠ Òp + ‹ V
Dt fl

MEC757 - Hiver 2025 - Séance 3 16

Retour sur la dérivée particulaire de la vitesse
I On peut écrire la dérivée particulaire d’une grandeur scalaire selon :

Df ˆf #» #»
(36) = + V · Òf
Dt ˆt
I Pour l’accélération particulaire, on voit écrit :
#» #»
DV ˆV #» #»
(37) = +ÒV · V
Dt ˆt
� expression valide uniquement en coordonnées cartésiennes �
I On utilisera l’expression suivante valide pour tous les systèmes de
coordonnées
3 24
#» 1 #» #»2
#» #»
DV ˆV #» V
(38) = +Ò ≠V ◊ Ò◊V
Dt ˆt 2

MEC757 - Hiver 2025 - Séance 3 17

Équations d’Euler

Hypothèses
H 1 Incompressible
H 2 Poids négligeable
H 3 non visqueux

#» #»
(39) Ò·V =0 Conservation de la masse � continuité

Conservation de la quantité de mouvement

3 24
#» 1 #» #»2
#» #»
DV ˆV #» V 1 #»
(40) = + Ò ≠ V ◊ Ò ◊ V = ≠ Òp
Dt ˆt 2 fl

� Les efforts de pression équilibrent

les variations temporelles et convectives de vitesse

MEC757 - Hiver 2025 - Séance 3 18

 Équations de Bernoulli

MEC757 - Hiver 2025 - Séance 3 19

Équation de Bernoulli
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Sur une ligne de courant

On note t le vecteur tangent à la ligne de courant

#»
#» #»
(41) V · t =V
3 24
#» 1 #» #»2
#» #»
DV V◆ #» V
ˆ◆
(42) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 24 1 #» #»2
#» V #» #» 1 #» #»
(43) · t ≠ V ◊ Ò ◊ V · t = ≠ Òp · t
#»
Ò
2 fl

MEC757 - Hiver 2025 - Séance 3 20

Équation de Bernoulli
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Sur une ligne de courant

On note t le vecteur tangent à la ligne de courant

#»
#» #»
(41) V · t =V
3 24
#» 1 #» #»2
#» #»
DV V◆ #» V
ˆ◆
(42) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 24 2⇠⇠
V #» 1 #» ⇠#» ⇠ dp
(43) ⇠◊ V · t = ≠
#»
d ≠V⇠ ◊ ⇠Ò
2 ⇠ fl

MEC757 - Hiver 2025 - Séance 3 20

Équation de Bernoulli
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Sur une ligne de courant

#» #»
(41) V · t =V
3 24
#» 1 #» #»2
#» #»
DV ˆ◆V◆ #» V
(42) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 4
1
(43) d flV 2 + p = 0
2

MEC757 - Hiver 2025 - Séance 3 20

Équation de Bernoulli
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Sur une ligne de courant

#» #»
(41) V · t =V
3 24
#» 1 #» #»2
#» #»
DV ˆ◆V◆ #» V
(42) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 4
1
(43) d flV 2 + p = 0
2

1
(44) fl V 2 + p = cst sur une ligne de courant
2
MEC757 - Hiver 2025 - Séance 3 20
Équation de Bernoulli - forme forte
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Irrotationnel
#» #»
(45) V · t =V
3 24
#» 1 #» #»2
#» #»
DV ˆ◆V◆ #» V
(46) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 24
#» V 1 #»
(47) Ò = ≠ Òp
2 fl

MEC757 - Hiver 2025 - Séance 3 21

Équation de Bernoulli - forme forte
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Irrotationnel
#» #»
(45) V · t =V
3 24
#» 1 #» #»2
#» #»
DV ˆ◆V◆ #» V
(46) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 4
#» 1
(47) Ò flV 2 + p = 0
2

MEC757 - Hiver 2025 - Séance 3 21

Équation de Bernoulli - forme forte
Hypothèses
H1 Incompressible
H2 Poids négligeable
H3 Non visqueux
H4 Régime permanent
H5 Irrotationnel
#» #»
(45) V · t =V
3 24
#» 1 #» #»2
#» #»
DV ˆ◆V◆ #» V
(46) = +Ò ≠V ◊ Ò◊V
Dt ◆ˆt 2
3 4
#» 1
(47) Ò flV 2 + p = 0
2

1
(48) fl V 2 + p = cst partout dans l’écoulement
2

MEC757 - Hiver 2025 - Séance 3 21

Écoulement Quasi-1D

I En régime permanent, pour un fluide incompressible

I Conservation de la masse

(49) A 1 V1 = A 2 V 2

I Pour un fluide non visqueux, irrotationnel � (48)

ı̂ 2 (p ≠ p )
ı 1 2
(50) V1 = ı 51 22 6
Ù
fl AA 12 ≠ 1

MEC757 - Hiver 2025 - Séance 3 22

Applications : soufflerie aérodynamique

I Soufflerie veine ouverte

MEC757 - Hiver 2025 - Séance 3 23

Applications : soufflerie aérodynamique
I Soufflerie veine ouverte NRC (Conseil national de recherches du
Canada)

MEC757 - Hiver 2025 - Séance 3 24

Applications : soufflerie aérodynamique
I Soufflerie veine fermée

MEC757 - Hiver 2025 - Séance 3 25

Applications : soufflerie aérodynamique
I Soufflerie veine fermée NRC

MEC757 - Hiver 2025 - Séance 3 26

Prises de pression statique et totale

I Si l’on suppose l’écoulement de fluide parfait, incompressible en

régime permanent et que l’on néglige les forces extérieures, on peut
appliquer la relation de Bernoulli sur une ligne de courant.
I Orifice perpendiculaire à l’écoulement : prise de pression statique

(51) pA = p Œ
I Orifice parallèle à l’écoulement : prise de pression totale
1 2
(52) p B = p Œ + flŒ V Œ
2
I On peut déterminer la vitesse VŒ avec les deux prises de pression
Û
2 (pB ≠ pA )
(53) VΠ=
fl

MEC757 - Hiver 2025 - Séance 3 27

Tube de Prandtl

Û
2 (pB ≠ pA )
(53) VΠ=
fl

I Pour connaı̂tre la vitesse VŒ , il faut déterminer pŒ , c’est le rôle des

orifices latéraux sur un tube de Prandtl.
I Les positions des prises de pression doivent s’affranchir de biais de
mesures

MEC757 - Hiver 2025 - Séance 3 28

Tube de Pitot et de Prandtl - cas compressible

I Si l’on suppose l’écoulement de fluide parfait, compressible,

isentropique flp“ = cst en régime permanent la relation de Bernoulli
devient
V2 “ p
(54) + = cst
2 “≠1fl
I On peut déterminer la vitesse avec deux prises de pression

(55)
Û 3 4
“≠1 pB pA Nécessite la connaissance des températures
VŒ = 2 ≠
“ flB flA d’arrêt : calibration de la probe

MEC757 - Hiver 2025 - Séance 3 29

Applications : Mesure de la vitesse sur avion
I Avion de chasse F-86H (1947) : Capteur immergé dans l’écoulement
sain, loin de toute perturbation.

MEC757 - Hiver 2025 - Séance 3 30

Applications : Mesure de la vitesse sur avion

I Pression statique et totale

prise indépendamment.
I Positionnement des prises de
pression où Cp ƒ 0
I Sensibilité à l’alignement
avec l’écoulement :
importance du design
MEC757 - Hiver 2025 - Séance 3 31
Applications : Mesure de la vitesse sur avion

et Formule 1 ...

I Pression statique et totale

prise indépendamment.
I Positionnement des prises de
pression où Cp ƒ 0
I Sensibilité à l’alignement
avec l’écoulement :
importance du design
MEC757 - Hiver 2025 - Séance 3 31
 Théorème de Kelvin

MEC757 - Hiver 2025 - Séance 3 32

Démonstration
I Pour un écoulement incompressible, non visqueux :
� équation d’Euler pour la quantité de mouvement :
#»
DV 1 #»
(35) = ≠ Òp
Dt fl
I La circulation est définie :
˛
#» #»
(56) =≠ V · dl
C

I La circulation transportée par l’écoulement :

#» #»
D D V #» #» Ddl
˛ ˛
#» #»
Ò ◊ Òf = 0
#»
=≠ · dl ≠ V ·
Dt C Dt C Dt
3 4 3 24
1 #» V
¨ ˛
#»
= Ò◊ Òp ds ≠ d
S fl C 2
=0

MEC757 - Hiver 2025 - Séance 3 33

Théorème de Kelvin
I Pour un écoulement incompressible, non visqueux, la circulation se
conserve au cours du mouvement.
D
(57) =0
Dt

MEC757 - Hiver 2025 - Séance 3 34

 Écoulements potentiels

MEC757 - Hiver 2025 - Séance 3 35

Définition du potentiel de vitesse
I Pour un écoulement irrotationnel, on peut poser Ï telle que :
#» #» #» #»
(58) V = ÒÏ et on a bien : Ò ◊ V = 0

I Ï est la fonction du potentiel de vitesse (c’est un scalaire)

I On a l’équivalence suivante (sans restriction existence et unicité) :
écoulement potentiel … écoulement irrotationnel
I Le modèle des écoulements potentiels est très utilisé en
aérodynamique, lorsque les effets visqueux peuvent être négligés (très
grand nombre de Reynolds), en 2D, en 3D en incompressible et en
compressible. Il donne une représentation fidèle de la structure
principale de l’écoulement.

MEC757 - Hiver 2025 - Séance 3 36

Définition de la fonction courant
I Pour un écoulement incompressible 2D, il existe une fonction Â
permettant de décrire l’écoulement :
I
2 2
u = ˆÂ #» #» ˆ  ˆ Â
(59) ˆy Ò·V = ≠ =0
v = ≠ ˆx ˆÂ
ˆxˆy ˆyˆx
I Â est la fonction de courant.
I Les iso-contours de  sont les lignes de courant.
#»
I C’est l’intensité d’un champ vectoriel  = Â(x, y) e#»z qui permet de
décrire le débit volumique passant dans un tube de courant, par unité
d’envergure.
#» 1 #» 2
(60) d = Ò · d #»
x = V ◊ d #»x · e#»z

MEC757 - Hiver 2025 - Séance 3 37

Écoulements potentiels

I Sous les hypothèses d’écoulement non visqueux plan incompressible et

irrotationnel
I La vitesse dérive d’un potentiel (irrotationnel) et elle est à divergence
nulle (incompressible)
#» #» #» 1 #» 2
(61) Ò·V = divV = div ÒÏ = Ï = 0

I D’autre part, on a également en écoulement plan incompressible et

irrotationnel,
3 4
#» #» ˆv ˆu #»
(62) Ò ◊ V = ≠ ez
ˆx ˆy
3 4 3 4 3 4
ˆv ˆu ˆ ˆÂ ˆ ˆÂ
(63) ≠ = ≠ ≠ =≠ Â=0
ˆx ˆy ˆx ˆx ˆy ˆy

MEC757 - Hiver 2025 - Séance 3 38

Principe de superposition
I Pour un écoulement non visqueux plan incompressible et irrotationnel
 et Ï sont solutions de l’équation de Laplace.
(64) Â=0 et Ï=0
I C’est une équation linéaire : la somme de solutions est solution
I Principe de superposition : un écoulement compliqué peut être
décomposé en une somme de problèmes élémentaires qui respectent
les mêmes hypothèses.
(65) Â = Â1 + Â2 + · · · + Ân
(66) Ï = Ï 1 + Ï2 + · · · + Ï n
I De fait, les champs de vitesses se somment, mais pas les champs de
pression �
#» #» #» #»
(67) V = V 1 + V 2 + ··· + V n
(68) p ”= p1 + p2 + · · · + pn

� Dans le cours, on va établir les solutions pour un certain nombre de

problèmes élémentaires et s’amuser à les superposer
MEC757 - Hiver 2025 - Séance 3 39
Sources élémentaires

Ces sources seront étudiées en TD

I Écoulement uniforme
I Source / puits 2D
I Tourbillon 2D
I Doublet Distribution du coefficient
de pression pour un
NACA0012 à 9° d'incidence

MEC757 - Hiver 2025 - Séance 3 40

 Conclusions

MEC757 - Hiver 2025 - Séance 3 41

Résolution d’un problème en aérodynamique

I À grand nombre de Reynolds, on Écoulement irrotationnel

peut distinguer la zone externe et la non visqueux incompressible
couche limite pour un écoulement
attaché et la couche limite est très méthode pression
mince. itérative

I Beaucoup d’applications
Écoulement
d’aérodynamique sont très bien couche limite
visqueux
décrites par les hypothèses
d’écoulement irrotationnel et non
visqueux.
I Pour obtenir le champ de pression,
nous étudierons les écoulements
plans potentiels incompressibles

MEC757 - Hiver 2025 - Séance 3 42

 Annexes

MEC757 - Hiver 2025 - Séance 3 43

In [3]: x = [Link]( (- [Link](-3,0.5)[::-1],[Link](-3,0.5)))
integ_cste = [0.25,0.5,1.0,2.0,3.0,5.0,8.0,10.0,12.]

In [4]: for i,A0 in enumerate(integ_cste):

[Link](x,A0/x,color=f'C{i}',label=fr'$A_0={A0}$')
[Link](x,-A0/x,linestyle='--',color=f'C{i}')
[Link](-10,10)
[Link]('x')
[Link]('y')
[Link](loc='center left',bbox_to_anchor=(1.0,0.5))
[Link]('image_hyperboles.png')
[Link]()
In [5]: x = [Link](-3,3,199)
y = [Link](-10,10,201)
X,Y = [Link](x,y)
U = -X
V = Y

In [6]: [Link](figsize=(5.0,3.0))
[Link](X,Y,U,V)
[Link](-10,10)
[Link]('x')
[Link]('y')
[Link]('image_streamlines.png')
[Link]()