Spectral features in ADAS Martin O’Mullane
ADAS-EU course, CEA, Cadarache, 3-October-2012
Special/Spectral features and tools for modelling them Not all spectra can be treated as isolated lines I
Quasi-continuum of high-Z ions.
I
Zeeman effect.
I
Balmer series.
I
He-like spectrum.
I
Stark split spectrum.
All are controlled by a limited set of inputs which enables the feature to be embedded into a parametrized fitting scheme. Spectral fitting tools for some of these features are being introduced into ADAS — many still in development so testers required
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Normal line emission Z
Rmax
PEC(r)Ne(r)Nz (r)ϕj→k (λ) dr
I(λ) = Rmin
The simplest line ‘feature’ is built-up from: I
Population calculation (adas208) for PEC
I
Measurement of electron density (and temperature) profiles
I
Ionisation balance along line of sight (Nz ) — either transport code or equilibrium balance (adas405)
I
Spectral profile — plasma and instrumentation (ϕj→k (λ))
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Feature photon emissivity coefficients The F − PEC is an extension of the familiar line emissivity coefficient. For an individual line: W
−13 1·10
10+
εj→k = Aj→k
M X
(exc)
Fjσ
NeNσ
σ=1
f−PEC (ph cm 3s−1 )
−14 8·10
Take a spectral interval, [λ0, λ1], subdivided into N pixels and a normalised emission profile, ϕj→k (λ) for each line.
−14 6·10 −14 4·10 −14 2·10
The envelope coefficient is:
0
feature
10000.0
(exc)[0,1]
1000.0 100.0 Te (eV)
photon
10.0 1.0 0.1
600
700
800 λ (Å)
900
1000
F PEC σ,i X j,k;λj→k ∈[0,1]
ADAS-EU course, CEA, Cadarache, 3-October-2012
emissivity
=
(exc) PEC σ,j→k
Z
λi+1
ϕj→k (λ) dλ, λi
www.adas.ac.uk
Generating spectral and power data — ADAS810 An adf42 driver file, made by adas8xx create drivers.pro, is used by ADAS810 to generate adf11/plt, adf15 and adf40 datasets.
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADAS801 — review or enter processing parameters
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
IDL (and fortran) routines to access F − PEC data
IDL> read_adf40,file=’/home/adas/adas/adf40/fpec40#w/fpec40#w_ic#w10.dat’, fulldata=al IDL> help, all, /st ESYM IZ0 IS IS1 NBLOCK NPIX WAVE_MIN WAVE_MAX NTE TE NDENS DENS FPEC TYPE
STRING LONG LONG LONG LONG LONG DOUBLE DOUBLE LONG DOUBLE LONG DOUBLE DOUBLE STRING
’W ’ 74 10 11 3 Array[3] Array[3] Array[3] Array[3] Array[14, 3] Array[3] Array[7, 3] Array[512, 14, 7, 3] Array[3]
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Interpolate onto our plasma conditions of interest IDL> wave=adas_vector(low=all.wave_min[0], high=all.wave_max[0], $ num=all.npix[0], /linear) IDL> file = ’/home/adas/adas/adf40/fpec40#w/fpec40#w_ic#w10.dat’ ILD> read_adf40, file=file, block=1, te=te, dens=dens, data=fpec_w10
For fortran call xxdata 40.for to read in the data and e3spln.for to interpolate these on the temperature and density pairs of interest.
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Hyperfine effects ADAS603 — Zeeman/Paschen-Back effect The spectral feature is the connected set of spectrum lines from magnetic field perturbation of multiplet components for a specified LS → L0S 0. I
λLSJM −L0S 0J 0M 0
I
Depends on magnetic field B.
I
Polarisation depends on whether M = M 0, M 0 + 1, M 0 − 1
I
Emissivities are proportional ro viewing angle.
Three IDL routines to retrieve the component data: I
adas603 get multiplets.pro lists available elements and transitions.
I
adas603 get hdlike.pro for H-like (H, D, T, He, C and Ne)
I
adas603 get zeeman.pro for other ions
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Hyperfine effects ADAS305 — Stark split beam emission A non-interactive code which calculates the hydrogen Balmer-α spectrum split into components due to the plasma magnetic field and arbitrary electric fields.
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Column emissivity (ph cm−2 s−1 A−1)
2.5•1013
2.0•1013
1.5•1013
1.0•1013
5.0•1012
0 6510
6520
ADAS-EU course, CEA, Cadarache, 3-October-2012
6530 Wavelength (A)
6540
6550
www.adas.ac.uk
Hyperfine effects — explore with ADAS605 A pedagogical tool to explore the response to feature inputs
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
Other features
I
CALCAT/Pickett code for diatomic band features (not in central ADAS) — ADAS series 9 is a more complete picture
I
Atomic data for He-like DR feature present in ADAS — specialized population code adas218 awaiting integration.
I
PPP code for series has been used in the past.
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADS602 — line fitting code
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADS603 — Zeeman feature line fitting code
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk
ADAS future development — afg/ffs afg — ADAS feature generator ffs — functions for spectral fitting
ADAS-EU course, CEA, Cadarache, 3-October-2012
www.adas.ac.uk