skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Hybrid Multiscale Finite Volume Method for Advection-Diffusion Equations Subject to Heterogeneous Reactive Boundary Conditions

Journal Article · · Multiscale Modeling & Simulation
DOI:https://doi.org/10.1137/15M1022537· OSTI ID:1344647

We present a hybrid scheme for the coupling of macro and microscale continuum models for reactive contaminant transport in fractured and porous media. The transport model considered is the advection-dispersion equation, subject to linear heterogeneous reactive boundary conditions. The Multiscale Finite Volume method (MsFV) is employed to define an approximation to the microscale concentration field defined in terms of macroscopic or \emph{global} degrees of freedom, together with local interpolator and corrector functions capturing microscopic spatial variability. The macroscopic mass balance relations for the MsFV global degrees of freedom are coupled with the macroscopic model, resulting in a global problem for the simultaneous time-stepping of all macroscopic degrees of freedom throughout the domain. In order to perform the hybrid coupling, the micro and macroscale models are applied over overlapping subdomains of the simulation domain, with the overlap denoted as the handshake subdomain $$\Omega^{hs}$$, over which continuity of concentration and transport fluxes between models is enforced. Continuity of concentration is enforced by posing a restriction relation between models over $$\Omega^{hs}$$. Continuity of fluxes is enforced by prolongating the macroscopic model fluxes across the boundary of $$\Omega^{hs}$$ to microscopic resolution. The microscopic interpolator and corrector functions are solutions to local microscopic advection-diffusion problems decoupled from the global degrees of freedom and from each other by virtue of the MsFV decoupling ansatz. The error introduced by the decoupling ansatz is reduced iteratively by the preconditioned GMRES algorithm, with the hybrid MsFV operator serving as the preconditioner.

Research Organization:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Organization:
USDOE
DOE Contract Number:
AC05-76RL01830
OSTI ID:
1344647
Report Number(s):
PNNL-SA-110744; KJ0401000
Journal Information:
Multiscale Modeling & Simulation, Vol. 14, Issue 4; ISSN 1540-3459
Publisher:
SIAM
Country of Publication:
United States
Language:
English

Similar Records

A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows
Journal Article · Fri Jul 12 00:00:00 EDT 2019 · Journal of Computational Physics · OSTI ID:1344647

Algorithms for propagating uncertainty across heterogeneous domains
Journal Article · Wed Dec 30 00:00:00 EST 2015 · SIAM Journal on Scientific Computing · OSTI ID:1344647

A hybrid pore‐scale and continuum‐scale model for solute diffusion, reaction, and biofilm development in porous media
Journal Article · Tue Mar 31 00:00:00 EDT 2015 · Water Resources Research · OSTI ID:1344647

Related Subjects