Deformation banding is one of the most common failure modes in geomaterials such as rock, concrete, and soil. It is well known that appearance of bands of intense localized deformation significantly reduces the load-carrying capacity of any structure that develops them. Furthermore, when dealing with fluid-saturated geomaterials, the interplay between the contraction/dilation of pores and development of pore fluid pressures is expected to influence not only the strength of the solid matrix but also its ability to block or transport such fluids. Accurate and thorough simulation of these phenomena (i.e., deformation banding and fluid flow) requires numerical models capable of capturing fine-scale mechanical processes such as mineral particle rolling and sliding in granular soils and the coupling between porosity and relative permeability. Until recently, these processes could not even be observed in the laboratory. Numerical models could only interpret material behavior as a macroscopic process and were, therefore, unable to model the very complex behavior of saturated geomaterials accurately.