Field distributions in vertebral bodies of the rat during electrical stimulation: a parametric study

The electrical field and current density distributions were found in the various tissues of a mathematical model of the experimental rat used to study systemic osteoporosis. The finite element method was used to solve the boundary value problem derived from Maxwell's equations using a quasistatic approximation for a 60 kHz external output signal applied via skin electrodes. A parametric study was done initially to determine the principle factors which effect the solution of the field in the vertebral bodies. Grid coarseness, model length, and intervertebral space width had little effect on the solution while trabecular bone and abdominal cavity conductivity values had strong effects. The two pair of transversely placed electrodes spaced by at least three vertebral bodies produced the most uniform field distributions and was used in the experimental rat model. The range of current density values in the trabecular bone was determined to be 3.0-5.0 microA/cm2 at the external output signal where evidence of a reversal of bone loss due to castration osteoporosis had been found in the experimental rat.