Two critical issues in Langevin simulation of gas flows

Jun Zhang

Research output: Chapter in Book/Report/Conference proceedingChapter (peer-reviewed)

Abstract

A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.
LanguageEnglish
Title of host publicationRarefied Gas Dynamics 2014
Subtitle of host publicationProceedings of the 29th International Symposium on Rarefied Gas Dynamics
Pages1087-1094
Number of pages8
Volume1328
DOIs
Publication statusPublished - Sep 2014
Event29th International Symposium on Rarefied Gas Dynamics - Hilton Xi’an Hotel, Xi'an, China
Duration: 13 Jul 201418 Jul 2014

Publication series

NameAIP Conference Proceedings
PublisherAmerican Institute of Physics
Volume1328

Conference

Conference29th International Symposium on Rarefied Gas Dynamics
CountryChina
CityXi'an
Period13/07/1418/07/14

Fingerprint

Gas Flow
Flow of gases
Monte Carlo methods
Direct Simulation Monte Carlo
Langevin Equation
Simulation
Simulation Methods
Monte Carlo method
Kinetic theory
Knudsen number
Couette Flow
Flow velocity
Transport properties
Boundary conditions
Collision
Rarefied Gas Flow
Molecules
Green's Formula
Hitting Time
Cell Size

Keywords

  • Langevin dynamics
  • Fokker-Planck equation
  • transport properties
  • DSMC method

Cite this

Zhang, J. (2014). Two critical issues in Langevin simulation of gas flows. In Rarefied Gas Dynamics 2014: Proceedings of the 29th International Symposium on Rarefied Gas Dynamics (Vol. 1328, pp. 1087-1094). (AIP Conference Proceedings; Vol. 1328). https://doi.org/10.1063/1.4902714
Zhang, Jun. / Two critical issues in Langevin simulation of gas flows. Rarefied Gas Dynamics 2014: Proceedings of the 29th International Symposium on Rarefied Gas Dynamics. Vol. 1328 2014. pp. 1087-1094 (AIP Conference Proceedings).
@inbook{41285d0b9cf846d9a13d927e7e3cc863,
title = "Two critical issues in Langevin simulation of gas flows",
abstract = "A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-B{\'e}nard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.",
keywords = "Langevin dynamics, Fokker-Planck equation, transport properties, DSMC method",
author = "Jun Zhang",
year = "2014",
month = "9",
doi = "10.1063/1.4902714",
language = "English",
isbn = "978-0-7354-1265-1",
volume = "1328",
series = "AIP Conference Proceedings",
publisher = "American Institute of Physics",
pages = "1087--1094",
booktitle = "Rarefied Gas Dynamics 2014",

}

Zhang, J 2014, Two critical issues in Langevin simulation of gas flows. in Rarefied Gas Dynamics 2014: Proceedings of the 29th International Symposium on Rarefied Gas Dynamics. vol. 1328, AIP Conference Proceedings, vol. 1328, pp. 1087-1094, 29th International Symposium on Rarefied Gas Dynamics, Xi'an, China, 13/07/14. https://doi.org/10.1063/1.4902714

Two critical issues in Langevin simulation of gas flows. / Zhang, Jun.

Rarefied Gas Dynamics 2014: Proceedings of the 29th International Symposium on Rarefied Gas Dynamics. Vol. 1328 2014. p. 1087-1094 (AIP Conference Proceedings; Vol. 1328).

Research output: Chapter in Book/Report/Conference proceedingChapter (peer-reviewed)

TY - CHAP

T1 - Two critical issues in Langevin simulation of gas flows

AU - Zhang, Jun

PY - 2014/9

Y1 - 2014/9

N2 - A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.

AB - A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.

KW - Langevin dynamics

KW - Fokker-Planck equation

KW - transport properties

KW - DSMC method

UR - http://rgd29.org/wordpress/wp-content/uploads/2014/12/RGD29-Proceedings.Table-of-Contents.pdf

U2 - 10.1063/1.4902714

DO - 10.1063/1.4902714

M3 - Chapter (peer-reviewed)

SN - 978-0-7354-1265-1

VL - 1328

T3 - AIP Conference Proceedings

SP - 1087

EP - 1094

BT - Rarefied Gas Dynamics 2014

ER -

Zhang J. Two critical issues in Langevin simulation of gas flows. In Rarefied Gas Dynamics 2014: Proceedings of the 29th International Symposium on Rarefied Gas Dynamics. Vol. 1328. 2014. p. 1087-1094. (AIP Conference Proceedings). https://doi.org/10.1063/1.4902714