Let $C$ be a nonempty closed subset of $\mathbb R^n$. It is known that any such set satisfies the following condition
(Unique CPP a.e). For almost every $x \in \mathbb R^n$, there exists a unique point in $c(x) \in C$ such that $\|x-c(x)\| = \mbox{dist}(x,C) := \inf_{c \in C}\|x-c\|$.
For example, see Theorem 1 of http://tzamfirescu.tricube.de/TZamfirescu-110.pdf.
In particular, if $C$ is a closed convex subset of $\mathbb R^n$, then it is a folklore fact that $C$ verifies he unique CPP everywhere.
For any $\epsilon > 0$, let $C^\epsilon := \{x \in \mathbb R^n \mid d(x,A) \le \epsilon\}$ be the $\epsilon$-expansion of $C$.
Question. Under what minimalistic conditions on $C$ does there exist a function $\epsilon:(0,\infty) \to [0,\infty)$, such that $\limsup_{L \to \infty} \epsilon(L) = 0$ and the mapping $u_C:x \mapsto (x-c(x))/\|x-c(x)\|$ is $L$-Lipschitz a.e on $\mathbb R^n \setminus C^{\epsilon(L)}$ for any sufficiently large $L>0$ ?
I'm particularly interested in the case of sets of the form $C := \cap_{i \in I} C_i$ with $C_i := \{x \in \mathbb R^n \mid f_i(x) \le 0\}$, where $f_i:\mathbb R^n \to \mathbb R$ are sufficiently smooth functions. Even the case where $I$ is a singleton $C$ is the sublevel set of a sufficiently smooth function is already interesting.
The case where $C$ is a closed convex subset of $\mathbb R^n$ has been solved here https://mathoverflow.net/a/412676/78539: it suffices to take $\varepsilon(L) = 2/L$.