[CIG-SHORT] Stress discontinuity at fault intersection in 3D

[Brad Aagaard]

Josimar,

You can output the stress tensor in the domain and use post-processing 
to compute the stress changes at whatever points you want and then 
resolve them into tractions for a specific orientation. We don't have 
utilities to do this, so you would need to write some (Python) code to 
do this post-processing.

Regards,
Brad

On 2/23/18 11:35 AM, Josimar Alves da Silva wrote:
> Hi Brad,
>
> I have a question regarding fault intersection in 3D. This is a follow
> up question from a previous recent discussion.
>
> When I look at the normal stress at one of my faults I notice a
> discontinuity in the stresses at the fault intersection. See attached
> figure for details.
>
> This is of course a result of the fault_edges definition as you
> explained in a recent post. My questions is:
>
> i) Is it possible to output normal and shear stresses on a plane inside
> the domain without having to define those planes as faults in PyLith ?
> This would perhaps bypass the fault intersection problem for cases where
> one only wants stresses on the fault. The goal is to avoid having to
> define fault buried edges that result in the attached stress discontinuity.
>
> thank you in advance,
> Josimar
>
>
>