Metodo dei volumi finiti

Il metodo numerico dei volumi finiti è un metodo utile nell'integrazione di equazioni differenziali alle derivate parziali. Tali equazioni vanno integrate in un volume sui cui confini sono imposte le condizioni al contorno.

L'interno di tale dominio viene quindi suddiviso in tanti volumi elementari, quindi tramite la forma integrale delle equazioni del problema considerato vengono scritte le relazioni che intercorrono tra i vari volumetti confinanti così da poter essere risolte per via numerica con l'ausilio del calcolatore. L'approssimazione risiede nel fatto che tali volumetti hanno dimensione finita e non infinitesima.

Esempio monodimensionale (1D)

modifica

Consideriamo il problema definito dalla seguente equazione alle derivate parziali:

 

in cui   rappresenta la variabile di stato e   rappresenta il flusso di  . In particolare, assumiamo   positivo o negativo a seconda della direzione del flusso. Se consideriamo l'equazione (1) relativa al flusso di materia attraverso una superficie di area costante, possiamo dividere l'intero dominio spaziale   in un certo numero di volumi finiti o celle, individuando con l'indice   il centro di ogni cella. Per una particolare cella   possiamo definire il volume medio di   al tempo   e  , come:

 

e il volume medio relativo al tempo  , come:

 

in cui   e   individuano le posizioni delle facce del flusso uscente ed entrante, relative alla   cella.

Integrando l'equazione (1) rispetto al tempo, otteniamo:

 

Per ottenere il volume medio di   al tempo  , integriamo   su tutto il volume della cella   e dividiamo il risultato per  , quindi

 

Assumiamo una certa regolarità di  , e che possiamo invertire l'ordine di integrazione. Essendo il flusso normale alla superficie di area unitaria della cella, siccome in una dimensione  , possiamo applicare il teorema della divergenza, sostituendo l'integrale di volume della divergenza con il valore di   assunto nelle facce   e   del volume finito, ovvero:

 

in cui   e  .

Possiamo inoltre derivare uno schema numerico semi-discreto per il problema seguente con il centro delle celle indicizzato con  , e utilizzando come indici per i flussi sulle facce  ; differenziando la (6) rispetto al tempo otteniamo:

 

dove i valori   dei flussi sulle facce possono essere ricavati da un'operazione di interpolazione o estrapolazione delle medie relative ad ogni cella. È da notare che l'equazione (7) è esatta per quanto riguarda i volumi medi, nel senso che a tale proposito non è stata introdotta, durante la trattazione svolta, nessuna approssimazione.

Voci correlate

modifica
Controllo di autoritàThesaurus BNCF 69870 · LCCN (ENsh95001595 · GND (DE4220855-5 · BNF (FRcb150249492 (data) · J9U (ENHE987007537003305171
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica