|
[Sponsors] |
Foam::error::printStack in interDyMFoam when changing the floatingObject case |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 6, 2014, 03:53 |
Foam::error::printStack in interDyMFoam when changing the floatingObject case
|
#1 |
New Member
Join Date: Jun 2014
Posts: 22
Rep Power: 11 |
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) |
|
August 6, 2014, 04:13 |
|
#2 |
New Member
Join Date: Jun 2014
Posts: 22
Rep Power: 11 |
I think I know the reason.
The Courant Number was too large. It begin to run after I reduce the delta T. |
|
August 7, 2014, 02:53 |
|
#3 |
Member
Join Date: May 2013
Posts: 51
Rep Power: 12 |
Two weeks ealier I had similar problem. I have found solution for this :
Code:
unset FOAM_SIGFPE seav |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Changing the grid on the same set-up | Katya | FLUENT | 7 | October 8, 2009 16:31 |
Free surface boudary conditions with SOLA-VOF | Fan | Main CFD Forum | 10 | September 9, 2006 12:24 |
Changing gambit file without change of case file?? | Asghari | FLUENT | 2 | August 28, 2006 13:48 |
Turbulent Flat Plate Validation Case | Jonas Larsson | Main CFD Forum | 0 | April 2, 2004 10:25 |
Body force - Does it work? | Jan Rusås | CFX | 5 | August 27, 2002 09:50 |