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

Merge branch 'issue/electricFieldGather' into 'development'

Browse source 
Incorrect tau in electrostatic pushers

See merge request JorgeGonz/fpakc!34
This commit is contained in:
Jorge Gonzalez 2023-01-13 14:38:09 +00:00
commit d36b2db342
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