[DataDriven]

Hi!

I'm working on a customized turbulence model in OF7. Somewhere in the calculations I have to add volSymmTensorField and volTensorField which I can't figure out.

On the forum I've seen one post about this, where it was suggested to multiply the volSymmTensorfield with tensor(1, 0, 0, 0, 1, 0, 0, 0, 1) if the computational overhead was not that important, which it isn't since they are very small 2D cases. See below:

  tauij_ = -2*nut() * S & tensor(1, 0, 0, 0, 1, 0, 0, 0, 1);

This is only part of the expression, later another volTensorField must be added, but this is already not allowed because S is symmetric and tauij is not.



The error message is:

 no match for ‘operator=’ (operand types are ‘Foam::volTensorField 
 {aka Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, 

Foam::volMesh>}’ and‘Foam::tmp<Foam::DimensionedField<Foam::SymmTensor<double>, 

Foam::volMesh> >’)

       tauij_ = -2*nut() * S & I

Another solution was to define a volTensorField and then set each individual component from the volSymmTensorField as below. But that's ugly and also I'm unsure whether indexing starts at 1 (I think the programmersGuide implies it starts at 1).

S.component(1) =  S_symm.component(1);
S.component(2) =  S_symm.component(2);
S.component(3) =  S_symm.component(3);
   
S.component(4) =  S_symm.component(2);
S.component(5) =  S_symm.component(4);
S.component(6) =  S_symm.component(5);
  
S.component(7) =  S_symm.component(3);
S.component(8) =  S_symm.component(5);
S.component(9) =  S_symm.component(6);

Can anybody help me with the proper transformation of a volSymmTensorField to volTensorField? Thank you in advance.



PS: these are the same right?

I = tensor(1, 0, 0, 0, 1, 0, 0, 0, 1)

[agustinvo]

Hi,


could you provide the equation you want to add in OpenFOAM? Maybe it is possible to find a way to implement it.