[Rody-]

Hi,

I am simulate a turbine cascade with rhoCentralFoam. I got 2 meshes, which have the same geometry, topology and boundary. Mesh 2 is just twice as fine as Mesh 1.

The simulation runs fine with the first Mesh, I got a nice solution.
Now I want to use the solution from Mesh 1 for Mesh 2 and countinue simulation with a finer mesh.

I use mapFields -consistent and everything is fine, i got a time case for the Mesh 2 and no error. But after I run the simulation, i got this error:

#0 Foam::error:PrintStack(Foam::Ostream&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::heThermo<Foam:PsiThermo, Foam:PureMixture<Foam::constTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam:PerfectGas<Foam: :specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heThermo<Foam:PsiThermo, Foam:PureMixture<Foam::constTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam:PerfectGas<Foam: :specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5 Foam::hePsiThermo<Foam:PsiThermo, Foam:PureMixture<Foam::constTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam:PerfectGas<Foam: :specie> >, Foam::sensibleInternalEnergy> > > >::hePsiThermo(Foam::fvMesh const&, Foam::word const&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6 Foam:PsiThermo::addfvMeshConstructorToTable<Foam:: hePsiThermo<Foam:PsiThermo, Foam:PureMixture<Foam::constTransport<Foam::specie s::thermo<Foam::hConstThermo<Foam:PerfectGas<Foam: :specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7 Foam::autoPtr<Foam:PsiThermo> Foam::basicThermo::New<Foam:PsiThermo>(Foam::fvMes h const&, Foam::word const&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#8 Foam:PsiThermo::New(Foam::fvMesh const&, Foam::word const&) in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#9
in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/sw/openfoam/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 28988 on node knoten-09 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

If I start the simulation with Mesh2 without using the solution from Mesh1, there isn't a error like this.

/thermophysicalProperties:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

mixture
{
specie
{
nMoles 1;
molWeight 28.9;
}
thermodynamics
{
Cp 1005;
Hf 0;
}
transport
{
mu 1.8e-5;
Pr 0.7;
}
}


// ************************************************** *********************** //

[fportela]

Have you inspected (in ParaView) Mesh 2 after the mapping?

If the mesh is twice as small, did you take this into account when prescribing the time step?

[Rody-]

After mapping Mesh 2 I inspected it with Tecplot 360 and it looks fine...
I also reduced the deltaT and it didnt work.

[shakirbsm]

Hi every one I am also getting the same problem in forwardStep with finer mesh. I am getting the same error message in both cases if I run in parallel or in serial. Can some body help?

[dmck]

hey guys


I've also got this error when i refined my mesh, I've managed the run the same simulation on a coarser mesh and it ran, I'm using sonicFoam

my simulation runs for 1600 iterations, it needs to run to 2000 iterations,

i've checked where the simulation fails when i have a pressure wave being reflected at the end of a pipe.

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5

i'd be very interested if anyone has any ideas on how to solve this issue

thanks

[Jetfire]

Hi,

I am simulating a compressor using rhoPimpleDyMFoam , I too get the same error after a few iterations. Did anyone figure out how to solve this ? please help

[dmck]

Hi Jetfire

I was able to solve my issue by modifying my mesh, basically my issue was with the connection between my pipe and my exit volume, when I tried to refine my mesh, OpenFoam was unable to tie the meshes together. so by modifying my mesh I was able to solve my issue.

I would suggest looking at your boundary conditions and the equations you are using to solve, to ensure you have selected the correct type of equation for your particular problem

good luck

[Jetfire]

Originally Posted by dmck

Hi Jetfire

I was able to solve my issue by modifying my mesh, basically my issue was with the connection between my pipe and my exit volume, when I tried to refine my mesh, OpenFoam was unable to tie the meshes together. so by modifying my mesh I was able to solve my issue.

I would suggest looking at your boundary conditions and the equations you are using to solve, to ensure you have selected the correct type of equation for your particular problem

good luck

Thanks for your reply
Is this error not with the thermoPhysicalProperties settings ?

[1] #0  Foam::error::printStack(Foam::Ostream&) in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigHandler(int) in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
[1] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
[1] #5  
[1]  in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam"
[1] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7  
[1]  in "/home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam"
[EAT-Standalone:05495] *** Process received signal ***
[EAT-Standalone:05495] Signal: Floating point exception (8)
[EAT-Standalone:05495] Signal code:  (-6)
[EAT-Standalone:05495] Failing at address: 0x3e800001577
[EAT-Standalone:05495] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f9fe423b4a0]
[EAT-Standalone:05495] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f9fe423b425]
[EAT-Standalone:05495] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f9fe423b4a0]
[EAT-Standalone:05495] [ 3] /home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x2ab) [0x7f9fe97a1f0b]
[EAT-Standalone:05495] [ 4] /home/eatin/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE7correctEv+0x32) [0x7f9fe97af5f2]
[EAT-Standalone:05495] [ 5] rhoPimpleDyMFoam() [0x41f217]
[EAT-Standalone:05495] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f9fe422676d]
[EAT-Standalone:05495] [ 7] rhoPimpleDyMFoam() [0x42660d]
[EAT-Standalone:05495] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 5495 on node EAT-Standalone exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

[dmck]

yea that is where the error is but if you look at the initial points on this thread they are to do with meshing,

in my case my mesh had created like zero geometry conditions. which the solver was unable to solve! by correcting my mesh I was able to resolve the issue.

there is some good information on selecting the thermoPhysicalProperties online, if you check out the slides on the link below - it could be that you haven't selected a incompatible set of parameters for your thermophysical model, the compatable links are shown on slide 6-8 on the link, or maybe you have specified a boundary condition wrongly or not included a boundary condition

as I haven't seen your files or your error message I'm purely speculating and these are just some of the things that I've done wrong in the past!!!

http://www.tfd.chalmers.se/~hani/kur...d_HN141009.pdf