Related papers: Coupling Fluid Plasma and Kinetic Neutral Models u…
Multi-fidelity Monte Carlo methods leverage low-fidelity and surrogate models for variance reduction to make tractable uncertainty quantification even when numerically simulating the physical systems of interest with high-fidelity models is…
Decoupling approach presents a novel solution/alternative to the highly time-consuming fluid-thermal-structural simulation procedures when thermal effects and resultant displacements on machine tools are analyzed. Using high dimensional…
We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of…
We propose a Multi-Cell Monte Carlo algorithm, or (MC)^2, for predicting stable phases in chemically complex crystalline systems. Free atomic transfer among cells is achieved via the application of the lever rule, where an assigned molar…
We present a promising coarse-graining strategy for linking micro- and mesoscales of soft matter systems. The approach is based on effective pairwise interaction potentials obtained from detailed atomistic molecular dynamics (MD)…
Phase separation plays a central role in the emergence of novel functionalities of correlated electron materials. The structure of the mixed-phase states depends strongly on the nonequilibrium phase-separation dynamics, which has so far yet…
In many applications, such as plasma edge simulation of a nuclear fusion reactor, a coupled PDE/kinetic description is required, which is usually solved with a coupled finite-volume/Monte-Carlo method. Different procedures have been…
This paper presents two novel ensemble domain decomposition methods for fast-solving the Stokes-Darcy coupled models with random hydraulic conductivity and body force. To address such random systems, we employ the Monte Carlo (MC) method to…
Kinetic Monte Carlo (KMC) is a powerful method for simulation of diffusion processes in various systems. The accuracy of the method, however, relies on the extent of details used for the parameterization of the model. Migration barriers are…
We develop a new Monte Carlo method that solves hyperbolic transport equations with stiff terms, characterized by a (small) scaling parameter. In particular, we focus on systems which lead to a reduced problem of parabolic type in the limit…
Adjoint methods form a class of importance sampling methods that are used to accelerate Monte Carlo (MC) simulations of transport equations. Ideally, adjoint methods allow for zero-variance MC estimators provided that the solution to an…
In this paper, we evaluate the performance of the multilevel Monte Carlo method (MLMC) for deterministic and uncertain hyperbolic systems, where randomness is introduced either in the modeling parameters or in the approximation algorithms.…
Neutrino-matter interactions play an important role in the post-merger evolution of neutron star-neutron star and black hole-neutron star mergers. Most notably, they determine the properties of the bright optical/infrared transients…
Restricted Boltzmann Machines are simple and powerful generative models that can encode any complex dataset. Despite all their advantages, in practice the trainings are often unstable and it is difficult to assess their quality because the…
In this paper, we focus on the construction of a hybrid scheme for the approximation of non-Maxwellian kinetic models with uncertainties. In the context of multiagent systems, the introduction of a kernel at the kinetic level is useful to…
Diffusion Monte Carlo (DMC) based on fixed-node approximation has enjoyed significant developments in the past decades and become one of the go-to methods when accurate ground state energy of molecules and materials is needed. The remaining…
One dimensional (1D) simulations of the flow and flooding of open channels are known to be inaccurate as the flow is multi-dimensional in nature, especially at the flooded regions. However, multi-dimensional simulations, even in two…
Kinetic descriptions of runaway electrons (RE) are usually based on Fokker-Planck models that determine the probability distribution function (PDF) of RE in 2-dimensional momentum space. Despite of the simplification involved, the…
Bounded plasmas are characterized by a rapid but smooth transition from quasi-neutrality in the volume to electron depletion close to the electrodes and chamber walls. The thin non-neutral region, the boundary sheath, comprises only a small…
Strongly correlated fermionic systems are of great interest in condensed matter physics and numerical methods are indispensable tools for their study. However, existing approaches such as exact diagonalization (ED) and stochastic quantum…