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.

An efficient semismooth Newton-AMG-based inexact primal-dual algorithm for generalized transport problems

A Fast Solver for Generalized Optimal Transport Problems Based on Dynamical System and Algebraic Multigrid

An Efficient Interior-Point Algorithm with New Non-Monotone Line Search Filter Method for Nonlinear Constrained Programming

A Globalized Inexact Semismooth Newton Method for Nonsmooth Fixed-point Equations involving Variational Inequalities

A Multiscale Semi-Smooth Newton Method for Optimal Transport

An Augmented Lagrangian Primal-Dual Semismooth Newton Method for Multi-Block Composite Optimization

A Sparse Smoothing Newton Method for Solving Discrete Optimal Transport Problems

A Stochastic Semismooth Newton Method for Nonsmooth Nonconvex Optimization.

Numerical method for feasible and approximately optimal solutions of multi-marginal optimal transport beyond discrete measures

A Newton-CG Augmented Lagrangian Method for Semidefinite Programming

A SemiSmooth Newton Method for Semidefinite Programs and Its Applications in Electronic Structure Calculations.

A Primal Majorized Semismooth Newton-CG Augmented Lagrangian Method for Large-Scale Linearly Constrained Convex Programming

An inexact regularized proximal Newton method for nonconvex and nonsmooth optimization

Structure Preserving Primal Dual Methods for Gradient Flows with Nonlinear Mobility Transport Distances

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

Nonlinearly preconditioned semismooth Newton algorithms for nonlinear nonsmooth systems

A nonsmooth primal-dual method with interwoven PDE constraint solver

A Semismooth Newton Based Augmented Lagrangian Method for Nonsmooth Optimization on Matrix Manifolds

Implementable Semismooth* Newton Methods for Generalized Equations are G-Semismooth Newton Methods

A projected semismooth Newton method for a class of nonconvex composite programs with strong prox-regularity

On the local convergence of the semismooth Newton method for composite optimization