[callumso]

Hi Foamers,

I am running a series of similar cases with SonicFoam. Most of them worked fine. But a handful of them just blew up all of a sudden when everything still looked rather normal at the previous time step.

Below is the message printed out for one of the case:

Time = 0.51623

Courant Number mean: 0.000144416 max: 4.330942
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
3422 3423
DILUPBiCG: Solving for Ux, Initial residual = 1.461462e-06, Final residual = 2.028677e-11, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 0.0002017862, Final residual = 9.653597e-12, No Iterations 12
DILUPBiCG: Solving for Uz, Initial residual = 0.0002056769, Final residual = 3.897662e-11, No Iterations 11
DILUPBiCG: Solving for e, Initial residual = 2.438534e-05, Final residual = 4.253788e-11, No Iterations 7
DILUPBiCG: Solving for p, Initial residual = 1.931868e-05, Final residual = 9.061885e-11, No Iterations 81
DILUPBiCG: Solving for p, Initial residual = 2.405003e-07, Final residual = 2.852968e-11, No Iterations 67
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.232529e-13, global = 1.979396e-15, cumulative = -1.438724e-11
3423 3423
DILUPBiCG: Solving for p, Initial residual = 8.452907e-07, Final residual = 3.624389e-11, No Iterations 72
DILUPBiCG: Solving for p, Initial residual = 1.477041e-07, Final residual = 6.175933e-11, No Iterations 51
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.81663e-13, global = -3.274056e-14, cumulative = -1.441998e-11
3423 3423
DILUPBiCG: Solving for p, Initial residual = 5.465708e-07, Final residual = 6.771115e-11, No Iterations 61
DILUPBiCG: Solving for p, Initial residual = 7.216307e-08, Final residual = 7.183256e-11, No Iterations 54
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.598445e-13, global = 1.141557e-14, cumulative = -1.440856e-11
3423 3423
DILUPBiCG: Solving for epsilon, Initial residual = 3.243357e-07, Final residual = 6.565833e-09, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 1.069664e-06, Final residual = 6.12174e-09, No Iterations 3
ExecutionTime = 38608.39 s ClockTime = 39120 s

Time = 0.51624

Courant Number mean: 0.0001443969 max: 4.323677
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
3423 3424
DILUPBiCG: Solving for Ux, Initial residual = 1.44806e-06, Final residual = 4.422659e-12, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 0.0002031021, Final residual = 9.751295e-11, No Iterations 10
DILUPBiCG: Solving for Uz, Initial residual = 0.0002066618, Final residual = 4.271786e-11, No Iterations 11
DILUPBiCG: Solving for e, Initial residual = 2.441768e-05, Final residual = 3.836755e-11, No Iterations 7
DILUPBiCG: Solving for p, Initial residual = 1.933993e-05, Final residual = 7.494482e-11, No Iterations 78
DILUPBiCG: Solving for p, Initial residual = 2.216237e-07, Final residual = 8.573657e-11, No Iterations 51
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.671348e-13, global = -2.6103e-14, cumulative = -1.443467e-11
3424 3424
DILUPBiCG: Solving for p, Initial residual = 9.442638e-07, Final residual = 7.056143e-11, No Iterations 67
DILUPBiCG: Solving for p, Initial residual = 1.311588e-07, Final residual = 9.744212e-11, No Iterations 47
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.58383e-13, global = -5.861799e-14, cumulative = -1.449328e-11
3424 3424
DILUPBiCG: Solving for p, Initial residual = 7.544892e-07, Final residual = 8.627015e-11, No Iterations 61
DILUPBiCG: Solving for p, Initial residual = 6.800743e-08, Final residual = 3.908017e-11, No Iterations 49
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.050143e-13, global = -2.038326e-15, cumulative = -1.449532e-11
3424 3424
DILUPBiCG: Solving for epsilon, Initial residual = 3.246423e-07, Final residual = 7.913151e-09, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 1.069163e-06, Final residual = 6.115701e-09, No Iterations 3
ExecutionTime = 38634.09 s ClockTime = 39146 s

Time = 0.51625

Courant Number mean: 0.0001443834 max: 4.323903
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
3424 3425
DILUPBiCG: Solving for Ux, Initial residual = 1.44179e-06, Final residual = 2.484585e-11, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 0.0002025631, Final residual = 1.016718e+17, No Iterations 1001
DILUPBiCG: Solving for Uz, Initial residual = 0.0002070403, Final residual = 4.698542e-11, No Iterations 11
DILUPBiCG: Solving for e, Initial residual = 2.445833e-05, Final residual = 4.306151e-11, No Iterations 7
DILUPBiCG: Solving for p, Initial residual = 0.9975288, Final residual = 2.605364, No Iterations 1001
DILUPBiCG: Solving for p, Initial residual = 0.5460555, Final residual = 16.64415, No Iterations 1001
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = -6.012919e+14, global = 1.057316e+11, cumulative = 1.057316e+11
3425 3425
DILUPBiCG: Solving for p, Initial residual = 0.9969977, Final residual = 0.03398452, No Iterations 1001
DILUPBiCG: Solving for p, Initial residual = 3.084765e-07, Final residual = 1.189066e-07, No Iterations 1001
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = -3.34374e+15, global = 3379692, cumulative = 1.05735e+11
3425 3425
DILUPBiCG: Solving for p, Initial residual = 0.9996619, Final residual = 1.179515, No Iterations 1001
DILUPBiCG: Solving for p, Initial residual = 5.310774e-14, Final residual = 5.310774e-14, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = -9.701425e+14, global = 1, cumulative = 1.05735e+11
3425 3425
DILUPBiCG: Solving for epsilon, Initial residual = 0.9997756, Final residual = 50231.57, No Iterations 1001
bounding epsilon, min: -1.416225e+72 max: 6.146947e+70 average: -3.553829e+66
DILUPBiCG: Solving for k, Initial residual = 0.01615606, Final residual = 8.162019e-16, No Iterations 2
bounding k, min: -3.304356e+66 max: 7.639573e+63 average: -5.087817e+60
ExecutionTime = 39004.48 s ClockTime = 39517 s

Time = 0.51626

Courant Number mean: 2.47646e+43 max: 1.355166e+53
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
3425 3426
[0] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/lib64/libc.so.6"
[1] #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/lib64/libc.so.6"
[0] #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #8 in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const
in "/home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #8 [1] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[1] #9

[1] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[1] #10 [0] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[0] #9

[1] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[1] #11 __libc_start_main[0] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[0] #10 in "/lib64/libc.so.6"
[1] #12

[0] in "/home/ttjs3/OpenFOAM/ttjs3-2.1.1/platforms/linux64GccDPOpt/bin/own_sonicFoam"
[0] #11 __libc_start_main[1] at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
[neso:23806] *** Process received signal ***
[neso:23806] Signal: Floating point exception (8)
[neso:23806] Signal code: (-6)
[neso:23806] Failing at address: 0x6ac900005cfe
[neso:23806] [ 0] /lib64/libc.so.6(+0x32bd0) [0x7f0bd44a5bd0]
[neso:23806] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7f0bd44a5b35]
[neso:23806] [ 2] /lib64/libc.so.6(+0x32bd0) [0x7f0bd44a5bd0]
[neso:23806] [ 3] /home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditioner15calcR eciprocalDERNS_5FieldIdEERKNS_9lduMatrixE+0x7f) [0x7f0bd541831f]
[neso:23806] [ 4] /home/ttjs3/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditionerC1ERKNS _9lduMatrix6solverERKNS_10dictionaryE+0x14e) [0x7f0bd541856e]
....


As it shows, at "Time = 0.51624" everything looked perfectly normal. But at "Time = 0.51625", the Uy equation just suddently ran into trouble. I have no clue how this happened. Below is the fvSolution file:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
"(p|U|e)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-10;
relTol 0;
}

rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0;
}

"(k|epsilon|R)"
{
$p;
tolerance 1e-8;
relTol 0;
}
}

PISO
{
nCorrectors 3;
nNonOrthogonalCorrectors 1;
}

And the fvSchemes file is:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

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

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phid,p) Gauss limitedLinear 1;
div(phi,K) Gauss limitedLinear 1;
div(phi,e) Gauss upwind;
div((muEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear limited 0.5;
laplacian(DkEff,k) Gauss linear limited 0.5;
laplacian(DREff,R) Gauss linear limited 0.5;
laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
laplacian((rho*(1|A(U))),p) Gauss linear limited 0.5;
laplacian(alphaEff,e) Gauss linear limited 0.5;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}

The controlDict file is:

application sonicFoam;

startFrom startTime;

startTime 0.5;

stopAt endTime;

endTime 0.65;

deltaT 0.00001;

writeControl adjustableRunTime;

writeInterval 0.00125;

purgeWrite 0;

writeFormat ascii;

writePrecision 7;

writeCompression off;

timeFormat general;

timePrecision 7;

runTimeModifiable true;

adjustTimeStep no;

maxCo 10;

maxDeltaT 1;

functions
(
probes1
{
type probes; // Type of functionObject
// Where to load it from (if not already in solver)
functionObjectLibs ("libsampling.so");
probeLocations // Locations to be probed. runTime modifiable!
(
(-0.1754 0.0 0.0)
(0.0 0.0 0.0)
(0.0455 0.0 0.0)
(0.2391 0.0 0.0)
(0.3831 0.0 0.0)
(0.5231 0.0 0.0)
(0.6241 0.0 0.0)
(0.7081 0.0 0.0)
(0.2391 0.015 0.0)
(0.2391 0.03 0.0)
(0.3831 0.015 0.0)
(0.3831 0.03 0.0)
);
// Fields to be probed. runTime modifiable!
fields
(
p
U
);
}
);

When I looked at the blown up solution, the entire field just got messed up and I can't tell which bit of the filed has gone haywired.

Is it possible to deduce something from the traces I gave above? Or are there some problems with the linear solvers? Thank you very much.

Regards,
Callum

[callumso]

Seems the problem can be solved by using four processors instead of two (or just try different number of partitions). something not unusual even for a commercial package. not sure this is worth reporting to the OpenFOAM team?