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

Compare commits

base: JGonzalez:a9cc3bba50a3f0cdecdd927b5714f27cde2cac0c
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:a1f94c1e4d8ae8d39aef9bec308d3d4f785415a4
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. "a9cc3bba50a3f0cdecdd927b5714f27cde2cac0c" and "a1f94c1e4d8ae8d39aef9bec308d3d4f785415a4" have entirely different histories.
a9cc3bba50 ... a1f94c1e4d

2 changed files with 5 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/modules/mesh/1DRad/moduleMesh1DRad.f90
View file

@ -103,7 +103,6 @@ MODULE moduleMesh1DRad
SUBROUTINE initEdge1DRad(self, n, p, ps) SUBROUTINE initEdge1DRad(self, n, p, ps)
USE moduleSpecies, only:nSpecies USE moduleSpecies, only:nSpecies
USE moduleErrors USE moduleErrors
use moduleConstParam, only: pi4
IMPLICIT NONE IMPLICIT NONE
CLASS(meshEdge1DRad), INTENT(out):: self CLASS(meshEdge1DRad), INTENT(out):: self
@ -122,7 +121,7 @@ MODULE moduleMesh1DRad
self%r = r1(1) self%r = r1(1)
self%surface = pi4 * self%r**2 self%surface = 1.D0
self%normal = (/ 1.D0, 0.D0, 0.D0 /) self%normal = (/ 1.D0, 0.D0, 0.D0 /)

10
src/modules/mesh/moduleMesh@boundaryEM.f90
View file

@ -159,7 +159,7 @@ submodule(moduleMesh) boundaryEM
! Inverse of non-dimensional capacitance ! Inverse of non-dimensional capacitance
self%invC = 1.0d0 / (capacitance / (eps_0 * L_ref)) self%invC = 1.0d0 / (capacitance / (eps_0 * L_ref))
self%charge = 0.0d0 self%charge = self%potential / self%invC
self%update => updateFloating self%update => updateFloating
@ -200,7 +200,7 @@ submodule(moduleMesh) boundaryEM
class(meshEdge), pointer:: edge class(meshEdge), pointer:: edge
real(8), allocatable:: mom_nodes(:) real(8), allocatable:: mom_nodes(:)
class(meshNode), pointer:: node class(meshNode), pointer:: node
real(8):: mom_center, edgeDensityCurrent, chargeTau real(8):: mom_center, edgeDensityCurrent
select type(self) select type(self)
type is(boundaryEMFloating) type is(boundaryEMFloating)
@ -230,14 +230,12 @@ submodule(moduleMesh) boundaryEM
end do end do
chargeTau = edgeDensityCurrent * edge%surface * tauMin
! Accumulate charge of speceis on surface ! Accumulate charge of speceis on surface
self%charge = self%charge + chargeTau self%charge = self%charge + edgeDensityCurrent * edge%surface * tauMin
end do end do
self%potential = self%potential + chargeTau * self%invC self%potential = self%charge * self%invC
end select end select