[mbcx4jc2]

Hi everyone,

I'm running into an error running simpleFoam (this problem also happens with icoFoam).

Casefile can be found at the following address:
https://drive.google.com/file/d/0B5R...it?usp=sharing

I'm getting the error pasted below. The log file also reveals that the solution is diverging which I suspect why this is happening. checkMesh looks good, boundary conditions appear sensible. Any ideas would be much appreciated!

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4
at simpleFoam.C:0
#5
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

[alexeym]

Hi,

Obviously you've got diverging solution. Look at all those 1001 in log file, it is default maximum number of iteration of the matrix solver after which it will give up even it doesn't get the solution within given tolerances.

I'm not quite get from the case files what you are trying to do but I've attached archive with the settings that can bear your setup. But the velocities in the narrow part of the mesh are around 4 km/s so I guess it doesn't make much sense.

I've attached only archive of the system folder as I've made no changes in other folders.

[mbcx4jc2]

Thanks alexeym!

It's running for me up to timestep 19.8 so obviously this is a big improvement. Basically, I'm just trying to measure permeability across this mesh by sampling pressure difference and flux through the two planes that are my inlets and outlets. Perhaps reducing the velocity bc at the inlet will reduce the absolute velocity in the narrow section but I imagine relatively speaking, the velocity would still be this high. Would reducing the initial velocity amount help the solution converge?

[alexeym]

Well,

You can try to estimate velocity in the narrow part from cross-section areas ratio.

As you are trying to run simulation in laminar case, you should also estimate Reynolds number to see if it is really so.

And answering your question: yes, it will help.

[jiaojiao]

Originally Posted by mbcx4jc2

Hi everyone,

I'm running into an error running simpleFoam (this problem also happens with icoFoam).

Casefile can be found at the following address:
https://drive.google.com/file/d/0B5R...it?usp=sharing

I'm getting the error pasted below. The log file also reveals that the solution is diverging which I suspect why this is happening. checkMesh looks good, boundary conditions appear sensible. Any ideas would be much appreciated!

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4
at simpleFoam.C:0
#5
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

I have a similar problem,I am running the MRFSimplefoam with of2.1 to simulate the 3D wind turbine, but it cause an error as I pasted

#0  Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#2  __restore_rt at sigaction.c:0
#3  Foam::writeFuns::insert(double, Foam::DynamicList<float, 0u, 2u, 1u>&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
#4  void Foam::writeFuns::write<double>(std::ostream&, bool, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh> const&, Foam::vtkMesh const&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
#5  void Foam::internalWriter::write<double>(Foam::volPointInterpolation const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
#6  main in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
#7  __libc_start_main in "/lib64/libc.so.6"
#8  Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
浮点数例外

can anybody help me? Where the problems are and how can I fix it? my fvshemes file is

ddtSchemes
{
    default         steadyState;
}

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

divSchemes
{
    default         Gauss upwind;
    div(phi,U)       Gauss limitedLinearV 1;//Gauss upwind;
    div(phi,k)       Gauss upwind;
    div(phi,epsilon)  Gauss upwind;//linear;
    div(phi,nuTilda) Gauss linearUpwind grad(nuTilda);
    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(DnuTildaEff,nuTilda) Gauss linear corrected;
    laplacian(1,p)  Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
}

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

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}

the log is

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}

No field sources present


SIMPLE: convergence criteria
    field p     tolerance 1e-05
    field U     tolerance 1e-05
    field nuTilda     tolerance 1e-05


Starting time loop

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0116175, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0478811, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0116854, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0747252, No Iterations 1
time step continuity errors : sum local = 1222.25, global = 625.783, cumulative = 625.783
smoothSolver:  Solving for epsilon, Initial residual = 0.999886, Final residual = 0.00670035, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.00694971, No Iterations 2
ExecutionTime = 39.72 s  ClockTime = 40 s

Time = 2

smoothSolver:  Solving for Ux, Initial residual = 0.255939, Final residual = 0.00653363, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.207718, Final residual = 0.00447832, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.253055, Final residual = 0.00686111, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.16482, No Iterations 1
time step continuity errors : sum local = 1507.88, global = 682.725, cumulative = 1308.51
smoothSolver:  Solving for epsilon, Initial residual = 0.36138, Final residual = 0.00191422, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.577132, Final residual = 0.00949441, No Iterations 2
ExecutionTime = 60.56 s  ClockTime = 66 s

Time = 3

smoothSolver:  Solving for Ux, Initial residual = 0.199729, Final residual = 0.00633206, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.27506, Final residual = 0.00436556, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.200162, Final residual = 0.00666538, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.486066, Final residual = 0.0946856, No Iterations 1
time step continuity errors : sum local = 824.775, global = 340.359, cumulative = 1648.87
smoothSolver:  Solving for epsilon, Initial residual = 0.298414, Final residual = 0.00162272, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.23253, Final residual = 0.00363581, No Iterations 2
ExecutionTime = 80.71 s  ClockTime = 86 s

Time = 4

smoothSolver:  Solving for Ux, Initial residual = 0.206463, Final residual = 0.00682647, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.447125, Final residual = 0.00887106, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.206458, Final residual = 0.00696421, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.443974, Final residual = 0.0504603, No Iterations 1
time step continuity errors : sum local = 661.598, global = -177.288, cumulative = 1471.58
smoothSolver:  Solving for epsilon, Initial residual = 0.432129, Final residual = 0.00656384, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.43187, Final residual = 0.00188758, No Iterations 2
ExecutionTime = 101.48 s  ClockTime = 112 s

Time = 5

smoothSolver:  Solving for Ux, Initial residual = 0.180274, Final residual = 0.00628132, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.279662, Final residual = 0.00753944, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.181797, Final residual = 0.00637471, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.615593, Final residual = 0.0888842, No Iterations 1
time step continuity errors : sum local = 1278.21, global = -675.965, cumulative = 795.615
smoothSolver:  Solving for epsilon, Initial residual = 0.310979, Final residual = 0.00371577, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.273199, Final residual = 0.00232058, No Iterations 2
ExecutionTime = 123.77 s  ClockTime = 134 s

Time = 6

smoothSolver:  Solving for Ux, Initial residual = 0.134419, Final residual = 0.00424664, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.13053, Final residual = 0.00380063, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.135395, Final residual = 0.00435766, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.784403, Final residual = 0.141384, No Iterations 1
time step continuity errors : sum local = 1529.09, global = -741.707, cumulative = 53.9077
smoothSolver:  Solving for epsilon, Initial residual = 0.231691, Final residual = 0.00282583, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.166224, Final residual = 0.00223461, No Iterations 2
ExecutionTime = 145.56 s  ClockTime = 161 s

Time = 7

smoothSolver:  Solving for Ux, Initial residual = 0.136914, Final residual = 0.00399246, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.105305, Final residual = 0.00289856, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.13624, Final residual = 0.00405558, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.721444, Final residual = 0.143873, No Iterations 1
time step continuity errors : sum local = 1409.36, global = -557.136, cumulative = -503.228
smoothSolver:  Solving for epsilon, Initial residual = 0.209118, Final residual = 0.0031759, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.156912, Final residual = 0.00272508, No Iterations 2
ExecutionTime = 167.51 s  ClockTime = 183 s

Time = 8

smoothSolver:  Solving for Ux, Initial residual = 0.142662, Final residual = 0.00416516, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.106274, Final residual = 0.0027194, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.142317, Final residual = 0.00423871, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.50711, Final residual = 0.0923087, No Iterations 1
time step continuity errors : sum local = 1146.77, global = -448.292, cumulative = -951.52
smoothSolver:  Solving for epsilon, Initial residual = 0.225749, Final residual = 0.00427787, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.174881, Final residual = 0.00347794, No Iterations 2
ExecutionTime = 189.73 s  ClockTime = 210 s

Time = 9

smoothSolver:  Solving for Ux, Initial residual = 0.132193, Final residual = 0.00415394, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0937998, Final residual = 0.00250987, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.133503, Final residual = 0.00419109, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.41403, Final residual = 0.0826546, No Iterations 1
time step continuity errors : sum local = 1231.54, global = -483.631, cumulative = -1435.15
smoothSolver:  Solving for epsilon, Initial residual = 0.238165, Final residual = 0.00537919, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.18193, Final residual = 0.00424883, No Iterations 2
ExecutionTime = 211.34 s  ClockTime = 232 s

Time = 10

smoothSolver:  Solving for Ux, Initial residual = 0.114206, Final residual = 0.00371409, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0774936, Final residual = 0.00227964, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.116838, Final residual = 0.00371665, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.390742, Final residual = 0.0831381, No Iterations 1
time step continuity errors : sum local = 1287.78, global = -275.8, cumulative = -1710.95
smoothSolver:  Solving for epsilon, Initial residual = 0.203471, Final residual = 0.00501127, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.169364, Final residual = 0.00457364, No Iterations 2
ExecutionTime = 233.55 s  ClockTime = 259 s

Time = 11

smoothSolver:  Solving for Ux, Initial residual = 0.103788, Final residual = 0.00356374, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0637353, Final residual = 0.00211766, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.104148, Final residual = 0.0034536, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.380292, Final residual = 0.0869304, No Iterations 1
time step continuity errors : sum local = 1302.78, global = -205.145, cumulative = -1916.1
smoothSolver:  Solving for epsilon, Initial residual = 0.160339, Final residual = 0.00424108, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.14831, Final residual = 0.00447062, No Iterations 2
ExecutionTime = 255.03 s  ClockTime = 281 s

Time = 12

smoothSolver:  Solving for Ux, Initial residual = 0.0950877, Final residual = 0.0033789, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0547081, Final residual = 0.0018934, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0942434, Final residual = 0.00320187, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.36055, Final residual = 0.0855903, No Iterations 1
time step continuity errors : sum local = 1272.49, global = -177.985, cumulative = -2094.08
smoothSolver:  Solving for epsilon, Initial residual = 0.141591, Final residual = 0.0039718, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.139064, Final residual = 0.00437892, No Iterations 2
ExecutionTime = 276.63 s  ClockTime = 307 s

Time = 13

smoothSolver:  Solving for Ux, Initial residual = 0.08845, Final residual = 0.00326364, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0481885, Final residual = 0.00174469, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0879515, Final residual = 0.00315651, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.346286, Final residual = 0.081073, No Iterations 1
time step continuity errors : sum local = 1200.51, global = -165.516, cumulative = -2259.6
smoothSolver:  Solving for epsilon, Initial residual = 0.13731, Final residual = 0.00392148, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.13561, Final residual = 0.0043371, No Iterations 2
ExecutionTime = 296.61 s  ClockTime = 327 s

Time = 14

smoothSolver:  Solving for Ux, Initial residual = 0.0826409, Final residual = 0.00302597, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0436394, Final residual = 0.0016168, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0823954, Final residual = 0.00302672, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.334584, Final residual = 0.0735457, No Iterations 1
time step continuity errors : sum local = 1086.94, global = -161.867, cumulative = -2421.46
smoothSolver:  Solving for epsilon, Initial residual = 0.13495, Final residual = 0.00391054, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.134138, Final residual = 0.00437452, No Iterations 2
ExecutionTime = 316.77 s  ClockTime = 354 s

Time = 15

smoothSolver:  Solving for Ux, Initial residual = 0.0779784, Final residual = 0.00287186, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0403581, Final residual = 0.00151513, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0772994, Final residual = 0.00287905, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.329316, Final residual = 0.0656906, No Iterations 1
time step continuity errors : sum local = 972.644, global = -160.998, cumulative = -2582.46
smoothSolver:  Solving for epsilon, Initial residual = 0.134529, Final residual = 0.0038749, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.139388, Final residual = 0.00462282, No Iterations 2
ExecutionTime = 338.2 s  ClockTime = 375 s

Time = 16

smoothSolver:  Solving for Ux, Initial residual = 0.075164, Final residual = 0.00285694, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.038141, Final residual = 0.00143983, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0761046, Final residual = 0.00311105, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.325418, Final residual = 0.0581975, No Iterations 1
time step continuity errors : sum local = 870.478, global = -83.0109, cumulative = -2665.47
smoothSolver:  Solving for epsilon, Initial residual = 0.439395, Final residual = 0.0130362, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.552394, Final residual = 0.0347422, No Iterations 2
ExecutionTime = 359.51 s  ClockTime = 402 s

Time = 17

smoothSolver:  Solving for Ux, Initial residual = 0.464726, Final residual = 0.00644288, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.129872, Final residual = 0.00153488, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.773062, Final residual = 0.0125177, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.720437, Final residual = 0.12145, No Iterations 1
time step continuity errors : sum local = 4831.8, global = -1203.49, cumulative = -3868.96
smoothSolver:  Solving for epsilon, Initial residual = 0.99852, Final residual = 0.0378969, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.99975, Final residual = 0.000130416, No Iterations 4
ExecutionTime = 381.33 s  ClockTime = 424 s

Time = 18

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.00994666, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0101358, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0100816, No Iterations 4
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.28906, No Iterations 2
time step continuity errors : sum local = 1.85649e+11, global = -7.84046e+10, cumulative = -7.84047e+10
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 0.0601755, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0695155, No Iterations 2
ExecutionTime = 403.73 s  ClockTime = 446 s

Time = 19

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.00701573, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.00701791, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0070213, No Iterations 4
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.203189, No Iterations 2
time step continuity errors : sum local = 2.55152e+32, global = 1.23564e+32, cumulative = 1.23564e+32
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 0.0624321, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0115553, No Iterations 2
ExecutionTime = 426 s  ClockTime = 468 s

Time = 20

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.00687755, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.00687755, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.00687755, No Iterations 4
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.246173, No Iterations 7
time step continuity errors : sum local = 3.63847e+82, global = 9.72572e+81, cumulative = 9.72572e+81
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 0.00747124, No Iterations 8
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 5.54463e-06, No Iterations 2
ExecutionTime = 454.86 s  ClockTime = 497 s

when it run to steps21,the errors take placed.

[norkistar]

your first step : time step already diverge, you can have a check about your mesh and your boundary condition settings

[jiaojiao]

i run checkMesh and it shows:
Create polyMesh for time = 0

Time = 0

Mesh stats
points: 644339
faces: 6657993
internal faces: 6626605
cells: 3261204
faces per cell: 4.07353
boundary patches: 4
point zones: 0
face zones: 1
cell zones: 2

Overall number of cells of each type:
hexahedra: 100740
prisms: 3984
wedges: 0
pyramids: 34318
tet wedges: 0
tetrahedra: 3122162
polyhedra: 0

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
blade 24146 21942 ok (closed singly connected)
inlet 1841 966 ok (non-closed singly connected)
outlet 1841 966 ok (non-closed singly connected)
symmetry 3560 3649 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-12.7106 -6.66098 -12.7139) (12.7102 27.7695 12.7057)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-2.3031e-15 -3.09845e-16 5.48772e-16) OK.
Max cell openness = 1.07338e-15 OK.
Max aspect ratio = 30.0527 OK.
Minimum face area = 6.8698e-07. Maximum face area = 0.772395. Face area magnitudes OK.
Min volume = 3.21357e-10. Max volume = 0.154152. Total volume = 17466.7. Cell volumes OK.
Mesh non-orthogonality Max: 70.2416 average: 19.945
*Number of severely non-orthogonal (> 70 degrees) faces: 3.
Non-orthogonality check OK.
<<Writing 3 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 3.90819 OK.
Coupled point location match (average 0) OK.

Mesh OK.

how ablout the mesh quality? if bad,does the problem lie in the volume range or the 3 non-orthogonal faces or something eles?

[jiaojiao]

Originally Posted by norkistar

your first step : time step already diverge, you can have a check about your mesh and your boundary condition settings

Hi,I now use Gambit to generate mesh, boundary condition doesnt change and it can run more than 150 steps with the first time steps continuty errors OK. I paste the log file,could you tell me where the problem is?the eorror messege is

#0  Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#2  __restore_rt at sigaction.c:0
#3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#4  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#8  Foam::fvMatrix<double>::solve() in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
#9  main in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
#10  __libc_start_main in "/lib64/libc.so.6"
#11  Foam::UOPstream::write(char) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
娴偣鏁颁緥澶

[jiaojiao]

It seems the log file didnt upload, I am so confused. I can only cut the part where error happens.
Time = 75

smoothSolver: Solving for Ux, Initial residual = 0.0602954, Final residual = 0.00373723, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0653118, Final residual = 0.00486974, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0597028, Final residual = 0.00373956, No Iterations 2
GAMG: Solving for p, Initial residual = 9.91772e-08, Final residual = 9.91772e-08, No Iterations 0
time step continuity errors : sum local = 0.508126, global = -0.00403281, cumulative = 0.0725249
smoothSolver: Solving for epsilon, Initial residual = 0.111438, Final residual = 2.47159e-07, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.0644941, Final residual = 0.00196909, No Iterations 2
ExecutionTime = 1101.44 s ClockTime = 1104 s

Time = 76

smoothSolver: Solving for Ux, Initial residual = 0.0594781, Final residual = 0.00349262, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.060454, Final residual = 0.00433672, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0589499, Final residual = 0.00349089, No Iterations 2
GAMG: Solving for p, Initial residual = 9.85788e-08, Final residual = 9.85788e-08, No Iterations 0
time step continuity errors : sum local = 0.536801, global = -0.00461372, cumulative = 0.0679112
smoothSolver: Solving for epsilon, Initial residual = 0.119377, Final residual = 1.54768e-07, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.0767105, Final residual = 0.00221075, No Iterations 2
ExecutionTime = 1115.63 s ClockTime = 1118 s

Time = 77

smoothSolver: Solving for Ux, Initial residual = 0.0592163, Final residual = 0.00329118, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0563358, Final residual = 0.00385813, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0588684, Final residual = 0.00329036, No Iterations 2
GAMG: Solving for p, Initial residual = 9.65561e-08, Final residual = 9.65561e-08, No Iterations 0
time step continuity errors : sum local = 0.565758, global = -0.00525311, cumulative = 0.0626581
smoothSolver: Solving for epsilon, Initial residual = 0.124077, Final residual = 9.12524e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.0872853, Final residual = 0.00261114, No Iterations 2
ExecutionTime = 1129.82 s ClockTime = 1133 s

Time = 78

smoothSolver: Solving for Ux, Initial residual = 0.0601651, Final residual = 0.00311845, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0534155, Final residual = 0.00342506, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0595368, Final residual = 0.0031283, No Iterations 2
GAMG: Solving for p, Initial residual = 9.28703e-08, Final residual = 9.28703e-08, No Iterations 0
time step continuity errors : sum local = 0.594105, global = -0.00596336, cumulative = 0.0566947
smoothSolver: Solving for epsilon, Initial residual = 0.130123, Final residual = 5.31763e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.0996443, Final residual = 0.00317308, No Iterations 2
bounding k, min: -27.2023 max: 6.08394e+06 average: 47.9263
ExecutionTime = 1144.26 s ClockTime = 1147 s

Time = 79

smoothSolver: Solving for Ux, Initial residual = 0.0614544, Final residual = 0.00297322, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0514696, Final residual = 0.00306002, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0607817, Final residual = 0.00298804, No Iterations 2
GAMG: Solving for p, Initial residual = 8.84393e-08, Final residual = 8.84393e-08, No Iterations 0
time step continuity errors : sum local = 0.620817, global = -0.00674595, cumulative = 0.0499488
smoothSolver: Solving for epsilon, Initial residual = 0.13737, Final residual = 2.93949e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.113938, Final residual = 0.00343359, No Iterations 2
bounding k, min: -216.755 max: 8.91232e+06 average: 61.9746
ExecutionTime = 1158.68 s ClockTime = 1161 s

Time = 80

smoothSolver: Solving for Ux, Initial residual = 0.063004, Final residual = 0.00283814, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0505971, Final residual = 0.00273304, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0624865, Final residual = 0.0028546, No Iterations 2
GAMG: Solving for p, Initial residual = 8.29655e-08, Final residual = 8.29655e-08, No Iterations 0
time step continuity errors : sum local = 0.645584, global = -0.00760253, cumulative = 0.0423462
smoothSolver: Solving for epsilon, Initial residual = 0.143916, Final residual = 1.66331e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.127062, Final residual = 0.00356757, No Iterations 2
bounding k, min: -19.5544 max: 1.31047e+07 average: 81.9192
ExecutionTime = 1173.05 s ClockTime = 1176 s

[Anamitra]

Hi I am using simpleIbFoam, each time I run it shows this error. The checkMesh looks ok. Can someone help with this issue ? Thanks in advance !!

anamitra@anamitra-HP-Pavilion-g6-Notebook-PC:~/IC_Files/Empty_Domain/simpleibfoam_vg$ simpleIbFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.1-7d8e040bf53d
Exec : simpleIbFoam
Date : Apr 21 2015
Time : 19:30:14
Host : anamitra-HP-Pavilion-g6-Notebook-PC
PID : 8630
CtrlDict : /home/anamitra/foam/foam-extend-3.1/etc/controlDict
Case : /home/anamitra/IC_Files/Empty_Domain/simpleibfoam_vg
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
External flow
Number of IB cells: 76
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
c1 10;
}

Create immersed boundary cell mask
Create immersed boundary face mask
Found immersed boundary patch 0 named ibsimplevg

Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 0.357069, Final residual = 0.000767327, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.999927, Final residual = 0.00105008, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.285562, Final residual = 0.000466337, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.984177, Final residual = 0.0909034, No Iterations 85
time step continuity errors : sum local = 0.175394, global = -0.00155391, cumulative = -0.00155391
Floating point exception (core dumped)
anamitra@anamitra-HP-Pavilion-g6-Notebook-PC:~/IC_Files/Empty_Domain/simpleibfoam_vg$

[kamakura117]

Anamitra:

It looks like you need to double check your BCs and ICs.

Also, it's taking a lot of iterations to converge your pressure, you might consider a different numerical method if your BC/ICs are okay.

[Anamitra]

Hi,

Thanks for the reply!! Could you suggest what changes are possible in the boundary conditions, I am still looking for it.

I am not very sure what factors can cause the floating point exception. I am solving a case of external flow around a cuboid (with Immersed boundary) placed around the centre of the flow domain.

Pressure looks like this:
internalField nonuniform List<scalar>
256000
(long list found from mapping fields)

boundaryField
{
ibsimplevg
{
type immersedBoundary;
refValue uniform 0;
refGradient uniform 0;
fixesValue no;
setDeadCellValue yes;
deadCellValue 0;
value nonuniform 0();
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
sides
{
type symmetryPlane;
}
wall_plate
{
type zeroGradient;
}
top
{
type zeroGradient;
}
}

[Kittipos]

Hi Foamer

I have an error with running rhoPorousSimpleFoam in my simple model.

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so"
#4  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so"
#5  
 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam"
#6  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam"
Floating point exception (core dumped)

Thank you