First attempt at Coulomb collisions #46

Merged
JorgeGonz merged 21 commits from feature/CoulombLinear into development 2023-07-16 14:47:59 +02:00
Showing only changes of commit a70788b95d - Show all commits

Fix an issue with compilation

Copy paste error. Now the code compiles.
Jorge Gonzalez 2023-03-08 19:42:15 +01:00

View file

@ -1067,7 +1067,9 @@ MODULE moduleMesh
END DO END DO
!Divide total momentum exchanged among all the particles of species j !Divide total momentum exchanged among all the particles of species j
normDeltaV = totalDeltaV_ij / REAL(cell%listPart_in(j)%amount) * (self%sp_i%weight*self%sp_i%m)/(self%sp_j%weight*self%sp_j%m) normDeltaV = totalDeltaV_ij / REAL(cell%listPart_in(j)%amount) * &
(coulombMatrix(k)%sp_i%weight*coulombMatrix(k)%sp_i%m) / &
(coulombMatrix(k)%sp_j%weight*coulombMatrix(k)%sp_j%m)
!Loop over particles of species_j !Loop over particles of species_j
partTemp => cell%listPart_in(j)%head partTemp => cell%listPart_in(j)%head
DO WHILE(ASSOCIATED(partTemp)) DO WHILE(ASSOCIATED(partTemp))
@ -1124,6 +1126,8 @@ MODULE moduleMesh
END DO END DO
DEALLOCATE(densityNodes, velocityNodes, temperatureNodes)
END DO END DO
!$OMP END DO !$OMP END DO