[<lorentz<]

Hello everyone,

I am trying to run a simulation with chtMultiRegionSimpleFoam using the Ellipting Blending Reynolds Stress Model (EBRSM). I am using OpenFOAM-v2306.
I started from the case in the tutorial section (EBRSM turbulence model with simpleFoam solver: $FOAM_TUTORIALS/verificationAndValidation/turbulenceModels/planeChannel) and I modified this case to solve it with chtMultiRegionSimpleFoam. The geometry is a rectangular domain with a volumetric momentum source, the BC are two walls at the same temperature in the upper and lower sides and a cyclic condition in the left and right side.
When I run the simulation with chtMultiRegionSimpleFoam it crashes during the first time step. here's the output:

Time = 1


Solving for fluid region region1
smoothSolver:  Solving for Ux, Initial residual = 0.9999999966456584, Final residual = 1.092316707381895e-11, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 3.15409836055221e-13, No Iterations 11
Pressure gradient source: uncorrected Ubar = 17.54514880973776, pressure gradient = 0.04728518145885395
smoothSolver:  Solving for h, Initial residual = 0.9999999417030885, Final residual = 5.504674013964634e-08, No Iterations 1000
Min/max T:373.1499998029909 373.1500011392363
GAMG:  Solving for p_rgh, Initial residual = 0.8116908744296613, Final residual = 9.470812287532965e-07, No Iterations 16
GAMG:  Solving for p_rgh, Initial residual = 5.8644174102195e-07, Final residual = 6.048763452135292e-08, No Iterations 3
Pressure gradient source: uncorrected Ubar = 17.5451678969471, pressure gradient = 0.02837495316405573
time step continuity errors : sum local = 1.822343471660741e-08, global = -7.188472161815702e-17, cumulative = -7.188472161815702e-17
Min/max rho:1 1
DICPBiCGStab:  Solving for f, Initial residual = 1, Final residual = 6.806232488646261e-20, No Iterations 3


--> FOAM FATAL ERROR: (openfoam-2306)
Incompatible dimensions for operation
    [epsilon[1 -1 -4 0 0 0 0] ] + [((Cstability*epsilon)|k)[0 2 -4 0 0 0 0] ]

    From void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::DimensionedField<Type, Foam::volMesh>&, const char*) [with Type = double]
    in file /usr/src/packages/BUILD/src/finiteVolume/lnInclude/fvMatrix.C at line 1897.

FOAM aborting

[stack trace]
=============
#1  Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2  Foam::error::exiting(int, bool) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#3  void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::DimensionedField<double, Foam::volMesh> const&, char const*) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
#4  Foam::tmp<Foam::fvMatrix<double> > Foam::operator+<double>(Foam::tmp<Foam::fvMatrix<double> > const&, Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#5  Foam::RASModels::EBRSM<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#6  ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
#7  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8  ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
=============
Aborted

It seems that there are incompatible dimensions while solving the epsilon equation.
I have already tried to switch to another turbulence model (realizableKE) and with that the simulation runs perfectly.
I've also tried to run another case with EBRSM turbulence using rhoSimpleFoam as a solver and the error is the same.
I am struggling to find out how to solve this problem, maybe the EBRSM turbulence model works only with incompressible solvers?
I have attached the case, any help is appreciated!

Regards

Lorenzo

[nnprashant]

Hi Lorenzo,
Were you able to find the root cause of the error. I have been trying to implement the same turbulence model and end up having the same error. The EBRSM model is compatible with compressible and incompressible flows.

Regards
Nagaraj