next up previous
Next: Soluzione agli elementi Up: Corso di Metodi Numerici Previous: Corso di Metodi Numerici

Equazione da risolvere

Il flusso dell'acqua in un mezzo poroso retto dall'equazione di Darcy:

 

dove:

Nel caso pigenerale il mezzo poroso da considerarsi non omogeneo e anisotropo. Il coefficiente di filtrazione diventa perciø un tensore funzione della posizione:

Se si assume la profondit dell'acquifero costante e pari a B il moto risulta bidimensionale in un piano orizzontale. In tal caso si puø definire un coefficiente di trasmissivit caratteristico dell'acquifero: . Dal principio di conservazione della massa di acqua si ricava l'equazione generale che regge il moto in un acquifero in pressione:

 

Ad essa vanno aggiunte le condizioni al contorno di Dirichlet o di Neumann:

  

Nella presente esercitazione si vuole risolvere numericamente il problema conservazione_massa in un dominio circolare di centro e raggio unitario, con trasmissivit costante ed unitaria. Sull'intero dominio vi una ricarica uniforme e costante pari a f(x,y)=4.

Il problema da risolvere dunque il seguente:

 

con condizioni al contorno di Dirichlet:

 

La soluzione analitica di tale problema data da:

L'equazione esercitazione viene risolta utilizzando un metodo agli elementi finiti (FEM) con elementi triangolari e funzioni base lineari.





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