-
Nonparametric density estimation for stationary processes under multiplicative measurement errors
Authors:
Duc Trong Dang,
Van Ha Hoang,
Phuc Hung Thai
Abstract:
This paper focuses on estimating the invariant density function $f_X$ of the strongly mixing stationary process $X_t$ in the multiplicative measurement errors model $Y_t = X_t U_t$, where $U_t$ is also a strongly mixing stationary process. We propose a novel approach to handle non-independent data, typical in real-world scenarios. For instance, data collected from various groups may exhibit interd…
▽ More
This paper focuses on estimating the invariant density function $f_X$ of the strongly mixing stationary process $X_t$ in the multiplicative measurement errors model $Y_t = X_t U_t$, where $U_t$ is also a strongly mixing stationary process. We propose a novel approach to handle non-independent data, typical in real-world scenarios. For instance, data collected from various groups may exhibit interdependencies within each group, resembling data generated from $m$-dependent stationary processes, a subset of stationary processes. This study extends the applicability of the model $Y_t = X_t U_t$ to diverse scientific domains dealing with complex dependent data. The paper outlines our estimation techniques, discusses convergence rates, establishes a lower bound on the minimax risk, and demonstrates the asymptotic normality of the estimator for $f_X$ under smooth error distributions. Through examples and simulations, we showcase the efficacy of our estimator. The paper concludes by providing proofs for the presented theoretical results.v
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
Bayesian inversion for Electrical Impedance Tomography by sparse interpolation
Authors:
Quang Huy Pham,
Viet Ha Hoang
Abstract:
We study the Electrical Impedance Tomography Bayesian inverse problem for recovering the conductivity given noisy measurements of the voltage on some boundary surface electrodes. The uncertain conductivity depends linearly on a countable number of uniformly distributed random parameters in a compact interval, with the coefficient functions in the linear expansion decaying at an algebraic rate. We…
▽ More
We study the Electrical Impedance Tomography Bayesian inverse problem for recovering the conductivity given noisy measurements of the voltage on some boundary surface electrodes. The uncertain conductivity depends linearly on a countable number of uniformly distributed random parameters in a compact interval, with the coefficient functions in the linear expansion decaying at an algebraic rate. We analyze the surrogate Markov Chain Monte Carlo (MCMC) approach for sampling the posterior probability measure, where the multivariate sparse adaptive interpolation, with interpolating points chosen according to a lower index set, is used for approximating the forward map. The forward equation is approximated once before running the MCMC for all the realizations, using interpolation on the finite element (FE) approximation at the parametric interpolating points. When evaluation of the solution is needed for a realization, we only need to compute a polynomial, thus cutting drastically the computation time. We contribute a rigorous error estimate for the MCMC convergence. In particular, we show that there is a nested sequence of interpolating lower index sets for which we can derive an interpolation error estimate in terms of the cardinality of these sets, uniformly for all the parameter realizations. An explicit convergence rate for the MCMC sampling of the posterior expectation of the conductivity is rigorously derived, in terms of the interpolating point number, the accuracy of the FE approximation of the forward equation, and the MCMC sample number. We perform numerical experiments using an adaptive greedy approach to construct the sets of interpolation points. We show the benefits of this approach over the simple MCMC where the forward equation is repeatedly solved for all the samples and the non-adaptive surrogate MCMC with an isotropic index set treating all the random parameters equally.
△ Less
Submitted 14 June, 2023;
originally announced June 2023.
-
Multilevel Markov Chain Monte Carlo for Bayesian Elliptic Inverse Problems with Besov Random Tree Priors
Authors:
Andreas Stein,
Viet Ha Hoang
Abstract:
We propose a multilevel Monte Carlo-FEM algorithm to solve elliptic Bayesian inverse problems with "Besov random tree prior". These priors are given by a wavelet series with stochastic coefficients, and certain terms in the expansion vanishing at random, according to the law of so-called Galton-Watson trees. This allows to incorporate random fractal structures and large deviations in the log-diffu…
▽ More
We propose a multilevel Monte Carlo-FEM algorithm to solve elliptic Bayesian inverse problems with "Besov random tree prior". These priors are given by a wavelet series with stochastic coefficients, and certain terms in the expansion vanishing at random, according to the law of so-called Galton-Watson trees. This allows to incorporate random fractal structures and large deviations in the log-diffusion, which occur naturally in many applications from geophysics or medical imaging. This framework entails two main difficulties: First, the associated diffusion coefficient does not satisfy a uniform ellipticity condition, which leads to non-integrable terms and thus divergence of standard multilevel estimators. Secondly, the associated space of parameters is Polish, but not a normed linear space. We address the first point by introducing cut-off functions in the estimator to compensate for the non-integrable terms, while the second issue is resolved by employing an independence Metropolis-Hastings sampler. The resulting algorithm converges in the mean-square sense with essentially optimal asymptotic complexity, and dimension-independent acceptance probabilities.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
Harmonic functions for singular quadrant walks
Authors:
Viet Hung Hoang,
Kilian Raschel,
Pierre Tarrago
Abstract:
We consider discrete (time and space) random walks confined to the quarter plane, with jumps only in directions $(i,j)$ with $i+j \geq 0$ and small negative jumps, i.e., $i,j \geq -1$. These walks are called singular, and were recently intensively studied from a combinatorial point of view. In this paper, we show how the compensation approach introduced in the 90ies by Adan, Wessels and Zijm may b…
▽ More
We consider discrete (time and space) random walks confined to the quarter plane, with jumps only in directions $(i,j)$ with $i+j \geq 0$ and small negative jumps, i.e., $i,j \geq -1$. These walks are called singular, and were recently intensively studied from a combinatorial point of view. In this paper, we show how the compensation approach introduced in the 90ies by Adan, Wessels and Zijm may be applied to compute positive harmonic functions with Dirichlet boundary conditions. In particular, in case the random walks have a drift with positive coordinates, we derive an explicit formula for the escape probability, which is the probability to tend to infinity without reaching the boundary axes. These formulas typically involve famous recurrent sequences, such as the Fibonacci numbers. As a second step, we propose a probabilistic interpretation of the previously constructed harmonic functions and prove that they allow to compute all positive harmonic functions of these singular walks. To that purpose, we derive the asymptotics of the Green functions in all directions of the quarter plane and use Martin boundary theory.
△ Less
Submitted 31 July, 2022;
originally announced August 2022.
-
Fast convergence of inertial dynamics with Hessian-driven damping under geometry assumptions
Authors:
Jean-François Aujol,
Charles Dossal,
Văn Hào Hoàng,
Hippolyte Labarrière,
Aude Rondepierre
Abstract:
First-order optimization algorithms can be considered as a discretization of ordinary differential equations (ODEs) \cite{su2014differential}. In this perspective, studying the properties of the corresponding trajectories may lead to convergence results which can be transfered to the numerical scheme. In this paper we analyse the following ODE introduced by Attouch et al. in \cite{attouch2016fast}…
▽ More
First-order optimization algorithms can be considered as a discretization of ordinary differential equations (ODEs) \cite{su2014differential}. In this perspective, studying the properties of the corresponding trajectories may lead to convergence results which can be transfered to the numerical scheme. In this paper we analyse the following ODE introduced by Attouch et al. in \cite{attouch2016fast}: \begin{equation*} \forall t\geqslant t_0,~\ddot{x}(t)+\fracα{t}\dot{x}(t)+βH_F(x(t))\dot{x}(t)+\nabla F(x(t))=0,\end{equation*} where $α>0$, $β>0$ and $H_F$ denotes the Hessian of $F$. This ODE can be derived to build numerical schemes which do not require $F$ to be twice differentiable as shown in \cite{attouch2020first,attouch2021convergence}. We provide strong convergence results on the error $F(x(t))-F^*$ and integrability properties on $\|\nabla F(x(t))\|$ under some geometry assumptions on $F$ such as quadratic growth around the set of minimizers. In particular, we show that the decay rate of the error for a strongly convex function is $O(t^{-α-\varepsilon})$ for any $\varepsilon>0$. These results are briefly illustrated at the end of the paper.
△ Less
Submitted 20 June, 2022; v1 submitted 14 June, 2022;
originally announced June 2022.
-
Functional Mixtures-of-Experts
Authors:
Faïcel Chamroukhi,
Nhat Thien Pham,
Van Hà Hoang,
Geoffrey J. McLachlan
Abstract:
We consider the statistical analysis of heterogeneous data for prediction in situations where the observations include functions, typically time series. We extend the modeling with Mixtures-of-Experts (ME), as a framework of choice in modeling heterogeneity in data for prediction with vectorial observations, to this functional data analysis context. We first present a new family of ME models, name…
▽ More
We consider the statistical analysis of heterogeneous data for prediction in situations where the observations include functions, typically time series. We extend the modeling with Mixtures-of-Experts (ME), as a framework of choice in modeling heterogeneity in data for prediction with vectorial observations, to this functional data analysis context. We first present a new family of ME models, named functional ME (FME) in which the predictors are potentially noisy observations, from entire functions. Furthermore, the data generating process of the predictor and the real response, is governed by a hidden discrete variable representing an unknown partition. Second, by imposing sparsity on derivatives of the underlying functional parameters via Lasso-like regularizations, we provide sparse and interpretable functional representations of the FME models called iFME. We develop dedicated expectation--maximization algorithms for Lasso-like (EM-Lasso) regularized maximum-likelihood parameter estimation strategies to fit the models. The proposed models and algorithms are studied in simulated scenarios and in applications to two real data sets, and the obtained results demonstrate their performance in accurately capturing complex nonlinear relationships and in clustering the heterogeneous regression data.
△ Less
Submitted 20 December, 2023; v1 submitted 4 February, 2022;
originally announced February 2022.
-
Discrete harmonic functions for non-symmetric Laplace operators in the quarter plane
Authors:
Viet Hung Hoang
Abstract:
We construct harmonic functions in the quarter plane for discrete Laplace operators. In particular, the functions are conditioned to vanish on the boundary and the Laplacians admit coefficients associated with transition probabilities of non-symmetric random walks. By solving a boundary value problem for generating functions of harmonic functions, we deduce explicit expressions for the generating…
▽ More
We construct harmonic functions in the quarter plane for discrete Laplace operators. In particular, the functions are conditioned to vanish on the boundary and the Laplacians admit coefficients associated with transition probabilities of non-symmetric random walks. By solving a boundary value problem for generating functions of harmonic functions, we deduce explicit expressions for the generating functions in terms of conformal mappings. These mappings are yielded from a conformal welding problem with quasisymmetric shift and contain information about the growth of harmonic functions. Further, we describe the set of harmonic functions as a vector space isomorphic to the space of formal power series.
△ Less
Submitted 18 October, 2022; v1 submitted 24 January, 2022;
originally announced January 2022.
-
Optimal finite elements for ergodic stochastic two-scale elliptic equations
Authors:
Viet Ha Hoang,
Chen Hui Pang,
Wee Chin Tan
Abstract:
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 solu…
▽ More
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.
△ Less
Submitted 12 January, 2022;
originally announced January 2022.
-
Bayesian Inversion of Log-normal Eikonal Equations
Authors:
Zhan Fei Yeo,
Viet Ha Hoang
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 January, 2023; v1 submitted 22 November, 2021;
originally announced November 2021.
-
Small lobe of comet 67P: Characterization of the Wosret region with ROSETTA-OSIRIS
Authors:
Sonia Fornasier,
Jules Bourdelle de Micas,
Pedro H. Hasselmann,
Van H. Hoang,
Maria Antonietta Barucci,
Holger Sierks
Abstract:
We investigated Wosret, a region located on the small lobe of the 67P/Churyumov-Gerasimenko comet subject to strong heating during the perihelion passage. This region includes Abydos, the final landing site of the Philae lander. We analyzed high-resolution images of the Wosret region acquired between 2015 and 2016 by the OSIRIS instrument on board the Rosetta spacecraft. We observed a few morpholo…
▽ More
We investigated Wosret, a region located on the small lobe of the 67P/Churyumov-Gerasimenko comet subject to strong heating during the perihelion passage. This region includes Abydos, the final landing site of the Philae lander. We analyzed high-resolution images of the Wosret region acquired between 2015 and 2016 by the OSIRIS instrument on board the Rosetta spacecraft. We observed a few morphological changes in Wosret, related to local dust coating removal with an estimated depth of $\sim$ 1 m, along with the formation of a cavity measuring 30 m in length and 6.5 m in depth, for a total estimated mass loss of 1.2 $\times$ 10$^6$ kg. The spectrophotometry of the region is typical of medium-red regions of comet 67P, with spectral slope values of 15-16 \%/(100 nm) in pre-perihelion data acquired at phase angle 60$^o$. Wosret has a spectral phase reddening of 0.0546 $\times 10^{-4}$ nm$^{-1} deg^{-1}$, which is about a factor of 2 lower than what was determined for the nucleus northern hemisphere regions, possibly indicating a reduced surface micro-roughness due to the lack of widespread dust coating. A few tiny bright spots are observed. Morphological features such as "goosebumps" or clods are widely present and larger in size than similar features located in the big lobe. Compared to Anhur and Khonsu, two southern hemisphere regions in the big lobe which are also exposed to high insolation during perihelion, Wosret exhibits fewer exposed volatiles and less morphological variations due to activity events. Our analysis indicates that the small lobe has different physical and mechanical properties than the big one and a lower volatile content, at least in its uppermost layers. These results support the hypothesis that comet 67P originated from the merging of two distinct bodies in the early Solar System.
△ Less
Submitted 2 July, 2021;
originally announced July 2021.
-
Constructing discrete harmonic functions in wedges
Authors:
Viet Hung Hoang,
Kilian Raschel,
Pierre Tarrago
Abstract:
We propose a systematic construction of signed harmonic functions for discrete Laplacian operators with Dirichlet conditions in the quarter plane. In particular, we prove that the set of harmonic functions is an algebra generated by a single element, which conjecturally corresponds to the unique positive harmonic function.
We propose a systematic construction of signed harmonic functions for discrete Laplacian operators with Dirichlet conditions in the quarter plane. In particular, we prove that the set of harmonic functions is an algebra generated by a single element, which conjecturally corresponds to the unique positive harmonic function.
△ Less
Submitted 18 January, 2023; v1 submitted 2 November, 2020;
originally announced December 2020.
-
Adaptive nonparametric estimation of a component density in a two-class mixture model
Authors:
Gaelle Chagny,
Antoine Channarond,
Van Ha Hoang,
Angelina Roche
Abstract:
A two-class mixture model, where the density of one of the components is known, is considered. We address the issue of the nonparametric adaptive estimation of the unknown probability density of the second component. We propose a randomly weighted kernel estimator with a fully data-driven bandwidth selection method, in the spirit of the Goldenshluger and Lepski method. An oracle-type inequality fo…
▽ More
A two-class mixture model, where the density of one of the components is known, is considered. We address the issue of the nonparametric adaptive estimation of the unknown probability density of the second component. We propose a randomly weighted kernel estimator with a fully data-driven bandwidth selection method, in the spirit of the Goldenshluger and Lepski method. An oracle-type inequality for the pointwise quadratic risk is derived as well as convergence rates over Holder smoothness classes. The theoretical results are illustrated by numerical simulations.
△ Less
Submitted 5 February, 2021; v1 submitted 30 July, 2020;
originally announced July 2020.
-
Homogenization of a multiscale multi-continuum system
Authors:
Jun Sur Richard Park,
Viet Ha Hoang
Abstract:
We study homogenization of a locally periodic two-scale dual-continuum system where each continuum interacts with the other. Equations for each continuum are written separately with interaction terms (exchange terms) added. The homogenization limit depends strongly on the scale of this continuum interaction term with respect to the microscopic scale. In J. S. R. Park and V. H. Hoang, {\it Hierarch…
▽ More
We study homogenization of a locally periodic two-scale dual-continuum system where each continuum interacts with the other. Equations for each continuum are written separately with interaction terms (exchange terms) added. The homogenization limit depends strongly on the scale of this continuum interaction term with respect to the microscopic scale. In J. S. R. Park and V. H. Hoang, {\it Hierarchical multiscale finite element method for multicontinuum media}, arXiv:1906.04635, we study in details the case where the interaction terms are scaled as $O(1/\ep^2)$ where $\ep$ is the microscale of the problem. We establish rigorously homogenization limit for this case where we show that in the homogenization limit, the dual-continuum structure disappears. In this paper, we consider the case where this term is scaled as $O(1/\ep)$. This case is far more interesting and difficult as the homogenized problem is a dual-continuum system which contains features that are not in the original two scale problem. In particular, the homogenized dual-continuum system contains extra convection terms and negative interaction coefficients while the interaction coefficient between the continua in the original two scale system obtains both positive and negative values. We prove rigorously the homogenization convergence. We also derive rigorously a homogenization convergence rate. Homogenization of dual-continuum system of this type has not been considered before.
△ Less
Submitted 15 October, 2019;
originally announced October 2019.
-
Multiscale simulations for upscaled multi-continuum flows
Authors:
Jun Sur Richard Park,
Siu Wun Cheung,
Tina Mai,
Viet Ha Hoang
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Hierarchical multiscale finite element method for multi-continuum media
Authors:
Jun Sur Richard Park,
Viet Ha Hoang
Abstract:
Simulation in media with multiple continua where each continuum interacts with every other is often challenging due to multiple scales and high contrast. One needs some types of model reduction. One of the approaches is multi-continuum technique, where every process in each continuum is modeled separately and an interaction term is added. Direct numerical simulation in multi scale media is usually…
▽ More
Simulation in media with multiple continua where each continuum interacts with every other is often challenging due to multiple scales and high contrast. One needs some types of model reduction. One of the approaches is multi-continuum technique, where every process in each continuum is modeled separately and an interaction term is added. Direct numerical simulation in multi scale media is usually not practicable. For this reason, one constructs the corresponding homogenized equations. The paper develops a hierarchical approach for solving cell problems at a dense network of macroscopic points with an essentially optimal computation cost. The method employs the fact that neighboring representative volume elements (RVEs) share similar features; and effective properties of the neighboring RVEs are close to each other. The hierarchical approach reduces computation cost by using different levels of resolution for cell problems at different macroscopic points. The method requires a hierarchy of macroscopic grid points and a corresponding nested approximation spaces with different levels of resolution. Each level of macroscopic points is assigned to an approximation finite element (FE) space which is used to solve the cell problems at the macroscopic points in that level. We prove rigorously that this hierarchical method achieves the same level of accuracy as that of the full solve where cell problems at every macroscopic point are solved using the FE spaces with the highest level of resolution, but at the essentially optimal computation cost. Numerical implementation that computes effective permeabilities of a two scale multicontinuum system via the numerical solutions of the cell problems supports the analytical results. Finally, we prove the homogenization convergence for our multiscale multi-continuum system.
△ Less
Submitted 11 June, 2019;
originally announced June 2019.
-
Linking surface morphology, composition, and activity on the nucleus of 67P/Churyumov-Gerasimenko
Authors:
S. Fornasier,
V. H. Hoang,
P. H. Hasselmann,
C. Feller,
M. A. Barucci,
J. D. P. Deshapriya,
H. Sierks,
G. Naletto,
P. L. Lamy,
R. Rodrigo,
D. Koschny,
B. Davidsson,
J. Agarwal,
C. Barbieri,
J. -L. Bertaux,
I. Bertini,
D. Bodewits,
G. Cremonese,
V. Da Deppo,
S. Debei,
M. De Cecco,
J. Deller,
S. Ferrari,
M. Fulle,
P. J. Gutierrez
, et al. (15 additional authors not shown)
Abstract:
The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October…
▽ More
The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion.
More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg.
We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions.
Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.
△ Less
Submitted 11 September, 2018;
originally announced September 2018.
-
Bayesian inverse problems for recovering coefficients of two scale elliptic equations
Authors:
Viet Ha Hoang,
Jia Hao Quek
Abstract:
We consider the Bayesian inverse homogenization problem of recovering the locally periodic two scale coefficient of a two scale elliptic equation, given limited noisy information on the solution. We consider both the uniform and the Gaussian prior probability measures. We use the two scale homogenized equation whose solution contains the solution of the homogenized equation which describes the mac…
▽ More
We consider the Bayesian inverse homogenization problem of recovering the locally periodic two scale coefficient of a two scale elliptic equation, given limited noisy information on the solution. We consider both the uniform and the Gaussian prior probability measures. We use the two scale homogenized equation whose solution contains the solution of the homogenized equation which describes the macroscopic behaviour, and the corrector which encodes the microscopic behaviour. We approximate the posterior probability by a probability measure determined by the solution of the two scale homogenized equation. We show that the Hellinger distance of these measures converges to zero when the microscale converges to zero, and establish an explicit convergence rate when the solution of the two scale homogenized equation is sufficiently regular. Sampling the posterior measure by Markov Chain Monte Carlo (MCMC) method, instead of solving the two scale equation using fine mesh for each proposal with extremely high cost, we can solve the macroscopic two scale homogenized equation. Although this equation is posed in a high dimensional tensorized domain, it can be solved with essentially optimal complexity by the sparse tensor product finite element method, which reduces the computational complexity of the MCMC sampling method substantially. We show numerically that observations on the macrosopic behaviour alone are not sufficient to infer the microstructure. We need also observations on the corrector. Solving the two scale homogenized equation, we get both the solution to the homogenized equation and the corrector. Thus our method is particularly suitable for sampling the posterior measure of two scale coefficients.
△ Less
Submitted 3 August, 2018; v1 submitted 17 July, 2018;
originally announced July 2018.
-
Nonparametric estimation of the fragmentation kernel based on a PDE stationary distribution approximation
Authors:
Van Ha Hoang,
Thanh Mai Pham Ngoc,
Vincent Rivoirard,
Viet Chi Tran
Abstract:
We consider a stochastic individual-based model in continuous time to describe a size-structured population for cell divisions. This model is motivated by the detection of cellular aging in biology. We address here the problem of nonparametric estimation of the kernel ruling the divisions based on the eigenvalue problem related to the asymptotic behavior in large population. This inverse problem i…
▽ More
We consider a stochastic individual-based model in continuous time to describe a size-structured population for cell divisions. This model is motivated by the detection of cellular aging in biology. We address here the problem of nonparametric estimation of the kernel ruling the divisions based on the eigenvalue problem related to the asymptotic behavior in large population. This inverse problem involves a multiplicative deconvolution operator. Using Fourier technics we derive a nonparametric estimator whose consistency is studied. The main difficulty comes from the non-standard equations connecting the Fourier transforms of the kernel and the parameters of the model. A numerical study is carried out and we pay special attention to the derivation of bandwidths by using resampling.
△ Less
Submitted 25 September, 2020; v1 submitted 25 October, 2017;
originally announced October 2017.
-
High dimensional finite elements for multiscale Maxwell wave equations
Authors:
Van Tiep Chu,
Viet Ha Hoang
Abstract:
We develop an essentially optimal numerical method for solving multiscale Maxwell wave equations in a domain $D\subset{\mathbb R}^d$. The problems depend on $n+1$ scales: one macroscopic scale and $n$ microscopic scales. Solving the macroscopic multiscale homogenized problem, we obtain the desired macroscopic and microscopic information. This problem depends on $n+1$ variables in ${\mathbb R}^d$,…
▽ More
We develop an essentially optimal numerical method for solving multiscale Maxwell wave equations in a domain $D\subset{\mathbb R}^d$. The problems depend on $n+1$ scales: one macroscopic scale and $n$ microscopic scales. Solving the macroscopic multiscale homogenized problem, we obtain the desired macroscopic and microscopic information. This problem depends on $n+1$ variables in ${\mathbb R}^d$, one for each scale that the original multiscale equation depends on, and is thus posed in a high dimensional tensorized domain. The straightforward full tensor product finite element (FE) method is exceedingly expensive. We develop the sparse tensor product FEs that solve this multiscale homogenized problem with essentially optimal number of degrees of freedom, that is essentially equal to that required for solving a macroscopic problem in a domain in ${\mathbb R}^d$ only, for obtaining a required level of accuracy. Numerical correctors are constructed from the FE solution. For two scale problems, we derive a rate of convergence for the numerical corrector in terms of the microscopic scale and the FE mesh width. Numerical examples confirm our analysis.
△ Less
Submitted 6 August, 2017;
originally announced August 2017.
-
Reconciling Bayesian and perimeter regularization for binary inversion
Authors:
Oliver R. A. Dunbar,
Matthew M. Dunlop,
Charles M. Elliott,
Viet Ha Hoang,
Andrew M. Stuart
Abstract:
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 re…
▽ More
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.
△ Less
Submitted 10 April, 2020; v1 submitted 6 June, 2017;
originally announced June 2017.
-
Homogenization of multiscale Maxwell wave equations
Authors:
Van Tiep Chu,
Viet Ha Hoang
Abstract:
We study homogenization of multiscale Maxwell wave equation that depends on $n$ separable microscopic scales in a domain $D\subset{\mathbb R}^d$ on a finite time interval $(0,T)$. Due to the non-compactness of the embedding of $H_0(\curl,D)$ in $L^2(D)^d$, homogenization of Maxwell wave equation can be significantly more complicated than that of scalar wave equations in the $H^1(D)$ setting, and r…
▽ More
We study homogenization of multiscale Maxwell wave equation that depends on $n$ separable microscopic scales in a domain $D\subset{\mathbb R}^d$ on a finite time interval $(0,T)$. Due to the non-compactness of the embedding of $H_0(\curl,D)$ in $L^2(D)^d$, homogenization of Maxwell wave equation can be significantly more complicated than that of scalar wave equations in the $H^1(D)$ setting, and requires analysis uniquely for Maxwell wave equations. We employ multiscale convergence. The homogenized Maxwell wave equation and the initial condition are deduced from the multiscale homogenized equation. When the coefficient of the second order time derivative in the multiscale equation depends on the microscopic scales, the derivation is significantly more complicated, comparing to scalar wave equations, due to the corrector terms for the solution $\ue$ of the multiscale equation in the $L^2(D)^d$ norm, which do not appear in the $H^1(D)$ setting. For two scale equations, we derive an explicit homogenization error estimate for the case where the solution $u_0$ of the homogenized equation belongs to $L^\infty((0,T);H^1(\curl,D))$. When $u_0$ only belongs to a weaker regularity space $L^\infty((0,T);H^s(\curl,D))$ for $0<s<1$, we contribute an approach to deduce a new homogenization error in this case, which depends on $s$. For general multiscale problems, a corrector is derived albeit without an explicit homogenization error estimate. These correctors and homogenization errors play an essential role in deriving numerical correctors for approximating the solutions to the multiscale problems numerically, as considered in our forthcoming publication.
△ Less
Submitted 21 May, 2017;
originally announced May 2017.
-
Bayesian Multiscale Finite Element Methods. Modeling missing subgrid information probabilistically
Authors:
Y. Efendiev,
W. T. Leung,
S. W. Cheung,
N. Guha,
V. H. Hoang,
B. Mallick
Abstract:
In this paper, we develop a Bayesian multiscale approach based on a multiscale finite element method. Because of scale disparity in many multiscale applications, computational models can not resolve all scales. Various subgrid models are proposed to represent un-resolved scales. Here, we consider a probabilistic approach for modeling un-resolved scales using the Multiscale Finite Element Method (c…
▽ More
In this paper, we develop a Bayesian multiscale approach based on a multiscale finite element method. Because of scale disparity in many multiscale applications, computational models can not resolve all scales. Various subgrid models are proposed to represent un-resolved scales. Here, we consider a probabilistic approach for modeling un-resolved scales using the Multiscale Finite Element Method (cf., [1, 2]). By representing dominant modes using the Generalized Multiscale Finite Element, we propose a Bayesian framework, which provides multiple inexpensive (computable) solutions for a deterministic problem. These approximate probabilistic solutions may not be very close to the exact solutions and, thus, many realizations are needed. In this way, we obtain a rigorous probabilistic description of approximate solutions. In the paper, we consider parabolic and wave equations in heterogeneous media. In each time interval, the domain is divided into subregions. Using residual information, we design appropriate prior and posterior distributions. The likelihood consists of the residual minimization. To sample from the resulting posterior distribution, we consider several sampling strategies. The sampling involves identifying important regions and important degrees of freedom beyond permanent basis functions, which are used in residual computation. Numerical results are presented. We consider two sampling algorithms. The first algorithm uses sequential sampling and is inexpensive. In the second algorithm, we perform full sampling using the Gibbs sampling algorithm, which is more accurate compared to the sequential sampling. The main novel ingredients of our approach consist of: defining appropriate permanent basis functions and the corresponding residual; setting up a proper posterior distribution; and sampling the posteriors.
△ Less
Submitted 9 February, 2017;
originally announced February 2017.
-
Adaptive wavelet multivariate regression with errors in variables
Authors:
Michaël Chichignoud,
Van Ha Hoang,
Thanh Mai Pham Ngoc,
Vincent Rivoirard
Abstract:
In the multidimensional setting, we consider the errors-in-variables model. We aim at estimating the unknown nonparametric multivariate regression function with errors in the covariates. We devise an adaptive estimator based on projection kernels on wavelets and a deconvolution operator. We propose an automatic and fully data driven procedure to select the wavelet level resolution. We obtain an or…
▽ More
In the multidimensional setting, we consider the errors-in-variables model. We aim at estimating the unknown nonparametric multivariate regression function with errors in the covariates. We devise an adaptive estimator based on projection kernels on wavelets and a deconvolution operator. We propose an automatic and fully data driven procedure to select the wavelet level resolution. We obtain an oracle inequality and optimal rates of convergence over anisotropic H{ö}lder classes. Our theoretical results are illustrated by some simulations.
△ Less
Submitted 12 January, 2016;
originally announced January 2016.
-
Homogenization error for two scale Maxwell equations
Authors:
Van Tiep Chu,
Viet Ha Hoang
Abstract:
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 finit…
▽ More
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$.
△ Less
Submitted 9 December, 2015;
originally announced December 2015.
-
Polynomial approximations of a class of stochastic multiscale elasticity problems
Authors:
Viet Ha Hoang,
Thanh Chung Nguyen,
Bingxing Xia
Abstract:
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…
▽ More
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.
△ Less
Submitted 28 November, 2015;
originally announced November 2015.
-
Estimating the Division Kernel of a Size-Structured Population
Authors:
Van Ha Hoang
Abstract:
We consider a size-structured population describing the cell divisions. The cell population is described by an empirical measure and we observe the divisions in the continuous time interval [0, T ]. We address here the problem of estimating the division kernel h (or fragmentation kernel) in case of complete data. An adaptive estimator of h is constructed based on a kernel function K with a fully d…
▽ More
We consider a size-structured population describing the cell divisions. The cell population is described by an empirical measure and we observe the divisions in the continuous time interval [0, T ]. We address here the problem of estimating the division kernel h (or fragmentation kernel) in case of complete data. An adaptive estimator of h is constructed based on a kernel function K with a fully data-driven bandwidth selection method. We obtain an oracle inequality and an exponential convergence rate, for which optimality is considered.
△ Less
Submitted 19 May, 2016; v1 submitted 9 September, 2015;
originally announced September 2015.
-
Determining White Noise Forcing From Eulerian Observations in the Navier Stokes Equation
Authors:
Viet Ha Hoang,
Kody J. H. Law,
Andrew M. Stuart
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 April, 2014; v1 submitted 19 March, 2013;
originally announced March 2013.
-
Complexity Analysis of Accelerated MCMC Methods for Bayesian Inversion
Authors:
Viet Ha Hoang,
Christoph Schwab,
Andrew M. Stuart
Abstract:
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…
▽ More
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.}
△ Less
Submitted 30 April, 2013; v1 submitted 10 July, 2012;
originally announced July 2012.
-
Bayesian inverse problems for Burgers and Hamilton-Jacobi equations with white-noise forcing
Authors:
Viet Ha Hoang
Abstract:
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. Inferenc…
▽ More
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.
△ Less
Submitted 10 May, 2011; v1 submitted 14 April, 2011;
originally announced April 2011.
-
Homogenization for advection-diffusion in a perforated domain
Authors:
P. H. Haynes,
V. H. Hoang,
J. R. Norris,
K. C. Zygalakis
Abstract:
The volume of a Wiener sausage constructed from a diffusion process with periodic, mean-zero, divergence-free velocity field, in dimension 3 or more, is shown to have a non-random and positive asymptotic rate of growth. This is used to establish the existence of a homogenized limit for such a diffusion when subject to Dirichlet conditions on the boundaries of a sparse and independent array of obst…
▽ More
The volume of a Wiener sausage constructed from a diffusion process with periodic, mean-zero, divergence-free velocity field, in dimension 3 or more, is shown to have a non-random and positive asymptotic rate of growth. This is used to establish the existence of a homogenized limit for such a diffusion when subject to Dirichlet conditions on the boundaries of a sparse and independent array of obstacles. There is a constant effective long-time loss rate at the obstacles. The dependence of this rate on the form and intensity of the obstacles and on the velocity field is investigated. A Monte Carlo algorithm for the computation of the volume growth rate of the sausage is introduced and some numerical results are presented for the Taylor--Green velocity field.
△ Less
Submitted 25 March, 2010; v1 submitted 21 March, 2010;
originally announced March 2010.