Strain rate partitioning

[ArneSpang]

Hi everyone,
compute_τij looks very useful. Is it possible to also make it return the strain rate of the individual rheological components? This would be important for shear heating and/or grain size reduction. Alternatively, returning effective viscosities of the individual components would also work. I had a look but at this point, I cannot follow the source code anymore.

[albert-de-montserrat]

Rather than modifying compute_τij, once you know the stress tensor/invariant after running compute_τij, you could add a function similar to compute_εII but using ntuple instead of nreduce, that should give you a tuple with the partitioned strain rate:

compute_elements_εII(v::CompositeRheology{T,N}, τII, args) where {T, N} = ntuple(i -> first(compute_εII( v.elements[i], τII, args)), Val(N))

[ArneSpang]

Thank you for the quick response! I will try this out

[boriskaus]

can this be closed?

[ArneSpang]

Yeah, the function was included around the last hackathon