Next-generation propulsion systems will likely take advantage of boundary layer ingestion (BLI), which requires simulation and design tools to perform optimization at an affordable cost while taking airframe interaction into account. A solution for this problem is the reduced-order body force model (BFM). It allows for computationally inexpensive flow analysis, and when combined with the adjoint method accelerates conceptual design optimization. In this thesis, Hall’s inviscid body force model is implemented into the direct and adjoint solvers of the open-source software SU2. Direct solver runs of the BFM show results consistent with what is expected from theory for a range of flow characteristics. Furthermore, fourteen-fold reductions in grid size as well as one to two orders of magnitude reduction in computational time are observed while retaining flow behavior. In the adjoint solver, registration of source terms is successful, which is confirmed by changes in the adjoint flow over the body force domain. Physical interpretations of the adjoint vector agree with theory. Total gradient computation is attempted but unsuccessful. It is therefore recommended for future work, as showing this capability is a large step towards faster turbomachinery conceptual design optimization.