8
$\begingroup$

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 on boundaries of $y$. I want to put Neumann Values. Can I do that?

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]
$\endgroup$
2
  • 6
    $\begingroup$ Something like NeumannValue[100, x == 0.2] + NeumannValue[200, y == 0.]? $\endgroup$ – Michael E2 Nov 29 '15 at 22:03
  • $\begingroup$ Exactly!! That was so easy (It´s my first time using Mathematica), Thank you very much, Greetings! $\endgroup$ – Edgar Olivos Nov 29 '15 at 22:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.