Implementation into the reynolds stress equation
 

Dear everyone!

The time evolution equation of the reynolds stress tensor R_ij = bar(u_i*u_j) is described as follows;

//-----------------------------------------------
tmp<fvSymmTensorMatrix> REqn
(
fvm::ddt(R_)
+ ...
//-----------------------------------------------

I would like to implement R_i1 = bar(u_i*u_1) and R_i2 = bar(u_i*u_2) into the equaion, how should I write?
(bar(u_i*u_1) = bar(u_i*u_j)*δ_j1 ? where δ_ij is the Kronecker delta)
Would you be able to help me?

aki