Calcul la vitesse d’un choc – équation de Burgers
Utilisez ce calculateur interactif pour déterminer la vitesse de propagation d’un choc dans l’équation de Burgers inviscide à partir des états gauche et droit d’une discontinuité initiale. L’outil applique la condition de Rankine-Hugoniot et affiche aussi une visualisation du profil de solution.
- Équation de Burgers inviscide
- Vitesse du choc
- Condition de Rankine-Hugoniot
- Visualisation avec Chart.js
Calculateur
Résultats
u_t + (u²/2)_x = 0s = (f(uL) - f(uR)) / (uL - uR) = ((uL²/2) - (uR²/2)) / (uL - uR) = (uL + uR) / 2
Guide expert : comprendre le calcul de la vitesse d’un choc avec l’équation de Burgers
Le calcul de la vitesse d’un choc avec l’équation de Burgers est un sujet classique en analyse des équations aux dérivées partielles, en mécanique des fluides simplifiée et en théorie des lois de conservation. Derrière cette expression se cache une idée très concrète : lorsqu’une grandeur physique se propage de manière non linéaire, des discontinuités peuvent apparaître. Ces discontinuités sont appelées chocs, et leur vitesse de déplacement ne se déduit pas d’une intuition purement géométrique, mais d’une relation de conservation rigoureuse.
L’équation de Burgers inviscide s’écrit sous la forme u_t + (u²/2)_x = 0. Elle est souvent utilisée comme modèle pédagogique parce qu’elle concentre les phénomènes essentiels des lois de conservation hyperboliques : transport non linéaire, croisement de caractéristiques, création de chocs et, dans d’autres configurations, formation d’éventails de raréfaction. En pratique, si l’état à gauche d’une discontinuité initiale est plus grand que l’état à droite, un choc se forme et sa vitesse se calcule grâce à la condition de Rankine-Hugoniot.
Pourquoi l’équation de Burgers est-elle si importante ?
Bien qu’elle soit plus simple que les équations d’Euler de la dynamique des gaz, l’équation de Burgers reproduit un mécanisme central des systèmes non linéaires : la vitesse locale dépend de l’état lui-même. Cela signifie que les zones où la solution est plus élevée peuvent se déplacer plus vite que les zones plus faibles. Petit à petit, le profil se raidit, les caractéristiques se croisent, puis une solution classique cesse d’exister. On remplace alors cette solution lisse par une solution faible qui admet une discontinuité mobile.
Ce cadre mathématique est extrêmement utile pour :
- introduire la notion de choc dans les lois de conservation ;
- comprendre la différence entre solution classique et solution faible ;
- illustrer la nécessité d’une condition d’entropie ;
- tester des schémas numériques de volumes finis ;
- interpréter des phénomènes de compression non linéaire.
La formule exacte de la vitesse du choc
Pour une loi de conservation scalaire sous la forme générale u_t + f(u)_x = 0, la vitesse d’un choc reliant un état gauche uL à un état droit uR est donnée par :
- on identifie le flux : f(u) ;
- on applique la relation de Rankine-Hugoniot ;
- on obtient s = (f(uL) – f(uR)) / (uL – uR).
Dans le cas particulier de Burgers inviscide, le flux vaut f(u) = u²/2. On obtient alors :
s = ((uL²/2) – (uR²/2)) / (uL – uR) = (uL + uR)/2
Cette formule a une interprétation élégante : la vitesse du choc est simplement la moyenne arithmétique des deux états. Par exemple, si uL = 3 et uR = 1, la vitesse du choc est s = 2. Si la discontinuité initiale se trouve en x0 = 0, la position du choc au temps t est :
x(t) = x0 + s t
Quand y a-t-il réellement un choc ?
Il ne suffit pas de brancher deux nombres dans la formule pour parler d’un choc physiquement admissible. Avec l’équation de Burgers, la structure de la solution dépend de l’ordre des états :
- si uL > uR, les caractéristiques convergent et un choc se forme ;
- si uL < uR, les caractéristiques divergent et la solution entropique correcte est une raréfaction ;
- si uL = uR, il n’y a ni choc ni raréfaction : la solution est constante.
C’est une nuance fondamentale. Mathématiquement, la relation de Rankine-Hugoniot fournit une vitesse potentielle de saut, mais la condition d’entropie de Lax décide si ce saut représente bien un choc admissible. Dans Burgers, cette condition revient à exiger que la vitesse des caractéristiques à gauche soit supérieure à la vitesse du choc, elle-même supérieure à la vitesse des caractéristiques à droite.
| Configuration | États | Structure obtenue | Vitesse ou loi |
|---|---|---|---|
| Compression | uL = 3, uR = 1 | Choc admissible | s = (3 + 1)/2 = 2.0 |
| Compression faible | uL = 1.8, uR = 1.2 | Choc admissible | s = 1.5 |
| Équilibre | uL = 2, uR = 2 | Solution constante | Pas de front mobile distinct |
| Détente | uL = 1, uR = 3 | Raréfaction | u(x,t) = (x-x0)/t dans l’éventail |
Interprétation physique du résultat
La vitesse du choc n’est pas choisie arbitrairement. Elle est imposée par le bilan de quantité conservée à travers la discontinuité. Si vous imaginez une petite boîte de contrôle entourant le front, la quantité qui entre et celle qui sort doivent équilibrer le mouvement du front lui-même. C’est exactement ce que traduit la condition de Rankine-Hugoniot.
Dans un langage plus concret, le choc se déplace à une vitesse intermédiaire entre les deux états, ce qui est cohérent avec l’idée d’une transition brutale entre une région plus rapide et une région plus lente. Le choc n’avance donc ni à la vitesse du plateau gauche, ni à celle du plateau droit, mais à une vitesse déterminée par le flux total.
Exemple de calcul détaillé
Prenons un problème de Riemann standard :
- état gauche : uL = 4 ;
- état droit : uR = 0.8 ;
- position initiale : x0 = -1 ;
- temps : t = 3.
Le flux de Burgers vaut f(u)=u²/2. Donc :
- f(uL)=4²/2=8
- f(uR)=0.8²/2=0.32
- s=(8-0.32)/(4-0.8)=7.68/3.2=2.4
La forme simplifiée donne exactement la même chose :
s=(uL+uR)/2=(4+0.8)/2=2.4
La position du choc au temps 3 est alors :
x(3)=x0 + s t = -1 + 2.4 × 3 = 6.2
Le calculateur ci-dessus automatise précisément cette chaîne logique. Il vérifie aussi si la configuration correspond à un choc, à une raréfaction ou à une solution uniforme.
Comparaison entre choc et raréfaction
Dans l’apprentissage des lois de conservation, l’erreur la plus fréquente consiste à appliquer mécaniquement la formule du choc sans vérifier la structure entropique réelle. Or le problème de Riemann pour Burgers possède deux comportements très différents. Le tableau suivant résume les distinctions essentielles.
| Critère | Choc | Raréfaction |
|---|---|---|
| Relation entre états | uL > uR | uL < uR |
| Comportement des caractéristiques | Elles se croisent | Elles s’écartent |
| Régularité de la solution | Discontinuité mobile | Profil continu auto-similaire |
| Loi principale | s = (uL + uR)/2 | u(x,t) = (x-x0)/t dans l’éventail |
| Exemple numérique | uL = 3, uR = 1, s = 2 | uL = 1, uR = 3, bords à x/t = 1 et 3 |
Ordres de grandeur utiles pour l’analyse numérique
En simulation numérique, la vitesse du choc influence directement le choix du pas de temps via une condition de type CFL. Voici quelques valeurs simples souvent utilisées comme cas tests pédagogiques :
- uL = 2, uR = 0 donne s = 1 ;
- uL = 5, uR = 1 donne s = 3 ;
- uL = 1.2, uR = -0.4 donne s = 0.4 ;
- uL = -1, uR = -3 donne s = -2, donc un choc se déplaçant vers la gauche.
On remarque que le choc peut se propager dans les deux sens selon le signe des états. La moyenne des deux valeurs donne immédiatement l’intuition de la direction du front.
Rôle de la viscosité et lien avec Burgers visqueuse
L’équation de Burgers visqueuse ajoute un terme diffusif : u_t + (u²/2)_x = nu u_xx. Quand la viscosité nu est petite mais non nulle, le choc n’est plus strictement discontinu ; il devient une couche de transition très raide mais lisse. Pourtant, dans la limite de viscosité nulle, la position de cette couche converge vers le choc de la version inviscide et sa vitesse reste cohérente avec Rankine-Hugoniot.
C’est pour cela que l’équation de Burgers sert aussi de pont entre :
- les modèles purement hyperboliques ;
- les problèmes avec diffusion ;
- la théorie asymptotique des couches internes ;
- la validation de méthodes numériques pour les fronts raides.
Sources académiques et institutionnelles recommandées
Pour approfondir le sujet, vous pouvez consulter les ressources suivantes :
- University of Maryland (.edu) – notes sur les lois de conservation et les chocs
- University of Pennsylvania (.edu) – notes de cours sur l’équation de Burgers
- NASA (.gov) – introduction institutionnelle aux ondes de choc et à la propagation compressible
Erreurs fréquentes dans le calcul de la vitesse d’un choc
- Confondre la vitesse du choc et la vitesse caractéristique. Pour Burgers, la vitesse caractéristique locale vaut u, alors que la vitesse du choc vaut (uL+uR)/2.
- Oublier la condition d’entropie. Si uL < uR, la solution correcte n’est pas un choc mais une raréfaction.
- Négliger la position initiale x0. La vitesse seule ne suffit pas ; pour connaître la position du front, il faut intégrer dans le temps.
- Choisir un domaine graphique trop petit. Le front peut sortir rapidement de la fenêtre d’observation si le temps ou les états sont grands.
- Interpréter la solution visqueuse comme identique au cas inviscide à toute échelle. La vitesse reste proche, mais l’épaisseur du front change.
Comment utiliser ce calculateur efficacement
Pour obtenir un résultat fiable et bien interprété, suivez cette méthode :
- entrez les valeurs de uL et uR ;
- fixez la position initiale x0 ;
- choisissez le temps t de visualisation ;
- vérifiez que le domaine graphique couvre la position attendue du front ;
- lancez le calcul ;
- interprétez le message : choc admissible, raréfaction ou état constant ;
- analysez enfin le graphique généré.
Conclusion
Le calcul de la vitesse d’un choc avec l’équation de Burgers repose sur une idée simple mais fondamentale : toute discontinuité admissible dans une loi de conservation doit respecter un bilan de flux. Pour Burgers inviscide, ce bilan conduit à une formule remarquablement compacte, s = (uL + uR)/2. Cette simplicité en fait un excellent laboratoire mathématique pour comprendre les chocs, les rarefactions, la condition d’entropie et les bases de l’analyse numérique hyperbolique.
Si vous cherchez un outil pratique pour tester des cas de figure, comparer plusieurs configurations ou illustrer un cours, le calculateur ci-dessus fournit une base robuste et visuelle. Il ne se contente pas de donner un nombre ; il vous aide à interpréter la nature de la solution et à visualiser l’évolution spatiale du problème.