Cauchy stress derived subfield poroelasticity example

Good afternoon everyone,

I want to have access to the Cauchy stress values (xx, yy, zz, xy) when I run step01_inflation example.

I’m reading through the PyLith documentation for poroelasticity and the user guide states the following:

This is what I have in the pylithapp.cfg file for both materials:

What else do I need to do in order to output the Cauchy stress values?

Best,
Eddie

I don’t believe you need to do anything else. What do you get in the HDF5 files? You can use h5dump -n FILENAME to see the structure.

Sorry for the delayed reply, we have been extra busy getting ready for the PyLith hackathon.

Hello Brad! I appreciate the response. My question is what do I need to do in order so that sigma_xx, sigma_yy, sigma_xy etc. appears in Paraview parameters section here:

Screenshot 2023-06-27 at 11.40.49 AM

The Cauchy stress is associated with a material rheology so it is only available in the material output files. You are looking at the output over the domain, which only includes the solution fields.