Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (master) Browse

fennecs/src/materials_mod.f90
dcbebf0c577cmaster

Last Change
Raw File

materials_mod.f90
View Options

!------------------------------------------------------------------------------
! EPFL/Swiss Plasma Cente
!------------------------------------------------------------------------------
!
! MODULE: materials
!
!> @author / date
!> S. Guinchard - EPFL/SPC / 11-25-22
!> P. Giroud-Garampon - EPFL/SPC / 13-02-2024
!
! DESCRIPTION:
!> Module storing data related to SEE and IIEE for different materials and gases.
!------------------------------------------------------------------------------
MODULE materials
USE constants
IMPLICIT NONE
!--------------------- Coefficients for the IIEE yield ---------------------
! The coefficients for Hydrogen are obtained by polynomial piecewise fit
! of the energy loss function dE/dx
! N.B. they fit dE/dx in MeV/cm
! The coefficients for other gases are obtained by piecewise fit
! of the IIEE yield
! 1. 304 STAINLESS STEEL
! 1. Hydrogen
REAL(KIND=db), DIMENSION(4), PARAMETER :: iiee_yield_H_SS_1 = &
(/0.0, 2.9778E02, 0.0, 0.0/)
REAL(KIND=db), DIMENSION(4), PARAMETER :: iiee_yield_H_SS_2 = &
(/1.834520315818557E02, 1.320304216355084E05, -8.583700602370013E06, 3.526140145560557E08/)
REAL(KIND=db), DIMENSION(4), PARAMETER :: iiee_yield_H_SS_3 = &
(/2.471679999999999E02, 9.695466666666670E04, -2.475200000000003E06, 2.325333333333340E07/)
REAL(KIND=db), DIMENSION(4), PARAMETER :: iiee_yield_H_SS_4 = &
(/2.533904454349683E03, -1.766016382825937E05, 8.202640024592019E06, -1.125320217235288E08/)
REAL(KIND=db), DIMENSION(4), PARAMETER :: iiee_yield_H_SS_5 = &
(/8.786057142856745E02, 2.856595238095524e+04, -1.834285714286391E05, 3.333333333338620E05/)
! 2. COPPER
! 1. Hydrogen
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Cu_1 = &
(/0.0, 2.812300E02, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Cu_2 = &
(/1.703970762187365E02, 1.273920899063999E05, -8.680492697427949E06, 3.535110244304671E08/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Cu_3 = &
(/2.618095714285664E02, 8.520814285714373E04, -2.094971428571471E06, 2.064000000000054E07/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Cu_4 = &
(/1.030650000000107E03, -6.231904761918026E03, 1.468571428571969E06, -2.506666666667393E07/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Cu_5 = &
(/7.117599999999986E02, 3.874166666666669E04, -5.099999999999992E05, 2.733333333333320E06/)
! 2. Helium (from https://doi.org/10.1103/PhysRevB.19.121)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Cu_1 = &
(/0.293714285714, 32.0952380952, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Cu_2 = &
(/0.125, 1.1154761904762E02, -9.166666666675E03, 4.761904761905E05/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Cu_3 = &
(/2.0E-02, 9.8E01, -2.0E03, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Cu_4 = &
(/6.6E-01, 2.6E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Cu_5 = &
(/6.6E-01, 2.6E01, 0.0, 0.0/)
! 3. Neon (from https://doi.org/10.1016/0039-6028(85)90136-0)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Cu_1 = &
(/0.225142857142857, 1.66612665684839, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Cu_2 = &
(/-0.10138381040292, 2.04712802207033E02, -2.12122342024929E04, 1.04199718472196E06/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Cu_3 = &
(/1.50158859504749, -1.99666130090673E02, 1.72321487767404E04, -3.60132754302814E05/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Cu_4 = &
(/1.04377749629525, 2.40390976067719E01, 0.0, 0.0/)
! 4. Argon (from https://doi.org/10.1016/0039-6028(79)90341-8)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Cu_1 = &
(/9.02857142857143E-02, 1.16100605299447E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Cu_2 = &
(/-4.555096430036E-02, 1.00774996605114E02, 7.7149611879632E01, -5.76830557431987E04/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Cu_3 = &
(/-2.5992776033708E-02, 1.07801648650678E02, -1.41855662651229E03, 4.91523169780497E03/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Cu_4 = &
(/-3.87791890832554E-01, 1.48442167310384E02, -3.1298539515514E03, 3.39747092112027E04/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Cu_5 = &
(/5.78008566722373E-01, 5.29272726756834E01, 0.0, 0.0/)
! 3. ALUMINUM
! 1. Hydrogen
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Al_1 = &
(/0.0, 2.100900000000000E02, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Al_2 = &
(/1.224483183711612E02, 9.731277035468460E04, -6.878546755544925E06, 2.795627748449915E08/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Al_3 = &
(/1.998525714285755E02, 6.229714285714197E04, -1.531771428571365E06, 1.567999999999851E07/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Al_4 = &
(/8.135440000000424E02, -1.333600000000502E04, 1.553600000000198E06, -2.624000000000260E07/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_H_Al_5 = &
(/4.052485714285687E02, 3.974642857142875E04, -6.631428571428614E05, 3.800000000000037E06/)
! 2. Helium (from https://doi.org/10.1103/PhysRevB.19.121)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Al_1 = &
(/0.190541702493551, 4.47291487532244E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Al_2 = &
(/1.76659521215925E-01, 5.34536358214894E01, -6.19216675102752E02, -2.93593275345582E04/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Al_3 = &
(/2.55E-01, 4.15E01, -5.00E02, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Al_4 = &
(/4.55E-01, 2.15E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_He_Al_5 = &
(/4.7E-01, 2.1E01, 0.0, 0.0/)
! 3. Neon (from https://doi.org/10.1103/PhysRevB.22.80)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Al_1 = &
(/0.14926913155632, 1.8132712892262E02, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Al_2 = &
(/0.290397887461225, 1.4159935798866E02, -1.19943419298762E04, 4.5403531271666E05/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Al_3 = &
(/0.741599478768533, 2.24368861780615E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Al_4 = &
(/0.731551934039544, 2.30818722295553E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ne_Al_5 = &
(/1.4059696031916, -8.36542116269722, 3.04658729274448E02, 0.0/)
! 4. Argon (https://doi.org/10.1016/0039-6028(79)90341-8)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Al_1 = &
(/6.80997420464317E-02, 3.90589439790723E01, 0.0, 0.0/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Al_2 = &
(/9.54897034753516E-02, -1.76674600320705E01, 1.54604387031526E04, -7.64344295477741E05/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Al_3 = &
(/-3.26838161662861E-01, 1.46366287920992E02, -5.18224285283153E03, 8.11411852216318E04/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Al_4 = &
(/5.51191348971975E-02, 7.50010618137324E01, -1.14344429560049E03, 9.72498529527904E03/)
REAL(KIND = db), DIMENSION(4), PARAMETER :: iiee_yield_Ar_Al_5 = &
(/5.25429038367765E-01, 3.37724550898206E01, 0.0, 0.0/)
!--------------- Coefficients for the SEE backscattering coefficient ----------------
! 3. ALUMINUM
REAL(KIND = db), PARAMETER :: back_coef_Al_1 = 0.0
REAL(KIND = db), DIMENSION(3), PARAMETER :: back_coef_Al_2 = &
(/-1.67E-03, 2.45E-03, -6.33E-06/)
REAL(KIND = db), PARAMETER :: back_coef_Al_3 = 0.235
!-------------------------- Parameters for the SEE yield ----------------------------
! 3. ALUMINUM
REAL(KIND = db), PARAMETER :: delta_m_Al = 1.034
REAL(KIND = db), PARAMETER :: E_m_Al = 275 ! [eV]
REAL(KIND = db), PARAMETER :: s_Al = 1.7475 ! fitting parameter for the yield curve
END MODULE materials
c4science · Help