Skip to content
GitLab
Closed
Issue created by Admin@OpenFOAM-adminMaintainer

reconstructPar in case using Turbulent DFSEM

I have created a test case using the new DFSEM turbulence synthesis BC. The solver runs. But when I try to reconstruct from the parallel run, the following floating point exception occurs:

Starting program: /home/hk354/OpenFOAM/OpenFOAM-plus/platforms/linux64GccDPInt32Opt/bin/reconstructPar -latestTime -fields \(k\ nut\ pMean\  pPrime2Mean\ U\ UMean\ UPrime2Mean\)
[Thread debugging using libthread_db enabled]
Using host libthread_db library "/lib/x86_64-linux-gnu/libthread_db.so.1".
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  plus                                  |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : plus-c9ef631fb806
Exec   : /home/hk354/OpenFOAM/OpenFOAM-plus/platforms/linux64GccDPInt32Opt/bin/reconstructPar -latestTime -fields (k nut pMean pPrime2Mean U UMean UPrime2Mean)
Date   : Jul 26 2016
Time   : 07:23:35
Host   : "einsiedel"
PID    : 15399
Case   : /home/hk354/SCRATCH/channel395DFSEM_hkmesh
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947

--> FOAM Warning : 
    From function void Foam::Time::setControls()
    in file db/Time/Time.C at line 230
    Increasing the timePrecision from 6 to 16 to support the formatting of the current time directory 2.994286488639947



Reconstructing fields for mesh region0

Time = 2.994286488639947

Reconstructing FV fields

    Reconstructing volScalarFields

        nut
        pMean
        pPrime2Mean
        k

    Reconstructing volVectorFields

        U
Turbulent DFSEM patch cycl_half0: interpolating field R from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field L from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field U from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field R from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field L from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field U from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor0/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field R from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor1/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field L from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor1/../constant/boundaryData/cycl_half0/0"
Turbulent DFSEM patch cycl_half0: interpolating field U from "/home/hk354/SCRATCH/channel395DFSEM_hkmesh/processor1/../constant/boundaryData/cycl_half0/0"

Program received signal SIGFPE, Arithmetic exception.
0x00007ffff650001a in Foam::turbulentDFSEMInletFvPatchVectorField::turbulentDFSEMInletFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) () from /home/hk354/OpenFOAM/OpenFOAM-plus/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
(gdb) backtrace
#0  0x00007ffff650001a in Foam::turbulentDFSEMInletFvPatchVectorField::turbulentDFSEMInletFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) () from /home/hk354/OpenFOAM/OpenFOAM-plus/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#1  0x00007ffff6507b62 in Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::turbulentDFSEMInletFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) () from /home/hk354/OpenFOAM/OpenFOAM-plus/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#2  0x00000000004a991d in Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&)
    ()
#3  0x00000000004a9d64 in Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) ()
#4  0x00000000004aa59b in Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) ()
#5  0x00000000004eacbd in Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() ()
#6  0x00000000004536c4 in Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) [clone .constprop.1296] ()
#7  0x00000000004ec86b in Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvFieldReconstructor::reconstructFvVolumeField<Foam::Vector<double> >(Foam::IOobject const&) const ()
#8  0x00000000004ecc8c in void Foam::fvFieldReconstructor::reconstructFvVolumeFields<Foam::Vector<double> >(Foam::IOobjectList const&, Foam::HashSet<Foam::word, Foam::string::hash> const&)
    ()
#9  0x000000000043f7c0 in main ()
(gdb)