[ozgunoglu]

Hello everyone,
I have a problem while using a custom solver together with kOmegaSST model.

Summary of the problem:
When I run with kOmegaSST turb. model in remote HPC (OFv1912), it gives an error at the end of the first iteration. The error only appears with kOmegaSST model, for example, no error appears with kOmega model. Meanwhile, this problem does not arise when I run exactly the same case/setup on my personal workstation.

The problem in detail:
I have a customized OpenFOAM solver, and I am trying to make a simulation using “kOmegaSST” RANS turbulence model. I am using the code both in my personal workstation and university cluster.

Running the case on my personal workstation.

No problem arises for both compilation of the solver and running the case.

Running the case in remote HPC

No problem when I compile the solver at remote HPC, but the error below appeared at the end of the first iteration

--> FOAM FATAL ERROR:   

Cannot dereference nullptr at index 2 in range [0,8)  
    From function const T& Foam::UPtrList<T>::operator[](Foam::label) const [with T = Foam::fvPatchField<double>; Foam::label = int]  
    in file /apps/Vera/software/OpenFOAM/v1912-foss-2019b/OpenFOAM-v1912/src/OpenFOAM/lnInclude/UPtrListI.H at line 204.  

FOAM aborting  

#0  Foam::error::printStack(Foam::Ostream&) at ??:?  
#1  Foam::error::abort() at ??:?  
#2  ? at /apps/Vera/software/OpenFOAM/v1912-foss-2019b/OpenFOAM-v1912/src/OpenFOAM/lnInclude/errorManip.H:93 (discriminator 4)  
#3  ? at /apps/Vera/software/OpenFOAM/v1912-foss-2019b/OpenFOAM-v1912/src/OpenFOAM/lnInclude/UPtrListI.H:207  
#4  ? at /apps/Vera/software/OpenFOAM/v1912-foss-2019b/OpenFOAM-v1912/src/OpenFOAM/lnInclude/scalarFieldField.C:92 (discriminator 2)  
#5  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at /apps/Vera/software/OpenFOAM/v1912-foss-2019b/OpenFOAM-v1912/src/OpenFOAM/lnInclude/GeometricScalarField.C:116  
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?  
#7  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::DbnsCavTurbulenceModel<Foam::equilibriumCavitationModel> > > > > >::F1(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?  
#8  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::DbnsCavTurbulenceModel<Foam::equilibriumCavitationModel> > > > > >::correct() at ??:?  
#9  ? at ~/OpenFOAM/applications/dbnsCavitatingFoam/dbnsCavitatingFoam.C:136  
#10  __libc_start_main in /lib64/libc.so.6  
#11  ? at ??:?  
Aborted (core dumped)

• The error only appears with “kOmegaSST” model. I have tried other RANS models and did not encounter any problem; the case is running just fine.
• I have tested different geometries/meshes and different initial-boundary conditions to understand if the problem is case dependent, but the error again appeared.
• I compiled the solver for different OpenFOAM versions (v18 and v20) and the problem is persistent. Moreover, different compilation methods (compile from source code or using pre-compiled versions) did not change the result

Hope that all is clear and someone would like to help.

Thanks!