scholarly journals Bounding of Flow and Transport Analysis in Heterogeneous Saturated Porous Media: A Minimum Energy Dissipation Principle for the Bounding and Scale-Up

Hydrology ◽  
2019 ◽  
Vol 6 (2) ◽  
pp. 33 ◽  
Author(s):  
Nelson ◽  
Williams

We apply minimum kinetic energy principles from classic mechanics to heterogeneous porous media flow equations to derive and evaluate rotational flow components to determine bounding homogenous representations. Kelvin characterized irrotational motions in terms of energy dissipation and showed that minimum dynamic energy dissipation occurs if the motion is irrotational; i.e., a homogeneous flow system. For porous media flow, reductions in rotational flow represent heterogeneity reductions. At the limit, a homogeneous system, flow is irrotational. Using these principles, we can find a homogenous system that bounds a more complex heterogeneous system. We present mathematics for using the minimum energy principle to describe flow in heterogeneous porous media along with reduced special cases with the necessary bounding and associated scale-up equations. The first, simple derivation involves no boundary differences and gives results based on direct Kelvin-type minimum energy principles. It provides bounding criteria, but yields only a single ultimate scale-up. We present an extended derivation that considers differing boundaries, which may occur between scale-up elements. This approach enables a piecewise less heterogeneous representation to bound the more heterogeneous system. It provides scale-up flexibility for individual model elements with differing sizes, and shapes and supports a more accurate representation of material properties. We include a case study to illustrate bounding with a single direct scale-up. The case study demonstrates rigorous bounding and provides insight on using bounding flow to help understand heterogeneous systems. This work provides a theoretical basis for developing bounding models of flow systems. This provides a means to justify bounding conditions and results.

1996 ◽  
Vol 175 (1-4) ◽  
pp. 429-452 ◽  
Author(s):  
R. Mackay ◽  
T.A. Cooper ◽  
A.V. Metcalfe ◽  
P.E. O'Connell

2007 ◽  
Vol 9 (4) ◽  
pp. 251-266 ◽  
Author(s):  
David J. Hill ◽  
Barbara S. Minsker ◽  
Albert J. Valocchi ◽  
Vladan Babovic ◽  
Maarten Keijzer

Due to the considerable computational demands of modeling solute transport in heterogeneous porous media, there is a need for upscaled models that do not require explicit resolution of the small-scale heterogeneity. This study investigates the development of upscaled solute transport models using genetic programming (GP), a domain-independent modeling tool that searches the space of mathematical equations for one or more equations that describe a set of training data. An upscaling methodology is developed that facilitates both the GP search and the implementation of the resulting models. A case study is performed that demonstrates this methodology by developing vertically averaged equations of solute transport in perfectly stratified aquifers. The solute flux models developed for the case study were analyzed for parsimony and physical meaning, resulting in an upscaled model of the enhanced spreading of the solute plume, due to aquifer heterogeneity, as a process that changes from predominantly advective to Fickian. This case study not only demonstrates the use and efficacy of GP as a tool for developing upscaled solute transport models, but it also provides insight into how to approach more realistic multi-dimensional problems with this methodology.


Author(s):  
Sérgio Felipe Ferreira Silva ◽  
Hanna Thaina Prates Arimatéia ◽  
Alexandre Santos Francisco ◽  
Weslley Luiz da Silva Assis

Multiscale methods are usually developed for solving second-order elliptic problems in which coefficients are of multiscale heterogeneous nature. The Multiscale Mixed Method (MuMM) was introduced aiming at the efficient and accurate approximation of large flow problems in highly heterogeneous porous media. In the MuMM numerical solver, first mixed multiscale basis functions are constructed, and next global domain decomposition iterations are performed to compute the discrete solution of the problems. However, this iterative procedure is a time-consuming step. In this paper, the authors improve the MuMM solver through the implementation of parallel computations in the step concerning the global iterative procedure. The parallel version of the solver employs the application programming interface Open Multi-Processing (OpenMP). The implementation with the OpenMP reduces significantly the computational effort to perform the domain decomposition iterations, as indicated by the numerical results.


Sign in / Sign up

Export Citation Format

Share Document