10 research outputs found
Importance Sampling Variance Reduction for the Fokker-Planck Rarefied Gas Particle Method
Models and methods that are able to accurately and efficiently predict the
flows of low-speed rarefied gases are in high demand, due to the increasing
ability to manufacture devices at micro and nano scales. One such model and
method is a Fokker-Planck approximation to the Boltzmann equation, which can be
solved numerically by a stochastic particle method. The stochastic nature of
this method leads to noisy estimates of the thermodynamic quantities one wishes
to sample when the signal is small in comparison to the thermal velocity of the
gas. Recently, Gorji et al have proposed a method which is able to greatly
reduce the variance of the estimators, by creating a correlated stochastic
process which acts as a control variate for the noisy estimates. However, there
are potential difficulties involved when the geometry of the problem is
complex, as the method requires the density to be solved for independently.
Importance sampling is a variance reduction technique that has already been
shown to successfully reduce the noise in direct simulation Monte Carlo
calculations. In this paper we propose an importance sampling method for the
Fokker-Planck stochastic particle scheme. The method requires minimal change to
the original algorithm, and dramatically reduces the variance of the estimates.
We test the importance sampling scheme on a homogeneous relaxation, planar
Couette flow and a lid-driven-cavity flow, and find that our method is able to
greatly reduce the noise of estimated quantities. Significantly, we find that
as the characteristic speed of the flow decreases, the variance of the noisy
estimators becomes independent of the characteristic speed
A kinetic FokkerâPlanck approach for modeling variable hard-sphere gas mixtures
Kinetic Fokker-Planck (FP) methods for modeling rarefied gas flows have received increasing attention over the last few years. However, formulating such models for realistic multi-species gases is still an open subject of research. Therefore, in this letter, we develop a kinetic FP model for describing gas mixtures with particles interacting according to the variable hard-sphere interaction potential. In accordance with the kinetic FP framework, a stochastic solution algorithm is employed in order to solve the model on a particle level. Different test cases are carried out, and the performance of the proposed method is compared with the direct simulation Monte Carlo algorithm
Treatment of long-range interactions arising in the Enskog-Vlasov description of dense fluids
The kinetic theory of rarefied gases and numerical schemes based on the Boltzmann equation, have evolved to the cornerstone of non-equilibrium gas dynamics. However, their counterparts in the dense regime remain rather exotic for practical non-continuum scenarios. This problem is partly due to the fact that long-range interactions arising from the attractive tail of molecular potentials, lead to a computationally demanding Vlasov integral. This study focuses on numerical remedies for efficient stochastic particle simulations based on the EnskogâVlasov kinetic equation. In particular, we devise a Poisson type elliptic equation which governs the underlying long-range interactions. The idea comes through fitting a Green function to the molecular potential, and hence deriving an elliptic equation for the associated fundamental solution. Through this transformation of the Vlasov integral, efficient Poisson type solvers can be readily employed in order to compute the mean field forces. Besides the technical aspects of different numerical schemes for treatment of the Vlasov integral, simulation results for evaporation of a liquid slab into the vacuum are presented. It is shown that the proposed formulation leads to accurate predictions with a reasonable computational cost
Controlling the bias error of Fokker- Planck methods for rarefied gas dynamics simulations
Direct simulation Monte-Carlo (DSMC) is the most established method for rarefied gas flow simulations. It is valid from continuum to near vacuum, but in cases involving small Knudsen numbers (Kn), it suffers from high computational cost. The Fokker-Planck (FP) method, on the other hand, is almost as accurate as DSMC for small to moderate Kn, but it does not have the computational drawback of DSMC, if Kn is small [P. Jenny, M. Torrilhon, and S. Heinz, âA solution algorithm for the fluid dynamic equations based on a stochastic model for molecular motion,â J. Comput. Phys. 229, 1077â1098 (2010) and H. Gorji, M. Torrilhon, and P. Jenny, âFokkerâPlanck model for computational studies of monatomic rarefied gas flows,â J. Fluid Mech. 680, 574â601 (2011)]. Especially attractive is the combination of the two approaches leading to the FP-DSMC method. Opposed to other hybrid methods, e.g., coupled DSMC/Navier-Stokes solvers, it is relatively straightforward to couple DSMC with the FP method since both are based on particle solution algorithms sharing the same data structure and having similar components. Regarding the numerical accuracy of such particle methods, one has to distinguish between spatial truncation errors, time stepping errors, statistical errors and bias errors. In this paper, the bias error of the FP method is analyzed in detail, and it is shown how it can be reduced without increasing the particle number to an exorbitant level. The effectiveness of the discussed bias error reduction scheme is demonstrated for uniform shear flow, for which an analytical reference solution was derived