Fail to reproduced Decay of homogeneous isotropic turbulence

[dimaZloy]

Hi!

I try to reproduce a simple test case decay of homogeneous isotropic turbulence from std OF / pisoFoam using QHDFoam.
In practice, I just updated files according to some tutorials from QHDFoam - its crashed from the beginning with a message that something is not implemented (it might be cyclic boundary conditions).

Also, I tried to run a simple 3d case with QHDFoam, which is crashed without any messages as well - it might be something wrong with the solver for 3d cases.

Thanks and regards,
boxIsoTurb-QHDFoam.zip

[mkraposhin]

Hey, the problem most likely is with the "cyclic" boundary conditions. The QGD/QHD approach extensively uses derivatives in tangential directions to calculate fluxes. The approximations of these derivatives are not available in OpenFOAM, therefore, we added them manually (fvsc namespace). The implementation is able to work with internal faces, faces on inter-processor boundaries, external faces with Dirichlet and Neuman B.C.

It looks like that for the cyclic B.C. the implementation needs to be updated.