thank you so much, it worked now, however raised then another question, let say I use Plot3D[Evaluate[u[t,x,y] /. bvpdisk], ....] it should give a similar answer also, and for some reason it does not ,,, could you please help me to understand the reason?