Calcul La Vitesse D Un Choc Equation Burgers

Calculez instantanément la vitesse d’une onde de choc pour l’équation de Burgers inviscide, identifiez le régime physique obtenu et visualisez le profil de solution à un instant donné grâce à un graphique interactif.

Calculateur interactif

Formule utilisée Equation de Burgers: u_t + (u²/2)_x = 0 Vitesse du choc (Rankine-Hugoniot): s = (f(uL) - f(uR)) / (uL - uR) = (uL + uR) / 2

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 cas d’école fondamental en mécanique des fluides, en mathématiques appliquées et en analyse numérique. Même si l’équation de Burgers est souvent présentée comme un modèle simplifié, elle concentre plusieurs phénomènes majeurs des lois de conservation hyperboliques : transport non linéaire, croisement des caractéristiques, apparition de discontinuités et sélection d’une solution entropique. En pratique, elle sert d’outil pédagogique puissant pour comprendre ce qui se passe dans des modèles plus complexes comme les équations d’Euler compressibles, les écoulements supersoniques, certains problèmes de trafic routier ou des schémas numériques de propagation d’ondes.

L’équation considérée ici est la forme conservative :

u_t + (u²/2)_x = 0

Dans cette écriture, le flux est f(u) = u² / 2. Lorsqu’on impose une donnée de Riemann, c’est-à-dire une condition initiale avec une valeur constante à gauche et une autre à droite, la dynamique dépend de la comparaison entre uL et uR. Si uL > uR, les caractéristiques convergent et un choc se forme ou persiste. Si uL < uR, les caractéristiques divergent et on observe une onde de détente auto-similaire. Cette distinction est capitale : beaucoup d’erreurs de calcul viennent du fait qu’on applique la formule d’un choc à une situation qui est en réalité une détente.

La formule exacte de la vitesse du choc

Pour une loi de conservation sous la forme u_t + f(u)_x = 0, la vitesse de propagation d’un choc satisfaisant les états constants uL et uR est donnée par la relation de Rankine-Hugoniot :

1. On identifie le flux f(u).
2. On calcule la différence de flux f(uL) – f(uR).
3. On divise par le saut d’état uL – uR.

Pour Burgers, cela donne :

s = (f(uL) – f(uR)) / (uL – uR) = ((uL²/2) – (uR²/2)) / (uL – uR) = (uL + uR) / 2

Cette formule est remarquablement simple : la vitesse du choc est la moyenne arithmétique des états gauche et droit. Si par exemple uL = 3 et uR = 1, alors la vitesse est s = 2. Si le choc est initialement en x0 = 0, sa position au temps t = 2 devient x = x0 + s t = 4.

Pourquoi cette formule est-elle physiquement cohérente ?

Dans l’équation de Burgers, la vitesse caractéristique vaut f'(u) = u. Autrement dit, l’information associée à une valeur locale u se propage à la vitesse u. Quand uL > uR, les caractéristiques provenant de la gauche vont plus vite que celles de la droite. Elles finissent donc par se rencontrer. Ce croisement rend la solution classique multivaluée si l’on persiste à suivre uniquement les caractéristiques. La théorie des solutions faibles remplace alors cette configuration impossible par une discontinuité mobile : le choc.

La vitesse s = (uL + uR) / 2 se situe exactement entre les deux vitesses caractéristiques associées aux états extrêmes. D’ailleurs, la condition d’entropie de Lax pour un choc de Burgers s’écrit :

• uL > s > uR

En remplaçant s par (uL + uR)/2, on voit immédiatement que cette condition est satisfaite dès que uL > uR. Cela confirme que la discontinuité retenue est physiquement admissible.

Méthode pratique pour calculer la vitesse d’un choc

Voici la procédure opérationnelle utilisée par le calculateur ci-dessus.

1. Entrer les deux états constants uL et uR.
2. Vérifier le régime :
   • si uL > uR, il s’agit d’un choc ;
   • si uL = uR, la solution reste uniforme ;
   • si uL < uR, il s’agit d’une détente.
3. En régime de choc, calculer s = (uL + uR)/2.
4. Calculer la position au temps t par x(t) = x0 + s t.
5. Tracer le profil en espace pour visualiser la rupture mobile.

Cette méthodologie est la base de très nombreux solveurs de Riemann utilisés dans les schémas de volumes finis. Même lorsqu’on travaille ensuite sur des systèmes plus complexes, cette logique de conservation locale et de vitesse de propagation reste centrale.

Exemples numériques de référence

Cas	uL	uR	Type de solution	Vitesse calculée	Position à t = 2 si x0 = 0
Benchmark A	3.0	1.0	Choc	2.0	4.0
Benchmark B	2.5	-0.5	Choc	1.0	2.0
Benchmark C	1.0	1.0	État constant	1.0	Pas de rupture à suivre
Benchmark D	0.5	2.0	Détente	Non applicable comme choc	Profil éventail

Le tableau précédent illustre un point essentiel : la même formule de moyenne peut être évaluée numériquement, mais elle n’a un sens de vitesse de choc admissible que lorsque la structure est réellement compressive, donc lorsque uL > uR. Dans le cas d’une détente, la physique de la solution est différente et le profil n’est pas une marche abrupte mobile, mais un éventail continu reliant les deux états.

Interprétation du profil : choc, état constant ou détente

1. Cas du choc

Si uL > uR, le profil à un temps donné est une marche : à gauche de la position du choc, la solution vaut uL, et à droite elle vaut uR. Le calculateur affiche alors la vitesse, la position et un graphique net avec changement brutal de niveau.

2. Cas uniforme

Si uL = uR, aucune onde particulière ne se forme. La solution reste constante partout. En termes numériques, c’est un cas utile pour tester la cohérence d’un code : aucun artefact ne devrait apparaître.

3. Cas de la détente

Si uL < uR, les caractéristiques divergent. La solution entropique prend alors la forme d’un éventail auto-similaire. Pour l’équation de Burgers, on écrit classiquement :

• u(x,t) = uL si (x – x0)/t < uL
• u(x,t) = (x – x0)/t si uL <= (x – x0)/t <= uR
• u(x,t) = uR si (x – x0)/t > uR

C’est précisément pour cette raison qu’un bon outil de calcul ne doit pas seulement fournir une valeur, mais aussi diagnostiquer le régime. Une réponse numérique isolée peut être trompeuse sans contexte analytique.

Applications concrètes et ordres de grandeur utiles

Bien que l’équation de Burgers soit plus simple que les modèles complets de la dynamique des gaz, elle est extrêmement utile pour représenter des mécanismes de concentration d’ondes, de formation de fronts et de transport non linéaire. On la rencontre dans :

• la modélisation pédagogique des ondes de choc en aérodynamique ;
• l’analyse de schémas numériques pour lois de conservation ;
• la propagation de fronts dans des systèmes diffusifs et convectifs ;
• des analogies avec le trafic routier et les congestions mobiles ;
• l’étude théorique des solutions faibles et des entropies.

Domaine	Indicateur mesuré	Valeur typique	Intérêt pour Burgers
Trafic autoroutier	Vitesse libre d’un véhicule léger	110 à 130 km/h	Illustre la propagation d’information et les fronts de densité
Écoulement supersonique	Mach de vol de chasse ou d’essai	1.2 à 2.0+	Met en évidence la formation de chocs compressifs
Calcul scientifique	Nombre de cellules d’un test 1D standard	100 à 1000	Ordre de grandeur courant pour visualiser diffusion numérique et capture de choc
Benchmark académique Burgers	Nombre de points de tracé interactif	100 à 200	Compromis efficace entre lisibilité et performance navigateur

Ces ordres de grandeur ne signifient pas que Burgers remplace un modèle réaliste complet. Ils montrent plutôt pourquoi ce cadre est si précieux pour l’intuition : il permet d’isoler le mécanisme clé sans tout le poids d’un système multiphysique plus difficile à résoudre.

Erreurs fréquentes lors du calcul

• Confondre vitesse caractéristique et vitesse de choc : pour Burgers, la vitesse caractéristique vaut u, alors que la vitesse de choc vaut (uL + uR)/2.
• Oublier la condition d’entropie : une discontinuité avec uL < uR n’est pas un choc admissible.
• Utiliser la forme non conservative sans précaution : près des discontinuités, la forme conservative est indispensable.
• Négliger la position initiale x0 : la vitesse seule ne suffit pas, il faut aussi localiser le front.
• Tracer un graphique sans bornes adaptées : si la fenêtre spatiale est trop petite, on peut croire à tort que le choc a disparu.

Pourquoi Burgers reste un modèle de référence en 2025

Dans l’enseignement supérieur, la recherche en méthodes numériques et les démonstrations de solveurs de Riemann, Burgers demeure incontournable. Sa valeur tient à quatre raisons. Premièrement, il présente une non-linéarité réelle. Deuxièmement, il produit des chocs exacts et des détentes exactes. Troisièmement, la relation de Rankine-Hugoniot y est simple et lisible. Quatrièmement, il permet de comparer facilement solution analytique et approximation numérique. C’est pour cela qu’on le retrouve encore dans les cours de calcul scientifique, de dynamique des fluides numérique et de théorie des équations aux dérivées partielles.

Ressources d’autorité pour approfondir

• NASA.gov – introduction aux ondes de choc normales
• MIT.edu – cours de méthodes numériques et EDP avec lois de conservation
• Wisc.edu – ressources universitaires sur les lois de conservation hyperboliques

Résumé opérationnel

Si vous devez retenir une seule idée, c’est celle-ci : pour l’équation de Burgers u_t + (u²/2)_x = 0, la vitesse d’un choc entre deux états constants est s = (uL + uR)/2, mais seulement dans le régime uL > uR. Ensuite, la position se calcule par x(t) = x0 + s t. Cette simplicité analytique fait de Burgers un laboratoire idéal pour comprendre la naissance, la propagation et la sélection physique des chocs dans les lois de conservation.

Le calculateur ci-dessus applique automatiquement cette logique, affiche le type de solution détecté et génère un graphique interactif pour vous permettre de vérifier immédiatement la cohérence du résultat.