Abstract

Continuous Galerkin formulations are appealing due to their low computational cost, whereas discontinuous Galerkin formulation facilitate adaptative mesh refinement and are more accurate in regions with jumps of physical parameters. Since many electromagnetic problems involve materials with different physical properties, this last point is very important. For this reason, in this article we have developed a combined cG–dG formulation for Maxwell’s problem that allows arbitrary finite element spaces with functins continuous in patches of finite elements and discontinuous on the interfaces of these patches. In particular, the second formulation we propose comes from a novel continuous Galerkin formulation that reduces the amount of stabilization introduced in the numerical system. In all cases, we have performed stability and convergence analyses of the methods. The outcome of this work is a new approach that keeps the low CPU cost of recent nodal continuous formulations with the ability to deal with coefficient jumps and adaptivity of discontinuous ones. All these methods have been tested using a problem with singular solution and another one with different materials, in order to prove that in fact the resulting formulations can properly deal with these problems.

Abstract

 A new mixed finite element approximation of Maxwell's problem is proposed, its main features being that it is based on a novel augmented formulation of the continuous problem and the introduction of a mesh dependent stabilizing term, which yields a very weak control on the divergence of the unknown. The method is shown to be stable and convergent in the natural H(curl 0; Ω) norm for this unknown. In particular, convergence also applies to singular solutions, for which classical nodal-based interpolations are known to suffer from spurious convergence upon mesh refinement.

Abstract

We introduce in this paper a variational subgrid scale model for the solution of the incompressible Navier-Stokes equations. With respect to classical multiscale-based stabilisation techniques, we retain the subgrid scale effects in the convective term and integrate the subgrid scale equation in time. The method is applied to the Navier-Stokes equations in an accelerating frame of reference and with Dirichlet (essential), Neumann (natural) and mixed boundary conditions. The concrete objective of the paper is to test a numerical algorithm for solving the non-linear subgrid scale equation and the introduction of the subgrid scale into the grid scale equation. The performance of the technique is demonstrated through the solution of two numerical examples: one to test the tracking of the subgrid scale in the convection term and the other to investigate the effects of considering the subgrid scale transient.