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.

Anisotropic Subvoxel-Smooth Conduction Model For Bioelectromagnetism Analysis

A Study of White Matter Anisotropic Conductivity on BEG Forward Solutions

Simulation Study of A Magnetocardiogram Based on A Virtual Heart Model: Effect of A Cardiac Equivalent Source and A Volume Conductor

An Efficient BiCGstab Solved Impedance Method for Induced Field Evaluation with a Hyperthermia Applicator.

A New Magneto-Cardiogram Study Using A Vector Model with A Virtual Heart and the Boundary Element Method

Towards the Cardiac Equivalent Source Models in Electrocardiography and Magnetocardiography: A Simulation Study

Calculation of EEG Problems with Anisotropic Conducting Media by the Finite Volume Method

Finite Volume Method for Investigating Anisotropic Conductivity in Eeg

A Survey on Brain Tissues Anisotropic Conductivity Model Based on Diffusion Tensor Imaging

Construction of boundary element models from magnetic resonance images for bioelectromagnetic inverse problems

Computation of Potential Distribution in Anisotropic Dielectric Media Sphere Model by Finite Element Method

Analytical Solution to the Three-Dimensional Head Eit Forward Problem: Calculation and Simulation with A 4-Layer Spherical Volume Conductor

Electromagnetic Modeling for FSS with Anisotropic Substrate by Using a Hybrid-Accelerated VSIE Method

3D forward modeling of controlled-source audio-frequency magnetotellurics in arbitrarily anisotropic media

Calculation of the realistic head model with anisotropic conducting media by the finite volume method

Conductivity Anisotropy Influence on Acoustic Sources for Magnetoacoustic Tomography with Magnetic Induction.

Calculation of EEG forward problem with anisotropy by the finite element method

Advancing Electromagnetic Field Modeling in Axially Anisotropic Media with Conformal Method in 3D Cylindrical Coordinate Systems

Research on the 3D mimetic finite volume method for loop-source TEM response in biaxial anisotropic formation

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

An Efficient Parallel Numerical Modeling of Bioheat Transfer in Realistic Tissue Structure