The suffusion phenomenon occurs when fine soil particles are detached by seepage flow and transported away from the matrix. This process is one of the main causes of failure of hydraulic structures and road embankments. This study aimed to build a numerical model for simulating the suffusion within a porous medium. This model combines a flow law and an erosion equation related to the evolution of soil porosity. In addition, the dispersion and the deposition kinetics of eroded particles were combined with detachability process. The equations describe the evolution of the instantaneous concentration of the fluidized solid and the variation of eroded mass. Sensitivity analysis allows highlighting the influence of the different parameters on the suffusion, particularly that deposition kinetics starts acting only below a given hydraulic gradient and beyond a sample length. The simulation results indicate that the suffusion process is strongly related to hydraulic conditions, physical soil characteristics and pore water chemicals. The adjustment of numerical results with experimental data from laboratory tests provides a good agreement. This model is devoted to investigate conditions leading hydraulic works to avoid suffering suffusion.