Software Package for Transcranial Magnetic Stimulation Coil and Coil Array Analysis and Design
Morales, L.; Wartman, W. A.; Ferreira, J.; Miles, A.; Daneshzand, M.; Lu, H. A.; Nummenmaa, A.; Deng, Z.-D.; Makaroff, S. N.
Show abstract
ObjectiveThis study aims to describe a MATLAB software package for transcranial magnetic stimulation (TMS) coil analysis and design. ApproachElectric and magnetic fields of the coils as well as their self- and mutual (for coil arrays) inductances are computed, with or without a magnetic core. Solid and stranded (Litz wire) conductors are also taken into consideration. The starting point is the centerline of a coil conductor(s), which is a 3D curve defined by the user. Then, a wire mesh and a computer aided design (CAD) mesh for the volume conductor of a given cross-section (circular, elliptical, or rectangular) are automatically generated. Self- and mutual inductances of the coil(s) are computed. Given the conductor current and its time derivative, electric and magnetic fields of the coil(s) are determined anywhere in space. Computations are performed with the fast multipole method (FMM), which is the most efficient way to evaluate the fields of many elementary current elements (current dipoles) comprising the current carrying conductor at a large number of observation points. This is the major underlying mathematical operation behind both inductance and field calculations. Main ResultsThe wire-based approach enables precise replication of even the most complex physical conductor geometries, while the FMM acceleration quickly evaluates large quantities of elementary current filaments. Agreement to within 0.74% was obtained between the inductances computed by the FMM method and ANSYS Maxwell 3D for the same coil model. Although not provided in this study, it is possible to evaluate non-linear magnetic cores in addition to the linear core exemplified. An experimental comparison was carried out against a physical MagVenture C-B60 coil; the measured and simulated inductances differed by only 1.25%, and nearly perfect correlation was found between the measured and computed E-field values at each observation point. SignificanceThe developed software package is applicable to any quasistatic inductor design, not necessarily to the TMS coils only.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- Evaluating computational efforts and physiological resolution of mathematical models of cardiac tissue 95%
- Nonlinear Dispersive Cell Model for Microdosimetry of Nanosecond Pulsed Electric Fields 93%
- A 1D Model Characterizing the Role of Spatiotemporal Contraction Distributions on Lymph Transport 92%
Similar papers in this journal
- Optimizing the Electric Field Strength in Multiple Targets for Multichannel Transcranial Electric Stimulation 96%
- Modeling implanted metals in electrical stimulation applications 94%
- Non-invasive stimulation with Temporal Interference: Optimization of the electric field deep in the brain with the use of a genetic algorithm 94%
Similar papers in this journal
- Advanced Biophysical Model to Capture Channel Variability for EQS Capacitive HBC 94%
- Bio-Physical Modeling of Galvanic Human Body Communication in Electro-Quasistatic Regime 94%
- Identifiability analysis and noninvasive online estimation of the first-order neural activation dynamics in the brain with closed-loop transcranial magnetic stimulation 94%
Similar papers in this journal
- Design and execution of a Verification, Validation, and Uncertainty Quantification plan for a numerical model of left ventricular flow after LVAD implantation 94%
- BlueRecording: A Pipeline for the efficient calculation of extracellular recordings in large-scale neural circuit models 94%
- Nano-scale solution of the Poisson-Nernst-Planck (PNP) equations in a fraction of two neighboring cells reveals the magnitude of intercellular electrochemical waves 93%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.