[letzel]

Dear Foamers,

in a modified LES solver, I calculate the turbulence intensity Iav from Uav and U2av as follows:

Uav  = ( Uav*(time-startAveraging) + U*dt ) / (time-startAveraging + dt);
for ( int i = 1; i <= 3; ++i)
{
  U2av.component(i) = ( U2av.component(i)*(time-startAveraging) + sqr(U.component(i))*dt ) / (time-startAveraging + dt);
}
for ( int i = 1; i <= 3; ++i)
{
  Iav.component(i) = sqrt( U2av.component(i) - sqr( Uav.component(i) ) ) / ( mag(Uav)+magUTinyOffset );  // avoid zero division
}

The term U2av.component(i) - sqr( Uav.component(i) ) may be negative, which then causes a fatal floating point error on calculation of its square root. I would like to set the term to zero in such cases, but don't know how to achieve this.

In Fortran, the WHERE function lets you conditionally manipulate arrays:

WHERE(array<0.0) array=0.0

Does OpenFOAM offer a similar functionality, or C++? If so, could you please explain how to use it?

Best regards,
Marcus