Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 426
    • Issues 426
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 8
    • Merge requests 8
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #2082
Closed
Open
Issue created Apr 30, 2021 by Hendrik Hetmann@Zymorui

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

Assignee
Assign to
Time tracking