Computing variational bounds for flow through random aggregates of spheres