Turbulent intensity field function: Floating point exception [divide by zero]
 

Hello to all CFD mates,

I wish to define a Turbulent Intensity (TI) scalar in Star-CCM+ so that I can study the turbulence better. However, I receive a floating point exception [divide by zero] error when I created a TI field function. The function's definition is: sqrt((2/3)*$TurbulentKineticEnergy)/(mag($$Velocity))

From the definition, it is clear that the zero velocity at the wall is the main issue that need to be considered when defining the function. I was thinking to redefine the function where it returns zero value to system when it is divided by any zero value, or maybe return zero value when TI equals to infinity. But I have no idea whether this can be done or not.

I am still seeking for a solution to modify the solution. At the same time, I hope anyone that has experiences in this type of issue can provide some advises that help me to solve the problem.

Many thanks in advance.

TI field function in STAR-CCM+
 

Hello STAR-CCM+ user,

In the properties window you have to tick mark the "ignore boundary values" option to avoid Divide by zero error message.

Hope this helps and answers your question.