Solve a pde on a 3d domain, also inside the domain.

Question

Hi,

I need to solve a system of PDE on a 3D domain but I need to see the solution also in the domain, not only on the surface. So I'd like to have a grid also in the inside of my domain.

it's possible only with FEniCS? What can i do to obtain that result?

Thanks
King regards
Valentina

Comments

Is it sufficient to save your solution to xdmf (or pvd) and visualize it in ParaView?

---

for example, if I write

mesh = BoxMesh (-2,-2,-2,2,2,2,10,10,10)

and I solve my PDE on this domain, both in pvd that in xdmf, I obtain the solution only on the surface of the cube, because it's "empty".
Also if I run

plot(mesh)

the box is empty, there's the mesh only of the surface.

---

But, in ParaView you can slice and clip the domain as you like?

---

No, I can't use slice filter.

---

That is very strange, I would look into that. Otherwise, what you can do is use cbcpost. Have a look at the documentation here, in particular section 3.1.2.2.

---

Generally, could I do what I described?
Thank you , I'm going to look it.

---

Yes, it's trivial in ParaView, and could be done using cbcpost.

---

It's strange but all 3D domain that I create are "Empty", also in paraview.

Sorry, how can I use cbcpost?
I done

To install simply run
python setup.py install --prefix=/path/to/your/favourite/site-packages

but i received a permission denied.

I've DOLFIN 1.5, is that a problem?

Thank you.

---

You probably need sudo.

---

great, thank you.