[aylalisa]

Dear Foamers,

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

I've tried:

Tu = sqrt(0.3 *mag(Rxx + Ryy + Rzz))/sqrt(sqr(U.component(0))+sqr(U.component(1))+sqr(U.component(2)));

and:

Tu = sqrt(0.3 *mag(Rxx + Ryy + Rzz))/sqrt(mag(sqr(U.component(0))+sqr(U.component(1))+sqr(U.component(2))));

error:

#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