Showing 1–50 of 58 results for author: Bordas, S P A

Intrinsic mixed-dimensional beam-shell-solid couplings in linear Cosserat continua via tangential differential calculus

Authors: Adam Sky, Jack S. Hale, Andreas Zilian, Stéphane P. A. Bordas, Patrizio Neff

Abstract: We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this app… ▽ More We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this approach circumvents the need for intermediate finite elements or mortar methods. We introduce the derivations of all models of various dimensions using tangential differential calculus. The coupling itself is then achieved by defining a mixed-dimensional action functional with consistent Sobolev trace operators. Finally, we present numerical examples involving a three-dimensional silicone-rubber block reinforced with a curved graphite shell on its lower surface, a three-dimensional silver block reinforced with a graphite plate and beams, and lastly, intersecting silver shells reinforced with graphite beams. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Gaussian process regression + deep neural network autoencoder for probabilistic surrogate modeling in nonlinear mechanics of solids

Authors: Saurabh Deshpande, Hussein Rappel, Mark Hobbs, Stéphane P. A. Bordas, Jakub Lengiewicz

Abstract: Many real-world applications demand accurate and fast predictions, as well as reliable uncertainty estimates. However, quantifying uncertainty on high-dimensional predictions is still a severely under-invested problem, especially when input-output relationships are non-linear. To handle this problem, the present work introduces an innovative approach that combines autoencoder deep neural networks… ▽ More Many real-world applications demand accurate and fast predictions, as well as reliable uncertainty estimates. However, quantifying uncertainty on high-dimensional predictions is still a severely under-invested problem, especially when input-output relationships are non-linear. To handle this problem, the present work introduces an innovative approach that combines autoencoder deep neural networks with the probabilistic regression capabilities of Gaussian processes. The autoencoder provides a low-dimensional representation of the solution space, while the Gaussian process is a Bayesian method that provides a probabilistic mapping between the low-dimensional inputs and outputs. We validate the proposed framework for its application to surrogate modeling of non-linear finite element simulations. Our findings highlight that the proposed framework is computationally efficient as well as accurate in predicting non-linear deformations of solid bodies subjected to external forces, all the while providing insightful uncertainty assessments. △ Less

Submitted 15 July, 2024; originally announced July 2024.

A short perspective on a posteriori error control and adaptive discretizations

Authors: Roland Becker, Stéphane P. A. Bordas, Franz Chouly, Pascal Omnes

Abstract: 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 s… ▽ More 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. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Journal ref: Advances in Applied Mechanics, volume 58, chapter 1, 2024

Graph Representation Learning Strategies for Omics Data: A Case Study on Parkinson's Disease

Authors: Elisa Gómez de Lope, Saurabh Deshpande, Ramón Viñas Torné, Pietro Liò, Enrico Glaab, Stéphane P. A. Bordas

Abstract: Omics data analysis is crucial for studying complex diseases, but its high dimensionality and heterogeneity challenge classical statistical and machine learning methods. Graph neural networks have emerged as promising alternatives, yet the optimal strategies for their design and optimization in real-world biomedical challenges remain unclear. This study evaluates various graph representation learn… ▽ More Omics data analysis is crucial for studying complex diseases, but its high dimensionality and heterogeneity challenge classical statistical and machine learning methods. Graph neural networks have emerged as promising alternatives, yet the optimal strategies for their design and optimization in real-world biomedical challenges remain unclear. This study evaluates various graph representation learning models for case-control classification using high-throughput biological data from Parkinson's disease and control samples. We compare topologies derived from sample similarity networks and molecular interaction networks, including protein-protein and metabolite-metabolite interactions (PPI, MMI). Graph Convolutional Network (GCNs), Chebyshev spectral graph convolution (ChebyNet), and Graph Attention Network (GAT), are evaluated alongside advanced architectures like graph transformers, the graph U-net, and simpler models like multilayer perceptron (MLP). These models are systematically applied to transcriptomics and metabolomics data independently. Our comparative analysis highlights the benefits and limitations of various architectures in extracting patterns from omics data, paving the way for more accurate and interpretable models in biomedical research. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Submitted to Machine Learning in Computational Biology 2024 as an extended abstract, 2 pages + 1 appendix

Discovering deposition process regimes: leveraging unsupervised learning for process insights, surrogate modeling, and sensitivity analysis

Authors: Geremy Loachamín Suntaxi, Paris Papavasileiou, Eleni D. Koronaki, Dimitrios G. Giovanis, Georgios Gakis, Ioannis G. Aviziotis, Martin Kathrein, Gabriele Pozzetti, Christoph Czettl, Stéphane P. A. Bordas, Andreas G. Boudouvis

Abstract: This work introduces a comprehensive approach utilizing data-driven methods to elucidate the deposition process regimes in Chemical Vapor Deposition (CVD) reactors and the interplay of physical mechanism that dominate in each one of them. Through this work, we address three key objectives. Firstly, our methodology relies on process outcomes, derived by a detailed CFD model, to identify clusters of… ▽ More This work introduces a comprehensive approach utilizing data-driven methods to elucidate the deposition process regimes in Chemical Vapor Deposition (CVD) reactors and the interplay of physical mechanism that dominate in each one of them. Through this work, we address three key objectives. Firstly, our methodology relies on process outcomes, derived by a detailed CFD model, to identify clusters of "outcomes" corresponding to distinct process regimes, wherein the relative influence of input variables undergoes notable shifts. This phenomenon is experimentally validated through Arrhenius plot analysis, affirming the efficacy of our approach. Secondly, we demonstrate the development of an efficient surrogate model, based on Polynomial Chaos Expansion (PCE), that maintains accuracy, facilitating streamlined computational analyses. Finally, as a result of PCE, sensitivity analysis is made possible by means of Sobol' indices, that quantify the impact of process inputs across identified regimes. The insights gained from our analysis contribute to the formulation of hypotheses regarding phenomena occurring beyond the transition regime. Notably, the significance of temperature even in the diffusion-limited regime, as evidenced by the Arrhenius plot, suggests activation of gas phase reactions at elevated temperatures. Importantly, our proposed methods yield insights that align with experimental observations and theoretical principles, aiding decision-making in process design and optimization. By circumventing the need for costly and time-consuming experiments, our approach offers a pragmatic pathway towards enhanced process efficiency. Moreover, this study underscores the potential of data-driven computational methods for innovating reactor design paradigms. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Integrating supervised and unsupervised learning approaches to unveil critical process inputs

Authors: Paris Papavasileiou, Dimitrios G. Giovanis, Gabriele Pozzetti, Martin Kathrein, Christoph Czettl, Ioannis G. Kevrekidis, Andreas G. Boudouvis, Stéphane P. A. Bordas, Eleni D. Koronaki

Abstract: This study introduces a machine learning framework tailored to large-scale industrial processes characterized by a plethora of numerical and categorical inputs. The framework aims to (i) discern critical parameters influencing the output and (ii) generate accurate out-of-sample qualitative and quantitative predictions of production outcomes. Specifically, we address the pivotal question of the sig… ▽ More This study introduces a machine learning framework tailored to large-scale industrial processes characterized by a plethora of numerical and categorical inputs. The framework aims to (i) discern critical parameters influencing the output and (ii) generate accurate out-of-sample qualitative and quantitative predictions of production outcomes. Specifically, we address the pivotal question of the significance of each input in shaping the process outcome, using an industrial Chemical Vapor Deposition (CVD) process as an example. The initial objective involves merging subject matter expertise and clustering techniques exclusively on the process output, here, coating thickness measurements at various positions in the reactor. This approach identifies groups of production runs that share similar qualitative characteristics, such as film mean thickness and standard deviation. In particular, the differences of the outcomes represented by the different clusters can be attributed to differences in specific inputs, indicating that these inputs are critical for the production outcome. Leveraging this insight, we subsequently implement supervised classification and regression methods using the identified critical process inputs. The proposed methodology proves to be valuable in scenarios with a multitude of inputs and insufficient data for the direct application of deep learning techniques, providing meaningful insights into the underlying processes. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Quantum-informed simulations for mechanics of materials: DFTB+MBD framework

Authors: Zhaoxiang Shen, Raúl I. Sosa, Stéphane P. A. Bordas, Alexandre Tkatchenko, Jakub Lengiewicz

Abstract: The macroscopic behaviors of materials are determined by interactions that occur at multiple lengths and time scales. Depending on the application, describing, predicting, and understanding these behaviors require models that rely on insights from electronic and atomic scales. In such cases, classical simplified approximations at those scales are insufficient, and quantum-based modeling is require… ▽ More The macroscopic behaviors of materials are determined by interactions that occur at multiple lengths and time scales. Depending on the application, describing, predicting, and understanding these behaviors require models that rely on insights from electronic and atomic scales. In such cases, classical simplified approximations at those scales are insufficient, and quantum-based modeling is required. In this paper, we study how quantum effects can modify the mechanical properties of systems relevant to materials engineering. We base our study on a high-fidelity modeling framework that combines two computationally efficient models rooted in quantum first principles: Density Functional Tight Binding (DFTB) and many-body dispersion (MBD). The MBD model is applied to accurately describe non-covalent van der Waals interactions. Through various benchmark applications, we demonstrate the capabilities of this framework and the limitations of simplified modeling. We provide an open-source repository containing all codes, datasets, and examples presented in this work. This repository serves as a practical toolkit that we hope will support the development of future research in effective large-scale and multiscale modeling with quantum-mechanical fidelity. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Enhancing Biomechanical Simulations Based on A Posteriori Error Estimates: The Potential of Dual Weighted Residual-Driven Adaptive Mesh Refinement

Authors: Huu Phuoc Bui, Michel Duprez, Pierre-Yves Rohan, Arnaud Lejeune, Stephane P. A. Bordas, Marek Bucki, Franz Chouly

Abstract: 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 probl… ▽ More 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. △ Less

Submitted 1 March, 2024; originally announced March 2024.

MSC Class: 65N30; 92C10; 74B20

Volume and Mass Conservation in Lagrangian Meshfree Methods

Authors: Pratik Suchde, Christian Leithäuser, Jörg Kuhnert, Stéphane P. A. Bordas

Abstract: Meshfree Lagrangian frameworks for free surface flow simulations do not conserve fluid volume. Meshfree particle methods like SPH are not mimetic, in the sense that discrete mass conservation does not imply discrete volume conservation. On the other hand, meshfree collocation methods typically do not use any notion of mass. As a result, they are neither mass conservative nor volume conservative at… ▽ More Meshfree Lagrangian frameworks for free surface flow simulations do not conserve fluid volume. Meshfree particle methods like SPH are not mimetic, in the sense that discrete mass conservation does not imply discrete volume conservation. On the other hand, meshfree collocation methods typically do not use any notion of mass. As a result, they are neither mass conservative nor volume conservative at the discrete level. In this paper, we give an overview of various sources of conservation errors across different meshfree methods. The present work focuses on one specific issue: unreliable volume and mass definitions. We introduce the concept of representative masses and densities, which are essential for accurate post-processing especially in meshfree collocation methods. Using these, we introduce an artificial compression or expansion in the fluid to rectify errors in volume conservation. Numerical experiments show that the introduced frameworks significantly improve volume conservation behaviour, even for complex industrial test cases such as automotive water crossing. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Breast simulation pipeline: from medical imaging to patient-specific simulations

Authors: Arnaud Mazier, Stéphane P. A. Bordas

Abstract: Breast-conserving surgery is the most acceptable operation for breast cancer removal from an invasive and psychological point of view. Before the surgical procedure, a preoperative MRI is performed in the prone configuration, while the surgery is achieved in the supine position. This leads to a considerable movement of the breast, including the tumor, between the two poses, complicating the surgeo… ▽ More Breast-conserving surgery is the most acceptable operation for breast cancer removal from an invasive and psychological point of view. Before the surgical procedure, a preoperative MRI is performed in the prone configuration, while the surgery is achieved in the supine position. This leads to a considerable movement of the breast, including the tumor, between the two poses, complicating the surgeon's task. In this work, a simulation pipeline allowing the computation of patient-specific geometry and the prediction of personalized breast material properties was put forward. Through image segmentation, a finite element model including the subject-specific geometry is established. By first computing an undeformed state of the breast, the geometrico-material model is calibrated by surface acquisition in the intra-operative stance. Using an elastic corotational formulation, the patient-specific mechanical properties of the breast and skin were identified to obtain the best estimates of the supine configuration. The final results are a Mean Absolute Error of 4.00mm for the mechanical parameters E_breast = 0.32 kPa and E_skin = 22.72 kPa, congruent with the current state-of-the-art. The Covariance Matrix Adaptation Evolution Strategy optimizer converges on average between 5 to 30 min depending on the initial parameters, reaching a simulation speed of 20s. To our knowledge, our model offers one of the best compromises between accuracy and speed. Satisfactory results were obtained for the estimation of breast deformation from preoperative to intra-operative configuration. Furthermore, we have demonstrated the clinical feasibility of such applications using a simulation framework that aims at the smallest disturbance of the actual surgical pipeline. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Level set topology optimization of metamaterial-based heat manipulators using isogeometric analysis

Authors: Chintan Jansari, Stéphane P. A. Bordas, Elena Atroshchenko

Abstract: We exploit level set topology optimization to find the optimal material distribution for metamaterial-based heat manipulators. The level set function, geometry, and solution field are parameterized using the non-uniform rational B-spline (NURBS) basis functions in order to take advantage of easy control of smoothness and continuity. In addition, NURBS approximations can produce conic geometries ex… ▽ More We exploit level set topology optimization to find the optimal material distribution for metamaterial-based heat manipulators. The level set function, geometry, and solution field are parameterized using the non-uniform rational B-spline (NURBS) basis functions in order to take advantage of easy control of smoothness and continuity. In addition, NURBS approximations can produce conic geometries exactly and provide higher efficiency for higher-order elements. The values of the level set function at the control points (called expansion coefficients) are utilized as design variables. For optimization, we use an advanced mathematical programming technique, Sequential Quadratic Programming (SQP). Taking into account a large number of design variables and the small number of constraints associated with our optimization problem, the adjoint method is utilized to calculate the required sensitivities with respect to the design variables. The efficiency and robustness of the proposed method are demonstrated by solving three numerical examples. We have also shown that the current method can handle different geometries and types of objective functions. In addition, regularization techniques such as Tikhonov regularization and volume regularization have been explored to reduce unnecessary complexity and increase the manufacturability of optimized topologies. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Multi-compartment poroelastic models of perfused biological soft tissues: implementation in FEniCSx

Authors: Thomas Lavigne, Stéphane Urcun, Pierre-Yves Rohan, Giuseppe Sciumè, Davide Baroli, Stéphane P. A. Bordas

Abstract: Soft biological tissues demonstrate strong time-dependent and strain-rate mechanical behavior, arising from their intrinsic visco-elasticity and fluid-solid interactions (especially at sufficiently large time scales). The time-dependent mechanical properties of soft tissues influence their physiological functions and are linked to several pathological processes. Poro-elastic modeling represents a… ▽ More Soft biological tissues demonstrate strong time-dependent and strain-rate mechanical behavior, arising from their intrinsic visco-elasticity and fluid-solid interactions (especially at sufficiently large time scales). The time-dependent mechanical properties of soft tissues influence their physiological functions and are linked to several pathological processes. Poro-elastic modeling represents a promising approach because it allows the integration of multiscale/multiphysics data to probe biologically relevant phenomena at a smaller scale and embeds the relevant mechanisms at the larger scale. The implementation of multi-phasic flow poro-elastic models however is a complex undertaking, requiring extensive knowledge. The open-source software FEniCSx Project provides a novel tool for the automated solution of partial differential equations by the finite element method. This paper aims to provide the required tools to model the mixed formulation of poro-elasticity, from the theory to the implementation, within FEniCSx. Several benchmark cases are studied. A column under confined compression conditions is compared to the Terzaghi analytical solution, using the L2-norm. An implementation of poro-hyper-elasticity is proposed. A bi-compartment column is compared to previously published results (Cast3m implementation). For all cases, accurate results are obtained in terms of a normalized Root Mean Square Error (RMSE). Furthermore, the FEniCSx computation is found three times faster than the legacy FEniCS one. The benefits of parallel computation are also highlighted. △ Less

Submitted 1 March, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: https://github.com/Th0masLavigne/Dolfinx_Porous_Media.git

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

Authors: Saurabh Deshpande, Raúl I. Sosa, Stéphane P. A. Bordas, Jakub Lengiewicz

Abstract: 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 bo… ▽ 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

Submitted 24 March, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Journal ref: Front. Mater. 10:1128954

Oncology and mechanics: landmark studies and promising clinical applications

Authors: Stéphane Urcun, Guillermo Lorenzo, Davide Baroli, Pierre-Yves Rohan, Giuseppe Sciumè, Wafa Skalli, Vincent Lubrano, Stéphane P. A. Bordas

Abstract: 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 i… ▽ 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

Submitted 16 November, 2022; originally announced November 2022.

Journal ref: Advances in Applied Mechanics, Volume 55, 2022, Pages 513-571

MAgNET: A Graph U-Net Architecture for Mesh-Based Simulations

Authors: Saurabh Deshpande, Stéphane P. A. Bordas, Jakub Lengiewicz

Abstract: In many cutting-edge applications, high-fidelity computational models prove to be too slow for practical use and are therefore replaced by much faster surrogate models. Recently, deep learning techniques have increasingly been utilized to accelerate such predictions. To enable learning on large-dimensional and complex data, specific neural network architectures have been developed, including convo… ▽ More In many cutting-edge applications, high-fidelity computational models prove to be too slow for practical use and are therefore replaced by much faster surrogate models. Recently, deep learning techniques have increasingly been utilized to accelerate such predictions. To enable learning on large-dimensional and complex data, specific neural network architectures have been developed, including convolutional and graph neural networks. In this work, we present a novel encoder-decoder geometric deep learning framework called MAgNET, which extends the well-known convolutional neural networks to accommodate arbitrary graph-structured data. MAgNET consists of innovative Multichannel Aggregation (MAg) layers and graph pooling/unpooling layers, forming a graph U-Net architecture that is analogous to convolutional U-Nets. We demonstrate the predictive capabilities of MAgNET in surrogate modeling for non-linear finite element simulations in the mechanics of solids. △ Less

Submitted 2 April, 2024; v1 submitted 1 November, 2022; originally announced November 2022.

Journal ref: Engineering Applications of Artificial Intelligence, Volume 133, Part B, 2024, 108055

Smart cloud collocation: geometry-aware adaptivity directly from CAD

Authors: Thibault Jacquemin, Pratik Suchde, Stéphane P. A. Bordas

Abstract: Computer Aided Design (CAD) is widely used in the creation and optimization of various industrial systems and processes. Transforming a CAD geometry into a computational discretization that be used to solve PDEs requires care and a deep knowledge of the selected computational method. In this article, we present a novel integrated collocation scheme based on smart clouds. It allows us to transform… ▽ More Computer Aided Design (CAD) is widely used in the creation and optimization of various industrial systems and processes. Transforming a CAD geometry into a computational discretization that be used to solve PDEs requires care and a deep knowledge of the selected computational method. In this article, we present a novel integrated collocation scheme based on smart clouds. It allows us to transform a CAD geometry into a complete point collocation model, aware of the base geometry, with minimum effort. For this process, only the geometry of the domain, in the form of a STEP file, and the boundary conditions are needed. We also introduce an adaptive refinement process for the resultant smart cloud using an \textit{a posteriori} error indication. The scheme can be applied to any 2D or 3D geometry, to any PDE and can be applied to most point collocation approaches. We illustrate this with the meshfree Generalized Finite Difference (GFD) method applied to steady linear elasticity problems. We further show that each step of this process, from the initial discretization to the refinement strategy, is connected and is affected by the approach selected in the previous step, thus requiring an integrated scheme where the whole solution process should be considered at once. △ Less

Submitted 25 October, 2022; originally announced October 2022.

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

Authors: Arnaud Mazier, Sidaty El Hadramy, Jean-Nicolas Brunet, Jack S. Hale, Stéphane Cotin, Stéphane P. A. Bordas

Abstract: 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. S… ▽ More 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. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Design of metamaterial-based heat manipulators by isogeometric shape optimization

Authors: Chintan Jansari, Stéphane P. A. Bordas, Elena Atroshchenko

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 18 May, 2022; originally announced May 2022.

A graph-based probabilistic geometric deep learning framework with online enforcement of physical constraints to predict the criticality of defects in porous materials

Authors: Vasilis Krokos, Stéphane P. A. Bordas, Pierre Kerfriden

Abstract: Stress prediction in porous materials and structures is challenging due to the high computational cost associated with direct numerical simulations. Convolutional Neural Network (CNN) based architectures have recently been proposed as surrogates to approximate and extrapolate the solution of such multiscale simulations. These methodologies are usually limited to 2D problems due to the high computa… ▽ More Stress prediction in porous materials and structures is challenging due to the high computational cost associated with direct numerical simulations. Convolutional Neural Network (CNN) based architectures have recently been proposed as surrogates to approximate and extrapolate the solution of such multiscale simulations. These methodologies are usually limited to 2D problems due to the high computational cost of 3D voxel based CNNs. We propose a novel geometric learning approach based on a Graph Neural Network (GNN) that efficiently deals with three-dimensional problems by performing convolutions over 2D surfaces only. Following our previous developments using pixel-based CNN, we train the GNN to automatically add local fine-scale stress corrections to an inexpensively computed coarse stress prediction in the porous structure of interest. Our method is Bayesian and generates densities of stress fields, from which credible intervals may be extracted. As a second scientific contribution, we propose to improve the extrapolation ability of our network by deploying a strategy of online physics-based corrections. Specifically, we condition the posterior predictions of our probabilistic predictions to satisfy partial equilibrium at the microscale, at the inference stage. This is done using an Ensemble Kalman algorithm, to ensure tractability of the Bayesian conditioning operation. We show that this innovative methodology allows us to alleviate the effect of undesirable biases observed in the outputs of the uncorrected GNN, and improves the accuracy of the predictions in general. △ Less

Submitted 5 November, 2023; v1 submitted 13 May, 2022; originally announced May 2022.

Comments: 68 pages; 52 figures

A Discrete Droplet Method for Modelling Thin Film Flows

Authors: Anand S Bharadwaj, Joerg Kuhnert, Stephane P. A. Bordas, Pratik Suchde

Abstract: In this paper, we present a new model to simulate the formation, evolution, and break up of a thin film of fluid flowing over a curved surface. Referred to as the discrete droplet method (DDM), the model captures the evolution of thin fluid films by tracking individual moving fluid droplets. In contrast to existing thin film models that solve a PDE to determine the film height, here, we compute th… ▽ More In this paper, we present a new model to simulate the formation, evolution, and break up of a thin film of fluid flowing over a curved surface. Referred to as the discrete droplet method (DDM), the model captures the evolution of thin fluid films by tracking individual moving fluid droplets. In contrast to existing thin film models that solve a PDE to determine the film height, here, we compute the film height by numerical integration based on the aggregation of droplets. The novelty of this approach in using droplets makes it suitable for simulating the formation of fluid films, and modelling thin film flows on partially wetted surfaces. The DDM is a Lagrangian approach, with a force balance on each droplet governing the motion, and derivatives approximated using a smoothed particle hydrodynamics (SPH) like approach. The proposed model is thoroughly validated by comparing results against analytical solutions, against the results of the shallow-water equations for thin film flow, and also against results from a full 3D resolved Navier Stokes model. We also present the use of the DDM on an industrial test case. The results highlight the effectiveness of the model for simulations of flows with thin films. △ Less

Submitted 11 May, 2022; originally announced May 2022.

An a posteriori error estimator for the spectral fractional power of the Laplacian

Authors: Raphaël Bulle, Olga Barrera, Stéphane P. A. Bordas, Franz Chouly, Jack S. Hale

Abstract: 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 a… ▽ More 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. △ Less

Submitted 10 March, 2023; v1 submitted 11 February, 2022; originally announced February 2022.

MSC Class: 65N15; 65N30

Journal ref: Computer Methods in Applied Mechanics and Engineering, Volume 407, 2023, 115943, ISSN 0045-7825

An adapted deflated conjugate gradient solver for robust extended/generalised finite element solutions of large scale, 3D crack propagation problems

Authors: Konstantinos Agathos, Tim Dodwell, Eleni Chatzi, Stephane P. A. Bordas

Abstract: An adapted deflation preconditioner is employed to accelerate the solution of linear systems resulting from the discretization of fracture mechanics problems with well-conditioned extended/generalized finite elements. The deflation space typically used for linear elasticity problems is enriched with additional vectors, accounting for the enrichment functions used, thus effectively removing low fre… ▽ More An adapted deflation preconditioner is employed to accelerate the solution of linear systems resulting from the discretization of fracture mechanics problems with well-conditioned extended/generalized finite elements. The deflation space typically used for linear elasticity problems is enriched with additional vectors, accounting for the enrichment functions used, thus effectively removing low frequency components of the error. To further improve performance, deflation is combined, in a multiplicative way, with a block-Jacobi preconditioner, which removes high frequency components of the error as well as linear dependencies introduced by enrichment. The resulting scheme is tested on a series of non-planar crack propagation problems and compared to alternative linear solvers in terms of performance. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Probabilistic Deep Learning for Real-Time Large Deformation Simulations

Authors: Saurabh Deshpande, Jakub Lengiewicz, Stéphane P. A. Bordas

Abstract: 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 respo… ▽ 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

Submitted 4 July, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Journal ref: Computer Methods in Applied Mechanics and Engineering, 2022, Volume 398

Neural-network acceleration of projection-based model-order-reduction for finite plasticity: Application to RVEs

Authors: S. Vijayaraghavan, L. Wu, L. Noels, S. P. A. Bordas, S. Natarajan, L. A. A. Beex

Abstract: Compared to conventional projection-based model-order-reduction, its neural-network acceleration has the advantage that the online simulations are equation-free, meaning that no system of equations needs to be solved iteratively. Consequently, no stiffness matrix needs to be constructed and the stress update needs to be computed only once per increment. In this contribution, a recurrent neural net… ▽ More Compared to conventional projection-based model-order-reduction, its neural-network acceleration has the advantage that the online simulations are equation-free, meaning that no system of equations needs to be solved iteratively. Consequently, no stiffness matrix needs to be constructed and the stress update needs to be computed only once per increment. In this contribution, a recurrent neural network is developed to accelerate a projection-based model-order-reduction of the elastoplastic mechanical behaviour of an RVE. In contrast to a neural network that merely emulates the relation between the macroscopic deformation (path) and the macroscopic stress, the neural network acceleration of projection-based model-order-reduction preserves all microstructural information, at the price of computing this information once per increment. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Machine learning in the social and health sciences

Authors: Anja K. Leist, Matthias Klee, Jung Hyun Kim, David H. Rehkopf, Stéphane P. A. Bordas, Graciela Muniz-Terrera, Sara Wade

Abstract: The uptake of machine learning (ML) approaches in the social and health sciences has been rather slow, and research using ML for social and health research questions remains fragmented. This may be due to the separate development of research in the computational/data versus social and health sciences as well as a lack of accessible overviews and adequate training in ML techniques for non data scie… ▽ More The uptake of machine learning (ML) approaches in the social and health sciences has been rather slow, and research using ML for social and health research questions remains fragmented. This may be due to the separate development of research in the computational/data versus social and health sciences as well as a lack of accessible overviews and adequate training in ML techniques for non data science researchers. This paper provides a meta-mapping of research questions in the social and health sciences to appropriate ML approaches, by incorporating the necessary requirements to statistical analysis in these disciplines. We map the established classification into description, prediction, and causal inference to common research goals, such as estimating prevalence of adverse health or social outcomes, predicting the risk of an event, and identifying risk factors or causes of adverse outcomes. This meta-mapping aims at overcoming disciplinary barriers and starting a fluid dialogue between researchers from the social and health sciences and methodologically trained researchers. Such mapping may also help to fully exploit the benefits of ML while considering domain-specific aspects relevant to the social and health sciences, and hopefully contribute to the acceleration of the uptake of ML applications to advance both basic and applied social and health sciences research. △ Less

Submitted 20 June, 2021; originally announced June 2021.

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

Authors: Paul Hauseux, Alberto Ambrosetti, Stéphane P. A. Bordas, Alexandre Tkatchenko

Abstract: 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 interact… ▽ More 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. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Circumferential Crack Modeling of Thin Cylindrical Shells in Modal Deformation

Authors: Ali Alijani, Olga Barrera, Stephane P. A. Bordas

Abstract: An innovative technique, called conversion, is introduced to model circumferential cracks in thin cylindrical shells. The semi-analytical finite element method is applied to investigate the modal deformation of the cylinder. An element including the crack is divided into three sub-elements with four nodes in which the stiffness matrix is enriched. The crack characteristics are included in the fini… ▽ More An innovative technique, called conversion, is introduced to model circumferential cracks in thin cylindrical shells. The semi-analytical finite element method is applied to investigate the modal deformation of the cylinder. An element including the crack is divided into three sub-elements with four nodes in which the stiffness matrix is enriched. The crack characteristics are included in the finite element method relations through conversion matrices and a rotational spring corresponding to the crack. Conversion matrices obtained by applying continuity conditions at the crack tip are used to transform displacements of the middle nodes to those of the main nodes. Moreover, another technique, called spring set, is represented based on a set of springs to model the crack as a separated element. Components of the stiffness matrix related to the separated element are incorporated while the geometric boundary conditions at the crack tip are satisfied. The effects of the circumferential mode number, the crack depth and the length of the cylinder on the critical buckling load are investigated. Experimental tests, ABAQUS modeling and results from literature are used to verify and validate the results and derived relations. In addition, the crack effect on the natural frequency is examined using the vibration analysis based on the conversion technique. △ Less

Submitted 15 May, 2021; originally announced May 2021.

Comments: 23 pages, 12 figures, and 5 tables

A rigged model of the breast for preoperative surgical planning

Authors: Arnaud Mazier, Sophie Ribes, Benjamin Gilles, Stéphane P. A Bordas, .

Abstract: In breast surgical practice, drawing is part of the preoperative planning procedure and is essential for a successful operation. In this study, we design a pipeline to assist surgeons with patient-specific breast surgical drawings. We use a deformable torso model containing the surgical patterns to match any breast surface scan. To be compatible with surgical timing, we build an articulated model… ▽ More In breast surgical practice, drawing is part of the preoperative planning procedure and is essential for a successful operation. In this study, we design a pipeline to assist surgeons with patient-specific breast surgical drawings. We use a deformable torso model containing the surgical patterns to match any breast surface scan. To be compatible with surgical timing, we build an articulated model through a skinning process coupled with shape deformers to enhance a fast registration process. On one hand, the scalable bones of the skinning account for pose and morphological variations of the patients. On the other hand, pre-designed artistic blendshapes create a linear space for guaranteeing anatomical variations. Then, we apply meaningful constraints to the model to find a trade-off between precision and speed. The experiments were conducted on 7 patients, in 2 different poses (prone and supine) with a breast size ranging from 36A and 42C (US/UK bra sizing). The acquisitions were obtained using the depth camera Structure Sensor, and the breast scans were acquired in less than 1 minute. The result is a registration method converging within a few seconds (3 maximum), reaching a Mean Absolute Error of 2.3 mm for mesh registration and 8.0 mm for breast anatomical landmarks. Compared to the existing literature, our model can be personalized and does not require any database. Finally, our registered model can be used to transfer surgical reference patterns onto any patient in any position. △ Less

Submitted 31 May, 2022; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 20 pages, 8 figures, submitted to Journal of Biomechanics

Journal ref: Journal of Biomechanics, Volume 128, 2021, 110645, ISSN 0021-9290,

On the finite element analysis of functionally graded sandwich curved beams via a new refined higher shear deformation theory

Authors: Mohamed-Ouejdi Belarbi, Mohammed Sid Ahmed Houari, Hicham Hirane1, Ahmed Amine Daikh, Stéphane Pierre Alain Bordas

Abstract: In the present paper, a new parabolic shear deformation beam theory is developed and applied to investigate the bending behavior of functionally graded (FG) sandwich curved beam. The present theory is exploited to satisfy parabolic variation of shear stress distribution along the thickness direction thereby obviating the use of any shear correction factors. The material properties of FG sandwich b… ▽ More In the present paper, a new parabolic shear deformation beam theory is developed and applied to investigate the bending behavior of functionally graded (FG) sandwich curved beam. The present theory is exploited to satisfy parabolic variation of shear stress distribution along the thickness direction thereby obviating the use of any shear correction factors. The material properties of FG sandwich beam change continuously from one surface to another according to a power-law function. Three common configurations of FG beams are used for the study, namely: (a) single layer FG beam; (b) sandwich beam with FG face sheets and homogeneous core and (c) sandwich beams with homogeneous face sheets and FG core. The governing equations derived herein are solved by employing the finite element method using a two-noded beam element, developed for this purpose. The robustness and reliability of the developed finite element model are demonstrated by comparing its results with those available by other researchers in existing literature. The comparison studies show that the proposed model is: (a) accurate and comparable with the literature; b) of fast rate of convergence to the reference solution; c) excellent in terms of numerical stability and d) valid for FG sandwich curved beams. Moreover, comprehensive numerical results are presented and discussed in detail to investigate the effects of volume fraction index, radius of curvature, material distributions, length-to-thickness ratio, face-to-core-thickness ratio, loadings and boundary conditions on the static response of FG curved sandwich beam. New referential results are reported which will be serve as a benchmark for future research. △ Less

Submitted 29 April, 2021; originally announced April 2021.

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

Authors: Arnaud Mazier, Alexandre Bilger, Antonio E. Forte, Igor Peterlik, Jack S. Hale, Stéphane P. A. Bordas, ., Institute of Computational Engineering, Department of Engineering, University of Luxembourg, Esch-sur-Alzette, Luxembourg., Harvard University, Cambridge, USA., Department of Electronics, Information, Bioengineering, Politecnico di Milano, Milan, Italy., Institute of Computer Science, Masaryk University, Czech Republic., Institute of Research , et al. (3 additional authors not shown)

Abstract: 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 config… ▽ More 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. △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 33 pages, 12 figures, submitted to Elsevier

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

Authors: Raphaël Bulle, Jack S. Hale, Alexei Lozinski, Stéphane P. A. Bordas, Franz Chouly

Abstract: 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 computationa… ▽ 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

Submitted 7 February, 2022; v1 submitted 8 February, 2021; originally announced February 2021.

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

Authors: Vasilis Krokos, Viet Bui Xuan, Stéphane P. A. Bordas, Philippe Young, Pierre Kerfriden

Abstract: 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… ▽ More 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. △ Less

Submitted 19 January, 2022; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 50 pages, 42 figures; corrected typos; clarifications

Taylor-series expansion based numerical methods: a primer, performance benchmarking and new approaches for problems with non-smooth solutions

Authors: Thibault Jacquemin, Satyendra Tomar, Konstantinos Agathos, Shoya Mohseni-Mofidi, Stéphane P. A. Bordas

Abstract: We provide a primer to numerical methods based on Taylor series expansions such as generalized finite difference methods and collocation methods. We provide a detailed benchmarking strategy for these methods as well as all data files including input files, boundary conditions, point distribution and solution fields, so as to facilitate future benchmarking of new methods. We review traditional meth… ▽ More We provide a primer to numerical methods based on Taylor series expansions such as generalized finite difference methods and collocation methods. We provide a detailed benchmarking strategy for these methods as well as all data files including input files, boundary conditions, point distribution and solution fields, so as to facilitate future benchmarking of new methods. We review traditional methods and recent ones which appeared in the last decade. We aim to help newcomers to the field understand the main characteristics of these methods and to provide sufficient information to both simplify implementation and benchmarking of new methods. Some of the examples are chosen within a subset of problems where collocation is traditionally known to perform sub-par, namely when the solution sought is non-smooth, i.e. contains discontinuities, singularities or sharp gradients. For such problems and other simpler ones with smooth solutions, we study in depth the influence of the weight function, correction function, and the number of nodes in a given support. We also propose new stabilization approaches to improve the accuracy of the numerical methods. In particular, we experiment with the use of a Voronoi diagram for weight computation, collocation method stabilization approaches, and support node selection for problems with singular solutions. With an appropriate selection of the above-mentioned parameters, the resulting collocation methods are compared to the moving least-squares method (and variations thereof), the radial basis function finite difference method and the finite element method. Extensive tests involving two and three dimensional problems indicate that the methods perform well in terms of efficiency (accuracy versus computational time), even for non-smooth solutions. △ Less

Submitted 24 September, 2020; originally announced September 2020.

Comments: Arch Computat Methods Eng (2019)

Distributed prediction of unsafe reconfiguration scenarios of modular robotic Programmable Matter

Authors: Benoît Piranda, Paweł Chodkiewicz, Paweł Hołobut, Stéphane P. A. Bordas, Julien Bourgeois, Jakub Lengiewicz

Abstract: We present a distributed framework for predicting whether a planned reconfiguration step of a modular robot will mechanically overload the structure, causing it to break or lose stability under its own weight. The algorithm is executed by the modular robot itself and based on a distributed iterative solution of mechanical equilibrium equations derived from a simplified model of the robot. The mode… ▽ More We present a distributed framework for predicting whether a planned reconfiguration step of a modular robot will mechanically overload the structure, causing it to break or lose stability under its own weight. The algorithm is executed by the modular robot itself and based on a distributed iterative solution of mechanical equilibrium equations derived from a simplified model of the robot. The model treats inter-modular connections as beams and assumes no-sliding contact between the modules and the ground. We also provide a procedure for simplified instability detection. The algorithm is verified in the Programmable Matter simulator VisibleSim, and in real-life experiments on the modular robotic system Blinky Blocks. △ Less

Submitted 3 May, 2021; v1 submitted 19 June, 2020; originally announced June 2020.

Journal ref: IEEE Transactions on Robotics, vol. 37, no. 6, pp. 2226-2233, Dec. 2021

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

Authors: Michel Duprez, Stéphane P. A. Bordas, Marek Bucki, Huu Phuoc Bui, Franz Chouly, Vanessa Lleras, Claudio Lobos, Alexei Lozinski, Pierre-Yves Rohan, Satyendra Tomar

Abstract: 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… ▽ More 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. △ Less

Submitted 18 June, 2018; originally announced June 2018.

Comments: 35 pages, 10 figures

Journal ref: Applied Mathematical Modelling Volume 77, Part 1, January 2020, Pages 709-723

A volume-averaged nodal projection method for the Reissner-Mindlin plate model

Authors: Alejandro Ortiz-Bernardin, Philip Köbrich, Jack S. Hale, Edgardo Olate-Sanzana, Stéphane P. A. Bordas, Sundararajan Natarajan

Abstract: We introduce a novel meshfree Galerkin method for the solution of Reissner-Mindlin plate problems that is written in terms of the primitive variables only (i.e., rotations and transverse displacement) and is devoid of shear-locking. The proposed approach uses linear maximum-entropy approximations and is built variationally on a two-field potential energy functional wherein the shear strain, writte… ▽ More We introduce a novel meshfree Galerkin method for the solution of Reissner-Mindlin plate problems that is written in terms of the primitive variables only (i.e., rotations and transverse displacement) and is devoid of shear-locking. The proposed approach uses linear maximum-entropy approximations and is built variationally on a two-field potential energy functional wherein the shear strain, written in terms of the primitive variables, is computed via a volume-averaged nodal projection operator that is constructed from the Kirchhoff constraint of the three-field mixed weak form. The stability of the method is rendered by adding bubble-like enrichment to the rotation degrees of freedom. Some benchmark problems are presented to demonstrate the accuracy and performance of the proposed method for a wide range of plate thicknesses. △ Less

Submitted 29 July, 2018; v1 submitted 8 March, 2018; originally announced March 2018.

Corotational Cut Finite Element Method for real-time surgical simulation: application to needle insertion simulation

Authors: Huu Phuoc Bui, Satyendra Tomar, Stéphane P. A. Bordas

Abstract: This paper describes the use of the corotational cut Finite Element Method (FEM) for real-time surgical simulation. Users only need to provide a background mesh which is not necessarily conforming to the boundaries/interfaces of the simulated object. The details of the surface, which can be directly obtained from binary images, are taken into account by a multilevel embedding algorithm applied to… ▽ More This paper describes the use of the corotational cut Finite Element Method (FEM) for real-time surgical simulation. Users only need to provide a background mesh which is not necessarily conforming to the boundaries/interfaces of the simulated object. The details of the surface, which can be directly obtained from binary images, are taken into account by a multilevel embedding algorithm applied to elements of the background mesh that cut by the surface. Boundary conditions can be implicitly imposed on the surface using Lagrange multipliers. The implementation is verified by convergence studies with optimal rates. The algorithm is applied to various needle insertion simulations (e.g. for biopsy or brachytherapy) into brain and liver to verify the reliability of method, and numerical results show that the present method can make the discretisation independent from geometric description, and can avoid the complexity of mesh generation of complex geometries while retaining the accuracy of the standard FEM. Using the proposed approach is very suitable for real-time and patient specific simulations as it improves the simulation accuracy by taking into account automatically and properly the simulated geometry. △ Less

Submitted 8 December, 2017; originally announced December 2017.

Comments: 25 pages, 25 figures

Skew-symmetric Nitsche's formulation in isogeometric analysis: Dirichlet and symmetry conditions, patch coupling and frictionless contact

Authors: Qingyuan Hu, Franz Chouly, Ping Hu, Gengdong Cheng, Stéphane Pierre Alain Bordas

Abstract: 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 formul… ▽ 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

Submitted 27 April, 2018; v1 submitted 28 November, 2017; originally announced November 2017.

Isogeometric analysis of thin Reissner-Mindlin plates and shells: locking phenomena and B-bar method

Authors: Qingyuan Hu, Yang Xia, Sundararajan Natarajan, Andreas Zilian, Ping Hu, Stéphane P. A. Bordas

Abstract: We propose a local type of B-bar formulation, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation allows the flexible utilization of basis functions of different order as the projec… ▽ More We propose a local type of B-bar formulation, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation allows the flexible utilization of basis functions of different order as the projection bases. The present formulation is much cheaper computationally than the classical $\bar{B}$ method. We show the numerical consistency of the scheme through numerical examples, moreover they show that the proposed formulation alleviates locking and yields good accuracy even for slenderness ratios of $1 \times 10^5$, and has the ability to capture deformations of thin shells using relatively coarse meshes. In addition it can be opined that the proposed method is less sensitive to locking and mesh distortion. △ Less

Submitted 22 February, 2018; v1 submitted 1 September, 2017; originally announced September 2017.

Comments: 32 pages, 21 figures

Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization

Authors: Gang Xu, Ming Li, Bernard Mourrain, Timon Rabczuk, Jinlan Xu, Stephane P. A. Bordas

Abstract: In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including Bézier extraction an… ▽ More In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including Bézier extraction and subdivision are performed on each boundary curve in order to generate a high-quality planar parameterization; then a robust planar domain partition framework is proposed to construct high-quality patch-meshing results with few singularities from the discrete boundary formed by connecting the end points of the resulting boundary segments. After the topology information generation of quadrilateral decomposition, the optimal placement of interior Bézier curves corresponding to the interior edges of the quadrangulation is constructed by a global optimization method to achieve a patch-partition with high quality. Finally, after the imposition of C1=G1-continuity constraints on the interface of neighboring Bézier patches with respect to each quad in the quadrangulation, the high-quality Bézier patch parameterization is obtained by a C1-constrained local optimization method to achieve uniform and orthogonal iso-parametric structures while keeping the continuity conditions between patches. The efficiency and robustness of the proposed method are demonstrated by several examples which are compared to results obtained by the skeleton-based parameterization approach. △ Less

Submitted 8 July, 2017; v1 submitted 3 July, 2017; originally announced July 2017.

Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximaTion (GIFT)

Authors: Elena Atroshchenko, Gang Xu, Satyendra Tomar, Stephane P. A. Bordas

Abstract: This paper presents an approach to generalize the concept of isogeometric analysis (IGA) by allowing different spaces for parameterization of the computational domain and for approximation of the solution field. The method inherits the main advantage of isogeometric analysis, i.e. preserves the original, exact CAD geometry (for example, given by NURBS), but allows pairing it with an approximation… ▽ More This paper presents an approach to generalize the concept of isogeometric analysis (IGA) by allowing different spaces for parameterization of the computational domain and for approximation of the solution field. The method inherits the main advantage of isogeometric analysis, i.e. preserves the original, exact CAD geometry (for example, given by NURBS), but allows pairing it with an approximation space which is more suitable/flexible for analysis, for example, T-splines, LR-splines, (truncated) hierarchical B-splines, and PHT-splines. This generalization offers the advantage of adaptive local refinement without the need to re-parameterize the domain, and therefore without weakening the link with the CAD model. We demonstrate the use of the method with different choices of the geometry and field splines, and show that, despite the failure of the standard patch test, the optimum convergence rate is achieved for non-nested spaces. △ Less

Submitted 9 July, 2017; v1 submitted 20 June, 2017; originally announced June 2017.

Comments: 39 pages, 24 figures, 10 tables, 2 algorithms

Controlling the Error on Target Motion through Real-time Mesh Adaptation: Applications to Deep Brain Stimulation

Authors: Huu Phuoc Bui, Satyendra Tomar, Hadrien Courtecuisse, Michel Audette, Stéphane Cotin, Stéphane P. A. Bordas

Abstract: We present an error-controlled mesh refinement procedure for needle insertion simulation and apply it to the simulation of electrode implantation for deep brain stimulation, including brain shift. Our approach enables to control the error in the computation of the displacement and stress fields around the needle tip and needle shaft by suitably refining the mesh, whilst maintaining a coarser mesh… ▽ More We present an error-controlled mesh refinement procedure for needle insertion simulation and apply it to the simulation of electrode implantation for deep brain stimulation, including brain shift. Our approach enables to control the error in the computation of the displacement and stress fields around the needle tip and needle shaft by suitably refining the mesh, whilst maintaining a coarser mesh in other parts of the domain. We demonstrate through academic and practical examples that our approach increases the accuracy of the displacement and stress fields around the needle without increasing the computational expense. This enables real-time simulations. The proposed methodology has direct implications to increase the accuracy and control the computational expense of the simulation of percutaneous procedures such as biopsy, brachytherapy, regional anesthesia, or cryotherapy and can be essential to the development of robotic guidance. △ Less

Submitted 30 September, 2017; v1 submitted 25 April, 2017; originally announced April 2017.

Comments: 21 pages, 14 figures

Isogeometric finite element analysis of time-harmonic exterior acoustic scattering problems

Authors: Tahsin Khajah, Xavier Antoine, Stéphane P. A. Bordas

Abstract: We present an isogeometric analysis of time-harmonic exterior acoustic problems. The infinite space is truncated by a fictitious boundary and (simple) absorbing boundary conditions are applied. The truncation error is included in the exact solution so that the reported error is an indicator of the performance of the isogeometric analysis, in particular of the related pollution error. Numerical res… ▽ More We present an isogeometric analysis of time-harmonic exterior acoustic problems. The infinite space is truncated by a fictitious boundary and (simple) absorbing boundary conditions are applied. The truncation error is included in the exact solution so that the reported error is an indicator of the performance of the isogeometric analysis, in particular of the related pollution error. Numerical results performed with high-order basis functions (third or fourth orders) showed no visible pollution error even for very high frequencies. This property combined with exact geometrical representation makes isogeometric analysis a very promising platform to solve high-frequency acoustic problems. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Bayesian inference for the stochastic identification of elastoplastic material parameters: Introduction, misconceptions and insights

Authors: Hussein Rappel, Lars A. A. Beex, Jack S. Hale, Stephane P. A. Bordas

Abstract: We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material parameters. For this purpose a single spring is considered, for which the stress-strain curves are artificially created. Besides offering a didactic introduction to BI, this paper proposes an approach… ▽ More We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material parameters. For this purpose a single spring is considered, for which the stress-strain curves are artificially created. Besides offering a didactic introduction to BI, this paper proposes an approach to incorporate statistical errors both in the measured stresses, and in the measured strains. It is assumed that the uncertainty is only due to measurement errors and the material is homogeneous. Furthermore, a number of possible misconceptions on BI are highlighted based on the purely elastic case. △ Less

Submitted 15 January, 2017; v1 submitted 8 June, 2016; originally announced June 2016.

Isogeometric analysis for functionally graded microplates based on modified couple stress theory

Authors: Hoang X. Nguyen, Tuan N. Nguyen, M. Abdel-Wahab, S. P. A. Bordas, H. Nguyen-Xuan, Thuc P. Vo

Abstract: Analysis of static bending, free vibration and buckling behaviours of functionally graded microplates is investigated in this study. The main idea is to use the isogeometric analysis in associated with novel four-variable refined plate theory and quasi-3D theory. More importantly, the modified couple stress theory with only one material length scale parameter is employed to effectively capture the… ▽ More Analysis of static bending, free vibration and buckling behaviours of functionally graded microplates is investigated in this study. The main idea is to use the isogeometric analysis in associated with novel four-variable refined plate theory and quasi-3D theory. More importantly, the modified couple stress theory with only one material length scale parameter is employed to effectively capture the size-dependent effects within the microplates. Meanwhile, the quasi-3D theory which is constructed from a novel seventh-order shear deformation refined plate theory with four unknowns is able to consider both shear deformations and thickness stretching effect without requiring shear correction factors. The NURBS-based isogeometric analysis is integrated to exactly describe the geometry and approximately calculate the unknown fields with higher-order derivative and continuity requirements. The convergence and verification show the validity and efficiency of this proposed computational approach in comparison with those existing in the literature. It is further applied to study the static bending, free vibration and buckling responses of rectangular and circular functionally graded microplates with various types of boundary conditions. A number of investigations are also conducted to illustrate the effects of the material length scale, material index, and length-to-thickness ratios on the responses of the microplates. △ Less

Submitted 2 April, 2016; originally announced April 2016.

Comments: 57 pages, 14 figures, 18 tables

On the equivalence between the cell-based smoothed finite element method and the virtual element method

Authors: Sundararajan Natarajan, Stéphane P. A. Bordas, Ean Tat Ooi

Abstract: We revisit the cell-based smoothed finite element method (SFEM) for quadrilateral elements and extend it to arbitrary polygons and polyhedrons in 2D and 3D, respectively. We highlight the similarity between the SFEM and the virtual element method (VEM). Based on the VEM, we propose a new stabilization approach to the SFEM when applied to arbitrary polygons and polyhedrons. The accuracy and the con… ▽ More We revisit the cell-based smoothed finite element method (SFEM) for quadrilateral elements and extend it to arbitrary polygons and polyhedrons in 2D and 3D, respectively. We highlight the similarity between the SFEM and the virtual element method (VEM). Based on the VEM, we propose a new stabilization approach to the SFEM when applied to arbitrary polygons and polyhedrons. The accuracy and the convergence properties of the SFEM are studied with a few benchmark problems in 2D and 3D linear elasticity. Later, the SFEM is combined with the scaled boundary finite element method to problems involving singularity within the framework of the linear elastic fracture mechanics in 2D. △ Less

Submitted 7 October, 2014; v1 submitted 7 July, 2014; originally announced July 2014.

Isogeometric finite element analysis of functionally graded plates using a refined plate theory

Authors: H. Nguyen-Xuan, Loc V. Tran, Chien H. Thai, S. Kulasegaram, S. P. A. Bordas

Abstract: We propose in this paper a novel inverse tangent transverse shear deformation formulation for functionally graded material (FGM) plates. The isogeometric finite element analysis (IGA) of static, free vibration and buckling problems of FGM plates is then addressed using a refined plate theory (RPT). The RPT enables us to describe the non-linear distribution of shear stresses through the plate thick… ▽ More We propose in this paper a novel inverse tangent transverse shear deformation formulation for functionally graded material (FGM) plates. The isogeometric finite element analysis (IGA) of static, free vibration and buckling problems of FGM plates is then addressed using a refined plate theory (RPT). The RPT enables us to describe the non-linear distribution of shear stresses through the plate thickness without any requirement of shear correction factors (SCF). IGA utilizes basis functions, namely B-splines or non-uniform rational B-splines (NURBS), which achieve easily the smoothness of any arbitrary order. It hence satisfies the C1 requirement of the RPT model. The present method approximates the displacement field of four degrees of freedom per each control point and retains the computational efficiency while ensuring the high accuracy in solution. △ Less

Submitted 7 October, 2013; originally announced October 2013.

Nitsche's method method for mixed dimensional analysis: conforming and non-conforming continuum-beam and continuum-plate coupling

Authors: Vinh Phu Nguyena, Pierre Kerfriden, Susanne Clausb Stephane P. A. Bordas

Abstract: A Nitche's method is presented to couple different mechanical models. They include coupling of a solid and a beam and of a solid and a plate. Both conforming and non-conforming formulations are presented. In a non-conforming for- mulation, the structure domain is overlapped by a refined solid model. Applications can be found in multi-dimensional analyses in which parts of a structure are modeled w… ▽ More A Nitche's method is presented to couple different mechanical models. They include coupling of a solid and a beam and of a solid and a plate. Both conforming and non-conforming formulations are presented. In a non-conforming for- mulation, the structure domain is overlapped by a refined solid model. Applications can be found in multi-dimensional analyses in which parts of a structure are modeled with solid elements and others are modeled using a coarser model with beam and/or plate elements. Discretisations are performed using both standard Lagrange elements and high order NURBS (Non Uniform Rational Bsplines) based isogeometric elements. We present various examples to demonstrate the performance of the method. △ Less

Submitted 13 August, 2013; originally announced August 2013.

Comments: 41 pages, 34 figures

An extended finite element method with smooth nodal stress

Authors: X. Peng, S. Kulasegaram, S. P. A. Bordas, S. C. Wu

Abstract: The enrichment formulation of double-interpolation finite element method (DFEM) is developed in this paper. DFEM is first proposed by Zheng \emph{et al} (2011) and it requires two stages of interpolation to construct the trial function. The first stage of interpolation is the same as the standard finite element interpolation. Then the interpolation is reproduced by an additional procedure using th… ▽ More The enrichment formulation of double-interpolation finite element method (DFEM) is developed in this paper. DFEM is first proposed by Zheng \emph{et al} (2011) and it requires two stages of interpolation to construct the trial function. The first stage of interpolation is the same as the standard finite element interpolation. Then the interpolation is reproduced by an additional procedure using the nodal values and nodal gradients which are derived from the first stage as interpolants. The re-constructed trial functions are now able to produce continuous nodal gradients, smooth nodal stress without post-processing and higher order basis without increasing the total degrees of freedom. Several benchmark numerical examples are performed to investigate accuracy and efficiency of DFEM and enriched DFEM. When compared with standard FEM, super-convergence rate and better accuracy are obtained by DFEM. For the numerical simulation of crack propagation, better accuracy is obtained in the evaluation of displacement norm, energy norm and the stress intensity factor. △ Less

Submitted 3 June, 2013; originally announced June 2013.

Space-time goal-oriented reduced basis approximation for linear wave equation

Authors: Khac Chi Hoang, Pierre Kerfriden, Stephane P. A. Bordas

Abstract: In this paper, we study numerically the linear damped second-order hyperbolic partial differential equation (PDE) with affine parameter dependence using a goal-oriented approach by finite element (FE) and reduced basis (RB) methods. The main contribution of this paper is the "goal-oriented" proper orthogonal decomposition (POD)-Greedy sampling procedure within the RB approximation context. First,… ▽ More In this paper, we study numerically the linear damped second-order hyperbolic partial differential equation (PDE) with affine parameter dependence using a goal-oriented approach by finite element (FE) and reduced basis (RB) methods. The main contribution of this paper is the "goal-oriented" proper orthogonal decomposition (POD)-Greedy sampling procedure within the RB approximation context. First, we introduce the RB recipe: Galerkin projection onto a space $Y_N$ spanned by solutions of the governing PDE at $N$ selected points in parameter space. This set of $N$ parameter points is constructed by the standard POD-Greedy sampling procedure already developed. Second, based on the affine parameter dependence, we make use of the offline-online computational procedures: in the offline stage, we generate the RB space; in the online stage, given a new parameter value, we calculate rapidly and accurately the space-time RB output of interest and its associated asymptotic error. The proposed goal-oriented POD-Greedy sampling procedure can now be implemented and will look for the parameter points such that it minimizes this (asymptotic) output error rather than the solution error (or, error indicator which is the dual norm of residual) as in the standard POD-Greedy procedure. Numerical results show that the new goal-oriented POD-Greedy sampling procedure improves significantly the accuracy of the space-time output computation in comparison with the standard POD-Greedy one. The method is thus ideally suited for repeated, rapid and reliable evaluation of input-output relationships within the space-time setting. △ Less

Submitted 14 September, 2013; v1 submitted 15 May, 2013; originally announced May 2013.