[maka]

I'm trying to calculate a tensor (d/dx_k)(s_ik * u_j). s_ik is a second order strain rate tensor (symm(fvc::grad(Uprime))). I read the programmer manual, but I have a question:

can I use divergence operator on third order tensor (div(s*u)) and how it is defined (on which index the summation convention works? is it the first?).

There is another way to do it but less compact:

(d/dx_k)(s_ik * u_j) = ((d/dx_k) s_ik) * u_j + s_ik * ((d/dx_k) u_j);

in the programmers guide: div(s) = (d/dx_i) s_ij

This transform the relation to (s_ij is symmetric):

(d/dx_k)(s_ik * u_j) = div(s) * u + s * fvc::grad(u); Am I right? Thanks.

Best regards,
Maaka

[MarcusRinzsch]

Hi, I realize this is nearly 20 years later but guess what, there's still virtually no documentation that I can find.

My situation is a little different in that I am calculating the Divergence numerically in MatLAB of a third order tensor in cylindrical coordinates and in my search of what formula to follow (instead of deriving my own and potentially making mistakes) I found this!
Have you found a solution? I'd really appreciate any response or direction. I have been looking for so long...

-Marcus