Variability Estimation in Resistive Switching Devices, a Numerical and Kinetic Monte Carlo Perspective