CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

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

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 11, 2019, 00:55
Smile LPT: Code gets stuck after "Added 1 new parcels"
  #1
Member
 
godfatherBond's Avatar
 
Maximus Arelius
Join Date: Jan 2017
Location: Morocco
Posts: 34
Rep Power: 7
godfatherBond is on a distinguished road
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
__________________
--
🃏Maximus🃏
godfatherBond is offline   Reply With Quote

Old   February 11, 2019, 04:59
Default
  #2
Member
 
ssa
Join Date: Sep 2018
Posts: 93
Rep Power: 6
ssa_cfd is on a distinguished road
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.!
ssa_cfd is offline   Reply With Quote

Old   February 11, 2019, 11:47
Thumbs up
  #3
Member
 
godfatherBond's Avatar
 
Maximus Arelius
Join Date: Jan 2017
Location: Morocco
Posts: 34
Rep Power: 7
godfatherBond is on a distinguished road
Quote:
Originally Posted by ssa_cfd View Post
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
__________________
--
🃏Maximus🃏

Last edited by godfatherBond; February 12, 2019 at 01:07.
godfatherBond is offline   Reply With Quote

Reply

Tags
lpt, openfoam 5.x, particle tracking, spray combustion, spray injection

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 13:30.