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

Compare commits

base: JGonzalez:086ce4f589aa05fae4b4866facf0223b5561024f
Branches Tags
JGonzalez:development
JGonzalez:feature/outflowAdaptive
JGonzalez:feature/floatingSurface
JGonzalez:feature/injectionQuasiNeutral
JGonzalez:improvements/mesh_restructure
JGonzalez:feature/see
JGonzalez:feature/CoulombLinear
JGonzalez:improvement/speciesArrays
JGonzalez:master
..
compare: JGonzalez:f4dab6aa71752c65e1c16be24ca6448f864f9ba9
Branches Tags
JGonzalez:feature/outflowAdaptive
JGonzalez:feature/floatingSurface
JGonzalez:feature/injectionQuasiNeutral
JGonzalez:improvements/mesh_restructure
JGonzalez:development
JGonzalez:feature/see
JGonzalez:feature/CoulombLinear
JGonzalez:improvement/speciesArrays
JGonzalez:master

No commits in common. "086ce4f589aa05fae4b4866facf0223b5561024f" and "f4dab6aa71752c65e1c16be24ca6448f864f9ba9" have entirely different histories.
086ce4f589 ... f4dab6aa71

1 changed files with 15 additions and 34 deletions
Download patch file Download diff file Expand all files Collapse all files

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

@ -546,11 +546,10 @@ submodule(moduleMesh) boundaryParticle
integer:: e, n, s integer:: e, n, s
integer, allocatable:: nodes(:) integer, allocatable:: nodes(:)
class(meshEdge), pointer:: edge class(meshEdge), pointer:: edge
real(8), allocatable:: den_nodes(:)
real(8), allocatable:: mom_nodes(:) real(8), allocatable:: mom_nodes(:)
class(meshNode), pointer:: node class(meshNode), pointer:: node
real(8):: den_center, mom_center, vel_center real(8):: mom_center
real(8):: vel_incident, vel_rest real(8):: mom_incident, mom_rest
real(8):: alpha real(8):: alpha
select type(self) select type(self)
@ -558,52 +557,34 @@ submodule(moduleMesh) boundaryParticle
do e = 1, size(self%edges) do e = 1, size(self%edges)
edge => self%edges(e)%obj edge => self%edges(e)%obj
vel_incident = 0.d0 mom_incident = 0.d0
vel_rest = 0.d0 mom_rest = 0.d0
nodes = edge%getNodes(edge%nNodes) nodes = edge%getNodes(edge%nNodes)
allocate(den_nodes(1:edge%nNodes))
allocate(mom_nodes(1:edge%nNodes)) allocate(mom_nodes(1:edge%nNodes))
do s = 1, nSpecies do s = 1, nSpecies
den_center = 0.0d0 do n = 1, edge%nNodes
mom_center = 0.0d0 node => mesh%nodes(nodes(n))%obj
vel_center = 0.0d0
select type(sp => species(s)%obj)
type is(speciesCharged)
do n = 1, edge%nNodes
node => mesh%nodes(nodes(n))%obj
den_nodes(n) = node%output(s)%den mom_nodes(n) = dot_product(node%output(s)%mom, edge%normal)
mom_nodes(n) = dot_product(node%output(s)%mom, edge%normal)
end do end do
den_center = edge%gatherF(edge%centerXi(), edge%nNodes, den_nodes) mom_center = qSpecies(s)*edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
if (den_center > 1.0d-10) then
mom_center = edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
vel_center = mom_center / den_center
else if (s == self%s_incident) then
vel_center = 0.0d0 mom_incident = mom_center
else
end if mom_rest = mom_rest + mom_center
end if
if (s == self%s_incident) then
vel_incident = vel_center
else
vel_rest = vel_rest + vel_center
end if
end select
end do end do
alpha = vel_rest - vel_incident alpha = mom_rest + mom_incident
! Apply correction with a factor of 0.1 to avoid fast changes ! Apply correction with a factor of 0.1 to avoid fast changes
self%velocity_shift(edge%n) = self%velocity_shift(edge%n) + 1.0d-2 * alpha self%velocity_shift(edge%n) = self%velocity_shift(edge%n) + 1.0d-2 * alpha
deallocate(den_nodes)
deallocate(mom_nodes) deallocate(mom_nodes)
end do end do