Condizione di Courant-Friedrichs-Lewy

In fluidodinamica numerica la condizione di Courant-Friedrichs-Lewy (spesso abbreviata con CFL) è una condizione necessaria per la convergenza numerica della soluzione di alcune equazioni alle derivate parziali (solitamente 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.

La condizione CFL è comunemente imposta per quei termini delle equazioni alle derivate parziali che rappresentano la convezione (parte iperbolica di una EDP). Per un caso unidimensionale la condizione CFL è scritta come:

$\frac {u \cdot \Delta\,t} {\Delta\,x} < C$

dove

u rappresenta la velocità [L/T]

Δt è l'intervallo temporale [T]

Δx è l'intervallo spaziale [L],

e la costante C dipende dalla particolare equazione che deve essere risolta, non da Δt o da Δx. Il numero

$\nu = \frac {u \cdot \Delta\,t} {\Delta\,x}$

è chiamato numero di Courant.

In un caso bidimensionale questo può scriversi come:

$\frac {u_ x \cdot \Delta\,t} {\Delta\,x} + \frac {u_ y \cdot \Delta\,t} {\Delta\,y} < C.$

La condizione CFL può diventare un grosso limite per il passo temporale Δt a causa del fatto che per certe equazioni alle derivate parziale del quarto ordine non lineari può diventare nella forma:

$\frac{\Delta t}{\Delta x^4} < C,$

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

Note

Portale Fisica

Portale Ingegneria

Portale Matematica