CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

calculation diverge after continue to run

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2010, 04:35
Default calculation diverge after continue to run
  #1
Member
 
jingjing
Join Date: Mar 2009
Location: shanghai,China
Posts: 30
Rep Power: 7
zhajingjing is on a distinguished road
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.



zhajingjing is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX doesn't continue calculation... mactech001 CFX 6 November 15, 2009 21:25
What's the best order to run this simulation in? siw CFX 1 November 4, 2009 19:42
serial run fine, but parallel run diverged phsieh2005 OpenFOAM Running, Solving & CFD 2 October 6, 2009 08:33
How to run and save tranisent case on Cray XD1 Leon FLUENT 0 October 3, 2006 21:59
Heat Flux Calculation under REPORTS Alberto Schroth FLUENT 0 May 16, 2000 08:19


All times are GMT -4. The time now is 22:57.