[petsc-users] Visualizing higher order finite element output in ParaView

[Jed Brown]

Anna Dalklint <anna.dalklint at solid.lth.se> writes:

> I want to save e.g. the discretized displacement field obtained from a quasi-static non-linear finite element simulation using 10 node tetrahedral elements (i.e. which has edge dofs). As mentioned, I use PetscSection to add the additional dofs on edges. I have also written my own Newton solver, i.e. I do not use SNES. In conclusion, what I want is to be able to save the discretized displacement field in each outer iteration of the Newton loop (where I increase the pseudo-time, i.e. scaling of the load). I would then preferably be able to load a stack of these files (call them u001, u002, u003… for each “load-step”) and step in “time” in ParaView.

Please use DMSetOutputSequenceNumber to record step number. You can either use one PetscViewer of type CGNS and call VecView in your loading loop or you can write a sequence of files by creating a new PetscViewer each time.