bond_percolation

openpnm.topotools.bond_percolation(ij, occupied_bonds)[source]

Calculates the site and bond occupancy status for a bond percolation process given a list of occupied bonds.

Parameters
  • ij (array_like) – An N x 2 array of [site_A, site_B] connections. A site is considered occupied if any of it’s connecting bonds are occupied.

  • occupied_bonds (bool) – A list indicating whether a bond is occupied or not

Returns

  • A tuple containing a list of site and bond labels, indicating which

  • cluster each belongs to. A value of -1 indicates uninvaded.

Notes

The connected_components function of scipy.sparse.csgraph will give ALL sites a cluster number whether they are occupied or not, so this function essentially adjusts the cluster numbers to represent a percolation process.