Source author record

Viet Ha Hoang

Viet Ha Hoang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.AP math.NA Numerical Analysis math.PR math.OC math.ST Methodology Statistics Theory

Catalog footprint

What is connected

9works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

Bayesian Inversion of Log-normal Eikonal Equations

We study the Bayesian inverse problem for inferring the log-normal slowness function of the eikonal equation given noisy observation data on its solution at a set of spatial points. We study approximation of the posterior probability measure by solving the truncated eikonal equation, which contains only a finite number of terms in the Karhunen-Loeve expansion of the slowness function, by the Fast Marching Method. The error of this approximation in the Hellinger metric is deduced in terms of the truncation level of the slowness and the grid size in the Fast Marching Method resolution. It is well known that the plain Markov Chain Monte Carlo procedure for sampling the posterior probability is highly expensive. We develop and justify the convergence of a Multilevel Markov Chain Monte Carlo method. Using the heap sort procedure in solving the forward eikonal equation by the Fast Marching Method, our Multilevel Markov Chain Monte Carlo method achieves a prescribed level of accuracy for approximating the posterior expectation of quantities of interest, requiring only an essentially optimal level of complexity. Numerical examples confirm the theoretical results.

preprint2022arXiv

Optimal finite elements for ergodic stochastic two-scale elliptic equations

We develop an essentially optimal finite element approach for solving ergodic stochastic two-scale elliptic equations whose two-scale coefficient may depend also on the slow variable. We solve the limiting stochastic two-scale homogenized equation obtained from the stochastic two-scale convergence in the mean (A. Bourgeat, A. Mikelic and S. Wright, J. reine angew. Math, Vol. 456, 1994), whose solution comprises of the solution to the homogenized equation and the corrector, by truncating the infinite domain of the fast variable and using the sparse tensor product finite elements. We show that the convergence rate in terms of the truncation level is equivalent to that for solving the cell problems in the same truncated domain. Solving this equation, we obtain the solution to the homogenized equation and the corrector at the same time, using only a number of degrees of freedom that is essentially equivalent to that required for solving one cell problem. Optimal complexity is obtained when the corrector possesses sufficient regularity with respect to both the fast and the slow variables. Although the regularity norm of the corrector depends on the size of the truncated domain, we show that the convergence rate of the approximation for the solution to the homogenized equation is independent of the size of the truncated domain. With the availability of an analytic corrector, we construct a numerical corrector for the solution of the original stochastic two-scale equation from the finite element solution to the truncated stochastic two-scale homogenized equation. Numerical examples of quasi-periodic two-scale equations, and a stochastic two-scale equation of the checker board type, whose coefficient is discontinuous, confirm the theoretical results.

preprint2020arXiv

Reconciling Bayesian and perimeter regularization for binary inversion

A central theme in classical algorithms for the reconstruction of discontinuous functions from observational data is perimeter regularization via the use of the total variation. On the other hand, sparse or noisy data often demands a probabilistic approach to the reconstruction of images, to enable uncertainty quantification; the Bayesian approach to inversion, which itself introduces a form of regularization, is a natural framework in which to carry this out. In this paper the link between Bayesian inversion methods and perimeter regularization is explored. In this paper two links are studied: (i) the maximum a posteriori (MAP) objective function of a suitably chosen Bayesian phase-field approach is shown to be closely related to a least squares plus perimeter regularization objective; (ii) sample paths of a suitably chosen Bayesian level set formulation are shown to possess finite perimeter and to have the ability to learn about the true perimeter.

preprint2019arXiv

Multiscale simulations for upscaled multi-continuum flows

We consider in this paper a challenging problem of simulating fluid flows, in complex multiscale media possessing multi-continuum background. As an effort to handle this obstacle, model reduction is employed. In \cite{rh2}, homogenization was nicely applied, to find effective coefficients and homogenized equations (for fluid flow pressures) of a dual-continuum system, with new convection terms and negative interaction coefficients. However, some degree of multiscale still remains. This motivates us to propose the generalized multiscale finite element method (GMsFEM), which is coupled with the dual-continuum homogenized equations, toward speeding up the simulation, improving the accuracy as well as clearly representing the interactions between the dual continua. In our paper, globally, each continuum is viewed as a system and connected to the other throughout the domain. We take into consideration the flow transfers between the dual continua and within each continuum itself. Such multiscale flow dynamics are modeled by the GMsFEM, which systematically generates either uncoupled or coupled multiscale basis (to carry the local characteristics to the global ones), via establishing local snapshots and spectral decomposition in the snapshot space. As a result, we will work with a system of two equations coupled with some interaction terms, and each equation describes one of the dual continua on the fine grid. Convergence analysis of the proposed GMsFEM is accompanied with the numerical results, which support the favorable outcomes.

preprint2015arXiv

Homogenization error for two scale Maxwell equations

For two scale elliptic equations in a domain $D$, standard homogenization errors are deduced with the assumption that the solution $u_0$ of the homogenized equation belongs to $H^2(D)$. For two scale Maxwell equations, the corresponding required regularity is $u_0\in H^1({\rm curl}, D)$. These regularity conditions normally do not hold in general polygonal domains, which are of interests for finite element discretization. The paper establishes homogenization errors when $u_0$ belongs to a weaker regularity space $H^{1+s}(D)$ for elliptic problems and $H^s({\rm curl},D)$ for Maxwell problems where $0<s<1$. Though we only present the results for two scale Maxwell equations when $u_0\in H^s({\rm curl},D)$ with $0<s<1$, the procedure works verbatim for elliptic equations when $u_0$ belongs to $H^{1+s}(D)$ with $0<s<1$.

preprint2015arXiv

Polynomial approximations of a class of stochastic multiscale elasticity problems

We consider a class of elasticity equations in ${\mathbb R}^d$ whose elastic moduli depend on $n$ separated microscopic scales, are random and expressed as a linear expansion of a countable sequence of random variables which are independently and identically uniformly distributed in a compact interval. The multiscale Hellinger-Reissner problem that allows for computing the stress directly, and the multiscale mixed problem for nearly incompressible isotropic materials are considered. The stochastic problems are studied via deterministic problems that depend on a countable number of real parameters. We study the multiscale homogenized problems that contain all the macroscopic and microscopic information, whose solutions are written as generalized polynomial chaos (gpc) expansions. We approximate these solutions by semidiscrete Galerkin approximating problems that project into the spaces of functions with only a finite number of $N$ gpc modes. We deduce bounds and summability properties for the solutions' gpc expansion coefficients, which imply explicit rates of convergence in terms of $N$ when the gpc modes used for the Galerkin approximation are chosen to correspond to the best $N$ terms in the gpc expansion. For the mixed problem for nearly incompressible materials, the rate of convergence for the best $N$ term approximation is independent of the Lamé constants' ratio. We establish parametric correctors in terms of the semidiscrete Galerkin approximations. For two scale problems, an explicit homogenization rate is deduced. Together with the best $N$ term rate, it provides an explicit convergence rate for the correctors of the parametric multiscale problems. For nearly incompressible materials, we obtain a homogenization rate that is independent of the ratio of the Lamé constants, so that the error for the corrector is also independent of this ratio.

preprint2014arXiv

Determining White Noise Forcing From Eulerian Observations in the Navier Stokes Equation

The Bayesian approach to inverse problems is of paramount importance in quantifying uncertainty about the input to and the state of a system of interest given noisy observations. Herein we consider the forward problem of the forced 2D Navier Stokes equation. The inverse problem is inference of the forcing, and possibly the initial condition, given noisy observations of the velocity field. We place a prior on the forcing which is in the form of a spatially correlated temporally white Gaussian process, and formulate the inverse problem for the posterior distribution. Given appropriate spatial regularity conditions, we show that the solution is a continuous function of the forcing. Hence, for appropriately chosen spatial regularity in the prior, the posterior distribution on the forcing is absolutely continuous with respect to the prior and is hence well-defined. Furthermore, the posterior distribution is a continuous function of the data. We complement this theoretical result with numerical simulation of the posterior distribution.

preprint2013arXiv

Complexity Analysis of Accelerated MCMC Methods for Bayesian Inversion

We study Bayesian inversion for a model elliptic PDE with unknown diffusion coefficient. We provide complexity analyses of several Markov Chain-Monte Carlo (MCMC) methods for the efficient numerical evaluation of expectations under the Bayesian posterior distribution, given data $δ$. Particular attention is given to bounds on the overall work required to achieve a prescribed error level $\varepsilon$. Specifically, we first bound the computational complexity of "plain" MCMC, based on combining MCMC sampling with linear complexity multilevel solvers for elliptic PDE. Our (new) work versus accuracy bounds show that the complexity of this approach can be quite prohibitive. Two strategies for reducing the computational complexity are then proposed and analyzed: first, a sparse, parametric and deterministic generalized polynomial chaos (gpc) "surrogate" representation of the forward response map of the PDE over the entire parameter space, and, second, a novel Multi-Level Markov Chain Monte Carlo (MLMCMC) strategy which utilizes sampling from a multilevel discretization of the posterior and of the forward PDE. For both of these strategies we derive asymptotic bounds on work versus accuracy, and hence asymptotic bounds on the computational complexity of the algorithms. In particular we provide sufficient conditions on the regularity of the unknown coefficients of the PDE, and on the approximation methods used, in order for the accelerations of MCMC resulting from these strategies to lead to complexity reductions over "plain" MCMC algorithms for Bayesian inversion of PDEs.}

preprint2011arXiv

Bayesian inverse problems for Burgers and Hamilton-Jacobi equations with white-noise forcing

The paper formulates Bayesian inverse problems for inference in a topological measure space given noisy observations. Conditions for the validity of the Bayes formula and the well-posedness of the posterior measure are studied. The abstract theory is then applied to Burgers and Hamilton-Jacobi equations on a semi-infinite time interval with forcing functions which are white noise in time. Inference is made on the white noise forcing, assuming the Wiener measure as the prior.