Let $\mathbb{D}^2$ be the closed $2$-dimensional unit disk, and let $g:\mathbb{D}^2 \to \mathbb{R}$ be a non-constant harmonic function (in particular smooth up to the boundary).
Does there exist a sequence of smooth harmonic functions $g_n$ on $\mathbb{D}^2 $, such that $g_n \to g$ in $W^{1,2}$ and $dg_n \neq 0$ everywhere on $ \text{int}(\mathbb{D}^2)$?
Since we can add additive constants to the $g_n$, we can arrange $\int_{\mathbb D^2} g_n=\int_{\mathbb D^2} g$, so the $W^{1,2}$ convergence of the $g_n$ is essentially equivalent to $dg_n \to dg$ in $L^2$. (via Poincare inequality).
Thinking on $dg$ as a vector field, I think that we can always approximate it with a non-zero vector field in $L^2$. However, the only procedure I know for doing that does not produce approximating vector fields which are gradients of harmonic functions (or gradients of anything, really).