Non-negative mixed finite element formulations for a tensorial diffusion equation

Abstract

We consider the tensorial diffusion equation, and address the discrete maximum-minimum principle of mixed finite element formulations. In particular, we address non-negative solutions (which is a special case of the maximum-minimum principle) of mixed finite element formulations. The discrete maximum-minimum principle is the discrete version of the maximum-minimum principle. In this paper we present two non-negative mixed finite element formulations for tensorial diffusion equations based on constrained optimization techniques (in particular, quadratic programming). These proposed mixed formulations produce non-negative numerical solutions on arbitrary meshes for low-order (i.e., linear, bilinear and trilinear) finite elements. The first formulation is based on the Raviart-Thomas spaces, and is obtained by adding a non-negative constraint to the variational statement of the Raviart-Thomas formulation. The second non-negative formulation based on the variational multiscale formulation. For the former formulation we comment on the affect of adding the non-negative constraint on the local mass balance property of the Raviart-Thomas formulation. We also study the performance of the active set strategy for solving the resulting constrained optimization problems. The overall performance of the proposed formulation is illustrated on three canonical test problems.

arXiv comment40 pages using amsart style file, and 15 figures

References in corpus (2)

• Well-Centered Triangulation
• Numerical method for Darcy flow derived using Discrete Exterior Calculus

Cited by in corpus (8)

• A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and non-negative constraint
• Variational inequality approach to enforce the non-negative constraint for advection-diffusion equations
• Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids
• A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
• Do current lattice Boltzmann methods for diffusion and diffusion-type equations respect maximum principles and the non-negative constraint?
• Large-scale Optimization-based Non-negative Computational Framework for Diffusion Equations: Parallel Implementation and Performance Studies
• On mesh restrictions to satisfy comparison principles, maximum principles, and the non-negative constraint: Recent developments and new results
• On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method