Variational inequality approach to enforcing the non-negative constraint for advection-diffusion equations
From MaRDI portal
Publication:2309846
DOI10.1016/J.CMA.2017.03.022zbMATH Open1439.35153arXiv1611.08758OpenAlexW2559474616MaRDI QIDQ2309846FDOQ2309846
Authors: J. Chang, K. B. Nakshatrala 
Publication date: 6 April 2020
Published in: Computer Methods in Applied Mechanics and Engineering (Search for Journal in Brave)
Abstract: Predictive simulations are crucial for the success of many subsurface applications, and it is highly desirable to obtain accurate non-negative solutions for transport equations in these numerical simulations. In this paper, we propose a computational framework based on the variational inequality (VI) which can also be used to enforce important mathematical properties (e.g., maximum principles) and physical constraints (e.g., the non-negative constraint). We demonstrate that this framework is not only applicable to diffusion equations but also to non-symmetric advection-diffusion equations. An attractive feature of the proposed framework is that it works with with any weak formulation for the advection-diffusion equations, including single-field formulations, which are computationally attractive. A particular emphasis is placed on the parallel and algorithmic performance of the VI approach across large-scale and heterogeneous problems. It is also shown that QP and VI are equivalent under certain conditions. State-of-the-art QP and VI solvers available from the PETSc library are used on a variety of steady-state 2D and 3D benchmarks, and a comparative study on the scalability between the QP and VI solvers is presented. We then extend the proposed framework to transient problems by simulating the miscible displacement of fluids in a heterogeneous porous medium and illustrate the importance of enforcing maximum principles for these types of coupled problems. Our numerical experiments indicate that VIs are indeed a viable approach for enforcing the maximum principles and the non-negative constraint in a large-scale computing environment. Also provided are Firedrake project files as well as a discussion on the computer implementation to help facilitate readers in understanding the proposed framework.
Full work available at URL: https://arxiv.org/abs/1611.08758
Recommendations
- Large-scale optimization-based non-negative computational framework for diffusion equations: parallel implementation and performance studies
- On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method
- A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations
- A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
- Enforcement of constraints and maximum principles in the variational multiscale method
quadratic programmingparallel computingvariational inequalitiesadvection-diffusion equationsmiscible displacementnon-negative solutions
Boundary value problems for second-order elliptic equations (35J25) Unilateral problems for linear elliptic equations and variational inequalities with linear elliptic operators (35J86) Numerical methods for variational inequalities and related problems (65K15) Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.) (74F10) Flows in porous media; filtration; seepage (76S05) Variational inequalities (global problems) in infinite-dimensional spaces (58E35)
Cites Work
- Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities
- Title not available (Why is that?)
- Newton's Method for Large Bound-Constrained Optimization Problems
- Automated solution of differential equations by the finite element method. The FEniCS book
- Firedrake, automating the finite element method by composing abstractions
- Unified form language: a domain-specific language for weak formulations of partial differential equations
- Title not available (Why is that?)
- Title not available (Why is that?)
- Elliptic partial differential equations of second order
- Title not available (Why is that?)
- Finite-Dimensional Variational Inequalities and Complementarity Problems
- Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations
- Solution of variational inequalities in mechanics
- A nonsmooth version of Newton's method
- Title not available (Why is that?)
- Numerical solution of saddle point problems
- Unified Hybridization of Discontinuous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems
- Title not available (Why is that?)
- Semismooth and Semiconvex Functions in Constrained Optimization
- An Introduction to Variational Inequalities and Their Applications
- Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
- Convergence Analysis of Some Algorithms for Solving Nonsmooth Equations
- Title not available (Why is that?)
- Variational inequalities
- A semismooth equation approach to the solution of nonlinear complementarity problems
- Minimization of functions having Lipschitz continuous first partial derivatives
- Chaotic streamlines in the ABC flows
- Title not available (Why is that?)
- Title not available (Why is that?)
- Title not available (Why is that?)
- Discontinuous Galerkin methods for flow and transport problems in porous media
- Applied functional analysis. Functional analysis, Sobolev spaces and elliptic differential equations
- Maximum principle and uniform convergence for the finite element method
- A special newton-type optimization method
- Title not available (Why is that?)
- A Note on Preconditioning for Indefinite Linear Systems
- Elliptic Equations: An Introductory Course
- The mimetic finite difference method for elliptic problems
- Title not available (Why is that?)
- Finance and variational inequalities
- Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes
- Linear elliptic differential systems and eigenvalue problems
- Finite volume monotone scheme for highly anisotropic diffusion operators on unstructured triangular meshes. (Schéma volumes finis monotone pour des opérateurs de diffusion fortement anisotropes sur des maillages de triangles non structurés).
- Miscible porous media displacements in the quarter five-spot configuration. Part 1. The homogeneous case
- Discontinuous Galerkin methods
- Title not available (Why is that?)
- Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids
- A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and the non-negative constraint
- On the performance of high‐order finite elements with respect to maximum principles and the nonnegative constraint for diffusion‐type equations
- A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations
- Non-negative mixed finite element formulations for a tensorial diffusion equation
- An explicit expression for the penalty parameter of the interior penalty method
- Numerical treatment of partial differential equations. Revised translation of the 3rd German edition of `Numerische Behandlung partieller Differentialgleichungen' by Martin Stynes.
- Chaotic and Resonant Streamlines in the ABC Flow
- Monotonicity recovering and accuracy preserving optimization methods for postprocessing finite element solutions
- Stability of semidiscrete formulations for parabolic problems at small time steps
- The semismooth algorithm for large scale complementarity problems
- On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method
- Flexible complementarity solvers for large-scale applications
- Large-scale optimization-based non-negative computational framework for diffusion equations: parallel implementation and performance studies
- Efficient Assembly of $H(\mathrm{div})$ and $H(\mathrm{curl})$ Conforming Finite Elements
- High order discontinuous Galerkin method for simulating miscible flooding in porous media
- Numerical solutions of the incompressible miscible displacement equations in heterogeneous media
- Discrete maximum principle for the weak Galerkin method for anisotropic diffusion problems
- A new nonlinear finite volume scheme preserving positivity for diffusion equations
- An algorithm for the optimization of finite element integration loops
- Automated generation and symbolic manipulation of tensor product finite elements
- A parallel edge orientation algorithm for quadrilateral meshes
Cited In (26)
- A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations
- Parallel generalized Lagrange-Newton method for fully coupled solution of PDE-constrained optimization problems with bound-constraints
- A minimum-type nonlinear complementarity simulator with constrained pressure residual (CPR) methods for wormhole propagation in carbonate acidization
- Scalable semismooth Newton methods with multilevel domain decomposition for subsurface flow and reactive transport in porous media
- High-order bounds-satisfying approximation of partial differential equations via finite element variational inequalities
- Parallel reservoir simulators for fully implicit complementarity formulation of multicomponent compressible flows
- A segregated approach for modeling the electrochemistry in the 3-D microstructure of li-ion batteries and its acceleration using block preconditioners
- Adaptive fully implicit simulator with multilevel Schwarz methods for gas reservoir flows in fractured porous media
- Parallel fully coupled methods for bound-preserving solution of subsurface flow and transport in porous media
- A full approximation scheme multilevel method for nonlinear variational inequalities
- Composable block solvers for the four-field double porosity/permeability model
- A Monotone Finite Volume Scheme with Fixed Stencils for 3D Heat Conduction Equation
- Application of regularization technique to variational adjoint method: A case for nonlinear convection-diffusion problem
- Variational inequality transport model on the sphere by the active-set reduced-space algorithm
- On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method
- A fully implicit constraint-preserving simulator for the black oil model of petroleum reservoirs
- A positivity-preserving finite volume scheme for heat conduction equation on generalized polyhedral meshes
- Lower bounds on the mix norm of passive scalars advected by incompressible enstrophy-constrained flows
- Comparative Study of Finite Element Methods Using the Time-Accuracy-Size(TAS) Spectrum Analysis
- A positivity-preserving finite volume scheme with least square interpolation for 3D anisotropic diffusion equation
- A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and the non-negative constraint
- A positivity-preserving pyramid scheme for anisotropic diffusion problems on general hexahedral meshes with nonplanar cell faces
- Nonlinearly preconditioned constraint-preserving algorithms for subsurface three-phase flow with capillarity
- A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
- A semi-Lagrangian mixed finite element method for advection-diffusion variational inequalities
- Large-scale optimization-based non-negative computational framework for diffusion equations: parallel implementation and performance studies
Uses Software
This page was built for publication: Variational inequality approach to enforcing the non-negative constraint for advection-diffusion equations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2309846)
