Condizione di Courant-Friedrichs-Lewy

In fluidodinamica numerica, la condizione di Courant-Friedrichs-Lewy, spesso abbreviata con CFL ed il cui nome è dovuto a Richard Courant, Kurt Friedrichs e Hans Lewy, è una condizione necessaria per la convergenza numerica della soluzione di alcune equazioni alle derivate parziali (di solito, equazioni di tipo iperbolico) ricavata nel 1928.

Questa condizione è sfruttata nell'impiego di schemi numerici espliciti temporali. Come conseguenza, il passo temporale deve essere più piccolo di un certo intervallo di tempo, altrimenti la simulazione produrrà risultati ampiamente scorretti. Per esempio, se un'onda attraversa una griglia di calcolo discreta, allora l'intervallo temporale deve essere più piccolo del tempo necessario all'onda per attraversare due punti adiacenti della griglia. Come corollario, se la distanza tra due punti adiacenti della griglia viene ridotta, il limite superiore dell'intervallo temporale sarà anch'esso diminuito. In sostanza, il dominio numerico (o discreto) di dipendenza deve includere il dominio analitico (o continuo) di dipendenza per poter assicurare che lo schema possa trovare l'informazione necessaria per creare la soluzione.

Descrizione

La condizione CFL è comunemente imposta per quei termini delle equazioni alle derivate parziali che rappresentano la convezione (o meglio, per i termini advettivi, relativi, cioè, ai moti orizzontali o prevalentemente orizzontali). Per un caso unidimensionale la condizione CFL è scritta come:

$C=\frac{u\cdot \mathrm{\Delta }t}{\mathrm{\Delta }x}<C_{max}$

dove $u$ rappresenta la velocità di flusso, $\mathrm{\Delta }t$ è l'intervallo temporale e $\mathrm{\Delta }x$ è l'intervallo spaziale. La costante $C_{max}$ dipende dalla tipologia di equazione che deve essere risolta e dal tipo di schema numerico utilizzato per la soluzione (esplicito o implicito). Se si utilizza uno schema esplicito allora $C_{max}$è dell'ordine di 1. Gli schemi impliciti invece sono meno sensibili alle instabilità numeriche, dunque sono tollerati valori di $C_{max}$più elevati.

Il numero adimensionale $C$ è chiamato numero di Courant.

In un caso bidimensionale la condizione CFL può scriversi come:

$C=\frac{u_x\cdot \mathrm{\Delta }t}{\mathrm{\Delta }x}+\frac{u_y\cdot \mathrm{\Delta }t}{\mathrm{\Delta }y}<C_{max}$

Nel caso $n$-dimensionale si ha:

$C=\mathrm{\Delta }t\sum _{i=1}^n\frac{u_{x_i}}{\mathrm{\Delta }{x}_i}\le C_{max}$

La condizione CFL può diventare un grosso limite per il passo temporale $\mathrm{\Delta }t$ a causa del fatto che per certe equazioni alle derivate parziali del quarto ordine non lineari può diventare nella forma:

$C=\frac{\mathrm{\Delta }t}{\mathrm{\Delta }{x}^4}<C_{max}$

e per questo motivo ci si sforza, in questi casi, di evitare questa condizione usando metodi numerici impliciti.

Bibliografia

• Template:De R. Courant, K.O. Friedrichs, H. Lewy, Ueber die partiellen Differenzgleichungen der mathematische Physik Math Ann. , 100 (1928) pp. 32–74
• Template:En S.K. Godunov, V.S. Ryaben'kii, The theory of difference schemes , North-Holland (1964)
• Template:En Courant, R.; Friedrichs, K.; and Lewy, H. On the Partial Difference Equations of Mathematical Physics. IBM J. 11, 215-234, 1967.

Voci correlate

• Convezione
• Equazione differenziale alle derivate parziali

Collegamenti esterni

• Template:Collegamenti esterni
• Template:SpringerEOM