[QMW]

Hello everyone!
I am a beginner of OpenFoam.In the modeling of the generation term of the turbulent kinetic energy equation in the turbulence model, I want to introduce a velocity divergence, so that the turbulent kinetic energy will increase sharply in theory when the over-excited wave flow ( because the severe compression of the shock wave will lead to a large velocity divergence ). However, I did not find this phenomenon in the verification, and I was very confused about this. At a moment, it was my code and the distribution nephogram of the turbulent kinetic energy. The model was Hifire-5, with a flow of 8 Mach and a shock wave at the head.

volScalarField::Internal G
(
this->GName(),

nut*
sqrt(2.0)*mag(symm(fvc::grad(this->U_)))*
sqrt(2.0*magSqr(skew(fvc::grad(this->U_))))
-(2/3)*(this->k_)*(fvc::div(this->U_))

);

[QMW]

Here is the picture of k distribution