This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large‐scale sampling‐based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right‐hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. We demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallel forward model problem solver, for large‐scale 3D MLMC simulations with up to 1.9·10^{9}unknowns.
 Award ID(s):
 2110773
 NSFPAR ID:
 10340868
 Date Published:
 Journal Name:
 HOUSEHOLDER SYMPOSIUM XXI
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Summary 
Abstract The reduction of a large‐scale symmetric linear discrete ill‐posed problem with multiple right‐hand sides to a smaller problem with a symmetric block tridiagonal matrix can easily be carried out by the application of a small number of steps of the symmetric block Lanczos method. We show that the subdiagonal blocks of the reduced problem converge to zero fairly rapidly with increasing block number. This quick convergence indicates that there is little advantage in expressing the solutions of discrete ill‐posed problems in terms of eigenvectors of the coefficient matrix when compared with using a basis of block Lanczos vectors, which are simpler and cheaper to compute. Similarly, for nonsymmetric linear discrete ill‐posed problems with multiple right‐hand sides, we show that the solution subspace defined by a few steps of the block Golub–Kahan bidiagonalization method usually can be applied instead of the solution subspace determined by the singular value decomposition of the coefficient matrix without significant, if any, reduction of the quality of the computed solution.

We combine datadriven reduced order models (ROM) with the Lippmann Schwinger integral equation to produce a direct nonlinear inversion method. The ROM is viewed as a Galerkin projection and is sparse due to Lanczos orthogonalization. Embedding into the continuous problem, a datadriven internal solution is produced. This internal solution is then used in the LippmannSchwinger equation, in a direct or iterative framework. The new approach also allows us to process nonsquare matrixvalued datatransfer functions, i.e., to remove the main limitation of the earlier versions of the ROM based inversion algorithms. We show numerical experiments for spectral domain data for which our inversion is far superior to the Born inversion.more » « less

null (Ed.)The 74 ka Youngest Toba Tuff (YTT) was discovered as cryptotephra in South African archaeological sites at Pinnacle Point (PP) 56N, Vleesbaai [1] and Klasies River on the Indian Ocean, and the Diepkloof Rock Shelter on the Atlantic coast nearly 750 km west of PP. The YTT eruption distributed tephra across eastern and southern Africa and provides a widespread isochron useful for dating archaeological deposits, testing age models, and precisely determining the timing of changes in human behavior. At PP, we demonstrated that the MIS 45 transition began just before the YTT eruption. Humans thrived both through the YTT event and the changing climate, and important changes in technology occurred just after the Toba eruption [1]. Controversy related to trapped charge age models at the Diepkloof rock shelter [2,3] were resolved by identifying YTT at a location in the stratigraphic section that confirmed the Jacobs et al. [2] model for the site, confirming that technological changes similar to those observed at PP occurred synchronously to those at Diepkloof, not substantially before as suggested by a prior published age model [3]. Processing samples with very low abundance cryptotephra, such as that found in South Africa, is a challenge and requires revision of standard laboratory techniques. Samples with high organic or clay content benefit by being treated with 10% HCl and 3% hydrogen peroxide (H2O2). This step does not degrade shard integrity or affect chemistry and is useful in separating shards from clay and organic particles allowing better recovery in heavy liquids. We also modified the heavy liquid density range from 1.95  2.55 to 2.2  2.5 g/cm3 to more effectively remove quartz and feldspar as well as biogenic silica. This density range captures shards ranging from rhyolite to dacite in composition. [1] Smith, E. I. et al. Nature 555, 511, 2018 [2] Jacobs, Z. & Roberts, R. G. Journal of Archaeological Science 63, 175192, 2015. [3] Tribolo, C. et al. Journal of Archaeological Science 40, 34013411 2013.more » « less

One of classical tasks of the network synthesis is to construct ROMs realized via ladder networks matching rational approximations of a targeted filter transfer function. The inverse scattering can be also viewed in the network synthesis framework. The key is continuum interpretation of the synthesized network in terms of the underlying medium properties, aka embedding. We describe such an embedding via finitedifference quadrature rules (FDQR), that can be viewed as extension of the concept of the Gaussian quadrature to finitedifference schemes. One of application of this approach is the solution of earlier intractable large scale inverse scattering problems. We also discuss an important open question in the FDQR related to Lothar’s earlier contributions, in particular, a possibility of finitedifference GaussKronrod rulesmore » « less