pimpleDyMFoam computation randomly stops

babapeti · January 27, 2017, 05:12

Dear All!

I have two problems with pimpledymfoam solver if I use adaptive mesh refinement (no moving mesh, just refinement and unrefinement) on OF v3.0+. First some words about the case I would like to simulate in 3D. It is a simple reaction-diffusion case in which there is a constant concentration gradient along the axis x. This gradient governs the diffusion of species A to species B and there is a reaction which converts B to A (this is why there is a constant gradient). Therefore a transient zone slowly propagates along the axis x. I planned that after testing this case, I put some difference in the physical properties of the species which will generate convection in the fluid. But now I am testing the simplest case. I have a working pimplefoam version of the solver on OFv2.2. In v2.2 I run the cases on 48 cpus but the computation is still very long. This is why I wanted to swich on AMR method.

The physics and the mesh geometry are quite simple (a cuboid filled with aqueous solution with no slip boundaries at the walls)

Now the two problems in short (I do not know if the problems are independent from each other or not):

Problem 1:
When I run the case the computation stops without any warning in the log file. CPUs still running on 100% but there are no new results written out anymore and no new notes written in the logfile.

I simply solved this with writing a small script which detects this stop and restarts the computation from the last time step which was saved.

The time when the computation "stops" looks quite random.
This problem only occurs when I run the case in parallel mode.

Problem 2 (and this is more weird than problem no.1.):
When I run the case the computation stops but in a different way. The computation still running, there are new notes in the logfile, the results are written out but the fields (concentration, U and the others I compute) do not evolve with time. Since I have a gradient in c they should change with time as long as I have species B in the domain.

The problem occurs both in parallel and single core mode. What is really interesting is that the number of iterations are really close to zero after this kind of stop (see post below).

If I turn off adaptive mesh refinement the problem does not occur! Unfortunately this is not true. In a recent case (single core, AMR off) the problem also occurs .

If I put some convection in the case, the stop can be delayed with about 100 seconds (in simulation time not in computation).
In single core mode the stop can be delayed with about 50 seconds.
The stop is not a kind of random stop. If I restart the computation from not the last result time step but from a step where I saw reasonable iteration numbers, the computation will crash exactly at the same time step.
The results before the crash are looking good, there aren't any sign of blowing up. (LAST 4 POINTS ARE WITH AMR ON)

At first I blamed it on the MPI but the second problem occurs also in single core mode. Maybe it is a memory problem, but how can it affect on the number of iterations. And if it is the case why it crashes at the same time step after restarting the computation?

I tried paralleltest and it ended up without error.

I am so confused. I tried a lot of scenarios but the problem appeared in all of them. Please if anyone has experiences like that and comments on any of the problems feel free to reply! If you need I can send the solver and the case to try on your machine.

See post #2 for attachments!

Thank You!

babapeti · January 27, 2017, 05:48

Screenshots:

Refinement ON case screenshots

Refinement OFF case screenshots

As you can see with refinement turned on after 8 seconds the evolution of the concentration field simply stops.

(I hope the links are working)

Pieces of the log after the computation stops (problem2 - number of iterations):

Code:

 Starting time loop

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.5

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 8.645913e-05, Final residual = 8.242017e-07, No Iterations 7
time step continuity errors : sum local = 3.510095e-15, global = 5.097001e-24, cumulative = 5.097001e-24
DICPCG:  Solving for p_rgh, Initial residual = 1.126512e-06, Final residual = 1.106687e-08, No Iterations 36
time step continuity errors : sum local = 5.623967e-17, global = 7.389175e-25, cumulative = 5.835918e-24
DICPCG:  Solving for p_rgh, Initial residual = 2.876526e-08, Final residual = 9.917772e-09, No Iterations 9
time step continuity errors : sum local = 4.773698e-17, global = -1.310264e-24, cumulative = 4.525655e-24
DICPCG:  Solving for p_rgh, Initial residual = 1.152389e-08, Final residual = 6.766281e-09, No Iterations 1
time step continuity errors : sum local = 3.09817e-17, global = -4.498018e-24, cumulative = 2.76371e-26
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 2.671926e-05, Final residual = 2.281364e-07, No Iterations 3
time step continuity errors : sum local = 9.695515e-16, global = 5.133446e-24, cumulative = 5.161083e-24
DICPCG:  Solving for p_rgh, Initial residual = 2.312907e-07, Final residual = 8.968321e-09, No Iterations 11
time step continuity errors : sum local = 4.047407e-17, global = -1.977283e-24, cumulative = 3.1838e-24
DICPCG:  Solving for p_rgh, Initial residual = 1.152137e-08, Final residual = 6.442881e-09, No Iterations 1
time step continuity errors : sum local = 2.763167e-17, global = 1.229716e-25, cumulative = 3.306772e-24
DICPCG:  Solving for p_rgh, Initial residual = 6.585413e-09, Final residual = 6.585413e-09, No Iterations 0
time step continuity errors : sum local = 2.793943e-17, global = 1.893967e-24, cumulative = 5.200739e-24
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 0.23 s  ClockTime = 0 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.501

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 7.627866e-05, Final residual = 6.238697e-07, No Iterations 8
time step continuity errors : sum local = 2.658767e-15, global = 5.216914e-25, cumulative = 5.722431e-24
DICPCG:  Solving for p_rgh, Initial residual = 6.55981e-07, Final residual = 9.589431e-09, No Iterations 41
time step continuity errors : sum local = 4.885325e-17, global = -3.220004e-24, cumulative = 2.502427e-24
DICPCG:  Solving for p_rgh, Initial residual = 2.288282e-08, Final residual = 8.055818e-09, No Iterations 9
time step continuity errors : sum local = 3.959288e-17, global = 5.196815e-24, cumulative = 7.699242e-24
DICPCG:  Solving for p_rgh, Initial residual = 9.468538e-09, Final residual = 9.468538e-09, No Iterations 0
time step continuity errors : sum local = 4.01714e-17, global = -1.526458e-24, cumulative = 6.172784e-24
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 2.757707e-08, Final residual = 9.388956e-09, No Iterations 1
time step continuity errors : sum local = 4.224741e-17, global = -2.932825e-24, cumulative = 3.239959e-24
DICPCG:  Solving for p_rgh, Initial residual = 1.007504e-08, Final residual = 7.046053e-09, No Iterations 1
time step continuity errors : sum local = 3.153773e-17, global = 4.653075e-25, cumulative = 3.705267e-24
DICPCG:  Solving for p_rgh, Initial residual = 7.435506e-09, Final residual = 7.435506e-09, No Iterations 0
time step continuity errors : sum local = 3.154602e-17, global = 4.208228e-24, cumulative = 7.913495e-24
DICPCG:  Solving for p_rgh, Initial residual = 7.435372e-09, Final residual = 7.435372e-09, No Iterations 0
time step continuity errors : sum local = 3.154545e-17, global = -7.747431e-25, cumulative = 7.138752e-24
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 0.27 s  ClockTime = 0 s

After some steps...

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.522

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 1.169379e-08, Final residual = 9.02939e-09, No Iterations 1
time step continuity errors : sum local = 3.846098e-17, global = 6.120112e-24, cumulative = 4.769134e-22
DICPCG:  Solving for p_rgh, Initial residual = 1.060744e-08, Final residual = 8.937469e-09, No Iterations 1
time step continuity errors : sum local = 3.809088e-17, global = 9.286871e-25, cumulative = 4.77842e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.971644e-09, Final residual = 8.971644e-09, No Iterations 0
time step continuity errors : sum local = 3.806326e-17, global = 3.323663e-25, cumulative = 4.781744e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.971396e-09, Final residual = 8.971396e-09, No Iterations 0
time step continuity errors : sum local = 3.806221e-17, global = 5.713057e-24, cumulative = 4.838875e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 9.130175e-09, Final residual = 9.130175e-09, No Iterations 0
time step continuity errors : sum local = 3.873585e-17, global = -1.05032e-24, cumulative = 4.828371e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.13285e-09, Final residual = 9.13285e-09, No Iterations 0
time step continuity errors : sum local = 3.874719e-17, global = 8.886698e-25, cumulative = 4.837258e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.133017e-09, Final residual = 9.133017e-09, No Iterations 0
time step continuity errors : sum local = 3.87479e-17, global = -3.168959e-24, cumulative = 4.805569e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.133027e-09, Final residual = 9.133027e-09, No Iterations 0
time step continuity errors : sum local = 3.874795e-17, global = 2.269943e-24, cumulative = 4.828268e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 2.01 s  ClockTime = 2 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.523

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 1.12552e-08, Final residual = 9.283247e-09, No Iterations 1
time step continuity errors : sum local = 3.943464e-17, global = 2.336049e-24, cumulative = 4.851628e-22
DICPCG:  Solving for p_rgh, Initial residual = 1.057525e-08, Final residual = 9.09871e-09, No Iterations 1
time step continuity errors : sum local = 3.869806e-17, global = 1.764931e-24, cumulative = 4.869278e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.120269e-09, Final residual = 9.120269e-09, No Iterations 0
time step continuity errors : sum local = 3.869381e-17, global = 5.424498e-24, cumulative = 4.923523e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.120252e-09, Final residual = 9.120252e-09, No Iterations 0
time step continuity errors : sum local = 3.869375e-17, global = 2.746371e-24, cumulative = 4.950986e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 9.281711e-09, Final residual = 9.281711e-09, No Iterations 0
time step continuity errors : sum local = 3.937875e-17, global = 2.006069e-24, cumulative = 4.971047e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.285808e-09, Final residual = 9.285808e-09, No Iterations 0
time step continuity errors : sum local = 3.939613e-17, global = 4.529846e-24, cumulative = 5.016346e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.286083e-09, Final residual = 9.286083e-09, No Iterations 0
time step continuity errors : sum local = 3.939731e-17, global = 1.809442e-24, cumulative = 5.03444e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.286099e-09, Final residual = 9.286099e-09, No Iterations 0
time step continuity errors : sum local = 3.939737e-17, global = 7.784517e-25, cumulative = 5.042225e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 2.05 s  ClockTime = 2 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.523

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 1.137321e-08, Final residual = 9.859569e-09, No Iterations 1
time step continuity errors : sum local = 4.184218e-17, global = 2.697257e-24, cumulative = 5.069197e-22
DICPCG:  Solving for p_rgh, Initial residual = 1.087583e-08, Final residual = 9.572986e-09, No Iterations 1
time step continuity errors : sum local = 4.067212e-17, global = 8.468252e-24, cumulative = 5.15388e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.586734e-09, Final residual = 9.586734e-09, No Iterations 0
time step continuity errors : sum local = 4.067285e-17, global = 3.660648e-24, cumulative = 5.190486e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.586707e-09, Final residual = 9.586707e-09, No Iterations 0
time step continuity errors : sum local = 4.067273e-17, global = 4.650171e-24, cumulative = 5.236988e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 9.739018e-09, Final residual = 9.739018e-09, No Iterations 0
time step continuity errors : sum local = 4.131893e-17, global = 2.450444e-24, cumulative = 5.261492e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.742975e-09, Final residual = 9.742975e-09, No Iterations 0
time step continuity errors : sum local = 4.133572e-17, global = 6.001968e-25, cumulative = 5.267494e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.743236e-09, Final residual = 9.743236e-09, No Iterations 0
time step continuity errors : sum local = 4.133683e-17, global = 2.285901e-24, cumulative = 5.290353e-22
DICPCG:  Solving for p_rgh, Initial residual = 9.74325e-09, Final residual = 9.74325e-09, No Iterations 0
time step continuity errors : sum local = 4.133688e-17, global = 1.390346e-24, cumulative = 5.304257e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 2.09 s  ClockTime = 2 s

After some steps...

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 1.199312e-08, Final residual = 7.862145e-09, No Iterations 1
time step continuity errors : sum local = 3.32797e-17, global = 6.093217e-24, cumulative = 7.055344e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.678425e-09, Final residual = 8.678425e-09, No Iterations 0
time step continuity errors : sum local = 3.681924e-17, global = 5.016008e-24, cumulative = 7.105504e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.711968e-09, Final residual = 8.711968e-09, No Iterations 0
time step continuity errors : sum local = 3.696156e-17, global = 3.737699e-24, cumulative = 7.142881e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.713971e-09, Final residual = 8.713971e-09, No Iterations 0
time step continuity errors : sum local = 3.697005e-17, global = 6.406955e-26, cumulative = 7.143522e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 8.789569e-09, Final residual = 8.789569e-09, No Iterations 0
time step continuity errors : sum local = 3.729078e-17, global = 1.891779e-24, cumulative = 7.16244e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.793263e-09, Final residual = 8.793263e-09, No Iterations 0
time step continuity errors : sum local = 3.730645e-17, global = 1.688955e-24, cumulative = 7.179329e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.793386e-09, Final residual = 8.793386e-09, No Iterations 0
time step continuity errors : sum local = 3.730698e-17, global = 1.305519e-24, cumulative = 7.192384e-22
DICPCG:  Solving for p_rgh, Initial residual = 8.793389e-09, Final residual = 8.793389e-09, No Iterations 0
time step continuity errors : sum local = 3.730699e-17, global = 1.281519e-24, cumulative = 7.2052e-22
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 2.45 s  ClockTime = 2 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.589

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 8.700032e-09, Final residual = 8.700032e-09, No Iterations 0
time step continuity errors : sum local = 3.691091e-17, global = 8.238318e-25, cumulative = 1.6815e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.771457e-09, Final residual = 8.771457e-09, No Iterations 0
time step continuity errors : sum local = 3.721394e-17, global = 2.235539e-24, cumulative = 1.683736e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.773125e-09, Final residual = 8.773125e-09, No Iterations 0
time step continuity errors : sum local = 3.722102e-17, global = 3.182128e-25, cumulative = 1.684054e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.773187e-09, Final residual = 8.773187e-09, No Iterations 0
time step continuity errors : sum local = 3.722128e-17, global = -2.201766e-24, cumulative = 1.681852e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 8.775236e-09, Final residual = 8.775236e-09, No Iterations 0
time step continuity errors : sum local = 3.722997e-17, global = -3.481434e-24, cumulative = 1.678371e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.775433e-09, Final residual = 8.775433e-09, No Iterations 0
time step continuity errors : sum local = 3.723081e-17, global = -8.705601e-26, cumulative = 1.678284e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.775446e-09, Final residual = 8.775446e-09, No Iterations 0
time step continuity errors : sum local = 3.723087e-17, global = -6.538259e-25, cumulative = 1.67763e-21
DICPCG:  Solving for p_rgh, Initial residual = 8.775447e-09, Final residual = 8.775447e-09, No Iterations 0
time step continuity errors : sum local = 3.723087e-17, global = -6.211939e-25, cumulative = 1.677009e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 7.27 s  ClockTime = 7 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.59

Selected 0 cells for refinement out of 7929.
Selected 0 split points out of a possible 597.
PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 1.258453e-08, Final residual = 4.407897e-09, No Iterations 1
time step continuity errors : sum local = 1.873894e-17, global = 4.938276e-24, cumulative = 1.681947e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.432957e-09, Final residual = 4.432957e-09, No Iterations 0
time step continuity errors : sum local = 1.880734e-17, global = 2.826262e-24, cumulative = 1.684773e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.433999e-09, Final residual = 4.433999e-09, No Iterations 0
time step continuity errors : sum local = 1.881176e-17, global = 4.118492e-25, cumulative = 1.685185e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.434048e-09, Final residual = 4.434048e-09, No Iterations 0
time step continuity errors : sum local = 1.881197e-17, global = -6.389413e-25, cumulative = 1.684546e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 4.435897e-09, Final residual = 4.435897e-09, No Iterations 0
time step continuity errors : sum local = 1.881981e-17, global = 5.886481e-24, cumulative = 1.690433e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.436121e-09, Final residual = 4.436121e-09, No Iterations 0
time step continuity errors : sum local = 1.882077e-17, global = 6.118302e-24, cumulative = 1.696551e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.436141e-09, Final residual = 4.436141e-09, No Iterations 0
time step continuity errors : sum local = 1.882085e-17, global = 2.889028e-24, cumulative = 1.69944e-21
DICPCG:  Solving for p_rgh, Initial residual = 4.436142e-09, Final residual = 4.436142e-09, No Iterations 0
time step continuity errors : sum local = 1.882086e-17, global = 2.228457e-24, cumulative = 1.701669e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 7.32 s  ClockTime = 7 s

Courant Number mean: 3.034126e-06 max: 0.0001558347
Time = 149.59

PIMPLE: iteration 1
DICPCG:  Solving for p_rgh, Initial residual = 7.927407e-09, Final residual = 7.927407e-09, No Iterations 0
time step continuity errors : sum local = 3.363296e-17, global = -5.385338e-24, cumulative = 1.696283e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.980661e-09, Final residual = 7.980661e-09, No Iterations 0
time step continuity errors : sum local = 3.38589e-17, global = 2.044312e-24, cumulative = 1.698328e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.982319e-09, Final residual = 7.982319e-09, No Iterations 0
time step continuity errors : sum local = 3.386593e-17, global = 1.431067e-24, cumulative = 1.699759e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.982383e-09, Final residual = 7.982383e-09, No Iterations 0
time step continuity errors : sum local = 3.38662e-17, global = 2.808914e-24, cumulative = 1.702568e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
PIMPLE: iteration 2
DICPCG:  Solving for p_rgh, Initial residual = 7.98412e-09, Final residual = 7.98412e-09, No Iterations 0
time step continuity errors : sum local = 3.387357e-17, global = 2.427221e-25, cumulative = 1.70281e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.984336e-09, Final residual = 7.984336e-09, No Iterations 0
time step continuity errors : sum local = 3.387449e-17, global = -7.879096e-25, cumulative = 1.702022e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.984352e-09, Final residual = 7.984352e-09, No Iterations 0
time step continuity errors : sum local = 3.387455e-17, global = 5.864928e-24, cumulative = 1.707887e-21
DICPCG:  Solving for p_rgh, Initial residual = 7.984352e-09, Final residual = 7.984352e-09, No Iterations 0
time step continuity errors : sum local = 3.387456e-17, global = 1.76403e-24, cumulative = 1.709651e-21
DILUPBiCG:  Solving for c, Initial residual = 9.999977e-07, Final residual = 9.999977e-07, No Iterations 0
ExecutionTime = 7.35 s  ClockTime = 7 s

If you need anything else which could help to find the solution please feel free to ask for it!

Thanks for the comments and hints!

Best Wishes!
Peter.

mo_na · January 31, 2017, 14:34

Hey Peter,

I am not familiar with pimpleDyMFoam or adaptive mesh refinement. But the only times where a sudden crash without any signs in the logfile happened to me where when I didn't have enough disk space to write the next time step. Probably thats not what happend here, but I thought its worth a try

Cheers,
Mona

babapeti · February 1, 2017, 03:10

Dear Mona,

Thanks for your reply but unfortunately I have plenty of space on the hard drive.

Thanks again,
Peter

javid629 · February 1, 2017, 18:55

hello, i am trying to simulate flow around a propeller, i create two mesh for stationary and rotary region and export them to openfoam,
then i merged them, after mergein i splte the final mesh for determain the rotary domain in dynamic mesh dict, when i try to run the case i saw the below error,
how i can solve this problem
thanks

Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone region0

PIMPLE: no residual control data found. Calculations will employ 3 corrector loops

Reading field p

--> FOAM FATAL ERROR:
Attempt to cast type patch to type lduInterface

From function To& Foam::refCast(From&) [with To = const Foam::lduInterface; From = const Foam::fvPatch]
in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

FOAM aborting

#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) at ??:?
#3 Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPa tchField(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#4 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#5 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam:

imensio nedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) at ??:?
#10 ? at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12 ? at ??:?
Aborted (core dumped)

consuto · January 24, 2018, 06:28

Quote:

Originally Posted by javid629

hello, i am trying to simulate flow around a propeller, i create two mesh for stationary and rotary region and export them to openfoam,
then i merged them, after mergein i splte the final mesh for determain the rotary domain in dynamic mesh dict, when i try to run the case i saw the below error,
how i can solve this problem
thanks

Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone region0

PIMPLE: no residual control data found. Calculations will employ 3 corrector loops

Reading field p

--> FOAM FATAL ERROR:
Attempt to cast type patch to type lduInterface

From function To& Foam::refCast(From&) [with To = const Foam::lduInterface; From = const Foam::fvPatch]
in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

FOAM aborting

#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) at ??:?
#3 Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPa tchField(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#4 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#5 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:

imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam:

imensio nedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) at ??:?
#10 ? at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12 ? at ??:?
Aborted (core dumped)

--------------------------------------------------
So is anyone who is able to find this and solve the problem?

January 27, 2017, 05:12	pimpleDyMFoam computation randomly stops	#1
babapeti New Member Peter Baba Join Date: Jan 2017 Posts: 3 Rep Power: 9	Dear All! I have two problems with pimpledymfoam solver if I use adaptive mesh refinement (no moving mesh, just refinement and unrefinement) on OF v3.0+. First some words about the case I would like to simulate in 3D. It is a simple reaction-diffusion case in which there is a constant concentration gradient along the axis x. This gradient governs the diffusion of species A to species B and there is a reaction which converts B to A (this is why there is a constant gradient). Therefore a transient zone slowly propagates along the axis x. I planned that after testing this case, I put some difference in the physical properties of the species which will generate convection in the fluid. But now I am testing the simplest case. I have a working pimplefoam version of the solver on OFv2.2. In v2.2 I run the cases on 48 cpus but the computation is still very long. This is why I wanted to swich on AMR method. The physics and the mesh geometry are quite simple (a cuboid filled with aqueous solution with no slip boundaries at the walls) Now the two problems in short (I do not know if the problems are independent from each other or not): Problem 1: When I run the case the computation stops without any warning in the log file. CPUs still running on 100% but there are no new results written out anymore and no new notes written in the logfile. I simply solved this with writing a small script which detects this stop and restarts the computation from the last time step which was saved. The time when the computation "stops" looks quite random. This problem only occurs when I run the case in parallel mode. Problem 2 (and this is more weird than problem no.1.): When I run the case the computation stops but in a different way. The computation still running, there are new notes in the logfile, the results are written out but the fields (concentration, U and the others I compute) do not evolve with time. Since I have a gradient in c they should change with time as long as I have species B in the domain. The problem occurs both in parallel and single core mode. What is really interesting is that the number of iterations are really close to zero after this kind of stop (see post below). If I turn off adaptive mesh refinement the problem does not occur! Unfortunately this is not true. In a recent case (single core, AMR off) the problem also occurs . If I put some convection in the case, the stop can be delayed with about 100 seconds (in simulation time not in computation). In single core mode the stop can be delayed with about 50 seconds. The stop is not a kind of random stop. If I restart the computation from not the last result time step but from a step where I saw reasonable iteration numbers, the computation will crash exactly at the same time step. The results before the crash are looking good, there aren't any sign of blowing up. (LAST 4 POINTS ARE WITH AMR ON) At first I blamed it on the MPI but the second problem occurs also in single core mode. Maybe it is a memory problem, but how can it affect on the number of iterations. And if it is the case why it crashes at the same time step after restarting the computation? I tried paralleltest and it ended up without error. I am so confused. I tried a lot of scenarios but the problem appeared in all of them. Please if anyone has experiences like that and comments on any of the problems feel free to reply! If you need I can send the solver and the case to try on your machine. See post #2 for attachments! Thank You! Last edited by babapeti; February 2, 2017 at 08:19.

February 1, 2017, 18:55	pimpleDyMFoam	#5
javid629 New Member Javid Join Date: Dec 2016 Posts: 2 Rep Power: 0	hello, i am trying to simulate flow around a propeller, i create two mesh for stationary and rotary region and export them to openfoam, then i merged them, after mergein i splte the final mesh for determain the rotary domain in dynamic mesh dict, when i try to run the case i saw the below error, how i can solve this problem thanks Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone region0 PIMPLE: no residual control data found. Calculations will employ 3 corrector loops Reading field p --> FOAM FATAL ERROR: Attempt to cast type patch to type lduInterface From function To& Foam::refCast(From&) [with To = const Foam::lduInterface; From = const Foam::fvPatch] in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) at ??:? #3 Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPa tchField(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #4 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #5 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam:imensio nedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:? #8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:? #9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) at ??:? #10 ? at ??:? #11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #12 ? at ??:? Aborted (core dumped) consuto likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Savonius pimpleDyMFoam	vitokad	OpenFOAM Pre-Processing	10	September 16, 2014 10:30
problem in the CFX12.1 parallel computation	BalanceChen	ANSYS	2	July 7, 2011 11:26
Discussion on the computation of the multi-stage axial compressor	BalanceChen	CFX	5	January 5, 2011 03:23
Cell face values computation un unstructured grids	Sergio Rossi	Main CFD Forum	2	May 28, 2006 11:04
Parallel LES computation stops with reason	vvqf	OpenFOAM Running, Solving & CFD	4	March 6, 2006 08:55

January 31, 2017, 14:34		#3
mo_na Member Mona Join Date: Mar 2016 Location: Berlin Posts: 49 Rep Power: 10	Hey Peter, I am not familiar with pimpleDyMFoam or adaptive mesh refinement. But the only times where a sudden crash without any signs in the logfile happened to me where when I didn't have enough disk space to write the next time step. Probably thats not what happend here, but I thought its worth a try Cheers, Mona

February 1, 2017, 03:10		#4
babapeti New Member Peter Baba Join Date: Jan 2017 Posts: 3 Rep Power: 9	Dear Mona, Thanks for your reply but unfortunately I have plenty of space on the hard drive. Thanks again, Peter