For Brownian motion on a convex polyhedral subset of a sphere or torus, the rate of convergence in distribution to uniformity is studied. The main result is a method to take a Markov coupling on the full sphere or torus and create a faster coupling on the convex polyhedral subset. Upper bounds on variation distance are computed, and applications are discussed.