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
- NSF-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
-
-
Abstract The classical serendipity and mixed finite element spaces suffer from poor approximation on nondegenerate, convex quadrilaterals. In this paper, we develop families of direct serendipity and direct mixed finite element spaces, which achieve optimal approximation properties and have minimal local dimension. The set of local shape functions for either the serendipity or mixed elements contains the full set of scalar or vector polynomials of degree r , respectively, defined directly on each element (i.e., not mapped from a reference element). Because there are not enough degrees of freedom for global $$H^1$$ H 1 or $$H(\text {div})$$ H ( div ) conformity, exactly two supplemental shape functions must be added to each element when $$r\ge 2$$ r ≥ 2 , and only one when $$r=1$$ r = 1 . The specific choice of supplemental functions gives rise to different families of direct elements. These new spaces are related through a de Rham complex. For index $$r\ge 1$$ r ≥ 1 , the new families of serendipity spaces $${\mathscr {DS}}_{r+1}$$ DS r + 1 are the precursors under the curl operator of our direct mixed finite element spaces, which can be constructed to have reduced or full $$H(\text {div})$$ H ( div ) approximation properties. One choice of direct serendipity supplements gives the precursor of the recently introduced Arbogast–Correa spaces (SIAM J Numer Anal 54:3332–3356, 2016. 10.1137/15M1013705 ). Other fully direct serendipity supplements can be defined without the use of mappings from reference elements, and these give rise in turn to fully direct mixed spaces. Our development is constructive, so we are able to give global bases for our spaces. Numerical results are presented to illustrate their properties.more » « less
-
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
-
more » « less
Abstract We construct new families of
direct serendipity anddirect mixed finite elements on general planar, strictly convex polygons that areH 1andH (div) conforming, respectively, and possess optimal order of accuracy for any order. They have a minimal number of degrees of freedom subject to the conformity and accuracy constraints. The name arises because the shape functions are defineddirectly on the physical elements, i.e., without using a mapping from a reference element. The finite element shape functions are defined to be the full spaces of scalar or vector polynomials plus a space of supplemental functions. The direct serendipity elements are the precursors of the direct mixed elements in a de Rham complex. The convergence properties of the finite elements are shown under a regularity assumption on the shapes of the polygons in the mesh, as well as some mild restrictions on the choices one can make in the construction of the supplemental functions. Numerical experiments on various meshes exhibit the performance of these new families of finite elements. -
A conforming discontinuous Galerkin finite element method for the Stokes problem on polytopal meshesmore » « less
Abstract A new discontinuous Galerkin finite element method for the Stokes equations is developed in the primary velocity‐pressure formulation. This method employs discontinuous polynomials for both velocity and pressure on general polygonal/polyhedral meshes. Most finite element methods with discontinuous approximation have one or more stabilizing terms for velocity and for pressure to guarantee stability and convergence. This new finite element method has the standard conforming finite element formulation, without any velocity or pressure stabilizers. Optimal‐order error estimates are established for the corresponding numerical approximation in various norms. The numerical examples are tested for low and high order elements up to the degree four in 2D and 3D spaces.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1142/S0218202523500112
