Resolva EDPs axissimétricas 

Resolva uma EDP sobre uma região axissimétrica.

In[1]:=

X Subscript[r, 1] = 1; Subscript[r, 2] = 2; \[CapitalOmega] = ImplicitRegion[ True, {{r, Subscript[r, 1], Subscript[r, 2]}, {z, 0, 1/2}}];

Especifique um valor de Neumann de 100 W/m² (watts por metro quadrado) na borda interna e uma temperatura de 10 °C na borda externa, com um coeficiente de condução térmica k de 10 W/m/C (watts por metro por grau).

In[2]:=

X k = 100; nv = -100; dc = 10; uif = NDSolveValue[{k/r D[ r D[u[r, z], r], r] + k D[u[r, z], {z, 2}] == NeumannValue[nv, r == 1], DirichletCondition[u[r, z] == 10, r == 2]}, u, {r, z} \[Element] \[CapitalOmega]];

Projete graficamente a solução.

In[3]:=

X dp = DensityPlot[uif[r, z], {r, z} \[Element] \[CapitalOmega], ColorFunction -> "Temperature"]

Out[3]=

Verifique a solução.

In[4]:=

X analyticalSolution[r_, z_] = (dc + nv*Subscript[r, 1]/k Log[r/Subscript[r, 2]]); Plot[uif[r, 0] - analyticalSolution[r, z], {r, Subscript[r, 1], Subscript[r, 2]}]

Out[4]=

Projete a solução em 3D.

mostre o input completo de Wolfram Languageoculte o input

In[5]:=

X Show[ RegionPlot3D[ 1 <= x^2 + y^2 <= 4 && 0 <= z <= 1/2, {x, -2, 2}, {y, -2, 2}, {z, 0, 1/2}, Boxed -> False, Axes -> False, PlotPoints -> 40, PlotStyle -> {Opacity[0.2]}, Mesh -> False], Graphics3D[{EdgeForm[Red], FaceForm[Gray], GraphicsComplex[{{1, 0, 0}, {2, 0, 0}, {2, 0, 1/2}, {1, 0, 1/2}}, {Texture[ Show[dp, Frame -> False, PlotRangePadding -> None]], Lighting -> {{"Ambient", White}}, Polygon[{{1, 2, 3, 4}}, VertexTextureCoordinates -> {{{0, 0}, {1, 0}, {1, 1}, {0, 1}}}]}]}] ]

Out[5]=