godfatherBond |
February 10, 2019 23:55 |
LPT: Code gets stuck after "Added 1 new parcels"
Hello Everyone,
I am trying to run a spray combustion simulation using LES. After running for ~5000 time steps (w/ dt=1e-06) the code gets stuck at :
Solving 3-D cloud fuelCloud1
Cloud: fuelCloud1 injector: model1
Added 1 new parcels
The code gets stuck at this for hours. The mesh quality is Ok as may be seen from this,
Code:
Mesh stats
points: 1218516
faces: 3632893
internal faces: 3611267
cells: 1207360
faces per cell: 6
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 1207360
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
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 Bounding box
inletJet 473 496 ok (non-closed singly connected) (-0.0022941398636 -0.0022941398636 0) (0.0022941398636 0.0022941398636 0)
inletCoflow 3300 3344 ok (non-closed singly connected) (-0.023 -0.023 0) (0.023 0.023 0)
outlet 3773 3796 ok (non-closed singly connected) (-0.046 -0.046 0.276) (0.046 0.046 0.276)
side 14080 14124 ok (non-closed singly connected) (-0.046 -0.046 0) (0.046 0.046 0.276)
Checking geometry...
Overall domain bounding box (-0.046 -0.046 0) (0.046 0.046 0.276)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-3.26238353783e-16 -9.91871009924e-17 1.64145405062e-16) OK.
Max cell openness = 3.11034785118e-16 OK.
Max aspect ratio = 11.3164429291 OK.
Minimum face area = 1.90157343094e-08. Maximum face area = 1.15407802292e-05. Face area magnitudes OK.
Min volume = 1.64010708419e-11. Max volume = 9.93740928183e-09. Total volume = 0.001362704. Cell volumes OK.
Mesh non-orthogonality Max: 42.553608524 average: 13.1870256687
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.48912315273 OK.
Coupled point location match (average 0) OK.
Face tets OK.
Min/max edge length = 0.000143749999992 0.0083636363637 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : min = 0.999999904452 average = 0.999999997621
All face flatness OK.
Cell determinant (wellposedness) : minimum: 0.178877632945 average: 2.13438959315
Cell determinant check OK.
Concave cell check OK.
Face interpolation weight : minimum: 0.24147048594 average: 0.496106306781
Face interpolation weight check OK.
Face volume ratio : minimum: 0.310136821802 average: 0.971367321108
Face volume ratio check OK.
Mesh OK.
End
The injection model used is the patch injection with
Code:
injectionModels
{
model1
{
type patchInjection;
patchName inletJet;
massTotal 0;
parcelBasisType fixed;
SOI 0.001;
duration 50;
parcelsPerSecond 1266683;
nParticle 10;
U0 (0 0 0);
flowRateProfile constant 1;
sizeDistribution
{
type general;
generalDistribution
{
distribution
(
( 3.8109384764E-06 0.0003484941 )
( 4.61043807331E-06 0.00107342111298 )
( 5.40993767023E-06 0.00461269121149 )
( 6.20943726714E-06 0.00981273225724 )
( 7.00893686406E-06 0.0169999950874 )
( 7.80843646097E-06 0.0253469930221 )
( 8.60793605788E-06 0.0355020732935 )
( 9.4074356548E-06 0.0430382816585 )
( 10.2069352517E-06 0.0456821381966 )
( 11.0064348486E-06 0.0487437280437 )
( 11.8059344455E-06 0.0494042002732 )
( 12.6054340425E-06 0.0491658370365 )
( 13.4049336394E-06 0.0502142074269 )
( 14.2044332363E-06 0.0466755561253 )
( 15.0039328332E-06 0.0451245975491 )
( 15.8034324301E-06 0.0425333745328 )
( 16.602932027E-06 0.0389389002971 )
( 17.4024316239E-06 0.0352220209195 )
( 18.2019312209E-06 0.0358097731086 )
( 19.0014308178E-06 0.0347394717635 )
( 19.8009304147E-06 0.0329766447077 )
( 20.6004300116E-06 0.0317908224798 )
( 21.3999296085E-06 0.0306212445313 )
( 22.1994292054E-06 0.0294889789261 )
( 22.9989288023E-06 0.026850356596 )
( 23.7984283993E-06 0.0245094694289 )
( 24.5979279962E-06 0.0229806679815 )
( 25.3974275931E-06 0.0215792478899 )
( 26.19692719E-06 0.0204036316226 )
( 26.9964267869E-06 0.0194293421011 )
( 27.7959263838E-06 0.0180227905273 )
( 28.5954259807E-06 0.0163912595075 )
( 29.3949255777E-06 0.0162337630891 )
( 30.1944251746E-06 0.015486247662 )
( 30.9939247715E-06 0.0150618921753 )
( 31.7934243684E-06 0.0143649152225 )
( 32.5929239653E-06 0.0133282181516 )
( 33.3924235622E-06 0.0127343999081 )
( 34.1919231591E-06 0.0120117516299 )
( 34.9914227561E-06 0.0110242003253 )
( 35.790922353E-06 0.0104723677058 )
( 36.5904219499E-06 0.0101565101053 )
( 37.3899215468E-06 0.0103536006335 )
( 38.1894211437E-06 0.0100271944756 )
( 38.9889207406E-06 0.00929628725347 )
( 39.7884203375E-06 0.00820976837546 )
( 40.5879199345E-06 0.00786692082195 )
( 41.3874195314E-06 0.00690523193374 )
( 42.1869191283E-06 0.00643648626995 )
( 42.9864187252E-06 0.00641398206238 )
( 43.7859183221E-06 0.00639147785481 )
( 44.585417919E-06 0.00558366956976 )
( 45.3849175159E-06 0.00490008428853 )
( 46.1844171128E-06 0.00454244382159 )
( 46.9839167098E-06 0.00427291014442 )
( 47.7834163067E-06 0.00441059652012 )
( 48.5829159036E-06 0.00370303351856 )
( 49.3824155005E-06 0.00372816886233 )
( 50.1819150974E-06 0.00377477657463 )
( 50.9814146943E-06 0.00374524141185 )
( 51.7809142912E-06 0.00338904631198 )
( 52.5804138882E-06 0.00297919419234 )
( 53.3799134851E-06 0.00261792278846 )
( 54.179413082E-06 0.0025452450932 )
( 54.9789126789E-06 0.00231808612267 )
( 55.7784122758E-06 0.00206519300722 )
( 56.5779118727E-06 0.00228040193365 )
( 57.3774114696E-06 0.00237071620751 )
( 58.1769110666E-06 0.00244165584872 )
( 58.9764106635E-06 0.00220079526334 )
( 59.7759102604E-06 0.00185840133667 )
( 60.5754098573E-06 0.0016519385678 )
( 61.3749094542E-06 0.00150762075531 )
( 62.1744090511E-06 0.00141185904696 )
( 62.973908648E-06 0.00158935457369 )
( 63.773408245E-06 0.00165107561334 )
( 64.5729078419E-06 0.00134556552197 )
( 65.3724074388E-06 0.00125110050926 )
( 66.1719070357E-06 0.00107725770879 )
( 66.9714066326E-06 0.000905239157519 )
( 67.7709062295E-06 0.00099533463522 )
( 68.5704058264E-06 0.00108543011292 )
( 69.3699054234E-06 0.000892307701056 )
( 70.1694050203E-06 0.000808784323647 )
( 70.9689046172E-06 0.000734195186446 )
( 71.7684042141E-06 0.000674855843045 )
( 72.567903811E-06 0.00061927141576 )
( 73.3674034079E-06 0.000639932102095 )
( 74.1669030048E-06 0.0006664632 )
( 74.9664026018E-06 0.0006664632 )
( 75.7659021987E-06 0.000616835891609 )
( 76.5654017956E-06 0.00049877617314 )
( 77.3649013925E-06 0.000399043111242 )
( 78.1644009894E-06 0.000385142526712 )
( 78.9639005863E-06 0.000371241942181 )
( 79.7634001832E-06 0.000357341357651 )
( 80.5628997802E-06 0.0003484941 )
( 81.3623993771E-06 0.0003484941 )
( 82.161898974E-06 0.0003484941 )
( 82.9613985709E-06 0.0003484941 )
);
}
}
}
}
The mesh does not have a wall and the particles are set to escape the boundaries.
Please share your views on this problem.
Thanks
P.S: I am using OpenFOAM 5.x.
RS
|