[tangd]

In the file kEpsilon.C, the source term for momentum equation is defined as:

- fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(fvc::grad(U)().T()))

But what does dev2() mean? What is the mathematical expression for the term dev2(fvc::grad(U)().T()))? Thanks!

[eugene]

grep dev2 OpenFOAM-1.3/src/OpenFOAM/lnInclude/*

edit OpenFOAM-1.3/src/OpenFOAM/lnInclude/TensorI.H

line 396-401

//- Return the deviatoric part of a tensor
template <class>
inline Tensor<cmpt> dev2(const Tensor<cmpt>& t)
{
return t - SphericalTensor<cmpt>::twoThirdsI*tr(t);
}

Hi,

I have a question conserning the operator "dev2". For example I am
interested in the following line:

fvc::div(muEff()*dev2(fvc::grad(U)().T()))

Can someone help me to translate it?

Thank you in andvance,
Rita