differential equations – Problem with “does not evaluate to a numeric scalar at the coordinate {2.75,3/2}, but it is outside domain


enter image description here

Problem with “does not evaluate to a numeric scalar at the coordinate {2.75,3/2}, but that point is outside domain. How can I fix it??

Remove("Global`*") ;
Needs("NDSolve`FEM`")
HeatTransferModelAxisymmetric(T_, {r_, z_}, k_, ρ_, Cp_,Velocity_, Source_) := Module({V, Q}, V = If(Velocity === "NoFlow", 0, Velocity.Inactive(Grad)(T, {r, z})); Q = If(Source === "NoSource", 0, Source);(1 - (r)^2 - ((1 - 0.5^2)/(Log(1/0.5)))*Log(1/r))*D(T, z) +1/r*D(-k*r*D(T, r), r) + D(-k*D(T, z), z) + V - Q)
op = HeatTransferModelAxisymmetric(T(r, z), {r, z}, k, ρ, Cp,"NoFlow", "NoSource")
parameters = {k -> 10, d -> 10};
Subscript(Γ, flux) = NeumannValue(40*(500 - T(r, z)), r == 2.5);
Subscript(Γ, temp) = DirichletCondition(T(r, z) == 1200, r == 1);
Subscript(Γ, enter) = DirichletCondition(T(r, z) == 800, z == 0);
Ω = Rectangle({1, 0}, {2.5, 3});
pde = {op == Subscript(Γ, flux),Subscript(Γ, temp),Subscript(Γ, enter)} /. parameters; Tfun =NDSolveValue(pde, T, {r, z} ∈ Ω);
MassTransferModelAxisymmetric(c_, {r_, z_}, d_, Velocity_, Source_) :=Module({V, Q},V = If(Velocity === "NoFlow", 0, Velocity.Inactive(Grad)(c, {r, z}));Q = If(Source === "NoSource", 0, Source);(1 - (r)^2 - ((1 - 0.5^2)/(Log(1/0.5)))*Log(1/r))*D(c, z) + 1/r*D(-d*r*D(c, r), r) + D(-d*D(c, z), z) + 4.67`*^14*Exp(-36635/(Tfun(r, z)))*c + V - Q)
op2 = MassTransferModelAxisymmetric(c(r, z), {r, z}, d, "NoFlow","NoSource")
Subscript(Γ, enter) = DirichletCondition(c(r, z) == 800, z == 0);
pde2 = {op2 == Subscript(Γ, enter)} /.parameters;cfun=NDSolveValue(pde2, c, {r, z} ∈ Ω);