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

1 Answer 1

2
$\begingroup$

NeumannValue type boundary conditions may overlap, while DirichletCondition may not overlap. So something like this:

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] + 
      NeumannValue[100 T[t, x, y], x == 0.2 && y <= 0.025], 
    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]

enter image description here

Note the kink at y==0.025. Probably it would be good if the mesh respected this discontinuity.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.