The stability of slopes consisting of unsaturated residual soils is strongly influenced by variations in local flow conditions, since in many cases much of their shear strength is conferred by the apparent cohesion caused by the suction phenomenon. With this in mind, the present study evaluated the possibility of creating a model where the final stability result is independent of the initial suction distribution through stability and flow analyzes performed in two and three-dimensional models. Flow and stability analyzes and modeling were performed using the SVDesigner, SVSlope and SVFlux-GE software components of the SVoffice package developed by SoilVision. It was tried to quantify the influence of the initial suction distribution in the cases studied. The soil parameters used as input data in the models were derived from studies developed by de Oliveira (2013) on the mechanisms of rupture of the mass movements that occurred in Nova Friburgo during the disaster that occurred in the Serrana Region of Rio de Janeiro in 2011, and the climatic data came from records in the region during and before the event. In the 2D models, it was possible, considering approximately eight months of climatic data, that different initial homogeneous suction conditions converged to a single stability result. It was also observed that the flow behavior differed between the 2D and 3D models, which had a significant impact on their stability analyzes.