Retrofitting seismically deficient beam-column joints (BCJs) in reinforced concrete (RC) structures is crucial to preventing collapse during high-intensity earthquakes. Fiber reinforced polymer (FRP) composites have proven effective for retrofitting these components, but debonding of the FRP from concrete surfaces remains a challenge. Recently, the grooving method (GM) has emerged as a promising solution to mitigate this issue, enhancing the bond between FRP and concrete. This paper presents a numerical investigation of BCJs retrofitted with FRP composites using the GM, focusing on various design parameters. The study utilized finite element models developed in ABAQUS, incorporating the concrete damage plasticity (CDP) model to simulate the nonlinear behavior of concrete. The interface between the concrete and FRP was modeled as a perfect bond, simulating the strong adhesion provided by the GM. The numerical results were validated against experimental data from two BCJ specimens, showing good agreement in terms of load-displacement behavior, peak loads, and failure modes. Key parameters studied include the beam longitudinal reinforcement ratio, column axial load ratio, and joint aspect ratio. The findings reveal that these parameters significantly affect the shear capacity of retrofitted joints, impacting the efficiency of the retrofitting.