“…Previous methods for tackling three dimensional Fickian diffusion have either relied on explicit finite difference calculation of uncomplicated geometries, which result in large simulation times due to the unstable nature of the discretisation method [6,15]. The alternating direction implicit finite difference method in three dimensions has also been used, but can suffer from oscillation and stability issues depending on the algorithm method used [31][32][33]. Boundary element simulation can be mathematically complex for this work, and finite element simulation can suffer from poor accuracy [7].…”