Two Nitsche-based mixed finite element discretizations for the seepage problem in Richards' equation

Federico Gatti,Andrea Bressan,Alessio Fumagalli,Domenico Gallipoli,Leonardo Maria Lalicata,Simone Pittaluga,Lorenzo Tamellini

2024-07-11

Abstract:This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in practice require switching between Neumann and Dirichlet boundary conditions on unknown portions on the boundary. Upon realizing the similarities of these conditions with unilateral contact problems in mechanics, we take inspiration from that literature to propose two approaches: the first method relies on a strongly consistent penalization term, whereas the second one is obtained by an hybridization approach, in which the value of the pressure on the surface is treated as a separate set of unknowns. The flow problem is discretized in mixed form with div-conforming elements so that the water mass is preserved. Numerical experiments show the validity of the proposed strategy in handling the seepage boundary conditions on geometries with increasing complexity.

Numerical Analysis

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the difficult boundary - condition problems encountered when dealing with groundwater seepage in the Richards equation, especially how to accurately determine and handle seepage faces in unsaturated soil. Specifically:

1. **Definition of seepage face**:
- A seepage face refers to the place where the saturated area in groundwater intersects with the ground surface, and groundwater flows out from these places and becomes surface water. For example, in the case of a slope, groundwater may seep out from the slope surface.

2. **Importance of the problem**:
- Correctly determining the seepage face is crucial for predicting key parameters (such as pore - water pressure, effective stress, and soil shear strength), which directly affects the safety assessment of slope stability, embankments, and other geotechnical engineering structures. Incorrect determination of the seepage face will lead to misjudgment of potential disasters (such as landslides and slope instability), thus posing a significant risk to infrastructure and life safety.

3. **Limitations of existing methods**:
- Existing finite - element algorithms usually determine the seepage face by switching between Neumann boundary conditions (flux - type) and Dirichlet boundary conditions (head - type) through iterative methods. However, these methods have problems such as imperfect mathematical theory, convergence and stability issues, and water - conservation problems.

4. **Proposed new methods**:
- The paper proposes two new mixed finite - element discretization strategies based on the Nitsche method to more accurately handle the boundary conditions of seepage faces. These two methods use strong - consistent penalty terms and hybridization methods respectively, and the latter treats the head values on the boundary as independent unknowns. These methods draw on the Karush - Kuhn - Tucker conditions in contact mechanics problems and their effectiveness on complex geometries has been verified through numerical experiments.

5. **Objective**:
- To provide a more robust and accurate method for simulating seepage faces in groundwater flow, thereby providing a more reliable basis for early warning and prevention of natural disasters such as landslide triggering.

Through these improvements, the paper hopes to more accurately capture the changes of seepage faces in groundwater flow simulation, and further improve the ability to predict and evaluate related geological disasters.