Annales de mathématiques pures et appliquées/Tome 14/Hydrodynamique, article 1

Recherches sur les lois générales du mouvement des fluides ; par M. F. Sarrus.

Annales de mathématiques pures et appliquées, Tome 14, 1823-1824 (p. 229-267).

HYDRODYNAMIQUE.

Recherches sur les lois générales du mouvement des fluides ;

Par M. F. Sarrus, docteur ès sciences.

≈≈≈≈≈≈≈≈≈

1. La théorie du mouvement des fluides est une des plus épineuses de la mécanique. On connaît depuis long-temps les équations différentielles de ce mouvement ; mais elles se montrent, pour la plupart, extrêmement rebelles, et ce n’est que dans quelques cas particuliers seulement que les géomètres modernes sont parvenus à les intégrer. En étudiant les découvertes dont ils ont enrichi cette partie de la science, nous avons cru reconnaître qu’ils avaient donné trop d’étendue à une proposition fort importante, et qu’en outre, la démonstration qu’ils en avaient donnée n’était pas à l’abri de toute objection. En cherchant à remplir cette lacune, nos recherches sur le même sujet se sont peu à peu étendues ; et nous sommes ainsi parvenu à quelques résultats nouveaux qui nous ont paru assez utiles ou du moins assez curieux pour mériter d’être connus. C’est l’ensemble de tout notre travail que nous nous proposons ici de mettre sous les yeux du lecteur.

2. Commençons par bien fixer le sens des notations dont nous nous proposons de faire usage dans la suite de cet écrit. Le temps sera constamment désigné par $t\,;$ et nous nommerons $x,y,z$ les coordonnées rectangulaires, pour l’époque $t,$ d’une molécule quelconque $M$ de la masse du fluide.

Toute quantité relative à cette molécule pourra, pour cette même époque $t,$ être considérée comme une fonction déterminée des variables $t,x,y,z,$ et de constantes qui devront être les mêmes pour toutes les autres molécules. C’est dans cette supposition, en en regardant les variables $t,x,y,z,$ comme absolument indépendantes, que seront prises, suivant les principes du calcul différentiel partiel, les dérivées

${\frac {\operatorname {d} u}{\operatorname {d} t}},\quad {\frac {\operatorname {d} u}{\operatorname {d} x}},\quad {\frac {\operatorname {d} u}{\operatorname {d} y}},\quad {\frac {\operatorname {d} u}{\operatorname {d} z}},$

quelle que puisse être d’ailleurs la fonction $u.$

3. Les coordonnées $x,y,z$ sont nécessairement fonctions du temps $t$ et des coordonnées primitives. Nommant ces dernières $a,b,c,$ toute quantité relative à la même molécule $M$ pourra être encore considérée comme une fonction déterminée de $t,a,b,c.$ En conséquence, pour éviter toute méprise, et ne pas confondre les valeurs de ${\frac {\operatorname {d} u}{\operatorname {d} t}},\ {\frac {\operatorname {d} ^{2}u}{\operatorname {d} t^{2}}},\ {\frac {\operatorname {d} ^{3}u}{\operatorname {d} t^{3}}}\ldots$ prises dans la première hypothèse avec ce qu’elles doivent être dans celle-ci, nous représenterons ces dernières respectivement par $\operatorname {d} 'u,\operatorname {d} '^{2}u,\operatorname {d} '^{3}u,\ldots .$ Quant aux coefficiens différentiels

${\frac {\operatorname {d} u}{\operatorname {d} a}},\quad {\frac {\operatorname {d} u}{\operatorname {d} b}},\quad {\frac {\operatorname {d} u}{\operatorname {d} c}},$

comme on ne saurait les confondre avec les dérivées d’une autre nature, nous n’emploîrons point de notations particulières pour les désigner.

On voit d’après cela que les vitesses de la molécule $M,$ parallèles aux coordonnées $x,y,z$ seront représentées respectivement par $\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z\,;$ que sa vitesse absolue aura conséquemment pour expression

${\sqrt {\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}}}$

et que les forces accélératrices correspondantes seront représentées par $\operatorname {d} '^{2}x+\operatorname {d} '^{2}y+\operatorname {d} '^{2}z.$

4. On peut observer, d’après ces conventions, que $\operatorname {d} 'u\operatorname {d} t$ n’est autre chose que la différentielle totale de $u,$ prise en regardant $x,y,z$ comme fonctions de $t,$ et que, par suite, on doit avoir identiquement

$\operatorname {d} 'u={\frac {\operatorname {d} u}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} u}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} u}{\operatorname {d} z}}\operatorname {d} 'z.$

On trouvera aussi

${\begin{aligned}{\frac {\operatorname {d} u}{\operatorname {d} a}}&={\frac {\operatorname {d} u}{\operatorname {d} x}}{\frac {\operatorname {d} x}{\operatorname {d} a}}+{\frac {\operatorname {d} u}{\operatorname {d} y}}{\frac {\operatorname {d} y}{\operatorname {d} a}}+{\frac {\operatorname {d} u}{\operatorname {d} z}}{\frac {\operatorname {d} z}{\operatorname {d} a}},\\\\{\frac {\operatorname {d} u}{\operatorname {d} b}}&={\frac {\operatorname {d} u}{\operatorname {d} x}}{\frac {\operatorname {d} x}{\operatorname {d} b}}+{\frac {\operatorname {d} u}{\operatorname {d} y}}{\frac {\operatorname {d} y}{\operatorname {d} b}}+{\frac {\operatorname {d} u}{\operatorname {d} z}}{\frac {\operatorname {d} z}{\operatorname {d} b}},\\\\{\frac {\operatorname {d} u}{\operatorname {d} c}}&={\frac {\operatorname {d} u}{\operatorname {d} x}}{\frac {\operatorname {d} x}{\operatorname {d} c}}+{\frac {\operatorname {d} u}{\operatorname {d} y}}{\frac {\operatorname {d} y}{\operatorname {d} c}}+{\frac {\operatorname {d} u}{\operatorname {d} z}}{\frac {\operatorname {d} z}{\operatorname {d} c}}\,;\end{aligned}}$

transformations dont nous ferons usage par la suite.

5. D’après les principes connus du calcul différentiel partiel, ou aura identiquement

${\begin{array}{rrr}{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} x}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} t}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} y}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} y\operatorname {d} t}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} z}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} z\operatorname {d} t}},\\\\{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} y}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} y\operatorname {d} x}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} z}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} z\operatorname {d} x}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} y\operatorname {d} z}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} z\operatorname {d} y}},\end{array}}$

et encore

${\begin{array}{rrr}\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} a}}={\frac {\operatorname {dd} 'u}{\operatorname {d} a}},&\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} b}}={\frac {\operatorname {dd} 'u}{\operatorname {d} b}},&\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} c}}={\frac {\operatorname {dd} 'u}{\operatorname {d} c}},\\\\{\frac {\operatorname {d} ^{2}u}{\operatorname {d} a\operatorname {d} b}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} b\operatorname {d} a}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} a\operatorname {d} c}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} c\operatorname {d} a}},&{\frac {\operatorname {d} ^{2}u}{\operatorname {d} b\operatorname {d} c}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} c\operatorname {d} b}}\,;\end{array}}$

mais il faudrait bien se garder de croire que l’on a de même

$\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} t}}={\frac {\operatorname {dd} 'u}{\operatorname {d} t}},\ \operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} x}}={\frac {\operatorname {dd} 'u}{\operatorname {d} x}},\ \operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} y}}={\frac {\operatorname {dd} 'u}{\operatorname {d} y}},\ \operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} z}}={\frac {\operatorname {dd} 'u}{\operatorname {d} z}},$

En effet, en vertu du n.^o 4, on a, par exemple,

$\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} t^{2}}}+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} t}}\operatorname {d} 'x+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} y\operatorname {d} t}}\operatorname {d} 'y+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} z\operatorname {d} t}}\operatorname {d} 'z,$

tandis que l’équation identique

$\operatorname {d} 'u={\frac {\operatorname {d} u}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} u}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} u}{\operatorname {d} z}}\operatorname {d} 'z,$

différentiée par rapport à $t,$ en regardant $x,y,z$ comme constans donne

${\frac {\operatorname {dd} 'u}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}u}{\operatorname {d} t^{2}}}+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} t\operatorname {d} z}}\operatorname {d} 'z$

$+{\frac {\operatorname {d} u}{\operatorname {d} x}}{\frac {\operatorname {dd} 'x}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} y}}{\frac {\operatorname {dd} 'y}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} z}}{\frac {\operatorname {dd} 'z}{\operatorname {d} t}},$

et par suite

${\frac {\operatorname {dd} 'u}{\operatorname {d} t}}=\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} x}}{\frac {\operatorname {dd} 'x}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} y}}{\frac {\operatorname {dd} 'y}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} z}}{\frac {\operatorname {dd} 'z}{\operatorname {d} t}},$

et non pas simplement

${\frac {\operatorname {dd} 'u}{\operatorname {d} t}}=\operatorname {d} '{\frac {\operatorname {d} u}{\operatorname {d} t}}.$

6. Nous supposerons que la caractéristique $\delta$ des variations est uniquement relative aux coordonnées primitives $a,b,c,$ et non au tems $t.$ Dès lors nous aurons

${\begin{aligned}\delta x&={\frac {\operatorname {d} x}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} x}{\operatorname {d} b}}\delta b+{\frac {\operatorname {d} x}{\operatorname {d} c}}\delta c,\\\\\delta y&={\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+{\frac {\operatorname {d} y}{\operatorname {d} c}}\delta c,\\\\\delta z&={\frac {\operatorname {d} z}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} z}{\operatorname {d} b}}\delta b+{\frac {\operatorname {d} z}{\operatorname {d} c}}\delta c,\\\\\end{aligned}}$

et en général

$\delta u={\frac {\operatorname {d} u}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} u}{\operatorname {d} y}}\delta y+{\frac {\operatorname {d} u}{\operatorname {d} z}}\delta z={\frac {\operatorname {d} u}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} u}{\operatorname {d} b}}\delta b+{\frac {\operatorname {d} u}{\operatorname {d} c}}\delta c\,;$

d’où il est facile de conclure

$\delta \operatorname {d} 'u=\operatorname {d} '\delta u,\qquad \delta {\frac {\operatorname {d} u}{\operatorname {d} t}}={\frac {\operatorname {d} \delta u}{\operatorname {d} t}},$

et toutes les autres relations qu’on en peut déduire par des différentiations répétées.

Si l’on avait

$\psi =\delta \phi +H,$

on en tirerait

$\operatorname {d} '\psi =\operatorname {d} '\delta \phi +\operatorname {d} 'H,$

et

${\frac {\operatorname {d} \psi }{\operatorname {d} t}}={\frac {\operatorname {d} \delta \phi +}{\operatorname {d} t}}+{\frac {\operatorname {d} H}{\operatorname {d} t}}.$

Dans les cas où $H$ sera une fonction de $a,b,c,$ sans $t,$ on aura $\operatorname {d} 'H=0,$ et par suite

$\operatorname {d} '\psi =\operatorname {d} '\delta \phi =\delta \operatorname {d} '\phi .$

Si, au contraire, $H$ est une fonction de $x,y,z,$ sans $t,$ on aura ${\frac {\operatorname {d} H}{\operatorname {d} t}}=0,$ et par suite

${\frac {\operatorname {d} \psi }{\operatorname {d} t}}={\frac {\operatorname {d} \delta \phi }{\operatorname {d} t}}=\delta {\frac {\operatorname {d} \phi }{\operatorname {d} t}}\,;$

et comme, dans ce cas ou les autres, $H$ peut ne pas être une variation exacte, ce qui pourtant serait nécessaire pour que $\psi$ fut une pareille variation, on en conclura que $\psi$ peut ne pas être une variation exacte, bien que l’une ou l’autre des dérivées $\operatorname {d} '\psi$ ou ${\frac {\operatorname {d} \psi }{\operatorname {d} t}}$ soit une telle variation.

C’est pour avoir négligé de faire cette observation, ou peut-être même pour avoir cru que le contraire était vrai, que les divers auteurs qui ont écrit sur l’hydrodynamique ont réputé, comme, généralement démontrée, une proposition qui n’est vraie que dans des cas particuliers.

7. Ces notations et remarques ainsi établies, soit $p$ la pression qu’éprouve la molécule $M,$ au bout du temps $t\,;$ soit $\Delta$ la densité correspondante de la même molécule ; soient, en outre, $X,Y,Z$ les trois forces parallèles aux axes des coordonnées auxquelles peuvent être réduites toutes celles qui la sollicitent, forces que nous réputerons positives ou négatives, suivant qu’elles tendront à augmenter ou à diminuer les coordonnées $x,y,z$ du point $M.$ Si les fluide était en équilibre, nous aurions, pour les équations de cet équilibre

$X-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}=0,\quad Y-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}=0,\quad Z-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}=0.$

Si, au contraire, le fluide est en mouvement, les forces accélératrices qui forment les premiers membres de ces trois équations ne seront pas nulles ; et comme ces mêmes forces accélératrices peuvent aussi être exprimées respectivement par $\operatorname {d} '^{2}x,\operatorname {d} '^{2}y,\operatorname {d} '^{2}z\,;$ il s’ensuit que les équations du mouvement sont

${\begin{aligned}X-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}&=\operatorname {d} '^{2}x,\\\\Y-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}&=\operatorname {d} '^{2}y,\\\\Z-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}&=\operatorname {d} '^{2}z.\end{aligned}}$

8. Dans la nature, les variations

$X\delta x+Y\delta y+Z\delta z,\qquad {\frac {\delta p}{\Delta }},$

sont exactes : et nous les supposerons telles dans tout ce qui va suivre. Posant donc

$X\delta x+Y\delta y+Z\delta z=\delta \nu \,;$

d’où

$X={\frac {\operatorname {d} \nu }{\operatorname {d} x}},\quad Y={\frac {\operatorname {d} \nu }{\operatorname {d} y}},\quad Z={\frac {\operatorname {d} \nu }{\operatorname {d} z}},$

les équations du mouvement deviendront

$\left.{\begin{aligned}{\frac {\operatorname {d} \nu }{\operatorname {d} x}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}&=\operatorname {d} '^{2}x,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} y}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}&=\operatorname {d} '^{2}y,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} z}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}&=\operatorname {d} '^{2}z.\end{aligned}}\right\}\quad (1)$

Soit prise successivement la somme de leurs produits d’abord par

${\frac {\operatorname {d} x}{\operatorname {d} a}},\quad {\frac {\operatorname {d} y}{\operatorname {d} a}},\quad {\frac {\operatorname {d} z}{\operatorname {d} a}},$

puis par

${\frac {\operatorname {d} x}{\operatorname {d} b}},\quad {\frac {\operatorname {d} y}{\operatorname {d} b}},\quad {\frac {\operatorname {d} z}{\operatorname {d} b}},$

puis enfin par

${\frac {\operatorname {d} x}{\operatorname {d} c}},\quad {\frac {\operatorname {d} y}{\operatorname {d} c}},\quad {\frac {\operatorname {d} z}{\operatorname {d} c}}\,;$

et réduisant à l’aide des relations du n.^o 4, on trouvera

$\left.{\begin{aligned}{\frac {\operatorname {d} \nu }{\operatorname {d} a}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} a}}&=\operatorname {d} '^{2}x.{\frac {\operatorname {d} x}{\operatorname {d} a}}+\operatorname {d} '^{2}y.{\frac {\operatorname {d} y}{\operatorname {d} a}}+\operatorname {d} '^{2}z.{\frac {\operatorname {d} z}{\operatorname {d} a}},\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} b}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} b}}&=\operatorname {d} '^{2}x.{\frac {\operatorname {d} x}{\operatorname {d} b}}+\operatorname {d} '^{2}y.{\frac {\operatorname {d} y}{\operatorname {d} b}}+\operatorname {d} '^{2}z.{\frac {\operatorname {d} z}{\operatorname {d} b}},\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} c}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} c}}&=\operatorname {d} '^{2}x.{\frac {\operatorname {d} x}{\operatorname {d} c}}+\operatorname {d} '^{2}y.{\frac {\operatorname {d} y}{\operatorname {d} c}}+\operatorname {d} '^{2}z.{\frac {\operatorname {d} z}{\operatorname {d} c}}\,;\end{aligned}}\right\}\quad (2)$

forme sous laquelle on devra employer les équations, lorsqu’on voudra déterminer $x,y,z,p,$ en fonction de $a,b,c,t\,;$ mais, si on veut déterminer $p,\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z,$ en fonction de $t,x,y,z,$ il faudra mettre les, équations (1) sous la forme

$\left.{\begin{aligned}{\frac {\operatorname {d} \nu }{\operatorname {d} x}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}&={\frac {\operatorname {dd} 'x}{\operatorname {d} t}}+{\frac {\operatorname {dd} 'x}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {dd} 'x}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {dd} 'x}{\operatorname {d} z}}\operatorname {d} 'z,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} y}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}&={\frac {\operatorname {dd} 'y}{\operatorname {d} t}}+{\frac {\operatorname {dd} 'y}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {dd} 'y}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {dd} 'y}{\operatorname {d} z}}\operatorname {d} 'z,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} z}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}&={\frac {\operatorname {dd} 'z}{\operatorname {d} t}}+{\frac {\operatorname {dd} 'z}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {dd} 'z}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {dd} 'z}{\operatorname {d} z}}\operatorname {d} 'z,\end{aligned}}\right\}\quad (3)$

qu’il sera facile de leur faire acquérir au moyen des relations du n.^o 4.

Les équations (2) ne sont pas en nombre suffisant pour déterminer les quatre inconnues $x,y,z,p,$ en fonction de $a,b,c,t\,;$ et il en est de même des équations (3) ; lorqu’on veut déterminer $p,\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z,$ en fonction de $t,x,y,z.$ Pour y suppléer, on y joint ordinairement celles qui peuvent résulter de la continuité du fluide.

9. La condition de continuité du fluide consiste en ce qu’un système qui, pour un instant donné, forme une masse finie et continue, a dû former jusqu’alors et devra former ensuite une pareille masse finie et continue.

On peut conclure de là que la masse d’un pareil système est invariable, et que sa surface est toujours formée des mêmes molécules, quelle que puisse être d’ailleurs sa forme primitive. Si donc le fluide est contenu dans un vase de figure quelconque, les molécules qui, pour un instant donné, sont contiguës aux parois du vase, n’ont pu auparavant et ne pourront ensuite que glisser contre ces parois.

Pour traduire ces conditions en langage analitique, considérons une partie quelconque (A) finie et continue du fluide en mouvement. Sa masse sera exprimée par l’intégrale

$\iiint \Delta \operatorname {d} x\operatorname {d} y\operatorname {d} z,$

prise entre les limites données par sa surface.

Cela posé, il résulte des principes connus sur la tranformation des intégrales que, si l’on fait, pour abréger,

$\left.{\begin{aligned}&{\frac {\operatorname {d} x}{\operatorname {d} a}}{\frac {\operatorname {d} y}{\operatorname {d} b}}{\frac {\operatorname {d} z}{\operatorname {d} c}}+{\frac {\operatorname {d} x}{\operatorname {d} b}}{\frac {\operatorname {d} y}{\operatorname {d} c}}{\frac {\operatorname {d} z}{\operatorname {d} a}}+{\frac {\operatorname {d} x}{\operatorname {d} c}}{\frac {\operatorname {d} y}{\operatorname {d} a}}{\frac {\operatorname {d} z}{\operatorname {d} b}}\\\\-&{\frac {\operatorname {d} x}{\operatorname {d} a}}{\frac {\operatorname {d} y}{\operatorname {d} c}}{\frac {\operatorname {d} z}{\operatorname {d} b}}-{\frac {\operatorname {d} x}{\operatorname {d} b}}{\frac {\operatorname {d} y}{\operatorname {d} a}}{\frac {\operatorname {d} z}{\operatorname {d} c}}-{\frac {\operatorname {d} x}{\operatorname {d} c}}{\frac {\operatorname {d} y}{\operatorname {d} b}}{\frac {\operatorname {d} z}{\operatorname {d} a}}\end{aligned}}\right\}=\theta ,$

l’intégrale précédente reviendra à

$\iiint \Delta \theta \operatorname {d} a\operatorname {d} b\operatorname {d} c,$

prise entre les limites données far la surface du fluide primitif (A), et comme, en désignant par $\Delta _{0}$ la valeur primitive de la densité $\Delta ,$ la masse primitive du même fluide est exprimée par

$\iiint \Delta _{0}\operatorname {d} a\operatorname {d} b\operatorname {d} c,$

on doit en conclure

$\iiint \theta \Delta \operatorname {d} a\operatorname {d} b\operatorname {d} c=\iiint \Delta _{0}\operatorname {d} a\operatorname {d} b\operatorname {d} c,$

et par suite

$\theta \Delta =\Delta _{0}\,;\qquad$(4)

en observant que les limites des intégrales qui composent l’équation précédente sont données par la surface primitive du fluide (A), laquelle surface est entièrement arbitraire.

Cette équation (4) sert à compléter le système d’équations (2) ; mais, pour pouvoir l’employer conjointement avec les équations (3), il faut la mettre sous une forme plus commode.

Comme l’origine des temps est entièrement arbitraire, si nous représentons par $\Delta _{1},x_{1},y_{1},z_{1}$ les valeurs de $\Delta ,x,y,z$ qui correspondent à la valeur $t+h$ de $t,$ nous devrons avoir encore

$\left\{{\begin{aligned}&{\frac {\operatorname {d} x_{1}}{\operatorname {d} x}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} y}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} z}}+{\frac {\operatorname {d} x_{1}}{\operatorname {d} y}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} z}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} x}}+{\frac {\operatorname {d} x_{1}}{\operatorname {d} z}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} x}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} y}}\\\\-&{\frac {\operatorname {d} x_{1}}{\operatorname {d} x}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} z}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} y}}-{\frac {\operatorname {d} x_{1}}{\operatorname {d} y}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} x}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} z}}-{\frac {\operatorname {d} x_{1}}{\operatorname {d} z}}{\frac {\operatorname {d} y_{1}}{\operatorname {d} y}}{\frac {\operatorname {d} z_{1}}{\operatorname {d} x}}\end{aligned}}\right\}\Delta _{1}=\Delta \,;$

mais, en développant suivant les puissances ascendantes de $h,$ nous avons

${\begin{aligned}\Delta _{1}&=\Delta +{\frac {h}{1}}\operatorname {d} '\Delta +{\frac {h^{2}}{1.2}}\operatorname {d} '^{2}\Delta +\ldots ,\\\\x_{1}&=x+{\frac {h}{1}}\operatorname {d} 'x+{\frac {h^{2}}{1.2}}\operatorname {d} '^{2}x+\ldots ,\\\\y_{1}&=y+{\frac {h}{1}}\operatorname {d} 'y+{\frac {h^{2}}{1.2}}\operatorname {d} '^{2}y+\ldots ,\\\\z_{1}&=z+{\frac {h}{1}}\operatorname {d} 'z+{\frac {h^{2}}{1.2}}\operatorname {d} '^{2}z+\ldots \,;\end{aligned}}$

substituant ces valeurs dans l’équation précédente, et ordonnant le résultat suivant les puissances de $h,$ nous trouverons

$\Delta +{\frac {h}{1}}\left\{\operatorname {d} '\Delta +\Delta \left({\frac {\operatorname {dd} 'x}{\operatorname {d} x}}+{\frac {\operatorname {dd} 'y}{\operatorname {d} y}}+{\frac {\operatorname {dd} 'z}{\operatorname {d} z}}\right)\right\}+{\frac {h^{2}}{1.2}}\left\{\ldots \right\}+\ldots =\Delta .$

Cette équation devant avoir identiquement lieu quelle que soit la valeur de $h,$ nous en conclurons

$0=\operatorname {d} '\Delta +\Delta \left({\frac {\operatorname {dd} 'x}{\operatorname {d} x}}+{\frac {\operatorname {dd} 'y}{\operatorname {d} y}}+{\frac {\operatorname {dd} 'z}{\operatorname {d} z}}\right).\qquad$(5)

Les coefficiens des puissances supérieures de $h$ nous donneraient bien une infinité d’autres équations ; mais elles ne seraient que les différentielles successives de l’équation (5), par rapport à la caractéristique $\operatorname {d} ',$ et conséquemment ne signifieraient rien de plus que cette équation, qu’on peut encore mettre sous la forme

$0={\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'y}{\operatorname {d} y}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'z}{\operatorname {d} z}},\qquad$(6)

en observant que

$\operatorname {d} '\Delta ={\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} \Delta }{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} \Delta }{\operatorname {d} z}}\operatorname {d} 'z.$

Lorsque le fluide est incompressible, la densité $\Delta$ de la molécule $M$ est invariable ; on a donc alors $\Delta =\Delta _{0}$ et $\operatorname {d} '\Delta =0.$ En conséquence, les équations (4) et (5) se réduisent à

$\left\{{\begin{aligned}&{\frac {\operatorname {d} x}{\operatorname {d} a}}{\frac {\operatorname {d} y}{\operatorname {d} b}}{\frac {\operatorname {d} z}{\operatorname {d} c}}+{\frac {\operatorname {d} x}{\operatorname {d} b}}{\frac {\operatorname {d} y}{\operatorname {d} c}}{\frac {\operatorname {d} z}{\operatorname {d} a}}+{\frac {\operatorname {d} x}{\operatorname {d} c}}{\frac {\operatorname {d} y}{\operatorname {d} a}}{\frac {\operatorname {d} z}{\operatorname {d} b}}\\\\-&{\frac {\operatorname {d} x}{\operatorname {d} a}}{\frac {\operatorname {d} y}{\operatorname {d} c}}{\frac {\operatorname {d} z}{\operatorname {d} b}}-{\frac {\operatorname {d} x}{\operatorname {d} b}}{\frac {\operatorname {d} y}{\operatorname {d} a}}{\frac {\operatorname {d} z}{\operatorname {d} c}}-{\frac {\operatorname {d} x}{\operatorname {d} c}}{\frac {\operatorname {d} y}{\operatorname {d} b}}{\frac {\operatorname {d} z}{\operatorname {d} a}}\end{aligned}}\right\}=1,\quad (7)$

$0={\frac {\operatorname {dd} 'x}{\operatorname {d} x}}+{\frac {\operatorname {dd} 'y}{\operatorname {d} y}}+{\frac {\operatorname {dd} 'z}{\operatorname {d} z}}\,;\qquad$(8)

en mettant au lieu de $\theta$ sa valeur.

De plus, dans le cas où le fluide ne sera point homogène, il faudra avoir égard à la condition que $\Delta$ doit se réduire à une fonction de $a,b,c$ sans $t,$ ou, ce qui revient au même, à l’équation

$0={\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} \Delta }{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} \Delta }{\operatorname {d} z}}\operatorname {d} 'z.$

Soit enfin $u=0$ l’équation en $t,x,y,z$ de la surface limite du même fluide (A) après le temps $t.$ Les valeurs primitives des coordonnées $x,y,z$ qui satisfont à cette équation devront être les coordonnées de la surface primitive de ce fluide (A), et par conséquent entièrement indépendantes de $t.$ Si donc on met dans cette équation $u=0$ au lieu de $x,y,z$ leurs valeurs en $a,b,c$ et $t,$ le temps $t$ devra disparaître du résultat ; condition qui servirait à déterminer la forme des fonctions arbitraires, dans les cas où l’on ferait usage des équations (2) et (4). Mais, si l’on faisait usage des équations (3) et (5), on observerait que, puisque $u$ doit se réduire à une fonction de $a,b,c$ sans $t,$ l’on doit avoir

$\operatorname {d} 'u={\frac {\operatorname {d} u}{\operatorname {d} t}}+{\frac {\operatorname {d} u}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} u}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} u}{\operatorname {d} z}}\operatorname {d} 'z=0.$

Telles sont les équations générales du mouvement des fluides, et les principes qui doivent servir à la détermination des fonctions arbitraires introduites par les intégrations. Nous aurions pu sans doute prendre ces équations soit dans les mécaniques céleste ou analitique, soit dans la mécanique de M. Poisson ; mais nous avons pensé qu’on trouverait plus commode de rencontrer au commencement de cet essai les considérations et les procédés nécessaires pour les obtenir.

10. Reprenons les équations (1), savoir

${\begin{aligned}{\frac {\operatorname {d} \nu }{\operatorname {d} x}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}&=\operatorname {d} '^{2}x,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} y}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}&=\operatorname {d} '^{2}y,\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} z}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}&=\operatorname {d} '^{2}z\,;\end{aligned}}$

prenant la somme de leurs produits respectifs par $\delta x,\delta y,\delta z,$ ajoutant les produits et réduisant, au moyen des relations données dans le n.^o 6, nous trouverons

$\delta \nu -{\frac {\delta p}{\Delta }}=\operatorname {d} '^{2}x\delta x+\operatorname {d} '^{2}y\delta y+\operatorname {d} '^{2}z\delta z\,;\qquad$(9)

faisant donc, pour abréger,

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\psi ,\qquad$(10)

nous aurons, en différentiant par rapport à la caractéristique $\operatorname {d} ',$

$\operatorname {d} '^{2}x\delta x+\operatorname {d} '^{2}y\delta y+\operatorname {d} '^{2}z\delta z+\operatorname {d} 'x\operatorname {d} '\delta x+\operatorname {d} 'y\operatorname {d} '\delta y+\operatorname {d} 'z\operatorname {d} '\delta z=\operatorname {d} '\psi ,$

ou encore

$\operatorname {d} '^{2}x\delta x+\operatorname {d} '^{2}y\delta y+\operatorname {d} '^{2}z\delta z=\operatorname {d} '\psi -{\frac {1}{2}}\delta \left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right)\,;$

en observant que

${\begin{aligned}&\operatorname {d} 'x\operatorname {d} '\delta x=\operatorname {d} 'x\delta \operatorname {d} 'x={\frac {1}{2}}\delta .\operatorname {d} 'x^{2},\\\\&\operatorname {d} 'y\operatorname {d} '\delta y=\operatorname {d} 'y\delta \operatorname {d} 'y={\frac {1}{2}}\delta .\operatorname {d} 'y^{2},\\\\&\operatorname {d} 'z\operatorname {d} '\delta z=\operatorname {d} 'z\delta \operatorname {d} 'z={\frac {1}{2}}\delta .\operatorname {d} 'z^{2}.\end{aligned}}$

Au moyen de cette dernière valeur de la quantité

$\operatorname {d} '^{2}x\delta x+\operatorname {d} '^{2}y\delta y+\operatorname {d} '^{2}z\delta z,$

l’équation (9) deviendra

$\delta \nu -{\frac {\delta p}{\Delta }}=\operatorname {d} '\psi -{\frac {1}{2}}\delta \left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right)\,;\qquad$(11)

de sorte que, si l’on détermine $\phi$ au moyen de l’équation

$\nu -\int {\frac {\delta p}{\Delta }}=\operatorname {d} '\phi -{\frac {1}{2}}\left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right),\qquad$(12)

ce qui est toujours possible, au moins par les quadratures, et ce qui donne

$\delta \nu -{\frac {\delta p}{\Delta }}=\delta \operatorname {d} '\phi -{\frac {1}{2}}\delta \left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right),$

on trouvera, en comparant cette dernière avec l’équation (11),

$\operatorname {d} '\psi =\delta \operatorname {d} '\phi =\operatorname {d} '\delta '\phi ,\quad$d’où$\quad \psi =\delta \phi +H\,;\qquad$(13)

en désignant par $H$ une fonction de $a,b,c$ sans $t,$ de la forme $A\delta a+B\delta b+C\delta c,$ qui d’ailleurs peut être quelconque, puisqu’elle est introduite par l’intégration, et qu’il faudra conséquemment déterminer au moyen de l’état primitif du fluide.

Cette équation (13) va nous conduire à quelques résultats importans.

11. Si l’on nomme $\psi _{1}$ et $\phi _{1}$ les valeurs de $\psi$ et $\phi$ qui correspondent à une valeur particulière quelconque $t_{1}$ du temps, nous devrons avoir de même

$\psi _{1}=\delta \phi _{1}+H\,;$

éliminant $H$ entre cette équation et l’équation (13), nous aurons

$\psi -\psi _{1}=\delta \phi -\delta \phi _{1}=\delta (\phi -\phi _{1}),$

qui, traduite en langage ordinaire, donne ce théorème :

Si, pour deux valeurs quelconques du temps, on calcule celles de $\psi$ qui leur correspondent, la différence des résultats sera toujours une variation exacte.

Ce résultat, entièrement nouveau, nous paraît curieux. Malheureusement il n’est utile qu’autant qu’on peut en déduire, comme cas particulier, que $\psi$ sera une variation exacte dans tous les temps si, pour un instant quelconque, cette fonction est une variation exacte, ou bien si elle est nulle ; théorème important, déjà connu depuis long-temps, mais qui n’avait pas été démontré jusqu’ici d’une manière aussi simple.

12. L’utilité de l’équation (13) ne se borne pas à la démonstration du précédent théorème. Si, en effet, on y met pour $\psi$ et $H$ les quantités qu’elles représentent, elle devient

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\delta \phi +A\delta a+B\delta b+C\delta c,\qquad$(14)

ou encore

$\left.{\begin{aligned}&\left(\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} a}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} a}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} a}}\right)\delta a\\\\+&\left(\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} b}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} b}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} b}}\right)\delta b\\\\+&\left(\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} c}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} c}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} c}}\right)\delta c\end{aligned}}\right\}=\delta \phi +A\delta a+B\delta b+C\delta c\,;$

d’où en observant que les variations $\delta a,\delta b,\delta c$ sont entièrement arbitraires, on tirera

$\left.{\begin{aligned}&\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} a}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} a}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} a}}={\frac {\operatorname {d} \phi }{\operatorname {d} a}}+A,\\\\&\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} b}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} b}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} b}}={\frac {\operatorname {d} \phi }{\operatorname {d} b}}+B,\\\\&\operatorname {d} 'x{\frac {\operatorname {d} x}{\operatorname {d} c}}+\operatorname {d} 'y{\frac {\operatorname {d} y}{\operatorname {d} c}}+\operatorname {d} 'z{\frac {\operatorname {d} z}{\operatorname {d} c}}={\frac {\operatorname {d} \phi }{\operatorname {d} c}}+C,\end{aligned}}\right\}\quad (15)$

qui sont les intégrales premières des équations (2). Il resterait encore à déterminer $x,y,z$ et $\phi ,$ au moyen de ces équations combinées avec l’équation (4) ; après quoi l’équation (12) ferait connaître $p\,;$ mais cette détermination surpasse les forces de l’analise. Toutefois ces équations pourront être utiles, lorsqu’on connaîtra le mouvement que prendrait le fluide, sans l’action de très-petites forces perturbatrices, comme il arrive dans la théorie des marées, où, sans l’action du soleil et de la lune, la mer prendrait un mouvement uniforme de rotation autour de l’axe de la terre ; mouvement qui n’est qu’extrêmement peu troublé par l’effet de cette action.

13. Nous avons identiquement

$\operatorname {d} \phi ={\frac {\operatorname {d} \phi }{\operatorname {d} t}}+{\frac {\operatorname {d} \phi }{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} \phi }{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} \phi }{\operatorname {d} z}}\operatorname {d} 'z,$

par suite de quoi l’équation (12) devient

$\nu -\int {\frac {\delta p}{\Delta }}={\frac {\operatorname {d} \phi }{\operatorname {d} t}}+{\frac {\operatorname {d} \phi }{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} \phi }{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} \phi }{\operatorname {d} z}}\operatorname {d} 'z-{\frac {1}{2}}\left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right)$ (16)

si donc l’on fait

$\operatorname {d} 'x={\frac {\operatorname {d} \phi }{\operatorname {d} x}}+P,\quad \operatorname {d} 'y={\frac {\operatorname {d} \phi }{\operatorname {d} y}}+Q,\quad \operatorname {d} 'z={\frac {\operatorname {d} \phi }{\operatorname {d} z}}+R,\qquad$(17)

et qu’on substitue ces valeurs dans l’équation (16), elle deviendra, réductions faites,

$\nu -\int {\frac {\delta p}{\Delta }}={\frac {\operatorname {d} \phi }{\operatorname {d} t}}+{\frac {1}{2}}\left\{\left({\frac {\operatorname {d} \phi }{\operatorname {d} x}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} y}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} z}}\right)^{2}\right\}-{\frac {1}{2}}\left(P^{2}+Q^{2}+R^{2}\right).$ (18)

au moyen des mêmes substitutions, l’équation (5) devient

${\begin{array}{cr}0={\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} x}}{\frac {\operatorname {d} \phi }{\operatorname {d} x}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} y}}{\frac {\operatorname {d} \phi }{\operatorname {d} y}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} z}}{\frac {\operatorname {d} \phi }{\operatorname {d} z}}+\Delta \left({\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z^{2}}}\right)&\\&(19)\\+P{\frac {\operatorname {d} \Delta }{\operatorname {d} x}}+Q{\frac {\operatorname {d} \Delta }{\operatorname {d} y}}+R{\frac {\operatorname {d} \Delta }{\operatorname {d} z}}+\Delta \left({\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} Q}{\operatorname {d} y}}+{\frac {\operatorname {d} R}{\operatorname {d} z}}\right)\,;&\\\end{array}}$

et l’équation (7), qui a lieu quand le fluide est incompressible, deviendra

$0={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z^{2}}}+{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} Q}{\operatorname {d} y}}+{\frac {\operatorname {d} R}{\operatorname {d} z}}\,;\qquad$(20)

d’où l’on voit que les équations du mouvement du fluide se simplifieraient d’une manière notable, si l’on pouvait déterminer les fonctions $P,Q,R.$ On a pour cela l’équation

$P\delta x+Q\delta y+R\delta z=A\delta a+B\delta b+C\delta c,$

que l’on obtient en mettant dans inéquation (14) les valeurs de $\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z,$ données par les équations (17). Il ne sera guère possible d’en tirer les valeurs de $P,Q,R,$ du moins en général. Si cependant ces fonctions ne dépendaient que de $x,y,z,$ et qu’on eût

$A=\operatorname {F} (a,b,c),\qquad B=\operatorname {F} '(a,b,c),\qquad C=\operatorname {F} ''(a,b,c),$

en désignant par $\operatorname {F} ,\operatorname {F} ',\operatorname {F} ''$ des fonctions quelconques de $a,b,c,$ déterminées au moyen de l’état initial du fluide, on en tirerait

$P=\operatorname {F} (x,y,z),\qquad Q=\operatorname {F} '(x,y,z),\qquad R=\operatorname {F} ''(x,y,z).$

14. Lorsque $\psi ,$ ou, ce qui revient au même, $\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z,$ est une variation exacte, l’équation (13) fait voir qu’il doit en être de même de $H.$ On peut donc alors représenter cette fonction par $\delta K,\ K$ étant une fonction de $a,b,c,$ sans $t,$ convenablement déterminée. Dès-lors, en observant que $\operatorname {d} 'K=0,$ les équations (12) et (13) pourront se mettre sous la forme

$\nu -\int {\frac {\delta p}{\Delta }}=\operatorname {d} '(\phi +K)-{\frac {1}{2}}\left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right),$

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\delta (\phi +k),$

et par conséquent, en mettant $\phi ,$ au lieu de $\phi +k,$

$\nu -\int {\frac {\delta p}{\Delta }}=\operatorname {d} '\phi -{\frac {1}{2}}\left(\operatorname {d} 'x^{2}+\operatorname {d} 'y^{2}+\operatorname {d} 'z^{2}\right),$

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\delta \phi ,$

qui sont absolument les mêmes que si l’en avait fait $H=0.$ Nous aurons donc alors

$P=0,\quad Q=0,\quad R=0\,;$

et par suite, les équations (17), (18), (19), (20) du n.^o  précédent, deviendront

$\operatorname {d} 'x={\frac {\operatorname {d} \phi }{\operatorname {d} x}},\quad \operatorname {d} 'y={\frac {\operatorname {d} \phi }{\operatorname {d} y}},\quad \operatorname {d} 'z={\frac {\operatorname {d} \phi }{\operatorname {d} z}}\,;\qquad$(21)

$\nu -\int {\frac {\delta p}{\Delta }}={\frac {\operatorname {d} \phi }{\operatorname {d} t}}+{\frac {1}{2}}\left\{\left({\frac {\operatorname {d} \phi }{\operatorname {d} x}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} y}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} z}}\right)^{2}\right\},\qquad$(22)

$0={\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} x}}{\frac {\operatorname {d} \phi }{\operatorname {d} x}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} y}}{\frac {\operatorname {d} \phi }{\operatorname {d} y}}+{\frac {\operatorname {d} \Delta }{\operatorname {d} z}}{\frac {\operatorname {d} \phi }{\operatorname {d} z}}+\Delta \left({\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z^{2}}}\right)$ (23)

$0={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y^{2}}}+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z^{2}}},\qquad$(24)

on peut observer de plus que le quarré de la vitesse absolue du fluide est

$\left({\frac {\operatorname {d} \phi }{\operatorname {d} x}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} y}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} z}}\right)^{2}.$

15. Lorsque les mouvemens du fluide sont très-petits, on peut négliger les termes affectés des puissances, et produits de puissances des vitesses $\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z.$ Dans ce cas les équations (3) deviennent

${\begin{aligned}{\frac {\operatorname {d} \nu }{\operatorname {d} x}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} x}}&={\frac {\operatorname {d} 'x}{\operatorname {d} t}},\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} y}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} y}}&={\frac {\operatorname {d} 'y}{\operatorname {d} t}},\\\\{\frac {\operatorname {d} \nu }{\operatorname {d} z}}-{\frac {1}{\Delta }}{\frac {\operatorname {d} p}{\operatorname {d} z}}&={\frac {\operatorname {d} 'z}{\operatorname {d} t}}\,;\end{aligned}}$

desquelles nous tirerons

$\delta \nu -{\frac {\delta p}{\Delta }}={\frac {\operatorname {d} 'x}{\operatorname {d} t}}\delta x+{\frac {\operatorname {d} 'y}{\operatorname {d} t}}\delta y+{\frac {\operatorname {d} 'z}{\operatorname {d} t}}\delta z={\frac {\operatorname {d} \psi }{\operatorname {d} t}}.$

C’est de cette dernière que les divers auteurs qui ont écrit sur l’hydrodynamique ont cru pouvoir conclure immédiatement que $\psi$ ou $\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z$ devait être une variation exacte ; conclusion dont nous avons déjà relevé l’inexactitude dans le n.^o 6. Nous allons encore le faire ici d’une autre manière.

Si l’on détermine $\phi$ au moyen de l’équation

$\nu -\int {\frac {\delta p}{\Delta }}={\frac {\operatorname {d} \phi }{\operatorname {d} t}},\qquad$(25)

ce qui est toujours possible, tout au moins par les quadratures, et ce qui donne

$\delta \nu -{\frac {\delta p}{\Delta }}=\delta {\frac {\operatorname {d} \phi }{\operatorname {d} t}},$

de laquelle on tirera

${\frac {\operatorname {dd} 'x}{\operatorname {d} t}}\delta x+{\frac {\operatorname {dd} 'y}{\operatorname {d} t}}\delta y+{\frac {\operatorname {dd} 'z}{\operatorname {d} t}}\delta z=\delta {\frac {\operatorname {d} \phi }{\operatorname {d} t}},$

et par suite

${\frac {\operatorname {dd} 'x}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} x}},\quad {\frac {\operatorname {dd} 'y}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} y}},\quad {\frac {\operatorname {dd} 'z}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} z}}\,;$

et de là, en intégrant par rapport à $t$

$\operatorname {d} 'x={\frac {\operatorname {d} \phi }{\operatorname {d} x}}+P,\quad \operatorname {d} 'y={\frac {\operatorname {d} \phi }{\operatorname {d} y}}+Q,\quad \operatorname {d} 'z={\frac {\operatorname {d} \phi }{\operatorname {d} z}}+R.\qquad$(26)

chacune des quantités $P,Q,R,$ étant une fonction arbitraire de $x,y,z,$ sans $t.$ Ce cas rentre donc dans celui que nous avons déjà considéré dans le n.^o 13.

Au moyen des valeurs de $\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z$ que nous venons de donner, on trouvera

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\delta \phi +P\delta x+Q\delta y+R\delta z,\qquad$(27)

d’où l’on voit que, pour que le premier membre de cette équation fût une variation exacte, il faudrait qu’il en fût de même de $P\delta x+Q\delta y+R\delta z\,;$ ce qui peut fort bien ne pas être, puisque $P,Q,R$ sont des fonctions entièrement arbitraires de $x,y,z.$

Si, par exemple, on fait

$\operatorname {d} 'x={\frac {\operatorname {d} \phi }{\operatorname {d} x}}+ny,\quad \operatorname {d} 'y={\frac {\operatorname {d} \phi }{\operatorname {d} y}}-nx,\quad \operatorname {d} 'z={\frac {\operatorname {d} \phi }{\operatorname {d} z}},$

$n$ étant un, nombre constant quelconque, on aura

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z=\delta \phi +n(y\delta x-x\delta y),$

dont le premier membre ne saurait être une variation exacte, tant que $n$ n’est point nul, quoique les mêmes valeurs donnent

${\frac {\operatorname {dd} 'x}{\operatorname {d} t}}\delta x+{\frac {\operatorname {dd} 'y}{\operatorname {d} t}}\delta y+{\frac {\operatorname {dd} 'z}{\operatorname {d} t}}\delta z=\delta {\frac {\operatorname {d} \phi }{\operatorname {d} t}}.$

16. Il semblerait résulter de ce que nous venons de dire, dans le dernier n.^o, que les résultats des recherches des grands géomètres de nos jours, sur les petites occillations des fluides, ont beaucoup moins de généralité qu’ils ne l’ont cru et annoncé. Heureusement il n’en est pas ainsi ; $\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z$ peut être réputée une variation exacte, toutes les fois que les molécules du fluide ne font que de très-petites occillations autour de leur position d’équilibre. Pour le prouver, faisons, en général,

$x=a+x_{1},\qquad y=b+y_{1},\qquad z=c+z_{1},$

en désignant par $x_{1},y_{1},z_{1},$ des fonctions de $x,y,z,$ nulles en même temps que $t$ et convenablement déterminées, dès-lors nous aurons

${\begin{aligned}&\operatorname {d} 'x=\operatorname {d} 'x_{1}={\frac {\operatorname {d} x}{\operatorname {d} t}}+{\frac {\operatorname {d} x_{1}}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} x_{1}}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} x_{1}}{\operatorname {d} z}}\operatorname {d} 'z,\\\\&\operatorname {d} 'y=\operatorname {d} 'y_{1}={\frac {\operatorname {d} y}{\operatorname {d} t}}+{\frac {\operatorname {d} y_{1}}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} y_{1}}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} y_{1}}{\operatorname {d} z}}\operatorname {d} 'z,\\\\&\operatorname {d} 'z=\operatorname {d} 'z_{1}={\frac {\operatorname {d} z}{\operatorname {d} t}}+{\frac {\operatorname {d} z_{1}}{\operatorname {d} x}}\operatorname {d} 'x+{\frac {\operatorname {d} z_{1}}{\operatorname {d} y}}\operatorname {d} 'y+{\frac {\operatorname {d} z_{1}}{\operatorname {d} z}}\operatorname {d} 'z,\end{aligned}}$

dans le cas actuel $x_{1},y_{1},z_{1},$ étant très-petits, nous pouvons négliger les termes de plus d’une dimension par rapport à ces quantités ; nous aurons donc

$\operatorname {d} 'x={\frac {\operatorname {d} x_{1}}{\operatorname {d} t}},\quad \operatorname {d} 'y={\frac {\operatorname {d} y_{1}}{\operatorname {d} t}},\quad \operatorname {d} 'z={\frac {\operatorname {d} z_{1}}{\operatorname {d} t}},$

et par suite, au moyen des équations (26)

${\frac {\operatorname {d} x_{1}}{\operatorname {d} t}}={\frac {\operatorname {d} \phi }{\operatorname {d} x}}+P,\quad {\frac {\operatorname {d} y_{1}}{\operatorname {d} t}}={\frac {\operatorname {d} \phi }{\operatorname {d} y}}+Q,\quad {\frac {\operatorname {d} z_{1}}{\operatorname {d} t}}={\frac {\operatorname {d} \phi }{\operatorname {d} z}}+R.$

d’où l’on tire, par l’élimination successive de $R,Q$ et $P,$

${\begin{aligned}&{\frac {\operatorname {d} ^{2}x_{1}}{\operatorname {d} t\operatorname {d} y}}-{\frac {\operatorname {d} ^{2}y_{1}}{\operatorname {d} t\operatorname {d} x}}={\frac {\operatorname {d} P}{\operatorname {d} y}}-{\frac {\operatorname {d} Q}{\operatorname {d} x}},\\\\&{\frac {\operatorname {d} ^{2}y_{1}}{\operatorname {d} t\operatorname {d} z}}-{\frac {\operatorname {d} ^{2}z_{1}}{\operatorname {d} t\operatorname {d} y}}={\frac {\operatorname {d} Q}{\operatorname {d} z}}-{\frac {\operatorname {d} R}{\operatorname {d} y}},\\\\&{\frac {\operatorname {d} ^{2}z_{1}}{\operatorname {d} t\operatorname {d} x}}-{\frac {\operatorname {d} ^{2}x_{1}}{\operatorname {d} t\operatorname {d} z}}={\frac {\operatorname {d} R}{\operatorname {d} x}}-{\frac {\operatorname {d} P}{\operatorname {d} z}}\,;\end{aligned}}$

et de là en intégrant par rapport à $t,$

$\left.{\begin{aligned}&{\frac {\operatorname {d} x_{1}}{\operatorname {d} y}}-{\frac {\operatorname {d} y_{1}}{\operatorname {d} x}}=t\left({\frac {\operatorname {d} P}{\operatorname {d} y}}-{\frac {\operatorname {d} Q}{\operatorname {d} x}}\right),\\\\&{\frac {\operatorname {d} y_{1}}{\operatorname {d} z}}-{\frac {\operatorname {d} z_{1}}{\operatorname {d} y}}=t\left({\frac {\operatorname {d} Q}{\operatorname {d} z}}-{\frac {\operatorname {d} R}{\operatorname {d} y}}\right),\\\\&{\frac {\operatorname {d} z_{1}}{\operatorname {d} x}}-{\frac {\operatorname {d} x_{1}}{\operatorname {d} z}}=t\left({\frac {\operatorname {d} R}{\operatorname {d} x}}-{\frac {\operatorname {d} P}{\operatorname {d} z}}\right).\end{aligned}}\right\}\quad (28)$

Nous n’ajoutons point de fonctions arbitraires de $x,y,z,$ parce que $x_{1},y_{1},z_{1},$ et par suite les premiers membres de ces équations, doivent être nuls en même temps que $t.$

Présentement, puisque $x_{1},y_{1},z_{1},$ doivent rester renfermés dans des limites très-resserrées, les premiers membres des équations (28) doivent jouir de la même propriété, ce qui ne pourra être qu’autant qu’on aura identiquement

${\frac {\operatorname {d} P}{\operatorname {d} y}}-{\frac {\operatorname {d} Q}{\operatorname {d} x}}=0,\quad {\frac {\operatorname {d} Q}{\operatorname {d} z}}-{\frac {\operatorname {d} R}{\operatorname {d} y}}=0,\quad {\frac {\operatorname {d} R}{\operatorname {d} x}}-{\frac {\operatorname {d} P}{\operatorname {d} z}}=0\,;$

qui sont précisément les équations nécessaires et suffisantes pour que

$P\delta x+Q\delta y+R\delta z,$

et par suite, en vertu de l’équation (27),

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z,$

soit une variation exacte, ce qui prouve ce que nous avons avancé au commencement du présent n.^o.

On voit même par là que, bien que les mouvemens soient très-petits, si la quantité

$\operatorname {d} 'x\delta x+\operatorname {d} 'y\delta y+\operatorname {d} 'z\delta z,$

n’est pas une variation exacte, le fluide éprouvera nécessairement un mouvement de translation ; mais il ne faudrait pas en conclure que la réciproque est également vraie.

17. Dans les applications usuelles de la théorie du mouvement des fluides, ce mouvement peut être regardé comme uniforme ; c’est-à-dire que les vitesses $\operatorname {d} 'x,\operatorname {d} 'y,\operatorname {d} 'z,$ sont fonctions de $x,y,z$ seulement ; de sorte que l’on a

${\frac {\operatorname {dd} 'x}{\operatorname {d} t}}=0,\quad {\frac {\operatorname {dd} 'y}{\operatorname {d} t}}=0,\quad {\frac {\operatorname {dd} 'z}{\operatorname {d} t}}=0.$

Si donc nous supposons que ce soit le cas traité dans le n.^o 14 ; nous aurons

${\begin{aligned}&{\frac {\operatorname {dd} 'x}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} x}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x\operatorname {d} t}}=0,\\\\&{\frac {\operatorname {dd} 'y}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} y}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y\operatorname {d} t}}=0,\\\\&{\frac {\operatorname {dd} 'z}{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} z}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z\operatorname {d} t}}=0\,;\end{aligned}}$

et par suite

$\delta {\frac {\operatorname {d} \phi }{\operatorname {d} t}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x\operatorname {d} t}}\delta x+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y\operatorname {d} t}}\delta y+{\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z\operatorname {d} t}}\delta z=0\,;$

d’où on conclura

${\frac {\operatorname {d} \phi }{\operatorname {d} t}}=T,$

en désignant par $T$ une fonction arbitraire de $t,$ qui doit être la même pour toutes les molécules du fluide. Faisant de plus, pour abréger,

$\left({\frac {\operatorname {d} \phi }{\operatorname {d} x}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} y}}\right)^{2}+\left({\frac {\operatorname {d} \phi }{\operatorname {d} z}}\right)^{2}=u^{2},$

$\int {\frac {\delta p}{\Delta }}=\gamma ,$

l’équation (22) deviendra

$\nu -\gamma =T+{\frac {u^{2}}{2}}$

En marquant de l’indice $1$ les diverses quantités qui, pour la même valeur $t$ du temps, sont relatives à une autre molécule quelconque, nous aurons de même

$\nu _{1}-\gamma _{1}=T+{\frac {u_{1}^{2}}{2}},$

et, par suite

$\left(\nu -\nu _{1}\right)-\left(\gamma -\gamma _{1}\right)={\frac {1}{2}}\left(u^{2}-u_{1}^{2}\right)\,;\qquad$(29)

équation qui établit une relation importante entre les vitesses absolues et simultanées $u$ et $u_{1}$ et les pressions correspondantes $p,p_{1},$ de deux molécules quelconques du fluide en mouvement.

Avant d’aller plus loin, nous ferons observer que, si le fluide se divise et éprouve un changement brusque aux points de division, dès-lors on n’a plus, pour ces points, les équations

${\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} x}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} x\operatorname {d} t}},\quad {\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} y}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} y\operatorname {d} t}},\quad {\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} t\operatorname {d} z}}={\frac {\operatorname {d} ^{2}\phi }{\operatorname {d} z\operatorname {d} t}}\,;$

sur lesquelles nous nous sommes fondés pour parvenir à l’équation (29). Si donc on veut que cette équation soit vraie dans tous les cas, il faut que les molécules auxquelles elle est relative soient prises toutes les deux dans la partie du fluide qui n’a point éprouvé de division.

18. Le plus souvent les forces accélératrices qui sollicitent le fluide en mouvement se réduisent à la pesanteur que l’on regarde comme une force constante et de direction parallèle à une droite fixe, dans toute l’étendue de la masse fluide. Représentant cette force par $g$ et prenant la coordonnée $z$ verticale, et dirigée de haut en bas, nous aurons

$\nu =gz,$

dont la substitution dans la formule (29) donnera

$g\left(z-z_{1}\right)-\left(\gamma -\gamma _{1}\right)={\frac {1}{2}}\left(u^{2}-u_{1}^{2}\right).$

Dans le cas on le fluide est homogène et incompressible, la densité $\Delta$ est absolument constante, et l’on a, par conséquent,

$\gamma =\int {\frac {\delta p}{\Delta }}={\frac {p}{\Delta }},$

en vertu de quoi la précédente équation devient

$y\left(z-z_{1}\right)-\left({\frac {p-p_{1}}{\Delta }}\right)={\frac {1}{2}}\left(u^{2}-u_{1}^{2}\right).\qquad$(30)

Si, pour tous les points de la surface libre du fluide, la pression doit être constante, on aura pour les divers points de cette surface

$g\left(z-z_{1}\right)={\frac {1}{2}}\left(u^{2}-u_{1}^{2}\right),$

et si les vitesses $u,u_{1}$ sont assez petites pour qu’on puisse en négliger les quarrés,

$g\left(z-z_{1}\right)=0\,;$

c’est-à-dire qu’alors cette surface est plane et horizontale. Il en serait encore de même, si tous les points de cette surface étaient doués de la même vitesse.

Il ne faut point perdre de vue, au surplus, que ces conclusions supposent que le mouvement du fluide est uniforme.

19. Lorsqu’un fluide pesant, homogène et incompressible, de l’eau, par exemple, est en équilibre dans un vase invariable de forme et de position. Si l’on pratique, dans la paroi de ce vase, une ouverture par laquelle le fluide puisse s’écouler, le mouvement, d’abord nul, augmentera avec rapidité ; et, si l’orifice d’écoulement est très-petit, après un intervalle de temps très-court, ce mouvement sera à très-peu près uniforme, et l’écoulement pourra conséquemment se calculer par la formule (30), savoir :

$g\left(z-z_{1}\right)-\left({\frac {p-p_{1}}{\Delta }}\right)={\frac {1}{2}}\left(u^{2}-u_{1}^{2}\right).$

À la surface supérieure du fluide, la vitesse $u$ est extrêmement petite, et l’on pourra en négliger le quarré. Si donc on suppose que $z_{1},p_{1},u_{1}$ soient relatifs à cette surface, l’équation précédente deviendra

$g\left(z-z_{1}\right)-\left({\frac {p-p_{1}}{\Delta }}\right)={\frac {u^{2}}{2}}.\qquad$(31)

Ordinairement, la pression est la même à l’orifice d’écoulement qu’à la surface du fluide. Si donc on veut que $z,p,u,$ appartiennent à un point quelconque de cet orifice, on aura $p=p_{1},$ et par suite

$g\left(z-z_{1}\right)={\frac {u^{2}}{2}}\,;\qquad$(32)

formule connue depuis très-long-temps ; mais à laquelle on n’était parvenu jusqu’ici qu’à l’aide d’hypothèses plus ou moins précaires.

Si le fluide éprouvait une plus grande pression à son niveau supérieur qu’à l’orifice d’écoulement, en désignant par $c$ cet excès de pression, on aurait, pour calculer la vitesse du fluide à l’orifice, la formule

$g\left(z-z_{1}\right)+{\frac {c}{\Delta }}={\frac {u^{2}}{2}}.$

20. La formule (31), mise sous la forme

$p=p_{1}+g\Delta \left(z-z_{1}\right)-{\frac {\Delta u^{2}}{2}},$

peut servir à calculer la pression que le fluide exerce contre un élément quelconque des parois du vase qui le renferme, ou de tout autre corps auquel le fluide serait contigu.

Supposons, par exemple, que l’on oppose perpendiculairement au choc de la veine fluide une figure plane quelconque, la vitesse $u$ du fluide contigu sera nulle, ou au moins très-petite : on aura donc

$p=p_{1}+g\Delta \left(z-z_{1}\right),$

pour la pression qu’éprouve chaque élément de la face de cette figure qui se trouve directement exposée au choc du fluide. De plus, si l’écoulement a lieu sous la pression ordinaire de l’air, l’autre face éprouvera une pression $p_{1},$ dont l’effet sera contraire à celui de la pression $p\,;$ il restera donc

$p-p_{1}=g\Delta \left(z-z_{1}\right),$

pour l’impulsion que reçoit, de la part du fluide, chaque élément de la figure en question.

21. Jusqu’ici nous n’avons pour ainsi dire considéré que les équations (1) et celles qui en dérivent. Nous allons présentement nous occuper de quelques résultats que l’on peut déduire des équations qui naissent de la condition de continuité du fluide, considérées isolément.

Reprenons l’équation (5). En multipliant tous ses termes par $\operatorname {d} x\operatorname {d} y\operatorname {d} z$ et intégrant, nous aurons

$0=\iiint \left({\frac {\operatorname {d} \Delta }{\operatorname {d} t}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'y}{\operatorname {d} y}}+{\frac {\operatorname {d} .\Delta \operatorname {d} 'z}{\operatorname {d} z}}\right)\operatorname {d} x\operatorname {d} y\operatorname {d} z,$

quelles que puissent être d’ailleurs les limites entre lesquelles il faille prendre les intégrales.

Présentement le terme

$\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}\operatorname {d} x\operatorname {d} y\operatorname {d} z,$

que renferme cette équation, intégré par rapport à $x,$ donne pour résultat

$\iint \Delta \operatorname {d} 'x\operatorname {d} y\operatorname {d} z\,;$

mais il reste à déterminer quelles sont les parties de ce résultat qui doivent êtres prises et celles qui doivent être laissées. Supposons pour cela que l’intégrale

$\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}\operatorname {d} x\operatorname {d} y\operatorname {d} z$

soit prise dans toute l’étendue d’une partie finie quelconque du fluide en mouvement, limitée par une surface rentrante que, pour la commodité du langage, nous désignerons par (O). L’équation de cette surface donnera pour $x,$ en fonction de $y$ et $z,$ différentes valeurs, les unes imaginaires, dont nous ferons abstraction, et les autres réelles, que nous représenterons par $x_{1},x_{2},x_{3},x_{4},\ldots \,;$ et que nous supposerons rangées dans un ordre tel que l’on ait

$x_{2}-x_{1}>0,\quad x_{3}-x_{2}>0,\quad x_{4}-x_{3}>0,\ldots \,;$

ce qui est toujours possible. Alors, toutes les molécules qui auront les mêmes coordonnées $y,z,$ mais pour lesquelles $x$ sera compris entre

$-{\frac {1}{0}}$ et $x_{1},\ x_{2}$ et $x_{3},\ x_{4}$ et $x_{3},\ldots$

seront nécessairement situées hors de la surface (O), tandis que celles pour lesquelles $x$ sera compris entre

$x_{1}$ et $x_{2},\ x_{3}$ et $x_{4},\ldots$

seront situées dans l’intérieur de cette même surface ; d’où l’on conclura qu’il ne faut prendre de l’intégrale

$\int \Delta \operatorname {d} 'x\operatorname {d} y\operatorname {d} z,$

que les parties qui sont comprises entre $x_{1}$ et $x_{2},\ x_{3}$ et $x_{4},\ldots$ et que, par suite, on a

$\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}\operatorname {d} x\operatorname {d} y\operatorname {d} z=\iint \left\{(\Delta \operatorname {d} 'x)_{2}-(\Delta \operatorname {d} 'x)_{1}\right\}\operatorname {d} y\operatorname {d} z$

$+\iint \left\{(\Delta \operatorname {d} 'x)_{4}-(\Delta \operatorname {d} 'x)_{3}\right\}\operatorname {d} y\operatorname {d} z+\ldots$

en représentant par $(\Delta \operatorname {d} 'x)_{1},(\Delta \operatorname {d} 'x)_{2},(\Delta \operatorname {d} 'x)_{3},\ldots$ les valeurs de $\Delta \operatorname {d} 'x$ qui correspondent à celles $x_{1},x_{2},x_{3},\ldots$ de $x.$

Le second membre de cette dernière équation peut être mis sous une forme beaucoup plus simple et plus facile à interpréter. Pour cela, par un point quelconque $(x,y,z)$ de la surface (O), menons-lui une normale. Prenons sur cette normale, dans l’intérieur de la surface (O), et à la distance $r$ de son pied, un second point $(\alpha ,\beta ,\gamma )$ et nous aurons

$\pm \left({\frac {x-\alpha }{r}}\right),$

pour la valeur du cosinus de l’angle que fait la normale avec l’axe des $x,$ angle qui est le même que celui que fait le plan tangent avec le plan des $yz$. Si donc nous représentons par $\operatorname {d} s^{2}$ l’élément de la surface (O) dont les coordonnées sont $x,y,z,$ et dont la projection sur le plan des $yz$ est $\operatorname {d} y\operatorname {d} z,$ nous aurons

$\operatorname {d} y\operatorname {d} z=\pm \left({\frac {x-\alpha }{r}}\right)\operatorname {d} s^{2},$

en observant de prendre le signe supérieur ou le signe inférieur, suivant que $x-\alpha$ sera positif ou négatif. Alors, en affectant des indices $1,2,3,\ldots$ les valeurs de

$\left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)$

qui correspondent aux valeurs $x_{1},x_{2},x_{3},\ldots ,$ nous aurons

${\begin{aligned}&(\Delta \operatorname {d} 'x)_{1}\operatorname {d} y\operatorname {d} z=-\left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)_{1},\\\\&(\Delta \operatorname {d} 'x)_{2}\operatorname {d} y\operatorname {d} z=+\left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)_{2},\end{aligned}}$

${\begin{aligned}&(\Delta \operatorname {d} 'x)_{3}\operatorname {d} y\operatorname {d} z=-\left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)_{3},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{aligned}}$

et par suite

$\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}\operatorname {d} x\operatorname {d} y\operatorname {d} z=$

$\iint \left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)_{1}+\iint \left(\Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\right)_{2}+\ldots$

ou, plus simplement

$\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}\operatorname {d} x\operatorname {d} y\operatorname {d} z=\iint \Delta \operatorname {d} 'x{\frac {x-\alpha }{r}}.\operatorname {d} s^{2}\,;$

pourvu que l’intégrale du second membre soit prise dans toute l’étendue de la surface (O).

Par des procédés entièrement analogues, on trouverait qu’entre les mêmes limites on a

${\begin{aligned}&\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'y}{\operatorname {d} y}}\operatorname {d} x\operatorname {d} y\operatorname {d} z=\iint \Delta \operatorname {d} 'y{\frac {y-\beta }{r}}.\operatorname {d} s^{2},\\\\&\iiint {\frac {\operatorname {d} .\Delta \operatorname {d} 'z}{\operatorname {d} z}}\operatorname {d} x\operatorname {d} y\operatorname {d} z=\iint \Delta \operatorname {d} 'z{\frac {z-\gamma }{r}}.\operatorname {d} s^{2}\,;\end{aligned}}$

substituant ces valeurs dans l’équation identique

$0=\iiint \left(\Delta +{\frac {\operatorname {d} .\Delta \operatorname {d} 'x}{\operatorname {d} x}}+\Delta +{\frac {\operatorname {d} .\Delta \operatorname {d} 'y}{\operatorname {d} y}}+\Delta +{\frac {\operatorname {d} .\Delta \operatorname {d} 'z}{\operatorname {d} z}}\right)\operatorname {d} x\operatorname {d} y\operatorname {d} z,$

nous trouverons définitivement

$0=\iiint {\frac {\operatorname {d} .\Delta }{\operatorname {d} t}}\operatorname {d} x\operatorname {d} y\operatorname {d} z$

$+\iint \left({\frac {x-\alpha }{r}}\operatorname {d} 'x+{\frac {y-\beta }{r}}\operatorname {d} 'y+{\frac {z-\gamma }{r}}\operatorname {d} 'z\right)\Delta \operatorname {d} s^{2},\qquad$(33)

ou plus simplement

$0=\iiint {\frac {\operatorname {d} .\Delta }{\operatorname {d} t}}\operatorname {d} x\operatorname {d} y\operatorname {d} z+\iint \Delta \operatorname {d} 'r\operatorname {d} s^{2},\qquad$(34)

en faisant, pour abréger,

$\operatorname {d} 'r={\frac {x-\alpha }{r}}\operatorname {d} 'x+{\frac {y-\beta }{r}}\operatorname {d} 'y+{\frac {z-\gamma }{r}}\operatorname {d} 'z\,;$

équation qui n’est que la différentielle de cette autre équation identique

$r^{2}=(x-\alpha )^{2}+(y-\beta )^{2}+(z-\gamma )^{2},$

prise en regardant $\alpha ,\beta ,\gamma$ comme constantes ; ce qui fait voir que $\operatorname {d} 'r$ exprime la vitesse du fluide décomposée suivant la normale $r\,;$ vitesse qui doit être regardée comme positive, lorsqu’elle tend à pousser contre la surface (O) le fluide situé dans l’intérieur de cette surface, tandis qu’elle sera réputée négative dans le cas contraire.

Si le fluide est homogène et incompressible, la densité $\Delta$ sera constante, et l’équation (34) deviendra

$0=\iint \operatorname {d} 'r\operatorname {d} s^{2},\qquad$(35)

qui serait encore vraie quand même le fluide serait simplement incompressible sans être homogène, comme on peut le déduire directement de l’équation (8), traitée de la même manière que nous venons de traiter l’équation (6).

Les équations (34) et (35) peuvent être utiles dans beaucoup de cas, notamment s’il s’agit de calculer les oscillations d’un fluide renfermé dans un vase.

22, Considérons, dans l’intérieur du fluide en mouvement, une ligne courbe rentrante quelconque $\mathrm {(P)}$ et une surface $\mathrm {(R),}$ assujettie aux seules conditions d’être limitée par cette ligne courbe, et de ne pas présenter de déchirures dans l’intervalle. Représentons toujours par $\operatorname {d} 'r$ la vitesse d’une molécule quelconque contiguë à cette surface, suivant la normale $r\,;$ regardant cette vitesse comme positive, lorsqu’elle tend à pousser le fluide contre la surface, et comme négative dans le cas contraire.

Cela posé, chaque élément de la surface $\mathrm {(R)}$ est contigu à deux molécules différentes, pour lesquelles la valeur absolue de $\operatorname {d} 'r$ sera identiquement la même. Quant au signe, il doit être différent, attendu que, quand une des deux molécules est poussée contre la surface, l’autre tend à s’en détacher. Si donc on prend l’intégrale

$\iint \operatorname {d} 'r\operatorname {d} s^{2}$

relativement aux molécules qui se trouvent situées d’un seul et même côté de la surface dans toute son étendue, la valeur absolue du résultat sera le même quel que soit le côté de cette surface qu’on aura choisi ; de sorte qu’en représentant cette valeur absolue par

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2},$

nous aurons

$\iint \operatorname {d} 'r\operatorname {d} s^{2}=\pm \iint (R)\operatorname {d} 'r\operatorname {d} s^{2},$

le signe $+$ ou $-$ devant être déterminé d’une manière convenable.

Présentement, soit $\mathrm {(R')}$ une autre surface, assujettie aux mêmes conditions que $\mathrm {(R)} \,;$ alors les surfaces $\mathrm {(R),(R')}$ formeront, par leur réunion, une troisième surface rentrante, pour laquelle nous aurons

$\iint \operatorname {d} 'r\operatorname {d} s^{2}=0.$

Si, comme nous le supposerons dans tout ce qui suit, le fluide est incompressible, l’intégrale qui compose le premier membre de cette dernière équation se décompose en deux parties, dont l’une relative à la surface $\mathrm {(R)}$ et l’autre relative à la surface $\mathrm {(R').}$ Mais, si les deux surfaces n’ont d’autres points communs que ceux de la courbe $\mathrm {(P)} ,$ tout le fluide intérieur sera situé d’un même côté de chacune de les surfaces, et les parties de l’intégrale

$\iint \operatorname {d} 'r\operatorname {d} s^{2},$

qui leur sont relatives seront exprimées, la première par

$\pm \iint (\mathrm {R} )\operatorname {d} 'r\operatorname {d} s^{2},$

et la seconde par

$\pm \iint (\mathrm {R} ')\operatorname {d} 'r\operatorname {d} s^{2},$

d’où nous conclurons que l’on doit avoir

$\iint (\mathrm {R} )\operatorname {d} 'r\operatorname {d} s^{2}=\iint (\mathrm {R} ')\operatorname {d} 'r\operatorname {d} s^{2}.\qquad$(36)

Si les surfaces $\mathrm {(R),(R')}$ se coupent en des points autres que ceux de la courbe $\mathrm {(P),}$ le fluide intérieur peut n’être pas situé en totalité, d’un même côté de chacune de ces surfaces ; alors le raisonnement qui nous a conduit à l’équation (36) cesse d’être applicable, et cette équation semble ne plus devoir être vraie. Cependant il est toujours possible de trouver une troisième surface $\mathrm {(R''),}$ assujettie aux mêmes conditions que $\mathrm {(R)}$ et $\mathrm {(R'),}$ et qui de plus ne les rencontre que suivant la courbe $\mathrm {(P).}$ Alors, on aura, par ce qui précède,

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2}=\iint (R'')\operatorname {d} 'r\operatorname {d} s^{2},$

$\iint (\mathrm {R'} )\operatorname {d} 'r\operatorname {d} s^{2}=\iint (R'')\operatorname {d} 'r\operatorname {d} s^{2},$

et par suite

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2}=\iint (R')\operatorname {d} 'r\operatorname {d} s^{2}\,;$

comme si les surfaces $\mathrm {(R),(R')}$ ne se rencontraient que suivant la courbe $\mathrm {(P).}$

On pourrait conclure de là que l’expression

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2}$

est entièrement indépendante de la nature de la surface $\mathrm {(R),}$ et quelle dépend seulement de la nature de la courbe $\mathrm {(P),}$ qui lui sert de limite et du temps $t.$

23. Pour donner une application de ce qui précède, supposons que le fluide soit contenu dans un vase invariable de forme et de position. Supposons encore que la courbe $\mathrm {(P),}$ tracée sur les parois du vase en fasse le tour, et qu’enfin la surface $\mathrm {(R)}$ se compose, 1.^o d’une partie plus ou moins étendue de ces parois ; 2.^o d’une section de surface quelconque, limitée par ces mêmes parois ; ce qui laisse la forme et la position de cette surface entièrement indépendante de la courbe $\mathrm {(P)}$ et de la forme du vase. En désignant par $\mathrm {(A)}$ la surface dont il s’agit, nous aurons

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2}=\iint (A)\operatorname {d} 'r\operatorname {d} s^{2}\,;\qquad$(37)

en observant que, pour toutes les molécules qui sont contiguës aux parois du vase, l’on a $\operatorname {d} 'r=0,$ et que, par suite, la partie de l’intégrale

$\iint (R)\operatorname {d} 'r\operatorname {d} s^{2},$

qui est relative à ces parois est identiquement nulle.

Présentement, d’après la conclusion qui termine le n.^o 22, le premier membre de l’équation (37) ne dépend que du temps $t$ et tout au plus de la courbe $\mathrm {(P)} \,;$ et, comme la section $\mathrm {(A)}$ est entièrement indépendante de cette courbe, nous en conclurons que l’expression

$\iint (A)\operatorname {d} 'r\operatorname {d} s^{2}$

ne dépend absolument que du temps $t$ et de la forme du vase, et nullement de la forme ou de la position de la surface à laquelle appartient la section $\mathrm {(A).}$

Si, pour fixer les idées, on suppose que cette section soit plane, et parallèle au plan des $xy,$ on aura

$\operatorname {d} 'r=\operatorname {d} 'z\quad$et$\quad \operatorname {d} s^{2}=\operatorname {d} x\operatorname {d} y,$

et par suite

$\iint \operatorname {d} 'z\operatorname {d} x\operatorname {d} y=\theta ,$

en désignant par $\theta$ une fonction de $t$ convenablement déterminée, mais qui doit être indépendante et de la position du plan coupant et de la direction des axes des coordonnées.

Supposons que le vase qui renferme le fluide soit très-étroit ; représentons par $x_{0},y_{0}$ les valeurs de $x,y$ pour un point quelconque de la section, et soit fait

$x=x_{0}+ih,\qquad y=y_{0}+ik,$

et par suite

$\operatorname {d} x=i\operatorname {d} h,\qquad \operatorname {d} y=i\operatorname {d} k,$

${\begin{aligned}\operatorname {d} 'z=(\operatorname {d} 'z)_{0}&+{\frac {i}{1}}\left\{h\left({\frac {\operatorname {dd} 'z}{\operatorname {d} x}}\right)_{0}+k\left({\frac {\operatorname {dd} 'z}{\operatorname {d} y}}\right)_{0}\right\}\\\\&+{\frac {i^{2}}{1.2}}\left\{h^{2}\left({\frac {\operatorname {d} ^{2}\operatorname {d} 'z}{\operatorname {d} x^{2}}}\right)_{0}+2hk\left({\frac {\operatorname {d} ^{2}\operatorname {d} 'z}{\operatorname {d} x\operatorname {d} y}}\right)_{0}+k^{2}\left({\frac {\operatorname {d} ^{2}\operatorname {d} 'z}{\operatorname {d} y^{2}}}\right)_{0}\right\}\\&+\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{aligned}}$

en dénotant respectivement par

$(\operatorname {d} 'z)_{0},\quad \left({\frac {\operatorname {dd} 'z}{\operatorname {d} x}}\right)_{0},\quad \left({\frac {\operatorname {dd} 'z}{\operatorname {d} y}}\right)_{0},\ldots$

ce que deviennent

$\operatorname {d} 'z,\quad {\frac {\operatorname {dd} 'z}{\operatorname {d} x}},\quad {\frac {\operatorname {dd} 'z}{\operatorname {d} y}},\ldots$

qui répondent au point $\left(x_{0},y_{0},z_{0}\right).$

Substituant ces valeurs dans l’équation (38), et négligeant les termes de plus de deux dimensions en $i,$ nous aurons

$(\operatorname {d} 'z)i^{2}\iint \operatorname {d} h\operatorname {d} k=\theta ,$

ou encore

$\varpi \operatorname {d} 'z=\theta ,\qquad$(39)

en observant que $i^{2}\iint \operatorname {d} h\operatorname {d} k$ exprime l’aire de la section que l’on considère et représentant cette aire par $\varpi .$

On voit par là que l’hypothèse du parallélisme des tranches peut être employée comme moyen approximatif ; et telle en est je crois la première démonstration générale. L’auteur de la mécanique analitique était déjà parvenu au même résultat ; mais seulement pour le cas où le vase, à très-peu près vertical, n’aurait que deux dimensions, et les procédés qu’il a mis en usage, pour parvenir à son but, seraient, à raison de leur complication, à peu près impraticables pour toute autre forme ou position du vase.

Nous devrions peut-être terminer cet essai par quelques applications particulières des formules qu’il a pour objet d’obtenir ; mais comme les seules applications qui ne soient point au-dessus de notre portée ont été déjà données par divers géomètres, beaucoup mieux que nous ne pourrions le faire nous-mêmes, nous croyons devoir nous borner à renvoyer à leurs ouvrages.