CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Foam::error::printStack in interDyMFoam when changing the floatingObject case (https://www.cfd-online.com/Forums/openfoam-solving/140012-foam-error-printstack-interdymfoam-when-changing-floatingobject-case.html)

chun August 6, 2014 03:53

Foam::error::printStack in interDyMFoam when changing the floatingObject case
 
Hello everyone,

I am trying to change the floatingObject case with interDyMFoam by openFoam 2.3.0.
The only change that I made was simply stick a turbine into the box. (I use snappyHexMesh to get the mesh done)

I can get the interDyMFoam running, but it get the error message very quick.

Code:

Interface Courant Number mean: 0.00845483025063 max: 46.9611791733
Courant Number mean: 0.179239126357 max: 85.6924669702
deltaT = 0.00012957563978
Time = 0.0112407

PIMPLE: iteration 1

Centre of mass: (0.500036791129 0.499459476047 0.410074192546)
Linear velocity: (0.00809747693356 0.0050168267431 -14.6817409152)
Angular velocity: (-0.761780378942 -0.0385987965828 -0.0271900008197)
Execution time for mesh.update() = 2.37 s
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 5.95903194681e-06, No Iterations 7
time step continuity errors : sum local = 2.9431910746e-09, global = 1.418932976e-09, cumulative = -0.0002136723294
smoothSolver:  Solving for alpha.water, Initial residual = 0.00109973974403, Final residual = 5.06941305157e-09, No Iterations 5
Phase-1 volume fraction = 0.44890634435  Min(alpha1) = -0.00134888936293  Max(alpha1) = 1.04137906511
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.44890634435  Min(alpha1) = -0.00167727959902  Max(alpha1) = 1.30118612275
GAMG:  Solving for p_rgh, Initial residual = 0.0233445876891, Final residual = 0.000127317324022, No Iterations 7
time step continuity errors : sum local = 3.73830976879e-06, global = -1.31450522516e-06, cumulative = -0.000214986834625
PIMPLE: iteration 2

Centre of mass: (0.500034286731 0.499471738212 0.411044501906)
Linear velocity: (-0.0141126916462 0.0961606865696 1.73393838823)
Angular velocity: (0.11494883232 -0.0435266761614 -0.012933016184)
Execution time for mesh.update() = 2.35 s
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 1.02556844088e-06, No Iterations 7
time step continuity errors : sum local = 2.56111418526e-10, global = 1.68491561144e-10, cumulative = -0.000214986666134
smoothSolver:  Solving for alpha.water, Initial residual = 0.000699357606737, Final residual = 2.36261277552e-09, No Iterations 3
Phase-1 volume fraction = 0.44890634933  Min(alpha1) = -0.00964830870267  Max(alpha1) = 1.05026732966
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.44890634933  Min(alpha1) = -0.0157291914946  Max(alpha1) = 1.05026732966
GAMG:  Solving for p_rgh, Initial residual = 0.000147432618839, Final residual = 8.16690647211e-07, No Iterations 5
time step continuity errors : sum local = 1.09685101299e-05, global = 3.68346223864e-06, cumulative = -0.000211303203895
PIMPLE: iteration 3

Centre of mass: (0.500002259257 0.499432967171 0.510236910937)
Linear velocity: (-0.256104781196 -0.224251359964 765.011355267)
Angular velocity: (25.4514739697 -0.815258408664 0.134150329908)
Execution time for mesh.update() = 2.36 s
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 1.25878373672e-06, No Iterations 7
time step continuity errors : sum local = 2.81810951347e-08, global = 1.19571609549e-08, cumulative = -0.000211291246734
smoothSolver:  Solving for alpha.water, Initial residual = 0.0526947575188, Final residual = 8.09904568507e-09, No Iterations 78
Phase-1 volume fraction = 0.448906355963  Min(alpha1) = -0.00146272212247  Max(alpha1) = 1.0000000104
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.448906355964  Min(alpha1) = -15.7709517047  Max(alpha1) = 14.1515141827
GAMG:  Solving for p_rgh, Initial residual = 0.00527085075396, Final residual = 1.85424905326e+29, No Iterations 1000
time step continuity errors : sum local = -2.76544388874e+29, global = 5.76168746684e+29, cumulative = 5.76168746684e+29
PIMPLE: iteration 4

Centre of mass: (0.500059814133 0.499197450026 0.485348024118)
Linear velocity: (5.45056704154e+30 3.14946176249e+30 1.03609251921e+32)
Angular velocity: (-4.11852589633e+31 3.4852527546e+31 -2.45766351592e+30)
Execution time for mesh.update() = 2.31 s
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 1.00003481069, No Iterations 101
time step continuity errors : sum local = -5.31265459993e+28, global = 3.83835740941e+29, cumulative = 9.60004487625e+29
smoothSolver:  Solving for alpha.water, Initial residual = 0.03657562748, Final residual = 9.24010741808e-09, No Iterations 62
Phase-1 volume fraction = 4.72845877569e+21  Min(alpha1) = -1.06464462909e-32  Max(alpha1) = 1.29130250216e+30
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 4.72845877571e+21  Min(alpha1) = -1.20322904843e+30  Max(alpha1) = 1.24206487342e+30
#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::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  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::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::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 
 at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 
 at ??:?

Floating point exception (core dumped)

Can anyone tell me what's the error message means? Thank you so much.

chun August 6, 2014 04:13

I think I know the reason.

The Courant Number was too large.
It begin to run after I reduce the delta T.

seav August 7, 2014 02:53

Two weeks ealier I had similar problem. I have found solution for this :

Code:

unset FOAM_SIGFPE
Try this before u start solver, maybe it will help.


seav :)


All times are GMT -4. The time now is 02:38.