CPL - Chalmers Publication Library
| Utbildning | Forskning | Styrkeområden | Om Chalmers | In English In English Ej inloggad.

Uncertainty estimation and a stopping rule in nonlinear gyrokinetic simulations

Michael Oberparleiter (Institutionen för rymd- och geovetenskap, Plasmafysik och fusionsenergi) ; Hans Nordman (Institutionen för rymd- och geovetenskap, Plasmafysik och fusionsenergi) ; Geert Verdoolaege ; Frank Jenko
Journal of Physics: Conference Series, Joint Varenna-Lausanne International Workshop on the Theory of Fusion Plasmas 2016 (1742-6588). Vol. 775 (2016), p. 012009.
[Konferensbidrag, refereegranskat]

We present a method to estimate the mean and uncertainty of fluctuating quantities, such as spatially averaged density and temperature fluctuations or radial fluxes, from initial value simulations of the Eulerian gyrokinetic code GENE[1, 2]. Since the time series are autocorrelated in time, the data is grouped into batches based on the autocorrelation time and their means form the sample for further statistical treatment, such as calculating the standard error of the mean. Based on this uncertainty estimate we develop a stopping rule for a nonlinear simulation: First, regression tests ensure that it has reached a stationary (quasisteady) state and data before this point is discarded. Then the previously described estimate is calculated. If the estimated relative error is below a prescribed threshold, the simulation is stopped. This scheme is applied to several previously performed GENE simulations ranging from simple benchmarks to modelling of JET and ASDEX discharges. It can be demonstrated that a number of simulations could be around 30% shorter if a maximal statistical relative uncertainty of 5% is desired for all monitored quantities.

Denna post skapades 2017-01-16. Senast ändrad 2017-10-10.
CPL Pubid: 247070


Läs direkt!

Lokal fulltext (fritt tillgänglig)

Länk till annan sajt (kan kräva inloggning)

Institutioner (Chalmers)

Institutionen för rymd- och geovetenskap, Plasmafysik och fusionsenergi (2013-2017)


Plasmafysik med fusion

Chalmers infrastruktur

C3SE/SNIC (Chalmers Centre for Computational Science and Engineering)