|
[Sponsors] |
LPT: Code gets stuck after "Added 1 new parcels" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 11, 2019, 00:55 |
LPT: Code gets stuck after "Added 1 new parcels"
|
#1 |
Member
Maximus Arelius
Join Date: Jan 2017
Location: Morocco
Posts: 36
Rep Power: 9 |
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 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 ) ); } } } } Please share your views on this problem. Thanks P.S: I am using OpenFOAM 5.x. RS
__________________
-- 🃏Maximus🃏 |
|
February 11, 2019, 04:59 |
|
#2 |
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.!
|
|
February 11, 2019, 11:47 |
|
#3 | |
Member
Maximus Arelius
Join Date: Jan 2017
Location: Morocco
Posts: 36
Rep Power: 9 |
Quote:
If somebody has other suggestions, please share. Any help is highly appericiated! Thanks
__________________
-- 🃏Maximus🃏 Last edited by godfatherBond; February 12, 2019 at 01:07. |
||
Tags |
lpt, openfoam 5.x, particle tracking, spray combustion, spray injection |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |