In a simple pde-problem, which I solve in a predefined mesh using NDSolveValue I would like to plot the result in the mesh NDSolve uses:
<< NDSolve`FEM`
netz = ToElementMesh[Rectangle[], "MeshElementType" -> TriangleElement, MaxCellMeasure -> {"Length" -> .1}]
u = NDSolveValue[{Laplacian[\[Psi][x, y], {x, y}] ==1 + \[Psi][x, y] + NeumannValue[0, x == 1] ,\[Psi][0, y] == 0, \[Psi][x, 0] == 0, \[Psi][x, 1] == 0}, \[Psi] ,Element[{x, y}, netz], Method -> "FiniteElement"]
If I plot the result u[x,y]
Plot3D[u[x, y], {x, 0, 1}, {y, 0, 1}, Mesh -> All]
I get a foursided mesh even though the solution mesh
u["ElementMesh"]["Wireframe"]
is triangular. How can I force Plot3D to show the right mesh?
ContourPlot
such asShow[{ContourPlot[u[x, y], {x, 0, 1}, {y, 0, 1}], netz["Wireframe"] }]
but this is not what you are asking about. $\endgroup$Plot3D[u[x, y] , Element[{x, y}, netz]]
would help $\endgroup$