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

SonicFoam crashes without any ominous signs

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 23, 2013, 07:47
Default SonicFoam crashes without any ominous signs
  #1
New Member
 
Jialin Su
Join Date: Mar 2013
Location: Loughborough
Posts: 29
Rep Power: 13
callumso is on a distinguished road
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 is offline   Reply With Quote

Old   April 25, 2013, 06:45
Default
  #2
New Member
 
Jialin Su
Join Date: Mar 2013
Location: Loughborough
Posts: 29
Rep Power: 13
callumso is on a distinguished road
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?
callumso is offline   Reply With Quote

Reply


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
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 1 September 24, 2019 01:15
Error trying to run steady-state sonicFoam dancfd OpenFOAM Running, Solving & CFD 2 February 12, 2013 04:15
flo-efd v11.0.0 crashes YoavF FloEFD, FloWorks & FloTHERM 3 June 21, 2012 13:37
Parallel runs with sonicDyMFoam crashes (works fine with sonicFoam) jnilsson OpenFOAM Running, Solving & CFD 0 March 9, 2012 07:45
Crash when using sonicFoam Horus OpenFOAM 1 June 16, 2010 13:57


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