Consider Cahn-Hilliard (see this) equation hich is known as the $H^{-1}$ gradient flow of Cahn-Hilliard energy functional, also it is easy to verify that this equation is mass preserving i.e. measure of field variable is preserved as time proceed.

To preserve mass exactly in numerical method, i want to project the gradient in to mass preserving space, by euclidean $L^2$ projection (it works for me in practice), but I do not know that does this action correct or no, more specifically, if I do $L^2$ projection of a function which is member of $H^{-1}$, is the resulted function after the $L^2$ projection is still a member of $H^{-1}$.

the projection space is convex and can be defined as follows,

$\mathcal{A} := [ u\in X\ |\ 0 \leq u \leq 1, \ \int_\Omega u\ dx = \mathtt{constant} ]$

assume $X$ as suitable Banach space.

if question is not clear please let me know to clarify.

PS: $H^{-1}$ is a member of $L^2$, right?