We compare the accuracy of five numerical schemes in modeling transport of nonreactive and reactive solutes in porous formations with heterogeneity increasing from low to very high. Two schemes, the Total Variation Diminishing (TVD) and the Eulerian-Lagrangian Method of Characteristics (MOC), are available in widely used packages. The other three schemes are the Random Walk Particle Tracking (RWPT), the Smoothed Particle Hydrodynamics (SPH) and a Streamline-Based (SB-FV) method, which we modified to improve its accuracy. The advective nature of the transport problem renders the numerical solution very challenging with the solutions provided by classic Eulerian methods that are plagued by numerical diffusion and spurious oscillations. Our analysis shows that TVD is severely affected by numerical diffusion, while the modified SB-FV method shows the tendency to underestimate dilution to an extent that increases with . In addition, we show that MOC is not mass-conservative, SPH is computationally demanding and cannot handle anisotropic dispersion, while RWPT develops spurious concentration fluctuations, which can be attenuated by increasing the number of particles at the expenses of an increase of the CPU time. Moreover, we investigate the effect of uniform and non-uniform local dispersion models on the overall plume dilution. These results help to consciously choose the numerical scheme according to investigation’s objectives and heterogeneity degree.

Numerical simulations of solute transport in highly heterogeneous formations: A comparison of alternative numerical schemes

Boso, Francesca;Bellin, Alberto;Dumbser, Michael
2013-01-01

Abstract

We compare the accuracy of five numerical schemes in modeling transport of nonreactive and reactive solutes in porous formations with heterogeneity increasing from low to very high. Two schemes, the Total Variation Diminishing (TVD) and the Eulerian-Lagrangian Method of Characteristics (MOC), are available in widely used packages. The other three schemes are the Random Walk Particle Tracking (RWPT), the Smoothed Particle Hydrodynamics (SPH) and a Streamline-Based (SB-FV) method, which we modified to improve its accuracy. The advective nature of the transport problem renders the numerical solution very challenging with the solutions provided by classic Eulerian methods that are plagued by numerical diffusion and spurious oscillations. Our analysis shows that TVD is severely affected by numerical diffusion, while the modified SB-FV method shows the tendency to underestimate dilution to an extent that increases with . In addition, we show that MOC is not mass-conservative, SPH is computationally demanding and cannot handle anisotropic dispersion, while RWPT develops spurious concentration fluctuations, which can be attenuated by increasing the number of particles at the expenses of an increase of the CPU time. Moreover, we investigate the effect of uniform and non-uniform local dispersion models on the overall plume dilution. These results help to consciously choose the numerical scheme according to investigation’s objectives and heterogeneity degree.
2013
Boso, Francesca; Bellin, Alberto; Dumbser, Michael
File in questo prodotto:
File Dimensione Formato  
13062717502616160.pdf

Solo gestori archivio

Tipologia: Versione editoriale (Publisher’s layout)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 616.17 kB
Formato Adobe PDF
616.17 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11572/93787
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 61
  • ???jsp.display-item.citation.isi??? 59
social impact