I am looking for a reference for the following result (if it is true, which I would expect):


Let $\Omega\subset\mathbb{R}^n$ be a bounded Lipschitz domain. Let $\Gamma\_0\subset\partial\Omega$ be sufficiently regular.


Let $V\_1:=\left\{\phi\in C^\infty\left(\overline{\Omega}\right):\phi\geq 0 ~\text{on}~ \Gamma\_0\right\}$.


Let $V\_2$ be the set of functions $u\in W^{1,p}(\Omega)$ with non-negative trace on $\Gamma\_0$.


Then $V\_1$ is dense in $V\_2$.


If we consider the case of "$\phi = 0$" instead of "$\phi \geq 0$", of course an analogous result for $\Gamma\_0=\partial\Omega$ is well known, and there is a result that it stays true if $\Gamma\_0\subset\partial\Omega$ is relatively open and has relative Lipschitz boundary.


It would already be very helpful to have a reference for the result with $\phi \geq 0$ and $\Gamma\_0=\partial\Omega$.



Density of smooth functions in Sobolev space, respecting nonnegative traces
More of a comment than an answer, but...


I don't know a reference, but perhaps you could prove it as follows. For $u \in W^{1,p}(\Omega)$ with nonnegative trace, write $u = u^+ - u^-$ which are both in $W^{1,p}(\Omega)$. Since $u^+$ is nonnegative everywhere, you should be able to approximate it by nonnegative smooth functions $\phi\_n$. And since $u^-$ has zero trace, by the result you quote, you can approximate it by smooth functions $\psi\_n$ vanishing on the boundary (maybe even compactly supported). Then $\phi\_n - \psi\_n$ are nonnegative on the boundary and converge to $u$.


