Home GeoForum Forums Geologia & Professione Idrogeologia & Idraulica FEMWATER

Originariamente inviato da s.zanarello:
Ciao Blizzard
per prima cosa una domanda: come fai a risolvere con modflow l'equazione di flusso tenendo conto della densità?

Comunque, secondo me anche risolvendo prima il flusso tenendo conto della densità e poi risolvendo il trasporto, tale metodo non è valido per simulare ad es. un'intrusione salina, in quanto si tratta di una soluzione disaccoppiata.
Per soluzione disaccoppiata intendo: in una prima fase viene determinato il campo di flusso (ad es con modflow), poi, sulla base del campo di flusso determinato, si risolvono le equazioni di trasporto. Tale approccio è disaccoppiato, in quanto presuppone che i fenomeni di trasporto (e quindi la densità del fluido), non influenzino il campo di moto, che viene risolto prima.
Un approccio accoppiato, invece, non vede questi due passi successivi (prima flusso, poi trasporto), ma ad ogni step temporale risolve contemporaneamente sia flusso che trasporto, in modo che anche il campo di flusso risulti influenzato dai fenomeni di trasporto. In questo caso, infatti, i fenomeni di deflusso e trasporto si influenzano a vicenda, e non possono quindi essere risolti separatamente. Ovviamente tale approccio numerico è decisamente più pesante dal punto di vista computazionale. Per maggior chiarezza ti invio per mail (indirizzo del tuo profilo) due pagine che spiegano sinteticamente, ma chiaramente la differenza tra approccio disaccoppiato ed accoppiato.

So per certo che modflow+mt3d non possono essere usati per un approccio accoppiato; l'usgs, infatti, ha creato il codice seawat per tali problemi. Seawat si basa proprio su modflow e mt3dms, ma i due codici sono integrati, per risolvere flusso e trasporto contemporaneamente (tenendo conto delle mutue influenze), e non prima uno e poi l'altro.
Spero di esserti stato utile.