![]() Bulle, Raphaël ![]() ![]() in Computers and Mathematics with Applications (2023), 131 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 ▲] Detailed reference viewed: 115 (13 UL)![]() Bulle, Raphaël ![]() ![]() E-print/Working paper (2022) We develop a novel a posteriori error estimator for the L2 error committed by the finite ele- ment discretization of the solution of the fractional Laplacian. Our a posteriori error estimator takes ... [more ▼] We develop a novel a posteriori error estimator for the L2 error committed by the finite ele- ment discretization of the solution of the fractional Laplacian. Our a posteriori error estimator takes advantage of the semi–discretization scheme using a rational approximation which allows 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. 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. [less ▲] Detailed reference viewed: 69 (3 UL)![]() Bulle, Raphaël ![]() ![]() in Applied Mathematics Letters (2020), 107 We provide a new argument proving the reliability of the Bank-Weiser estimator for Lagrange piecewise linear finite elements in both dimension two and three. The extension to dimension three constitutes ... [more ▼] We provide a new argument proving the reliability of the Bank-Weiser estimator for Lagrange piecewise linear finite elements in both dimension two and three. The extension to dimension three constitutes the main novelty of our study. In addition, we present a numerical comparison of the Bank-Weiser and residual estimators for a three-dimensional test case. [less ▲] Detailed reference viewed: 165 (41 UL)![]() Bulle, Raphaël ![]() ![]() Scientific Conference (2020, July) Detailed reference viewed: 104 (8 UL)![]() ; ; Sensale, Marco ![]() 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 ▲] Detailed reference viewed: 42 (1 UL)![]() ; ; Sensale, Marco ![]() Scientific Conference (2019, December 09) The human skin behaves as an elastic membrane initially prestressed but not uniformly. The presence of anatomical sites favorable to the appearance of some tumors, a keloid in our case, while other sites ... [more ▼] The human skin behaves as an elastic membrane initially prestressed but not uniformly. The presence of anatomical sites favorable to the appearance of some tumors, a keloid in our case, while other sites never develop them attests to the importance of the mechanical environment of the tissue. Thus, a mechanical characterization of the tumored skin is necessary to understand the keloid expansion from a mechanical point of view. Our case study consists in modeling a bi-material structure composed of a keloid skin surrounded by healthy skin located on upper left arm of a young female. From the experimental measurements in vivo, by combining force sensor, displacement sensor and Digital Image Correlation techniques, we perform a mechanical analysis to characterize the mechanical stress fields over the entire area and on the interface ‘healthy skin/keloid skin’. Since the mechanical behavior of the tumorous skin is unknown, many physical models can be implemented and assessed very easily inside the specific digital software to fit with the real data. Once a set of mechanical parameters for both the healthy skin and the keloid skin are identified, the stress fields around the keloid are calculated. Next steps consist in determining matching preferential directions in order to define as precisely as possible the specifications of a device for preventing the growth of keloids. [less ▲] Detailed reference viewed: 48 (3 UL)![]() Hauseux, Paul ![]() ![]() ![]() Scientific Conference (2018, July 23) We present probabilistic approaches aiming at the selection of the best constitutive model and to identify their parameters from experimental data. These parameters are always associated with some degree ... [more ▼] We present probabilistic approaches aiming at the selection of the best constitutive model and to identify their parameters from experimental data. These parameters are always associated with some degree of uncertainty. It is therefore important to study how this statistical uncertainty in parameters propagates to a safety-critical quantity of interest in the output of a model. Efficient Monte Carlo methods based on variance reduction techniques (Sensitivity Derivatives Monte Carlo methods [Hauseux et al. 2017] and MultiLevel Monte Carlo [Giles 2015] methods) are employed to propagate this uncertainty for both random variables and random fields. Inverse and forward problems are strongly connected. In a bayesian setting [Matthies et al. 2017], developing methods that reduce the number of evaluations of the forward model to an absolute minimum to achieve convergence is crucial for tractable computations. Numerical results in the context of soft tissue biomechanics are presented and discussed. [less ▲] Detailed reference viewed: 267 (8 UL)![]() ; ; et al 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 ▲] Detailed reference viewed: 191 (1 UL)![]() Sensale, Marco ![]() Scientific Conference (2017, August) Detailed reference viewed: 122 (14 UL) |
||