16,473 research outputs found
Robust scaling in fusion science: case study for the L-H power threshold
In regression analysis for deriving scaling laws in the context of fusion studies, standard regression methods are usually applied, of which ordinary least squares (OLS) is the most popular. However, concerns have been raised with respect to several assumptions underlying OLS in its application to fusion data. More sophisticated statistical techniques are available, but they are not widely used in the fusion community and, moreover, the predictions by scaling laws may vary significantly depending on the particular regression technique. Therefore we have developed a new regression method, which we call geodesic least squares regression (GLS), that is robust in the presence of significant uncertainty on both the data and the regression model. The method is based on probabilistic modeling of all variables involved in the scaling expression, using adequate probability distributions and a natural similarity measure between them (geodesic distance). In this work we revisit the scaling law for the power threshold for the L-to-H transition in tokamaks, using data from the multi-machine ITPA databases. Depending on model assumptions, OLS can yield different predictions of the power threshold for ITER. In contrast, GLS regression delivers consistent results. Consequently, given the ubiquity and importance of scaling laws and parametric dependence studies in fusion research, GLS regression is proposed as a robust and easily implemented alternative to classic regression techniques
Modeling and interpolation of the ambient magnetic field by Gaussian processes
Anomalies in the ambient magnetic field can be used as features in indoor
positioning and navigation. By using Maxwell's equations, we derive and present
a Bayesian non-parametric probabilistic modeling approach for interpolation and
extrapolation of the magnetic field. We model the magnetic field components
jointly by imposing a Gaussian process (GP) prior on the latent scalar
potential of the magnetic field. By rewriting the GP model in terms of a
Hilbert space representation, we circumvent the computational pitfalls
associated with GP modeling and provide a computationally efficient and
physically justified modeling tool for the ambient magnetic field. The model
allows for sequential updating of the estimate and time-dependent changes in
the magnetic field. The model is shown to work well in practice in different
applications: we demonstrate mapping of the magnetic field both with an
inexpensive Raspberry Pi powered robot and on foot using a standard smartphone.Comment: 17 pages, 12 figures, to appear in IEEE Transactions on Robotic
Surrogate Accelerated Bayesian Inversion for the Determination of the Thermal Diffusivity of a Material
Determination of the thermal properties of a material is an important task in
many scientific and engineering applications. How a material behaves when
subjected to high or fluctuating temperatures can be critical to the safety and
longevity of a system's essential components. The laser flash experiment is a
well-established technique for indirectly measuring the thermal diffusivity,
and hence the thermal conductivity, of a material. In previous works,
optimization schemes have been used to find estimates of the thermal
conductivity and other quantities of interest which best fit a given model to
experimental data. Adopting a Bayesian approach allows for prior beliefs about
uncertain model inputs to be conditioned on experimental data to determine a
posterior distribution, but probing this distribution using sampling techniques
such as Markov chain Monte Carlo methods can be incredibly computationally
intensive. This difficulty is especially true for forward models consisting of
time-dependent partial differential equations. We pose the problem of
determining the thermal conductivity of a material via the laser flash
experiment as a Bayesian inverse problem in which the laser intensity is also
treated as uncertain. We introduce a parametric surrogate model that takes the
form of a stochastic Galerkin finite element approximation, also known as a
generalized polynomial chaos expansion, and show how it can be used to sample
efficiently from the approximate posterior distribution. This approach gives
access not only to the sought-after estimate of the thermal conductivity but
also important information about its relationship to the laser intensity, and
information for uncertainty quantification. We also investigate the effects of
the spatial profile of the laser on the estimated posterior distribution for
the thermal conductivity
- …