NAMar 20
H(curl)-based approximation of the Stokes problem with slip boundary conditionsWietse M. Boon, Ralf Hiptmair, Wouter Tonnon et al.
Reformulating the incompressible Stokes equations with the velocity sought in H(curl) has recently emerged as a promising approach for the design of helicity-preserving schemes in magnetohydrodynamics and pressure-robust finite element methods on polygonal meshes. A key challenge in this setting, however, is the treatment of Navier slip boundary conditions. In this paper, we overcome this difficulty by recasting the slip condition as a Robin boundary condition and proving well-posedness of the resulting continuous problem. We further identify the geometric and regularity assumptions on the domain and the exact solution under which the classical Stokes solution is recovered. Finally, we study a conforming finite element Galerkin discretization, establishing stability and a priori error estimates. Numerical experiments validate the optimal convergence rates predicted by the theory.
NAMay 8
On structure-preserving and pointwise conservative continuous DG schemes for hyperbolic systemsRémi Abgrall, Michael Dumbser, Pierre-Henri Maire et al.
We present a new class of structure-preserving semi-discrete continuous-discontinuous Galerkin (CG-DG) finite element schemes for linear and nonlinear hyperbolic systems of partial differential equations on unstructured simplex meshes that automatically satisfy the following properties: i) the new schemes are not only cellwise conservative, but also locally pointwise conservative everywhere, hence they satisfy the integral form of the conservation law on arbitrary control volumes that do not have to coincide with the mesh at all; ii) the new methods naturally satisfy the two basic vector calculus identities $\nabla \cdot \nabla \times \mathbf{A}$ and $\nabla \times \nabla Z$ exactly pointwise locally and globally everywhere on the discrete level; iii) for linear symmetric hyperbolic systems the schemes are naturally energy conservative for the square energy, i.e. nonlinearly stable in the $L^2$ norm. The key ingredient of the new CG-DG schemes is the use of two different but compatible approximation spaces: the classical DG space $\mathcal{U}_h^N$ of discontinuous piecewise polynomials of degree up to $N$ and a classical finite element space $\mathcal{W}_h^{N+1}$ of globally continuous piecewise polynomials of degree $N+1$. In the new CG-DG schemes, the discrete solution $\mathbf{u}_h$ is sought in $\mathcal{U}_h^N$, while a suitable discrete flux field $\tilde{\mathbf{f}}_h$ is computed in $\mathcal{W}_h^{N+1}$. For $N=0$ our new schemes are directly related to cell-centered finite volume schemes with suitable vertex-based fluxes. All claimed properties of the schemes are first mathematically proven and are then also verified via suitable numerical tests. We show applications of our approach to three linear and nonlinear hyperbolic systems.
NAApr 7
Robust H(curl)-based finite element methods for the incompressible MHD systemLourenço Beirão da Veiga, Sergio Gómez, Ilaria Perugia et al.
We propose and analyze a class of finite element methods for the time-dependent incompressible magnetohydrodynamics system based on $H(\mathrm{curl})$-conforming discretizations for both the velocity and the magnetic field. This choice is guided by the aim of developing methods that are also suitable for the types of solutions arising in problems posed on nonconvex domains. Within this framework, we introduce three stabilized formulations, and study how the stabilization mechanisms employed influence their structural properties. In particular, we focus on suitability for nonconvex polyhedral domains, the need for Lagrange multipliers for the magnetic field, pressure-robustness, and quasi-robustness with respect to both the fluid and magnetic Reynolds numbers. The proposed formulations are further assessed through numerical experiments, highlighting their practical performance.