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.

Finite Element Refinements for Inverse Electrocardiography: Hybrid Shaped Elements and High-Order Element Truncation

Noninvasive Electroardiographic Imaging: Application Of Hybrid Methods For Solving The Electrocardiography Inverse Problem

Solving the ECG forward problem by means of standard h- and h-hierarchical adaptive linear boundary element method: comparison with two refinement schemes.

Solving the Electrocardiography Inverse Problem by Using an Optimal Algorithm Based on the Total Least Squares Theory

Truncated total least squares: a new regularization method for the solution of ECG inverse problems.

Combination of the LSQR Method and a Genetic Algorithm for Solving the Electrocardiography Inverse Problem

Combining Regularization Frameworks for Solving the Electrocardiography Inverse Problem.

The Application of Subspace Preconditioned LSQR Algorithm for Solving the Electrocardiography Inverse Problem

Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

On epicardial potential reconstruction using regularization schemes with the L1-norm data term.

[EM Algorithm for the Inverse Problem of Electrocardiography].

High-order spectral/ hp element discretisation for reaction-diffusion problems on surfaces: Application to cardiac electrophysiology

Simulation of Active Cardiac Electromechanical Dynamics

How Much Geometrical Detail Do We Need In Cardiac Electrophysiological Imaging? A Generic Heart-Torso Representation For Fast Subject-Specific Customization

Meshfree Implementation of Individualized Active Cardiac Dynamics.

Conforming Discretizations of Boundary Element Solutions of the Electroencephalography Forward Problem

The Importance of FEM Model for Linearized EIT Image Reconstruction.

Finite element-based space-time total variation-type regularization of the inverse problem in electrocardiographic imaging

A hybrid finite element/finite difference method for reconstruction of dielectric properties of conductive objects

A Two-Dimensional Eight-Node Quadrilateral Inverse Element for Shape Sensing and Structural Health Monitoring

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