Calculating averages with respect to probability measures on submanifolds (level-sets of a certain function) is crucial in computational statistical mechanics. Inspired by the idea of encoding constraints via stiff drift terms (so called soft constraints), in this talk I will present sampling schemes on submanifolds using various diffusion processes. This talk is based on joint work with Lara Neureither and Wei Zhang.