Abstract:Objective In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems. Approach We propose, describe, and systematically investigate an AMR method using the Boundary Element Method with Fast Multipole Acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy. The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a "silver-standard" solution found by subsequent 4:1 global refinement of the adaptively-refined model. Main Results Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models. Significance This error has especially important implications for TES dosing prediction – where the stimulation strength plays a central role – and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.

Open Meshed Anatomy: Towards a comprehensive finite element hexahedral mesh derived from open atlases

Three-Dimensional Mesh Generation for Human Heart Model

A Finite Element Mesh Generation Algorithm and Application to Human Head Modeling

Hexahedral mesh of anatomical atlas for construction of computational human brain models: Applications to modeling biomechanics and bioelectric field propagation

Fast and high-quality tetrahedral mesh generation from neuroanatomical scans

Image-based Meshing of Patient-Specific Data:Converting Medical Scans into Highly Accurate Computational Models

Development of High-Quality Hexahedral Human Brain Meshes Using Feature-Based Multi-Block Approach

Facial Hexahedral Mesh Transferring by Volumetric Mapping Based on Harmonic Fields

Meshfree Implementation of Individualized Active Cardiac Dynamics.

Delaunay Tetrahedralization of the Heart Based on Integration of Open Source Codes

A Meshfree Representation for Cardiac Medical Image Computing.

3D Anatomical Shape Atlas Construction Using Mesh Quality Preserved Deformable Models

Multi-compartment human head modeling: generating adaptive tetrahedral mesh with GPU acceleration

Constrained tetrahedral mesh generation of human organs on segmented volume

Improving Model-Based Functional Near-Infrared Spectroscopy Analysis Using Mesh-Based Anatomical and Light-Transport Models

Techniques for the Generation of 3D Finite Element Meshes of Human Organs

A fast and robust patient specific Finite Element mesh registration technique: application to 60 clinical cases

Volumetric Mesh Construction from Scattered Prior Data: Application to Cardiac MR Image Analysis.

An adaptive H-refinement method for the boundary element fast multipole method for quasi-static electromagnetic modeling

The generation of tetrahedral mesh models for neuroanatomical MRI

Subject specific finite element mesh generation of the pelvis from biplanar x-ray images: application to 120 clinical cases