A new framework of high-order unfitted finite element methods using ALE maps for moving-domain problems

Wenhao Lu,Chuwen Ma,Weiying Zheng

2024-04-24

Abstract:As a sequel to our previous work [C. Ma, Q. Zhang and W. Zheng, SIAM J. Numer. Anal., 60 (2022)], [C. Ma and W. Zheng, J. Comput. Phys. 469 (2022)], this paper presents a generic framework of arbitrary Lagrangian-Eulerian unfitted finite element (ALE-UFE) methods for partial differential equations (PDEs) on time-varying domains. The ALE-UFE method has a great potential in developing high-order unfitted finite element methods. The usefulness of the method is demonstrated by a variety of moving-domain problems, including a linear problem with explicit velocity of the boundary (or interface), a PDE-domain coupled problem, and a problem whose domain has a topological change. Numerical experiments show that optimal convergence is achieved by both third- and fourth-order methods on domains with smooth boundaries, but is deteriorated to the second order when the domain has topological changes.

Numerical Analysis

What problem does this paper attempt to address?

The paper attempts to address the problem of solving partial differential equations (PDEs) on time-varying domains using high-order non-conforming finite element methods. Specifically:

1. **Proposing a new framework**: The paper introduces a high-order non-conforming finite element method (ALE-UFE) based on Arbitrary Lagrangian-Eulerian (ALE) mapping for solving partial differential equations on time-varying domains.

2. **Handling moving domain problems**: For moving domain problems, traditional numerical methods often require constant mesh updates, which can lead to additional computational overhead. The method proposed in this paper aims to overcome this challenge, especially when dealing with large deformations or topological changes.

3. **Validating the method's effectiveness**: The effectiveness of the method is validated through various moving domain problems, including linear problems with explicit boundary velocities, PDE-domain coupling problems, and problems with topological changes.

4. **Convergence analysis**: Numerical experiments show that on domains with smooth boundaries, the method can achieve optimal convergence rates of 3rd and 4th order; however, in the presence of topological changes, the convergence rate degrades to 2nd order.

In summary, the main objective of the paper is to develop an efficient and robust numerical method suitable for time-varying domain problems and to demonstrate its effectiveness and superiority in various complex scenarios.