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

Small change to ionization boundary

Browse source 
Very small change in formats to better account the ionization processes.
Minimum impact in results.
This commit is contained in:
Jorge Gonzalez 2022-12-12 17:36:10 +01:00
commit 2ee58b6b77
1 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/modules/mesh/moduleMeshBoundary.f90
View file

@ -116,7 +116,7 @@ MODULE moduleMeshBoundary
CLASS(meshEdge), INTENT(inout):: edge CLASS(meshEdge), INTENT(inout):: edge
CLASS(particle), INTENT(inout):: part CLASS(particle), INTENT(inout):: part
REAL(8):: vRel, eRel, mRel !relative velocity, energy and mass REAL(8):: vRel, eRel, mRel !relative velocity, energy and mass
REAL(8):: nIonizations !Number of ionizations based on eRel INTEGER:: nIonizations !Number of ionizations based on eRel
REAL(8):: pIonization !Probability of ionization of each event REAL(8):: pIonization !Probability of ionization of each event
INTEGER:: p INTEGER:: p
REAL(8):: v0(1:3) !random velocity of neutral REAL(8):: v0(1:3) !random velocity of neutral
@ -130,11 +130,11 @@ MODULE moduleMeshBoundary
eRel = mRel*vRel**2*5.D-1 eRel = mRel*vRel**2*5.D-1
!Maximum number of possible ionizations based on relative energy !Maximum number of possible ionizations based on relative energy
nIonizations = eRel/bound%eThreshold nIonizations = FLOOR(eRel/bound%eThreshold)
DO p = 1, FLOOR(nIonizations) DO p = 1, nIonizations
!Get probability of ionization !Get probability of ionization
pIonization = 1.D0 - DEXP(-bound%n0*bound%crossSection%get(eRel)*vRel*bound%effectiveTime/nIonizations) pIonization = 1.D0 - DEXP(-bound%n0*bound%crossSection%get(eRel)*vRel*bound%effectiveTime/REAL(nIonizations))
!If a random number is below the probability of ionization, create new pair of ion-electron !If a random number is below the probability of ionization, create new pair of ion-electron
IF (random() < pIonization) THEN IF (random() < pIonization) THEN