The parameters \(\tau _j\)  in the stabilization terms must be chosen in an appropriate way. With too small parameters, the effect is negligible, while with too large parameters, degrees of freedom can be constrained too strongly, potentially degenerating the accuracy of the approximation [1]}.
On Dirichlet boundaries, the required value to provide coercivity of the bilinear form correlates inversely with the Nitsche parameter.
On Neumann boundaries, coercivity of the bilinear form does not depend on the stabilization term.
Nevertheless, stabilization can still be applied to the Neumann boundary in order to repair the conditioning problems discussed in Section REF , for which a scaling with \(h^{2j + 1}\)  instead of \(h^{2j - 1}\)  in (REF ) suffices [2]}.
