Search

Error control by means of a posteriori error estimators or indica-tors and adaptive discretizations, such as adaptive mesh refinement, have emerged in the late seventies. Since then, numerous theoretical developments and improvements have been made, as well as the first attempts to introduce them into real-life industrial applications. The present introductory chapter provides an overview of the subject, highlights some of the achievements to date and discusses possible perspectives.

We discuss how slip conditions for the Stokes equation can be handled using Nitsche method, for a stabilized finite element discretization. Emphasis is made on the interplay between stabilization and Nitsche terms. Well-posedness of the discrete problem and optimal convergence rates, in natural norm for the velocity and the pressure, are established, and illustrated with various numerical experiments. The proposed method fits naturally in the context of a finite element implementation while being accurate, and allows an increased flexibility in the choice of the finite element pairs., Comment: 16 pages, 10 figures

The Finite Element Method (FEM) is a well-established procedure for computing approximate solutions to deterministic engineering problems described by partial differential equations. FEM produces discrete approximations of the solution with a discretisation error that can be an be quantified with \emph{a posteriori} error estimates. The practical relevance of error estimates for biomechanics problems, especially for soft tissue where the response is governed by large strains, is rarely addressed. In this contribution, we propose an implementation of \emph{a posteriori} error estimates targeting a user-defined quantity of interest, using the Dual Weighted Residual (DWR) technique tailored to biomechanics. The proposed method considers a general setting that encompasses three-dimensional geometries and model non-linearities, which appear in hyperelastic soft tissues. We take advantage of the automatic differentiation capabilities embedded in modern finite element software, which allows the error estimates to be computed generically for a large class of models and constitutive laws. First we validate our methodology using experimental measurements from silicone samples, and then illustrate its applicability for patient-specific computations of pressure ulcers on a human heel.

We develop a novel a posteriori error estimator for the $L^2$ error committed by the finite element discretization of the solution of the fractional Laplacian. Our a posteriori error estimator takes advantage of the semi-discretization scheme using rational approximations which allow to reformulate the fractional problem into a family of non-fractional parametric problems. The estimator involves applying the implicit Bank-Weiser error estimation strategy to each parametric non-fractional problem and reconstructing the fractional error through the same rational approximation used to compute the solution to the original fractional problem. In addition we propose an algorithm to adapt both the finite element mesh and the rational scheme in order to balance the discretization errors. We provide several numerical examples in both two and three-dimensions demonstrating the effectivity of our estimator for varying fractional powers and its ability to drive an adaptive mesh refinement strategy.

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.

Errors in biomechanics simulations arise from modeling and discretization. Modeling errors are due to the choice of the mathematical model whilst discretization errors measure the impact of the choice of the numerical method on the accuracy of the approximated solution to this specific mathematical model. A major source of discretization errors is mesh generation from medical images, that remains one of the major bottlenecks in the development of reliable, accurate, automatic and efficient personalized, clinically-relevant Finite Element (FE) models in biomechanics. The impact of mesh quality and density on the accuracy of the FE solution can be quantified with \emph{a posteriori} error estimates. Yet, to our knowledge, the relevance of such error estimates for practical biomechanics problems has seldom been addressed, see [25]. In this contribution, we propose an implementation of some a posteriori error estimates to quantify the discretization errors and to optimize the mesh. More precisely, we focus on error estimation for a user-defined quantity of interest with the Dual Weighted Residual (DWR) technique. We test its applicability and relevance in two situations, corresponding to computations for a tongue and an artery, using a simplified setting, i.e., plane linearized elasticity with contractility of the soft-tissue modeled as a pre-stress. Our results demonstrate the feasibility of such methodology to estimate the actual solution errors and to reduce them economically through mesh refinement., Comment: 35 pages, 10 figures

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.

This paper is devoted to the partial null controllability issue of parabolic linear systems with n equations. Given a bounded domain in R N, we study the effect of m localized controls in a nonempty open subset only controlling p components of the solution (p, m < n). The first main result of this paper is a necessary and sufficient condition when the coupling and control matrices are constant. The second result provides, in a first step, a sufficient condition of partial null controllability when the matrices only depend on time. In a second step, through an example of partially controlled 2x2 parabolic system, we will provide positive and negative results on partial null controllability when the coefficients are space dependent.

We propose a new procedure for clustering nucleotide sequences based on the "Laplacian Eigenmaps" and Gaussian Mixture modelling. This proposal is then applied to a set of 100 DNA sequences from the mitochondrially encoded NADH dehydrogenase 3 (ND3) gene of a collection of Platyhelminthes and Nematoda species. The resulting clusters are then shown to be consistent with the gene phylogenetic tree computed using a maximum likelihood approach. This comparison shows in particular that the clustering produced by the methodology combining Laplacian Eigenmaps with Gaussian Mixture models is coherent with the phylogeny as well as with the NCBI taxonomy. We also developed a Python package for this procedure which is available online.

We present a convergence analysis of the penalty method applied to unilateral contact problems in two and three space dimensions. We first consider, under various regularity assumptions on the exact solution to the unilateral contact problem, the convergence of the continuous penalty solution as the penalty parameter $\varepsilon$ vanishes. Then, the analysis of the finite element discretized penalty method is carried out. Denoting by $h$ the discretization parameter, we show that the error terms we consider give the same estimates as in the case of the constrained problem when the penalty parameter is such that $\varepsilon = h$.

In this work a local projection stabilization method is proposed to solve a fictitious domain problem. The method adds a suitable fluctuation term to the formulation thus rendering the natural space for the Lagrange multiplier stable. Stability and convergence are proved and these results are illustrated by a numerical experiment., Comment: Submitted Preprint

We introduce and analyze a Nitsche-based domain decomposition method for the solution of hypersingular integral equations. This method allows for discretizations with non-matching grids without the necessity of a Lagrangian multiplier, as opposed to the traditional mortar method. We prove its almost quasi-optimal convergence and underline the theory by a numerical experiment., Comment: 21 pages, 5 figures

In the presented study, a numerical model which predicts the flow-induced collapse within the pharyngeal airway is validated using in vitro measurements. Theoretical simplifications were considered to limit the computation time. Systematic comparisons between simulations and measurements were performed on an in vitro replica, which reflects asymmetries of the geometry and of the tissue properties at the base of the tongue and in pathological conditions (strong initial obstruction). First, partial obstruction is observed and predicted. Moreover, the prediction accuracy of the numerical model is of 4.2% concerning the deformation (mean quadratic error on the constriction area). It shows the ability of the assumptions and method to predict accurately and quickly a fluid-structure interaction.

This paper gives a review of the experience provided by our group in terms of computer aided maxillo-facial surgery.

One of the most important issue in soft tissue modeling is to assess the quality of the simulations. A validation protocol is presented based on two CT scans of the patient acquired before and after cranio-maxillofacial surgery. The actual bones repositioning realized during the intervention are accurately measured and reproduced. A evaluation of the soft tissue deformation is then computed using a finite element model of the face. The simulations are therefore compared, qualitatively and quantitatively, with the actual outcome of the surgery. This protocol enable to rigorously evaluate different modeling methods, and to assess the clinical relevance of soft tissue simulation in maxillofacial surgery.

Obstructive sleep apnea syndrome is now considered as a major health care topic. An in-vitro setup which reproduces and simplifies upper airway geometry has been the basis to study the fluid/walls interaction that leads to an apnea. It consists of a rigid pipe (the pharynx) in contact with a deformable latex cylinder filled with water (the tongue). Air flows out of the rigid pipe and induces pressure forces on the cylinder. We present a numerical model of this setup: a finite element model of the latex cylinder is in interaction with a fluid model. Simulation of an hypopnea (partial collapsus of the airway) has been possible and in agreement with observations from the in-vitro setup. The same phenomenon has been simulated on a soft palate model obtained from a patient sagittal radiography. These first results encourage us to improve the model so as it could reproduce the complete apnea phenomenon, and be used for a planification purpose in sleep apnea surgery.