The stress in static granular media under gravity

A fundamental open problem in the mechanics of granular media is the determination of the stress in the static state. It is known that the static stress depends strongly on how the grain assembly is created and the nature of confining boundaries. Non-trivial spatial variations have been observed even in simple geometries, posing long-standing challenges to continuum modelling. In this paper, we create gravity-deposited grain packings computationally and devise a method to visualise the paths of load transmission, which we call force lines. We show that the force lines reflect the flow during deposition, thereby encoding preparation history. We then show that the force lines coincide with ensemble averaged biased random walks in the particle contact network; this identification yields a closure relation for the stress, which together with the static momentum balances fully determines the stress field. The model makes accurate predictions for the stress in piles and silos, even for unusual deposition methods, thereby showing promise for more general scenarios.