Resuelva PDEs axisimétricos 

Resuelva un PDE en una región axisimé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 un valor de Neumann de 100 W/m²(vatios por metro cuadrado) en el borde interior y una temperatura de 10 °C (degrees) en el borde exterior, con un coeficiente de conducción térmica k de 10 W/m/C (vatios por metro por grado).

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(*[W/m^2]*), r == 1], DirichletCondition[u[r, z] == 10, r == 2]}, u, {r, z} \[Element] \[CapitalOmega]];

Grafique la solución.

In[3]:=

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

Out[3]=

Verifique la solución.

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]=

Grafique la solución en 3D.

muestre la entrada completa de Wolfram Languageoculte la entrada

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]=