Researcher profile

Stéphane P. A. Bordas

Stéphane P. A. Bordas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

A Bayesian multiscale CNN framework to predict local stress fields in structures with microscale features

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.

preprint2022arXiv

Design of metamaterial-based heat manipulators by isogeometric shape optimization

There has been a growing interest in controlled heat flux manipulation to increase the efficiency of thermal apparatus. Heat manipulators control and manipulate heat flow. A key to the effective performance of these heat manipulators is their thermal design. Such designs can be achieved by a periodic assembly of unit cells (known as metamaterials or meta-structure), whose geometry and material properties can be optimized for a specific objective. In this work, we focus on thermal metamaterial-based heat manipulators such as thermal concentrator (which concentrates the heat flux in a specified region of the domain). The main scope of the current work is to optimize the shape of the heat manipulators using Particle Swarm Optimization (PSO) method. The geometry is defined using NURBS basis functions due to the higher smoothness and continuity and the thermal boundary value problem is solved using Isogeometric Analysis (IGA). Often, nodes as design variables (as in Lagrange finite element method) generate the serrate shapes of boundaries which need to be smoothened later. For the NURBS-based boundary with the control points as design variables, the required smoothness can be predefined through knot vectors and smoothening in the post-processing can be avoided. The optimized shape generated by PSO is compared with the other shape exploited in the literature. The effects of the number of design variables, the thermal conductivity of the materials used, as well as some of the geometry parameters on the optimum shapes are also demonstrated.

preprint2022arXiv

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

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.

preprint2022arXiv

Probabilistic Deep Learning for Real-Time Large Deformation Simulations

For many novel applications, such as patient-specific computer-aided surgery, conventional solution techniques of the underlying nonlinear problems are usually computationally too expensive and are lacking information about how certain can we be about their predictions. In the present work, we propose a highly efficient deep-learning surrogate framework that is able to accurately predict the response of bodies undergoing large deformations in real-time. The surrogate model has a convolutional neural network architecture, called U-Net, which is trained with force-displacement data obtained with the finite element method. We propose deterministic and probabilistic versions of the framework. The probabilistic framework utilizes the Variational Bayes Inference approach and is able to capture all the uncertainties present in the data as well as in the deep-learning model. Based on several benchmark examples, we show the predictive capabilities of the framework and discuss its possible limitations

preprint2022arXiv

SOniCS: Develop intuition on biomechanical systems through interactive error controlled simulations

This new approach allows the user to experiment with model choices easily and quickly without requiring in-depth expertise, as constitutive models can be modified by one line of code only. This ease in building new models makes SOniCS ideal to develop surrogate, reduced order models and to train machine learning algorithms for uncertainty quantification or to enable patient-specific simulations. SOniCS is thus not only a tool that facilitates the development of surgical training simulations but also, and perhaps more importantly, paves the way to increase the intuition of users or otherwise non-intuitive behaviors of (bio)mechanical systems. The plugin uses new developments of the FEniCSx project enabling automatic generation with FFCx of finite element tensors such as the local residual vector and Jacobian matrix. We validate our approach with numerical simulations such as manufactured solutions, cantilever beams, and benchmarks provided by FEBio. We reach machine precision accuracy and demonstrate the use of the plugin for a real-time haptic simulation involving a surgical tool controlled by the user in contact with a hyperelastic liver. We include complete examples showing the use of our plugin for simulations involving Saint Venant-Kirchhoff, Neo-Hookean, Mooney-Rivlin, and Holzapfel Ogden anisotropic models as supplementary material.

preprint2021arXiv

Colossal Atomic Force Response in van der Waals Materials Arising From Electronic Correlations

Understanding static and dynamic phenomena in complex materials at different length scales requires reliably accounting for van der Waals (vdW) interactions, which stem from long-range electronic correlations. While the important role of many-body vdW interactions has been extensively documented when it comes to the stability of materials, much less is known about the coupling between vdW interactions and atomic forces. Here we analyze the Hessian force response matrix for a single and two vdW-coupled atomic chains to show that a many-body description of vdW interactions yields atomic force response magnitudes that exceed the expected pairwise decay by 3-5 orders of magnitude for a wide range of separations between the perturbed and the observed atom. Similar findings are confirmed for graphene and carbon nanotubes. This colossal force enhancement suggests implications for phonon spectra, free energies, interfacial adhesion, and collective dynamics in materials with many interacting atoms.

preprint2021arXiv

Inverse deformation analysis: an experimental and numerical assessment using the FEniCS Project

In this paper, we develop a framework for solving inverse deformation problems using the FEniCS Project finite element software. We validate our approach with experimental imaging data acquired from a soft silicone beam under gravity. In contrast with inverse iterative algorithms that require multiple solutions of a standard elasticity problem, the proposed method can compute the undeformed configuration by solving only one modified elasticity problem. This modified problem has a complexity comparable to the standard one. The framework is implemented within an open-source pipeline enabling the direct and inverse deformation simulation directly from imaging data. We use the high-level Unified Form Language (UFL) of the FEniCS Project to express the finite element model in variational form and to automatically derive the consistent Jacobian. Consequently, the design of the pipeline is flexible: for example, it allows the modification of the constitutive models by changing a single line of code. We include a complete working example showing the inverse deformation of a beam deformed by gravity as supplementary material.

preprint2018arXiv

Quantifying discretization errors for soft-tissue simulation in computer assisted surgery: a preliminary study

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.