arrow
Volume 20, Issue 5
A Probabilistic Automatic Steady State Detection Method for the Direct Simulation Monte Carlo

A. Karchani, O. Ejtehadi & R. S. Myong

Commun. Comput. Phys., 20 (2016), pp. 1183-1209.

Published online: 2018-04

Export citation
  • Abstract

The statistical error associated with sampling in the DSMC method can be categorized as type I and II, which are caused by the incorrect rejection and acceptance of the null hypothesis, respectively. In this study, robust global and local automatic steady state detection methods were developed based on an ingenious method based purely on the statistics and kinetics of particles. The key concept is built upon probabilistic automatic reset sampling (PARS) to minimize the type II error caused by incorrect acceptance of the samples that do not belong to the steady state. The global steady state method is based on a relative standard variation of collisional invariants, while the local steady state method is based on local variations in the distribution function of particles at each cell. In order to verify the capability of the new methods, two benchmark cases – the one-dimensional shear-driven Couette flow and the two-dimensional high speed flow past a vertical wall – were extensively investigated. Owing to the combined effects of the automatic detection and local reset sampling, the local steady state detection method yielded a substantial gain of 30-36% in computational cost for the problem studied. Moreover, the local reset feature outperformed the automatic detection feature in overall computational savings.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-20-1183, author = {}, title = {A Probabilistic Automatic Steady State Detection Method for the Direct Simulation Monte Carlo}, journal = {Communications in Computational Physics}, year = {2018}, volume = {20}, number = {5}, pages = {1183--1209}, abstract = {

The statistical error associated with sampling in the DSMC method can be categorized as type I and II, which are caused by the incorrect rejection and acceptance of the null hypothesis, respectively. In this study, robust global and local automatic steady state detection methods were developed based on an ingenious method based purely on the statistics and kinetics of particles. The key concept is built upon probabilistic automatic reset sampling (PARS) to minimize the type II error caused by incorrect acceptance of the samples that do not belong to the steady state. The global steady state method is based on a relative standard variation of collisional invariants, while the local steady state method is based on local variations in the distribution function of particles at each cell. In order to verify the capability of the new methods, two benchmark cases – the one-dimensional shear-driven Couette flow and the two-dimensional high speed flow past a vertical wall – were extensively investigated. Owing to the combined effects of the automatic detection and local reset sampling, the local steady state detection method yielded a substantial gain of 30-36% in computational cost for the problem studied. Moreover, the local reset feature outperformed the automatic detection feature in overall computational savings.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.080815.240316a}, url = {http://global-sci.org/intro/article_detail/cicp/11187.html} }
TY - JOUR T1 - A Probabilistic Automatic Steady State Detection Method for the Direct Simulation Monte Carlo JO - Communications in Computational Physics VL - 5 SP - 1183 EP - 1209 PY - 2018 DA - 2018/04 SN - 20 DO - http://doi.org/10.4208/cicp.080815.240316a UR - https://global-sci.org/intro/article_detail/cicp/11187.html KW - AB -

The statistical error associated with sampling in the DSMC method can be categorized as type I and II, which are caused by the incorrect rejection and acceptance of the null hypothesis, respectively. In this study, robust global and local automatic steady state detection methods were developed based on an ingenious method based purely on the statistics and kinetics of particles. The key concept is built upon probabilistic automatic reset sampling (PARS) to minimize the type II error caused by incorrect acceptance of the samples that do not belong to the steady state. The global steady state method is based on a relative standard variation of collisional invariants, while the local steady state method is based on local variations in the distribution function of particles at each cell. In order to verify the capability of the new methods, two benchmark cases – the one-dimensional shear-driven Couette flow and the two-dimensional high speed flow past a vertical wall – were extensively investigated. Owing to the combined effects of the automatic detection and local reset sampling, the local steady state detection method yielded a substantial gain of 30-36% in computational cost for the problem studied. Moreover, the local reset feature outperformed the automatic detection feature in overall computational savings.

A. Karchani, O. Ejtehadi & R. S. Myong. (2020). A Probabilistic Automatic Steady State Detection Method for the Direct Simulation Monte Carlo. Communications in Computational Physics. 20 (5). 1183-1209. doi:10.4208/cicp.080815.240316a
Copy to clipboard
The citation has been copied to your clipboard