[lhcamilo]

Hi everyone,

I have been working on a solver for a reynolds stresses transport equation, which is aimed at solving each component of the said equation.

one of these components is the divergence of a third order tensor, namely: d_k( u_k u_i u_j).

My attempts to declare this third order tensor have resulted in compiling errors, which lead me to search around for implemetations of third order tensors, but I could not find that.

In fact, after googling the issue a bit I came across a 2010 post iwhich suggested that OF did not have support for such tensors, but rather each component would have to be coded in a list of scalars.

I have sort of formulated a work around but before I went there, I committed to that I was wondering if this had already been implemented in OF in some way. (Note, I had been using OF 2.0.x but I am currently compiling 2.3.x ...)

For my work around, how would I go about taking the derivative in a single direction? My work around would involve doing something like :


d_x {u_x (u_i u_j)} + d_y {u_y (u_i u_j)} + d_z {u_z (u_i u_j)}

Thanks in advance for your help