Valores de Neumann 

Compute o efeito de arrefecimento de um duto com um líquido de arrefecimento em uma seção transversal axissimétrica de um tubo.

Valores de Neumann prescrevem o fluxo ao longo da borda do contorno.

In[1]:=

X \[CapitalOmega] = ImplicitRegion[ y >= 0 && (x - 3/2)^2 + y^2 >= 1/25 && x^2 + y^2 >= 1 && x^2 + y^2 <= 4 && y <= x Tan[\[Pi]/8], {x, y}];

Visualize a região de simulação.

mostre o input completo de Wolfram Languageoculte o input

In[2]:=

X Show[Graphics[{Circle[{0, 0}, 2], Circle[{0, 0}, 1], Circle[#, 1/5] & /@ {{-1, -1}, {-1, 1}, {1, -1}, {1, 1}, {3/2, 0}, {-3/2, 0}, {0, 3/2}, {0, -3/2}}}], RegionPlot[\[CapitalOmega]], ImageSize -> 300]

Out[2]=

Especifique as condições de contorno de Dirichlet para definir as temperaturas no interior e exterior do tubo.

In[3]:=

X Subscript[\[CapitalGamma], D] = {DirichletCondition[u[x, y] == 200., x^2 + y^2 == 1], DirichletCondition[u[x, y] == 15., x^2 + y^2 == 4]};

Especifique um fluxo para dentro do duto de arrefecimento no tubo.

In[4]:=

X Subscript[\[CapitalGamma], N] = NeumannValue[-1000., (x - 3/2)^2 + y^2 == 1/25];

Resolva a equação de Laplace.

In[5]:=

X uif = NDSolveValue[{-10 \!\( \*SubsuperscriptBox[\(\[Del]\), \({x, y}\), \(2\)]\(u[x, y]\)\) == Subscript[\[CapitalGamma], N], Subscript[\[CapitalGamma], D]}, u, {x, y} \[Element] \[CapitalOmega]];

Projete graficamente a solução.

In[6]:=

X cp = ContourPlot[uif[x, y], {x, y} \[Element] \[CapitalOmega], Mesh -> None, ColorFunction -> "TemperatureMap", Contours -> 21, AspectRatio -> Automatic, Frame -> False, ImageSize -> 300]

Out[6]=

mostre o input completo de Wolfram Languageoculte o input

In[7]:=

X Show[{Graphics[{Circle[{0, 0}, 2], Circle[{0, 0}, 1], Circle[#, 1/5] & /@ {{-1, -1}, {-1, 1}, {1, -1}, {1, 1}, {3/2, 0}, {-3/2, 0}, {0, 3/2}, {0, -3/2}}}, PlotRange -> {{0.7, 2.2}, {-0.5, 1}}], cp}]

Out[7]=