“…Notable examples include the solutions of Obdam and Veling [1987], Janković and Barnes [1999a], and Barnes and Janković [1999]. A number of special closed-form solutions, including flow through circular and elliptical inclusions, are documented in Strack [1989] and there are a number of one-dimensional solutions available for both steady state and transient flow in piecewise-homogeneous domains [e.g., Yeh and Kuo, 2010;Liang and Zhang, 2013;Hu and Jiao, 2014]. Simulation of flow through highly heterogeneous conductivity fields using analytical approaches is possible via the superposition of solutions for piecewise constant heterogeneity [e.g., Janković et al, 2003]; this approach is very successful, enabling the solution to problems with tens of thousands of inclusions.…”