I want to solve numerically the heat exchange of a fin in transient state and two-dimentional, here is my attempt but instead of put DirichletCondition in Boundaries of "y" I want to put Neumann Values, I can do that?

Thanks for your help!!

Clear["Global`*"]
tf = 0.02;
A = 90 10^-2;
Tc = 10;
Tf = 5;
S1 = NDSolve[{A (D[T[t, x, y], x, x] + D[T[t, x, y], y, y]) - 
     D[T[t, x, y], t] == NeumannValue[100, x == 0.2], 
   DirichletCondition[T[t, x, y] == Tc, x == 0], 
   DirichletCondition[T[t, x, y] == Tf, y == 0], 
   DirichletCondition[T[t, x, y] == Tf, y == 0.05], T[0, x, y] == Tf},
   T, {t, 0, tf}, {x, 0, 0.2}, {y, 0, 0.05}]
G1 = Plot3D[Evaluate[T[tf, x, y] /. S1], {x, 0, 0.2}, {y, 0, 0.05}, 
  PlotRange -> All, AxesLabel -> Automatic]