Browse ORBi

- What it is and what it isn't
- Green Road / Gold Road?
- Ready to Publish. Now What?
- How can I support the OA movement?
- Where can I learn more?

ORBi

Taylor-Series Expansion Based Numerical Methods: A Primer, Performance Benchmarking and New Approaches for Problems with Non-smooth Solutions Jacquemin, Thibault Augustin Marie ; Tomar, Satyendra ; Agathos, Konstantinos et al in Archives of Computational Methods in Engineering (2019) 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 ... [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 ▲] Detailed reference viewed: 79 (11 UL)Weak and strong from meshless methods for linear elastic problem under fretting contact conditions ; ; et al in Tribology International (2019), 138 We present numerical computation of stresses under fretting fatigue conditions derived from closed form expressions. The Navier-Cauchy equations, that govern the problem, are solved with strong and weak ... [more ▼] We present numerical computation of stresses under fretting fatigue conditions derived from closed form expressions. The Navier-Cauchy equations, that govern the problem, are solved with strong and weak form meshless numerical methods. The results are compared to the solution obtained from well-established commercial package ABAQUS, which is based on finite element method (FEM). The results show that the weak form meshless solution exhibits similar behavior as the FEM solution, while, in this particular case, strong form meshless solution performs better in capturing the peak in the surface stress. This is of particular interest in fretting fatigue, since it directly influences crack initiation. The results are presented in terms of von Mises stress contour plots, surface stress profiles, and the convergence plots for all three methods involved in the study. [less ▲] Detailed reference viewed: 61 (1 UL)ADVANCES IN GEOMETRY INDEPENDENT APPROXIMATIONS ; ; Bordas, Stéphane et al Scientific Conference (2019, April 11) We present recent advances in geometry independent field approximations. The GIFT approach is a generalisation of isogeometric analysis where the approximation used to describe the field variables no ... [more ▼] We present recent advances in geometry independent field approximations. The GIFT approach is a generalisation of isogeometric analysis where the approximation used to describe the field variables no-longer has to be identical to the approximation used to describe the geometry of the domain. As such, the geometry can be described using usual CAD representations, e.g. NURBS, which are the most common in the CAD area, whilst local refinement and meshes approximations can be used to describe the field variables, enabling local adaptivity. We show in which cases the approach passes the patch test and present applications to various mechanics, fracture and multi-physics problems. [less ▲] Detailed reference viewed: 316 (24 UL)Corotational cut finite element method for real-time surgical simulation: Application to needle insertion simulation Bui, Huu Phuoc ; Tomar, Satyendra ; Bordas, Stéphane in Computer Methods in Applied Mechanics and Engineering (2018), 345 We present the corotational cut Finite Element Method (FEM) for real-time surgical simulation. The only requirement of the proposed method is a background mesh, which is not necessarily conforming to the ... [more ▼] We present the corotational cut Finite Element Method (FEM) for real-time surgical simulation. The only requirement of the proposed method is 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 which is applied to elements of the background mesh that are cut by the surface. Dirichlet boundary conditions can be implicitly imposed on the surface using Lagrange multipliers, whereas traction or Neumann boundary conditions, which is/are applied on parts of the surface, can be distributed to the background nodes using shape functions. The implementation is verified by convergences studies, of the geometry and of numerical solutions, which exhibit optimal rates. To verify the reliability of the method, it is applied to various needle insertion simulations (e.g. for biopsy or brachytherapy) into brain and liver models. The numerical results show that, while retaining the accuracy of the standard FEM, the proposed method can (1) make the discretisation independent from geometric description, (2) avoid the complexity of mesh generation for complex geometries, and (3) provide computational speed suitable for real-time simulations. Thereby, the proposed method is very suitable for patient-specific simulations as it improves the simulation accuracy by automatically, and properly, taking the simulated geometry into account, while keeping the low computational cost. [less ▲] Detailed reference viewed: 52 (0 UL)Adaptive Isogeometric analysis for plate vibrations: An efficient approach of local refinement based on hierarchical a posteriori error estimation ; ; Tomar, Satyendra et al in Computer Methods in Applied Mechanics and Engineering (2018), 342 This paper presents a novel methodology of local adaptivity for the frequency-domain analysis of the vibrations of Reissner–Mindlin plates. The adaptive discretization is based on the recently developed ... [more ▼] This paper presents a novel methodology of local adaptivity for the frequency-domain analysis of the vibrations of Reissner–Mindlin plates. The adaptive discretization is based on the recently developed Geometry Independent Field approximaTion (GIFT) framework, which may be seen as a generalization of the Iso-Geometric Analysis (IGA).Within the GIFT framework, we describe the geometry of the structure exactly with NURBS (Non-Uniform Rational B-Splines), whilst independently employing Polynomial splines over Hierarchical T-meshes (PHT)-splines to represent the solution field. The proposed strategy of local adaptivity, wherein a posteriori error estimators are computed based on inexpensive hierarchical h-refinement, aims to control the discretization error within a frequency band. The approach sweeps from lower to higher frequencies, refining the mesh appropriately so that each of the free vibration mode within the targeted frequency band is sufficiently resolved. Through several numerical examples, we show that the GIFT framework is a powerful and versatile tool to perform local adaptivity in structural dynamics. We also show that the proposed adaptive local h-refinement scheme allows us to achieve significantly faster convergence rates than a uniform h-refinement. [less ▲] Detailed reference viewed: 159 (5 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximaTion (GIFT) ; Tomar, Satyendra ; et al in International Journal for Numerical Methods in Engineering (2018) 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 ... [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 ▲] Detailed reference viewed: 216 (8 UL)Real-time Error Control for Surgical Simulation: Application to Percutaneous Interventions Bui, Huu Phuoc ; Tomar, Satyendra ; et al Presentation (2017, August) Detailed reference viewed: 101 (11 UL)Real-time Error Control for Surgical Simulation Bui, Huu Phuoc ; Tomar, Satyendra ; Bordas, Stéphane E-print/Working paper (2017) Real-time simulations are becoming increasingly common for various applications, from geometric design to medical simulation. Two of the main factors concurrently involved in defining the accuracy of ... [more ▼] Real-time simulations are becoming increasingly common for various applications, from geometric design to medical simulation. Two of the main factors concurrently involved in defining the accuracy of surgical simulations are: the modeling error and the discretization error. Most work in the area has been looking at the above sources of error as a compounded, lumped, overall error. Little or no work has been done to discriminate between modeling error (e.g. needle-tissue interaction, choice of constitutive models) and discretization error (use of approximation methods like FEM). However, it is impossible to validate the complete surgical simulation approach and, more importantly, to understand the sources of error, without evaluating both the discretization error and the modeling error. Our objective is thus to devise a robust and fast approach to measure the discretization error via a posteriori error estimates, which are then used for local remeshing in surgical simulations. To ensure that the approach can be used in clinical practice, the method should be robust enough to deal, as realistically as possible, with the interaction of surgical tools with the organ, and fast enough for real-time simulations. The approach should also lead to an improved convergence so that an economical mesh is obtained at each time step. The final goal is to achieve optimal convergence and the most economical mesh, which will be studied in our future work. [less ▲] Detailed reference viewed: 588 (26 UL)Real-time error control for surgical simulation Bui, Huu Phuoc ; Tomar, Satyendra ; et al Poster (2016, December 12) Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use ... [more ▼] Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational elasticity and a frictional needle/tissue interaction model based on friction. The problem is solved using finite elements within SOFA. The refinement strategy relies upon a hexahedron-based finite element method, combined with a posteriori error estimation driven local $h$-refinement, for simulating soft tissue deformation. Results: We control the local and global error level in the mechanical fields (e.g. displacement or stresses) during the simulation. We show the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. Conclusions: Error control guarantees that a tolerable error level is not exceeded during the simulations. Local mesh refinement accelerates simulations. Significance: Our work provides a first step to discriminate between discretization error and modeling error by providing a robust quantification of discretization error during simulations. [less ▲] Detailed reference viewed: 243 (20 UL)Simulating topological changes in real time for surgical assistance Bordas, Stéphane ; ; et al Speeches/Talks (2016) Detailed reference viewed: 442 (38 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis Tomar, Satyendra ; ; et al Presentation (2016, June 07) In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry ... [more ▼] In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is employed for the numerical solution of the problem over the domain. However, there are situations where this tight integration is a bane rather than a boon. Such situations arise where, e.g., (1) the geometry of the domain is simple enough to be represented by low order NURBS, whereas the unknown (exact) solution of the problem is sufficiently regular, and thus, the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry, (2) the constraint of using the same space for the geometry and the numerical solution is particularly undesirable, such as in the shape and topology optimization, and (3) the solution of the problem has low regularity but for the curved boundary of the domain one can employ higher order NURBS. Therefore, we propose to weaken this constraint. An extensive study of patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution will be discussed. It will be shown, with concrete reasoning, that why patch test fails in certain cases, and that those cases should be avoided in practice. Thereafter, selective numerical examples will be presented to address some of the above-mentioned situations, and it will be shown that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces, and thus, improved accuracy of the numerical solution. [less ▲] Detailed reference viewed: 173 (9 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis Bordas, Stéphane ; Tomar, Satyendra ; et al Scientific Conference (2016, June 05) In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry ... [more ▼] In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is employed for the numerical solution of the problem over the domain. However, there are situations where this tight integration is a bane rather than a boon. Such situations arise where, e.g., (1) the geometry of the domain is simple enough to be represented by low order NURBS, whereas the unknown (exact) solution of the problem is sufficiently regular, and thus, the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry, (2) the constraint of using the same space for the geometry and the numerical solution is particularly undesirable, such as in the shape and topology optimization, and (3) the solution of the problem has low regularity but for the curved boundary of the domain one can employ higher order NURBS. Therefore, we propose to weaken this constraint. An extensive study of patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution will be discussed. It will be shown, with concrete reasoning, that why patch test fails in certain cases, and that those cases should be avoided in practice. Thereafter, selective numerical examples will be presented to address some of the above-mentioned situations, and it will be shown that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces, and thus, improved accuracy of the numerical solution. Powered by [less ▲] Detailed reference viewed: 150 (5 UL)Generalizing the isogeometric concept: weakening the tight coupling between geometry and simulation in IGA Tomar, Satyendra ; ; et al Presentation (2016, June 02) In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the ... [more ▼] In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is also employed for the numerical solution of the problem over the domain. However, in certain situations, such as, when the geometry of the domain can be represented by low order NURBS but the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry; or in the shape and topology optimization where the constraint of using the same space for the geometry and the numerical solution is not favorable, this tight coupling is disadvantageous. Therefore, we study the effect of decoupling the spaces for the geometry representation and the numerical solution, though still using the prevalent functions in CAD/CAGD. To begin with, we perform the patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution. This shows that certain cases, perhaps intuitive, should be avoided in practice because patch test fails. The above-mentioned situations are further explored with some numerical examples, which shows that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces. [1] J. Cottrell, T.J.R. Hughes, and Y. Bazilevs. Isogeometric Analysis: Toward Integration of CAD and FEA, volume 80. Wiley, Chichester, 2009. [2] T.J.R. Hughes, J. Cottrell, and Y. Bazilevs. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 194:4135–4195, 2005. [less ▲] Detailed reference viewed: 180 (11 UL)Linear smoothing over arbitrary polytopes for compressible and nearly incompressible linear elasticity ; Tomar, Satyendra ; Bordas, Stéphane et al Scientific Conference (2016, June) We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. In this approach, a volume-averaged nodal projection operator is ... [more ▼] We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. In this approach, a volume-averaged nodal projection operator is constructed to project the dilatational strain into an approximation space of equal or lower-order than the approximation space for the displacement field, resulting in a locking-free method. The formulation uses the usual Wachspress interpolants over arbitrary polytopes and the stability of the method is ensured by the addition of bubble like functions. The smoothed strains are evaluated based on the linear smoothing procedure. This further softens the bilinear form allowing the procedure to search for a solution satisfying the divergence- free condition. The divergence-free condition of the proposed approach is verified through systematic numerical study. The formulation delivers optimal convergence rates in the energy and L2-norms. Inf-sup tests are presented to demonstrated the stability of the formulation. [less ▲] Detailed reference viewed: 236 (3 UL)Generalizing the isogeometric concept: weakening the tight coupling between geometry and simulation in IGA Bordas, Stéphane ; Tomar, Satyendra ; et al Scientific Conference (2016, May 30) In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the ... [more ▼] In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is also employed for the numerical solution of the problem over the domain. However, in certain situations, such as, when the geometry of the domain can be represented by low order NURBS but the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry; or in the shape and topology optimization where the constraint of using the same space for the geometry and the numerical solution is not favorable, this tight coupling is disadvantageous. Therefore, we study the effect of decoupling the spaces for the geometry representation and the numerical solution, though still using the prevalent functions in CAD/CAGD. To begin with, we perform the patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution. This shows that certain cases, perhaps intuitive, should be avoided in practice because patch test fails. The above-mentioned situations are further explored with some numerical examples, which shows that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces. [less ▲] Detailed reference viewed: 145 (3 UL)On the convergence of stresses in fretting fatigue ; Bordas, Stéphane ; Tomar, Satyendra et al in Materials (2016), 9(8), Fretting is a phenomenon that occurs at the contacts of surfaces that are subjected to oscillatory relative movement of small amplitudes. Depending on service conditions, fretting may significantly reduce ... [more ▼] Fretting is a phenomenon that occurs at the contacts of surfaces that are subjected to oscillatory relative movement of small amplitudes. Depending on service conditions, fretting may significantly reduce the service life of a component due to fretting fatigue. In this regard, the analysis of stresses at contact is of great importance for predicting the lifetime of components. However, due to the complexity of the fretting phenomenon, analytical solutions are available for very selective situations and finite element (FE) analysis has become an attractive tool to evaluate stresses and to study fretting problems. Recent laboratory studies in fretting fatigue suggested the presence of stress singularities in the stick-slip zone. In this paper, we constructed finite element models, with different element sizes, in order to verify the existence of stress singularity under fretting conditions. Based on our results, we did not find any singularity for the considered loading conditions and coefficients of friction. Since no singularity was found, the present paper also provides some comments regarding the convergence rate. Our analyses showed that the convergence rate in stress components depends on coefficient of friction, implying that this rate also depends on the loading condition. It was also observed that errors can be relatively high for cases with a high coefficient of friction, suggesting the importance of mesh refinement in these situations. Although the accuracy of the FE analysis is very important for satisfactory predictions, most of the studies in the literature rarely provide information regarding the level of error in simulations. Thus, some recommendations of mesh sizes for those who wish to perform FE analysis of fretting problems are provided for different levels of accuracy. [less ▲] Detailed reference viewed: 105 (2 UL)Controlling the Error on Target Motion through Real-time Mesh Adaptation: Applications to Deep Brain Stimulation Bui, Huu Phuoc ; Tomar, Satyendra ; et al E-print/Working paper (n.d.) 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 ... [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 ▲] Detailed reference viewed: 425 (34 UL)Real-time error controlled adaptive mesh refinement in surgical simulation: Application to needle insertion simulation ; Tomar, Satyendra ; et al in IEEE Transactions on Biomedical Engineering (n.d.) This paper presents the first real-time discretisation-error-driven adaptive finite element approach for corotational elasticity problems involving strain localisation. We propose a hexahedron-based ... [more ▼] This paper presents the first real-time discretisation-error-driven adaptive finite element approach for corotational elasticity problems involving strain localisation. We propose a hexahedron-based finite element method combined with local oct-tree $h$-refinement, driven by a posteriori error estimation, for simulating soft tissue deformation. This enables to control the local error and global error level in the mechanical fields during the simulation. The local error level is used to refine the mesh only where it is needed, while maintaining a coarser mesh elsewhere. We investigate the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. [less ▲] Detailed reference viewed: 645 (63 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximaTion (GIFT) ; ; Tomar, Satyendra et al E-print/Working paper (n.d.) 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 ... [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 ▲] Detailed reference viewed: 274 (13 UL) |
||