A novel approach to uncertainty quantification in codes simulating fusion plasma, deterministic sampling (DS), is evaluated. This method uses a few carefully selected samples and can be used to propagate input parameter uncertainties through calculations where other sampling methods are unmanageable due to time constraints. The primary analysis is performed on the transport code TRANSP, but another faster code is also tested where a comparison with Monte Carlo sampling is made. The tests, performed with two pulses at the Joint European Torus (JET), show that even lower order DS will give a reliable estimation of the standard deviation of the calculated neutron rate. However, a higher order DS can give information about higher output moments, such as skewness an kurtosis. The TRANSP-simulated neutron rates of both examined pulses are found to have an uncertainty with an upward skewness, meaning input parameter uncertainties are can better explain an underestimation of the neutron rate than an overestimation. This information can, for example, be lead to a better benchmarking comparison between the measured and calculated neutron rates.

Propagating transport-code input parameter uncertainties with deterministic sampling / Sahlberg, A; Hellesen, C; Eriksson, J; Conroy, S; Ericsson, G; King, D; Subba, F. - In: PLASMA PHYSICS AND CONTROLLED FUSION. - ISSN 0741-3335. - 60:12(2018). [10.1088/1361-6587/aae80b]

Propagating transport-code input parameter uncertainties with deterministic sampling

Subba, F
2018

Abstract

A novel approach to uncertainty quantification in codes simulating fusion plasma, deterministic sampling (DS), is evaluated. This method uses a few carefully selected samples and can be used to propagate input parameter uncertainties through calculations where other sampling methods are unmanageable due to time constraints. The primary analysis is performed on the transport code TRANSP, but another faster code is also tested where a comparison with Monte Carlo sampling is made. The tests, performed with two pulses at the Joint European Torus (JET), show that even lower order DS will give a reliable estimation of the standard deviation of the calculated neutron rate. However, a higher order DS can give information about higher output moments, such as skewness an kurtosis. The TRANSP-simulated neutron rates of both examined pulses are found to have an uncertainty with an upward skewness, meaning input parameter uncertainties are can better explain an underestimation of the neutron rate than an overestimation. This information can, for example, be lead to a better benchmarking comparison between the measured and calculated neutron rates.
File in questo prodotto:
File Dimensione Formato  
61291_uq_article.pdf

accesso aperto

Tipologia: 1. Preprint / submitted version [pre- review]
Licenza: Pubblico - Tutti i diritti riservati
Dimensione 3.67 MB
Formato Adobe PDF
3.67 MB Adobe PDF Visualizza/Apri
Sahlberg_2018_Plasma_Phys._Control._Fusion_60_125010.pdf

accesso riservato

Tipologia: 2a Post-print versione editoriale / Version of Record
Licenza: Non Pubblico - Accesso privato/ristretto
Dimensione 1.16 MB
Formato Adobe PDF
1.16 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11583/2986756