Mack (@sac-bsa) and I are trying to track down a singularity in the viscoelastic bending beam benchmark. It is easy to see in the stress_xx field at the end time and at other time steps (e.g., step 36).
The benchmark is on a grid with square cells of side h.
ParaView’s rendering appears to make one cell into 4 square cells of side h/2. The singularity appears as a triangle in the upper left hand corner of 1 cell with vertices (0, h/2) (0, h), and (h/2, h). The actual location of the singularity appears to be in the cell that is directly in the center of the computational domain.
Where is the data that ASPECT outputs to the *.vtu file located?
On the corners of the cell?
At the midpoint of the cell?
At the centers of the cell edges?
Some combination of the above?