[n00bie]

Hi everybody,

I am currently trying to run a simulation with rhoSimpleFoam.
I used the tutorial case "angledDuctExplicitFixedCoeff",
replaced the files in the polyMesh-folder with my own mesh and adjusted the boundary conditions.
When I tried to run rhoSimpleFoam on the case, I got the following error message:

#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::heThermo<Foam:: psiThermo, Foam:: pureMixture<Foam::sutherlandTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam: psiThermo::addfvMeshConstructorToTable<Foam::hePsi Thermo<Foam: psiThermo, Foam: pureMixture<Foam::sutherlandTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam: psiThermo> Foam::basicThermo::New<Foam:: psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam:: psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7
at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
at ??:?
Floating point exception (core dumped)

I have already searched the forums and I know, that to display the lines I would need to use the debug-version of OpenFOAM.
However, since I am on a university computer, I would rather not change anything in the program files of OpenFOAM.
I am using OpenFOAM 2.3.1 by the way.

From the error message I gather, that it is a problem with the thermophysical properties.
I did not modify them in any way, they are still the same as in the tutorial case.
I am not sure how to change them since I can't figure out, where they are used.

If you need to have a look at any of my files let me know and I will post them as soon as possible,
but so far I don't really know what you are going to need.

By the way I don't think it is a problem with k or epsilon being zero.

I hope you can help me.

Sincerely,
Sabrina

[n00bie]

Actually I think it may be a problem with my mesh.
I executed a checkMesh -allTopology -allGeometry and it failed 2 checks.
So I assume, the error was caused by this and I will need to
modify my mesh in order to run rhoSimpleFoam properly.

[enginpower]

Hi Sabrina,
I'm stuck to a similar problem, but my mesh is all ok.

Could you provide some useful instructions, you found about the debug, to understand which line of the code gets the error?
I would be really grateful if you or someone else want/s to share these instructions here or in my post

Thanks

[n00bie]

I think most of the advice to help with this error I found in this thread:
Foam::error::PrintStack

[n00bie]

It took me a while to find it again, but here are the instructions on the openfoamwiki.

https://openfoamwiki.net/index.php/FAQ/Troubleshooting

[ingarcia1703]

THe problem appears due to the boundary conditions. It may be due to the Pressure outlet, the units are established in Pa and you should take into account a total pressure, not a p/rho (m2/s2) pressure as it is taken in incompressible flow (simpleFoam).