How to access vertices of ASPECT mesh with deformation (e.g. free surface case)

[jbnaliboff]

Hi Zhikui,

Thank you for posting to the forum!

The standard ASPECT visualization post processing plugins all write out the deformed mesh at each time step.

One option would be to take the standard ASPECT output data files (HDF5, VTK, or ASCII) and process + rewrite that data in the OpenFOAM format.

Here is a post on how to load the VTK data into numpy arrays:

Loading VTK data into Python

Hi all, This topic has come up a number of times in recent conversations, so I thought I would share the general workflow for importing VTK data directly into python (numpy arrays). Below is a snippet of python code that should largely work out of the box for 3D cartesian models, where all one has to do is provide the full path and name of the .pvtu in a string. Please note that lines in bold are comments, and you would need to add a # symbol before them in python. I excluded them as they prod…

Would this approach work or is there a reason that this needs to be done within ASPECT, rather than as an external post processing step?

Cheers,
John