Computation of the dynamic compression effects in spine discs using integral methods

Abstract

The computational modeling using integral methods of dynamic loading and its effects on the nutrients transport in spine discs is addressed in this work. The numerical simulation and analysis was carried out using the Boundary Element Method (BEM) and a 3D model (axisymmetric) of the disc. The boundary model was discretized using linear interpolated elements and a multi-region approach. Concentration and production of three nutrients as lactate, oxygen and glucose were obtained. The maximum lactate concentration was observed very close to the interface between the nucleus and the inner annulus. A relatively simple model discretized with 130 boundary elements yielded very similar results to these coming from more complex FEM-based models. The numerical efforts in the domain and boundary discretizations were optimized using the BEM. Our results are in good agreement with those obtained using with finite element-based models. As expected, the dynamic loading increased the oxygen–glucose consumption and the lactate production, thus leading to a poor oxygen–glucose concentration at the nucleus of the disc. All of that is a favorable environment for a disc degeneration mechanism to be developed.