View Issue Details
ID  Project  Category  View Status  Date Submitted  Last Update 

0003234  OpenFOAM  Bug  public  20190509 21:04  20190521 09:42 
Reporter  rmn_990  Assigned To  will  
Priority  low  Severity  minor  Reproducibility  have not tried 
Status  closed  Resolution  unable to reproduce  
Product Version  6  
Fixed in Version  
Summary  0003234: MPPICFoam/DPMFoam Wrong number of parcels on stick boundary  
Description  By using MPPICFoam, I injected 550,000 parcels (nParticle=1) inside a pipe and I considered sticking boundary conditions for the pipe wall. What happens is that the number of parcels which are deposited on the pipe wall is higher than what is injected. Flow is laminar, incompressible with gravity in parallel computation. Attached you can find checkMesh report, KinematicCloudProperties file, and part of the log file which can show the number of stuck particles. I tried doing this simulation with a lower number of parcels but I faced this problem again.  
Tags  No tags attached.  

checkMesh (3,280 bytes)
/**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: v1812   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ Build : v1812 OPENFOAM=1812 Arch : "LSB;label=32;scalar=64" Exec : checkMesh Date : May 09 2019 Time : 15:57:16 Host : raminatomslab PID : 27520 I/O : uncollated Case : /home/ramin/Desktop/temp2_animation nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring runtime modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Time = 0 Mesh stats points: 316025 internal points: 290311 faces: 2489345 internal faces: 2439265 cells: 1148898 faces per cell: 4.28986 boundary patches: 5 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 0 prisms: 333018 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 815880 polyhedra: 0 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology INLET 1258 1014 ok (nonclosed singly connected) WALL_IN 6584 3388 ok (nonclosed singly connected) LUMEN_WALL 34444 17318 ok (nonclosed singly connected) OUTLET 1248 1009 ok (nonclosed singly connected) WALL_OUT 6546 3369 ok (nonclosed singly connected) Checking faceZone topology for multiply connected surfaces... No faceZones found. Checking basic cellZone addressing... No cellZones found. Checking geometry... Overall domain bounding box (0.00635 0.00635 0) (0.00635 0.00635 0.08001) Mesh has 3 geometric (nonempty/wedge) directions (1 1 1) Mesh has 3 solution (nonempty) directions (1 1 1) Boundary openness (5.05755e18 1.44969e15 7.95046e16) OK. Max cell openness = 5.22612e16 OK. Max aspect ratio = 12.4196 OK. Minimum face area = 1.87741e09. Maximum face area = 1.76629e07. Face area magnitudes OK. Min volume = 1.23317e13. Max volume = 2.39655e11. Total volume = 9.16435e06. Cell volumes OK. Mesh nonorthogonality Max: 62.513 average: 14.6154 Nonorthogonality check OK. Face pyramids OK. Max skewness = 2.49864 OK. Coupled point location match (average 0) OK. Mesh OK. End
kinematicCloudProperties (3,682 bytes)
/** C++ **\ =========  \\ / F ield  OpenFOAM: The Open Source CFD Toolbox \\ / O peration  Website: https://openfoam.org \\ / A nd  Version: 6 \\/ M anipulation  \**/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object particleProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solution { active true; coupled true; transient yes; cellValueSourceCorrection off; maxCo 1.0; interpolationSchemes { rho cell; U cellPoint; mu cell; } averagingMethod dual; integrationSchemes { U Euler; } sourceTerms { schemes { U semiImplicit 0.00078; } } } constantProperties { rho0 1050; alphaMax 0.9; } subModels { particleForces { ErgunWenYuDrag { alphac alphac; } gravity; } injectionModels { model1 { type patchInjection; parcelBasisType fixed; SOI 0; patch INLET; duration 5; nParticle 1; parcelsPerSecond 110000; massTotal 2.84e11; U0 (0 0 0); flowRateProfile constant 2.5e8; sizeDistribution { type fixedValue; fixedValueDistribution { value 0.000001; } } } } dispersionModel none; patchInteractionModel localInteraction; localInteractionCoeffs { patches ( WALL_IN { type rebound; e 0.97; mu 0.09; } WALL_OUT { type rebound; e 0.97; mu 0.09; } LUMEN_WALL { type stick; } OUTLET { type escape; } INLET { type rebound; e 0.97; mu 0.09; } ); } heatTransferModel none; surfaceFilmModel none; packingModel implicit; explicitCoeffs { particleStressModel { type HarrisCrighton; alphaPacked 0.6; pSolid 10.0; beta 2.0; eps 1.0e7; } correctionLimitingMethod { type absolute; e 0.9; } } implicitCoeffs { alphaMin 0.0001; rhoMin 1.0; applyLimiting true; applyGravity false; particleStressModel { type HarrisCrighton; alphaPacked 0.6; pSolid 5.0; beta 2.0; eps 1.0e2; } } dampingModel none; isotropyModel stochastic; stochasticCoeffs { timeScaleModel { type isotropic; alphaPacked 0.6; e 0.9; } } stochasticCollisionModel none; radiation off; } cloudFunctions {} // ************************************************************************* //
log file (18,237 bytes)
Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.74047e05, Final residual = 4.93955e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.33195e16 1.07411e16 1.65086e15) Linear momentum = 1.65971e15 Linear kinetic energy = 8.10248e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2264940, 1.18592e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 26.3872 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000773465, Final residual = 2.90678e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742042, Final residual = 2.83409e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655442, Final residual = 2.45896e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000363683, Final residual = 1.23966e06, No Iterations 2 time step continuity errors : sum local = 3.08911e08, global = 1.97904e11, cumulative = 2.03826e07 GAMGPCG: Solving for p, Initial residual = 0.00015826, Final residual = 6.95268e07, No Iterations 2 time step continuity errors : sum local = 1.73254e08, global = 2.88311e11, cumulative = 2.03855e07 ExecutionTime = 155431 s ClockTime = 156154 s Courant Number mean: 0.18969 max: 10.7764 Time = 45.353 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107993, Final residual = 4.95123e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.33643e16 1.06987e16 1.65083e15) Linear momentum = 1.65968e15 Linear kinetic energy = 8.10281e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2265432, 1.18618e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000773749, Final residual = 2.90797e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742236, Final residual = 2.83527e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655512, Final residual = 2.45934e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000363957, Final residual = 1.27501e06, No Iterations 2 time step continuity errors : sum local = 3.17721e08, global = 1.9626e11, cumulative = 2.03835e07 GAMGPCG: Solving for p, Initial residual = 0.00015847, Final residual = 7.15167e07, No Iterations 2 time step continuity errors : sum local = 1.78213e08, global = 2.8688e11, cumulative = 2.03864e07 ExecutionTime = 155440 s ClockTime = 156163 s Courant Number mean: 0.189689 max: 10.7761 Time = 45.354 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.73532e05, Final residual = 4.9377e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.34123e16 1.06896e16 1.65081e15) Linear momentum = 1.6597e15 Linear kinetic energy = 8.103e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2265924, 1.18644e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774036, Final residual = 2.90918e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742425, Final residual = 2.8364e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655578, Final residual = 2.45973e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364179, Final residual = 1.32287e06, No Iterations 2 time step continuity errors : sum local = 3.29648e08, global = 1.92535e11, cumulative = 2.03844e07 GAMGPCG: Solving for p, Initial residual = 0.00015865, Final residual = 7.39101e07, No Iterations 2 time step continuity errors : sum local = 1.84178e08, global = 2.84279e11, cumulative = 2.03873e07 ExecutionTime = 155448 s ClockTime = 156171 s Courant Number mean: 0.189687 max: 10.7759 Time = 45.355 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000108002, Final residual = 4.95222e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.34462e16 1.07276e16 1.65075e15) Linear momentum = 1.65969e15 Linear kinetic energy = 8.10319e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2266416, 1.1867e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774323, Final residual = 2.91041e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742607, Final residual = 2.83745e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655639, Final residual = 2.46011e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364374, Final residual = 1.36658e06, No Iterations 2 time step continuity errors : sum local = 3.40541e08, global = 1.88494e11, cumulative = 2.03854e07 GAMGPCG: Solving for p, Initial residual = 0.000158807, Final residual = 7.60462e07, No Iterations 2 time step continuity errors : sum local = 1.89501e08, global = 2.82329e11, cumulative = 2.03882e07 ExecutionTime = 155457 s ClockTime = 156180 s Courant Number mean: 0.189686 max: 10.7757 Time = 45.356 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.73883e05, Final residual = 4.93987e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.34128e16 1.08115e16 1.65053e15) Linear momentum = 1.65949e15 Linear kinetic energy = 8.10337e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2266911, 1.18695e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.00077461, Final residual = 2.9117e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.00074278, Final residual = 2.83845e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655697, Final residual = 2.46053e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364577, Final residual = 1.40622e06, No Iterations 2 time step continuity errors : sum local = 3.50419e08, global = 1.84088e11, cumulative = 2.03864e07 GAMGPCG: Solving for p, Initial residual = 0.000158996, Final residual = 7.80148e07, No Iterations 2 time step continuity errors : sum local = 1.94407e08, global = 2.82342e11, cumulative = 2.03892e07 ExecutionTime = 155465 s ClockTime = 156189 s Courant Number mean: 0.189684 max: 10.7754 Time = 45.357 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107947, Final residual = 4.95018e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.33038e16 1.08835e16 1.65063e15) Linear momentum = 1.65955e15 Linear kinetic energy = 8.10373e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2267405, 1.18721e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774901, Final residual = 2.91308e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742951, Final residual = 2.8394e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.00065576, Final residual = 2.461e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364802, Final residual = 1.44512e06, No Iterations 2 time step continuity errors : sum local = 3.60116e08, global = 1.79768e11, cumulative = 2.03874e07 GAMGPCG: Solving for p, Initial residual = 0.000159202, Final residual = 7.96791e07, No Iterations 2 time step continuity errors : sum local = 1.98555e08, global = 2.8587e11, cumulative = 2.03903e07 ExecutionTime = 155474 s ClockTime = 156197 s Courant Number mean: 0.189683 max: 10.7752 Time = 45.358 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.74196e05, Final residual = 4.93861e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.32042e16 1.09311e16 1.65066e15) Linear momentum = 1.65954e15 Linear kinetic energy = 8.10417e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2267899, 1.18747e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000775202, Final residual = 2.91455e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000743115, Final residual = 2.84027e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655825, Final residual = 2.46151e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000365067, Final residual = 1.46651e06, No Iterations 2 time step continuity errors : sum local = 3.65447e08, global = 1.77551e11, cumulative = 2.03885e07 GAMGPCG: Solving for p, Initial residual = 0.000159396, Final residual = 8.04721e07, No Iterations 2 time step continuity errors : sum local = 2.00532e08, global = 2.93408e11, cumulative = 2.03914e07 ExecutionTime = 155483 s ClockTime = 156207 s Courant Number mean: 0.189681 max: 10.7749 Time = 45.359 Evolving kinematicCloud Solving 3D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107982, Final residual = 4.95282e05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e10 Linear momentum = (1.32957e16 1.09242e16 1.65081e15) Linear momentum = 1.65975e15 Linear kinetic energy = 8.1046e18 model1: number of parcels added = 550000 mass introduced = 2.87979e10 Parcel fate (number, mass) : patch WALL_IN  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch WALL_OUT  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch LUMEN_WALL  escape = 0, 0  stick = 2268393, 1.18773e09 Parcel fate (number, mass) : patch OUTLET  escape = 0, 0  stick = 0, 0 Parcel fate (number, mass) : patch INLET  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 9.96681e05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000775507, Final residual = 2.91605e05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.00074327, Final residual = 2.84107e05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655884, Final residual = 2.46198e05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000365264, Final residual = 1.48518e06, No Iterations 2 time step continuity errors : sum local = 3.701e08, global = 1.77174e11, cumulative = 2.03896e07 GAMGPCG: Solving for p, Initial residual = 0.000159529, Final residual = 8.09602e07, No Iterations 2 time step continuity errors : sum local = 2.01749e08, global = 3.02367e11, cumulative = 2.03927e07 ExecutionTime = 155492 s ClockTime = 156216 s Courant Number mean: 0.18968 max: 10.7747 Time = 45.36 

I also tried icoUncoupledKinematicParcelFoam and it provides a correct number of parcels on the stick boundary condition. In my case just MPPICFoam should be used and I shouldn't use an uncoupled solver. 

Can you upload a case that reproduces the issue, please? 

Thanks. Attached you can find it. 

https://1drv.ms/u/s!AkqWyjWvlskVlhXALhLGEEo5G4en 

That case took you 43 hours to run on 40 cores. That is too big to debug. Have you got anything smaller which demonstrates the problem? 

You are right. I have uploaded another case with log file and some time step outputs. https://1drv.ms/u/s!AkqWyjWvlskVlhdQ_T10HMAZCyG5 sorry if it is still not very fast. Because this error can be seen after some time. 

I did this simulation with different particle size and number but I faced this bug again. icoUncoupledKinematicParcelFoam works without this bug and the number of particles on surfaces with stick bc is true. 

If you turn the packing model off in constant/kinematicCloudProperties, does that fix it? I.e., replace "packingModel implicit;" with "packingModel none;". 

No, it will not be fixed because I did this simulation with packingModel implicit in OP5 and everything worked well. 

Will, I tried what you told about packingModel and it doesn't work on OP6. 

log_PackingOFF.txt (17,524 bytes)
Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 237.993 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388221 Current mass in system = 2.13436e10 Linear momentum = (5.81069e16 8.67692e16 3.05555e14) Linear momentum = 3.05733e14 Linear kinetic energy = 6.09623e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161783, 8.89448e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161084034, 8.85609e08 Parcel fate: patch OUTLET (number, mass)  escape = 161783, 8.8945e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 2.16902e05 Min dense number of parcels = 5.10842 Continous phase1 volume fraction = 1 Min(alphac) = 0.999978 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000673368, Final residual = 0.000305183, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000641985, Final residual = 0.000304481, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.06681e06, Final residual = 1.41077e07, No Iterations 3 time step continuity errors : sum local = 8.07579e11, global = 8.512e12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.5909e06, Final residual = 2.24164e07, No Iterations 5 time step continuity errors : sum local = 1.2832e10, global = 1.3728e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.94458e06, Final residual = 9.61243e09, No Iterations 448 time step continuity errors : sum local = 5.50252e12, global = 1.35174e14, cumulative = 1 ExecutionTime = 46640.1 s ClockTime = 46849 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388198 Current mass in system = 2.13423e10 Linear momentum = (5.79615e16 8.63523e16 3.05466e14) Linear momentum = 3.05643e14 Linear kinetic energy = 6.09596e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161806, 8.89575e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161130283, 8.85863e08 Parcel fate: patch OUTLET (number, mass)  escape = 161806, 8.89576e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 1.80752e05 Min dense number of parcels = 4.27346 Continous phase1 volume fraction = 1 Min(alphac) = 0.999982 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000694296, Final residual = 0.000306, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000660099, Final residual = 0.000306549, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.02897e06, Final residual = 1.44372e07, No Iterations 3 time step continuity errors : sum local = 8.26447e11, global = 9.14134e12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.6074e06, Final residual = 2.16544e07, No Iterations 6 time step continuity errors : sum local = 1.23959e10, global = 1.47832e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97092e06, Final residual = 9.45472e09, No Iterations 445 time step continuity errors : sum local = 5.41227e12, global = 2.9251e14, cumulative = 1 ExecutionTime = 46641.6 s ClockTime = 46851 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238.007 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388188 Current mass in system = 2.13417e10 Linear momentum = (5.87336e16 8.58988e16 3.05461e14) Linear momentum = 3.05638e14 Linear kinetic energy = 6.09576e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161816, 8.8963e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161176536, 8.86117e08 Parcel fate: patch OUTLET (number, mass)  escape = 161816, 8.89631e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 2.00488e05 Min dense number of parcels = 4.74008 Continous phase1 volume fraction = 1 Min(alphac) = 0.99998 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000693061, Final residual = 0.000306107, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.00065908, Final residual = 0.000306787, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03002e06, Final residual = 1.50115e07, No Iterations 3 time step continuity errors : sum local = 8.59314e11, global = 9.76491e12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.57637e06, Final residual = 2.19112e07, No Iterations 6 time step continuity errors : sum local = 1.25428e10, global = 1.56968e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.9644e06, Final residual = 9.96242e09, No Iterations 444 time step continuity errors : sum local = 5.70286e12, global = 2.69179e14, cumulative = 1 ExecutionTime = 46643.4 s ClockTime = 46852 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238.014 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388173 Current mass in system = 2.13409e10 Linear momentum = (5.97394e16 8.50179e16 3.05426e14) Linear momentum = 3.05602e14 Linear kinetic energy = 6.09535e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161831, 8.89712e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161222806, 8.86372e08 Parcel fate: patch OUTLET (number, mass)  escape = 161831, 8.89713e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 1.78598e05 Min dense number of parcels = 4.22254 Continous phase1 volume fraction = 1 Min(alphac) = 0.999982 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.00069305, Final residual = 0.000305972, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000659156, Final residual = 0.000306633, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.00649e06, Final residual = 1.29918e07, No Iterations 4 time step continuity errors : sum local = 7.43704e11, global = 1.04235e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.61046e06, Final residual = 2.18902e07, No Iterations 6 time step continuity errors : sum local = 1.25309e10, global = 1.66328e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97581e06, Final residual = 9.74588e09, No Iterations 446 time step continuity errors : sum local = 5.57895e12, global = 1.40598e14, cumulative = 1 ExecutionTime = 46644.6 s ClockTime = 46854 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238.021 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388161 Current mass in system = 2.13403e10 Linear momentum = (5.95773e16 8.48949e16 3.05405e14) Linear momentum = 3.05581e14 Linear kinetic energy = 6.09507e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161843, 8.89778e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161269081, 8.86626e08 Parcel fate: patch OUTLET (number, mass)  escape = 161843, 8.89779e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 1.65841e05 Min dense number of parcels = 3.92093 Continous phase1 volume fraction = 1 Min(alphac) = 0.999983 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000692808, Final residual = 0.000305801, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000658822, Final residual = 0.000306701, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03043e06, Final residual = 1.30005e07, No Iterations 4 time step continuity errors : sum local = 7.44199e11, global = 1.10879e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.61072e06, Final residual = 2.23184e07, No Iterations 6 time step continuity errors : sum local = 1.27758e10, global = 1.77773e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.96983e06, Final residual = 9.30023e09, No Iterations 448 time step continuity errors : sum local = 5.3238e12, global = 4.33324e15, cumulative = 1 ExecutionTime = 46645.9 s ClockTime = 46855 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238.028 Evolving kinematicCloud Solving 3D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388144 Current mass in system = 2.13393e10 Linear momentum = (5.97055e16 8.45076e16 3.0535e14) Linear momentum = 3.05526e14 Linear kinetic energy = 6.09494e16 Injector model1:  parcels added = 550004  mass introduced = 3.0238e10 Parcel fate: system (number, mass)  escape = 161860, 8.89872e11 Parcel fate: patch WALL_IN (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch WALL_OUT (number, mass)  escape = 0, 0  stick = 0, 0 Parcel fate: patch LUMEN_WALL (number, mass)  escape = 0, 0  stick = 161315354, 8.8688e08 Parcel fate: patch OUTLET (number, mass)  escape = 161860, 8.89873e11  stick = 0, 0 Parcel fate: patch INLET (number, mass)  escape = 0, 0  stick = 0, 0 Min cell volume fraction = 0 Max cell volume fraction = 1.75191e05 Min dense number of parcels = 4.12606 Continous phase1 volume fraction = 1 Min(alphac) = 0.999982 Max(alphac) = 1 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000693007, Final residual = 0.000305759, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.000659283, Final residual = 0.000306587, No Iterations 100 Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03613e06, Final residual = 1.31899e07, No Iterations 4 time step continuity errors : sum local = 7.55043e11, global = 1.18886e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.62919e06, Final residual = 2.30573e07, No Iterations 7 time step continuity errors : sum local = 1.3199e10, global = 1.898e11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97176e06, Final residual = 9.3065e09, No Iterations 449 time step continuity errors : sum local = 5.32743e12, global = 7.14281e15, cumulative = 1 ExecutionTime = 46647.1 s ClockTime = 46856 s Courant Number mean: 0.115384 max: 1.0014 Interface Courant Number mean: 0 max: 0 deltaT = 0.00689655 Time = 238.034 

I cannot see why the reporting should be incorrect, and I cannot reproduce the problem. I have run your case from a time of 30 on 16 cores and over the first few second of physical time the number of stuck parcels reported matches the saved data. I don't doubt that there is a problem, I just can't reproduce or efficiently investigate it with the information you have given me. If you want this investigated further you will need to supply a *smaller* case which demonstrates the problem, so that it can be run and tested and modified in a much shorter time frame. If you can produce such as case, then please open a new report. In the meantime you will have to rely on the actual written out "active" field to determine how many parcels have stuck to the wall. 
Date Modified  Username  Field  Change 

20190509 21:04  rmn_990  New Issue  
20190509 21:04  rmn_990  File Added: checkMesh  
20190509 21:04  rmn_990  File Added: kinematicCloudProperties  
20190509 21:04  rmn_990  File Added: log file  
20190509 21:21  rmn_990  Note Added: 0010456  
20190510 08:54  will  Note Added: 0010457  
20190510 15:11  rmn_990  Note Added: 0010459  
20190510 15:12  rmn_990  Note Added: 0010460  
20190513 08:59  will  Note Added: 0010462  
20190513 18:27  rmn_990  Note Added: 0010464  
20190515 04:26  rmn_990  Note Added: 0010465  
20190516 15:16  will  Note Added: 0010468  
20190516 15:16  will  Note Edited: 0010468  View Revisions 
20190516 17:29  rmn_990  Note Added: 0010476  
20190518 18:33  rmn_990  Note Added: 0010479  
20190518 18:36  rmn_990  File Added: log_PackingOFF.txt  
20190521 09:41  will  Priority  high => low 
20190521 09:41  will  Severity  major => minor 
20190521 09:41  will  Description Updated  View Revisions 
20190521 09:42  will  Assigned To  => will 
20190521 09:42  will  Status  new => closed 
20190521 09:42  will  Resolution  open => unable to reproduce 
20190521 09:42  will  Note Added: 0010480 