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

pimpelDyMfoam suddenly blow up (velocity)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 12, 2016, 05:59
Default pimpelDyMfoam suddenly blow up (velocity)
  #1
New Member
 
LU
Join Date: Oct 2015
Posts: 11
Rep Power: 7
louvel is on a distinguished road
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
louvel is offline   Reply With Quote

Reply

Tags
ami patches, diverged, pimpledymfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Problem with Velocity Poisson Equation and Vector Potential Poisson Equation mykkujinu2201 Main CFD Forum 1 August 12, 2017 13:15
Superimposing 2 Velocity fields For Model inputs in Comsol f.black Main CFD Forum 0 August 2, 2017 09:35
atmBoundaryLayerInletVelocity - Velocity Profile not continuous through domain sdfij6354 OpenFOAM Running, Solving & CFD 3 July 26, 2017 16:16
Velocity in Porous medium : HELP! HELP! HELP! Kali Sanjay Phoenics 0 November 6, 2006 06:10
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 05:13


All times are GMT -4. The time now is 07:47.