  The discrete element method constitutes a general class of modeling
techniques to simulate the microscopic behavior (i.e. at the particle scale) of
granular/soil materials. We present a contact dynamics method, accounting for
the cohesive nature of fine powders and soils. A modification of the model
adjusted to capture the essential physical processes underlying the dynamics of
generation and collapse of loose systems is able to simulate "quicksand"
behavior of a collapsing soil material, in particular of a specific type, which
we call "living quicksand". We investigate the penetration behavior of an
object for varying density of the material. We also investigate the dynamics of
the penetration process, by measuring the relation between the driving force
and the resulting velocity of the intruder, leading to a "power law" behavior
with exponent 1/2, i.e. a quadratic velocity dependence of the drag force on
the intruder.
