CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Divition by zero error in rhoSimpleFoam with SA

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 15, 2017, 17:25
Default Divition by zero error in rhoSimpleFoam with SA
  #1
New Member
 
Andrea Matiz C
Join Date: Feb 2017
Posts: 9
Rep Power: 2
andreamc is on a distinguished road
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:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#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() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa 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
andreamc is offline   Reply With Quote

Old   June 1, 2017, 17:36
Default
  #2
New Member
 
Arturo Cajal
Join Date: May 2017
Posts: 2
Rep Power: 0
Cajal is on a distinguished road
Hello Andrea,

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

Thanks.
Cajal is offline   Reply With Quote

Old   June 2, 2017, 14:01
Default
  #3
New Member
 
Andrea Matiz C
Join Date: Feb 2017
Posts: 9
Rep Power: 2
andreamc is on a distinguished road
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!
andreamc is offline   Reply With Quote

Old   June 2, 2017, 14:16
Default
  #4
New Member
 
Arturo Cajal
Join Date: May 2017
Posts: 2
Rep Power: 0
Cajal is on a distinguished road
Andrea,

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

Best.
Cajal is offline   Reply With Quote

Reply

Tags
bump, rhosimplefoam, spalartallmaras

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure instability with rhoSimpleFoam philipp. OpenFOAM 13 October 30, 2016 04:39
rhoSimpleFoam angledDuctExplicitFixedCoeff tutorial fails in parallel donQi OpenFOAM Running, Solving & CFD 1 February 22, 2016 20:47
Switching from simpleFoam to rhoSimpleFoam sebastian OpenFOAM 11 January 7, 2015 05:32
rhoSimpleFoam. patchField error. 123 OpenFOAM Running, Solving & CFD 4 June 6, 2014 15:22
transsonic nozzle with rhoSimpleFoam Unseen OpenFOAM Running, Solving & CFD 7 April 16, 2014 03:38


All times are GMT -4. The time now is 15:02.