nonlinear solver very unstable
mixed dirichlet bc not clear