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]