In this work, a known Eddy Current (EC) model is adapted to characterize subsurface defects in carbon fiber-reinforced polymer (CFRP) plates intended for the civil aerospace industry. The considered defects include delaminations, microcracks, porosity, fiber breakage, and the simultaneous presence of these defects. Each defect is modeled as an additive variation in the material’s electrical conductivity tensor, allowing for a detailed mathematical representation of the defect’s influence on the CFRP’s electromagnetic behavior. The additivity of the variations in the conductivity tensor is justified by the assumption that the defects are not visible to the naked eye, implying that the material does not require non-destructive testing. The adapted EC model admits a unique and stable solution by verifying that all analytical steps are satisfied. To reconstruct 2D maps of the magnetic flux density amplitude, a FEM formulation is adopted, based on the energy functional because it ensures a stable and consistent numerical formulation given its coercivity. Moreover, the numerical approach allows precise and reliable numerical solutions, enhancing the capability to detect and quantify defects. The numerical results show that the obtained 2D maps are entirely superimposable on those highlighting the distribution of mechanical stress states known in the literature, offering a clear advantage in terms of detection costs. This approach provides an effective and economical solution for the non-destructive inspection of CFRP, ensuring accurate and timely defect diagnosis for maintaining structural integrity.