Hard carbon (HC) has been predominantly used as a typical anode material of sodium-ion batteries (SIBs) but its sodiation mechanism has been debated. In this work, we investigate the adsorption of Na atoms on defective graphene under propylene carbonate (PC) and water solvent as well as vacuum conditions to clarify the sodiation mechanism of HC. Within the joint density functional theory framework, we use the nonlinear polarizable continuum model for PC and the charge-asymmetric nonlocally-determined local electric solvation model for water. Our calculations reveal that the centre of each point defect such as mono-vacancy (MV), di-vacancy (DV) and Stone-Wales is a preferable adsorption site and the electrolyte enhances the Na adsorption through implicit interaction. Furthermore, we calculate the formation energies of multiple Na atom arrangements on the defective graphene and estimate the electrode potential versus Na/Na+, verifying that the multiple Na adsorption on the MV and DV defective graphene under the PC electrolyte conditions is related to the slope region of the discharge curve in HC. This reveals new prospects for optimizing anodes and electrolytes for high performance SIBs.