Abstract:Large-scale modeling and predictive simulations of the subsurface flow and reactive transport system in porous media is significantly challenging, due to the high nonlinearity of the governing equations and the strong heterogeneity of material coefficients. The design of novel mathematical models and state-of-the-art methods for the flow simulation through porous media typically needs to satisfy the so-called bound-preserving property, i.e., the computed solution should stay inside a physically meaningful range. This paper presents a robust, scalable numerical framework based on the variational inequality formulation and the semismooth Newton method in a fully implicit manner, to model and simulate highly nonlinear flows without violating the boundedness requirement of the solution. Rigorous theoretical analysis for the variational inequality formulation of the problem is provided for facilitating the design of algorithms. Specifically, our approach further enhances the numerical formulation by utilizing a family of multilevel monolithic overlapping Schwarz methods for efficiently preconditioning, and the parallel implementation of the bound-preserving solver is based on the fast and robust domain decomposition technique. Numerical experiments are presented to demonstrate the efficiency and parallel scalability of the solution strategy for both standard benchmarks as well as realistic flow problems involving strong heterogeneity and high nonlinearity. We also show that the proposed framework is more robust and efficient than the commonly used inexact Newton algorithm in terms of the bound-preserving property.

Numerical study of a domain decomposition method for a two-scale linear transport equation

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

Numerical solution of 2-D incompressible viscous flow with high order compact scheme using domain decomposition and matched method

Numerical studies of a class of reaction–diffusion equations with Stefan conditions

Numerical methods for multiscale transport equations and application to two-phase porous media flow

Space-time domain decomposition for advection-diffusion problems in mixed formulations

Numerical Optimal Transport from 1D to 2D using a Non-local Monge-Ampère Equation

Asymptotic analysis of domain decomposition for optimal transport

Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations

On the implementation of linear finite element method for nonlocal diffusion model over 2D domain

The Neumann numerical boundary condition for transport equations

High-Order Numerical Method for 1D Non-local Diffusive Equation

Domain and Subband decomposition approach for 2 D simulation of quantum transport phenomena ∗

Lattice-boltzmann simulation of two-dimensional super-diffusion

Scalable Semismooth Newton Methods with Multilevel Domain Decomposition for Subsurface Flow and Reactive Transport in Porous Media

A conservative Eulerian finite element method for transport and diffusion in moving domains

Numerical Solution of the Nonlocal Diffusion Equation on the Real Line.

On Convergence of a Domain Decomposition Method for a Scalar Conservation Law

A high order conservative semi-Lagrangian discontinuous Galerkin method for two-dimensional transport simulations

A New Continuum Approach for Nonlinear Kinetic Simulation and Transport Analysis

A numerical study on solving a fractional time-space diffusion equation via the finite difference method