title: Finite-element time-domain algorithms for modeling linear Debye and Lorentz dielectric dispersions at low frequencies.
abstract: We present what we believe to be the first algorithms that use a simple scalar-potential formulation to model linear Debye and Lorentz dielectric dispersions at low frequencies in the context of finite-element time-domain (FETD) numerical solutions of electric potential. The new algorithms, which permit treatment of multiple-pole dielectric relaxations, are based on the auxiliary differential equation method and are unconditionally stable. We validate the algorithms by comparison with the results of a previously reported method based on the Fourier transform. The new algorithms should be useful in calculating the transient response of biological materials subject to impulsive excitation. Potential applications include FETD modeling of electromyography, functional electrical stimulation, defibrillation, and effects of lightning and impulsive electric shock.
authors.affiliation: Rehabilitation Institute of Chicago, IL 60611, USA. n-stoykov@northwestern.edu
authors.forename: Nikolay S
authors.initials: NS
authors.lastname: Stoykov
authors.forename: Todd A
authors.initials: TA
authors.lastname: Kuiken
authors.forename: Madeleine M
authors.initials: MM
authors.lastname: Lowery
authors.forename: Allen
authors.initials: A
authors.lastname: Taflove
mesh_terms: D000465:Algorithms; D003198:Computer Simulation; D004553:Electric Conductivity; D004574:Electromagnetic Fields; D020342:Finite Element Analysis; D008954:Models, Biological; D008959:Models, Neurological; D011874:Radiometry / Q000379:methods*; D015203:Reproducibility of Results; D012542:Scattering, Radiation; D012680:Sensitivity and Specificity