Abstract Finite element methods developed for unfitted meshes have been widely applied to various interface problems. However, many of them resort to non-conforming spaces for approximation, which is a critical obstacle for the extension to $\textbf{H}(\text{curl})$ equations. This essential issue stems from the underlying Sobolev space $\textbf{H}^s(\text{curl};\,\Omega)$ , and even the widely used penalty methodology may not yield the optimal convergence rate. One promising approach to circumvent this issue is to use a conforming test function space, which motivates us to develop a Petrov–Galerkin immersed finite element (PG-IFE) method for $\textbf{H}(\text{curl})$ -elliptic interface problems. We establish the Nédélec-type IFE spaces and develop some important properties including their edge degrees of freedom, an exact sequence relating to the $H^1$ IFE space and optimal approximation capabilities. We analyse the inf-sup condition under certain assumptions and show the optimal convergence rate, which is also validated by numerical experiments.
more »
« less
Immersed virtual element methods for electromagnetic interface problems in three dimensions
Finite element methods for electromagnetic problems modeled by Maxwell-type equations are highly sensitive to the conformity of approximation spaces, and non-conforming methods may cause loss of convergence. This fact leads to an essential obstacle for almost all the interface-unfitted mesh methods in the literature regarding the application to electromagnetic interface problems, as they are based on non-conforming spaces. In this work, a novel immersed virtual element method for solving a three-dimensional (3D) H(curl) interface problem is developed, and the motivation is to combine the conformity of virtual element spaces and robust approximation capabilities of immersed finite element spaces. The proposed method is able to achieve optimal convergence. To develop a systematic framework, the [Formula: see text], H(curl) and H(div) interface problems and their corresponding problem-orientated immersed virtual element spaces are considered all together. In addition, the de Rham complex will be established based on which the Hiptmair–Xu (HX) preconditioner can be used to develop a fast solver for the H(curl) interface problem.
more »
« less
- PAR ID:
- 10423062
- Date Published:
- Journal Name:
- Mathematical Models and Methods in Applied Sciences
- Volume:
- 33
- Issue:
- 03
- ISSN:
- 0218-2025
- Page Range / eLocation ID:
- 455 to 503
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We consider the mathematical analysis and numerical approximation of a system of nonlinear partial differential equations that arises in models that have relevance to steady isochoric flows of colloidal suspensions. The symmetric velocity gradient is assumed to be a monotone nonlinear function of the deviatoric part of the Cauchy stress tensor. We prove the existence of a weak solution to the problem, and under the additional assumption that the nonlinearity involved in the constitutive relation is Lipschitz continuous we also prove uniqueness of the weak solution. We then construct mixed finite element approximations of the system using both conforming and nonconforming finite element spaces. For both of these we prove the convergence of the method to the unique weak solution of the problem, and in the case of the conforming method we provide a bound on the error between the analytical solution and its finite element approximation in terms of the best approximation error from the finite element spaces. We propose first a Lions–Mercier type iterative method and next a classical fixed-point algorithm to solve the finite-dimensional problems resulting from the finite element discretisation of the system of nonlinear partial differential equations under consideration and present numerical experiments that illustrate the practical performance of the proposed numerical method.more » « less
-
Surface Stokes and Navier–Stokes equations are used to model fluid flow on surfaces. They have attracted significant recent attention in the numerical analysis literature because approximation of their solutions poses significant challenges not encountered in the Euclidean context. One challenge comes from the need to simultaneously enforce tangentiality and $H^1$ conformity (continuity) of discrete vector fields used to approximate solutions in the velocity-pressure formulation. Existing methods in the literature all enforce one of these two constraints weakly either by penalization or by use of Lagrange multipliers. Missing so far is a robust and systematic construction of surface Stokes finite element spaces which employ nodal degrees of freedom, including MINI, Taylor–Hood, Scott–Vogelius, and other composite elements which can lead to divergence-conforming or pressure-robust discretizations. In this paper we construct surface MINI spaces whose velocity fields are tangential. They are not $H^1$-conforming, but do lie in $H(div)$ and do not require penalization to achieve optimal convergence rates. We prove stability and optimal-order energy-norm convergence of the method and demonstrate optimal-order convergence of the velocity field in $L_2$ via numerical experiments. The core advance in the paper is the construction of nodal degrees of freedom for the velocity field. This technique also may be used to construct surface counterparts to many other standard Euclidean Stokes spaces, and we accordingly present numerical experiments indicating optimal-order convergence of nonconforming tangential surface Taylor–Hood elements.more » « less
-
We construct direct serendipity finite elements on general cuboidal hexahedra, which are 𝐻1-conforming and optimally approximate to any order. The new finite elements are direct in that the shape functions are directly defined on the physical element. Moreover, they are serendipity by possessing a minimal number of degrees of freedom satisfying the conformity requirement. Their shape function spaces consist of polynomials plus (generally nonpolynomial) supplemental functions, where the polynomials are included for the approximation property and supplements are added to achieve 𝐻1 -conformity. The finite elements are fully constructive. The shape function spaces of higher order 𝑟 ≥ 3 are developed first, and then the lower order spaces are constructed as subspaces of the third order space. Under a shape regularity assumption, and a mild restriction on the choice of supplemental functions, we develop the convergence properties of the new direct serendipity finite elements. Numerical results with different choices of supplements are compared on two mesh sequences, one regularly distorted and the other one randomly distorted. They all possess a convergence rate that aligns with the theory, while a slight difference lies in their performance.more » « less
-
We construct direct serendipity finite elements on general cuboidal hexahedra, which are H1-conforming and optimally approximate to any order. The new finite elements are direct in that the shape functions are directly defined on the physical element. Moreover, they are serendipity by possessing a minimal number of degrees of freedom satisfying the conformity requirement. Their shape function spaces consist of polynomials plus (generally nonpolynomial) supplemental functions, where the polynomials are included for the approximation property and supplements are added to achieve H1-conformity. The finite elements are fully constructive. The shape function spaces of higher order 𝑟 ≥ 3 are developed first, and then the lower order spaces are constructed as subspaces of the third order space. Under a shape regularity assumption, and a mild restriction on the choice of supplemental functions, we develop the convergence properties of the new direct serendipity finite elements. Numerical results with different choices of supplements are compared on two mesh sequences, one regularly distorted and the other one randomly distorted. They all possess a convergence rate that aligns with the theory, while a slight difference lies in their performance.more » « less