On Monte Carlo Hybrid Methods for Linear Algebra

Author(s):  
Diego Davila ◽  
Vassil Alexandrov ◽  
Oscar A. Esquivel-Flores
1998 ◽  
Vol 120 (3) ◽  
pp. 547-560 ◽  
Author(s):  
J. R. Howell

The use of the Monte Carlo method in radiative heat transfer is reviewed. The review covers surface-surface, enclosure, and participating media problems. Discussion is included of research on the fundamentals of the method and on applications to surface-surface interchange in enclosures, exchange between surfaces with roughness characteristics, determination of configuration factors, inverse design, transfer through packed beds and fiber layers, participating media, scattering, hybrid methods, spectrally dependent problems including media with line structure, effects of using parallel algorithms, practical applications, and extensions of the method. Conclusions are presented on needed future work and the place of Monte Carlo techniques in radiative heat transfer computations.


2017 ◽  
Vol 145 (7) ◽  
pp. 2533-2553 ◽  
Author(s):  
Andreas S. Stordal ◽  
Hans A. Karlsen

In high-dimensional dynamic systems, standard Monte Carlo techniques that asymptotically reproduce the posterior distribution are computationally too expensive. Alternative sampling strategies are usually applied and among these the ensemble Kalman filter (EnKF) is perhaps the most popular. However, the EnKF suffers from severe bias if the model under consideration is far from linear. Another class of sequential Monte Carlo methods is kernel-based Gaussian mixture filters, which reduce the bias but maintain the robustness of the EnKF. Although many hybrid methods have been introduced in recent years, not many have been analyzed theoretically. Here it is shown that the recently proposed adaptive Gaussian mixture filter can be formulated in a rigorous Bayesian framework and that the algorithm can be generalized to a broader class of interpolated kernel filters. Two parameters—the bandwidth of the kernel and a weight interpolation factor—determine the filter performance. The new formulation of the filter includes particle filters, EnKF, and kernel-based Gaussian mixture filters as special cases. Techniques from particle filter literature are used to calculate the asymptotic bias of the filter as a function of the parameters and to derive a central limit theorem. The asymptotic theory is then used to determine the parameters as a function of the sample size in a robust way such that the error norm vanishes asymptotically, whereas the normalized error is sample independent and bounded. The parameter choice is tested on the Lorenz 63 model, where it is shown that the error is smaller or equal to the EnKF and the optimal particle filter for a varying sample size.


1998 ◽  
Vol 334 (1-2) ◽  
pp. 209-213
Author(s):  
Jun-ichiro Takano ◽  
Osamu Takai ◽  
Yoshiaki Kogure ◽  
Masao Doyama

Author(s):  
Ivan Dimov ◽  
Vassil Alexandrov ◽  
Rumyana Papancheva ◽  
Christian Weihrauch

Sign in / Sign up

Export Citation Format

Share Document