Related papers: Multilevel Asymptotic-Preserving Monte Carlo for P…
The stochastic particle method based on Bhatnagar-Gross-Krook (BGK) or ellipsoidal statistical BGK (ESBGK) model approximates the pairwise collisions in the Boltzmann equation using a relaxation process. Therefore, it is more efficient to…
Ideal gases can be modeled by the Boltzmann equation from statistical physics. Instead of trying to track the position and velocity of a large number of gas molecules, it is possible to describe the particles with a particle distribution…
This article is concerned with the multilevel Monte Carlo (MLMC) methods for approximating expectations of some functions of the solution to the Heston 3/2-model from mathematical finance, which takes values in $(0, \infty)$ and possesses…
In this article we develop a new sequential Monte Carlo (SMC) method for multilevel (ML) Monte Carlo estimation. In particular, the method can be used to estimate expectations with respect to a target probability distribution over an…
In this paper, we generalize the Quantum Kinetic Monte Carlo (QKMC) method of the Schrodinger equation, which was first proposed by [Z. Cai and J. Lu. SIAM J. Sci. Comput., 40(3):B706-B722, 2018] to the Lindblad equation. This algorithm…
In some recent works [G. Dimarco, L. Pareschi, Hybrid multiscale methods I. Hyperbolic Relaxation Problems, Comm. Math. Sci., 1, (2006), pp. 155-177], [G. Dimarco, L. Pareschi, Hybrid multiscale methods II. Kinetic equations, SIAM…
Partial differential equations (p.d.e) equipped of spatial derivatives of fractional order capture anomalous transport behaviors observed in diverse fields of Science. A number of numerical methods approximate their solutions in dimension…
In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations.…
Estimating risk measures such as large loss probabilities and Value-at-Risk is fundamental in financial risk management and often relies on computationally intensive nested Monte Carlo methods. While Multi-Level Monte Carlo (MLMC)…
We consider the computational efficiency of Monte Carlo (MC) and Multilevel Monte Carlo (MLMC) methods applied to partial differential equations with random coefficients. These arise, for example, in groundwater flow modelling, where a…
In this paper, a physics-oriented stochastic kinetic scheme will be developed that includes random inputs from both flow and electromagnetic fields via a hybridization of stochastic Galerkin and collocation methods. Based on the BGK-type…
In this comprehensive and detailed study, vacancy-mediated self-diffusion of A- and B-elements in 'triple-defect' B2-ordered ASB(1-S) binaries is simulated by means of a kinetic Monte Carlo (KMC) algorithm involving atomic jumps to…
Monte Carlo simulations of quantum field theories on a lattice become increasingly expensive as the continuum limit is approached since the cost per independent sample grows with a high power of the inverse lattice spacing. Simulations on…
We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that are modeled using the Richards' equation. We propose a stochastic extension for the empirical models that are…
We present a method to model the interaction and the dynamics of atoms excited to Rydberg states. We show a way to solve the optical Bloch equations for laser excitation of the frozen gas in good agreement with the experiment. A second…
Monte Carlo simulation is an unbiased numerical tool for studying classical and quantum many-body systems. One of its bottlenecks is the lack of general and efficient update algorithm for large size systems close to phase transition or with…
The design of particle simulation methods for collisional plasma physics has always represented a challenge due to the unbounded total collisional cross section, which prevents a natural extension of the classical Direct Simulation Monte…
The Multilevel Monte Carlo (MLMC) approach usually works well when estimating the expected value of a quantity which is a Lipschitz function of intermediate quantities, but if it is a discontinuous function it can lead to a much slower…
The Direct Simulation Monte Carlo method is applied to solve the Boltzmann equation in the steady planar Couette flow for Maxwell molecules and hard spheres. Nonequilibrium boundary conditions based on the solution of the…
In this article we consider static Bayesian parameter estimation for partially observed diffusions that are discretely observed. We work under the assumption that one must resort to discretizing the underlying diffusion process, for…