Pierre Kerfriden

Pierre Kerfriden, Jean-Charles Passieux, Stephane Pierre-Alain Bordas

This paper proposes a novel technique to reduce the computational burden associated with the simulation of localised failure. The proposed methodology affords the simulation of damage initiation and propagation whilst concentrating the computational effort where it is most needed, i.e. in the localisation zones. To do so, a local/global technique is devised where the global (slave) problem (far from the zones undergoing severe damage and cracking) is solved for in a reduced space computed by the classical Proper Orthogonal Decomposition, while the local (master) degrees of freedom (associated with the part of the structure where most of the damage is taking place) are fully resolved. Both domains are coupled through a local/global technique. This method circumvents the difficulties associated with model order reduction for the simulation of highly non-linear mechanical failure and offers an alternative or complementary approach to the development of multiscale fracture simulators.

Susanne Claus, Pierre Kerfriden

In this article, we present a cut finite element method for two-phase Navier-Stokes flows. The main feature of the method is the formulation of a unified continuous interior penalty stabilisation approach for, on the one hand, stabilising advection and the pressure-velocity coupling and, on the other hand, stabilising the cut region. The accuracy of the algorithm is enhanced by the development of extended fictitious domains to guarantee a well defined velocity from previous time steps in the current geometry. Finally, the robustness of the moving-interface algorithm is further improved by the introduction of a curvature smoothing technique that reduces spurious velocities. The algorithm is shown to perform remarkably well for low capillary number flows, and is a first step towards flexible and robust CutFEM algorithms for the simulation of microfluidic devices.

Octavio A. González-Estrada, Juan José Ródenas, Stéphane P. A. Bordas et al.

Purpose: This paper aims at assessing the effect of (1) the statical admissibility of the recovered solution; (2) the ability of the recovered solution to represent the singular solution; on the accuracy, local and global effectivity of recovery-based error estimators for enriched finite element methods (e.g. the extended finite element method, XFEM). Design/methodology/approach: We study the performance of two recovery techniques. The first is a recently developed superconvergent patch recovery procedure with equilibration and enrichment (SPR-CX). The second is known as the extended moving least squares recovery (XMLS), which enriches the recovered solutions but does not enforce equilibrium constraints. Both are extended recovery techniques as the polynomial basis used in the recovery process is enriched with singular terms for a better description of the singular nature of the solution. Findings: Numerical results comparing the convergence and the effectivity index of both techniques with those obtained without the enrichment enhancement clearly show the need for the use of extended recovery techniques in Zienkiewicz-Zhu type error estimators for this class of problems. The results also reveal significant improvements in the effectivities yielded by statically admissible recovered solutions. Originality/value: This work shows that both extended recovery procedures and statical admissibility are key to an accurate assessment of the quality of enriched finite element approximations.

Olivier Allix, Pierre Kerfriden, Pierre Gosselet

In quasi-static nonlinear time-dependent analysis, the choice of the time discretization is a complex issue. The most basic strategy consists in determining a value of the load increment that ensures the convergence of the solution with respect to time on the base of preliminary simulations. In more advanced applications, the load increments can be controlled for instance by prescribing the number of iterations of the nonlinear resolution procedure, or by using an arc-length algorithm. These techniques usually introduce a parameter whose correct value is not easy to obtain. In this paper, an alternative procedure is proposed. It is based on the continuous control of the residual of the reference problem over time, whose measure is easy to interpret. This idea is applied in the framework of a multiscale domain decomposition strategy in order to perform 3D delamination analysis.

Peng Yu, Cosmin Anitescu, Satyendra Tomar et al.

This paper presents a novel methodology of local adaptivity for the frequency-domain analysis of the vibrations of Reissner-Mindlin plates. The adaptive discretization is based on the recently developed Geometry Independent Field approximaTion (GIFT) framework, which may be seen as a generalisation of the Iso-Geometric Analysis (IGA). Within the GIFT framework, we describe the geometry of the structure exactly with NURBS (Non-Uniform Rational B-Splines), whilst independently employing Polynomial splines over Hierarchical T-meshes (PHT)-splines to represent the solution field. The proposed strategy of local adaptivity, wherein a posteriori error estimators are computed based on inexpensive hierarchical $h-$refinement, aims to control the discretisation error within a frequency band. The approach sweeps from lower to higher frequencies, refining the mesh appropriately so that each of the free vibration mode within the targeted frequency band is sufficiently resolved. Through several numerical examples, we show that the GIFT framework is a powerful and versatile tool to perform local adaptivity in structural dynamics. We also show that the proposed adaptive local $h-$refinement scheme allows us to achieve significantly faster convergence rates than when using a uniform $h-$refinement.

Susanne Claus, Pierre Kerfriden

In this paper, we propose a novel unfitted finite element method for the simulation of multiple body contact. The computational mesh is generated independently of the geometry of the interacting solids, which can be arbitrarily complex. The key novelty of the approach is the combination of elements of the CutFEM technology, namely the enrichment of the solution field via the definition of overlapping fictitious domains with a dedicated penalty-type regularisation of discrete operators, and the LaTIn hybrid-mixed formulation of complex interface conditions. Furthermore, the novel P1-P1 discretisation scheme that we propose for the unfitted LaTIn solver is shown to be stable, robust and optimally convergent with mesh refinement. Finally, the paper introduces a high-performance 3D level-set/CutFEM framework for the versatile and robust solution of contact problems involving multiple bodies of complex geometries, with more than two bodies interacting at a single point.

Aurèle Goetz, Ali Riza Durmaz, Martin Müller et al.

Materials' microstructures are signatures of their alloying composition and processing history. Therefore, microstructures exist in a wide variety. As materials become increasingly complex to comply with engineering demands, advanced computer vision (CV) approaches such as deep learning (DL) inevitably gain relevance for quantifying microstrucutures' constituents from micrographs. While DL can outperform classical CV techniques for many tasks, shortcomings are poor data efficiency and generalizability across datasets. This is inherently in conflict with the expense associated with annotating materials data through experts and extensive materials diversity. To tackle poor domain generalizability and the lack of labeled data simultaneously, we propose to apply a sub-class of transfer learning methods called unsupervised domain adaptation (UDA). These algorithms address the task of finding domain-invariant features when supplied with annotated source data and unannotated target data, such that performance on the latter distribution is optimized despite the absence of annotations. Exemplarily, this study is conducted on a lath-shaped bainite segmentation task in complex phase steel micrographs. Here, the domains to bridge are selected to be different metallographic specimen preparations (surface etchings) and distinct imaging modalities. We show that a state-of-the-art UDA approach surpasses the naïve application of source domain trained models on the target domain (generalization baseline) to a large extent. This holds true independent of the domain shift, despite using little data, and even when the baseline models were pre-trained or employed data augmentation. Through UDA, mIoU was improved over generalization baselines from 82.2%, 61.0%, 49.7% to 84.7%, 67.3%, 73.3% on three target datasets, respectively. This underlines this techniques' potential to cope with materials variance.

Vasilis Krokos, Viet Bui Xuan, Stéphane P. A. Bordas et al.

Multiscale computational modelling is challenging due to the high computational cost of direct numerical simulation by finite elements. To address this issue, concurrent multiscale methods use the solution of cheaper macroscale surrogates as boundary conditions to microscale sliding windows. The microscale problems remain a numerically challenging operation both in terms of implementation and cost. In this work we propose to replace the local microscale solution by an Encoder-Decoder Convolutional Neural Network that will generate fine-scale stress corrections to coarse predictions around unresolved microscale features, without prior parametrisation of local microscale problems. We deploy a Bayesian approach providing credible intervals to evaluate the uncertainty of the predictions, which is then used to investigate the merits of a selective learning framework. We will demonstrate the capability of the approach to predict equivalent stress fields in porous structures using linearised and finite strain elasticity theories.

Susanne Claus, Samuel Bigot, Pierre Kerfriden

In this article, we develop a cut finite element method for one-phase Stefan problems, with applications in laser manufacturing. The geometry of the workpiece is represented implicitly via a level set function. Material above the melting/vaporisation temperature is represented by a fictitious gas phase. The moving interface between the workpiece and the fictitious gas phase may cut arbitrarily through the elements of the finite element mesh, which remains fixed throughout the simulation, thereby circumventing the need for cumbersome re-meshing operations. The primal/dual formulation of the linear one-phase Stefan problem is recast into a primal non-linear formulation using a Nitsche-type approach, which avoids the difficulty of constructing inf-sup stable primal/dual pairs. Through the careful derivation of stabilisation terms, we show that the proposed Stefan-Signorini-Nitsche CutFEM method remains stable independently of the cut location. In addition, we obtain optimal convergence with respect to space and time refinement. Several 2D and 3D examples are proposed, highlighting the robustness and flexibility of the algorithm, together with its relevance to the field of micro-manufacturing.

Navid Rahimi, Pierre Kerfriden, Frank C Langbein et al.

Simplifying the geometry of a CAD model using defeaturing techniques enables more efficient discretisation and subsequent simulation for engineering analysis problems. Understanding the effect this simplification has on the solution helps to decide whether the simplification is suitable for a specific simulation problem. It can also help to understand the functional effect of a geometry feature. The effect of the simplification is quantified by a user-defined quantity of interest which is assumed to be (approximately) linear in the solution. A bound on the difference between the quantity of interest of the original and simplified solutions based on the energy norm is derived. The approach is presented in the context of electrostatics problems, but can be applied in general to a range of elliptic partial differential equations. Numerical results on the efficiency of the bound are provided for electrostatics problems with simplifications involving changes inside the problem domain as well as changes to the boundaries.