This talk will be devoted to the usage of new discretization schemes on polyhedral meshes in an industrial context. These discretizations called CDO [1, 2] (Compatible Discrete Operator) or Hybrid High Order [3,4] (HHO) schemes have been recently implemented in Code Saturne [5]. Code Saturne is an open-source code developed at EDF R&D aiming at simulating single-phase flows. First, the advantages of robust polyhedral discretizations will be recalled. Then, the underpinning principles of CDO schemes will be presented as well as some applications: diffusion equations, transport problems, groundwater flows or the discretization of the Stokes equations. High Performance Computing (HPC) aspects will be also discussed as it is an essential feature in an industrial context either to address complex and large computational domains or to get a quick answer. Some highlights on the main outlooks will be given to conclude.[-]

We present a lowest order Serendipity Virtual Element method, and show its use for the numerical solution of linear magneto-static problems in three dimensions. The method can be applied to very general decompositions of the computational domain (as is natural for Virtual Element Methods) and uses as unknowns the (constant) tangential component of the magnetic eld H on each edge, and the vertex values of the Lagrange multiplier p (used to enforce the solenoidality of the magnetic induction B = µH). In this respect the method can be seen as the natural generalization of the lowest order Edge Finite Element Method (the so-called ”first kind N´ed´elec” elements) to polyhedra of almost arbitrary shape, and as we show on some numerical examples it exhibits very good accuracy (for being a lowest order element) and excellent robustness with respect to distortions. Hints on a whole family of elements will also be given.[-]

Some convergence properties for the approximation of second order elliptic problems with a variety of boundary conditions (homogeneous Dirichlet, homogeneous or non-homogeneous Neumann or Fourier boundary conditions), using a given discretisation method, can be obtained when this method is plugged into the Gradient Discretisation Method (GDM) framework.
Instead of defining one GDM framework for each of these boundary conditions, we show that these properties can be stated using the same abstract tools for all the above boundary conditions. Then these tools enable the application of the GDM to a larger class of elliptic problems.[-]

Recently there has been a surge in interest in cut, or unfitted, finite element methods. In this class of methods typically the computational mesh is independent of the geometry. Interfaces and boundaries are allowed to cut through the mesh in a very general fashion. Constraints on the boundaries such as boundary or transmission conditions are typically imposed weakly using Nitsche's method. In this talk we will discuss how these ideas can be combined in a fruitful way with the idea of hybridization, where additional degrees of freedom are added on the interfaces to further improve the decoupling of the systems, allowing for static condensation of interior unknowns. In the first part of the talk we will discuss how hybridization can be combined with the classical cut finite element method, using standard H1 -conforming finite elements in each subdomain, leading to a robust method allowing for the integration of polytopal geometries, where the subdomains are independent of the underlying mesh. This leads to a framework where it is easy to integrate multiscale features such as strongly varying coefficients, or multidimensional coupling, as in flow in fractured domains. Some examples of such applications will be given. In the second part of the talk we will focus on the Hybridized High Order Method (HHO) and show how cut techniques can be introduced in this context. The HHO is a recently introduced nonconforming method that allows for arbitrary order discretization of diffusive problems on polytopal meshes. HHO methods have hybrid unknowns, made of polynomials in the mesh elements and on the faces, without any continuity requirement. They rely on high-order local reconstructions, which are used to build consistent Galerkin contributions and appropriate stabilization terms designed to preserve the high-order approximation properties of the local reconstructions. Here we will show how cut element techniques can be introduced as a tool for the handling of (possibly curved) interfaces or boundaries that are allowed to cut through the polytopal mesh. In this context the cut element method plays the role of a local interface model, where the associated degrees of freedom are eliminated in the static condensation step. Issues of robustness and accuracy will be discussed and illustrated by some numerical examples.[-]

The aim of the talk is to introduce a nonlinear Discrete Duality Finite Volume scheme to approximate the solutions of drift-diffusion equations. The scheme is built to preserve at the discrete level even on severely distorted meshes the energy / energy dissipation relation. This relation is of paramount importance to capture the long-time behavior of the problem in an accurate way. To enforce it, the linear convection diffusion equation is rewritten in a nonlinear form before being discretized. This is a joint work with Clément Cancès (Lille) and Stella Krell (Nice).[-]