JGonzalez/fpakc
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Issue with reflecting particles solved.

Browse source
This commit is contained in:
Jorge Gonzalez 2020-10-23 13:41:44 +02:00
commit a04b46a0f0
2 changed files with 12 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

10
src/modules/moduleMeshCyl.f95
View file

@ -484,6 +484,16 @@ MODULE moduleMeshCyl
CLASS IS (meshEdgeCyl) CLASS IS (meshEdgeCyl)
!Particle encountered an edge, execute boundary !Particle encountered an edge, execute boundary
CALL nextElement%fBoundary(part) CALL nextElement%fBoundary(part)
!If particle is still inside the domain, call findCell
IF (part%n_in) THEN
IF(PRESENT(oldCell)) THEN
CALL self%findCell(part, oldCell)
ELSE
CALL self%findCell(part)
END IF
END IF
CLASS DEFAULT CLASS DEFAULT
WRITE(*,*) "ERROR, CHECK findCellCylQuad" WRITE(*,*) "ERROR, CHECK findCellCylQuad"

2
src/modules/moduleMeshCylBoundary.f95
View file

@ -49,6 +49,8 @@ MODULE moduleMeshCylBoundary
part%v(1) = cosT*vpp(1) + sinT*vpp(2) part%v(1) = cosT*vpp(1) + sinT*vpp(2)
part%v(2) = -sinT*vpp(1) + cosT*vpp(2) part%v(2) = -sinT*vpp(1) + cosT*vpp(2)
part%n_in = .TRUE.
END SUBROUTINE reflection END SUBROUTINE reflection
!Absoption in a surface !Absoption in a surface