next up previous
Next: Il codice di Up: Il codice di Previous: Soluzione del sistema

Calcolo dei bilanci di massa

La correttezza della soluzione molto spesso controllata verificando che il principio di conservazione della massa venga rispettato. Per questo si deve calcolare la massa di acqua che entra nel dominio e verificare che questa sia uguale alla massa d'acqua che esce dal dominio stesso.

La massa d'acqua entrata o uscita dai nodi di Dirichlet viene calcolata semplicemente moltiplicando per il vettore soluzione la matrice di rigidezza , come calcolata prima dell'imposizione delle condizioni al contorno.

Per fare ciø bisogna salvare i coefficienti diagonali della matrice prima che essi siano cambiati dall'imposizione delle condizioni al contorno.

La massa d'acqua entrata o uscita dai nodi di Neumann data direttamente dal valore della portata in quei nodi.

Il bilancio di massa viene eseguito sommando algebricamente i valori ottenuti per ciascun nodo di Dirichlet e di Neumann. Il risultato deve essere teoricamente nullo. In pratica il bilancio di massa si considera corretto se tale risultato è uguale ad una frazione piccola della massa entrante o uscente.



Andrea Pellizzon
Thu Jun 5 11:00:04 MET DST 1997