The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem
Autor: | Andreas Nüßing, Christian Engwer, Carsten H. Wolters, Heinrich Brinck |
---|---|
Rok vydání: | 2016 |
Předmět: |
Adult
Male FOS: Computer and information sciences J.3 Computer science Pipeline (computing) Finite Element Analysis Biomedical Engineering G.1.8 010103 numerical & computational mathematics 35J25 35J75 35Q92 65N30 68U20 92C50 01 natural sciences Models Biological Computational Engineering Finance and Science (cs.CE) 03 medical and health sciences 0302 clinical medicine Level set G.1.10 Discontinuous Galerkin method FOS: Mathematics Humans Polygon mesh I.6.0 Mathematics - Numerical Analysis 0101 mathematics Computer Science - Computational Engineering Finance and Science ComputingMethodologies_COMPUTERGRAPHICS Conservation law Models Statistical Triangulation (social science) Brain Electroencephalography Numerical Analysis (math.NA) Finite element method Quantitative Biology - Neurons and Cognition FOS: Biological sciences Neurons and Cognition (q-bio.NC) Hexahedron Algorithm Head 030217 neurology & neurosurgery |
Zdroj: | IEEE Transactions on Biomedical Engineering |
DOI: | 10.48550/arxiv.1601.07810 |
Popis: | Objective: The purpose of this study is to introduce and evaluate the unfitted discontinuous Galerkin finite element method (UDG-FEM) for solving the electroencephalography (EEG) forward problem. Methods: This new approach for source analysis does not use a geometry conforming volume triangulation, but instead uses a structured mesh that does not resolve the geometry. The geometry is described using level set functions and is incorporated implicitly in its mathematical formulation. As no triangulation is necessary, the complexity of a simulation pipeline and the need for manual interaction for patient specific simulations can be reduced and is comparable with that of the FEM for hexahedral meshes. In addition, it maintains conservation laws on a discrete level. Here, we present the theory for UDG-FEM forward modeling, its verification using quasi-analytical solutions in multi-layer sphere models and an evaluation in a comparison with a discontinuous Galerkin (DG-FEM) method on hexahedral and on conforming tetrahedral meshes. We furthermore apply the UDG-FEM forward approach in a realistic head model simulation study. Results: The given results show convergence and indicate a good overall accuracy of the UDG-FEM approach. UDG-FEM performs comparable or even better than DG-FEM on a conforming tetrahedral mesh while providing a less complex simulation pipeline. When compared to DG-FEM on hexahedral meshes, an overall better accuracy is achieved. Conclusion: The UDG-FEM approach is an accurate, flexible and promising method to solve the EEG forward problem. Significance: This study shows the first application of the UDG-FEM approach to the EEG forward problem. Comment: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible |
Databáze: | OpenAIRE |
Externí odkaz: |