Novelle approach to simulating spinal cord stimulation during tSCS using CT images and FEM

Abstract Transcutaneous spinal cord stimulation (tSCS) offers non-invasive relief for chronic pain and improves motor function in spinal cord injured (SCI) patients. However, its mechanisms are not currently fully understood, and patientspecific factors, such as Body mass index (BMI) and age, complicate treatments. This paper aims to understand tSCS better by developing a novel Finite Element Model (FEM) of the human body using CT scans. Three subjects (sex, male, female, male. Age: 26, 27, 64. BMI: 38.9, 24.1, 28.4) underwent a CT scan, performed on a Cannon Aquilion Prime (Slice thickness [mm]: 0.8, Voxel size [mm 3 ]: 0.564, 0.328, 0.527), which imaged the trunk of the body, from top of the abdomen to the bottom of the pelvis. The images were then used in Materialize Mimics Research 21.0 to create 3D images of individual organs, skin, fat, muscles, skeleton, and spinal cord. After pre-processing in Autodesk Meshmixer, the models were converted into solid CAD objects in Ansys SpaceClaim R2021 and combined into a single abdominal model. Ansys Maxwell R2021 was then used for simulations. Five different two-electrode configurations were tested in the prototype phase with a simplified model setup. The positive electrodes were placed over the (Thoracic) T10, T12, (Lumbar) L2 and L4 vertebrae sequentially, with the negative electrode over (Sacral) S2. The simulations took on average 47.4 hours. A marked decline in electrical current penetration depth was observed as the electrodes were placed closer together which was consistent with known current distribution patterns. A preliminary validation test was also performed using a lamb’s thigh. Two electrodes were placed a known distance apart and a stimulation was given, needle electrodes were then inserted in a grid-like pattern to obtain voltage values. The setup was recreated and simulated in Ansys Maxwell. The resulting average percent difference was 44.93 ± 33.3 % (5Vpp Square wave) and 35.02 ± 23.38 % (5V DC). In both instances, the highest difference was at the edges of the electrodes and the lowest difference in the midpoint between electrodes. Later versions of FEM models incorporated more organs and had improved on previous mesh generation flaws, but encountered new mesh generation errors which could not be rectified before the conclusion of the master’s project. Despite complications, this project has provided a pipeline for creating similar models and shown their usability. Future work will involve overcoming the current mesh generation errors, reducing calculation time, and performing a thorough validation test.