“…Since ǫ ≪ 1, there is an inner region near each disk where the concentrations are at steady-state to leading order, and an outer region where the concentrations are timedependent to leading order. We will find the solutions in both regions using matched asymptotics in a similar manner to the single-disk problem descibed in our previous article [18]. In the inner region near disk j, the coordinate system is (r j , z, θ j ) as defined above, and variables will be C O, j (r j , z, θ j , t) and C R, j (r j , z, θ j , t).…”