We present a model-order reduction technique that overcomes the computational burden associated with the application of Monte Carlo methods to the solution of the groundwater flow equation with random hydraulic conductivity. The method is based on the Galerkin projection of the high-dimensional model equations onto a subspace, approximated by a small number of pseudo-optimally chosen basis functions (principal components). To obtain an efficient reduced-order model, we develop an offline algorithm for the computation of the parameter-independent principal components. Our algorithm combines a greedy algorithm for the snapshot selection in the parameter space and an optimal distribution of the snapshots in time. Moreover, we introduce a residual-based estimation of the error associated with the reduced model. This estimation allows a considerable reduction of the number of full system model solutions required for the computation of principal components. We demonstrate the robustness of our methodology by way of numerical examples, comparing the empirical statistics of the ensemble of the numerical solutions obtained using the traditional Monte Carlo method and our reduced model. The numerical results show that our methodology significantly reduces the computational requirements (CPU time and storage) for the solution of the Monte Carlo simulation, ensuring a good approximation of the mean and variance of the head. The analysis of the empirical probability density functions at the observation wells suggests that our reduced model produces good results and is most accurate in the regions with large drawdown. ©2013. American Geophysical Union. All Rights Reserved.

A reduced-order model for groundwater flow equation with random hydraulic conductivity: Application to Monte Carlo methods

Pasetto D.;
2013-01-01

Abstract

We present a model-order reduction technique that overcomes the computational burden associated with the application of Monte Carlo methods to the solution of the groundwater flow equation with random hydraulic conductivity. The method is based on the Galerkin projection of the high-dimensional model equations onto a subspace, approximated by a small number of pseudo-optimally chosen basis functions (principal components). To obtain an efficient reduced-order model, we develop an offline algorithm for the computation of the parameter-independent principal components. Our algorithm combines a greedy algorithm for the snapshot selection in the parameter space and an optimal distribution of the snapshots in time. Moreover, we introduce a residual-based estimation of the error associated with the reduced model. This estimation allows a considerable reduction of the number of full system model solutions required for the computation of principal components. We demonstrate the robustness of our methodology by way of numerical examples, comparing the empirical statistics of the ensemble of the numerical solutions obtained using the traditional Monte Carlo method and our reduced model. The numerical results show that our methodology significantly reduces the computational requirements (CPU time and storage) for the solution of the Monte Carlo simulation, ensuring a good approximation of the mean and variance of the head. The analysis of the empirical probability density functions at the observation wells suggests that our reduced model produces good results and is most accurate in the regions with large drawdown. ©2013. American Geophysical Union. All Rights Reserved.
2013
49
File in questo prodotto:
File Dimensione Formato  
2013_Pasetto_et_al_wrr.pdf

accesso aperto

Descrizione: articolo principale
Tipologia: Versione dell'editore
Licenza: Creative commons
Dimensione 476.92 kB
Formato Adobe PDF
476.92 kB Adobe PDF Visualizza/Apri

I documenti in ARCA 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/10278/3722960
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 30
  • ???jsp.display-item.citation.isi??? 28
social impact