When simulating complex flows, there are some physical situations that exhibit large fluctuations in particle density such as: planetary reentry, ablation due to arcing, rocket exhaust plumes, etc. When simulating these events, a high level of physical accuracy can be achieved with kinetic methods otherwise known as particle methods. However, this high level of physical accuracy requires large amounts of computation time. If the simulated flow is in collisional equilibrium, then less computationally intensive continuum methods, otherwise known as fluid methods, can be utilized. Hybrid Particle-Continuum (HPC) codes attempt to blend particle and fluid solutions in order to reduce computation time for transitional flows that exhibit both continuum and rarefied flow in a single domain. This thesis details the development of an HPC code in OpenFoam for Cal Poly's Aerospace Engineering department. The primary benchmark for the solver, named hybridFoam, was to simulate a 1D sod-shock simulation. This primary goal was achieved and a collection of test simulations were conducted to map out the solvers current capabilities and identify where future development efforts should focus.