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

Incorrect tau in electrostatic pushers

Browse source 
During the improve performant step, an error in the electrostatic
pushers was introduced, resulting in these using the minimum time step
and not the species time step when calculating the acceleration.
This commit is contained in:
Jorge Gonzalez 2023-01-13 15:32:25 +01:00
commit cca57d5ab0
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/modules/solver/pusher/modulePusher.f90
View file

@ -24,7 +24,7 @@ MODULE modulePusher
!Get the electric field at particle position !Get the electric field at particle position
qmEFT = mesh%cells(part%cell)%obj%gatherElectricField(part%Xi) qmEFT = mesh%cells(part%cell)%obj%gatherElectricField(part%Xi)
qmEFt = qmEFt*part%species%qm*tauMin qmEFt = qmEFt*part%species%qm*tauIn
!Update velocity !Update velocity
part%v = part%v + qmEFt part%v = part%v + qmEFt
@ -127,7 +127,7 @@ MODULE modulePusher
part_temp = part part_temp = part
!Get electric field at particle position !Get electric field at particle position
qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi) qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi)
qmEFt = qmEFt*part_temp%species%qm*tauMin qmEFt = qmEFt*part_temp%species%qm*tauIn
!z !z
part_temp%v(1) = part%v(1) + qmEFt(1) part_temp%v(1) = part%v(1) + qmEFt(1)
part_temp%r(1) = part%r(1) + part_temp%v(1)*tauIn part_temp%r(1) = part%r(1) + part_temp%v(1)*tauIn
@ -203,7 +203,7 @@ MODULE modulePusher
part_temp = part part_temp = part
!Get electric field at particle position !Get electric field at particle position
qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi) qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi)
qmEFt = qmEFt*part_temp%species%qm*tauMin qmEFt = qmEFt*part_temp%species%qm*tauIn
!r,theta !r,theta
v_p_oh_star(1) = part%v(1) + qmEFt(1) v_p_oh_star(1) = part%v(1) + qmEFt(1)
x_new = part%r(1) + v_p_oh_star(1)*tauIn x_new = part%r(1) + v_p_oh_star(1)*tauIn