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/)
-   -   Divition by zero error in rhoSimpleFoam with SA (https://www.cfd-online.com/Forums/openfoam-solving/184970-divition-zero-error-rhosimplefoam-sa.html)

andreamc March 15, 2017 16:25

Divition by zero error in rhoSimpleFoam with SA
 
Hello

I am trying to simulate the bump-in-channel verification case from the TRM web by NASA with compressible, steady solver rhoSimpleFoam and SpalartAllmaras with OF 3.1. I've succeed at simulating the case for the three coarse grids. However when simulating exactly the same case but with a finer mesh it crashes. The error I get is:

#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::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? at ??:?
Excepción de coma flotante (`core' generado)

I know this error comes out when there is a mathematical inconsistency like division by zero. So I've checked all the boundary conditions but nothing I've changed has worked so far and they are the same BC I used in the previous simulation with coarse meshes. I've changed also schemes but this has not show any improvement in the simulation.
Last I think from the error I get that the problem is in the thermophysicalProperties file but I've also checked that and nothing :(

My thermophysicalProperties file looks like follows:

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

mixture
{
specie
{
nMoles 1;
molWeight 28.9;
}
thermodynamics
{
Cp 1007;//1004.5;
Hf 0;// 2.544e+06;
}
transport
{
As 1.4792e-06;
Ts 116;
// mu 1.8e-05;
// Pr 0.7;
}
}


In conclusion, I've tried to change everything from BC, schemes, solver, and thermoProperties and nothing has worked so far.
Could anyone help me to understand this?

Ps. I tried with incompressible simpleFoam and it runs, I haven't check the results but it doesn't crashes.
Thanks in advance

Cajal June 1, 2017 17:36

Hello Andrea,

I'm having the same issue. Any update on this?

Thanks.

andreamc June 2, 2017 14:01

Hi Arturo

I was able to solve my problem by changing the boundary condition alphat in the wall to alphatJayatilllekeWallFunction. Additionally I had to increase the number of outercorrectors.
That solved my problem, hope it works for you!

Cajal June 2, 2017 14:16

Andrea,

Thanks for your kindly reply, I will try to change the parameters you mention.

Best.


All times are GMT -4. The time now is 03:25.