interDyMFoam "time step continuity errors"
Hello all,
I am having problems with my "waveDyMFoam" which I modified from "interDyMFoam". I am trying to run a simple 2D case where a box is floating on the water under the effect of waves. The video of the case can be seen here:
https://www.youtube.com/watch?featur...&v=YKbj_7JMRl8
Also the case can be downloaded from here:
https://sites.google.com/site/jordim...edirects=0&d=1
At first the case was not running because of some missing schemes in the "fvSchemes" file. So, I added the following lines under "divSchemes" section:
div((muEff*dev(T(grad(U))))) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
When I try to run the case, I get the following errors:
Quote:
sixDoFRigidBodyMotion constraints converged in 1 iterations
Constraint force: (0 0 0)
Constraint moment: (0 0 0)
Centre of mass: (10.3657 0.00896592 0.05)
Linear velocity: (0.169457 0.0262615 3.6077e17)
Angular velocity: (3.19807e28 1.37522e26 5.1003e12)
GAMG: Solving for cellDisplacementx, Initial residual = 0.000224045, Final residual = 8.12868e06, No Iterations 2
GAMG: Solving for cellDisplacementy, Initial residual = 0.000224045, Final residual = 8.12868e06, No Iterations 2
Execution time for mesh.update() = 0.23 s
time step continuity errors : sum local = 5.05277e05, global = 3.01844e19, cumulative = 5.40294e08
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 8.97864e08, No Iterations 17
GAMG: Solving for pcorr, Initial residual = 0.250457, Final residual = 5.28244e08, No Iterations 15
time step continuity errors : sum local = 8.53279e05, global = 1.02447e11, cumulative = 5.40397e08
MULES: Solving for alpha1
Phase1 volume fraction = 0.666764 Min(alpha1) = 1.93586e23 Max(alpha1) = 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0130186, Final residual = 1.17418e10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0946691, Final residual = 7.58036e11, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.261841, Final residual = 5.56369e08, No Iterations 12
GAMG: Solving for p_rgh, Initial residual = 0.000997295, Final residual = 4.38728e08, No Iterations 7
time step continuity errors : sum local = 7.35319e21, global = 4.08938e21, cumulative = 5.40397e08
GAMG: Solving for p_rgh, Initial residual = 0.000600283, Final residual = 5.20561e08, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 0.00035003, Final residual = 4.53092e08, No Iterations 5
time step continuity errors : sum local = 7.67551e21, global = 6.55771e21, cumulative = 5.40397e08
GAMG: Solving for p_rgh, Initial residual = 0.000171914, Final residual = 7.8525e08, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 5.06367e05, Final residual = 6.19844e09, No Iterations 6
time step continuity errors : sum local = 1.04979e21, global = 1.00984e21, cumulative = 5.40397e08
ExecutionTime = 47953.3 s ClockTime = 48046 s
Interface Courant Number mean: 2.38452e06 max: 0.00400263
Courant Number mean: 0.000200405 max: 0.425284
deltaT = 2.62873e13
> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 14 to 16 to distinguish between timeNames at time 12.8093
Time = 12.80933457039102
sixDoFRigidBodyMotion constraints converged in 1 iterations
Constraint force: (0 0 0)
Constraint moment: (0 0 0)
Centre of mass: (10.3657 0.00896592 0.05)
Linear velocity: (0.169262 0.026284 3.60775e17)
Angular velocity: (3.19814e28 1.37523e26 5.10294e12)
GAMG: Solving for cellDisplacementx, Initial residual = 0.000221903, Final residual = 8.00376e06, No Iterations 2
GAMG: Solving for cellDisplacementy, Initial residual = 0.000221903, Final residual = 8.00376e06, No Iterations 2
Execution time for mesh.update() = 0.23 s
time step continuity errors : sum local = 5.01594e05, global = 1.26375e19, cumulative = 5.40397e08
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 8.65576e08, No Iterations 18
GAMG: Solving for pcorr, Initial residual = 0.181331, Final residual = 5.39143e08, No Iterations 15
time step continuity errors : sum local = 7.78955e05, global = 1.19518e11, cumulative = 5.40516e08
MULES: Solving for alpha1
Phase1 volume fraction = 0.666764 Min(alpha1) = 1.88721e23 Max(alpha1) = 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0129951, Final residual = 6.3627e11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0946374, Final residual = 9.16735e11, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.205321, Final residual = 7.70156e08, No Iterations 11
GAMG: Solving for p_rgh, Initial residual = 0.000734827, Final residual = 3.76525e08, No Iterations 7
time step continuity errors : sum local = 3.98678e21, global = 2.59775e22, cumulative = 5.40516e08
GAMG: Solving for p_rgh, Initial residual = 0.00041682, Final residual = 7.59192e08, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.000186038, Final residual = 4.30575e08, No Iterations 5
time step continuity errors : sum local = 4.24988e21, global = 3.26612e21, cumulative = 5.40516e08
GAMG: Solving for p_rgh, Initial residual = 0.000101912, Final residual = 4.99173e08, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 3.9727e05, Final residual = 8.3358e09, No Iterations 6
time step continuity errors : sum local = 8.22266e22, global = 6.35014e22, cumulative = 5.40516e08
ExecutionTime = 47954.4 s ClockTime = 48047 s
Interface Courant Number mean: 2.30322e06 max: 0.0032262
Courant Number mean: 0.000207624 max: 0.544706
deltaT = 1.20649e13
Time = 12.80933457039114
sixDoFRigidBodyMotion constraints converged in 1 iterations
Constraint force: (0 0 0)
Constraint moment: (0 0 0)
Centre of mass: (10.3657 0.00896592 0.05)
Linear velocity: (0.169412 0.0262994 3.60779e17)
Angular velocity: (3.19822e28 1.37522e26 5.10071e12)
GAMG: Solving for cellDisplacementx, Initial residual = 0.00021926, Final residual = 7.84774e06, No Iterations 2
GAMG: Solving for cellDisplacementy, Initial residual = 0.00021926, Final residual = 7.84774e06, No Iterations 2
Execution time for mesh.update() = 0.23 s
time step continuity errors : sum local = 3.57512e05, global = 7.09289e19, cumulative = 5.40516e08
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 7.51115e08, No Iterations 20
GAMG: Solving for pcorr, Initial residual = 0.121942, Final residual = 9.90107e08, No Iterations 14
time step continuity errors : sum local = 8.06169e05, global = 3.4904e11, cumulative = 5.40865e08
#0 Foam::error::printStack(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_64linuxgnu/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6
in "/home/meta/OpenFOAM/root2.2.0/platforms/linux64GccDPOpt/bin/waveDyMFoam"
#7 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6"
#8
in "/home/meta/OpenFOAM/root2.2.0/platforms/linux64GccDPOpt/bin/waveDyMFoam"
Floating point exception (core dumped)

And my analysis crashes around the 12th second. I checked my mesh and everything seemed ok to me.
Do you have any idea what may be the cause of that problem?
Thank you very much for your time and help.
Kilroy
