Analyse vectorielle

De testwiki
Aller à la navigation Aller à la recherche

L'analyse vectorielle est une branche des mathématiques qui étudie les champs de scalaires et de vecteurs suffisamment réguliers des espaces euclidiens, c'est-à-dire les applications différentiables d'un ouvert d'un espace euclidien Modèle:Formule à valeurs respectivement dans et dans Modèle:Formule. Du point de vue du mathématicien, l'analyse vectorielle est donc une branche de la géométrie différentielle. Cette dernière inclut l'analyse tensorielle qui apporte des outils plus puissants et une analyse plus concise entre autres des champs de vecteurs.

Mais l'importance de l'analyse vectorielle provient de son utilisation intensive en physique et dans les sciences de l'ingénieur. C'est de ce point de vue qu'elle sera présentée dans l'article, et c'est pourquoi elle sera le plus souvent limitée au cas où E=3 est l'espace usuel à trois dimensions. Dans ce cadre, un champ de vecteurs associe à chaque point de l'espace un vecteur (à trois composantes réelles), tandis qu'un champ de scalaires y associe un réel. Par exemple, dans le cas de l'eau d'un lac, la donnée de sa température en chaque point forme un champ de scalaires, et celle de sa vitesse en chaque point, un champ de vecteurs (pour une approche plus théorique, voir Géométrie différentielle).

Le calcul vectoriel et l'analyse vectorielle furent développés à la fin du Modèle:S- par J. Willard Gibbs et Oliver Heaviside à partir de la théorie des quaternions (due à Hamilton) ; la plupart des notations et de la terminologie furent établies par Gibbs et Edwin Bidwell Wilson dans leur livre de 1901, Modèle:Langue (Analyse vectorielle).

Principaux opérateurs différentiels linéaires

Le gradient, la divergence et le rotationnel sont les trois principaux opérateurs différentiels linéaires du premier ordre. Cela signifie qu'ils ne font intervenir que des dérivées partielles (ou différentielles) premières des champs, à la différence, par exemple, du laplacien qui fait intervenir des dérivées partielles du second ordre.

On les rencontre en particulier :

Opérateur formel nabla

Modèle:Article détaillé L'opérateur nabla tire son nom d'une lyre antique qui avait la même forme de triangle pointant vers le bas. Il s'agit d'un opérateur formel de 3 défini en coordonnées cartésiennes par :

=(xyz).

On écrit aussi pour souligner que formellement, l'opérateur nabla a les caractéristiques d'un vecteur. Il ne contient certes pas de valeurs scalaires, mais on utilise ses éléments constitutifs (que l'on peut voir comme des opérations en attente d'argument — des opérateurs différentiels) très exactement comme on aurait utilisé les valeurs scalaires composant un vecteur.

La notation nabla fournit un moyen commode pour exprimer les opérateurs vectoriels en coordonnées cartésiennes ; dans d'autres systèmes de coordonnées, elle est encore utilisable au prix de précautions supplémentaires ; pour plus de précisions, et des interprétations plus théoriques (en particulier la relation avec la dérivée covariante), voir les articles détaillés Nabla et Connexion de Koszul.

Opérateur différentiel gradient

Modèle:Article détaillé Le gradient est un opérateur linéaire qui s'applique à un champ de scalaires et décrit un champ de vecteurs qui représente la variation de la valeur du champ scalaire dans l'espace Modèle:Mvar considéré. Concrètement, le gradient indique la direction de la plus grande variation du champ scalaire et l'intensité de cette variation. Par exemple, sur une surface définie par l'équation réduite Modèle:Math, où f:E=2 est une application différentiable, le gradient de l'altitude est dirigé selon la composante horizontale de la ligne de plus grande pente (montante), et sa norme augmente avec la pente.

Pour tout vecteur Modèle:Mvar,
df(a)(h)=(gradaf)h,
Modèle:Math désigne la valeur sur le vecteur Modèle:Mvar de la différentielle de la fonction Modèle:Mvar au point Modèle:Mvar.
C'est donc la définition de l'application linéaire tangente du champ scalaire Modèle:Math en Modèle:Formule.
  • On en déduit facilement[1] que :
En dimension 3, le vecteur normal en a=(xa,ya,za) à la surface d'équation Modèle:Formule, où Modèle:Mvar est une constante, est donné par :
gradaf.
  • Il en résulte immédiatement[2] que :
Pour tout vecteur Modèle:Mvar, la dérivée directionnelle fv(a) de la fonction Modèle:Mvar suivant le vecteur Modèle:Mvar au point Modèle:Mvar est donnée par :
(gradaf)v=(f(a))v=df(a)(v).
  • En dimension 3 et coordonnées cartésiennes (en base orthonormée), le champ de gradients vérifie :
grad f=f=(fxfyfz).
Cette relation peut servir, dans le cas particulier où elle s'applique, de définition du gradient. Elle se généralise naturellement en dimension quelconque, en ajoutant des composantes au nabla.

Application linéaire tangente d'un champ de vecteurs

Soit Modèle:Formule le point translaté de Modèle:Formule par la translation de vecteur h ; alors :

F(M)F(M)=(F^)Mh+o(h)

définit l'opérateur linéaire noté par un chapeau pour signifier que sa représentation dans une base est une matrice carrée [3-3], application linéaire tangente du champ de vecteurs Modèle:Formule.

Le déterminant de cet opérateur est le Jacobien de la transformation qui à Modèle:Formule associe Modèle:Formule.

Sa trace définira la divergence du champ de vecteurs Modèle:Formule.

Cela permettra de donner du rotationnel du champ de vecteurs Modèle:Formule une définition intrinsèque.

On pourra vérifier que symboliquement :

(F^)Mh=(h)F

Opérateur divergence

Modèle:Article détaillé La divergence s'applique à un champ de tenseurs d'ordre Modèle:Math et le transforme en un champ de tenseurs d'ordre Modèle:Math. Pratiquement, la divergence d'un champ de vecteurs exprime sa tendance à fluer localement hors d'un petit volume entourant le point Modèle:Formule où est calculée la divergence.

En dimension 3 et en coordonnées cartésiennes, si F est un tenseur d'ordre 1, alors c'est un vecteur et on peut définir la divergence par la relation

divF=F=Fxx+Fyy+Fzz

F=(Fx,Fy,Fz) désigne le champ de vecteurs auquel est appliqué l'opérateur divergence. La divergence peut être vue, formellement, comme le produit scalaire de l'opérateur nabla par le vecteur « générique » du champ auquel elle est appliquée, ce qui justifie la notation . Bien entendu, cette définition se généralise naturellement en dimension quelconque.

La définition indépendante du choix de la base est :

divF=Tr(F^).

Une autre définition possible, plus générale mais plus difficile à formaliser, consiste à définir la divergence d'un champ de vecteurs en un point comme le flux local du champ autour de ce point.

Opérateur rotationnel

Modèle:Article détaillé

Le rotationnel transforme un champ de vecteurs en un autre champ de vecteurs. Plus difficile à se représenter aussi précisément que le gradient et la divergence, il exprime la tendance qu'a un champ à tourner autour d'un point : sa circulation locale sur un petit lacet entourant le point M est non nulle. Par exemple :

  • dans une tornade, le vent tourne autour de l'œil du cyclone et le champ de vecteurs vitesse du vent a un rotationnel non nul autour de l'œil. Le rotationnel de ce champ de vitesse (autrement dit le champ de vorticité ou encore champ tourbillon) est d'autant plus intense que l'on est proche de l'œil ;
  • le rotationnel du champ des vitesses V(M)=Ω0OM d'un solide qui tourne à vitesse constante Ω0 est constant, dirigé selon l'axe de rotation et orienté de telle sorte que la rotation ait lieu, par rapport à lui, dans le sens direct, et vaut simplement 2Ω0.

Dans un espace à 3 dimensions et en coordonnées cartésiennes, on peut définir le rotationnel par la relation

rot F=F=(Fz/yFy/zFx/zFz/xFy/xFx/y)

F=(Fx,Fy,Fz) désigne le champ de vecteurs auquel est appliqué l'opérateur rotationnel. L'analogie formelle avec un produit vectoriel justifie la notation .

Cela peut aussi s'écrire, par abus de notation (c'est aussi une astuce mnémotechnique), à l'aide d'un déterminant :

rot F=|ijkxyzFxFyFz|

(i,j,k) désigne la base canonique. Cette dernière expression est un peu plus compliquée que la précédente, mais elle se généralise facilement à d'autres systèmes de coordonnées.

  • Une définition intrinsèque (parmi d'autres) du rotationnel est la suivante :

À partir du champ F, on peut construire le champ X0F (où X0 est un vecteur uniforme) dont la divergence est une forme linéaire de X0 et donc exprimable par un produit scalaire KX0, où K est l'opposé du rotationnel de F :

div(X0F)=rotFX0.

Une autre définition possible, plus générale mais plus difficile à formaliser, consiste à définir le rotationnel d'un champ de vecteurs en un point comme la circulation locale du champ autour de ce point.

Opérateurs d'ordre supérieur

Opérateur laplacien

Le plus utilisé des opérateurs d'ordre 2 est le laplacien, du nom du mathématicien Pierre-Simon de Laplace. Le laplacien d'un champ est égal à la somme des dérivées secondes de ce champ par rapport à chacune des variables.

En dimension 3 et en coordonnées cartésiennes, il s'écrit :

Δ=2=2x2+2y2+2z2.

Cette définition a un sens aussi bien pour un champ de scalaires que pour un champ de vecteurs. On parle respectivement de laplacien scalaire et de laplacien vectoriel. Le laplacien scalaire d'un champ de scalaires est un champ de scalaires alors que le laplacien vectoriel d'un champ de vecteurs est un champ de vecteurs. Pour distinguer ce dernier, on le note parfois Δ (afin que les novices n'oublient pas qu'il s'agit de l'opérateur grad divrot rot) ; la notation Δ est plutôt à déconseiller.

L'autre notation du laplacien qui apparaît ci-dessus, 2, invite à le considérer, formellement, comme le carré scalaire de l'opérateur nabla «  ».

Le laplacien apparaît dans l'écriture de plusieurs équations aux dérivées partielles qui jouent un rôle essentiel en physique.

Opérateur laplacien vectoriel

Le laplacien d'un champ de vecteurs A est un vecteur défini par le laplacien scalaire de chacune des composantes du champ vectoriel, ainsi en coordonnées cartésiennes, il est défini par :

ΔA=2A=(.)A=(2x2+2y2+2z2).[AxAyAz]=[2Axx2+2Axy2+2Axz22Ayx2+2Ayy2+2Ayz22Azx2+2Azy2+2Azz2]=[ΔAxΔAyΔAz]

Le laplacien vectoriel est présent :

Quelques formules différentielles

Attention : les formules suivantes sont valables à condition que certaines hypothèses soient vérifiées (la fonction scalaire dans la première formule doit être 𝒞2(Ω), où Ω, par exemple. De même, si f désigne la fonction vectorielle concernée dans la seconde formule, il faut vérifier f𝒞2(Ω), Ωn.)

  • rot(grad)=0
  • div(rot)=0
  • Δ=div(grad) (appliqué à un scalaire)

Formules dites de Leibniz pour les produits

  • grad(X0B)=(X0grad)B+X0rot(B) (où X0 est un vecteur uniforme) et évidemment :
  • grad(AB)=(Agrad)B+ArotB+(Bgrad)A+BrotA
  • grad(FF)=2(Fgrad)F+2F(rot(F)) (dite de Bernoulli, en mécanique des fluides)
  • div(X0B)=X0rot(B) (où X0 est un vecteur uniforme, définition intrinsèque du rotationnel)
  • div(AB)=ArotB+BrotA
  • rot(X0B)=X0divB(X0grad)B (où X0 est un vecteur uniforme, par définition de l'application linéaire tangente)
  • rot(AB)=A divB(Agrad)BB divA+(Bgrad)A


  • grad(fg)=fgrad(g)+ggrad(f) (symétrique en f et g)
  • div(ρV)=ρdivV+grad(ρ)V
  • rot(ρV)=ρrotV+grad(ρ)V
  • Δ(fg)=fΔg+2grad(f)grad(g)+gΔf
  • div(fgrad(g)ggrad(f))=fΔggΔf

Quelques formules utiles

  • Soient Modèle:Formule et Modèle:Formule deux champs scalaires, il existe un champ de vecteurs A(M) tel que :
    rotA=gradfgradg
  • Le champ central OM=r joue un rôle très important en physique. Aussi convient-il de mémoriser ces quelques évidences :
    • son application linéaire tangente est la matrice identité (cf. la définition !),
    • donc divr=3 et rot(X0r)=2X0 (où X0 est un vecteur uniforme) et rot(r)=0
  • D'autre part X0=grad(X0r) (où X0 est un vecteur uniforme).
  • Et aussi : gradf(r)=f(r)u avec u=rren particulier grad(r2)=2r (évident car d(rr)=d(r2))
  • Δf(r)=f(r)+2rf(r) , sauf en r=0
  • Le champ newtonien, soit rr3, est très souvent étudié car c'est le seul champ central à divergence nulle (évident si l'on pense en termes de flux) sauf pour Modèle:Formule, où elle vaut 4πδ(r) ; ce résultat est le théorème de Gauss pour l'angle solide). Il en résulte que Δ(1/r)=4πδ(r). Donc
    Δ(X0/r)=4πX0δ(r)
    (où X0 est un vecteur uniforme) qui se décompose en :
    grad(div)(X0/r)=4πX0δ(r)(1/3)

(où X0est un vecteur uniforme), et

rot(rot)(X0/r)=+4πX0δ(r)(2/3)

(où

X0

est un vecteur uniforme) ce qui est moins évident (cf. moment magnétique).

Expressions des opérateurs en différentes coordonnées

Coordonnées cylindriques

gradf=frur+1rfθuθ+fzuz
divA=1rr(rAr)+1rAθθ+Azz
rot(A)=(1rAzθAθz)ur+(ArzAzr)uθ+1r(r(rAθ)Arθ)uz
Δf=1rr(rfr)+1r22fθ2+2fz2

Coordonnées sphériques

gradf=frur+1rfθuθ+1rsinθfφuφ
divA=1r2r(r2Ar)+1rsinθθ(sinθAθ)+1rsinθAφφ
rot(A)=1rsinθ(θ(sinθAφ)Aθφ)ur+(1rsinθArφ1rr(rAφ))uθ+1r(r(rAθ)Arθ)uφ
Δf=1r2r(r2fr)+1r2sinθθ(sinθfθ)+1r2sin2θ2fφ2

Annexes

Modèle:Autres projets

Bibliographie

Références

Modèle:Références

Articles connexes

Modèle:Palette

Modèle:Portail

  1. Pour tout vecteur Modèle:Mvar, f(a+h)f(a)=(gradaf)h+o(h).
    Or, pour tout vecteur Modèle:Mvar tel que Modèle:Math ∈ surface, Modèle:Math ;
    donc Modèle:Math ; donc (gradaf)h=o(h).
  2. Il y a plusieurs définitions non équivalentes entre elles d'une dérivée directionnelle.
    Explication seulement pour tout vecteur unitaire Modèle:Mvar (fixé) :
    Pour tout vecteur Modèle:Mvar, f(a+h)f(a)=(gradaf)h+o(h).
    Or, pour tout réel Modèle:Mvar, (gradaf)(tu)=((gradaf)u)t ;
    donc f(a+tu)f(a)=((gradaf)u)t+o(|t|).