Nous avons déjà rencontré la notion de champ gradient d’une fonction \(f\) de deux ou trois variables, à la Section 2.6. Il s’agit d’une nouvelle fonction \(\vnabla f\) qui à chaque point \(\vr\) associe un vecteur \(\vnabla f(\vr)\text{.}\) Au-delà des champs gradients et de leur utilité pour le calcul des dérivées directionnelles, la notion de champ de vecteurs est très utile pour modéliser plusieurs phénomènes. Pensons par exemple à la force \(\vF\) qui s’exerce sur une particule qui se déplace. La force \(\vF\) peut varier en fonction de la position de la particule (pensons à la proximité d’une charge électrique). De même, imaginons un fluide qui s’écoule dans une région du plan ou de l’espace. À chaque point de cette région, le fluide possède une vitesse spécifique; celle-ci variera en fonction d’éventuels obstacles rencontrés par le fluide et de la géométrie du milieu d’écoulement. Ceci motive la définition suivante.
Définition4.1.1.
Soit \(\cD \subseteq \R^n\) une région donnée. Un champ de vecteurs sur \(\cD\) est une fonction \(\vF : \cD \to \R^n\) qui à chaque point \(\vr\) de \(\cD\) associe un vecteur \(\vF(\vr)\text{.}\)
Les cas qui nous intéressent davantage sont :
Lorsque \(n=2\text{,}\)\(\cD\) est une région du plan. Dans ce cas, le champ \(\vF\) s’écrit au moyen de ses composantes \(\vF = P\vi + Q\vj\text{,}\) c’est-à-dire \(\vF(\vr) = P(\vr)\vi + Q(\vr)\vj\text{,}\) ou encore \(\vF(x,y) = P(x,y)\vi + Q(x,y)\vj\text{.}\)
Lorsque \(n=3\text{,}\)\(\cD\) est une région de l’espace. Dans ce cas, on écrit \(\vF = P\vi + Q\vj + R\vk \text{,}\) c’est-à-dire \(\vF(\vr) = P(\vr)\vi + Q(\vr)\vj + R(\vr)\vk\text{,}\) ou encore \(\vF(x,y,z) = P(x,y,z)\vi + Q(x,y,z)\vj + R(x,y,z)\vk\text{.}\)
On représente les champs de vecteurs au moyen de diagrammes. Cette représentation est surtout utile en deux dimensions, elle consiste à considérer un certain nombre de points sur la région \(\cD\) et d’en dessiner, au point \(\vr \in \cD\text{,}\) le vecteur \(\vF(\vr)\) ou un multiple scalaire de celui-ci (la constante sera la même pour tous les points).
Exemple4.1.2.
Soit \(\vF(x,y) = (1-y^2)\vi + 0 \vj\) défini pour \(-1 \leqslant y \leqslant 1\text{.}\)
Notons que la composante verticale du champ est nulle, de sorte qu’en tout point les flèches représentant ce champ seront horizontales. De plus, lorsque \(y = \pm 1\text{,}\) on aura \(\vF(x,y) = \vZero\text{.}\) Par ailleurs, les composantes du champ ne dépendent pas de la variable \(x\text{.}\) Ceci entraîne que le champ est invariant par translation horizontale. Notons aussi que la composante horizontale du champ, \(1-y^2\text{,}\) est toujours positive (les flèches pointent toujours vers la droite) et maximale lorsque \(y=0\text{.}\)
Ce champ peut être utilisé pour modéliser l’écoulement d’un fluide le long d’un canal se trouvant sur la bande horizontale \(-1 \leqslant y \leqslant 1\text{.}\) La vitesse du fluide est plus grande au milieu du canal, le frottement avec la paroi entraîne un ralentissement à mesure qu’on s’éloigne du milieu.
Exemple4.1.4.La source ponctuelle.
Imaginons que
tout l’espace est rempli avec un fluide incompressible, disons de l’eau;
d’une façon à déterminer, on trouve un moyen de produire plus d’eau à l’origine, disons \(4\pi m\) litres par seconde;
ceci force l’eau à couler vers l’extérieur, créant un “flux” symétrique en s’éloignant de l’origine.
Voyons quel est le champ de vecteurs résultant, \(\vv(x,y,z)\text{.}\) Comme le flux est symétrique, la vitesse au point \((x,y,z)\)
doit pointer en s’éloignant de l’origine, c’est-à-dire que le vecteur de direction du champ \(\vv(x,y,z)\) doit être le vecteur unitaire radial :
Nous devons simplement déterminer la fonction \(v(r)\text{.}\) Fixons \(r \geqslant0\) et concentrons-nous sur la sphère \(x^2+y^2+z^2=r^2\text{.}\) Elle est esquissée dans la figure ci-bas.
Pendant un très bref laps de temps de \(dt\) secondes, \(4\pi m\,dt\) litres d’eau sont créés à l’origine (le point rouge au centre). Puisque l’eau est incompressible, \(4\pi m\,dt\) litres d’eau doivent sortir à travers la sphère pendant cet intervalle de temps, afin de faire de la place pour la nouvelle eau.
En tout point de la sphère, l’eau s’écoule vers l’extérieur avec vitesse \(v(r)\text{.}\) Ainsi, pendant ce laps de temps, l’eau près de la surface de la sphère bouge d’une distance \(v(r)\,dt\text{.}\) En particulier, l’eau qui se trouvait dans la croûte mince \(r-v(r)\,dt \leqslant \sqrt{x^2+y^2+z^2} \leqslant r\ \ \) au début sort à travers la sphère \(\sqrt{x^2+y^2+z^2} = r\) pendant cet intervalle de temps. La croûte est esquissée en gris dans la figure. Le volume d’eau dans la partie grise est essentiellement l’aire de \(4\pi r^2\) multipliée par l’épaisseur, c’est-à-dire \(v(r)\,dt\text{.}\) Ainsi, l’égalité de volumes donne
La représentation de ce champ est obtenue comme suit. Pour chaque \((x,y)\text{,}\) le vecteur \(\frac{m}{r(x,y)}\,\hat\vr(x,y)\) est dessiné avec son origine en \((x,y)\text{.}\) Notons que \(\frac{m}{r(x,y)}\,\hat\vr(x,y)\)
pointe en direction opposée à l’origine,
a une norme \(\frac{m}{r(x,y)}\) qui
dépend seulement de \(r=|(x,y)|\text{,}\) et
est très longue lorsque \((x,y)\) est près de l’origine, et
a une norme décroissante \(\frac{1}{r}\) lorsque \(r=|(x,y)|\) grandit.
Voici une esquisse de quelques-uns de ces vecteurs.
Notons que, lorsque \(|(x,y)|\rightarrow 0\text{,}\) la norme de la vitesse \(|\vv(x,y)|\rightarrow\infty\text{.}\) Ceci est une conséquence de notre simplification qui consiste à penser que l’eau est produite en un seul point.
Exemple4.1.6.Le tourbillon.
Dans cet exemple, nous allons esquisser le champ de vecteurs
où \(\Om\) est une constante positive. Nous allons donner une procédure efficace pour produire une esquisse grossière, elle sera suffisante pour avoir une idée réaliste du contexte. Cette procédure se généralise à d’autres champs de vecteurs.
Pour commencer, concentrons-nous sur la composante horizontale \(\vi\cdot\vv(x,y)\) de notre champ. Voyons dans quelle région du plan \(Oxy\) celle-ci est positive, négative ou nulle.
Ceci partitionne naturellement le plan \(Oxy\) en 9 morceaux, selon que chacune de ces composantes est positive, négative ou nulle :
\(\vi\cdot\vv \geqslant0\) et \(\vj\cdot\vv \geqslant0\) sur \(\big\{\ (x,y)\in\R^2\ \big|\ y \lt 0,\ x \geqslant0\ \big\}\)
\(\vi\cdot\vv \geqslant0\) et \(\vj\cdot\vv=0\) sur \(\big\{\ (x,y)\in\R^2\ \big|\ y \lt 0,\ x=0\ \big\}\)
\(\vi\cdot\vv \geqslant0\) et \(\vj\cdot\vv \lt 0\) sur \(\big\{\ (x,y)\in\R^2\ \big|\ y \lt 0,\ x \lt 0\ \big\}\)
\(\vi\cdot\vv=0\) et \(\vj\cdot\vv \geqslant0\) sur \(\big\{\ (x,y)\in\R^2\ \big|\ y=0,\ x \geqslant0\ \big\}\)
et ainsi de suite.
Pensons maintenant à \(\vv(x,y)\) comme la vitesse au point \((x,y)\) d’un fluide.
Regardons le premier point ci-haut. Il nous dit que, sur la première des 9 régions, c’est-à-dire, sur \(\big\{\ (x,y)\in\R^2\ \big|\ y \lt 0,\ x \geqslant0\ \big\}\text{,}\) qui est le quatrième quadrant, la vitesse horizontale \(\vi\cdot\vv \geqslant0\text{.}\) Ceci veut dire que le fluide bouge vers la droite. Nous signifions ceci dans notre esquisse par une flèche qui va horizontalement vers la droite (il s’agit de la flèche bleue dans la figure ci-bas). La composante verticale \(\vj\cdot\vv \geqslant0\text{,}\) ce qui signifie que le fluide bouge aussi vers le haut. On signifie ceci par la flèche verticale vers le haut (c’est la flèche rouge).
Par la suite, concentrons-nous sur le deuxième point dans la liste ci-haut. Il nous dit que, sur la deuxième région, c’est-à-dire sur \(\big\{\ (x,y)\in\R^2\ \big|\ y \lt 0,\ x=0\ \big\}\text{,}\) qui est la moitié inférieure de l’axe \(Oy\text{,}\) la composante horizontale \(\vi\cdot\vv \geqslant0\text{,}\) ce qui veut dire que le fluide bouge vers la droite. On indique ceci par une flèche horizontale vers la droite dont la source est un point sur la moitié inférieure de l’axe \(Oy\) (c’est la deuxième flèche bleue dans la figure ci-après). La composante verticale est \(\vj\cdot\vv=0\text{,}\) ce qui veut dire que le fluide n’a pas de vitesse verticale. On indique ceci par l’absence de flèche verticale ayant pour source un point sur la partie inférieure de l’axe des \(y\text{.}\)
On poursuit de la sorte.
Lorsque nous aurons considéré les 9 régions, nous aurons l’esquisse suivante.
De ceci, on voit que, par exemple, dans le premier quadrant,
le fluide bouge vers la gauche et vers le haut,
le fluide croise l’axe \(x\) verticalement (donc, près de l’axe \(x\text{,}\) les flèches seront presque verticales), et
le fluide croise l’axe des \(y\) horizontalement (donc, près de l’axe \(y\text{,}\) les flèches seront presque horizontales), et
il y a un point, à savoir \((0,0)\text{,}\) où le champ de vitesses est exactement \(\vZero\text{.}\) De plus, \(\vv(x,y)=\Omega(-y\vi+x\vj)\) est plus petit lorsque \((x,y)\) est plus près de \((0,0)\text{,}\) et \(\vv(x,y)\) a une grande norme lorsque \((x,y)\) est loin de \((0,0)\text{.}\)
Avec toute cette information, nous pouvons esquisser notre champ de vecteurs.
Ceci montre le champ qui tourbillonne autour de l’origine, d’où son nom.
Si \(f\) est une fonction de deux ou trois variables, alors son gradient \(\vnabla f\) définit un champ de vecteurs (là où il existe !). Nous avons déjà considéré le gradient à la Section 2.6, mais d’un point de vue un peu différent. Supposons que nous ayons une fonction de deux variables, \(f(x,y)\text{.}\) Son champ gradient est le champ de vecteurs \(\vnabla f = f_x \vi + f_y \vj\) ou, plus explicitement, \(\vnabla f(x,y) = f_x(x,y)\vi + f_y(x,y)\vj\text{.}\) De la même façon, si l’on a une fonction de trois variables \(f\text{,}\) son gradient est \(\vnabla f = f_x\vi + f_y\vj + f_z \vk\text{,}\) ou encore \(\vnabla f(x,y,z) = f_x(x,y,z)\vi + f_y(x,y,z)\vj + f_z(x,y,z)\vk\text{.}\)
Exemple4.1.9.
Soit \(f(x,y) = x\sin y + y\text{.}\) Son champ gradient est \(\vnabla f(x,y) = \sin y\, \vi + (x\cos y + 1)\vj\text{.}\) La figure ci-après, à gauche, montre le champ gradient. Dans la figure de droite, on trouve le champ gradient \(\vnabla f\) ainsi que les courbes de niveau de \(f\text{.}\)
Le champ gradient de \(f(x,y) = x\sin y + y\)
Le champ gradient de \(f(x,y) = x\sin y + y\) ainsi que ses courbes de niveau
On remarque que la norme du gradient est plus grande là où les courbes de niveau de \(f\) sont le plus près les unes des autres. Ceci résulte du fait que la norme du gradient est la dérivée directionnelle de la fonction dans la direction du gradient. Plus cette dérivée est grande, plus la valeur de \(f\) change vite.
Dans les sections qui suivent, nous devront faire face au problème de décider si un champ \(\vF\) est est un champ gradient ou pas. Nous reviendrons sur ce point. Contentons-nous de donner un exemple.
Exemple4.1.10.
Montrer qu’il n’existe pas de fonction \(f\) continue telle que \(\vnabla f(x,y) = xy\, \vi + (y^2+1)\, \vj = P(x,y)\,\vi + Q(x,y)\,\vj\text{.}\)
Solution.
Si \(f\) était une telle fonction, ses dérivées partielles seraient \(P\) et \(Q\text{,}\) qui sont de classe \(C^1\text{.}\) On devrait donc avoir \(f\) de classe \(C^2\text{.}\) Ainsi, en vertu du théorème de Clairaut,
Or, \(\pdiff{P}{y} = 2y\) et \(\pdiff{Q}{x} = y\text{.}\) Ainsi, une telle fonction \(f\) ne peut pas exister.
En trois dimensions, le concept est le même. La représentation graphique est toujours de mise, même si elle pourrait être un peu plus difficile à interpréter. Dans ce cas, il est commun d’utiliser aussi des couleurs pour les flèches.
Exemple4.1.11.
On considère le champ de vecteurs \(\vF = y\vi + \vj +(z^2 -1) \vk\text{.}\)
La composante \(y\) du champ est constante, les flèches pointent toujours dans la direction des \(y\) positifs. Les deux autres composantes dépendent de la position.
Sous-sectionLes lignes de courant
Étant donné un champ de vecteurs \(\vF\text{,}\) une ligne de courant de \(\vF\) est une courbe paramétrée \(\cC\) telle que, en tout point de \(\cC\text{,}\) le vecteur vitesse est donné par \(\vF\) évalué en ce point. Si le champ est le champ de vitesses d’un fluide et que l’on dépose un objet qui flotte en un point donné, la ligne de courant passant par ce point représente la trajectoire que cet objet suivrait, mené par le courant.
Plus précisément, si \(\vF = p\vi + Q \vj\) est donné et \(\cC\) est une ligne de courant paramétrée par \(\vr(t) = x(t) \vi + y(t) \vj\text{,}\) alors le vecteur vitesse est \(\vv(t) = x'(t) \vi + y'(t) \vj\) et il doit vérifier
Remarquons que différentes valeurs de \(C_1\) et \(C_2\) donnent différentes lignes de courant. L’interprétation est claire : non seulement il faut connaître la vitesse \(\vF\text{,}\) mais il faut aussi connaître un point sur la trajectoire, ce qui permet de déterminer \(C_1\) et \(C_2\text{.}\)
Dans cet exemple, nous avons \(x(t) \cdot y(t) = C_1 C_2 \text{,}\) de sorte que les lignes de courant se trouvent sur des hyperboles d’équations \(xy = k\text{.}\) Notons que l’équation \(xy = 0\) a pour courbe la réunion des deux axes. Chaque demi-axe est à lui tout seul une ligne de courant.
Voyons maintenant comment trouver les lignes de courant d’un champ donné.
Exemple4.1.16.
Trouver les lignes de courant du champ de vecteurs \(\vF(x,y) = -2y\,\vi + x\, \vj\text{.}\)
Solution.
Supposons que nous avons une ligne de courant \(\vr(t) = x(t) \,\vi + y(t) \, \vj\text{.}\) Puisqu’il s’agit d’une ligne de courant, nous devons avoir
On remarque que ce vecteur n’est pas égal à \(-2y(u)\vi + x(u) \vj\text{,}\) il lui est seulement proportionnel. Ce paramétrage ne vérifie donc pas la définition des lignes de courant. Cependant, posons plutôt \(\vr(t) = \sqrt{2C}\cos (t\sqrt{2})\, \vi + \sqrt{C}\sin (t\sqrt{2})\, \vj\text{.}\) On vérifie aisément qu’il s’agit d’un paramétrage de l’ellipse et, par ailleurs, la condition de la définition des lignes de courant est satisfaite.
En général il est difficile de déterminer les lignes de courant d’un champ de vecteurs. Il s’agit de quelque chose qui est plutôt abordé dans les cours d’équations différentielles ou de mathématiques avancées.
ExercicesExercices
Pour se pratiquer.
1.
Ci-dessous se trouve une esquisse du champ de vecteurs \(\vv(x,y)\text{.}\)
Trouver les régions où les coordonnées \(x\) et \(y\) sont positives, négatives et nulles :
\begin{equation*}
\vv(x,y)\cdot \vi \begin{cases}
\gt 0 & \mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
=0 &\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
\lt 0&\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $}
\end{cases}
\end{equation*}
\begin{equation*}
\vv(x,y) \cdot \vj
\begin{cases}
\gt 0 & \text{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
=0 & \text{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
\lt 0&\text{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $}
\end{cases}
\end{equation*}
Supposer que \(\vv(x,y)\) se comporte comme prévu aux points qui ne sont pas visibles. Autrement dit, les échantillons sont représentatifs d’une fonction vectorielle continue et lisse. Supposer également que les repères sur les axes correspondent à des distances unitaires.
Réponse.
\begin{equation*}
\vv(x,y)\cdot \vi \left.\begin{cases}
\gt 0 & \text{ lorsque } \fbox{$x \gt 0$}\\
=0 &\text{ lorsque } \fbox{$x=0$}\\
\lt 0&\text{ lorsque } \fbox{$x \lt 0$}
\end{cases}\right\}
\end{equation*}
et
\begin{equation*}
\vv(x,y)\cdot \vj \left.\begin{cases}
\gt 0 & \text{ lorsque } \fbox{$-2 \lt x \lt 2$}\\
=0 &\mbox{ lorsque } \fbox{$x \in \{-2,2\}$}\\
\lt 0&\mbox{ lorsque } \fbox{$x \lt -2$ ou $x \gt 2$}
\end{cases}\right\}
\end{equation*}
au moins pour \((x,y)\) montré dans l’esquisse.
2.
Ci-dessous se trouve une esquisse du champ de vecteurs \(\vv(x,y)\text{.}\)
Trouver les régions où les coordonnées \(x\) et \(y\) sont positives, négatives et nulles :
\begin{equation*}
\vv(x,y)\cdot \vi \begin{cases}
\gt 0 & \mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $}\\
=0 &\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
\lt 0&\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $}
\end{cases}
\end{equation*}
\begin{equation*}
\vv(x,y)\cdot \vj \begin{cases}
\gt 0 & \mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
=0 &\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $} \\
\lt 0&\mbox{ lorsque } \fbox{$\vphantom{L}\qquad\qquad $}
\end{cases}
\end{equation*}
Supposer que \(\vv(x,y)\) se comporte comme prévu aux points qui ne sont pas visibles. Supposer également que les repères sur les axes correspondent à des distances unitaires.
Réponse.
\begin{equation*}
\vv(x,y)\cdot \vi \left.\begin{cases}
\gt 0 & \mbox{ lorsque } \fbox{$y \gt -x$}\\
=0 &\mbox{ lorsque } \fbox{$y=-x$}\\
\lt 0&\mbox{ lorsque } \fbox{$y \lt -x$}
\end{cases}\right\}
\end{equation*}
et
\begin{equation*}
\vv(x,y)\cdot \vj \left.\begin{cases}
\gt 0 & \mbox{ lorsque } \fbox{$y \lt x$}\\
=0 &\mbox{ lorsque } \fbox{$y=x$}\\
\lt 0&\mbox{ lorsque } \fbox{$y \gt x$}
\end{cases}\right\}
\end{equation*}
au moins pour \((x,y)\) montré dans l’esquisse.
3.
Une plateforme comportant de nombreux petits tapis roulants est alignée sur un plan de coordonnées. Chaque bande transporteuse déplace un objet situé au-dessus d’elle dans la direction de l’origine, et une bande transporteuse en position \((x,y)\) provoque le déplacement d’un objet situé au-dessus d’elle à une vitesse \(y\text{.}\) Supposer que les objets n’interfèrent pas les uns avec les autres.
Donner une formule vectorielle pour la vitesse d’un objet en position \((x,y)\text{.}\)
Réponse.
\(\vv(x,y)=\frac{-y}{\sqrt{x^2+y^2}}(x,y)\)
4.
Trouver les lignes de courant du champ de vecteurs \(\vF(x,y) = x^2\, \vi + y^2\, \vj\text{.}\)
Pour mieux comprendre.
5.
Soit \(\vF = P\,\vi + Q\,\vj\) le champ de vecteurs à deux dimensions esquissé ci-dessous.
Déterminer les signes de \(P\text{,}\)\(Q\text{,}\)\(\pdiff{Q}{x}\) et \(\pdiff{Q}{y}\) au point \(A\text{.}\)
Réponse.
\(P \gt 0\)
\(Q \gt 0\)
\(\pdiff{Q}{x} \lt 0\)
\(\pdiff{Q}{y} \gt 0\)
6.
Imaginer que le champ de vecteurs \(\vv(x,y) = x\,\vi+y\,\vj\) est le champ de vitesses d’un fluide en mouvement.
Au temps \(0\text{,}\) vous laissez tomber une brindille dans le fluide au point \((1,1)\text{.}\) Quelle est la position approximative de la brindille au moment \(t=0.01\text{?}\)
Au temps \(0\text{,}\) vous laissez tomber une brindille dans le fluide au point \((0,0)\text{.}\) Quelle est la position de la brindille au moment \(t=0.01\text{?}\)
Au temps \(0\text{,}\) vous laissez tomber une brindille dans le fluide au point \((0,0)\text{.}\) Quelle est la position de la brindille au moment \(t=10\text{?}\)
Réponse.
(a) \((1.01,\,1.01)\)
(b) \((0,\,0)\)
(c) \((0,\,0)\)
7.
Imaginer que le champ de vecteurs \(\vv(x,y) = 2x\,\vi -\vj\) est le champ de vitesses d’un fluide en mouvement. Au temps \(0\text{,}\) vous laissez tomber une brindille dans le fluide au point \((0,0)\text{.}\) Quelle est la position de la brindille au temps \(t=10\) ?
Réponse.
\((0\,,\,-10)\)
8.
Des abeilles amicales volent vers votre visage depuis toutes les directions. La vitesse de chaque abeille est inversement proportionnelle à sa distance par rapport à votre visage. Trouver un champ de vecteurs pour la vitesse de l’essaim.
Réponse.
Si votre visage est à l’origine, alors \(\vv(x,y,z)=-\frac{\alpha}{\sqrt{x^2+y^2+z^2}}(x,y,z)\) pour une certaine constante positive \(\alpha\text{.}\)
9.
Esquisser le champ de vecteurs \(\vv(x,y)=(x^2,y)\text{.}\)