4,514 research outputs found
Coordinate Transformation and Polynomial Chaos for the Bayesian Inference of a Gaussian Process with Parametrized Prior Covariance Function
This paper addresses model dimensionality reduction for Bayesian inference
based on prior Gaussian fields with uncertainty in the covariance function
hyper-parameters. The dimensionality reduction is traditionally achieved using
the Karhunen-\Loeve expansion of a prior Gaussian process assuming covariance
function with fixed hyper-parameters, despite the fact that these are uncertain
in nature. The posterior distribution of the Karhunen-Lo\`{e}ve coordinates is
then inferred using available observations. The resulting inferred field is
therefore dependent on the assumed hyper-parameters. Here, we seek to
efficiently estimate both the field and covariance hyper-parameters using
Bayesian inference. To this end, a generalized Karhunen-Lo\`{e}ve expansion is
derived using a coordinate transformation to account for the dependence with
respect to the covariance hyper-parameters. Polynomial Chaos expansions are
employed for the acceleration of the Bayesian inference using similar
coordinate transformations, enabling us to avoid expanding explicitly the
solution dependence on the uncertain hyper-parameters. We demonstrate the
feasibility of the proposed method on a transient diffusion equation by
inferring spatially-varying log-diffusivity fields from noisy data. The
inferred profiles were found closer to the true profiles when including the
hyper-parameters' uncertainty in the inference formulation.Comment: 34 pages, 17 figure
Numerical approximation of poroelasticity with random coefficients using Polynomial Chaos and Hybrid High-Order methods
In this work, we consider the Biot problem with uncertain poroelastic
coefficients. The uncertainty is modelled using a finite set of parameters with
prescribed probability distribution. We present the variational formulation of
the stochastic partial differential system and establish its well-posedness. We
then discuss the approximation of the parameter-dependent problem by
non-intrusive techniques based on Polynomial Chaos decompositions. We
specifically focus on sparse spectral projection methods, which essentially
amount to performing an ensemble of deterministic model simulations to estimate
the expansion coefficients. The deterministic solver is based on a Hybrid
High-Order discretization supporting general polyhedral meshes and arbitrary
approximation orders. We numerically investigate the convergence of the
probability error of the Polynomial Chaos approximation with respect to the
level of the sparse grid. Finally, we assess the propagation of the input
uncertainty onto the solution considering an injection-extraction problem.Comment: 30 pages, 15 Figure
A non-adapted sparse approximation of PDEs with stochastic inputs
We propose a method for the approximation of solutions of PDEs with
stochastic coefficients based on the direct, i.e., non-adapted, sampling of
solutions. This sampling can be done by using any legacy code for the
deterministic problem as a black box. The method converges in probability (with
probabilistic error bounds) as a consequence of sparsity and a concentration of
measure phenomenon on the empirical correlation between samples. We show that
the method is well suited for truly high-dimensional problems (with slow decay
in the spectrum)
Motion Planning of Uncertain Ordinary Differential Equation Systems
This work presents a novel motion planning framework, rooted in nonlinear programming theory, that treats uncertain fully and under-actuated dynamical systems described by ordinary differential equations. Uncertainty in multibody dynamical systems comes from various sources, such as: system parameters, initial conditions, sensor and actuator noise, and external forcing. Treatment of uncertainty in design is of paramount practical importance because all real-life systems are affected by it, and poor robustness and suboptimal performance result if it’s not accounted for in a given design. In this work uncertainties are modeled using Generalized Polynomial Chaos and are solved quantitatively using a least-square collocation method. The computational efficiency of this approach enables the inclusion of uncertainty statistics in the nonlinear programming optimization process. As such, the proposed framework allows the user to pose, and answer, new design questions related to uncertain dynamical systems.
Specifically, the new framework is explained in the context of forward, inverse, and hybrid dynamics formulations. The forward dynamics formulation, applicable to both fully and under-actuated systems, prescribes deterministic actuator inputs which yield uncertain state trajectories. The inverse dynamics formulation is the dual to the forward dynamic, and is only applicable to fully-actuated systems; deterministic state trajectories are prescribed and yield uncertain actuator inputs. The inverse dynamics formulation is more computationally efficient as it requires only algebraic evaluations and completely avoids numerical integration. Finally, the hybrid dynamics formulation is applicable to under-actuated systems where it leverages the benefits of inverse dynamics for actuated joints and forward dynamics for unactuated joints; it prescribes actuated state and unactuated input trajectories which yield uncertain unactuated states and actuated inputs.
The benefits of the ability to quantify uncertainty when planning the motion of multibody dynamic systems are illustrated through several case-studies. The resulting designs determine optimal motion plans—subject to deterministic and statistical constraints—for all possible systems within the probability space
Open TURNS: An industrial software for uncertainty quantification in simulation
The needs to assess robust performances for complex systems and to answer
tighter regulatory processes (security, safety, environmental control, and
health impacts, etc.) have led to the emergence of a new industrial simulation
challenge: to take uncertainties into account when dealing with complex
numerical simulation frameworks. Therefore, a generic methodology has emerged
from the joint effort of several industrial companies and academic
institutions. EDF R&D, Airbus Group and Phimeca Engineering started a
collaboration at the beginning of 2005, joined by IMACS in 2014, for the
development of an Open Source software platform dedicated to uncertainty
propagation by probabilistic methods, named OpenTURNS for Open source Treatment
of Uncertainty, Risk 'N Statistics. OpenTURNS addresses the specific industrial
challenges attached to uncertainties, which are transparency, genericity,
modularity and multi-accessibility. This paper focuses on OpenTURNS and
presents its main features: openTURNS is an open source software under the LGPL
license, that presents itself as a C++ library and a Python TUI, and which
works under Linux and Windows environment. All the methodological tools are
described in the different sections of this paper: uncertainty quantification,
uncertainty propagation, sensitivity analysis and metamodeling. A section also
explains the generic wrappers way to link openTURNS to any external code. The
paper illustrates as much as possible the methodological tools on an
educational example that simulates the height of a river and compares it to the
height of a dyke that protects industrial facilities. At last, it gives an
overview of the main developments planned for the next few years
Surrogate and reduced-order modeling: a comparison of approaches for large-scale statistical inverse problems [Chapter 7]
Solution of statistical inverse problems via the frequentist or Bayesian approaches described in earlier chapters can be a computationally intensive endeavor, particularly when faced with large-scale forward models characteristic of many engineering and science applications. High computational cost arises in several ways. First, thousands or millions of forward simulations may be required to evaluate estimators of interest or to characterize a posterior distribution. In the large-scale setting, performing so many forward simulations is often computationally intractable. Second, sampling may be complicated by the large dimensionality of the input space--as when the inputs are fields represented with spatial discretizations of high dimension--and by nonlinear forward dynamics that lead to multimodal, skewed, and/or strongly correlated posteriors. In this chapter, we present an overview of surrogate and reduced order modeling methods that address these computational challenges. For illustration, we consider a Bayesian formulation of the inverse problem. Though some of the methods we review exploit prior information, they largely focus on simplifying or accelerating evaluations of a stochastic model for the data, and thus are also applicable in a frequentist context.Sandia National Laboratories (Laboratory Directed Research and Development (LDRD) program)United States. Dept. of Energy (Contract DE-AC04-94AL85000)Singapore-MIT Alliance Computational Engineering ProgrammeUnited States. Dept. of Energy (Award Number DE-FG02-08ER25858 )United States. Dept. of Energy (Award Number DESC00025217
- …