A partition of unity approach to fluid mechanics and fluid-structure interaction

For problems involving large deformations of thin structures, simulating fluid-structure interaction (FSI) remains challenging largely due to the need to balance computational feasibility, efficiency, and solution accuracy. Overlapping domain techniques have been introduced as a way to combine the fluid-solid mesh conformity, seen in moving-mesh methods, without the need for mesh smoothing or re-meshing, which is a core characteristic of fixed mesh approaches. In this work, we introduce a novel overlapping domain method based on a partition of unity approach. Unified function spaces are defined as a weighted sum of fields given on two overlapping meshes. The method is shown to achieve optimal convergence rates and to be stable for steady-state Stokes, Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for FSI in the case of a 2D mock aortic valve simulation. These initial results point to the potential applicability of the method to a wide range of FSI applications, enabling boundary layer refinement and large deformations without the need for re-meshing or user-defined stabilization.Comment: 34 pages, 15 figur

Structure-preserving mesh coupling based on the Buffa-Christiansen complex

The state of the art for mesh coupling at nonconforming interfaces is presented and reviewed. Mesh coupling is frequently applied to the modeling and simulation of motion in electromagnetic actuators and machines. The paper exploits Whitney elements to present the main ideas. Both interpolation- and projection-based methods are considered. In addition to accuracy and efficiency, we emphasize the question whether the schemes preserve the structure of the de Rham complex, which underlies Maxwell's equations. As a new contribution, a structure-preserving projection method is presented, in which Lagrange multiplier spaces are chosen from the Buffa-Christiansen complex. Its performance is compared with a straightforward interpolation based on Whitney and de Rham maps, and with Galerkin projection.Comment: 17 pages, 7 figures. Some figures are omitted due to a restricted copyright. Full paper to appear in Mathematics of Computatio

Hybrid coupling of CG and HDG discretizations based on Nitsche’s method

This is a post-peer-review, pre-copyedit version of an article published in Computational mechanics. The final authenticated version is available online at: http://dx.doi.org/10.1007/s00466-019-01770-8A strategy to couple continuous Galerkin (CG) and hybridizable discontinuous Galerkin (HDG) discretizations based only on the HDG hybrid variable is presented for linear thermal and elastic problems. The hybrid CG-HDG coupling exploits the definition of the numerical flux and the trace of the solution on the mesh faces to impose the transmission conditions between the CG and HDG subdomains. The con- tinuity of the solution is imposed in the CG problem via Nitsche’s method, whereas the equilibrium of the flux at the interface is naturally enforced as a Neumann con- dition in the HDG global problem. The proposed strategy does not affect the core structure of CG and HDG discretizations. In fact, the resulting formulation leads to a minimally-intrusive coupling, suitable to be integrated in existing CG and HDG libraries. Numerical experiments in two and three dimensions show optimal global convergence of the stress and superconvergence of the displacement field, locking-free approximation, as well as the potential to treat structural problems of engineering interest featuring multiple materials with compressible and nearly incompressible behaviors.Peer ReviewedPostprint (author's final draft

Computational Electromagnetism and Acoustics

It is a moot point to stress the significance of accurate and fast numerical methods for the simulation of electromagnetic fields and sound propagation for modern technology. This has triggered a surge of research in mathematical modeling and numerical analysis aimed to devise and improve methods for computational electromagnetism and acoustics. Numerical techniques for solving the initial boundary value problems underlying both computational electromagnetics and acoustics comprise a wide array of different approaches ranging from integral equation methods to finite differences. Their development faces a few typical challenges: highly oscillatory solutions, control of numerical dispersion, infinite computational domains, ill-conditioned discrete operators, lack of strong ellipticity, hysteresis phenomena, to name only a few. Profound mathematical analysis is indispensable for tackling these issues. Many outstanding contributions at this Oberwolfach conference on Computational Electromagnetism and Acoustics strikingly confirmed the immense recent progress made in the field. To name only a few highlights: there have been breakthroughs in the application and understanding of phase modulation and extraction approaches for the discretization of boundary integral equations at high frequencies. Much has been achieved in the development and analysis of discontinuous Galerkin methods. New insight have been gained into the construction and relationships of absorbing boundary conditions also for periodic media. Considerable progress has been made in the design of stable and space-time adaptive discretization techniques for wave propagation. New ideas have emerged for the fast and robust iterative solution for discrete quasi-static electromagnetic boundary value problems

Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations

International audienceWe present equilibrated flux a posteriori error estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed finite element discretizations of the two-dimensional Poisson problem. Relying on the equilibration by mixed finite element solution of patchwise Neumann problems, the estimates are guaranteed, locally computable, locally efficient, and robust with respect to polynomial degree. Maximal local overestimation is guaranteed as well. Numerical experiments suggest asymptotic exactness for the incomplete interior penalty discontinuous Galerkin scheme

Vertex-based Compatible Discrete Operator schemes on polyhedral meshes for advection-diffusion equations

International audienceWe devise and analyze vertex-based, Péclet-robust, lowest-order schemes for advection-diffusion equations that support polyhedral meshes. The schemes are formulated using Compatible Discrete Operators (CDO), namely primal and dual discrete differential operators, a discrete contraction operator for advection, and a discrete Hodge operator for diffusion. Moreover, discrete boundary operators are devised to weakly enforce Dirichlet boundary conditions. The analysis sheds new light on the theory of Friedrichs' operators at the purely algebraic level. Moreover, an extension of the stability analysis hinging on inf-sup conditions is presented to incorporate divergence-free velocity fields under some assumptions. Error bounds and convergence rates for smooth solutions are derived, and numerical results are presented on three-dimensional polyhedral meshes