81 research outputs found
Failure Probability Estimation and Detection of Failure Surfaces via Adaptive Sequential Decomposition of the Design Domain
We propose an algorithm for an optimal adaptive selection of points from the
design domain of input random variables that are needed for an accurate
estimation of failure probability and the determination of the boundary between
safe and failure domains. The method is particularly useful when each
evaluation of the performance function g(x) is very expensive and the function
can be characterized as either highly nonlinear, noisy, or even discrete-state
(e.g., binary). In such cases, only a limited number of calls is feasible, and
gradients of g(x) cannot be used. The input design domain is progressively
segmented by expanding and adaptively refining mesh-like lock-free geometrical
structure. The proposed triangulation-based approach effectively combines the
features of simulation and approximation methods. The algorithm performs two
independent tasks: (i) the estimation of probabilities through an ingenious
combination of deterministic cubature rules and the application of the
divergence theorem and (ii) the sequential extension of the experimental design
with new points. The sequential selection of points from the design domain for
future evaluation of g(x) is carried out through a new learning function, which
maximizes instantaneous information gain in terms of the probability
classification that corresponds to the local region. The extension may be
halted at any time, e.g., when sufficiently accurate estimations are obtained.
Due to the use of the exact geometric representation in the input domain, the
algorithm is most effective for problems of a low dimension, not exceeding
eight. The method can handle random vectors with correlated non-Gaussian
marginals. The estimation accuracy can be improved by employing a smooth
surrogate model. Finally, we define new factors of global sensitivity to
failure based on the entire failure surface weighted by the density of the
input random vector.Comment: 42 pages, 24 figure
A Dimension-Adaptive Multi-Index Monte Carlo Method Applied to a Model of a Heat Exchanger
We present an adaptive version of the Multi-Index Monte Carlo method,
introduced by Haji-Ali, Nobile and Tempone (2016), for simulating PDEs with
coefficients that are random fields. A classical technique for sampling from
these random fields is the Karhunen-Lo\`eve expansion. Our adaptive algorithm
is based on the adaptive algorithm used in sparse grid cubature as introduced
by Gerstner and Griebel (2003), and automatically chooses the number of terms
needed in this expansion, as well as the required spatial discretizations of
the PDE model. We apply the method to a simplified model of a heat exchanger
with random insulator material, where the stochastic characteristics are
modeled as a lognormal random field, and we show consistent computational
savings
Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds
In this paper, we study the efficient numerical integration of functions with
sharp gradients and cusps. An adaptive integration algorithm is presented that
systematically improves the accuracy of the integration of a set of functions.
The algorithm is based on a divide and conquer strategy and is independent of
the location of the sharp gradient or cusp. The error analysis reveals that for
a function (derivative-discontinuity at a point), a rate of convergence
of is obtained in . Two applications of the adaptive integration
scheme are studied. First, we use the adaptive quadratures for the integration
of the regularized Heaviside function---a strongly localized function that is
used for modeling sharp gradients. Then, the adaptive quadratures are employed
in the enriched finite element solution of the all-electron Coulomb problem in
crystalline diamond. The source term and enrichment functions of this problem
have sharp gradients and cusps at the nuclei. We show that the optimal rate of
convergence is obtained with only a marginal increase in the number of
integration points with respect to the pure finite element solution with the
same number of elements. The adaptive integration scheme is simple, robust, and
directly applicable to any generalized finite element method employing
enrichments with sharp local variations or cusps in -dimensional
parallelepiped elements.Comment: 22 page
Edge Detection by Adaptive Splitting II. The Three-Dimensional Case
In Llanas and Lantarón, J. Sci. Comput. 46, 485–518 (2011) we proposed an algorithm (EDAS-d) to approximate the jump discontinuity set of functions defined on subsets of ℝ d . This procedure is based on adaptive splitting of the domain of the function guided by the value of an average integral. The above study was limited to the 1D and 2D versions of the algorithm. In this paper we address the three-dimensional problem. We prove an integral inequality (in the case d=3) which constitutes the basis of EDAS-3. We have performed detailed computational experiments demonstrating effective edge detection in 3D function models with different interface topologies. EDAS-1 and EDAS-2 appealing properties are extensible to the 3D cas
A Node Elimination Algorithm for Cubature of High-Dimensional Polytopes
Node elimination is a numerical approach to obtain cubature rules for the
approximation of multivariate integrals. Beginning with a known cubature rule,
nodes are selected for elimination, and a new, more efficient rule is
constructed by iteratively solving the moment equations. This paper introduces
a new criterion for selecting which nodes to eliminate that is based on a
linearization of the moment equation. In addition, a penalized iterative solver
is introduced, that ensures that weights are positive and nodes are inside the
integration domain. A strategy for constructing an initial quadrature rule for
various polytopes in several space dimensions is described. High efficiency
rules are presented for two, three and four dimensional polytopes. The new
rules are compared with rules that are obtained by combining tensor products of
one dimensional quadrature rules and domain transformations, as well as with
known analytically constructed cubature rules.Comment: 18 pages, 6 figure
Numerical Integration in S-PLUS or R: A Survey
This paper reviews current quadrature methods for approximate calculation of integrals within S-Plus or R. Starting with the general framework, Gaussian quadrature will be discussed first, followed by adaptive rules and Monte Carlo methods. Finally, a comparison of the methods presented is given. The aim of this survey paper is to help readers, not expert in computing, to apply numerical integration methods and to realize that numerical analysis is an art, not a science.
Scattering of general incident beams by diffraction gratings
The paper is devoted to the electromagnetic scattering of arbitrary time-harmonic fields by periodic structures. The Floquet-Fourier transform converts the full space Maxwell problem to a two-parameter family of diffraction problems with quasiperiodic incidence waves, for which conventional grating methods become applicable. The inverse transform is given by integrating with respect to the parameters over a infinite strip in ℝ². For the computation of the scattered fields we propose an algorithm, which extends known adaptive methods for the approximate calculation of multiple integrals. The novel adaptive approach provides autonomously the expansion of the incident field into quasiperiodic waves in order to approximate the scattered fields within a prescribed error tolerance. Some application examples are numerically examined
- …