direct image analysis methods for surgical simulation and
play

Direct image-analysis methods for surgical simulation and mixed - PowerPoint PPT Presentation

Direct image-analysis methods for surgical simulation and mixed meshfree methods Jack S. Hale Marie Curie Postdoctoral Fellow Research Unit in Engineering Science Universit du Luxembourg Collaborators: Stephane P. A. Bordas, Universite du


  1. Direct image-analysis methods for surgical simulation and mixed meshfree methods Jack S. Hale Marie Curie Postdoctoral Fellow Research Unit in Engineering Science Université du Luxembourg Collaborators: Stephane P. A. Bordas, Universite du Luxembourg/Cardiff University Pierre Kerfriden, Cardiff University Juan J. Rodenas Garcia, Universidad Politecnica de Valencia � Pedro M. Baiz, Imperial College London Alejandro Ortiz-Bernardin, Universidad de Chile Christian J. Cyron, Yale University Université Pierre-et-Marie-Curie Institut Jean le Rond d’Alembert 26th May 1400

  2. An Overview

  3. The problems and my research Mixed Meshfree variational methods forms Numerical Direct image- locking in Mesh distortion in analysis transition incompressible hyperelastic elasticity, plates and for surgical problems shells, Cosserat simulators. elasticity Enriched implicit Patch projection techniques boundary methods

  4. Project 1: Direct image to mesh analysis for surgical simulation

  5. Enriched and implicit High Performance boundary methods Computing Surgical Simulation Model-order Real-time simulation reduction

  6. RealTCut and SHACRA Inria

  7. How can we move from an image… COLONIX, OSIRIX

  8. …or perhaps a series of images… Source: COLONIX, OSIRIX

  9. Pipelines to analysis Traditional Implicit Boundary } } Acquire images Acquire images Implicit Segment images Segment images { NURBS Mesh Surfaces Implicit Explicit Mesh Volume Perform analysis Perform analysis

  10. Nested Octree Discretisation Geometry O g O d M

  11. How to transfer geometric information back to the discretisation? Enrichment E [ V p g V p d � h d ( O d ) h g ( O g )] V p g V p d p d > p g h d ( O d ) h g ( O g ) h d > h g p g = 1

  12. Octree Level 5/Level 3

  13. Surface 6a2e86c

  14. Features • 2D and 3D problems on the same code-path. • parallel hybrid MPI/TBB assembly and solution. • fast and robust computational geometry using CGAL. • automatic Delaunay tessellation of integration subdomains. • completely separate representation of discretisation and geometry via nested octree data structures. • constructs to represent soft and hard segmentations of image data. • implicit representation via level-sets, inside-outside functions. • independent hpe-type adaptivity on discretisation and geometry. • fast refinement and coarsening operations.

  15. Outlook • We are developing a cartesian grid implicit boundary/enriched finite element method toolkit within deal.ii. • By uncoupling discretisation and geometry we hope to produce a method particularly suited to image-based analysis. • Many challenges ahead, particularly with imposition of Dirichlet boundary conditions on hyperelastic soft-tissue model.

  16. Project 2: Mesh-based and Meshfree Partitions of Unity

  17. Partition of Unity N � u h ( x ) = φ i u i i =1 N � φ i = 1 i =1

  18. �� � �� �������� �� �������� ������� ������ ���� ����� ��������� ����������� ����� ���������� ��� ������ �� ���� �������� ���� � ������ ������ ����� ��� ������� ������� ����� ���� ��� ���� �� ��������������� �������� ��� ��� ������������� ������������� �� ��� ����� ��������� ���������� ���� ��� �������� ����� �� ��� �������� �� ��� ������� �� ��� �������� ��� �� ��� ������ �� ��� ������� ������ ������ ����� �� ���������� ��� ����� ��� ��� ��������������� �������� ���� Finite elements 1 . 0 0 . 8 0 . 6 φ 0 . 4 0 . 2 0 . 0 0 . 0 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 x ������ ���� ����� ��������� ����������� ����� � � ����� ������� ������ �� ���� �������� ���� �

  19. �� ��������� ������ ������� ���� � �������� ������� � �� � ������� ������� �� ��� ����� ������ ���� �� ��� ����� ������� ������ ����� ��������� ��� ����� �� ��� ��������� ������� ���� ������ ���� ���������� ��������� �� ����� ������������ ��������� ����� ������ ������ �� ��� ���������� ��� �� ���� ����� �������� �� ��� ��� �� ������� ��� ���������� ������ ��� ������ ��� ���� ������ ������ ������ ������� ��� ��������� �� ��� ������� � ������ ������� ��� � ������������ Mesh-based Partition of Unity

  20. �� � ������������ ������ ������� ���� � �������� ������� � �� � ������� ������� �� ��� ����� ������ ���� �� ��� ����� ������� ������ ����� ��������� ��� ����� �� ��� ��������� ������� ���� ������ ���� ���������� ��������� �� ����� ������������ ��������� ����� ������ ������ �� ��� ���������� ��� �� ���� ����� �������� �� ��� ��� �� ������� ��� ���������� ������ ��� ������ ��������� ��� ���� ������ ������ ������ ������� ��� ��������� �� ��� ������� � ������ ������� ��� Finite Element Method Reference Mesh

  21. �� ������ ������ ���� �������� ��������� �� ����� ������������ ��������� ��������� ��� ������� �������� ������� ����������� �� ��� ���� ��� �� ���� ����� �������� �� � ��� �� ������� ��� ���������� ���� ������� ��� ����� ���� �� ����� � ���� ���� ������ ������ ������� ��� ��������� �� ��� �������� � ������ ������� ��� � ������������ Meshfree-based Partition of Unity

  22. Meshfree partition of unity is constructed in global space Good: Not so good: improved approximation expensive. • • properties. extra flexibility also means extra • high-order continuity. complexity: • less dofs for given accuracy. integration. � • • less sensitive to poor quality adaptivity/error measures. • • node distributions. � mathematical proof. • eases mesh generation. stability. � • • boundary conditions. � • �

Recommend


More recommend