[agustinvo]

Hello,


during my simulations I have to add an extra term in the Reynolds stress tensor R. This new term is only applicable when using eddy viscosity models, so I was planning to add it using the fvOptions.correct() function to enable this modification whenever is needed.


As far as I found, I can import R in my code using:

const volSymmTensorField& R = turbModel.R();

but obviously later on I have to correct it, therefore I can't call it using const. How can I copy this R field and apply on it an fvOptions properly?.


Regards,
Agustín

[HPE]

Out of my head:

tmp<volSymmTensorField> tR = turbModel.R();
auto& R = tR.ref();

There might be something similar in one of those fvOptions: link. Have a look, please.

[agustinvo]

Thanks! Now it works.


EDIT: well, it works as far as the turbulence models belongs to the eddyViscosity class. If I try to use a Reynolds stress model, I get the next error message:

Attempt to acquire non-const reference to const object from a tmp<N4Foam14GeometricFieldINS_10SymmTensorIdEENS_12fvPatchFieldENS_7volMeshEEE> 

    From function T& Foam::tmp<T>::ref() const [with T = Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>] 
    in file /home/user/OpenFOAM/OpenFOAM-7/src/OpenFOAM/lnInclude/tmpI.H at line 187.

I suppose it happens because R is now a model variable and not a model output. The code should be able to take the R value (it doesn't matter its origin) and apply, if needed, the fvOption. I will check this again.




EDIT2: I managed to implement this feature. The Reynolds stress needed to be called like:

tmp<volSymmTensorField> tR = turbModel.R();
volSymmTensorField& R = const_cast<volSymmTensorField&>(tR());

and it runs now. The reference couldn't be constant, otherwise you cannot modify the R variable. Using this code I'm able to modify the R field if needed.