View Issue Details

IDProjectCategoryView StatusLast Update
0001113OpenFOAMBugpublic2014-01-03 09:38
ReportermartinB Assigned Touser2 
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSOpenSUSEOS Version13.1
Summary0001113: Problems with cellZoneInjection
DescriptionWhen using cellZoneInjection as injection model, some of the parcels do not rebound at the walls. Instead they leave the computational domain. This behaviour does not occur, if cellZoneInjection is replace by coneInjection or another injection method.
Two animated GIF files are attached to visualize the problem.
Steps To ReproduceAttached is a case based on the hopper tutorial. Parcels are first injected with cellZoneInjection and they leave the domain. Next the parcels from first time step (0.005/lagrangian) are copied to another case. These parcels are rebounding correctly.
Run the Allrun script to reproduce the problem.
Additional InformationparaFoam tends to crash when postprocessing the bad case. ParaView works fine.
TagsNo tags attached.

Activities

martinB

2013-12-20 21:37

reporter  

bugReport.tar.gz (3,307 bytes)

martinB

2013-12-20 21:37

reporter  

failure.gif (911,862 bytes)

martinB

2013-12-20 21:38

reporter  

success.gif (954,311 bytes)

martinB

2013-12-21 13:05

reporter   ~0002703

The problem seems to be located at line 101 of CellZoneInjection.C:

If
injectorTetPts.append(cellTetIs[tetI].faceBasePt());

is replaced by
injectorTetPts.append(1);

the injector and the particle rebound work fine.

will

2013-12-31 10:01

manager   ~0002725

Last edited: 2013-12-31 10:02

I've reproduced the problem, and the proposed fix works. However, based on my understanding of the tet decomposition, I think line 101 of CellZoneInjection.C should read:

injectorTetPts.append(cellTetIs[tetI].facePtA());

user2

2014-01-03 09:38

  ~0002739

Thanks for the report - the tetPt() should have been used instead of the faceBasePt(). Resolved by commit 95610f9

Issue History

Date Modified Username Field Change
2013-12-20 21:37 martinB New Issue
2013-12-20 21:37 martinB File Added: bugReport.tar.gz
2013-12-20 21:37 martinB File Added: failure.gif
2013-12-20 21:38 martinB File Added: success.gif
2013-12-21 13:05 martinB Note Added: 0002703
2013-12-31 10:01 will Note Added: 0002725
2013-12-31 10:01 will Note Edited: 0002725
2013-12-31 10:02 will Note Edited: 0002725
2014-01-03 09:38 user2 Note Added: 0002739
2014-01-03 09:38 user2 Status new => resolved
2014-01-03 09:38 user2 Fixed in Version => 2.2.x
2014-01-03 09:38 user2 Resolution open => fixed
2014-01-03 09:38 user2 Assigned To => user2