Comment calculer nombre Courant CFL stabilité CFD numérique

Calcule nombre Courant CFL C = u*dt/dx critère stabilité schemas numériques explicites CFD computational fluid dynamics. C<=1 stable; C>1 divergence. Richard Courant + Friedrichs + Lewy 1928. CNRS + Polytechnique CFD.

Paramètres de calcul
Ce calculateur vous a-t-il été utile ?

Comment calculer le nombre de Courant (C ou CFL Courant-Friedrichs-Lewy) ratio adimensionnel fondamental analyse numérique et dynamique fluides computationnelle CFD mesurant la stabilité des schemas numériques explicites pour équations hyperboliques et de transport: C = u * dt / dx. u vitesse caracteristique m/s (air CFD compressible 100-1000; eau incompressible 0.1-100; transport substances 0.001-1; climat atmospherique 30 m/s vent; oceanographie 0.1-1; meteorologie 10-50); dt pas temporel s (CFD aerospace 1e-7 a 1e-3; meteorologie 60-3600); dx taille maillage m (CFD aerodynamique 0.0001-0.01; meteo 1000-50000).

Critère CFL: C <= 1 condition nécessaire stabilité schemas explicites; C <= 0.5 marge securite industrielle; C <= 0.1 haute précision.

Pas maximum: dt_max = dx/u pour C=1 + dt_max = dx/(|u|+a) compressible vitesse son a + multidimensionnel C = u*dt/dx + v*dt/dy + w*dt/dz <= 1.

Schemas numériques: FTBS upwind first-order stable C<=1; FTCS centre toujours instable advection; Lax-Friedrichs 1954 dissipatif + stable C<=1; Lax-Wendroff 1960 second-order + stable C<=1; MacCormack 1969 predictor-corrector + stable C<=1; Crank-Nicolson 1947 implicit inconditionnellement stable; Beam-Warming 1976 implicit.

Extensions CFL: 2D advection C = u*dt/dx + v*dt/dy <= 1; 3D Navier-Stokes + diffusion Fo = nu*dt/dx^2 <= 0.5; compressible Euler C = (|u|+a)*dt/dx.

Équations couvertes CFL: advection ut + a*ux = 0; wave équation utt = c^2*uxx; Burgers; Navier-Stokes; shallow water + WRF MPAS; oceanographique POP MOM NEMO; sismique; plasma MHD. Histoire: Richard Courant (Lublinitz 1888 - New Rochelle 1972) Gottingen 1908-1933 + NYU 1936-1958 fondateur Courant Institute Mathematical Sciences NYU 1965; 1928 Math Ann 100:32 avec Friedrichs et Lewy classique CFL condition + classical référence numerical analysis + 1924-1937 'Methoden der mathematischen Physik' Springer 2 vol classique avec David Hilbert; pionnier numerical PDEs + Courant Institute world-leading + Gauss Medal 1965 + National Medal Science 1965.

Kurt Friedrichs (1901-1982) NYU 1937-1982; 1928 CFL + 1948 'Supersonic Flow and Shock Waves' classique + Steele Prize 1976.

Hans Lewy (1904-1988) UC Berkeley 1935-1988; 1928 CFL + Lewy example 1957 + Wolf Prize 1984.

John von Neumann (1903-1957) Princeton-Los Alamos; analyse stabilité 1944-1947 + ENIAC numerical weather prediction 1950 + game theory + computer architecture + National Medal 1956.

Peter Lax (Budapest 1926 - NYC 2025) NYU Courant 1958-1995; 1954 Lax-Friedrichs + 1960 Lax-Wendroff classique + 1964 Lax theorem + Wolf Prize 1987 + Abel Prize 2005. Heritage: LeVeque 1992 + 2002 Cambridge 558 + Toro 2009 Springer 724 + Hirsch 2007 + Anderson 1995 + Ferziger-Peric 2020 + Strikwerda 2004.

France recherche CFD + numérique: CNRS INSU INSIS Paris 33000 + INRIA + LJLL Sorbonne + CMAP Polytechnique + LMD; INRIA Paris + Saclay + Grenoble + Sophia + Rennes 4500 world-leading + ANR; CEA Saclay DRF 20000; ONERA Palaiseau Toulouse 2000 + CFD aerospace world-leading + ELSA elsA; IFP Énergies Nouvelles Rueil 1700 + reservoir CFD; BRGM Orleans 1100; Meteo France 2700 + ARPEGE AROME numerical weather prediction; SHOM Brest 500; CSTB Marne-la-Vallee 900; CERFACS Toulouse 100 + climate ocean; GENCI Paris 50 supercomputing France.

France universites CFD + numérique: Sorbonne Paris 6 LJLL Laboratoire Jacques-Louis Lions 50000 + Yvon Maday + Frederic Hecht FreeFEM world-leading + Frederic Coquel + Olivier Pironneau; Ecole Polytechnique CMAP Palaiseau 3000 + Gregoire Allaire + Habib Ammari + Antonin Chambolle; ENS Paris 4000; Universite Paris-Saclay LMO 65000; ENPC Marne-la-Vallee 2000; Mines ParisTech CEMEF Sophia 1500; Universite Aix-Marseille 80000; INSA Lyon 6500; INSA Rouen 1700; Universite Cote Azur Acumes Sophia 35000; Universite Bordeaux 50000; Universite Toulouse III IMT 32000; Universite Strasbourg 50000.

France industrie CFD aerospace + automotive: Airbus Toulouse Hambourg 134000 + A350 A380 + CFD aerodynamique world-leading; Safran Paris 92000 + jet engines CFD + Snecma Turbomeca; Dassault Aviation Paris 12500 + Rafale Falcon + CFD; Thales Alenia Space France Cannes 8500; ArianeGroup 8500 + Ariane 6 CFD; Stellantis France PSA 80000 + Peugeot Citroen CFD; Renault Paris 153000 + CFD; Michelin Clermont-Ferrand 132000 + tire CFD; Faurecia Nanterre 114000; Valeo Paris 109000 + automotive electronics.

France industrie CFD énergie + petrole: Total Énergies La Defense 105000 + reservoir CFD + offshore; EDF 100000 + CFD nuclear + EDF R&D Chatou; Engie Courbevoie 100000 + offshore wind CFD; Naval Group Paris 16000 + sous-marins SNLE CFD; Framatome 14000 + nuclear CFD.

France industrie meteo + supercomputing: Meteo France 2700 + ARPEGE AROME ALADIN; GENCI Paris 50 + Jean Zay supercomputer; Atos Eviden Paris 109000 + supercomputing; CERFACS Toulouse 100.

Formule Courant CFL

Définition: C = u*dt/dx. u vitesse + dt pas temporel + dx maillage + adimensionnel.

Pas maximum: dt_max = dx/u + compressible dt_max = dx/(|u|+a).

Multidimensionnel: C = sum(u_i*dt/dx_i) <= 1.

Diffusion Fourier: Fo = nu*dt/dx^2 <= 0.5.

Valeurs typiques: aerospace C=1 limite; eau C=1; climat C=0.36 stable.

Histoire et contexte

Richard Courant (Lublinitz 1888 - New Rochelle 1972) Gottingen + NYU 1936-1958 fondateur Courant Institute 1965; 1928 Math Ann 100:32 avec Friedrichs Lewy classique CFL + 'Methoden mathematischen Physik' Hilbert; pionnier numerical PDEs + Courant Institute world-leading + Gauss Medal 1965 + National Medal Science 1965.

Kurt Friedrichs (1901-1982) NYU; 1928 CFL + 1948 supersonic + Steele Prize 1976.

Hans Lewy (1904-1988) Berkeley + Wolf Prize 1984.

John von Neumann (1903-1957) Princeton-Los Alamos + ENIAC weather 1950 + Manhattan + game theory + National Medal 1956.

Peter Lax (1926-2025) NYU Courant + Lax-Friedrichs 1954 + Lax-Wendroff 1960 + Wolf 1987 + Abel 2005. Heritage: LeVeque + Toro + Hirsch + Anderson + Strikwerda.

France recherche: CNRS 33000; INRIA 4500 world-leading; CEA Saclay 20000; ONERA 2000 CFD aerospace world-leading + elsA; IFP Énergies 1700; Meteo France 2700 ARPEGE AROME; CSTB 900; CERFACS Toulouse 100; GENCI 50 Jean Zay supercomputer.

Universites: Sorbonne LJLL 50000 + Yvon Maday + Frederic Hecht FreeFEM world-leading; Polytechnique CMAP 3000 + Allaire + Pironneau; Paris-Saclay 65000; ENPC 2000; Mines ParisTech CEMEF 1500; Aix-Marseille 80000; INSA Lyon 6500; Cote Azur Acumes 35000; Bordeaux 50000; Toulouse III IMT 32000.

Industrie CFD: Airbus 134000 A350 world-leading; Safran 92000 jet engines; Dassault Aviation 12500; ArianeGroup 8500 Ariane 6; Stellantis France 80000; Renault 153000; Michelin 132000; Valeo 109000.

Énergie: Total 105000 reservoir; EDF 100000 nuclear; Engie 100000 offshore wind; Naval Group 16000; Framatome 14000.

Meteo: Meteo France 2700 ARPEGE; Atos Eviden 109000 supercomputing.

Tableau Courant

  • C = u*dt/dx
  • C <= 0.5: stable marge
  • C <= 1: limite CFL
  • C > 1: instable
  • dt_max = dx/u
  • Fo <= 0.5: diffusion

5 exemples numériques

  1. Aerospace: u=100 + dt=1e-5 + dx=0.001 + C=1 + Safran ONERA elsA.
  2. Eau canal: u=1 + dt=0.01 + C=1 + Sorbonne LJLL FreeFEM.
  3. Instable: u=100 + dt=0.001 + dx=0.05 + C=2.
  4. Industriel: u=10 + dt=0.001 + C=0.1 + Renault.
  5. Climat: u=30 + dt=60 + dx=5000 + C=0.36 + Meteo France ARPEGE.

Limites

Linear stability + Boundary GKS + Implicit + Adaptive.

Erreurs courantes

  • Confondre C et Fo.
  • Compressible vitesse son.
  • 2D 3D extension.
  • Implicit damping.
  • Boundary GKS.

Applications

  1. Aerospace: Airbus + Safran world-leading + ONERA elsA.
  2. Automotive: Stellantis + Renault + Michelin.
  3. Naval: Naval Group sous-marins.
  4. Climat: Meteo France ARPEGE + AROME.
  5. Oceanographie: NEMO Ifremer.
  6. Petrole: Total reservoir + IFPEN.
  7. Sismique: BRGM.
  8. Plasma: CEA WEST.
  9. Acoustique: ONERA.
  10. Combustion: Safran + ONERA.

Mis à jour en juillet 2026

Questions fréquentes

Que mesure le nombre Courant CFL?

Nombre Courant C = u*dt/dx mesure le ratio adimensionnel qui determine la stabilité des schemas numériques explicites pour équations hyperboliques et de transport en CFD computational fluid dynamics. C1 schemas explicites divergent. dt_max = dx/u pour explicite + Crank-Nicolson 1947 implicit inconditionnellement stable. Richard Courant + Friedrichs + Lewy 1928 paper Math Ann 100:32 founding paper.

Qu'est-ce que la condition CFL Courant-Friedrichs-Lewy 1928?

Condition CFL Courant-Friedrichs-Lewy 1928 paper Math Ann 100:32: pour schemas explicites finite differences sur équations hyperboliques le domaine de dependance numérique doit contenir le domaine de dependance physique + sinon instabilite + divergence exponentielle. Nécessaire mais pas suffisante en général. Pionnier Courant Institute NYU + classical référence. Strikwerda 2004 SIAM 434 pages classique.

Difference entre C et nombre Fourier Fo?

C = u*dt/dx critère advection (transport) + Fo = nu*dt/dx^2 critère diffusion (viscosité). Tous deux nécessaires pour Navier-Stokes: C <= 1 explicit advection + Fo <= 0.5 explicit diffusion. Combinee C + Fo <= 1 condition nécessaire. Strikwerda 2004 distinction classique + Anderson 'CFD' 1995 McGraw-Hill 547 pages classique textbook + classical référence numerical PDEs.

Schemas explicites vs implicites avantages?

Schemas explicites (FTBS upwind + Lax-Wendroff + MacCormack): faible cout par time step + faciles a implementer + parallelisables + mais contraints C<=1; schemas implicites (Crank-Nicolson 1947 + Beam-Warming 1976): inconditionnellement stables + permettent dt grand + mais cout eleve par time step + système linéaire chaque pas. Trade-off + classical référence Anderson 1995 + LeVeque 2002 Cambridge 558 pages.

Qu'est-ce que l'analyse stabilité von Neumann?

Analyse stabilité von Neumann 1944-1947 Manhattan Project Los Alamos: technique d'analyse spectrale via decomposition de Fourier pour determiner la stabilité des schemas numériques linéaires + amplification factor |g| <= 1 + condition nécessaire. John von Neumann 1903-1957 Hongrie-USA Princeton-Los Alamos + ENIAC numerical weather prediction 1950 avec Charney Fjortoft + classical référence Strikwerda 2004 + Anderson 1995.

Théorème de Lax que dit-il?

Théorème de Lax 1964 fondamental numerical analysis: pour schemas numériques d'équations aux dérivées partielles bien posees consistance + stabilité + convergence sont equivalentes. Peter Lax 1926-2025 Hongrie-USA NYU Courant Institute + Wolf Prize 1987 + Abel Prize 2005 + 1954 Lax-Friedrichs scheme + 1960 Lax-Wendroff classical référence + LeVeque 2002 Cambridge + Strikwerda 2004 SIAM classique.

Lax-Wendroff vs upwind lequel choisir?

Lax-Wendroff 1960 second-order precis + dispersif oscillations pres discontinuites chocs + adapte smooth solutions; upwind first-order 1952 stable + dissipatif + monotonique + adapte shock waves sans oscillations + moins precis. Schemas TVD Total Variation Diminishing Harten 1983 combinent: high-order smooth + first-order shock + Toro 'Riemann Solvers' 2009 Springer 724 pages classique + LeVeque 2002 classical référence.

Commentaires (0)

Pas encore de commentaire. Soyez le premier !