Maths TD7

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 5

Master Dynamique terrestre et risques naturels Math ematiques pour g eologues

Op erateurs di erentiels
On etudie en g eosciences des fonctions scalaires des coordonn ees despace, comme la temp erature, ou bien des vecteurs dont les trois composantes sont des fonctions des coordonn ees, comme la pesanteur ou le champ magn etique. Lorsque ces fonctions ont des d eriv ees partielles, on peut d enir dautres scalaires ou vecteurs qui restent les m emes pour tout r ef erentiel, ce quon appelle des invariants de la fonction ou du champ de vecteurs. Ils en fournissent des propri et es intrins` eques et locales. Pour une fonction, les invariants qui nous seront utiles sont le gradient (un vecteur) et le laplacien (un scalaire). Pour un champ de vecteurs ce sont le rotationnel (un vecteur), la divergence (un scalaire) et le laplacien vectoriel (un vecteur).

Produit scalaire et vectoriel

Soit deux vecteurs a et b ayant pour composantes dans un r ef erentiel cart esien ax , ay , az et bx , by , bz respectivement. On d enit : le produit scalaire : a.b = ax b x + ay by + az bz ay b z az b y le produit vectoriel : a b = az bx ax bz ax b y ay b x

Notion de circulation dun champ de vecteurs

On appelle travail de A ` a B du vecteur a le long dune courbe (C ), dont un segment innit esimal est le vecteur dl, la somme des produits scalaires innit esimaux
B

a.dl.
A

Lorsque (C ) est une courbe ferm ee (A = B ), la position de A sur (C ) na pas dimportance, seulement le sens de parcours. Le travail est alors appel e la circulation de a sur (C ), a.dl.

Notion de ux dun champ de vecteurs

Consid erons une surface continue (S ). Soit n sa normale unitaire, toujours issue du m eme c ot e. On appelle ux du vecteur a ` a travers (S ) le scalaire : = a.ndS

Cette notion sera utilis ee pour introduire la divergence.

Le gradient

La forme di erentielle totale dune fonction f (x, y, z ), o` u x, y et z sont les trois variables de lespace, est f f f df = dx + dy + dz, x y z

qui peut s ecrire sous la forme dun produit scalaire de deux vecteurs u et dl avec f dx x f u = y dl = dy f dz
z

Le champ vectoriel u sexprime par un op erateur nomm e gradient que lon note : grad(f ) = f =
f x f y f z

Ce vecteur nest autre quune extension de la classique d eriv ee dune fonction ` a un espace de dimension sup erieure. Il indique donc la pente locale de la fonction, le vecteur obtenu etant dirig e le long de la plus grande pente au champ f . f est orthogonal aux isosurfaces de f .

Expression en coordonn ees cyclindriques


Sur la base (er , e , ez ) on a f =
f r 1 f r f z

Expression en coordonn ees sph eriques


Sur la base (er , e , e ) on a f =
f r 1 f r f 1 r. sin

La divergence
ux uy uz + + . x y z

La divergence dun champ vectoriel u est un scalaire d eni par : div (u) = .u =

An de d enir le sens physique de la divergence consid erons un volume rectangulaire de c ot es dx, dy et dz .


z

Ux dz dy x

dx y

Le ux de u sortant de la face de droite dans la direction x est ux (x + dx, y, z )dydz . De m eme le ux de u entrant par la face de gauche dans la direction x est ux (x, y, z )dydz . Le bilan de ux entre ces deux surfaces est donc ux ux dxdydz = dV x x Le m eme raisonnement peut etre fait dans la direction y et dans la direction z . Le bilan de ux au travers des faces du volume peut donc s ecrire [ux (x + dx, y, z ) ux (x, y, z )]dydz = d = uz ux uy + + x y z dV = (.u)dV.

L equation ci-dessus permet dexprimer ce quest la divergence : la divergence est donc le bilan de ux dun champ de vecteurs par unit e de volume. Sous forme int egrale on obtient le th eor` eme de la divergence (Green-Ostrogradsky) : u.dS =
int

div (u)dV

Le ux ` a travers une surface ferm ee S dun champ u est egale ` a lint egrale sur tout le volume d elimit e par S de la divergence de u. Si u est un champ de vitesse, alors la divergende de u mesure laccroissement total de volume par unit e de temps et par unit e de volume. Si en un point A la divergence est positive (n egative) alors A est un point dexpension (de compression). Si la divergence de u est nulle en tout point dune r egion D alors le corps ayant le champ de vitesse u est incompressible dans cette r egion.

Expression en coordonn ees cyclindriques


.u = uz 1 (r ur ) 1 u + + . r r r z

Expression en coordonn ees sph eriques


.u = 1 (r2 ur ) 1 (sin u ) 1 u + + r2 r r sin r sin

Le rotationnel

Le rotationnel dun champ vectoriel u est un vecteur d eni par uy uz y z uz x . rotu = u = u z x uy ux x y On interpr` ete le rotationnel u comme laxe dun tourbillon en m ecanique des uides, le champ u repr esentant la vitesse du liquide. Le rotationnel dun gradient est nul : rot(grad(f ) = f = 0 Soit u un champ vectoriel et (C ) un contour ferm e. On peut montrer que le ux du rotationnel de u` a travers une surface sappuyant sur (C ) est egale ` a la circulation le long de (C ) du champ u. Ce th eor` eme est appel e th eor` eme du rotationnel (Stokes) u.dl = u.dS

Expression en coordonn ees cyclindriques

1 r u 1 uz r z uz ur z r (r u ) ur ) (r r

u =

Expression en coordonn ees sph eriques

1 r sin 1 r sin 1 r (sin u ) u (r sin u ur r (r u ) ur r

u =

Le laplacien

Le dernier op erateur que nous utiliserons est le laplacien. Le laplacien est d eni comme la divergence du gradient. On distingue le laplacien scalaire laplacien(f ) = div (grad(f )) = f = 2 f = De m eme on d enit le laplacien vectoriel comme
2 ux x2 2 uy x2 2 uz x2

2f 2f 2f + + . x2 y 2 z 2

+ + +

u = 2 u =

2 ux y 2 2 uy y 2 2 uz y 2

+ + +

2 ux z 2 2 uy z 2 2 uz z 2

Le laplacien dune fonction mesure la di erence entre la valeur de la fonction en un point et sa moyenne autour de ce point. Ainsi le laplacien est nul ou tr` es petit lorsque la fonction varie sans a coups. `

Expression en coordonn ees cyclindriques


2 f = 1 r r r f r + 1 2f 2f + . r2 2 z 2

Expression en coordonn ees sph eriques


2 f = 1 r2 r r2 f r + 1 r2 sin sin f + 1 2f . r2 sin2 2

Relations fondamentales
div (grad) = laplacien div (rot) = 0 rot(rot) = grad(div ) laplacien rot(grad) = 0

Exercices
1. On consid` ere un champ v purement divergent en 2D (a) D eterminer le syst` eme d equations di erentielles v eri e par les composantes du champ v (b) Donner un exemple simple de vecteur v (c) Tracer le champ v 2. Idem mais pour un champ purement rotationnel.

Vous aimerez peut-être aussi