Beyond finite elements: A comprehensive, patient-specific neurosurgical simulation utilizing a meshless method

    Research output: Contribution to journalArticle

    22 Citations (Scopus)

    Abstract

    To be useful in clinical (surgical) simulations, a method must use fully nonlinear (both geometric and material) formulations to deal with large (finite) deformations of tissues. The method must produce meaningful results in a short time on consumer hardware and not require significant manual work while discretizing the problem domain. In this paper, we showcase the Meshless Total Lagrangian Explicit Dynamics Method (MTLED) which meets these requirements, and use it for computing brain deformations during surgery. The problem geometry is based on patient-specific MRI data and includes the parenchyma, tumor, ventricles and skull. Nodes are distributed automatically through the domain rendering the normally difficult problem of creating a patient-specific computational grid a trivial exercise. Integration is performed over a simple, regular background grid which does not need to conform to the geometry boundaries. Appropriate nonlinear material formulation is used. Loading is performed by displacing the parenchyma surface nodes near the craniotomy and a finite frictionless sliding contact is enforced between the skull (rigid) and parenchyma. The meshless simulation results are compared to both intraoperative MRIs and Finite Element Analysis results for multiple 2D sections. We also calculate Hausdorff distances between the computed deformed surfaces of the ventricles and those observed intraoperatively. The difference between previously validated Finite Element results and the meshless results presented here is less than 0.2 mm. The results are within the limits of neurosurgical and imaging equipment accuracy (∼1 mm) and demonstrate the method's ability to fulfill all of the important requirements for surgical simulation. Copyright © 2012 Elsevier Ltd. All rights reserved
    Original languageEnglish
    Pages (from-to)2698–2701
    JournalJournal of Biomechanics
    Volume45
    Issue number15
    Early online date27 Aug 2012
    DOIs
    Publication statusPublished - 11 Oct 2012

    Fingerprint

    Magnetic resonance imaging
    Geometry
    Skull
    Surgery
    Tumors
    Brain
    Tissue
    Hardware
    Imaging techniques
    Finite element method
    Finite Element Analysis
    Craniotomy
    Exercise
    Equipment and Supplies
    Neoplasms

    Cite this

    @article{a0019a206d054ccc9c32995410df49fa,
    title = "Beyond finite elements: A comprehensive, patient-specific neurosurgical simulation utilizing a meshless method",
    abstract = "To be useful in clinical (surgical) simulations, a method must use fully nonlinear (both geometric and material) formulations to deal with large (finite) deformations of tissues. The method must produce meaningful results in a short time on consumer hardware and not require significant manual work while discretizing the problem domain. In this paper, we showcase the Meshless Total Lagrangian Explicit Dynamics Method (MTLED) which meets these requirements, and use it for computing brain deformations during surgery. The problem geometry is based on patient-specific MRI data and includes the parenchyma, tumor, ventricles and skull. Nodes are distributed automatically through the domain rendering the normally difficult problem of creating a patient-specific computational grid a trivial exercise. Integration is performed over a simple, regular background grid which does not need to conform to the geometry boundaries. Appropriate nonlinear material formulation is used. Loading is performed by displacing the parenchyma surface nodes near the craniotomy and a finite frictionless sliding contact is enforced between the skull (rigid) and parenchyma. The meshless simulation results are compared to both intraoperative MRIs and Finite Element Analysis results for multiple 2D sections. We also calculate Hausdorff distances between the computed deformed surfaces of the ventricles and those observed intraoperatively. The difference between previously validated Finite Element results and the meshless results presented here is less than 0.2 mm. The results are within the limits of neurosurgical and imaging equipment accuracy (∼1 mm) and demonstrate the method's ability to fulfill all of the important requirements for surgical simulation. Copyright {\circledC} 2012 Elsevier Ltd. All rights reserved",
    author = "Karol Miller and A. Horton and Grand Joldes and Adam Wittek",
    year = "2012",
    month = "10",
    day = "11",
    doi = "10.1016/j.jbiomech.2012.07.031",
    language = "English",
    volume = "45",
    pages = "2698–2701",
    journal = "Journal of Biomechanics",
    issn = "0021-9290",
    publisher = "Elsevier",
    number = "15",

    }

    Beyond finite elements: A comprehensive, patient-specific neurosurgical simulation utilizing a meshless method. / Miller, Karol; Horton, A.; Joldes, Grand; Wittek, Adam.

    In: Journal of Biomechanics, Vol. 45, No. 15, 11.10.2012, p. 2698–2701.

    Research output: Contribution to journalArticle

    TY - JOUR

    T1 - Beyond finite elements: A comprehensive, patient-specific neurosurgical simulation utilizing a meshless method

    AU - Miller, Karol

    AU - Horton, A.

    AU - Joldes, Grand

    AU - Wittek, Adam

    PY - 2012/10/11

    Y1 - 2012/10/11

    N2 - To be useful in clinical (surgical) simulations, a method must use fully nonlinear (both geometric and material) formulations to deal with large (finite) deformations of tissues. The method must produce meaningful results in a short time on consumer hardware and not require significant manual work while discretizing the problem domain. In this paper, we showcase the Meshless Total Lagrangian Explicit Dynamics Method (MTLED) which meets these requirements, and use it for computing brain deformations during surgery. The problem geometry is based on patient-specific MRI data and includes the parenchyma, tumor, ventricles and skull. Nodes are distributed automatically through the domain rendering the normally difficult problem of creating a patient-specific computational grid a trivial exercise. Integration is performed over a simple, regular background grid which does not need to conform to the geometry boundaries. Appropriate nonlinear material formulation is used. Loading is performed by displacing the parenchyma surface nodes near the craniotomy and a finite frictionless sliding contact is enforced between the skull (rigid) and parenchyma. The meshless simulation results are compared to both intraoperative MRIs and Finite Element Analysis results for multiple 2D sections. We also calculate Hausdorff distances between the computed deformed surfaces of the ventricles and those observed intraoperatively. The difference between previously validated Finite Element results and the meshless results presented here is less than 0.2 mm. The results are within the limits of neurosurgical and imaging equipment accuracy (∼1 mm) and demonstrate the method's ability to fulfill all of the important requirements for surgical simulation. Copyright © 2012 Elsevier Ltd. All rights reserved

    AB - To be useful in clinical (surgical) simulations, a method must use fully nonlinear (both geometric and material) formulations to deal with large (finite) deformations of tissues. The method must produce meaningful results in a short time on consumer hardware and not require significant manual work while discretizing the problem domain. In this paper, we showcase the Meshless Total Lagrangian Explicit Dynamics Method (MTLED) which meets these requirements, and use it for computing brain deformations during surgery. The problem geometry is based on patient-specific MRI data and includes the parenchyma, tumor, ventricles and skull. Nodes are distributed automatically through the domain rendering the normally difficult problem of creating a patient-specific computational grid a trivial exercise. Integration is performed over a simple, regular background grid which does not need to conform to the geometry boundaries. Appropriate nonlinear material formulation is used. Loading is performed by displacing the parenchyma surface nodes near the craniotomy and a finite frictionless sliding contact is enforced between the skull (rigid) and parenchyma. The meshless simulation results are compared to both intraoperative MRIs and Finite Element Analysis results for multiple 2D sections. We also calculate Hausdorff distances between the computed deformed surfaces of the ventricles and those observed intraoperatively. The difference between previously validated Finite Element results and the meshless results presented here is less than 0.2 mm. The results are within the limits of neurosurgical and imaging equipment accuracy (∼1 mm) and demonstrate the method's ability to fulfill all of the important requirements for surgical simulation. Copyright © 2012 Elsevier Ltd. All rights reserved

    U2 - 10.1016/j.jbiomech.2012.07.031

    DO - 10.1016/j.jbiomech.2012.07.031

    M3 - Article

    VL - 45

    SP - 2698

    EP - 2701

    JO - Journal of Biomechanics

    JF - Journal of Biomechanics

    SN - 0021-9290

    IS - 15

    ER -