Finite element analysis of electric field distribution during direct current stimulation of the spinal cord: Implications for device design

Spinal cord injury (SCI) arises from damage to the spinal cord, often caused by trauma or disease. The resulting sensorimotor dysfunction is variable and dependent on the extent of the injury. Despite years of research, curative options for SCI remain limited. However, recent advancements in electric field stimulated axonal regrowth have shown promise for neuronal regeneration. One roadblock in the development of therapeutic treatments based on this is a lack of understanding of the exogenous electric field distribution in the injured tissue, and in particular, how this is influenced by electrode geometry and placement. To better understand this electric field, and provide a means by which it can be optimized, we have developed a finite element model of such spinal cord treatment. We investigate the impact of variations in electrode geometry, spinal cord size, and applied current magnitude as well as looking at several injury models in relation to clinically observed outcomes. Through this, we show that electrode shape has little effect on the induced electric field, that the placement of these electrodes has a noticeable influence on the field distribution, and that the magnitude of this field is governed by both the applied current and the spinal cord morphology. We also show that the injury modality influences the induced field distribution and that a stronger understanding of the injury will help decide treatment parameters. This work provides guidance in the design of electrodes for future clinical application in direct current electric field stimulation for axonal regeneration.