pimpelDyMfoam suddenly blow up (velocity)
 

Hello Everyone

I am using PimpeDymesh with AMI to simulate a rotating wheel, and the time step is 2e-5 , all the schemes are first order, and until t=0.82893 everything is fun, but from the next step t=0.82895, it suddenly blow up.

anyone had similar problem or some tipps

checkMesh:

Checking geometry...
Overall domain bounding box (-3 -3 -0.125228) (8 3 3.00002)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-7.5053e-15 5.60214e-16 -1.12528e-14) OK.
Max cell openness = 3.79016e-16 OK.
Max aspect ratio = 16.1504 OK.
Minimum face area = 8.61023e-09. Maximum face area = 0.0194853. Face area magnitudes OK.
Min volume = 2.94566e-11. Max volume = 0.00265095. Total volume = 206.227. Cell volumes OK.
Mesh non-orthogonality Max: 59.9619 average: 4.93528
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.97977 OK.
Coupled point location match (average 0) OK.

Mesh OK.


moveDynamicMesh -checkAMI

log file of pimpleDymFoam
Time = 0.0068
solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 0.0068 transformation: ((0.0059926 0 0.000961795) (0.376604 (0 -0.926374 0)))
AMI: Creating addressing and weights between 34240 source faces and 34240 target faces
AMI: Patch source sum(weights) min/max/average = 0.422807, 1.15043, 0.998492
AMI: Patch target sum(weights) min/max/average = 0.334496, 1.09189, 0.998241
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Mesh topology OK.
Boundary openness (-4.20596e-16 6.0293e-16 -1.87871e-13) OK.
Max cell openness = 3.79016e-16 OK.
Max aspect ratio = 16.1504 OK.
Minimum face area = 8.61023e-09. Maximum face area = 0.0194853. Face area magnitudes OK.
Min volume = 2.94566e-11. Max volume = 0.00265095. Total volume = 206.227. Cell volumes OK.
Mesh non-orthogonality Max: 59.9619 average: 4.93528
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.97977 OK.
Mesh geometry OK.
Mesh OK.
Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2
ExecutionTime = 21617.4 s ClockTime = 21842 s

the mesh looks good, and the schemes and solver work very well for other simple mesh as well, so I think it should be some problem with the mesh...but check mesh and movedynamikmesh are all fine.



Time = 0.82895

solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 0.82895 transformation: ((-0.000274797 0 0.000571264) (0.995308 (0 0.0967573 0)))
AMI: Creating addressing and weights between 34240 source faces and 34240 target faces
AMI: Patch source sum(weights) min/max/average = 0.618387, 1.11881, 0.999584
AMI: Patch target sum(weights) min/max/average = 0.621077, 1.17541, 0.99958
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 3.8942e-06, Final residual = 7.78869e-12, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 4.89506e-05, Final residual = 2.59989e-12, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 3.73426e-05, Final residual = 1.55752e-11, No Iterations 4
GAMGPCG: Solving for p, Initial residual = 0.265868, Final residual = 0.000483203, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.00707617, Final residual = 1.47055e-05, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.000799827, Final residual = 1.27174e-06, No Iterations 2
time step continuity errors : sum local = 5.815e-13, global = 1.20908e-14, cumulative = 4.69174e-11
GAMGPCG: Solving for p, Initial residual = 0.173067, Final residual = 0.000382646, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.00499796, Final residual = 8.93821e-06, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.000573963, Final residual = 9.29454e-07, No Iterations 2
time step continuity errors : sum local = 4.24875e-13, global = 3.03523e-15, cumulative = 4.69204e-11
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 2.89581e-06, Final residual = 1.36012e-12, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 3.61649e-05, Final residual = 1.69233e+27, No Iterations 1001
DILUPBiCG: Solving for Uz, Initial residual = 2.7915e-05, Final residual = 1.50301e-11, No Iterations 5
GAMGPCG: Solving for p, Initial residual = 1, Final residual = 0.00120953, No Iterations 3
GAMGPCG: Solving for p, Initial residual = 0.0880561, Final residual = 0.000323297, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.00996636, Final residual = 6.33094e-05, No Iterations 2
time step continuity errors : sum local = 5.83761e+20, global = -5.3204e+19, cumulative = -5.3204e+19
GAMGPCG: Solving for p, Initial residual = 0.513317, Final residual = 0.00056665, No Iterations 3
GAMGPCG: Solving for p, Initial residual = 0.292931, Final residual = 0.000974897, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.0395915, Final residual = 4.18906e-07, No Iterations 5
time step continuity errors : sum local = 8.56991e+17, global = -8.10343e+16, cumulative = -5.3285e+19
DILUPBiCG: Solving for nuTilda, Initial residual = 0.373567, Final residual = 0.0108753, No Iterations 3
ExecutionTime = 171134 s ClockTime = 172533 s