Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAdaptive stratified sampling for non-smooth problems
72
0
0.0
(
0
)
Added by
Per Pettersson
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Science and engineering problems subject to uncertainty are frequently both computationally expensive and feature nonsmooth parameter dependence, making standard Monte Carlo too slow, and excluding efficient use of accelerated uncertainty quantification methods relying on strict smoothness assumptions. To remedy these challenges, we propose an adaptive stratification method suitable for nonsmooth problems and with significantly reduced variance compared to Monte Carlo sampling. The stratification is iteratively refined and samples are added sequentially to satisfy an allocation criterion combining the benefits of proportional and optimal sampling. Theoretical estimates are provided for the expected performance and probability of failure to correctly estimate essential statistics. We devise a practical adaptive stratification method with strata of the same kind of geometrical shapes, cost-effective refinement satisfying a greedy variance reduction criterion. Numerical experiments corroborate the theoretical findings and exhibit speedups of up to three orders of magnitude compared to standard Monte Carlo sampling.
rate research
Read More
We propose an efficient numerical strategy for simulating fluid flow through porous media with highly oscillatory characteristics. Specifically, we consider non-linear diffusion models. This scheme is based on the classical homogenization theory and uses a locally mass-conservative formulation. In addition, we discuss some properties of the standard non-linear solvers and use an error estimator to perform a local mesh refinement. The main idea is to compute the effective parameters in such a way that the computational complexity is reduced without affecting the accuracy. We perform some numerical examples to illustrate the behaviour of the adaptive scheme and of the non-linear solvers. Finally, we discuss the advantages of the implementation of the numerical homogenization in a periodic media and the applicability of the same scheme in non-periodic test cases such as SPE10th project.
The analysis of linear ill-posed problems often is carried out in function spaces using tools from functional analysis. However, the numerical solution of these problems typically is computed by first discretizing the problem and then applying tools from (finite-dimensional) linear algebra. The present paper explores the feasibility of applying the Chebfun package to solve ill-posed problems. This approach allows a user to work with functions instead of matrices. The solution process therefore is much closer to the analysis of ill-posed problems than standard linear algebra-based solution methods.
In this paper, we study temporal splitting algorithms for multiscale problems. The exact fine-grid spatial problems typically require some reduction in degrees of freedom. Multiscale algorithms are designed to represent the fine-scale details on a coarse grid and, thus, reduce the problems size. When solving time-dependent problems, one can take advantage of the multiscale decomposition of the solution and perform temporal splitting by solving smaller-dimensional problems, which is studied in the paper. In the proposed approach, we consider the temporal splitting based on various low dimensional spatial approximations. Because a multiscale spatial splitting gives a good decomposition of the solution space, one can achieve an efficient implicit-explicit temporal discretization. We present a recently developed theoretical result in our earlier work and adopt it in this paper for multiscale problems. Numerical results are presented to demonstrate the efficiency of the proposed splitting algorithm.
We consider the inverse source problems with multi-frequency sparse near field measurements. In contrast to the existing near field operator based on the integral over the space variable, a multi-frequency near field operator is introduced based on the integral over the frequency variable. A factorization of this multi-frequency near field operator is further given and analysed. Motivated by such a factorization, we introduce a multi-frequency sampling method to reconstruct the source support. Its theoretical foundation is then derived from the properties of the factorized operators and a properly chosen point spread function. Numerical examples are provided to illustrate the multi-frequency sampling method with sparse near field measurements. Finally we briefly discuss how to extend the near field case to the far field case.
Consider the electromagnetic scattering of a time-harmonic plane wave by an open cavity which is embedded in a perfectly electrically conducting infinite ground plane. This paper is concerned with the numerical solutions of the transverse electric and magnetic polarizations of the open cavity scattering problems. In each polarization, the scattering problem is reduced equivalently into a boundary value problem of the two-dimensional Helmholtz equation in a bounded domain by using the transparent boundary condition (TBC). An a posteriori estimate based adaptive finite element method with the perfectly matched layer (PML) technique is developed to solve the reduced problem. The estimate takes account both of the finite element approximation error and the PML truncation error, where the latter is shown to decay exponentially with respect to the PML medium parameter and the thickness of the PML layer. Numerical experiments are presented and compared with the adaptive finite element TBC method for both polarizations to illustrate the competitive behavior of the proposed method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
