This paper studies the optimal stopping problem in the presence of model uncertainty (ambiguity). We develop a numerically implementable method to solve this problem in a general setting, allowing for general time-consistent ambiguity-averse preferences and general payoff processes driven by jump diffusions. Our method consists of three steps. First, we construct a suitable Doob martingale associated with the solution to the optimal stopping problem using backward stochastic calculus. Second, we employ this martingale to construct an approximated upper bound to the solution using duality. Third, we introduce backward-forward simulation to obtain a genuine upper bound to the solution, which converges to the true solution asymptotically. We a...