HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems

Computational Engineering

This Workshop treated a variety of finite element methods and applications in computational engineering and expanded their mathematical foundation in engineering analysis. Among the 53 participants were mathematicians and engineers with focus on mixed and nonstandard finite element schemes and their applications

Computational Engineering

The focus of this Computational Engineering Workshop was on the mathematical foundation of state-of-the-art and emerging finite element methods in engineering analysis. The 52 participants included mathematicians and engineers with shared interest on discontinuous Galerkin or Petrov-Galerkin methods and other generalized nonconforming or mixed finite element methods

A variational framework for mathematically nonsmooth problems in solid and structure mechanics

This dissertation presents a new paradigm for addressing multi-physics problems with interfaces in the field of Additive Manufacturing and the modeling of fibrous composite materials. The unique process of adding the material layer by layer in the AM techniques raises the issue about the stability of the interfaces between the layers and along the boundaries of multi-constituent materials. A stabilized interface formulation is developed to model debonding in monotonic loading, fatigue effects in cyclic loading, and thermal effects at interfaces which severely impact the functional life of those materials and structures. The formulation is based on embedding Discontinuous Galkerin (DG) ideas in a Continuous Galerkin (CG) framework. Starting from a mixed method incorporating the Lagrange multiplier along the interface, a pure displacement formulation is derived using the Variational Multiscale Method (VMS). From a mathematical and computational perspective, the key factor influencing the accuracy and robustness of the interface formulation is the design of the numerical flux and the penalty or stability terms. Analytical expressions that are free from user-defined parameters are naturally derived for the numerical flux and stability tensor which are functions of the evolving geometric and material nonlinearity. The proposed framework is extended for debonding at finite strains across general bimaterial interfaces. An interfacial gap function is introduced that evolves subject to constraints imposed by opening and/or sliding interfaces. An internal variable formalism is derived together with the notion of irreversibility of damage results in a set of evolution equations for the gap function that seamlessly tracks interface debonding by treating damage and friction in a unified way. Tension debonding, compression damage, and frictional sliding are accommodated, and return mapping algorithms in the presence of evolving strong discontinuities are developed. This derivation variationally embeds the interfacial kinematic models that are crucial to capturing the physical and mathematical properties involving large strains and damage. The framework is extended for monolithic coupling of thermomechanical fields in the class of problems that have embedded weak and strong discontinuities in the mechanical and thermal fields. Since the derived expressions are a function of the mechanical and thermal fields, the resulting stabilized formulation contains numerical flux and stability tensors that provide an avenue to variationally embed interfacial kinetic and kinematic models for more robust representation of interfacial physics. Representative numerical tests involving large strains and rotations, damage phenomena, and thermal effects are performed to confirm the robustness and accuracy of the method. Comparison of the results with both experimental and numerical results from literature are presented.Ope

The Hybrid High-Order Method for Polytopal Meshes: Design, Analysis, and Applications

International audienceHybrid High-Order (HHO) methods are new generation numerical methods for models based on Partial Differential Equations with features that set them apart from traditional ones. These include: the support of polytopal meshes including non star-shaped elements and hanging nodes; the possibility to have arbitrary approximation orders in any space dimension; an enhanced compliance with the physics; a reduced computational cost thanks to compact stencil and static condensation. This monograph provides an introduction to the design and analysis of HHO methods for diffusive problems on general meshes, along with a panel of applications to advanced models in computational mechanics. The first part of the monograph lays the foundation of the method considering linear scalar second-order models, including scalar diffusion, possibly heterogeneous and anisotropic, and diffusion-advection-reaction. The second part addresses applications to more complex models from the engineering sciences: non-linear Leray-Lions problems, elasticity and incompressible fluid flows

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described