Essaye de refaire ton calcul avec f(r)=r. Tu verras que ton erreur se situe dans le fait que tu n’as dérivé que dans une seule direction et pas trois (d’où le facteur trois). Le fait que ça se factorise vient des symétries par rapport à x,y et z.
Dis autrement, et peut-être de manière plus intuitive,
div(rf(r))=f(r)divr+rgradf(r)
Tu retrouves bien ta dérivée (via le gradient), et le 3 vient du fait que la dérivée de la fonction r -> r vaut 1 par rapport à x, à y et à z. Le fameux « tu n’as dérivé que dans une seule direction et pas trois » d’Holosmos vient de là.
Tu doit utiliser les coordonnées sphériques. La tu donnes la formule de divergence pour le cylindrique (ce qui doit te conduire à une mauvaise réponse ^^).
Je te donne une interprétation de la formule (à la physicienne) :
La divergence donne en chaque point de l’espace une valeur qui représente la somme de ce qui "rentre" et de ce qui "sort" d’une petit élément de volume. En effet on l’applique sur un champ vectoriel. Pour se le représenter on peut par exemple imaginer une rivière et le champ vectoriel représente la vitesse de l’eau. En chaque petite élément de volume de rivière il y a de l’eau qui rentre et de l’eau qui sort. On voit que si l’eau à une vitesse plus importante à l’entrée du volume qu’à la sortie alors nécessairement la quantité d’eau dans ce volume doit augmenter. Si on calcul la divergence de la vitesse à cet endroit elle sera négative. Inversement si la vitesse est plus faible l’entre qu’a la sortie alors l’eau quitte le volume => divergence positive.
En cartésien ces petites éléments de volume sont décrit par des cubes de coté dx, dy et dz. En cylindrique ou sphérique c’est plus vraiment des cubes, certains cotés sont des surfaces de cylindre (ou de sphère). Du coup ça complique un peu le calcul.
Par exemple en coordonnée cylindrique les éléments de volume sont comme ceci :
Image utilisateur
Pour retrouver "avec les mains" l’opérateur divergence il faut faire la somme de ce qu’il rentre et sort de ces "cubes". On doit donc faire la somme de ce qu’il rentre par "la face 1" moins ce qu’il sort par la face opposé (disons "la face 4") plus ce qu’il rentre par la face 2 moins ce qu’il sort par la face opposé : la face 5 plus ce qu’il rentre par la face 3 moins ce qu’il sors par la face opposé : la face 6
Le premier truc à voir est que uniquement la composante orthogonal du champ de vecteur ( ie de la "vitesse") à la face du "cube" élémentaire participe à faire rentrer ou sortir du "fluide". Ainsi seul la composante Fx participera à remplir le cube à travers les deux faces orthogonales à x^, seul la composante Fy participera à remplir le cube à travers les deux faces orthogonales à y^ et bien sur seul la composante Fz participera à remplir le cube à travers les deux faces orthogonal à z^.
Et la c’est pareil en cylindrique ( ou sphérique) :
seul la composante Fr participera à remplir le cube à travers les deux faces orthogonal à r^, etc.
Essayons de retrouver le terme en r de la divergence en cylindrique (tu pourras le faire en sphérique et voir que tu obtiens bien le bon résultat ) :
Il suffit de sommer ce qui rentre/sort par les deux faces orthogonal à r^ afin d’avoir le "terme en r".
Ce qui rentre par la face c’est le produit de la composante orthogonal à la face de ta fonction vectoriel multiplié par la surface de ta face (on peut l’exprimer sous forme de produit scalaire mais c’est pas important ici).
Il y a un truc auquel il faut faire attention : c’est le volume du cube élémentaire. En cartésien le volume du cube est constant car la dimension de ses arêtes est constante : V=dxdydz mais en cylindrique (ou sphérique) ce n’est plus le cas. C’est pour cela que la formule est un peu plus compliqué.
En effet on a en cylindrique un "cube" elementaire qui a un volume dV=dr×rdθ×dz. En d’autre terme plus on regarde loin de l’axe Oz (r grand) plus le cube à un volume important. Ainsi la divergence doit être normalisée par le volume du cube. C’est assez logique puisque ce qui nous intéresse c’est "ce qui rentre" par élément de volume identique ("par metre cube" par exemple). Sans cette normalisation la valeur obtenu n’a plus aucun sens et correspondra proche de l’axe Oz (r petit) à "ce qui rentre" dans un tous petit cube et correspondra loin de l’axe Oz (r grand) à "ce qui rentre" dans une cube gigantesque… ça n’aurait alors pas de sens de comparer ces deux valeurs. Normaliser par le volume du cube permet de rendre ces valeurs comparable.
Ainsi on divise la somme de ce qui rentre et sort dans les petits cube par leur petit volume et on divise la somme de ce qui rentre et ce qui sort des grands cube par leur grand volume afin d’avoir des valeurs qui puissent être comparées.
Concrètement on écrit cela comme : divr^=rdrdθdzFr(r+dr,θ,z)(r+dr)dθdz−Fr(r,θ,z)rdθdz
Et si on re-arrange un peut l’expression on voit : r1(rdrFr(r+dr,θ,z)−Fr(r,θ,z)+Fr(r+dr,θ,z))
Et bien sur on evalue cette expression pour des cube infinitésimaux cad dr tend vers zero :
edit : mais je sais pas si c’est très clair. En me relisant j’ai comme un doute…
Tu peux essayer le même raisonnement pour les autres faces de l’élément de volume en cylindrique : les faces orthogonales a θ^ et z^ et pour les faces de l’élément de volume en sphériques.
Tu devrais retrouver l’expression de la divergence dans ces coordonnées !
edit2 : correction des erreurs de signe et des d ronds
Tu doit utiliser les coordonnées sphériques. La tu donnes la formule de divergence pour le cylindrique (ce qui doit te conduire à une mauvaise réponse ^^).
Je ne suis pas certain, je pense que quel que soit le repère utilisé, le champ vectoriel F reste le même.
Après, j’obtiens 2f(r)+rf′(r) comme résultat avec les coordonnées cylindriques, je ne comprends pas pourquoi cela ne marche pas…
Plus sérieusement, c’est à cause des faces "d’entrées" et de "sorties" et des volumes élémentaire dont la surface/le volume n’évoluent pas de la même maniere avec r en sphérique et en cylindrique:
En cylindrique on a :
Face proche de l’axe : dS1=rdθdz
Face loin de l’axe : dS2=(r+dr)dθdz
Volume : dV=dr×rdθ×dz=rdrdθdz
En sphérique on a (tu trouveras sur google image des dessins similaire a celui que j’ai montré pour les coordonnée sphérique pour t’aider à réfléchir)
Face proche de l’origine : dS1=rdθrsin(θ)dϕ=r2sin(θ)dθdϕ
Face loin de l’origine : dS2=(r+dr)2sin(θ)dθdϕ
Volume : dV=dr×rdθ×dz=r2sin(θ)drdθdϕ
Si tu refais le calcul que j’ai fais tu verras que tu n’obtiens pas le même résultat en sphérique qu’en cylindrique. Et c’est bien confirmé par la page wiki. (ouf !)
Également : quand je parle de somme sur les faces dans mon post précédant, ce sont des sommes algébrique c’est pour cela que j’ai compté le premier terme en positif et le second en négatif.
En y réfléchissant un peu c’est logique : pour une vitesse toujours orienté positivement selon l’axe r^ dans une face "ca rentre" (donc +) et sur la face opposé "ca sort" (donc − ).
Du coup, je pense que le calcul est quand même possible en coordonnées cylindriques.
Le calcul est possible dans n’importe quel système de coordonnées arbitraire (enfin, tant que tu peux calculer son Jacobien), c’est juste que tu vas te fatiguer pour rien à recalculer un changement de repère.
@Ludwig : comme l’as dit adri1 le calcul est évidement possible dans n’importe quelle systeme de coordonnée et heureusement (sinon il y a un problème) !
Mais en fait, ce que tu fais en passant par le cylindrique c’est ce compliquer la tache (mais très bien pour apprendre à manier les dérivées ). Dans ton exercice on te donne directement une expression de r comme le r des coordonnées sphérique. Il y a juste a appliquer la formule des coordonnée sphériques.
Mais comme tu le dis :
Et bon… faute de pouvoir interpréter cette formule, j’ai terminé mon exercice sans faire d’histoires.
Au début on trouve forcement ces formules un peu obscure. Elle ne sont pas si évidente à interpréter. Et on a tendance a faire comme tu as fait au début et penser que ça se généralise comme ceci : divf=∂r∂fr+∂θ∂fθ+∂ϕ∂fϕ ?
D’où ma réponse qui explique pourquoi ce n’est pas si simple : il faut prendre en compte la déformation lié au système de coordonnée que tu utilises.
Et egalement pour essayer de donner un peu de substance à cet opérateur qu’est la divergence.
Partant de la formule de la divergence pour un vecteur exprimé en coordonnées cartésiennes, comment peut-on facilement retrouver l’expression de la divergence en coordonnées sphériques ?
Connectez-vous pour pouvoir poster un message.
Connexion
Pas encore membre ?
Créez un compte en une minute pour profiter pleinement de toutes les fonctionnalités de Zeste de Savoir. Ici, tout est gratuit et sans publicité.
Créer un compte