Finite-element time-domain algorithms for modeling linear Debye and Lorentz dielectric dispersions at low frequencies

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.

[1]  T. Preston,et al.  Finite Elements for Electrical Engineers , 1984 .

[2]  R. W. Lau,et al.  The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. , 1996, Physics in medicine and biology.

[3]  William D. Hurt,et al.  Multiterm Debye Dispersion Relations for Permittivity of Muscle , 1985, IEEE Transactions on Biomedical Engineering.

[4]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[5]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[6]  K. Foster,et al.  Dielectric properties of tissues and biological materials: a critical review. , 1989, Critical reviews in biomedical engineering.

[7]  D. B. Heppner,et al.  Considerations of quasi-stationarity in electrophysiological systems. , 1967, The Bulletin of mathematical biophysics.

[8]  D. Coppola,et al.  Electrically mediated delivery of vector plasmid DNA elicits an antitumor effect , 2002, Gene Therapy.

[9]  Karl Johan Åström,et al.  Transient Response Analysis , 1985 .

[10]  Todd A. Kuiken,et al.  Frequency- and time-domain FEM models of EMG: capacitive effects and aspects of dispersion , 2002, IEEE Transactions on Biomedical Engineering.

[11]  A. Taflove,et al.  Finite element modeling of electromagnetic signal propagation in a phantom arm , 2001, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[12]  S. O. Nelson,et al.  Low-frequency dielectric properties of biological tissues : A review with some new insights , 1998 .

[13]  R. W. Lau,et al.  The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. , 1996, Physics in medicine and biology.

[14]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[15]  H-T Huang,et al.  Electroporation-mediated pain-killer gene therapy for mononeuropathic rats , 2002, Gene Therapy.

[16]  H. Nakamura,et al.  Ectopic bone formation by electroporatic transfer of bone morphogenetic protein-4 gene. , 2002, Bone.