Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2009 B. Schrefler, E. O˜ nate and M. Papadrakakis (Eds) c

CIMNE, Barcelona, 2009

A ROBUST AND EFFICIENT UNCERTAINTY QUANTIFICATION METHOD FOR COUPLED FLUID-STRUCTURE INTERACTION PROBLEMS JEROEN A.S. WITTEVEEN∗ AND HESTER BIJL† Faculty of Aerospace Engineering, Delft University of Technology Kluyverweg 1, 2629HS Delft, The Netherlands ∗ e-mail: † e-mail:

[email protected], web page: http://www.jeroenwitteveen.com [email protected], web page: http://www.lr.tudelft.nl/aerodynamics

Key words: Uncertainty quantification, Fluid-structure interaction, Polynomial Chaos, Extrema diminishing, Error bounds. Abstract. A robust and efficient uncertainty quantification method is presented for resolving the effect of uncertainty on the behavior of multi-physics systems. The extrema diminishing method in probability space maintains a bounded error due to the interpolation of deterministic samples at constant phase in a transonic airfoil flutter problem.

1

INTRODUCTION

Numerical errors in multi-physics simulations start to reach acceptable engineering accuracy levels due to the increasing availability of computational resources. Nowadays, uncertainties in multi-scale models, physical parameter variations, and lack of knowledge of initial and boundary conditions have a larger effect on computational predictions than discretization errors. It is, therefore, vital to take these uncertainties into account in coupled problems to obtain reliable computational predictions for reducing safety factors by robust design optimization. Classical Monte Carlo uncertainty quantification for modeling random parameters is computationally intensive compared to the more efficient Polynomial Chaos method. The global polynomial approximation employed in the Polynomial Chaos formulation can, however, give unreliable results for discontinuous responses. Polynomial Chaos methods also require a fast increasing number of samples with time to maintain a constant accuracy in unsteady problems. In this paper, a robust extrema diminishing Polynomial Chaos method is presented based on Newton-Cotes quadrature in a simplex elements discretization of probability space [1, 2]. The method results in a bounded error in time-dependent problems by 1

Jeroen A.S. Witteveen and Hester Bijl

performing the uncertainty quantification interpolation of oscillatory samples efficiently at constant phase. The application to a transonic airfoil flutter problem shows that the asymptotic pitch angle standard deviation is 16 times larger than the initial condition. 2

MATHEMATICAL UNCERTAINTY QUANTIFICATION PROBLEM

Consider a dynamical system subject to na uncorrelated second-order random input parameters a(ω) = {a1 (ω), .., ana (ω)} ∈ A with parameter space A ∈ Rna , which governs an oscillatory response u(x, t, a) L(x, t, a; u(x, t, a)) = S(x, t, a),

(1)

with operator L and source term S defined on domain D × T × A. The spatial and temporal dimensions are defined as x ∈ D and t ∈ T , with D ⊂ Rd , d = {1, 2, 3}, and T = [0, tmax ]. A realization of the set of outcomes Ω of the probability space (Ω, F , P ) is denoted by ω ∈ Ω, with F ⊂ 2Ω the σ-algebra of events and P a probability measure. Here we consider a non-intrusive uncertainty quantification method l which reuses an existing deterministic solver for fluid-structure interaction simulation of (1). Nonintrusive uncertainty quantification method l is a combination of a sampling method g and an interpolation method h. Sampling method g defines the ns sampling points ak ≡ a(ωk ) and returns the deterministic samples v(x, t) = {v1 (x, t), .., vns (x, t)} with vk (x, t) ≡ u(x, t, ak ). Interpolation method h constructs an interpolation surface w(x, t, a) through the ns samples v(x, t) as a weighted approximation of u(x, t, a). 3

ROBUST EXTREMA DIMINISHING METHOD

The multi-element Polynomial Chaos method l based on Newton-Cotes quadrature points in simplex elements [1] evaluates the statistical moment integral by dividing parameter space A into ne non-overlapping simplex elements Aj ⊂ A. A piecewise polynomial approximation w(x, t, a) is then constructed based on ns deterministic solutions vj,k (x, t) = u(x, t, aj,k ) for the values of the random parameters aj,k that correspond to the n ˜ s Newton-Cotes quadrature points of degree d in the elements Aj n ˜s ne X ne Z X X cj,k vj,k (x, t)i , (2) w(x, t, a)ifa (a)da = µwi (x, t) = j=1

Aj

j=1 k=1

where cj,k are Polynomial Chaos Newton-Cotes quadrature weights. Here, second degree Newton-Cotes quadrature is considered in combination with adaptive mesh refinement in probability space, see Figure 1. It is proven in [2] that the resulting approach satisfies the extrema diminishing (ED) robustness concept in probability space min(w(a)) ≥ min(u(a)) ∧ max(w(a)) ≤ max(u(a)) ∀u(a). A

A

A

A

(3)

The ED property leads to the advantage that no non-zero probabilities of unphysical realizations can be predicted due to overshoots at discontinuities in the response surface.

2

Jeroen A.S. Witteveen and Hester Bijl

(a) Element

(b) Initial grid

(c) Adapted grid

Figure 1: Discretization of two-dimensional parameter space A using 2-simplex elements and seconddegree Newton-Cotes quadrature points given by the dots.

4

EFFICIENT INTERPOLATION AT CONSTANT PHASE

Assume that solving equation (1) for realizations of the random parameters ak results in oscillatory samples vk (t) = u(ak ), of which the phase vφk (t) = φ(t, ak ) is a well-defined monotonically increasing function of time t for k = 1, .., ns . In order to interpolate the samples v(t) = {v1 (t), .., vna (t)} at constant phase [1], they are transformed into functions ˆ (vφ (t)) according to vˆk (vφk (t)) = vk (t) for k = 1, .., ns , see Figure 2. The inof their phase v ˆ (vφ (t)) is transformed back to an approximation terpolation w(w ˆ φ (t, a), a) of the samples v in the time domain w(t, a) = w(w ˆ φ (t, a), a). This uncertainty quantification formulation for oscillatory responses is proven to achieve a bounded error εˆ(ϕ, a) = w(ϕ, ˆ a) − uˆ(ϕ, a) as function of phase ϕ for periodic responses according to εˆ(ϕ, a) < δ

∀ϕ ∈ R, a ∈ A,

(4)

∀ϕ ∈ [0, 1], a ∈ A.

(5)

where δ is defined by εˆ(ϕ, a) < δ,

The error ε(t, a) = w(t, a) − u(t, a) is also bounded in time under certain conditions, see [2]. The phases vφ (t) are extracted from the samples using a trial-and-error algorithm based on the local extrema of the time series v(t). v^ k

vk

phase φ

time t

(a) samples vk (t)

(b) samples vˆk (φ)

Figure 2: Oscillatory samples as function of time and phase.

3

Jeroen A.S. Witteveen and Hester Bijl

(a) t = 2.5

(b) Stochastic grid

Figure 3: Response surface of angle of attack α(ω) as function of random natural frequency ratio ω ¯ (ω) and free stream velocity U∞ (ω) for the transonic airfoil flutter problem.

5

STOCHASTIC TRANSONIC AIRFOIL FLUTTER

The stochastic post flutter behavior of an elastically mounted airfoil in Euler flow is analyzed at a bifurcation parameter value U ∗ of 130% of the deterministic linear bifurcation point. The response of a structural model of a pitch-plunge airfoil with cubic nonlinear spring stiffness at a Mach number of M∞ = 0.8 is considered. The randomness in the ratio of natural frequencies ω ¯ (ω) and free stream velocity U∞ (ω) is given by a uniform and beta distribution with a coefficient of variation of 10% and 1%, respectively. The response surface approximation of the angle of attack α(t, ω) as function of the random parameters ω ¯ (ω) and U∞ (ω) given in Figure 3a shows a highly oscillatory response surface at t = 2.5. The resulting asymptotic standard deviation of σα = 1.6o is a factor 16 larger than the initial angle of attack α(0) = 0.1o . This result is obtained using the time-independent grid with ns = 9 samples and ne = 2 elements in probability space shown in Figure 3b. 6

CONCLUSIONS

The presented extrema diminishing uncertainty quantification method with bounded error predicts a 16 times larger asymptotic standard deviation compared to the initial condition due to the effect of physical uncertainties in a transonic airfoil flutter problem. REFERENCES [1] J.A.S. Witteveen and H. Bijl. An alternative unsteady adaptive stochastic finite elements formulation based on interpolation at constant phase. Comput. Method Appl. M. 198, 578–591, (2008). [2] J.A.S. Witteveen and H. Bijl. A TVD uncertainty quantification method with bounded error applied to transonic airfoil flutter. Commun. Comput. Phys. 6, 406– 432, (2009).

4

a robust and efficient uncertainty quantification method ...

∗e-mail: [email protected], web page: http://www.jeroenwitteveen.com. †e-mail: ... Numerical errors in multi-physics simulations start to reach acceptable ...

348KB Sizes 0 Downloads 318 Views

Recommend Documents

efficient uncertainty quantification in unsteady ...
of the probability distribution and statistical moments µui (x,t) of the output u(x, t, a), .... subject to a symmetric unimodal beta distribution with β1 = β2 = 2 with a ...

Efficient Uncertainty Quantification in Computational Fluid-Structure ...
Sep 21, 2007 - Abstract. Uncertainty quantification in complex flow and fluid-structure interaction simulations requires efficient uncertainty quantification meth-.

Uncertainty Quantification for Stochastic Subspace ...
Uncertainty Quantification for Stochastic Subspace Identification on. Multi-Setup Measurements. Michael Döhler, Xuan-Binh Lam and Laurent Mevel. Abstract— ...

Uncertainty Quantification for Laminar-Turbulent ... - Jeroen Witteveen
Jan 7, 2011 - 8. Considerable effort has been spent in the past two decades to develop transition models for engineering applications to predict transitional ...

Uncertainty Quantification for Laminar-Turbulent ...
entirely on empirical correlations obtained from existing data sets for simple flow configurations. ... Postdoctoral Fellow, Center for Turbulence Research, Building 500, Stanford University, Stanford, CA ... 4 - 7 January 2011, Orlando, Florida.

Uncertainty Quantification for Stochastic Damage ...
finite element model of the structure. Damage localization using both finite element information and modal parameters estimated from ambient vibration data collected from sensors is possible by the Stochastic Dynamic Damage Location Vector (SDDLV) ap

Uncertainty Quantification for Multi-Frequency ...
reliable predictions, which can be utilized in robust design optimization and reducing .... on a 60 × 15 × 30m domain using an unstructured hexahedral mesh. ... Results for the time evolution of the mean µL(t) and the standard deviation σL(t) of 

Uncertainty quantification and error estimation in ... - Semantic Scholar
non-reactive simulations, Annual Research Briefs, Center for Turbulence Research, Stanford University (2010) 57–68. 9M. Smart, N. Hass, A. Paull, Flight data ...

Uncertainty quantification and error estimation in scramjet simulation
is overwhelmed by abundant uncertainties and errors. This limited accuracy of the numerical prediction of in-flight performance seriously hampers scramjet design. The objective of the Predictive Science Academic Alliance Program (PSAAP) at Stanford U

Uncertainty quantification and error estimation in scramjet simulation
to improve the current predictive simulation capabilities for scramjet engine flows. ..... solution is obtained by the automatic differentiation software package ...

Uncertainty quantification and error estimation in scramjet simulation
V.A. Aleatoric flight conditions. The flow conditions of the HyShot II flight are uncertain due to the failure of the radar tracking system during the experiment. Therefore, the free stream conditions were inferred from pressure measurements in the u

Enabling Robust and Efficient Distributed ...
relatively recent P2P-based storage services that allow data to be stored and retrieved among peers [3]. ... recently, for cloud computing services as well [2], [18]. ...... [45] R. O'Dell and R. Wattenhofer, “Information dissemination in highly ..

A Simple and Efficient Sampling Method for Estimating ...
Jul 24, 2008 - Storage and Retrieval; H.3.4 Systems and Software: Perfor- mance Evaluation ...... In TREC Video Retrieval Evaluation Online. Proceedings ...

Quantification of uncertainty in nonlinear soil models at ...
Recent earth- quakes in Japan ... al Research Institute for Earth Science and Disaster. Prevention ..... not predicting a large degree of nonlinear soil be- havior.

Uncertainty Quantification in MD simulations of ...
that, based on the target application, only certain molecules or ions can pass ...... accordance with Figure 8 (b) which showed that the noise level in the Cl− ..... Figure 12 shows that the predictive samples form a “cloud” demonstrating that

Quantification of uncertainty in nonlinear soil models at ...
Jun 17, 2013 - DEEPSOIL. – ABAQUS. • Equivalent-linear models: Within SHAKE, the following modulus-reduction and damping relationships are tested: – Zhang et al. (2005). – Darendeli (2001). • Nonlinear models: - DEEPSOIL (Hashash et al., 20

Uncertainty quantification in MD simulations of ...
the accuracy with which the potential function can reproduce the atomic .... σ) and Coulombic charge, in multiples of the electron charge |e|, for each atom type.

Uncertainty Quantification in Fluid-Structure Interaction ...
the probability distribution and statistical moments µui (x,t) of the output u(x, t, a), which ...... beta distribution with a coefficient of variation of 1% around a mean of ...

Uncertainty Quantification for the Trailing-Edge Noise of ...
on the restricted domain with the above extracted velocity profiles, directly yields .... from LWT RANS computations and (dash-dot) uncertainty bounds around inlet ..... G., Wang, M. & Roger, M. 2003 Analysis of flow conditions in free-jet experi-.

MORE EFFICIENT TESTS ROBUST TO ...
MacKinnon and White (1985) considered a number of possible forms of the estimator. They showed that the HCCME's, ... estimator appears to be serious and disappears very slowly as the number of observations increases. ...... with a few additional inst

Efficient and Robust Music Identification with Weighted Finite-State ...
be used to give a compact representation of all song snippets for a large .... song during the course of three iterations of acoustic model training. mpx stands for ...... transducers in speech recognition,” Computer Speech and Language, vol.

Efficient and Robust Music Identification with Weighted Finite-State ...
a database of over 15 000 songs while achieving an identification ... our system is robust to several different types of noise and ..... can be distributed along a path. ...... [9] D. Pye, “Content-based methods for the management of digital music,