  Dual three-dimensional networks of structural and transport elements were
combined to model the effect of fracture on mass transport in quasi-brittle
geomaterials. Element connectivity of the structural network, representing
elasticity and fracture, was defined by the Delaunay tessellation of a random
set of points. The connectivity of transport elements within the transport
network was defined by the Voronoi tessellation of the same set of points. A
new discretisation strategy for domain boundaries was developed to apply
boundary conditions for the coupled analyses. The properties of transport
elements were chosen to evolve with the crack opening values of neighbouring
structural elements. Through benchmark comparisons involving non-stationary
transport and fracture, the proposed dual network approach was shown to be
objective with respect to element size and orientation.
