View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0003234 | OpenFOAM | Bug | public | 2019-05-09 21:04 | 2019-05-21 09:42 |
Reporter | rmn_990 | Assigned To | will | ||
Priority | low | Severity | minor | Reproducibility | have not tried |
Status | closed | Resolution | unable to reproduce | ||
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 : ramin-atomslab PID : 27520 I/O : uncollated Case : /home/ramin/Desktop/temp2_animation nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied 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 (non-closed singly connected) WALL_IN 6584 3388 ok (non-closed singly connected) LUMEN_WALL 34444 17318 ok (non-closed singly connected) OUTLET 1248 1009 ok (non-closed singly connected) WALL_OUT 6546 3369 ok (non-closed 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 (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (-5.05755e-18 -1.44969e-15 -7.95046e-16) OK. Max cell openness = 5.22612e-16 OK. Max aspect ratio = 12.4196 OK. Minimum face area = 1.87741e-09. Maximum face area = 1.76629e-07. Face area magnitudes OK. Min volume = 1.23317e-13. Max volume = 2.39655e-11. Total volume = 9.16435e-06. Cell volumes OK. Mesh non-orthogonality Max: 62.513 average: 14.6154 Non-orthogonality 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.84e-11; U0 (0 0 0); flowRateProfile constant 2.5e-8; 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.0e-7; } 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.0e-2; } } 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 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.74047e-05, Final residual = 4.93955e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.33195e-16 -1.07411e-16 1.65086e-15) |Linear momentum| = 1.65971e-15 Linear kinetic energy = 8.10248e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18592e-09 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.96681e-05 Min dense number of parcels = 26.3872 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000773465, Final residual = 2.90678e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742042, Final residual = 2.83409e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655442, Final residual = 2.45896e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000363683, Final residual = 1.23966e-06, No Iterations 2 time step continuity errors : sum local = 3.08911e-08, global = 1.97904e-11, cumulative = -2.03826e-07 GAMGPCG: Solving for p, Initial residual = 0.00015826, Final residual = 6.95268e-07, No Iterations 2 time step continuity errors : sum local = 1.73254e-08, global = -2.88311e-11, cumulative = -2.03855e-07 ExecutionTime = 155431 s ClockTime = 156154 s Courant Number mean: 0.18969 max: 10.7764 Time = 45.353 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107993, Final residual = 4.95123e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.33643e-16 -1.06987e-16 1.65083e-15) |Linear momentum| = 1.65968e-15 Linear kinetic energy = 8.10281e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18618e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000773749, Final residual = 2.90797e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742236, Final residual = 2.83527e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655512, Final residual = 2.45934e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000363957, Final residual = 1.27501e-06, No Iterations 2 time step continuity errors : sum local = 3.17721e-08, global = 1.9626e-11, cumulative = -2.03835e-07 GAMGPCG: Solving for p, Initial residual = 0.00015847, Final residual = 7.15167e-07, No Iterations 2 time step continuity errors : sum local = 1.78213e-08, global = -2.8688e-11, cumulative = -2.03864e-07 ExecutionTime = 155440 s ClockTime = 156163 s Courant Number mean: 0.189689 max: 10.7761 Time = 45.354 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.73532e-05, Final residual = 4.9377e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.34123e-16 -1.06896e-16 1.65081e-15) |Linear momentum| = 1.6597e-15 Linear kinetic energy = 8.103e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18644e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774036, Final residual = 2.90918e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742425, Final residual = 2.8364e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655578, Final residual = 2.45973e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364179, Final residual = 1.32287e-06, No Iterations 2 time step continuity errors : sum local = 3.29648e-08, global = 1.92535e-11, cumulative = -2.03844e-07 GAMGPCG: Solving for p, Initial residual = 0.00015865, Final residual = 7.39101e-07, No Iterations 2 time step continuity errors : sum local = 1.84178e-08, global = -2.84279e-11, cumulative = -2.03873e-07 ExecutionTime = 155448 s ClockTime = 156171 s Courant Number mean: 0.189687 max: 10.7759 Time = 45.355 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000108002, Final residual = 4.95222e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.34462e-16 -1.07276e-16 1.65075e-15) |Linear momentum| = 1.65969e-15 Linear kinetic energy = 8.10319e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.1867e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774323, Final residual = 2.91041e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742607, Final residual = 2.83745e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655639, Final residual = 2.46011e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364374, Final residual = 1.36658e-06, No Iterations 2 time step continuity errors : sum local = 3.40541e-08, global = 1.88494e-11, cumulative = -2.03854e-07 GAMGPCG: Solving for p, Initial residual = 0.000158807, Final residual = 7.60462e-07, No Iterations 2 time step continuity errors : sum local = 1.89501e-08, global = -2.82329e-11, cumulative = -2.03882e-07 ExecutionTime = 155457 s ClockTime = 156180 s Courant Number mean: 0.189686 max: 10.7757 Time = 45.356 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.73883e-05, Final residual = 4.93987e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.34128e-16 -1.08115e-16 1.65053e-15) |Linear momentum| = 1.65949e-15 Linear kinetic energy = 8.10337e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18695e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.00077461, Final residual = 2.9117e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.00074278, Final residual = 2.83845e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655697, Final residual = 2.46053e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364577, Final residual = 1.40622e-06, No Iterations 2 time step continuity errors : sum local = 3.50419e-08, global = 1.84088e-11, cumulative = -2.03864e-07 GAMGPCG: Solving for p, Initial residual = 0.000158996, Final residual = 7.80148e-07, No Iterations 2 time step continuity errors : sum local = 1.94407e-08, global = -2.82342e-11, cumulative = -2.03892e-07 ExecutionTime = 155465 s ClockTime = 156189 s Courant Number mean: 0.189684 max: 10.7754 Time = 45.357 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107947, Final residual = 4.95018e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.33038e-16 -1.08835e-16 1.65063e-15) |Linear momentum| = 1.65955e-15 Linear kinetic energy = 8.10373e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18721e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000774901, Final residual = 2.91308e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000742951, Final residual = 2.8394e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.00065576, Final residual = 2.461e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000364802, Final residual = 1.44512e-06, No Iterations 2 time step continuity errors : sum local = 3.60116e-08, global = 1.79768e-11, cumulative = -2.03874e-07 GAMGPCG: Solving for p, Initial residual = 0.000159202, Final residual = 7.96791e-07, No Iterations 2 time step continuity errors : sum local = 1.98555e-08, global = -2.8587e-11, cumulative = -2.03903e-07 ExecutionTime = 155474 s ClockTime = 156197 s Courant Number mean: 0.189683 max: 10.7752 Time = 45.358 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 9.74196e-05, Final residual = 4.93861e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.32042e-16 -1.09311e-16 1.65066e-15) |Linear momentum| = 1.65954e-15 Linear kinetic energy = 8.10417e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18747e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000775202, Final residual = 2.91455e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.000743115, Final residual = 2.84027e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655825, Final residual = 2.46151e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000365067, Final residual = 1.46651e-06, No Iterations 2 time step continuity errors : sum local = 3.65447e-08, global = 1.77551e-11, cumulative = -2.03885e-07 GAMGPCG: Solving for p, Initial residual = 0.000159396, Final residual = 8.04721e-07, No Iterations 2 time step continuity errors : sum local = 2.00532e-08, global = -2.93408e-11, cumulative = -2.03914e-07 ExecutionTime = 155483 s ClockTime = 156207 s Courant Number mean: 0.189681 max: 10.7749 Time = 45.359 Evolving kinematicCloud Solving 3-D cloud kinematicCloud GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.000107982, Final residual = 4.95282e-05, No Iterations 1000 Cloud: kinematicCloud Current number of parcels = 550000 Current mass in system = 2.87979e-10 Linear momentum = (1.32957e-16 -1.09242e-16 1.65081e-15) |Linear momentum| = 1.65975e-15 Linear kinetic energy = 8.1046e-18 model1: number of parcels added = 550000 mass introduced = 2.87979e-10 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.18773e-09 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.96681e-05 Min dense number of parcels = 28.6195 DILUPBiCGStab: Solving for U.airx, Initial residual = 0.000775507, Final residual = 2.91605e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airy, Initial residual = 0.00074327, Final residual = 2.84107e-05, No Iterations 1 DILUPBiCGStab: Solving for U.airz, Initial residual = 0.000655884, Final residual = 2.46198e-05, No Iterations 1 GAMGPCG: Solving for p, Initial residual = 0.000365264, Final residual = 1.48518e-06, No Iterations 2 time step continuity errors : sum local = 3.701e-08, global = 1.77174e-11, cumulative = -2.03896e-07 GAMGPCG: Solving for p, Initial residual = 0.000159529, Final residual = 8.09602e-07, No Iterations 2 time step continuity errors : sum local = 2.01749e-08, global = -3.02367e-11, cumulative = -2.03927e-07 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388221 Current mass in system = 2.13436e-10 Linear momentum = (5.81069e-16 -8.67692e-16 3.05555e-14) |Linear momentum| = 3.05733e-14 Linear kinetic energy = 6.09623e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161783, 8.89448e-11 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.85609e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161783, 8.8945e-11 - 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.16902e-05 Min dense number of parcels = 5.10842 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.06681e-06, Final residual = 1.41077e-07, No Iterations 3 time step continuity errors : sum local = 8.07579e-11, global = -8.512e-12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.5909e-06, Final residual = 2.24164e-07, No Iterations 5 time step continuity errors : sum local = 1.2832e-10, global = -1.3728e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.94458e-06, Final residual = 9.61243e-09, No Iterations 448 time step continuity errors : sum local = 5.50252e-12, global = -1.35174e-14, 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388198 Current mass in system = 2.13423e-10 Linear momentum = (5.79615e-16 -8.63523e-16 3.05466e-14) |Linear momentum| = 3.05643e-14 Linear kinetic energy = 6.09596e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161806, 8.89575e-11 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.85863e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161806, 8.89576e-11 - 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.80752e-05 Min dense number of parcels = 4.27346 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.02897e-06, Final residual = 1.44372e-07, No Iterations 3 time step continuity errors : sum local = 8.26447e-11, global = 9.14134e-12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.6074e-06, Final residual = 2.16544e-07, No Iterations 6 time step continuity errors : sum local = 1.23959e-10, global = 1.47832e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97092e-06, Final residual = 9.45472e-09, No Iterations 445 time step continuity errors : sum local = 5.41227e-12, global = 2.9251e-14, 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388188 Current mass in system = 2.13417e-10 Linear momentum = (5.87336e-16 -8.58988e-16 3.05461e-14) |Linear momentum| = 3.05638e-14 Linear kinetic energy = 6.09576e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161816, 8.8963e-11 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.86117e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161816, 8.89631e-11 - 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.00488e-05 Min dense number of parcels = 4.74008 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03002e-06, Final residual = 1.50115e-07, No Iterations 3 time step continuity errors : sum local = 8.59314e-11, global = -9.76491e-12, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.57637e-06, Final residual = 2.19112e-07, No Iterations 6 time step continuity errors : sum local = 1.25428e-10, global = -1.56968e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.9644e-06, Final residual = 9.96242e-09, No Iterations 444 time step continuity errors : sum local = 5.70286e-12, global = -2.69179e-14, 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388173 Current mass in system = 2.13409e-10 Linear momentum = (5.97394e-16 -8.50179e-16 3.05426e-14) |Linear momentum| = 3.05602e-14 Linear kinetic energy = 6.09535e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161831, 8.89712e-11 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.86372e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161831, 8.89713e-11 - 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.78598e-05 Min dense number of parcels = 4.22254 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.00649e-06, Final residual = 1.29918e-07, No Iterations 4 time step continuity errors : sum local = 7.43704e-11, global = 1.04235e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.61046e-06, Final residual = 2.18902e-07, No Iterations 6 time step continuity errors : sum local = 1.25309e-10, global = 1.66328e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97581e-06, Final residual = 9.74588e-09, No Iterations 446 time step continuity errors : sum local = 5.57895e-12, global = 1.40598e-14, 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388161 Current mass in system = 2.13403e-10 Linear momentum = (5.95773e-16 -8.48949e-16 3.05405e-14) |Linear momentum| = 3.05581e-14 Linear kinetic energy = 6.09507e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161843, 8.89778e-11 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.86626e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161843, 8.89779e-11 - 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.65841e-05 Min dense number of parcels = 3.92093 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03043e-06, Final residual = 1.30005e-07, No Iterations 4 time step continuity errors : sum local = 7.44199e-11, global = -1.10879e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.61072e-06, Final residual = 2.23184e-07, No Iterations 6 time step continuity errors : sum local = 1.27758e-10, global = -1.77773e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.96983e-06, Final residual = 9.30023e-09, No Iterations 448 time step continuity errors : sum local = 5.3238e-12, global = 4.33324e-15, 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 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 388144 Current mass in system = 2.13393e-10 Linear momentum = (5.97055e-16 -8.45076e-16 3.0535e-14) |Linear momentum| = 3.05526e-14 Linear kinetic energy = 6.09494e-16 Injector model1: - parcels added = 550004 - mass introduced = 3.0238e-10 Parcel fate: system (number, mass) - escape = 161860, 8.89872e-11 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.8688e-08 Parcel fate: patch OUTLET (number, mass) - escape = 161860, 8.89873e-11 - 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.75191e-05 Min dense number of parcels = 4.12606 Continous phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 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 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 3.03613e-06, Final residual = 1.31899e-07, No Iterations 4 time step continuity errors : sum local = 7.55043e-11, global = 1.18886e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 4.62919e-06, Final residual = 2.30573e-07, No Iterations 7 time step continuity errors : sum local = 1.3199e-10, global = 1.898e-11, cumulative = 1 DICPCG: Solving for p_rgh, Initial residual = 1.97176e-06, Final residual = 9.3065e-09, No Iterations 449 time step continuity errors : sum local = 5.32743e-12, global = -7.14281e-15, 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 |
---|---|---|---|
2019-05-09 21:04 | rmn_990 | New Issue | |
2019-05-09 21:04 | rmn_990 | File Added: checkMesh | |
2019-05-09 21:04 | rmn_990 | File Added: kinematicCloudProperties | |
2019-05-09 21:04 | rmn_990 | File Added: log file | |
2019-05-09 21:21 | rmn_990 | Note Added: 0010456 | |
2019-05-10 08:54 | will | Note Added: 0010457 | |
2019-05-10 15:11 | rmn_990 | Note Added: 0010459 | |
2019-05-10 15:12 | rmn_990 | Note Added: 0010460 | |
2019-05-13 08:59 | will | Note Added: 0010462 | |
2019-05-13 18:27 | rmn_990 | Note Added: 0010464 | |
2019-05-15 04:26 | rmn_990 | Note Added: 0010465 | |
2019-05-16 15:16 | will | Note Added: 0010468 | |
2019-05-16 15:16 | will | Note Edited: 0010468 | |
2019-05-16 17:29 | rmn_990 | Note Added: 0010476 | |
2019-05-18 18:33 | rmn_990 | Note Added: 0010479 | |
2019-05-18 18:36 | rmn_990 | File Added: log_PackingOFF.txt | |
2019-05-21 09:41 | will | Priority | high => low |
2019-05-21 09:41 | will | Severity | major => minor |
2019-05-21 09:41 | will | Description Updated | |
2019-05-21 09:42 | will | Assigned To | => will |
2019-05-21 09:42 | will | Status | new => closed |
2019-05-21 09:42 | will | Resolution | open => unable to reproduce |
2019-05-21 09:42 | will | Note Added: 0010480 |