Let $B_+\subset\mathbb R^n$ be the intersection of the unit ball with the upper halfspace and $u$ be a harmonic function in $B_+$ that vanishes identically in the flat part of the boundary of $B_+$.
We can easily see that $u$ can be extended to $B$ by $-u$ in $B_-$ in a harmonic fashion by using the mean value property (Schwarz reflection principle). This allows us, for example, to use interior estimates near this part of the boundary.
I wonder whether there are other similar ideas or techniques when $u$ solves weakly a PDE like $$ \operatorname{div}(A(x)\nabla u(x))=0, \quad x\in B_+ $$ where $A$ is a symmetric, uniformly elliptic matrix with Lipschitz coefficients. To be specific, assume that $A$ is only defined on $\overline B_+$ (so that we can redefine $A$ on $B_-$ anyway we want but keeping the global Lipschitz condition on the coefficients) and $u$ also vanishes identically in the flat part of the boundary.
What I have tried: if the coefficients of $A$ have higher regularity I think that we can find a biLipschitz map $T$ from the upper half space to the lower half space and which is the identity on the hyperplane $\mathbb R^{n-1}$ such that $$ \operatorname{div}(A_T\nabla (u \circ T)) = 0, \quad x \in T^{-1}(B_+) $$ where $$ A_T = \vert JT \vert DT^{-1} A \circ T (DT^{-1})^t $$ and then we can join both domains $B_+$ and $T^{-1}(B_+)$ and both PDEs and solutions in a way that $A$ is still globally Lipschitz.