Skip to content

DEShybrid leads to crash for single and mixed precision versions of OpenFOAM

Summary

OpenFOAM crashes if using the DEShybrid blending scheme with single (SP) or mixed precision (SPDP) compilation.

Steps to reproduce

Compile OpenFOAM in single or mixed precision mode. Run any case with DEShybrid blending scheme for div(phi,U).

Example case

pitzDaily with SA-DDES and DEShybrid blending scheme DEShybrid_pitzDaily.tar.gz

What is the current bug behaviour?

OpenFOAM crashes.

What is the expected correct behavior?

No crash :-)

Relevant logs and/or images

Starting time loop

fieldAverage fieldAverage1:
    Restarting averaging for fields:
        U: starting averaging at time 0
        p: starting averaging at time 0

Sampled surface:
    nearWall -> vtk

Time = 1e-05

Courant Number mean: 8.74894e-05 max: 0.0315889
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::pow(Foam::Field<float>&, Foam::UList<float> const&, float const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<float> const&) at ??:?
#5  Foam::tmp<Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> const&, float const&) at ??:?
#6  Foam::DEShybrid<Foam::Vector<float> >::calcBlendingFactor(Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<float, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#7  Foam::DEShybrid<Foam::Vector<float> >::blendingFactor(Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#8  Foam::DEShybrid<Foam::Vector<float> >::weights(Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#9  Foam::fv::gaussConvectionScheme<Foam::Vector<float> >::fvmDiv(Foam::GeometricField<float, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#10  Foam::fv::boundedConvectionScheme<Foam::Vector<float> >::fvmDiv(Foam::GeometricField<float, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<float>, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#11  ? at ??:?
#12  ? at ??:?
#13  __libc_start_main in /lib64/libc.so.6
#14  ? at /home/abuild/rpmbuild/BUILD/glibc-2.26/csu/../sysdeps/x86_64/start.S:122
Gleitkomma-Ausnahme (Speicherabzug geschrieben)

Environment information

  • OpenFOAM version : v1912
  • Operating system : openSUSE 15.1
  • Hardware info :
  • Compiler : gcc7

Possible fixes