CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Post-Processing

error - Foam::sigFpe::sigHandler(int) at ??:?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 21, 2014, 12:43
Default error - Foam::sigFpe::sigHandler(int) at ??:?
  #1
Senior Member
 
Join Date: Nov 2012
Location: Bavaria
Posts: 144
Rep Power: 4
aylalisa is on a distinguished road
Dear Foamers,

I try to compute the turbulence intensity (utiliy) and receive an error:

I've tried:
Code:
Tu = sqrt(0.3 *mag(Rxx + Ryy + Rzz))/sqrt(sqr(U.component(0))+sqr(U.component(1))+sqr(U.component(2)));
and:
Code:
Tu = sqrt(0.3 *mag(Rxx + Ryy + Rzz))/sqrt(mag(sqr(U.component(0))+sqr(U.component(1))+sqr(U.component(2))));
error:
Code:
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
 #1  Foam::sigFpe::sigHandler(int) at ??:?
 #2   in "/lib64/libc.so.6"
 #3  
 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
 #4  
 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
 #5  
  at ??:?
 #6  __libc_start_main in "/lib64/libc.so.6"
 #7  
  at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
Could somebody help? What's wrong with the division?

Best regards,
Ayla
aylalisa is offline   Reply With Quote

Old   August 25, 2014, 22:56
Default
  #2
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 268
Rep Power: 9
kmooney is on a distinguished road
Hi Ayla,

Can you be 100% sure that U is non-zero throughout the entire domain?

If you want to make the evaluation more fpe-proof add a

Code:
+VSMALL
somewhere in the denominator.

I hope that helps!
Cheers,
Kyle
kmooney is offline   Reply With Quote

Old   August 26, 2014, 13:45
Default
  #3
Senior Member
 
Join Date: Nov 2012
Location: Bavaria
Posts: 144
Rep Power: 4
aylalisa is on a distinguished road
Hi Kyle,

thanks for your answer!
I think the code is wrong anyway.

I've changed it accordingly:
Quote:
Tu = (sqrt((1/3) * (Rxx + Ryy + Rzz))/stabilise(mag(U),SMALL*UUnits));
I guess entries of R are already divided by rho and units are m2/s2. How could I figure out?

The trace elements of R and the values of mag(U) are good but Tu = 0, why?
Whereas in ParaFoam I can compute Tu with the same formula and Tu is not 0.

Greetings!
Ayla
aylalisa is offline   Reply With Quote

Reply

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



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