Nitsche’s method for two and three dimensional NURBS patch coupling

in Computational Mechanics (in press)

We present a Nitche’s method to couple non-conforming two and three-dimensional NURBS (Non Uniform Rational B-splines) patches in the context of isogeometric analysis (IGA). We present results for linear ... [more ▼]

We present a Nitche’s method to couple non-conforming two and three-dimensional NURBS (Non Uniform Rational B-splines) patches in the context of isogeometric analysis (IGA). We present results for linear elastostatics in two and and three-dimensions. The method can deal with surface-surface or volume-volume coupling, and we show how it can be used to handle heterogeneities such as inclusions. We also present preliminary results on modal analysis. This simple coupling method has the potential to increase the applicability of NURBS-based isogeometric analysis for practical applications. [less ▲]

Convolution, aggregation and attention based deep neural networks for accelerating simulations in mechanics

in Frontiers in Materials (2023)

Deep learning surrogate models are being increasingly used in accelerating scientific simulations as a replacement for costly conventional numerical techniques. However, their use remains a significant ... [more ▼]

Deep learning surrogate models are being increasingly used in accelerating scientific simulations as a replacement for costly conventional numerical techniques. However, their use remains a significant challenge when dealing with real-world complex examples. In this work, we demonstrate three types of neural network architectures for efficient learning of highly non-linear deformations of solid bodies. The first two architectures are based on the recently proposed CNN U-NET and MAgNET (graph U-NET) frameworks which have shown promising performance for learning on mesh-based data. The third architecture is Perceiver IO, a very recent architecture that belongs to the family of attention-based neural networks–a class that has revolutionised diverse engineering fields and is still unexplored in computational mechanics. We study and compare the performance of all three networks on two benchmark examples, and show their capabilities to accurately predict the non-linear mechanical responses of soft bodies. [less ▲]

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

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 ▲]

Probabilistic Deep Learning for Real-Time Large Deformation Simulations

in Computer Methods in Applied Mechanics and Engineering (2022), 398(0045-7825), 115307

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 ... [more ▼]

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. [less ▲]

SOniCS: Interfacing SOFA and FEniCS for advanced constitutive models

Scientific Conference (2022, August)

The Simulation Open Framework Architecture (SOFA) is a software environment for building simulations with a particular focus on real-time medical applications, e.g. surgery. Its scope is far broader than ... [more ▼]

The Simulation Open Framework Architecture (SOFA) is a software environment for building simulations with a particular focus on real-time medical applications, e.g. surgery. Its scope is far broader than the FEniCS Project, encompassing e.g. rigid body dynamics, interfacing with haptic devices, contact and visualisation. Naturally, it also includes some finite element models of soft tissue mechanics, but these capabilities are currently ‘pre-baked’ and limited to a few simple constitutive models. The goal of this work is to incorporate state-of-the-art code generation tools from the FEniCS Project into SOFA in order to hugely increase SOFA’s capabilities in terms of soft tissue mechanics. To this end we have developed a new SOFA plugin named SOniCS. For adding a new material model in SOniCS, the user describes its strain energy density function using UFL (Unified Form Language) syntax. Then, using FFCx (FEniCSx Form Compiler) we generate the C code associated with the kernels corresponding to the automatically differentiated cell-local residual and stiffness forms. Finally, we assemble these kernels in SOFA into global tensors and solve the resulting non-linear systems of equations. The result is that it is now possible to straightforwardly implement complex material models such as the Holzapfel-Ogden anisotropic model into SOFA, and to use them alongside SOFA’s existing strong feature set in medical simulation. [less ▲]

Oncology and mechanics: landmark studies and promising clinical applications

in Advances in Applied Mechanics (2022), 55

Clinical management of cancer has continuously evolved for several decades. Biochemical, molecular and genomics approaches have brought and still bring numerous insights into cancerous diseases. It is now ... [more ▼]

Clinical management of cancer has continuously evolved for several decades. Biochemical, molecular and genomics approaches have brought and still bring numerous insights into cancerous diseases. It is now accepted that some phenomena, allowed by favorable biological conditions, emerge via mechanical signaling at the cellular scale and via mechanical forces at the macroscale. Mechanical phenomena in cancer have been studied in-depth over the last decades, and their clinical applications are starting to be understood. If numerous models and experimental setups have been proposed, only a few have led to clinical applications. The objective of this contribution is to propose to review a large scope of mechanical findings which have consequences on the clinical management of cancer. This review is mainly addressed to doctoral candidates in mechanics and applied mathematics who are faced with the challenge of the mechanics-based modeling of cancer with the aim of clinical applications. We show that the collaboration of the biological and mechanical approaches has led to promising advances in terms of modeling, experimental design and therapeutic targets. Additionally, a specific focus is brought on imaging-informed mechanics-based models, which we believe can further the development of new therapeutic targets and the advent of personalized medicine. We study in detail several successful workflows on patient-specific targeted therapies based on mechanistic modeling. [less ▲]

A CutFEM Method for a Mechanistic Modelling of Astrocytic Metabolism in 3D Physiological Morphologies

Scientific Conference (2022, June 07)

Investigating neurodegenerative diseases can be done complementary through biological and computational experiments. A good computational approach describing a simplification of the reality and focusing ... [more ▼]

Investigating neurodegenerative diseases can be done complementary through biological and computational experiments. A good computational approach describing a simplification of the reality and focusing only on some features of the problem can help getting insights on the field. The question addressed in our work is the role of astrocytes in neurodegeneration. These cells have two interesting characteristics that we want to investigate in our model: first, their role as metabolic mediator between neurons and blood vessels and second, their peculiar morphology. In fact, metabolic dysfunctions and morphological changes have been noticed in astrocyte affected by neuropathology. Computationally the main difficulty arising from solving a metabolic model into cellular shape comes from the complexity of the domain. The shape of astrocytes are very ramified, with thin branches and sharp edges. As shown in our previous work \cite{Farina}, a \cutfem{} \cite{Burman} approach is a suitable tool to deal with this issue. In our latest work we use real human three-dimensional astrocyte morphologies obtained via microscopy \cite{Salamanca} as domain to solve our system. The performed simulations highlight the effect of morphological changes on the system output. Suggesting that our model can be crucial in understanding the morphological-dependency in neuropathologies and that the spatial component cannot be neglected. [less ▲]

Predicting depression in old age: Combining life course data with machine learning

Scientific Conference (2022, June 03)

Colossal Enhancement of Atomic Force Response in van der Waals Materials Arising from Many-Body Electronic Correlations

in Physical Review Letters (2022)

An adaptive multiscale quasicontinuum approach for mechanical simulations of elastoplastic periodic lattices

in Mechanics Research Communications (2022), 126

The quasicontinuum method is a multiscale method that combines locally supported coarse-grained domains, with small regions in which the microstructural model is fully resolved. This contribution proposes ... [more ▼]

The quasicontinuum method is a multiscale method that combines locally supported coarse-grained domains, with small regions in which the microstructural model is fully resolved. This contribution proposes the first adaptive formulation of the method for microstructural elastoplasticity. The microstructural model uses an elastoplastic beam description. The indicator for refinement is the occurrence of plastic deformation, such that plasticity can only occur in fully resolved regions. An illustrative numerical example of a scratch test of an elastoplastic Kelvin lattice demonstrates the capabilities of the resulting framework. [less ▲]