View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0000493 | OpenFOAM | Bug | public | 2012-04-03 13:17 | 2013-08-13 16:06 |
| Reporter | Assigned To | ||||
| Priority | high | Severity | major | Reproducibility | always |
| Status | resolved | Resolution | no change required | ||
| Platform | Linux | OS | Ubuntu | OS Version | 10.04 |
| Summary | 0000493: icoUncoupledKinematicParcelFoam spinning particles | ||||
| Description | While trying to simulate the charging of a hopper, a few particles colliding the wall, don't follow their expected path, that would be descending to the bottom, but they stick to the wall and ascend to the top of the hopper itself. The hopper geometry is far different from the tutorial one, indeed it is a round shaped tank, and particles are introduced from above by patchInjection. Can this behaviour be due to the fact that OF works with cartesian coordinates, while this particular geometry should need cylindrical coordinates? In addition, the same "kinematicCloudProperties" parameters have been used to simulate a slice of the hopper, and the strange behaviour reported above doesn't occur. I'm including the entire case if it can be of any help. | ||||
| Tags | No tags attached. | ||||
|
2012-04-03 13:42
|
kinematicCloudProperties (3,767 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object reactingCloud1Properties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solution
{
active true;
coupled false;
transient yes;
cellValueSourceCorrection off;
interpolationSchemes
{
rho cell;
U cellPoint;
mu cell;
}
integrationSchemes
{
U Euler;
}
}
constantProperties
{
parcelTypeId 1;
rhoMin 1e-15;
minParticleMass 1e-15;
rho0 1050;
youngsModulus 5.4e9;
poissonsRatio 0.35;
constantVolume false;
}
subModels
{
particleForces
{
sphereDrag;
gravity;
}
injectionModel patchInjection;
dispersionModel none;
patchInteractionModel standardWallInteraction;
heatTransferModel none;
surfaceFilmModel none;
collisionModel pairCollision;
radiation off;
patchInjectionCoeffs
{
SOI 0;
massTotal 0;
parcelBasisType fixed;
nParticle 1;
patchName inletCentral;
duration 4;
parcelsPerSecond 130;
U0 (-3.53553391 -3.53553391 0);
flowRateProfile constant 35;
sizeDistribution
{
type normal;
normalDistribution
{
minValue 0.037;
maxValue 0.1;
expectation 0.051;
variance 0.01355;
}
}
}
pairCollisionCoeffs
{
// Maximum possible particle diameter expected at any time
maxInteractionDistance 0.1;
writeReferredParticleCloud no;
pairModel pairSpringSliderDashpot;
pairSpringSliderDashpotCoeffs
{
useEquivalentSize no;
alpha 1; //1.4;
b 1.5;
mu 0.58;
cohesionEnergyDensity 0;
collisionResolutionSteps 12;
};
wallModel wallLocalSpringSliderDashpot;
wallLocalSpringSliderDashpotCoeffs
{
useEquivalentSize no;
collisionResolutionSteps 12;
walls
{
youngsModulus 3.4e9;
poissonsRatio 0.25;
alpha 1.4;
b 1.5;
mu 1.5;
cohesionEnergyDensity 0;
}
// frontAndBack
// {
// youngsModulus 200e9;
// poissonsRatio 0.25;
// alpha 1.4;
// b 1.5;
// mu 1.5;
// cohesionEnergyDensity 0;
// }
};
}
standardWallInteractionCoeffs
{
type rebound;
e 0.1;
}
}
cloudFunctions
{}
// ************************************************************************* //
|
|
2012-04-03 13:42
|
blockMeshDict (16,229 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.001;
coordinateSystem cylindrical;
vertices
(
(-2222.8486093203 0 -201.959734561733) //0
(-1936.95973456173 0 487.848609320304) //1
(-1247.1513906797 0 201.959734561733) //2
(-1533.04026543827 0 -487.848609320304) //3
(-3406.36339143178 1683.88593354921 -626.91667603538)//4
(-2518.91667603538 1683.07374924591 1514.36339143178)//5
(-377.636608568224 1681.11406645079 626.91667603538)//6
(-1265.08332396462 1681.92625075409 -1514.36339143178)//7
(-3923.2307501067 2591.16351269417 -753.332759885101)//8
(-2856.8327598851 2585.20749447004 1819.7307501067)//9
(-283.769249893299 2570.83648730583 753.332759885101)//10
(-1350.1672401149 2576.79250552996 -1819.7307501067)//11
(-3923.2307501067 5402.42966774603 -753.332759885101) //12
(-2856.8327598851 5398.36874622957 1819.7307501067) //13
(-283.769249893299 5388.57033225397 753.332759885101) //14
(-1350.1672401149 5392.63125377043 -1819.7307501067) //15
(-2426.74274053272 7059.12226411401 -373.510759090023) //16
(-1898.01075909002 6817.90352603646 902.242740532721) //17
(-622.257259467279 6235.87773588599 373.510759090023) //18
(-1150.98924090998 6477.09647396354 -902.242740532721) //19
(-649.13 -794.4 -137.248) //20 a
(-733.714692949178 -953.277075053695 -320) //21 b
(-952.055126121641 -1096.69230602708 -320) //22 c
(-1220.13 -1090.96 -147.1586) //23 d
(-1220.13 -1090.96 147.1586) //24 e
(-952.055126121641 -1096.69230602708 320) //25 f
(-733.714692949178 -953.277075053695 320) //26 g
(-649.13 -794.4 137.248) //27 h
(-2222.76815754937 0 202.153961954245) //28
(-1532.84603804576 0 487.768157549368) //29
(-1247.23184245063 0 -202.153961954245) //30
(-1937.15396195425 0 -487.768157549368) //31
(-3406.11365612594 1683.8857049934 627.519590398654) //32
(-1264.48040960135 1681.92569897157 1514.11365612594) //33
(-377.886343874061 1681.1142950066 -627.519590398654) //34
(-2519.51959039865 1683.07430102843 -1514.11365612594) //35
(-3922.93065633899 2591.16183661829 754.05725032945) //36
(-1349.44274967055 2576.78845912484 1819.43065633899) //37
(-284.069343661013 2570.83816338171 -754.05725032945) //38
(-2857.55725032945 2585.21154087516 -1819.43065633899) //39
(-3922.93065633899 5402.42852496702 754.05725032945) //40
(-1349.44274967055 5392.62849485785 1819.43065633899) //41
(-284.069343661013 5388.57147503298 -754.05725032945) //42
(-2857.55725032945 5398.37150514215 -1819.43065633899) //43
(-2426.59395070615 7059.05438304106 373.869969508513) //44
(-1150.63003049149 6476.93259455602 902.093950706152) //45
(-622.406049293849 6235.94561695894 -373.869969508513) //46
(-1898.36996950851 6818.06740544398 -902.093950706152) //47
(-1755.39548825747 6752.83941630179 95.6708579876307) //48
(-1620.09646307434 6691.11313292332 230.969883170751) //49
(-1428.75474709908 6603.81898600373 230.969883170751) //50
(-1293.45572191596 6542.09270262526 95.6708579876307) //51
(-1293.45572191596 6542.09270262526 -95.6708579876307)//52
(-1428.75474709908 6603.81898600373 -230.969883170751) //53
(-1620.09646307434 6691.11313292332 -230.969883170751) //54
(-1755.39548825747 6752.83941630179 -95.6708579876307) //55
(-2934.84153253085 5398.66580933942 344.41508875547) //56
(-2447.76504187162 5396.81098645039 831.491579414706) //57
(-1758.93486436067 5394.1878707706 831.491579414706) //58
(-1271.85837370144 5392.33304788157 344.41508875547) //59
(-1271.85837370144 5392.33304788157 -344.41508875547)//60
(-1758.93486436067 5394.1878707706 -831.491579414706)//61
(-2447.76504187162 5396.81098645039 -831.491579414706)//62
(-2934.84153253085 5398.66580933942 -344.41508875547)//63
(-2934.84153253085 2585.64318703114 344.41508875547) //64
(-2447.76504187162 2582.92278012723 831.491579414706) //65
(-1758.93486436067 2579.07554379689 831.491579414706) //66
(-1271.85837370144 2576.35513689298 344.41508875547) //67
(-1271.85837370144 2576.35513689298 -344.41508875547)//68
(-1758.93486436067 2579.07554379689 -831.491579414706)//69
(-2447.76504187162 2582.92278012723 -831.491579414706)//70
(-2934.84153253085 2585.64318703114 -344.41508875547)//71
(-2584.78478185934 1683.13403122195 287.012573962892)//72<--
(-2178.88770630997 1682.76255738832 692.909649512255)//73
(-1604.86255838419 1682.23721405587 692.909649512255)//74
(-1198.96548283483 1681.86574022224 287.012573962892)//75
(-1198.96548283483 1681.86574022224 -287.012573962892)//76
(-1604.86255838419 1682.23721405587 -692.909649512255)//77
(-2178.88770630997 1682.76255738832 -692.909649512255)//78
(-2584.78478185934 1683.13403122195 -287.012573962892)//79
(-1938.2132713048 0 84.190355029115)//80
(-1819.15012914365 0 203.253497190261)//81
(-1650.76941908542 0 203.253497190261)//82
(-1531.70627692427 0 84.190355029115)//83
(-1531.70627692427 0 -84.190355029115)//84
(-1650.76941908542 0 -203.253497190261)//85
(-1819.15012914365 0 -203.253497190261)//86
(-1938.2132713048 0 -84.190355029115)//87
(-952.055126121641 -1096.69230602708 38.2683431950523)//88
(-952.055126121641 -1096.69230602708 92.3879532683006)//89
(-733.714692949178 -953.277075053695 92.3879532683006)//90
(-733.714692949178 -953.277075053695 38.2683431950523)// 91
(-733.714692949178 -953.277075053695 -38.2683431950523)//92
(-733.714692949178 -953.277075053695 -92.3879532683006)//93
(-952.055126121641 -1096.69230602708 -92.3879532683006)//94
(-952.055126121641 -1096.69230602708 -38.2683431950523)//95
);
blocks
(
hex (0 28 80 87 4 32 72 79) (10 10 10) simpleGrading (1 1 1)
hex (28 1 81 80 32 5 73 72) (10 10 10) simpleGrading (1 1 1)
hex (1 29 82 81 5 33 74 73) (10 10 10) simpleGrading (1 1 1)
hex (29 2 83 82 33 6 75 74) (10 10 10) simpleGrading (1 1 1)
hex (2 30 84 83 6 34 76 75) (10 10 10) simpleGrading (1 1 1)
hex (30 3 85 84 34 7 77 76) (10 10 10) simpleGrading (1 1 1)
hex (3 31 86 85 7 35 78 77) (10 10 10) simpleGrading (1 1 1)
hex (31 0 87 86 35 4 79 78) (10 10 10) simpleGrading (1 1 1)
hex (80 81 82 83 72 73 74 75) (10 10 10) simpleGrading (1 1 1)
hex (80 83 84 87 72 75 76 79) (10 10 10) simpleGrading (1 1 1)
hex (87 84 85 86 79 76 77 78) (10 10 10) simpleGrading (1 1 1)
hex (4 32 72 79 8 36 64 71) (10 10 10) simpleGrading (1 1 1)
hex (32 5 73 72 36 9 65 64) (10 10 10) simpleGrading (1 1 1)
hex (5 33 74 73 9 37 66 65) (10 10 10) simpleGrading (1 1 1)
hex (33 6 75 74 37 10 67 66) (10 10 10) simpleGrading (1 1 1)
hex (6 34 76 75 10 38 68 67) (10 10 10) simpleGrading (1 1 1)
hex (34 7 77 76 38 11 69 68) (10 10 10) simpleGrading (1 1 1)
hex (7 35 78 77 11 39 70 69) (10 10 10) simpleGrading (1 1 1)
hex (35 4 79 78 39 8 71 70) (10 10 10) simpleGrading (1 1 1)
hex (72 73 74 75 64 65 66 67) (10 10 10) simpleGrading (1 1 1)
hex (72 75 76 79 64 67 68 71) (10 10 10) simpleGrading (1 1 1)
hex (79 76 77 78 71 68 69 70) (10 10 10) simpleGrading (1 1 1)
hex (64 65 66 67 56 57 58 59) (10 10 10) simpleGrading (1 1 1)
hex (64 67 68 71 56 59 60 63) (10 10 10) simpleGrading (1 1 1)
hex (71 68 69 70 63 60 61 62) (10 10 10) simpleGrading (1 1 1)
hex (8 36 64 71 12 40 56 63) (10 10 10) simpleGrading (1 1 1)
hex (36 9 65 64 40 13 57 56) (10 10 10) simpleGrading (1 1 1)
hex (9 37 66 65 13 41 58 57) (10 10 10) simpleGrading (1 1 1)
hex (37 10 67 66 41 14 59 58) (10 10 10) simpleGrading (1 1 1)
hex (10 38 68 67 14 42 60 59) (10 10 10) simpleGrading (1 1 1)
hex (38 11 69 68 42 15 61 60) (10 10 10) simpleGrading (1 1 1)
hex (11 39 70 69 15 43 62 61) (10 10 10) simpleGrading (1 1 1)
hex (39 8 71 70 43 12 63 62) (10 10 10) simpleGrading (1 1 1)
hex (56 57 58 59 48 49 50 51) (10 10 10) simpleGrading (1 1 1)
hex (56 59 60 63 48 51 52 55) (10 10 10) simpleGrading (1 1 1)
hex (63 60 61 62 55 52 53 54) (10 10 10) simpleGrading (1 1 1)
hex (12 40 56 63 16 44 48 55) (10 10 10) simpleGrading (1 1 1)
hex (40 13 57 56 44 17 49 48) (10 10 10) simpleGrading (1 1 1)
hex (13 41 58 57 17 45 50 49) (10 10 10) simpleGrading (1 1 1)
hex (41 14 59 58 45 18 51 50) (10 10 10) simpleGrading (1 1 1)
hex (14 42 60 59 18 46 52 51) (10 10 10) simpleGrading (1 1 1)
hex (42 15 61 60 46 19 53 52) (10 10 10) simpleGrading (1 1 1)
hex (15 43 62 61 19 47 54 53) (10 10 10) simpleGrading (1 1 1)
hex (43 12 63 62 47 16 55 54) (10 10 10) simpleGrading (1 1 1)
hex (88 89 90 91 80 81 82 83) (10 10 10) simpleGrading (1 1 1)//<---
hex (95 88 91 92 87 80 83 84) (10 10 10) simpleGrading (1 1 1)
hex (95 92 93 94 87 84 85 86) (10 10 10) simpleGrading (1 1 1)
hex (23 24 88 95 0 28 80 87) (10 10 10) simpleGrading (1 1 1)
hex (24 25 89 88 28 1 81 80) (10 10 10) simpleGrading (1 1 1)
hex (25 26 90 89 1 29 82 81) (10 10 10) simpleGrading (1 1 1)
hex (26 27 91 90 29 2 83 82) (10 10 10) simpleGrading (1 1 1)
hex (27 20 92 91 2 30 84 83) (10 10 10) simpleGrading (1 1 1)
hex (20 21 93 92 30 3 85 84) (10 10 10) simpleGrading (1 1 1)
hex (21 22 94 93 3 31 86 85) (10 10 10) simpleGrading (1 1 1)
hex (22 23 95 94 31 0 87 86) (10 10 10) simpleGrading (1 1 1)//<---
//hex () (10 10 10) simpleGrading (1 1 1)
);
//mergePatchPairs
//(
//);
edges
(
arc 0 28 (-2222.808393166 0 202.056852289)
arc 28 1 (-1937.056852289 0 487.808393166)
arc 1 29 (-1532.943147711 0 487.808393166)
arc 29 2 (-1247.191606834 0 202.056852289)
arc 2 30 (-1247.191606834 0 -202.056852289)
arc 30 3 (-1532.943147711 0 -487.808393166)
arc 3 31 (-1937.056852289 0 -487.808393166)
arc 31 0 (-2222.808393166 0 -202.056852289)
arc 4 32 (-3406.2385537852 1683.88581929877 627.218145648314)
arc 32 5 (-2519.21814564831 1683.07402514855 1514.2385537852)
arc 5 33 (-1264.78185435169 1681.92597485145 1514.2385537852)
arc 33 6 (-377.761446214801 1681.11418070123 627.218145648314)
arc 6 34 (-377.761446214801 1681.11418070123 -627.218145648314)
arc 34 7 (-1264.78185435169 1681.92597485145 -1514.2385537852)
arc 7 35 (-2519.21814564831 1683.07402514855 -1514.2385537852)
arc 35 4 (-3406.2385537852 1683.88581929877 -627.218145648314)
arc 8 36 (-3923.08073928002 2591.16267485762 753.695020045366)
arc 36 9 (-2857.19502004537 2585.20951775603 1819.58073928002)
arc 9 37 (-1349.80497995463 2576.79048224397 1819.58073928002)
arc 37 10 (-283.919260719983 2570.83732514238 753.695020045366)
arc 10 38 (-283.919260719983 2570.83732514238 -753.695020045366)
arc 38 11 (-1349.80497995463 2576.79048224397 -1819.58073928002)
arc 11 39 (-2857.19502004537 2585.20951775603 -1819.58073928002)
arc 39 8 (-3923.08073928002 2591.16267485762 -753.695020045366)
arc 12 40 (-3923.08073928002 5402.42909649383 753.695020045366)
arc 40 13 (-2857.19502004537 5398.37012574275 1819.58073928002)
arc 13 41 (-1349.80497995463 5392.62987425725 1819.58073928002)
arc 41 14 (-283.919260719983 5388.57090350617 753.695020045366)
arc 14 42 (-283.919260719983 5388.57090350617 -753.695020045366)
arc 42 15 (-1349.80497995463 5392.62987425725 -1819.58073928002)
arc 15 43 (-2857.19502004537 5398.37012574275 -1819.58073928002)
arc 43 12 (-3923.08073928002 5402.42909649383 -753.695020045366)
arc 16 44 (-2426.66836347811 7059.08833172504 373.69037175077)
arc 44 17 (-1898.19037170566 6817.98546911917 902.168363496795)
arc 17 45 (-1150.80962829434 6477.01453088083 902.168363496795)
arc 45 18 (-622.331636503205 6235.91166826644 373.690371705661)
arc 18 46 (-622.331636503205 6235.91166826644 -373.690371705661)
arc 46 19 (-1150.80962829434 6477.01453088083 -902.168363496795)
arc 19 47 (-1898.19037170566 6817.98546911917 -902.168363496795)
arc 47 16 (-2426.66836347811 7059.08833172504 -373.69037175077)
//arc 25 26 (-696.21 -900.74 320)
//arc 22 21 (-696.21 -900.74 -320)
//arc 89 90 (-696.21 -900.74 92.3879532683006)
//arc 94 93 (-696.21 -900.74 -92.3879532683006)
//arc 88 91 (-696.21 -900.74 38.2683431950523)
//arc 95 92 (-696.21 -900.74 -38.2683431950523)
//arc 27 91 (-696.21 -900.74 83.5832500830023)
//arc 20 92 (-696.21 -900.74 -83.5832500830023)
//arc 24 88 (-696.21 -900.74 -83.5832500830023)
//arc 23 95 (-696.21 -900.74 83.5832500830023)
//arc 90 26 (-696.21 -696.21 138.065147194396)
);
boundary
(
walls
{
type wall;
faces
(
// (0 28 80 87) //<-------|e
// (28 1 81 80) // |l
// (1 29 82 81) // |i
// (29 2 83 82) // |m
// (2 30 84 83) // |i
// (30 3 85 84) // |n
// (3 31 86 85) // |a
// (31 0 87 86) // |r
// (80 81 82 83)// |e
// (80 83 84 87)// |dopo inserimento ottagonale
// (87 84 85 86)//<-------
(88 89 90 91)
(95 88 91 92)
(95 92 93 94)
(23 24 88 95)
(24 25 89 88)
(25 26 90 89)
(26 27 91 90)
(27 20 92 91)
(20 21 93 92)
(21 22 94 93)
(22 23 95 94)
(23 24 28 0)
(24 25 1 28)
(25 26 29 1)
(26 27 2 29)
(27 20 30 2)
(20 21 3 30)
(21 22 31 3)
(22 23 0 31)
(0 28 32 4)
(28 1 5 32)
(1 29 33 5)
(29 2 6 33)
(2 30 34 6)
(30 3 7 34)
(3 31 35 7)
(31 0 4 35)
(4 32 36 8)
(32 5 9 36)
(5 33 37 9)
(33 6 10 37)
(6 34 38 10)
(34 7 11 38)
(7 35 39 11)
(35 4 8 39)
(8 36 40 12)
(36 9 13 40)
(9 37 41 13)
(37 10 14 41)
(10 38 42 14)
(38 11 15 42)
(11 39 43 15)
(39 8 12 43)
(12 40 44 16)
(40 13 17 44)
(13 41 45 17)
(41 14 18 45)
(14 42 46 18)
(42 15 19 46)
(15 43 47 19)
(43 12 16 47)
(16 44 48 55)
(44 17 49 48)
(17 45 50 49)
(45 18 51 50)
(18 46 52 51)
(46 19 53 52)
(19 47 54 53)
(47 16 55 54)
);
}
inletCentral
{
type patch;
faces
(
(48 49 50 51)
(48 51 52 55)
(55 52 53 54)
);
}
dummy
{
type empty;
faces
(
//
);
}
// outlet
//{
//type patch;
//faces
//(
//(24 25 26 27)
//(24 27 20 23)
//(23 20 21 22)
//);
//}
);
// ************************************************************************* //
|
|
2012-04-03 13:43
|
U (1,500 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type fixedValue;
value uniform (0 0 0);
}
inletCentral
{
type zeroGradient;
value uniform (0 0 0);
}
dummy
{
type empty;
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
connect1
{
type empty;
}
connect2
{
type empty;
}
}
// ************************************************************************* //
|
|
|
What you are observing here is a by-product of employing the patchInjection model to introduce your particles - this injection model is not necessarily suitable for DEM calculations since it allows particles to be injected within the bounds of existing particles - this leads to excessive repulsive forces causing the particles to 'shoot' apart at high/non-physical speeds, and follow unexpected trajectories. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2012-04-03 13:17 |
|
New Issue | |
| 2012-04-03 13:42 |
|
File Added: kinematicCloudProperties | |
| 2012-04-03 13:42 |
|
File Added: blockMeshDict | |
| 2012-04-03 13:43 |
|
File Added: U | |
| 2013-08-13 16:06 |
|
Note Added: 0002395 | |
| 2013-08-13 16:06 |
|
Status | new => resolved |
| 2013-08-13 16:06 |
|
Resolution | open => no change required |
| 2013-08-13 16:06 |
|
Assigned To | => user2 |
