When I add the option "space charge conservation" included into the module "semiconductor", and assign a specific domain (characterized by a insulating material) to this option, I can compute the whole study. It should work, but I do not why, it finally converges to a solution where D·n = 0 at the boundary of this insulating material with a highly doped semiconductor (to a metal grade), where D·n=0 should not be a suitable boundary condition.