Over 10% of the global populatin does not have access to high quality drinking water caused by the introduction of contaminants into aquifers. Colloids, collections of nanoparticles, have been explored as a possible solution to remediate polluted groundwater. To maximize the effectiveness and efficiency of this process, mathematical models have been developed to exambine the movement of the colloids. This research study aims to conduct analysis and implement more efficient algorithms for simulations of a previously developed model. The governing differential equation system consists of the MNM1D (Micro- and Nanoparticle transport Model) in porous media in one-dimensional geometry coupled with the Darcy-Forchheimer fluid model modified to include elctromagnetic effects. The resulting model accounts for attachment and detachment phenomena, both of the linear and Langmuirian type, as well as changes to hydrochemical parameters such as maximum colloidal particle concentration in the porous medium. In this work, the computational efficiency of the schemes was improved by implementing a parallel infrastructure for benchmarck problems with a flexible algorithm that is efficient, robust, and stable. These improvements in the reliability and efficiency of simulations of nanoparticle transport in porous media will contribute to the creation of an efficient method to counteract contaminants in groundwater and ultimately increase availability of clean drinking water.