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.

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

Three-Dimensional Mesh Generation for Human Heart Model

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

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

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

Finite Element Modeling of Human Head from MRI

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

Facial Hexahedral Mesh Transferring by Volumetric Mapping Based on Harmonic Fields

Fast and high-quality tetrahedral mesh generation from neuroanatomical scans

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

Numerical Simulation of Fluid-Structure Interaction Relations Between Skull and Brain Based on ALE and Overlapping Mesh Methods

Meshfree Implementation of Individualized Active Cardiac Dynamics.

Accuracy and efficiency of finite element head models: The role of finite element formulation and material laws

A Semi-Automatic Method of Generating Subject-Specific Pediatric Head Finite Element Models for Impact Dynamic Responses to Head Injury

Sheet Operation Based Block Decomposition of Solid Models for Hex Meshing.

A Fast Methodology for Generating Skeletal FEM with Detailed Human Geometric Features Based on CPD and RBF Algorithms

Development and Preliminary Validation of a Parametric Pediatric Head Finite Element Model for Population-Based Impact Simulations

Automated Block Decomposition of Solid Models Based on Sheet Operations

The generation of tetrahedral mesh models for neuroanatomical MRI

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

Development of a Three-dimensional Finite Element Model of the Human Brain Based on CT