# twoPhaseEulerFoam - floating point exception (nutb)

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

 April 9, 2009, 05:07 twoPhaseEulerFoam - floating point exception (nutb) #1 New Member   Join Date: Mar 2009 Posts: 28 Rep Power: 9 Hi, I just started with the twoPhaseEulerFoam Solver. I got while "calculating nutb" a floating point exception. I saw that nutb is turbulent kinematic viscosity of phase b. Where could be the problem? Br

April 9, 2009, 10:41
#2
Member

Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 83
Rep Power: 9
Quote:
 Originally Posted by Hectux Hi, I just started with the twoPhaseEulerFoam Solver. I got while "calculating nutb" a floating point exception. I saw that nutb is turbulent kinematic viscosity of phase b. Where could be the problem? Br
I guess you're using 1.5? If you look at the createFields.H in the solvers's folder and search for the line "Calculating field nutb" you will find:

Info<< "Calculating field nutb\n" << endl;
volScalarField nutb
(
IOobject
(
"nutb",
runTime.timeName(),
mesh,
IOobject::AUTO_WRITE
),
Cmu*sqr(k)/epsilon
);

Most floating point exceptions are divide-by-zeros. As you can see the square of the turbulent kinetic energy, k, is divided by epsilon. I guess you have a zero in your epsilon field? You can eliminate the floating point exception by using a small value such as 1e-13 instead.

However, in 1.5.x version this problem is eliminated by using:

Cmu*sqr(k)/max(epsilon, dimensionedScalar("smallEps",epsilon.dimensions(), 1e-6))

Where the value of epsilon in the denominator is limited to larger than 1e-6 and thus removing the possibility of a floating point exception caused by a zero in the epsilon field.

So another option is to update to 1.5.x or to modify the line yourself and recompile the solver.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Ricky Wong FLUENT 11 September 8, 2016 02:46 maryam CFX 3 May 18, 2010 07:44 Alex FLUENT 2 April 21, 2009 01:29 Riyaz Main CFD Forum 0 November 14, 2008 07:30 touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27

All times are GMT -4. The time now is 20:34.