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

Hello Andrea,

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

Thanks.

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!

Andrea,

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

Best.