solver does not converge without at least 1 Dirichlet BC