Les fonctions de courant Ψ /

potentielle Φ
Exemples de méthodes numériques
pour résoudre les EDP’S
S9
 La fonction de courant ψ
En deux dimensions, une méthode pour satisfaire
l'équation de continuité consiste à remplacer les deux
composantes inconnues de la vitesse u et v par une
seule fonction inconnue, ce qui réduit le nombre
d'inconnues et simplifie donc l'analyse d'un problème
d'écoulement de fluide idéal.
La fonction de courant ψ est l'équation qui représente
toutes les équations des lignes de courant.
En deux dimensions, elle est une fonction de x et y, et
pour l'équation de chaque ligne de courant, elle est
égale à une constante spécifique ψ(x,y) = C.
Composantes de la vitesse.
 Par définition, la vitesse d'une
particule de fluide est toujours
tangente à la ligne de courant le
long de laquelle elle se déplace.
Par conséquent, nous pouvons
relier les composantes de vitesse
u et v à la pente de la tangente
par proportion. Comme le
montre la figure , ou:
Maintenant, si nous prenons la dérivée
totale de l'équation de la ligne de courant
ψ(x,y) = C, qui décrit la ligne de courant

les deux composantes de la vitesse peuvent

être reliées à ψ. Donc par comparaison on
trouve :
et v
Par conséquent, si nous connaissons l'équation
d'une ligne de courant,ψ(x,y)=C , les composantes
de la vitesse d'une particule qui se déplace le long
de cette ligne peuvent être obtenues à l'aide de
ces équations.
 En obtenant les composantes de vitesse de cette
manière, nous pouvons montrer que pour un
écoulement permanent, la fonction de courant ψ
satisfait automatiquement l'équation de continuité.

et v

)=0
Nous verrons plus loin que dans certains
problèmes, il est commode d'exprimer la
fonction de courant Ψ et les composantes de
la vitesse en fonction de leurs coordonnées
polaires, r et θ.
Sans preuve, si ψ(r,θ)= C est donné, alors les
composantes radiales et transversales de la
vitesse sont:
,
 Débit volumétrique.
 La fonction d'écoulement peut également être
utilisée pour déterminer le débit volumétrique
entre deux lignes d'écoulement quelconques. Par
exemple, considérons le volume de contrôle
différentiel triangulaire de la figure qui est situé à
l'intérieur du tube entre les deux lignes de courant
ψ et ψ + dψ.
 Comme il s'agit d'un écoulement bidimensionnel,
nous considérerons l'écoulement dq à travers cet
élément en tant que mesure du débit par unité de
profondeur z, c'est-à-dire en m2/s.
 On suppose que l’écoulement est incompressible
et permanent :
Q (entrant)=Q (sortant)

q=
 LA FONCTION POTENTIELLE
 Une autre façon de relier les composantes de la vitesse à
une fonction unique est d'utiliser la fonction potentielle Φ où
Φ = Φ(x,y).
 Les composantes de la vitesse sont déterminées à partir de
Φ(x,y). à l'aide des équations suivantes.

 La vitesse résultante est donc:

 La fonction potentielle Φ ne décrit que l'écoulement

irrotationnel.
( )- ( ))=0
 Classification des équations différentielles
partielles

• La caractérisation dépend des racines des termes d'ordre supérieur

(ici d'ordre 2) :
• > 0 alors l'équation est dite hyperbolique.
• = 0  alors l'équation est dite parabolique.
• < 0  alors l'équation est dite elliptique.
• Remarque : si a, b et c dépendent eux-mêmes de x et y, les
équations peuvent être de types différents, en fonction de la
position dans l'espace x-y. Dans ce cas, les équations sont de type
mixte..
 Problèmes elliptiques
 Les équations elliptiques sont caractéristiques des problèmes d'équilibr
ce qui inclut de nombreux écoulements à l'état permanent (mais pas
tous).
 Les exemples sont les flux potentiels, la distribution de la température a
régime permanent dans une tige de matériau solide et les distributions
de contraintes à l'équilibre dans des objets solides soumis à des charge
appliquées.
 Pour les écoulements potentiels, la vitesse est exprimée en termes de
potentiel de vitesse : . L'écoulement étant incompressible, =0,
ce qui donne =0 .
 Cette équation est également connue sous le nom d'équation de
Laplace. également connue sous le nom d'équation de Laplace :
=0
𝟐 𝟐
𝟐 + 𝟐 =0
 Problèmes paraboliques
 Les équations paraboliques décrivent les problèmes de marche. Il s'agit
notamment de problèmes dépendant du temps qui impliquent des
quantités importantes de dissipation.
 Les exemples sont les écoulements visqueux instationnaire et la conduction
thermique instationnaire.
 L'écoulement de la couche limite visqueuse en régime permanent est
également parabolique (marche le long de la ligne de courant, pas dans
le temps).
 Un exemple est la distribution transitoire de la température dans une tige
de refroidissement :
Les équations hyperboliques
 Les équations hyperboliques sont typiques des
problèmes de marche avec une dissipation
négligeable.
 L'équation des ondes est un exemple :

• Il s'agit de décrire le déplacement transversal d'une

corde lors de vibrations de faible amplitude. Si y est le
déplacement, x la coordonnée le long de la corde, et a
l'amplitude initiale.
• Les ondes de choc dans les écoulements transsoniques
et supersoniques sont des exemples d'écoulements régis
par des équations hyperboliques.
 Méthodes de différences finies
 Après avoir développé un modèle
mathématique (c'est-à-dire les
équations directrices, les conditions
aux limites, la géométrie, etc.), nous
devons résoudre ces équations - cette
étape s'appelle la discrétisation.
 Nous apprendrons la méthode des
différences finies (FDM)
 L'objectif de la méthode des
différences finies est d'obtenir une
approximation des dérivées
 Différence en avant
Différence en arrière:
Différence centrée :
Différence centrée :
Exemple 1
Convertir les dérivées en une relation algébrique
entre des points voisins
La discrétisation temporelle

Différence en avant
Différence en arrière
Les erreurs numériques :

Les normes
 Stabilité

Un schéma de
différences finies est dit
stable au sens naïf du
terme lorsque les erreurs
d'arrondi ou de
troncature n'augmentent
pas indéfiniment.
 Analyse de la stabilité des EDP
 En général, il est très difficile d'effectuer une analyse de
stabilité perspicace et en forme fermée des EDP non
linéaires multidimensionnelles.
 Nous avons généralement recours à des équations
modèles - telles que les équations d'advection et de
diffusion d'advection - pour tester les idées sur la
stabilité.
 Un outil utile pour les EDP linéaires est connu sous le
nom d'analyse de stabilité de von-Neumann.

=
 La consistance de la solution numérique
 Une approximation par différences finies d'une EDP
linéaire continue est dite consistante si l'erreur de
troncature s'approche de zéro lorsque le pas de temps
et l'espacement de la grille s'approchent tous deux de
zéro
 Un schéma numérique est consistant si son opérateur
discret (avec les différences finies) converge vers
l'opérateur continu (avec les dérivées) de l'EDP pour
(erreur de troncature nulle).
 Théorème d'équivalence de Lax-Richtmyer: Pour une
EDP linéaire bien posée, un schéma de différences finies
consistant est convergent si, et seulement si, il est stable.
Analyse dimensionnelle/modélisation
physique/Les similitudes
S10
Introduction
 Dans les chapitres précédents, nous avons présenté de
nombreuses équations importantes de la mécanique des
fluides et illustré leur application à la résolution de certains
problèmes pratiques d'ingénierie.
 Dans tous ces cas, nous avons pu obtenir une solution
algébrique sous la forme d'une équation décrivant
l'écoulement.
 Dans certains cas, cependant, un problème peut impliquer
un écoulement très compliqué, et la combinaison de
variables telles que la vitesse, la pression, la masse volumique,
la viscosité, etc. qui le décrivent peut ne pas être entièrement
comprise.
 Dans ce cas, l'écoulement est souvent étudié par le biais
d'une expérience.
Malheureusement, le travail expérimental peut être
coûteux et prendre du temps, et il devient donc
important de pouvoir minimiser la quantité de données
expérimentales ,
 La meilleure façon d'y parvenir est d'effectuer une
analyse dimensionnelle de toutes les variables
pertinentes.
 Plus précisément, l'analyse dimensionnelle est une
branche des mathématiques qui permet d'organiser
toutes ces variables en ensemble de groupes sans
dimension.
Une fois ces groupes obtenus, nous pouvons les utiliser
pour obtenir le maximum d'informations à partir d'un
nombre minimum d'expériences.
 Les nombres adimensionnels
La méthode d'application de l'analyse dimensionnelle, a
été développée à l'origine par Lord Rayleigh.
Elle a ensuite été améliorée par Edgar Buckingham, sa
méthode consiste à combiner les variables qui décrivent
l'écoulement en un ensemble de ratios ou de "nombres"
sans dimension.
 Ces nombres sont souvent représentés par un rapport
des forces qui agissent sur le fluide à l'intérieur de
l'écoulement.
Et comme ces rapports apparaissent fréquemment dans
l'étude de la mécanique des fluides expérimentale, nous
présenterons ici certains des plus importants d'entre eux.
LE THÉORÈME DE BUCKINGHAM Pi
Le théorème Pi de Buckingham stipule que si un
phénomène d'écoulement dépend de n variables
physiques, telles que la vitesse, la pression et la
viscosité, et si ces variables physiques comportent
m dimensions primaires, telles que M, L et T, alors
l'analyse dimensionnelle permet d'organiser les n
variables en (n - m) nombres indépendants sans
dimension ou groupements.
 Chacun de ces groupements est appelé un terme
Π (Pi), car ce symbole est utilisé en mathématiques
pour symboliser un produit.
Une fois que la relation fonctionnelle entre les termes Π est
établie, elle peut être étudiée expérimentalement pour
voir comment elle est liée au comportement du flux à
l'aide de modèles.
Les groupements qui ont le plus d'influence sont retenus,
et ceux qui n'ont qu'un faible effet sur le flux sont rejetés.
En fin de compte, ce processus aboutira à une équation
empirique, dont les coefficients et exposants inconnus
seront ensuite déterminés par des expériences
supplémentaires.
La preuve du théorème de Buckingham Pi est assez
longue et peut être trouvée dans les livres relatifs à
l'analyse dimensionnelle. Nous ne nous intéresserons ici
qu'à son application.
La traînée de cette voiture est influencée par la masse
volumique ρ, la viscosité et la vitesse de l'air V, ainsi que
par la surface projetée S de la voiture dans l'écoulement.
(© Takeshi Takahara/Science Source)
 Exemple 1:

Etablissez le nombre de Reynolds

pour un fluide s'écoulant dans la
conduite , à l'aide de l'analyse
dimensionnelle, en sachant que
l'écoulement est fonction de la
masse volumique ρ et de la viscosité
du fluide, ainsi que de sa vitesse V
et du diamètre D de la conduite.
 SIMILITUDE

En raison de la complexité de
l'écoulement, l'effet du vent sur
les immeubles de grande
hauteur est souvent étudié à
l'aide d'un modèle réduit en
soufflerie.
 Les modèles réduits sont souvent utilisés pour étudier
l'écoulement tridimensionnel autour d'un objet à la
géométrie complexe, tel qu'un bâtiment, une automobile
ou un avion.
 En effet, il peut être assez difficile de décrire l'écoulement
à l'aide d'une solution analytique ou numérique.
 Même si l'écoulement peut être décrit par une analyse
numérique, dans les cas complexes, il convient de le valider
par une étude expérimentale correspondante, en utilisant
un modèle pour vérifier les résultats.
 Cette démarche est nécessaire pour la simple raison que
les hypothèses formulées dans le cadre d'une étude
numérique peuvent ne pas refléter fidèlement la situation
réelle, compte tenu de la complexité de l'écoulement.
 La similitude
La similitude est un processus mathématique qui
permet de s'assurer que c'est bien le cas. Elle exige
que le modèle et le flux qui l'entoure maintiennent non
seulement une similitude géométrique avec le
prototype, mais aussi une similitude cinématique et
dynamique.
Similitude géométrique: Si le modèle et l'écoulement
du fluide sont géométriquement similaires au
prototype, toutes les dimensions linéaires du modèle
doivent être dans les mêmes proportions que celles du
prototype, et tous ses angles doivent être identiques.
Le degré de réalisation de la
similitude géométrique dépend
du type de problème et de la
précision requise pour sa
résolution. Par exemple, une
similitude géométrique exacte
exige également que la rugosité
de la surface du modèle soit
proportionnelle à celle du
prototype.
 Dans certains cas, cependant, cela n'est pas possible, car la taille
réduite d'un modèle peut exiger que sa surface soit excessivement
lisse.
 De même, dans certaines modélisations, l'échelle verticale du modèle
peut devoir être exagérée pour produire les caractéristiques
d'écoulement appropriées.
 C'est le cas des études de rivières, dont le fond peut être difficile à
mettre à l'échelle.
 Similitude cinématique. Les principales dimensions définissant la
similitude cinématique sont la longueur et le temps. Par exemple, la
vitesse du fluide aux points correspondants entre le modèle de jet et
son prototype doit avoir une magnitude proportionnelle et être dans la
même direction.
 Puisque la vitesse dépend de la longueur et du temps, V = L/T, alors:
Similitude dynamique.
Pour que les lignes de courant autour du
prototype et de son modèle soient semblables, il
faut que les forces agissant sur les particules de
fluide correspondantes soient proportionnelles
dans les deux cas. Comme indiqué
précédemment, la force d'inertie Fi est
généralement considérée comme la force la plus
importante qui influence l'écoulement des fluides
autour d'un objet. C'est pourquoi il est d'usage
d'utiliser cette force, ainsi que chacune des
autres forces F qui influencent l'écoulement, pour
produire les rapports de similitude dynamique
entre le modèle et le prototype. Nous pouvons
exprimer symboliquement chaque rapport de

force comme suit =

 Écoulement permanent dans une
conduite.
 Si l'écoulement de l'eau dans une conduite doit être étudié à
l'aide d'un modèle, l'expérience a montré que les forces d'inertie
et visqueuses sont importantes.
 Par conséquent, pour l'écoulement dans une conduite, le nombre
de Reynolds doit être le même pour le modèle et le prototype afin
d'obtenir une similitude dynamique. Pour ce rapport, le diamètre D
de la conduite devient la "longueur caractéristique" L, et donc

Souvent, le même fluide est utilisé pour le modèle et le

prototype, de sorte que les propriétés ρ et sont identiques. Si
c'est le cas, alors
Écoulement a surface libre :
Pour l'écoulement dans un canal ouvert,
les forces causées par la tension
superficielle et
la compressibilité peuvent être
négligées sans introduire d'erreurs
appréciables, et si la longueur du canal
est courte , les pertes par frottement
seront négligeables, et les effets
visqueux peuvent donc également être
négligés. Par conséquent, l'écoulement
dans un canal ouvert est principalement
régi par la force d'inertie et la gravité.
Par conséquent, le nombre de Froude
peut être utilisé pour établir la similitude
entre le prototype et son modèle.
Pour ce nombre, la profondeur du fluide, h,
dans le canal est choisie comme "longueur
caractéristique" L, et comme g est le même
dans les deux cas, nous avons:
odèle hydraulique du barrage de Bluestone Lake sur la New River près
Hinton, en Virginie-Occidentale. L'échelle du modèle est de 1:65 à la
s verticalement et horizontalement, et le nombre de Reynolds, bien que
s inférieur à la valeur du prototype, est suffisamment élevé pour que le
t soit turbulent. (Avec l'aimable autorisation du U.S. Army Corps of
gineers Waterways Experiment Station).
Exemple 2
Un modèle de la voiture doit être construit à
l'échelle 1/4 et doit être testé à 20°C dans un tunnel
d'eau. Déterminer la vitesse requise de l'eau si la
voiture réelle se déplace à 30 m/s dans l'air à cette
même température.
la viscosité crée la force prédominante et nous
devons donc utiliser la similitude du nombre de
Reynolds
Un avion de chasse testé dans un voiture testée dans un canal d'eau
tunnel d'eau Image de la NASA

Sous-marin testé dans un canal de

vent . image de la NASA
paramètre Description Dimension
d'échelle
L longueur
caractéristique
V Vitesse
caractéristique
f Fréquence
caractéristique
différence de
pression de
référence
g l'accélération
gravitationnelle
 Equations de Navier Stokes adimensionnelles

St] ]
Exemple 1:
l'eau s'écoule sur un objet a 20 °c on donne ρ=980 kg/ et
µ= kg/m.s
 Le diamètre de l'objet est : D=0.6 m
 La vitesse loin de l’obstacle V=4 m/s
 La fréquence de vibration de l’objet est f=4Hz
 la différence de pression entre la pression de stagnation et
la pression statique du courant libre est d'environ 8000 Pa

Comparer l’ordre de grandeur de chaque terme dans les

équations de Stokes de Navier adimensionnelles