Objective To investigate the impact of different coil packing densities on intra-secular hemodynamics during embolization of intracranial aneurysm, so as to provide the theoretical evidence for endovascular treatment of intracranical aneurysms. Methods Based on digital subtraction angiography image data, the vessel numerical models with different packing densities were constructed using computational fluid dynamics software. The hemodynamic parameters, including the wall shear stress (WSS) and flow steamline and pressure, were obtained from the numerical simulation using Ansys CFX software. Results The aneurysm volume and low WSS area in the aneurysm sac were decreased after the first session of embolization compared with before treatment; the vortex in the aneurysm was decreased and the aneurysmal flow velocity had no significant changes. Meanwhile, the WSS at the aneurysm dome increased to 2.13 Pa. Continuous embolization greatly decreased the overall aneurysmal flow velocity and the aneurysmal flow velocity at the distal part of the residual aneurysm neck. The average wall pressure of the inflow tract had no noticeable changes before and after embolization. Conclusion The intra-saccular hemodynamics during embolization of intracranial aneurysm can be influenced by the coil packing density and coil distribution in the aneurysm.