[t.teschner]

I am trying to get the DeardorffDiffStress model to work and can't seem to get the compressbile version to run. Running the model in incompressible mode works just fine, but when switching to compressible flow, I am getting the following error:

Incompatible dimensions for operation
    [R[1 -1 -3 0 0 0 0] ] - [R[0 2 -3 0 0 0 0] ]

Makes sort of sense, I am using

[0 2 -2 0 0 0 0]

as my dimensions when setting up the incompressible version for the Reynolds stresses, so multiplying by density should then remove this issue, so setting R to

[1 -1 -2 0 0 0 0]

then results in this error:

--> FOAM FATAL ERROR: (openfoam-2212)
Different dimensions for '(a = b)'
     dimensions : [0 2 -1 0 0 0 0] != [0.5 0.5 -1 0 0 0 0]


    From bool Foam::checkDims(const char*, const Foam::dimensionSet&, const Foam::dimensionSet&)
    in file dimensionSet/dimensionSet.C at line 56.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#1  Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2  Foam::error::exiting(int, bool) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#3  Foam::dimensionSet::operator=(Foam::dimensionSet const&) const in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator=(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam
#5  Foam::LESModels::DeardorffDiffStress<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut() in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#6  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam
#7  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam

So at least I know the issue is in correctNut() but that does not really help. The solver does not start (i.e. no momentum equation is solved, it crashes before then, so probably the issue is happening once the laplacian term is evaluating the viscosity), just wondering where these dimensions of

[0 2 -1 0 0 0 0] != [0.5 0.5 -1 0 0 0 0]

do come from. The first set seems to be something of a kinematic viscosity, but no idea where the second one is comming from (or what operation is producing kg^0.5 * m^0.5 s^-1 ... seems to be equivalent of sqrt(R)).


Any idea what to change to get this to work for compressible flows?