Volume 31, Issue 2

We consider flux-corrected finite element discretizations of 3D convection-dominated transport problems and assess the computational efficiency of algorithms based on such approximations. The methods under investigation include flux-corrected transport schemes and monolithic limiters. We discretize in space using a continuous Galerkin method and ℙ 1 or ℚ 1 finite elements. Time integration is performed using the Crank–Nicolson method or an explicit strong stability preserving Runge–Kutta method. Nonlinear systems are solved using a fixed-point iteration method, which requires solution of large linear systems at each iteration or time step. The great variety of options in the choice of discretization methods and solver components calls for a dedicated comparative study of existing approaches. To perform such a study, we define new 3D test problems for time dependent and stationary convection–diffusion–reaction equations. The results of our numerical experiments illustrate how the limiting technique, time discretization and solver impact on the overall performance.

This paper constructs and analyzes a boundary correction finite element method for the Stokes problem based on the Scott–Vogelius pair on Clough–Tocher splits. The velocity space consists of continuous piecewise polynomials of degree k , and the pressure space consists of piecewise polynomials of degree ( k – 1) without continuity constraints. A Lagrange multiplier space that consists of continuous piecewise polynomials with respect to the boundary partition is introduced to enforce boundary conditions and to mitigate the lack of pressure-robustness. We prove several inf-sup conditions, leading to the well-posedness of the method. In addition, we show that the method converges with optimal order and the velocity approximation is divergence-free.

In this paper a time-explicit weak Galerkin finite element method is introduced and analyzed for parabolic equations. The main idea relies on the inclusion of a stabilization term in the temporal direction in addition to the usual static stabilization in the weak Galerkin framework. Both semi-discrete and fully-discrete schemes in time are presented, as well as their stability and error analysis. Numerical results are reported for this new explicit weak Galerkin finite element method.

We introduce and analyse a fully discrete approximation for a mathematical model for the solidification and liquidation of materials of negligible specific heat. The model is a two-sided Mullins–Sekerka problem. The discretization uses finite elements in space and an independent parameterization of the moving free boundary. We prove unconditional stability and exact volume conservation for the introduced scheme. Several numerical simulations, including for nearly crystalline surface energies, demonstrate the practicality and accuracy of the presented numerical method.