[nick0520]

Hi everyone!
I'am a new user of the forum. Firstly, I would like to show you my appreciation for this website since it has helped me a lot sincee the beginning of my journey in openFoam/CFD.

I am student currently working on a thesis about flow simulations around an f1-2022-like car.
I am running the simulations on the university cluster and I have completed my first simulation lately. I am trying to visualize the results in paraview locally, on my computer, so I ran foamToVTK on the cluster but it does not complete all the time steps.
The simulations has 15 time steps but foamToVTK stops at the fourth one.

I have no clue why is doing this, so if anyone is willing to help me I would highly appreciate it.

My guess is that is something related to the dimension of the files (" void vtkWriteOps double" - double precision?).
In controlDict I am using the following:

writeFormat ascii;

writePrecision 12;

writeCompression compressed;

What do you think?

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : dev-69858a80ec41
Exec : foamToVTK
Date : Jun 11 2022
Time : 04:30:12
Host : "cm004.hpc.nyu.edu"
PID : 1092320
I/O : uncollated
Case : /scratch/nm3584/openfoam-test/formula4/formula3
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Time: 0
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_0.vtk"
Original cells:8841237 points:9911068 Additional cells:9288785 additional points:795265

Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/frontAndBack/frontAndBack_0.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/inlet/inlet_0.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/outlet/outlet_0.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/lowerWall/lowerWall_0.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/upperWall/upperWall_0.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/F1-75/F1-75_0.vtk"
Time: 100
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/frontAndBack/frontAndBack_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/inlet/inlet_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/outlet/outlet_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/lowerWall/lowerWall_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/upperWall/upperWall_100.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/F1-75/F1-75_100.vtk"
surfScalarFields : phi
Time: 200
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/frontAndBack/frontAndBack_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/inlet/inlet_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/outlet/outlet_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/lowerWall/lowerWall_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/upperWall/upperWall_200.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/F1-75/F1-75_200.vtk"
surfScalarFields : phi
Time: 300
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/frontAndBack/frontAndBack_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/inlet/inlet_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/outlet/outlet_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/lowerWall/lowerWall_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/upperWall/upperWall_300.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/F1-75/F1-75_300.vtk"
surfScalarFields : phi
Time: 400
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/frontAndBack/frontAndBack_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/inlet/inlet_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/outlet/outlet_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/lowerWall/lowerWall_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/upperWall/upperWall_400.vtk"
Patch : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/F1-75/F1-75_400.vtk"
surfScalarFields : phi
Time: 500
volScalarFields : p nut k omega
volVectorFields : U

Internal : "/scratch/nm3584/openfoam-test/formula4/formula3/VTK/formula3_500.vtk"
#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::vtkWriteOps::insert(double, Foam:: DynamicList<float, 0u, 2u, 1u>&) at ??:?
#4 void Foam::vtkWriteOps::write<double>(std:: ostream&, bool, Foam:: DimensionedField<double, Foam::volMesh> const&, Foam::vtkMesh const&) in "/opt/openfoam-dev/platforms/linux64GccDPInt32Opt/bin/foamToVTK"
#5 ? in "/opt/openfoam-dev/platforms/linux64GccDPInt32Opt/bin/foamToVTK"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/opt/openfoam-dev/platforms/linux64GccDPInt32Opt/bin/foamToVTK"
/scratch/work/public/apps/openfoam/20220221/run-openfoam-dev-paraview56_20220221.bash: line 27: 1092301 Floating point exceptionsingularity exec $nv --overlay /scratch/work/public/singularity/openfoam-dev-paraview56_20220221-depends.sqf:ro /scratch/work/public/singularity/openfoam-dev-paraview56_20220221.sif /bin/bash -c "
source /ext3/openmpi-env.sh
source /opt/openfoam-dev/etc/bashrc
export PATH=.:\$PATH
export LD_LIBRARY_PATH=/opt/slurm/lib64:\$LD_LIBRARY_PATH
$args
"

[olesen]

Floating point exception sounds like there might be some interpolation happening that provokes a divide-by-zero error (just a wild guess based on the error message).


If you have access to any recent www.openfoam.com version (current latest is v2112), it would be interesting to see if the problem also exists there. In this version, the VTK conversion is handled using a completely different infrastructure (also supports vtu,vtp,vtm output).


Another thing to try might be foamToEnsight. The ensight case files can also be opened in paraview.

[nick0520]

Originally Posted by olesen

Floating point exception sounds like there might be some interpolation happening that provokes a divide-by-zero error (just a wild guess based on the error message).

Thanks for your quick reply!

Originally Posted by olesen

If you have access to any recent www.openfoam.com version (current latest is v2112), it would be interesting to see if the problem also exists there. In this version, the VTK conversion is handled using a completely different infrastructure (also supports vtu,vtp,vtm output).

I have access to what is installed in the cluster, from my understanding is a recent version; If I'm not wrong it was the latest release in late April 22'.

Originally Posted by olesen

Another thing to try might be foamToEnsight. The ensight case files can also be opened in paraview.

I'll give it a try and let you know, thanks!

[olesen]

Originally Posted by nick0520

I have access to what is installed in the cluster, from my understanding is a recent version; If I'm not wrong it was the latest release in late April 22'.

I was referring to a release of the www.openfoam.com version (OpenCFD Ltd.) whereas your headers indicate you have the www.openfoam.org version (OpenFOAM Foundation Ltd.) - not quite the same thing. You can also visit my blog post (https://olesenm.github.io/2020/11/11...de-provenance/) for some background.

[nick0520]

Originally Posted by olesen

I was referring to a release of the www.openfoam.com version (OpenCFD Ltd.) whereas your headers indicate you have the www.openfoam.org version (OpenFOAM Foundation Ltd.) - not quite the same thing. You can also visit my blog post (https://olesenm.github.io/2020/11/11...de-provenance/) for some background.

I apologize for my late reply.

I was able to solve my issue, most likely caused by some "bugs" (holes) in my geometry.

[rideste]

Originally Posted by nick0520

I apologize for my late reply.

I was able to solve my issue, most likely caused by some "bugs" (holes) in my geometry.

Hi nicola, i have a geometry with a hole and foamToVTK give me the same error.
How did you solve this bug?