Bayesian Inference for Microwave Diagnostics at Joint European Torus

• The major goal of this thesis is to perform Bayesian inference jointly for electron kinetic profiles, a toroidal magnetic field correction, diagnostic sensitivities, and wall reflection properties given noisy measurements of four microwave diagnostics at the tokamak JET. Besides the measurements, this kind of inference considers prior knowledge, like profile length-scales, to update objectively the information about the physics parameters. Such a probabilistic update, i.e. a posterior probability distribution or posterior, states the plausibility of parameter combinations and captures parameter uncertainties and correlations. The already existing Bayesian framework Minerva was used to carry out the inference.Contrary to standard approaches, plasma physics and diagnostic models and physics parameters are used to make predictions which are compared objectively with the data provided by one reflectometer and three electron cyclotron emission (ECE) diagnostics; two broadband Martin-Puplett interferometers and one heterodyne radiometer. In addition, these models represent reality more closely, for instance: (i) the smoothnesses of temperature and density profiles are modelled and estimated by length-scale parameters, (ii) the models SPECE and ECEPT, predicting each broadband ECE spectra, take into account relativistic and density effects like the finite optical thickness, (iii) the working principle and measurement uncertainty of each diagnostic are considered. As an example for the latter point, the two interferometers supply convoluted ECE spectra over several harmonic ranges up to 500 GHz and in mainly ordinary and extra-ordinary wave mode polarisations. The uncertainties on these spectra originate mostly from absolute calibrations for which a dedicated, robust and reliable procedure had to be established during this thesis.Further intermediate achievements of this thesis are for example: (i) the analytical derivation of the generalised square-exponential covariance function, enabling the estimation of multiple length-scales for electron temperature and density profiles in the plasma core and edge domains, (ii) the already existing ray-tracer SPECE, which predicts accurately but slowly ECE spectra, was parallelised successfully by a client-server approach, (iii) the derivation of the model ECEPT to predict quickly and sufficiently accurately broadband ECE spectra for an Ohmically heated plasma, and (iv) the extension of the multi-reflection model to allow different properties for the high-field and low-field side walls of a fusion device, affecting substantially the predictions of ECE spectra at frequencies for which the plasma has a low optical thickness.For a low temperature and low density plasma, the joint inference of more than 200 parameters was carried out for a given flux surface geometry, using either SPECE or ECEPT as predictor for the measured ECE spectra from 50 GHz to 280 GHz (first three harmonic ranges at least) with contributions from the ordinary and extra-ordinary wave-modes. For the ECEPT case, the shape of the joint posterior could be explored numerically. Related findings are, for example, smooth electron temperature and density profiles with values at the centre of 1.5 keV and 1.75x1e19 m^-3 and at the separatrix of 80 eV and 2x1e18 m^-3 with uncertainties of the order of 10 eV and few 1e17 m^-3. Furthermore, the smoothness of each core and edge profile originates in the inferred length-scale. These results were confirmed by the most likely parameter combination of the posterior when the more accurate predictor SPECE is used. The only exceptions are the correction to the magnetic field, which increased from 1.4% to 2%, and the reflectivity of the ITER-like wall, which elevated from 0.72 to 0.92.With ECEPT in use, multiple posterior correlation features could be revealed of which many could be explained. For instance, each profile smoothness, due to the estimated length-scale, implies considerable correlationsbetween nearby locations. In addition, a global correlation feature is evidentfor edge and core density profiles, likely caused by the measurement principle of the reflectometer.The second part of this thesis focusses on Bayesian inference about spectra in the field of Fourier transform spectroscopy. An example application has been carried out for the calibration data measured with one of the Martin-Puplett interferometers at JET. Compared to standard analysis techniques, more information could be extracted about the spectra from the so-called double-, single- and zero-sided data domains. The spectra have been modelled as Brownian bridge processes a priori, by which the global trend of the data located in double- and single-sided domains is modelled. This enables the estimation of the posterior uncertainties of the spectra due to non-probed data domains, especially the zero-sided domain. Such consideration is not made in any form by conventional techniques, because no method has been developed to estimate the seemingly lost information contained in the zero-sided domain. The lower and upper limits of the spectra have been estimated by (33+/-1.7) GHz and (913+/-2.9) GHz. These limits are rated to be more likely by a factor of 1e277 than the limits at 0 GHz and 3747 GHz (Nyquist frequency) assumed by conventional analysis approaches but never checked for plausibility. 

• Det huvudsakliga målet med denna avhandling är att utföra Bayesiansk inferens gemensamt för elektronernas kinetiska profiler, en toroidal magnetfältskorrigering, diagnostiska känsligheter och väggreflektionsegenskaper givet brusiga mätningar från fyra mikrovågsdiagnostiker vid tokamaken JET. Förutom mätningarna tar denna typ av inferens hänsyn till tidigare kunskap, t.ex. profilernas längdskalor, för att objektivt uppdatera informationen om fysikaliska parametrar. En sådan probabilistisk uppdatering, dvs. en efterföljande sannolikhetsfördelning eller posterior, anger sannolikheten för parameterkombinationer och fångar upp parameterosäkerheter och korrelationer. Det redan existerande bayesianska ramverket Minerva användes för att genomföra inferensen. I motsats till standardmetoder används modeller för plasmafysik och diagnostik samt fysikaliska parametrar för att göra förutsägelser som jämförs objektivt med data från en reflektometer och tre ECE-diagnostiker (elektroncyklotronemission); två bredbandiga Martin-Puplett-interferometrar och en heterodynradiometer. Dessutom ger dessa modeller en nära beskrivning av verkligheten, till exempel: (ii) modellerna SPECE och ECEPT, som förutsäger varje bredbandigt ECE-spektra, tar hänsyn till relativistiska effekter och densitetseffekter som den ändliga optiska tjockleken, (iii) arbetsprincipen och mätosäkerheten för varje diagnostik beaktas. Som exempel på den sistnämnda punkten kan nämnas att de två interferometrarna levererar sammanflätade ECE-spektra över flera harmoniska områden upp till 500 GHz och i huvudsakligen ordinära och extraordinära vågmoder. Osäkerheterna i dessa spektrum beror främst på absoluta kalibreringar för vilka ett särskilt, robust och tillförlitligt förfarande behövde fastställas under detta avhandlingsarbete.Andra resultat i denna avhandling är till exempel: (i) den analytiska härledningen av den generaliserade kvadratisk-exponentiella kovariansfunktionen, som gör det möjligt att uppskatta flera längdskalor för elektronernas temperatur- och densitetsprofiler i plasmats kärn- och kantområden, (ii) den redan existerande strålspåraren SPECE, som förutsäger exakta ECE-spektra, men är tidskrävande, parallelliserades framgångsrikt med hjälp av en klient-servermetod, (iii) härledning av modellen ECEPT för att snabbt och tillräckligt exakt förutsäga bredbandiga ECE-spektra för ett ohmskt uppvärmt plasma, och (iv) utvidgning av multireflektionsmodellen för att tillåta olika egenskaper för högfälts- och lågfälts-sidoväggarna i en fusionsanläggning, vilket i hög grad påverkar förutsägelserna av ECE-spektra vid frekvenser för vilka plasman har en liten optisk tjocklek.För ett plasma med låg temperatur och låg densitet genomfördes en gemensam uppskattning av mer än 200 parametrar för en given geometri av flödesytan, med antingen SPECE eller ECEPT som prediktor för de uppmätta ECE-spektra från 50 GHz till 280 GHz (åtminstone de tre första harmoniska områdena) med bidrag från vanliga och extraordinära vågmoder. För ECEPT-fallet kan formen på den gemensamma posterior undersökas numeriskt. Relaterade resultat är till exempel jämna profiler för elektrontemperatur och elektrontäthet med värden i mitten av 1,5 keV och 1,75×1e19 m^−3 och vid separatrisen av 80 eV och 2×1e18 m^−3 med osäkerheter i storleksordningen 10 eV och några 1e17 m^−3. Dessutom har jämnheten hos varje kärn- och kantprofil sitt ursprung i den härledda längdskalan. Dessa resultat bekräftades av den mest sannolika parameterkombinationen av posterior när den mer exakta prediktorn SPECE används. De enda undantagen är korrigeringen av magnetfältet, som ökade från 1,4% till 2%, och reflektiviteten hos den ITER-liknande väggen, som ökade från 0,72 till 0,92. Med ECEPT kunde flera posterior korrelationsegenskaper upptäckas, varav vissa kunde förklaras. Till exempel innebär varje profils jämnhet, på grund av den uppskattade längdskalan, betydande korrelationer mellan närliggande platser. Dessutom är det en global korrelation för profilers kant- och kärntäthet, vilket troligen beror på reflektometerns mätprincip. Den andra delen av denna avhandling fokuserar på Bayesiansk inferens av spektrum inom området Fouriertransformspektroskopi. Ett exempel på tillämpning som utförts är kalibreringsdata som uppmätts med en av Martin-Puplett-interferometrarna vid JET. Jämfört med vanliga analysmetoder kunde mer information om spektra utvinnas från de så kallade dubbel-, enkel- och nollsidiga datadomänerna. Spektra har på förhand modellerats som "Brownian-bridge" processer, varigenom den globala trenden för de data som finns i dubbel- och enkelsidiga områden modelleras. Detta gör det möjligt att uppskatta posterior osäkerheterna i spektra på grund av ej uppmätta datadomäner, särskilt den nollsidiga domänen. Sådana överväganden görs inte i någon form med konventionella tekniker, eftersom ingen metod har utvecklats för att uppskatta den till synes förlorade informationen i den nollsidiga domänen. Spektras nedre och övre gränser har uppskattats till (33 ± 1,7) GHz och (913 ± 2,9) GHz. Dessa gränser bedöms vara mer sannolika med en faktor 1e277 än gränserna 0 GHz och 3747 GHz (Nyquist-frekvens) som antas i konventionella analysmetoder men vars rimlighet aldrig har kontrollerats.

Ämnesord

NATURVETENSKAP  -- Fysik -- Fusion, plasma och rymdfysik (hsv//swe)

NATURAL SCIENCES  -- Physical Sciences -- Fusion, Plasma and Space Physics (hsv//eng)

Nyckelord

ECE

Bayesian

Fourier transform

spectroscopy

joint inference

length-scale

Gaussian process

microwave diagnostic

JET

Electrical Engineering

Elektro- och systemteknik

Publikations- och innehållstyp

vet (ämneskategori)

dok (ämneskategori)