CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   OFv1606+ - Error in running the channel395DFSEM tutorial (https://www.cfd-online.com/Forums/openfoam-solving/187780-ofv1606-error-running-channel395dfsem-tutorial.html)

rob3rt 0ng May 15, 2017 06:41

OFv1606+ - Error in running the channel395DFSEM tutorial
 
3 Attachment(s)
Hi,

I already have the version 2.4.x on my laptop, and I've successfully compiled the OFv1606+ (without the 64-bit integer bit version) as shown in the attached logMake file. I also attached my .bashrc file.

Upon executing Allrun command for channel395DFSEM tutorial, a dynamicCode folder pops up therein contains four sub-folders of codeStreamTemplate.C and a folder of platforms/linux64GccDPInt32Opt/lib/libcodeStream_xxxxx. Please see attached.

And when I run reconstructPar -latestTime, I have got this error:
Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  v1606+                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : v1606+Exec  : reconstructPar -latest
TimeDate  : May 15 2017
Time  : 19:03:27
Host  : "roberto-Precision-M6700"
PID    : 7903
Case  : /home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM
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 timeReconstructing fields for mesh region0
Time = 12.1
Reconstructing FV fields
    Reconstructing volScalarFields
        nut        k_0        pMean        p        pPrime2Mean        k        Q1        yPlus
    Reconstructing volVectorFields
        vorticity1        U_0
Turbulent DFSEM patch inlet: interpolating field R from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor0/../constant/boundaryData/inlet/0"
Turbulent DFSEM patch inlet: interpolating field L from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor0/../constant/boundaryData/inlet/0"
Turbulent DFSEM patch inlet: interpolating field U from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor0/../constant/boundaryData/inlet/0"
Turbulent DFSEM patch inlet: interpolating field R from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor1/../constant/boundaryData/inlet/0"
Turbulent DFSEM patch inlet: interpolating field L from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor1/../constant/boundaryData/inlet/0"
Turbulent DFSEM patch inlet: interpolating field U from "/home/robert/OpenFOAM/roberto/tutorial-v1606+/channel395DFSEM/processor1/../constant/boundaryData/inlet/0"
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::turbulentDFSEMInletFvPatchVectorField::turbulentDFSEMInletFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#4  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::turbulentDFSEMInletFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#5  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#6  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#7  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#8  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#9  ? at reconstructPar.C:?
#10  ? at ??:?
#11  ? at ??:?
#12  ? at ??:?
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? at ??:?Floating point exception

Has anyone encountered this problem before?

I'm guessing that it has something to do with the installation.

Thanks and regards,
Robert

rob3rt 0ng May 19, 2017 21:40

2 Attachment(s)
Please see attached is the self-generated dynamicCode folder.

The simulation runs okay. But when trying to reconstruct I get an error basically saying Floating Point Exception due to the scalar and vector operations done by the turbulentInletDFSEMFvField.

Can anyone help with this since maybe this is a pretty common issue for people who installed multiple FOAM versions.

Thanks and regards
Robert

[Moderator note: Moved from https://www.cfd-online.com/Forums/op...roblems-2.html and edited to remove the duplicate content]

wyldckat May 27, 2017 10:10

Quick answer: I've tested this now and this is already fixed in OpenFOAM+ v1612+. You will need to upgrade to v1612+ if you really want this feature.

In addition, you cannot use the case simulated with v1606+ and reconstruct it with v1612+, at least it didn't work for me.

rob3rt 0ng May 27, 2017 11:10

Thanks for the reply, Bruno.

I've just got the OF v1606+ version compiled on the HPC today (but havent been able to test anything yet). So are you saying that the turbulentInletDFSEM BC can't be used in v1606+ regardless the machine architecture?

Regards,
Robert

wyldckat May 27, 2017 12:00

Quick answer: From what I can figure out, the only problem is related to a limitation regarding reconstruction of the fields for the parallel case. In other words, you can still run the cases with v1606+, but you won't be able to reconstruct the results with reconstructPar. You can only do reconstruction with version v1612+.

The architecture issue that there was with 32 and 64-bit labels has nothing to do with this issue.


All times are GMT -4. The time now is 11:32.