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

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: 84 (14 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: 129 (18 UL)Simulating topological changes in real time for surgical assistance Bordas, Stéphane ; ; et al Speeches/Talks (2016) Detailed reference viewed: 310 (37 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: 105 (7 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: 82 (4 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: 80 (10 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: 131 (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: 76 (3 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: 368 (56 UL)Controlling the Error on Target Motion through Real-time Mesh Adaptation: Applications to Deep Brain Stimulation Bui, Huu Phuoc ; Tomar, Satyendra ; et al in Submitted to Frontiers (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: 38 (16 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 one approach to generalize the concept of isogeometric analysis (IGA) by allowing different spaces for parameterization of the computational domain and for approximation of the ... [more ▼] This paper presents one 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: 17 (5 UL) |
||