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

Steady state current doesn't seem to go to 0, but it still not perfect

Browse source
This commit is contained in:
Jorge Gonzalez 2026-03-13 23:00:10 +01:00
commit a1386b0b9c
1 changed files with 6 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/modules/mesh/moduleMesh@boundaryParticle.f90
View file

@ -479,7 +479,7 @@ submodule(moduleMesh) boundaryParticle
type is(boundaryQuasiNeutrality)
write(fileID, '(A,",",A)') '"Edge id"', '"alpha"'
do e = 1, size(self%edges)
write(fileID, '('//fmtColInt//','//fmtReal//')') self%edges(e)%obj%n, self%alpha(self%edges(e)%obj%n)
write(fileID, '('//fmtColInt//','//fmtColReal//')') self%edges(e)%obj%n, self%alpha(self%edges(e)%obj%n)
end do
@ -501,16 +501,12 @@ submodule(moduleMesh) boundaryParticle
v_cut = 2.d0 * 40.d3/v_ref !This will be the drag velocity of the ions in the future
if (dot_product(part%v,-edge%normal) > v_cut) then
! print *,part%v(1), v_cut
! part%v = part%v + v_cut*edge%normal
part%v(1) = part%v(1) - v_cut
! print *,part%v(1)
call genericReflection(edge, part)
! print *,part%v(1)
! print *
call genericReflection(edge, part)
part%v(1) = part%v(1) + v_cut
if (dot_product(part%v,edge%normal) <= 0.d0) then
else
call genericTransparent(edge, part)
end if