# calculation diverge after continue to run

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

 April 28, 2010, 04:35 calculation diverge after continue to run #1 Member   jingjing Join Date: Mar 2009 Location: shanghai,China Posts: 30 Rep Power: 10 Sponsored Links Hi，everyone Recently I find that when I continue to run my case(eg, I calculated my case for 20s last time, now I modified the ebdTime to be 40,so I continue running), the calculation diverge. I try several cases, they all fail like this. This phenomenon never happened before. Does anyone meet this problem? you see, // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 25 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: displacementLaplacian water depth 0.6 wave type linearWave wave omega 6.283185 wave height 0.02 wave number 4.084595 w=1.836603 piston stroke 0.005444834 Selecting motion diffusion: inverseDistance Reading environmentalProperties Reading field pd Reading field gamma Reading field U Reading/calculating face flux field phi Reading transportProperties Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting RAS turbulence model laminar time step continuity errors : sum local = 1.444482e-08, global = 1.435022e-08, cumulative = 1.435022e-08 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 7.714203e-08, No Iterations 140 DICPCG: Solving for pcorr, Initial residual = 3.427171e-09, Final residual = 3.427171e-09, No Iterations 0 time step continuity errors : sum local = 1.114302e-15, global = 2.682517e-17, cumulative = 1.435022e-08 Courant Number mean: 0.01505891 max: 0.1339163 Starting time loop Courant Number mean: 0.01505891 max: 0.1339163 Time = 25.002 linear Wave piston position is 4.293759e-07 DICPCG: Solving for cellDisplacementx, Initial residual = 0.9988885, Final residual = 9.096984e-08, No Iterations 383 DICPCG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 1.37 s time step continuity errors : sum local = 1.434755e-08, global = -1.434755e-08, cumulative = 2.671165e-12 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.846891e-08, No Iterations 139 DICPCG: Solving for pcorr, Initial residual = 4.348726e-09, Final residual = 4.348726e-09, No Iterations 0 time step continuity errors : sum local = 1.41279e-15, global = -3.167157e-17, cumulative = 2.671133e-12 MULES: Solving for gamma Liquid phase volume fraction = 0.5000462 Min(gamma) = -2.304434e-67 Max(gamma) = 1.000002 MULES: Solving for gamma Liquid phase volume fraction = 0.5000011 Min(gamma) = -2.297683e-67 Max(gamma) = 1.000004 DICPCG: Solving for pd, Initial residual = 0.1916157, Final residual = 8.835227e-08, No Iterations 84 DICPCG: Solving for pd, Initial residual = 3.197623e-09, Final residual = 3.197623e-09, No Iterations 0 time step continuity errors : sum local = 4.759465e-10, global = 7.660377e-11, cumulative = 7.92749e-11 DICPCG: Solving for pd, Initial residual = 0.0006697087, Final residual = 8.359662e-08, No Iterations 49 DICPCG: Solving for pd, Initial residual = 8.633154e-08, Final residual = 8.633154e-08, No Iterations 0 time step continuity errors : sum local = 1.244683e-08, global = -3.657849e-09, cumulative = -3.578575e-09 ExecutionTime = 4 s ClockTime = 4 s Courant Number mean: 0.06487583 max: 12.65412 Time = 25.004 linear Wave piston position is 1.718486e-06 DICPCG: Solving for cellDisplacementx, Initial residual = 0.8184182, Final residual = 9.771286e-08, No Iterations 382 DICPCG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 1.29 s time step continuity errors : sum local = 1.244683e-08, global = -3.65785e-09, cumulative = -7.236424e-09 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.506461e-08, No Iterations 155 DICPCG: Solving for pcorr, Initial residual = 1.184242e-08, Final residual = 1.184242e-08, No Iterations 0 time step continuity errors : sum local = 1.183273e-15, global = -6.732176e-17, cumulative = -7.236424e-09 MULES: Solving for gamma Liquid phase volume fraction = 0.4999555 Min(gamma) = -2.855421 Max(gamma) = 1.000028 MULES: Solving for gamma Liquid phase volume fraction = 0.4999103 Min(gamma) = -11.62382 Max(gamma) = 12.92754 DICPCG: Solving for pd, Initial residual = 0.01535727, Final residual = 9.311828e-08, No Iterations 69 DICPCG: Solving for pd, Initial residual = 9.157189e-08, Final residual = 9.157189e-08, No Iterations 0 time step continuity errors : sum local = 1.405371e-08, global = -2.960173e-09, cumulative = -1.01966e-08 DICPCG: Solving for pd, Initial residual = 0.01140352, Final residual = 8.428359e-08, No Iterations 47 DICPCG: Solving for pd, Initial residual = 8.423795e-08, Final residual = 8.423795e-08, No Iterations 0 time step continuity errors : sum local = 1.305664e-08, global = 1.663829e-09, cumulative = -8.532768e-09 ExecutionTime = 6.32 s ClockTime = 7 s Courant Number mean: 0.01547009 max: 0.9084346 Time = 25.006 linear Wave piston position is 3.867127e-06 DICPCG: Solving for cellDisplacementx, Initial residual = 0.6518518, Final residual = 7.768823e-08, No Iterations 382 DICPCG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 1.28 s time step continuity errors : sum local = 1.305664e-08, global = 1.663829e-09, cumulative = -6.868939e-09 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.751944e-08, No Iterations 154 DICPCG: Solving for pcorr, Initial residual = 2.320168e-08, Final residual = 2.320168e-08, No Iterations 0 time step continuity errors : sum local = 1.273276e-15, global = 5.851911e-17, cumulative = -6.868939e-09 MULES: Solving for gamma Liquid phase volume fraction = 0.4999103 Min(gamma) = -7.181216 Max(gamma) = 4.842551 MULES: Solving for gamma Liquid phase volume fraction = 0.4999103 Min(gamma) = -6.619966 Max(gamma) = 2.834347 DICPCG: Solving for pd, Initial residual = 0.07596338, Final residual = 8.243381e-08, No Iterations 67 DICPCG: Solving for pd, Initial residual = 9.909663e-07, Final residual = 7.865661e-08, No Iterations 13 time step continuity errors : sum local = 1.07695e-09, global = 3.32348e-10, cumulative = -6.536591e-09 DICPCG: Solving for pd, Initial residual = 0.0440518, Final residual = 9.535207e-08, No Iterations 53 DICPCG: Solving for pd, Initial residual = 8.70003e-08, Final residual = 8.70003e-08, No Iterations 0 time step continuity errors : sum local = 1.295907e-09, global = -4.037227e-10, cumulative = -6.940313e-09 ExecutionTime = 8.72 s ClockTime = 9 s Courant Number mean: 0.0212218 max: 7.507377 Time = 25.008 linear Wave piston position is 6.87496e-06 DICPCG: Solving for cellDisplacementx, Initial residual = 0.5380525, Final residual = 8.313747e-08, No Iterations 381 DICPCG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 1.27 s time step continuity errors : sum local = 1.295907e-09, global = -4.037227e-10, cumulative = -7.344036e-09 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.740195e-08, No Iterations 151 DICPCG: Solving for pcorr, Initial residual = 1.216642e-08, Final residual = 1.216642e-08, No Iterations 0 time step continuity errors : sum local = 1.262291e-16, global = -3.557281e-18, cumulative = -7.344036e-09 MULES: Solving for gamma Liquid phase volume fraction = 0.4999104 Min(gamma) = -6.494791 Max(gamma) = 4.573024 MULES: Solving for gamma Liquid phase volume fraction = 0.4999104 Min(gamma) = -27.25873 Max(gamma) = 25.77633 DICPCG: Solving for pd, Initial residual = 0.1978681, Final residual = 8.198099e-08, No Iterations 77 DICPCG: Solving for pd, Initial residual = 3.041559e-09, Final residual = 3.041559e-09, No Iterations 0 time step continuity errors : sum local = 1.577804e-09, global = 3.218275e-10, cumulative = -7.022209e-09 DICPCG: Solving for pd, Initial residual = 0.09315195, Final residual = 8.495972e-08, No Iterations 83 DICPCG: Solving for pd, Initial residual = 9.250593e-10, Final residual = 9.250593e-10, No Iterations 0 time step continuity errors : sum local = 4.835948e-08, global = 1.095943e-08, cumulative = 3.937219e-09 ExecutionTime = 11.14 s ClockTime = 12 s Courant Number mean: 27.73857 max: 38630.37 Time = 25.01 linear Wave piston position is 1.074151e-05 DICPCG: Solving for cellDisplacementx, Initial residual = 0.4571736, Final residual = 9.227115e-08, No Iterations 380 DICPCG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 1.28 s time step continuity errors : sum local = 4.835948e-08, global = 1.095943e-08, cumulative = 1.489665e-08 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.711118e-08, No Iterations 156 DICPCG: Solving for pcorr, Initial residual = 1.439181e-08, Final residual = 1.439181e-08, No Iterations 0 time step continuity errors : sum local = 4.842107e-15, global = 3.515976e-16, cumulative = 1.489665e-08 MULES: Solving for gamma Liquid phase volume fraction = 0.4999105 Min(gamma) = -49507.96 Max(gamma) = 51611.79 MULES: Solving for gamma ... By the way, my case is a wave tank that generate stable linear wave, so after initial time, the run time information for every time step should looks likes no-change.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post mactech001 CFX 6 November 15, 2009 22:25 siw CFX 1 November 4, 2009 20:42 phsieh2005 OpenFOAM Running, Solving & CFD 2 October 6, 2009 08:33 Leon FLUENT 0 October 3, 2006 21:59 Alberto Schroth FLUENT 0 May 16, 2000 08:19