LPT: Code gets stuck after "Added 1 new parcels"

godfatherBond · February 11, 2019, 00:55

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

ssa_cfd · February 11, 2019, 04:59

I had this problem in manual injection. I think it is unable to handle when there are more particles. Try to reduce the number of parcelspersecond. In my case it worked.!

godfatherBond · February 11, 2019, 11:47

Quote:

Originally Posted by ssa_cfd

I had this problem in manual injection. I think it is unable to handle when there are more particles. Try to reduce the number of parcelspersecond. In my case it worked.!

Thanks for replying! I will try it asap but I am using ~1.27 million parcels per second now as you may see in the file above, is it too large?

If somebody has other suggestions, please share.

Any help is highly appericiated!
Thanks

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
own Code vs. commercial code	Bernhard Mueck	Main CFD Forum	10	February 16, 2000 11:07
State of the art in CFD technology	Juan Carlos GARCIA SALAS	Main CFD Forum	39	November 1, 1999 15:34
Use of an hydrodynamic model integrated in the vertical by the PHOENICS code	Hafsia Zouhaier	Main CFD Forum	2	March 15, 1999 09:23
public CFD Code development	Heinz Wilkening	Main CFD Forum	38	March 5, 1999 12:44
What kind of Cmmercial CFD code you feel well?	Lans	Main CFD Forum	13	October 27, 1998 11:20

February 11, 2019, 04:59		#2
ssa_cfd Member ssa Join Date: Sep 2018 Posts: 93 Rep Power: 8	I had this problem in manual injection. I think it is unable to handle when there are more particles. Try to reduce the number of parcelspersecond. In my case it worked.!