Intervertebral discs (IVD) provide flexibility to the back and ensure functional distributions of the spinal loads. They are avascular, and internal diffusion-dependent metabolic transport is vital to supply nutrients to disc cells1, but interactions with personalized IVD shapes and mechanics remain poorly explored. Poromechanical finite element models of seven personalized lumbar IVD geometries, with mean heights ranging from 8 to 16 mm were coupled with a reactive oxygen, glucose and lactate transport model linked with tissue deformations and osmosis . In previous studies, reduced formulations of the divergence of the solute flux (∇ .