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 | 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.
https://hdl.handle.net/11583/2986756