Hierarchical a posteriori error estimation of Bank-Weiser type in the FEniCS Project

E-print/Working paper (2021)

In the seminal paper of Bank and Weiser [Math. Comp., 44 (1985), pp.283-301] a new a posteriori estimator was introduced. This estimator requires the solution of a local Neumann problem on every cell of ... [more ▼]

In the seminal paper of Bank and Weiser [Math. Comp., 44 (1985), pp.283-301] a new a posteriori estimator was introduced. This estimator requires the solution of a local Neumann problem on every cell of the finite element mesh. Despite the promise of Bank-Weiser type estimators, namely locality, computational efficiency, and asymptotic sharpness, they have seen little use in practical computational problems. The focus of this contribution is to describe a novel implementation of hierarchical estimators of the Bank-Weiser type in a modern high-level finite element software with automatic code generation capabilities. We show how to use the estimator to drive (goal-oriented) adaptive mesh refinement and to mixed approximations of the nearly-incompressible elasticity problems. We provide comparisons with various other used estimators. An open-source implementation based on the FEniCS Project finite element software is provided as supplementary material. [less ▲]

Practical aspects of the Bank-Weiser estimator implementation and Biomechanics applications.

Scientific Conference (2020, July)

Parameter identification problem in bimaterial human skin and sensitivity analysis : Uncertainties in biomechanics of skin

Scientific Conference (2019, December 09)

The proposed paper concerns the prediction of the numerical response of a biomechanical structure submitted to an unknown external loading state. The methodology is based on homogeneous and then ... [more ▼]

The proposed paper concerns the prediction of the numerical response of a biomechanical structure submitted to an unknown external loading state. The methodology is based on homogeneous and then heterogeneous structures such as healthy or pathological cutaneous tissues that can be mechanically tested in vivo under a patchy knowledge of boundary conditions. Experimental data corresponding to the extension of a piece of skin located between two pads with displacement enslavement, represent input data to the numerical model. Data are reaction force on one pad and displacement field between the two pads and all around. The numerical model consists of a representation of the bi-material domain geometry with neo-hookean behaviors. The boundary conditions and loadings of the experimental extension test are imposed. The materials parameters have been identified by inverse method starting from a constrained cost function minimizing the difference between the calculated displacements field and experimental displacements field obtained by digital image correlation and taking into account the reaction force as a constraint. An analysis of the model sensitivity to material parameters is presented. [less ▲]

Skew-symmetric Nitsche’s formulation in isogeometric analysis: Dirichlet and symmetry conditions, patch coupling and frictionless contact

in Computer Methods in Applied Mechanics and Engineering (2018), 341

A simple skew-symmetric Nitsche’s formulation is introduced into the framework of isogeometric analysis (IGA) to deal with various problems in small strain elasticity: essential boundary conditions ... [more ▼]

A simple skew-symmetric Nitsche’s formulation is introduced into the framework of isogeometric analysis (IGA) to deal with various problems in small strain elasticity: essential boundary conditions, symmetry conditions for Kirchhoff plates, patch coupling in statics and in modal analysis as well as Signorini contact conditions. For linear boundary or interface conditions, the skew-symmetric formulation is parameter-free. For contact conditions, it remains stable and accurate for a wide range of the stabilization parameter. Several numerical tests are performed to illustrate its accuracy, stability and convergence performance. We investigate particularly the effects introduced by Nitsche’s coupling, including the convergence performance and condition numbers in statics as well as the extra “outlier” frequencies and corresponding eigenmodes in structural dynamics. We present the Hertz test, the block test, and a 3D self-contact example showing that the skew-symmetric Nitsche’s formulation is a suitable approach to simulate contact problems in IGA. [less ▲]