CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

simpleFoam Floating point error (Core Dumped)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 24, 2013, 03:23
Post simpleFoam Floating point error (Core Dumped)
  #1
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 13
sam.ho is on a distinguished road
Hi
I am simulating a flow valve with following parameters
Density = 830 [kg/m3]
Viscosity = 8.26e-06 [kg/m.s]
Mass Flow rate = 145 lph (litres per hr)
Re=1.87202 X 10^5
The valve has one inlet and 4 outlets.
The boundary conditions are as follow

U file
[dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
INLET
{
type flowRateInletVelocity;
volumetricFlowRate constant 40.6e-6;
value uniform (0 0 0);
}

OUTLET_1
{
type inletOutlet;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
OUTLET_2
{
type inletOutlet;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
OUTLET_3
{
type inletOutlet;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
OUTLET_4
{
type inletOutlet;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
WALL_VALVE
{
type fixedValue;
value uniform (0 0 0);
}
WALL_SPRING
{
type fixedValue;
value uniform (0 0 0);
}
WALL_HOUSING
{
type fixedValue;
value uniform (0 0 0);
}
WALL_SPRING_SUPPORT
{
type fixedValue;
value uniform (0 0 0);
}
}

// ************************************************** *********************** //]

p file
[CODE]
[dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET_1
{
type fixedValue;
value $internalField;
}
OUTLET_2
{
type fixedValue;
value $internalField;
}
OUTLET_3
{
type fixedValue;
value $internalField;
}
OUTLET_4
{
type fixedValue;
value $internalField;
}

WALL_VALVE
{
type zeroGradient;
}
WALL_SPRING
{
type zeroGradient;
}
WALL_HOUSING
{
type zeroGradient;
}
WALL_SPRING_SUPPORT
{
type zeroGradient;
}


}

// ************************************************** *********************** //]

k file
[CODE]
[dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 0.0001133785;

boundaryField
{
INLET
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.0322;
value uniform 0.0001133785;
}

OUTLET_1
{
type inletOutlet;
inletValue uniform 0.0001133785;
value uniform 0.0001133785;
}
OUTLET_2
{
type inletOutlet;
inletValue uniform 0.0001133785;
value uniform 0.0001133785;
}
OUTLET_3
{
type inletOutlet;
inletValue uniform 0.0001133785;
value uniform 0.0001133785;
}
OUTLET_4
{
type inletOutlet;
inletValue uniform 0.0001133785;
value uniform 0.0001133785;
}

WALL_VALVE
{
type kqRWallFunction;
value uniform 0.0001133785;
}
WALL_SPRING
{
type kqRWallFunction;
value uniform 0.0001133785;
}
WALL_HOUSING
{
type kqRWallFunction;
value uniform 0.0001133785;
}
WALL_SPRING_SUPPORT
{
type kqRWallFunction;
value uniform 0.0001133785;
}

}


// ************************************************** *********************** //]
Epsilon file
[CODE]
[dimensions [ 0 2 -3 0 0 0 0 ];

internalField uniform 0.01533;

boundaryField
{

INLET
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.001; // This is based on diameter
value uniform 0.0153;
}

OUTLET_1
{
type inletOutlet;
inletValue uniform 0.0153;
value uniform 0.0153;
}
OUTLET_2
{
type inletOutlet;
inletValue uniform 0.0153;
value uniform 0.0153;
}
OUTLET_3
{
type inletOutlet;
inletValue uniform 0.0153;
value uniform 0.0153;
}
OUTLET_4
{
type inletOutlet;
inletValue uniform 0.0153;
value uniform 0.0153;
}
WALL_VALVE
{
type epsilonWallFunction;
value uniform 0.0153;
}

WALL_SPRING
{
type epsilonWallFunction;
value uniform 0.0153;
}

WALL_HOUSING
{
type epsilonWallFunction;
value uniform 0.0153;
}

WALL_SPRING_SUPPORT
{
type epsilonWallFunction;
value uniform 0.0153;
}

}


// ************************************************** *********************** //]

controlDict file
[CODE]
[application simpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 100;

deltaT 0.1;

writeControl timeStep;

writeInterval 1;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 8;

runTimeModifiable true;
]

fvSchemes file
[CODE]
[ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,R) bounded Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) bounded Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}
// ***************]

fvSolutions
[CODE][solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 0.01;
}

U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 3;
pRefCell 0;
pRefValue 0;
residualControl
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
}
}

relaxationFactors
{
fields
{
p 0.2;
}
equations
{
U 0.2;
k 0.2;
epsilon 0.2;
R 0.2;
nuTilda 0.2;
}
}


// ************************************************** *********************** //]

Snap shots of mesh and Velocity profile is attached. When i run simulation it ran up to 5.3 after that we got floating point error (core Dumped) message.Preceding the bounding epsilon problem.
Please Help me to get this simulation.
Attached Images
File Type: jpg Screenshot from 2013-12-24 11:30:00.jpg (73.3 KB, 18 views)
File Type: jpg Screenshot from 2013-12-24 12:09:41.jpg (57.4 KB, 17 views)
sam.ho is offline   Reply With Quote

Old   December 25, 2013, 20:21
Default
  #2
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
Hi Sam.ho,
This post is not the solution to your problem. But an advise that might help others to help you. Always test your setup with simple model test cases before proceeding towards the solution of your actual problem. It could be easier to catch the problems in a simple test problem, where you don't have so many boundary conditions and complex setup options.
Always post the error with files "(i.e. setup (+case) or test case)" with which you encounter the problem. If you can't share the full problem make a small test case demonstrating the problem. That might help others to find the problem and help you. Or may be you will solve the problem yourself in y process to make a test case.

Best Regards,
Ali
owayz is offline   Reply With Quote

Old   December 31, 2013, 06:57
Default
  #3
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 13
sam.ho is on a distinguished road
Hi owayz,

I have simulated simple cases like flow past circular cylinder and square cylinder and i got satisfactory results.
But if i apply the similar concept here I am getting error ?

Can anybody suggest a solution ?
sam.ho is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
What is causing this simpleFoam floating point exception? pajot OpenFOAM Running, Solving & CFD 15 April 20, 2018 03:39
[snappyHexMesh] How to define to right point for locationInMesh Mirage12 OpenFOAM Meshing & Mesh Conversion 7 March 13, 2016 14:07
Floating point exception (core dumped) Punt3r OpenFOAM Running, Solving & CFD 1 September 3, 2013 08:51
Floating point error and divergence detected aannjj FLUENT 0 July 2, 2013 03:44
Floating point error (core dumped) sihaqqi OpenFOAM Post-Processing 1 April 17, 2013 06:57


All times are GMT -4. The time now is 08:29.