Darvishi, V.; Chan, E. Y. K.; Duckworth, H.; Parker, T. D.; Sharp, D. J.; Ghajari, M.

Converting medical images into anatomically detailed, subject-specific finite element (FE) models is a long-standing bottleneck in brain computational modelling. These models are used to predict brain tissue deformation, e.g. in traumatic brain injury, particle diffusion in brain drug delivery, and other biophysical phenomena across neurological disorders. However, existing model creation workflows depend on manual image segmentation, proprietary meshing software, and labour-intensive repair of meningeal and interface structures, limiting reproducibility and cohort analysis. Here we present PARS, a fully automated, open-source pipeline that converts a T1-weighted MRI scan into a simulation-ready FE head model. PARS combines anatomical parcellation with tissue maps and uses iterative neighbourhood-based reclassification, yielding a gap-free whole-head label volume. The volume is directly converted into a hexahedral mesh, augmented with algorithmically reconstructed falx, tentorium, pia and dura mater, and refined by Laplacian smoothing under a node-locking scheme that controls element quality and the explicit-solver stable timestep. We evaluated PARS on 23 subjects spanning cranial volumes of 832 to 1,329 cubic centimeter, at 1.0, 1.5 and 2.0 mm MRI resolutions. At 1 mm, meshes achieved a median Scaled Jacobian of 0.976, and total intracranial volume error of ~0.54; quality remained high at 1.5 mm (SJ of 0.933) and 2 mm (SJ of 0.921). Model creation runtime ranged from 9 to 38 minutes per subject. Models generated by PARS have been validated against cadaveric brain displacement data and demonstrated utility across traumatic brain injury and normal pressure hydrocephalus research. PARS provides an open-access, reproducible resource that substantially lowers the barriers to subject-specific brain modelling.