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

Division by zero in Xoodles

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 18, 2007, 07:16
Default Description: While starting t
  #1
Member
 
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 86
Rep Power: 8
hannes is on a distinguished road
Description:
While starting the Xoodles/pitzDaily3D tutorial, Xoodles fails and reports a floating point error.
The reason seems to be in XiFoam/bEqn.H in line 79. There, muu is divided by (rhou*epsilon). It seems, that epsilon may be zero in some cases.
Change the statement to:

volScalarField tauEta = sqrt(thermo->muu()/(
rhou*epsilon + dimensionedScalar("1e-6", rhou.dimensions()*epsilon.dimensions(), 1e-6)
));


Solver/Application:
Xoodles, XiFoam, engineFoam

Source file:
$FOAM_APP/solvers/combustion/XiFoam/bEqn.H

Testcase:
$FOAM_TUTORIALS/Xoodles/pitzDaily3D

Platform:
x86_64

Version:
1.4 (also 1.3)

Notes:
hannes is offline   Reply With Quote

Old   April 18, 2007, 08:13
Default Thanks for the bug report, the
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
Thanks for the bug report, the problem is that in the initial k field the internal and wall boundary values are set to zero from which a zero value of epsilon is calculated. The best solution to this problem is to initialise the k field to a small but non-zero value and also select a non-zero value for the walls.

Henry
henry is offline   Reply With Quote

Old   August 4, 2008, 09:34
Default Hi! It seems that this prob
  #3
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,834
Rep Power: 38
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi!

It seems that this problem is fixed in 1.5 for the Xoodles/pitzDaily (default value set to 1e-8) but NOT in the Xoodles/pitzDaily3D-tutorial-case (values are still 0)

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   August 4, 2008, 11:04
Default Thanks for the report, I will
  #4
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
Thanks for the report, I will fix it in pitzDaily3D also.

H
henry 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
About flow division sunnysun OpenFOAM Running, Solving & CFD 5 March 2, 2009 12:34
Division by zero vitke OpenFOAM Running, Solving & CFD 5 September 1, 2008 05:35
PitzDaily tutorial with Xoodles christianvhoersten OpenFOAM Running, Solving & CFD 0 January 11, 2008 10:17
How to run the cases Xoodles%5cpitzDaily3D and Xoodles%5cpitzDaily cfdfans OpenFOAM Running, Solving & CFD 4 October 18, 2007 03:47
Soot model Xoodles deepblue17 OpenFOAM Running, Solving & CFD 0 May 2, 2006 11:57


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