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)
modificaConsideriamo 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
modificaControllo di autorità | Thesaurus BNCF 69870 · LCCN (EN) sh95001595 · GND (DE) 4220855-5 · BNF (FR) cb150249492 (data) · J9U (EN, HE) 987007537003305171 |
---|