 |
 |
 |
|
|
|
[Sponsors] | ||||
March 7, 2014, 20:31
|
floating point exception
|
#1 |
|
New Member
J. Brindley
Join Date: Mar 2014
Posts: 11
Rep Power: 12  |
Not sure what's going on. Can't get simpleFoam to run
Code:
Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.081410493, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.061596729, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.091684187, No Iterations 3 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0059667553, No Iterations 5 GAMG: Solving for p, Initial residual = 0.40650934, Final residual = 0.0012374181, No Iterations 3 GAMG: Solving for p, Initial residual = 0.036894522, Final residual = 0.00027957842, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0062564286, Final residual = 4.8770264e-05, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0013606213, Final residual = 1.2695641e-05, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00035198339, Final residual = 3.117196e-06, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00010876533, Final residual = 1.0124916e-06, No Iterations 3 GAMG: Solving for p, Initial residual = 3.7107364e-05, Final residual = 3.2700493e-07, No Iterations 3 GAMG: Solving for p, Initial residual = 1.4011258e-05, Final residual = 1.2603484e-07, No Iterations 3 GAMG: Solving for p, Initial residual = 5.5796704e-06, Final residual = 3.7850784e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 2.3340101e-06, Final residual = 6.442118e+119, No Iterations 1000 time step continuity errors : sum local = 1.2632859e+123, global = 4.4873194e+101, cumulative = 4.4873194e+101 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:? #9 Foam::incompressible::RASModels::kOmegaSST::correct() at ??:? #10 at ??:? #11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #12 at ??:? Floating point exception (core dumped) |
|
|
|
|
|
March 7, 2014, 21:06
|
|
#2 |
|
New Member
J. Brindley
Join Date: Mar 2014
Posts: 11
Rep Power: 12 |
[
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : checkMesh
Date : Mar 07 2014
Time : 20:04:09
Host : "AVA-411296"
PID : 26818
Case : /home/kathyb/OpenFOAM/kathyb-2.3.0/run/SplitTipTesting/AOA/simplesolvefixed(donotuse)
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points: 30927
faces: 72559
internal faces: 60772
cells: 21783
faces per cell: 6.1208741
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 15749
prisms: 1429
wedges: 0
pyramids: 0
tet wedges: 128
tetrahedra: 2
polyhedra: 4475
Breakdown of polyhedra by number of faces:
faces number of cells
4 1186
5 814
6 740
7 218
8 90
9 678
10 86
11 47
12 401
13 3
14 6
15 150
18 56
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
inlet 0 0 ok (empty)
outlet 0 0 ok (empty)
walls 40 67 ok (non-closed singly connected)
frontBack 0 0 ok (empty)
Foil 11747 13554 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-0.020220271 -0.010812593 -0.00023150215) (0.020260984 0.010818352 0.15183358)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-5.8212676e-17 4.5742309e-17 1.775036e-17) OK.
Max cell openness = 2.9738872e-16 OK.
Max aspect ratio = 21.38308 OK.
Minimum face area = 6.6009305e-09. Maximum face area = 1.3716404e-05. Face area magnitudes OK.
Min volume = 1.1250083e-12. Max volume = 2.2923998e-08. Total volume = 2.8560307e-05. Cell volumes OK.
Mesh non-orthogonality Max: 62.160689 average: 13.151348
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.6576963 OK.
Coupled point location match (average 0) OK.
Mesh OK.
End
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 0.02;
relTol 0;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0.1;
}
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0.1;
}
omega
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0.1;
}
cellDisplacement
{
solver GAMG;
tolerance 1e-5;
relTol 0;
smoother GaussSeridel;
cacheAgglomeration true;
nCellsInCoursestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
pReffCell 0;
pRefValue 0;
residualControl
{
p 1e-4;
U 1e-4;
}
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
omega 0.7;
}
}
cache
{
grad(U);
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
grad(k) cellLimited Gauss linear 1;
grad(omega) cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss upwind;;
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited corrected 0.5;
laplacian(DomegaEff,omega) Gauss linear limited 0.33;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.33;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //
|
|
|
|
|
|
|
March 8, 2014, 17:56
|
|
#3 |
|
Senior Member
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 299
Rep Power: 23 |
Hi olddirtybrindly,
Do you get the same crash when you run the case with pisofoam? Cheers, L |
|
|
|
|
|
 |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Floating point exception with pimpleDyMFoam | ebah6 | OpenFOAM Running, Solving & CFD | 9 | November 1, 2017 06:58 |
| Inlet Velocity Profile BC - Floating Point exception during solution initialization | Janshi | STAR-CCM+ | 4 | March 14, 2012 11:21 |
| simpleFoam Floating point exception error -help | sudhasran | OpenFOAM Running, Solving & CFD | 3 | March 12, 2012 17:23 |
| Pipe flow in settlingFoam floating point exception | jochemvandenbosch | OpenFOAM Running, Solving & CFD | 4 | February 16, 2012 04:24 |
| block-structured mesh for t-junction | Robert@cfd | ANSYS Meshing & Geometry | 20 | November 11, 2011 05:59 |
Linear Mode
