Modeling subsurface water resource systems involving heterogeneous porous media using the variational multiscale formulation

• Published in: Journal of hydrology (Amsterdam) Vol. 428-429; pp. 1 - 14
• Main Authors: Turner, D.Z., Nakshatrala, K.B., Martinez, M.J., Notz, P.K.
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier B.V 27.03.2012; Elsevier
• Subjects: Earth sciences; Earth, ocean, space; engineering; Exact sciences and technology; Flow through porous media; Hydrogeology; Hydrology; Hydrology. Hydrogeology; Local mass balance; Mixed finite element formulations; permeability; porous media; Raviart–Thomas spaces; Spatial heterogeneity; subsurface flow; Variational multiscale formulation; viscosity; water resources; Spatial heterogeneity; Flow through porous media; Local mass balance; Variational multiscale formulation; Raviart–Thomas spaces; Mixed finite element formulations; experimental studies; mass balance; permeability; ground water; subsurface flow; Raviart-Thomas spaces; water resources; boundary conditions; advection; models; hydrodynamics; transport; velocity; finite element analysis; materials properties; subsurface; heterogeneity; interpolation; viscosity; heterogeneous materials; flow field; porous media
• ISSN: 0022-1694; 1879-2707
• DOI: 10.1016/j.jhydrol.2011.12.024

[Display omitted] ► Variational multiscale methods perform poorly for aquifer injection problems. ► Highly pressure dependent viscosity models show reduced local mass balance error. ► We present estimates of the error produced by variational multiscale methods. ► These results are useful for judging accuracy of models that inform policy decisions. This work compares two popular mixed finite element formulations used to model subsurface flow and transport in heterogeneous porous media, namely, the lowest order Raviart–Thomas and the variational multiscale stabilized formulations. Comparison is made based on performance for several problems of engineering relevance that involve highly heterogeneous material properties (permeability ratios of up to 1×105), open flow boundary conditions (pressure driven flows), and large scale domains in two dimensions. Numerical experiments are performed to show the degree to which mass conservation is violated when a flow field computed using either element is used as the advection velocity in a transport model. The results reveal that the equal-order interpolation under the variational multiscale formulation shows considerable mass production or loss for problems that involve flow tangential to layers of differing permeability. But the violation of local mass balance is marginal for problems in which flow is orthogonal to the layers of differing permeability. For problems involving pressure dependent viscosity, we show that models with a high degree of pressure dependence exhibit improved performance for the variational multiscale method. The results are useful in establishing rudimentary estimates of the error produced by using the variational mutliscale formulation for several different types of problems related to subsurface water resource systems.