Adaptive Mesh Refinement for Variational Inequalities

TL;DR

This paper introduces two adaptive mesh refinement techniques, VCES and UDO, for variational inequalities, improving free boundary resolution and solver efficiency in finite element methods.

Contribution

It presents novel AMR strategies specifically designed for variational inequalities, with theoretical analysis and implementation in Firedrake.

Findings

Both methods effectively refine meshes around free boundaries.

Numerical results show improved convergence rates and computational efficiency.

The techniques are compatible with existing finite element solvers.

Abstract

Variational inequalities play a pivotal role in a wide array of scientific and engineering applications. This project presents two techniques for adaptive mesh refinement (AMR) in the context of variational inequalities, with a specific focus on the classical obstacle problem. We propose two distinct AMR strategies: Variable Coefficient Elliptic Smoothing (VCES) and Unstructured Dilation Operator (UDO). VCES uses a nodal active set indicator function as the initial iterate to a time-dependent heat equation problem. Solving a single step of this problem has the effect of smoothing the indicator about the free boundary. We threshold this smoothed indicator function to identify elements near the free boundary. Key parameters such as timestep and threshold values significantly influence the efficacy of this method. The second strategy, UDO, focuses on the discrete identification of elements adjacent to the free boundary, employing a graph-based approach to mark neighboring elements for refinement. This technique resembles the dilation morphological operation in image processing, but tailored for unstructured meshes. We also examine the theory of variational inequalities, the convergence behavior of finite element solutions, and implementation in the Firedrake finite element library. Convergence analysis reveals that accurate free boundary estimation is pivotal for solver performance. Numerical experiments demonstrate the effectiveness of the proposed methods in dynamically enhancing mesh resolution around free boundaries, thereby improving the convergence rates and computational efficiency of variational inequality solvers. Our approach integrates seamlessly with existing Firedrake numerical solvers, and it is promising for solving more complex free boundary problems.