Intrathecal drug and gene vector delivery is a procedure to release a solute within the cerebrospinal fluid. This procedure is currently used in clinical practice and shows promise for treatment of several central nervous system pathologies. However, intrathecal delivery protocols and systems are not yet optimized. The aim of this study was to investigate the effects of injection parameters on solute distribution within the cervical subarachnoid space using a numerical platform. We developed a numerical model based on a patient-specific three dimensional geometry of the cervical subarachnoid space with idealized dorsal and ventral nerve roots and denticulate ligament anatomy. We considered the drug as massless particles within the flow field and with similar properties as the CSF, and we analyzed the effects of anatomy, catheter position, angle and injection flow rate on solute distribution within the cerebrospinal fluid by performing a series of numerical simulations. Results were compared quantitatively in terms of drug peak concentration, spread, accumulation rate and appearance instant over 15 seconds following the injection. Results indicated that solute distribution within the cervical spine was altered by all parameters investigated within the time range analyzed following the injection. The presence of spinal cord nerve roots and denticulate ligaments increased drug spread by 60% compared to simulations without these anatomical features. Catheter position and angle were both found to alter spread rate up to 86%, and catheter flow rate altered drug peak concentration up to 78%. The presented numerical platform fills a first gap towards the realization of a tool to parametrically assess and optimize intrathecal drug and gene vector delivery protocols and systems. Further investigation is needed to analyze drug spread over a longer clinically relevant time frame.