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

Do not put initial charge, just change the voltaje over time

Browse source
This commit is contained in:
Jorge Gonzalez 2026-04-10 22:23:24 +02:00
commit bec3a1045e
1 changed files with 6 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

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 = self%potential / self%invC self%charge = 0.0d0
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 real(8):: mom_center, edgeDensityCurrent, chargeTau
select type(self) select type(self)
type is(boundaryEMFloating) type is(boundaryEMFloating)
@ -230,12 +230,14 @@ 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 + edgeDensityCurrent * edge%surface * tauMin self%charge = self%charge + chargeTau
end do end do
self%potential = self%charge * self%invC self%potential = self%potential + chargeTau * self%invC
end select end select