Hostname: page-component-745bb68f8f-5r2nc Total loading time: 0 Render date: 2025-01-23T12:13:29.558Z Has data issue: false hasContentIssue false
  • English
  • Français

A Particle Fokker-Planck Algorithm with Multiscale Temporal Discretization for Rarefied and Continuum Gas Flows

Published online by Cambridge University Press:  21 June 2017

Fei Fei ,
Zhaohui Liu ,
Jun Zhang  and
Chuguang Zheng
Fei Fei*
Affiliation:
School of Aerospace Engineering, Huazhong University of Science and Technology, Wuhan 430074, P.R. China State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, P.R. China
Zhaohui Liu*
Affiliation:
State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, P.R. China
Jun Zhang*
Affiliation:
School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, United Kingdom School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, P.R. China
Chuguang Zheng*
Affiliation:
State Key Laboratory of Coal Combustion, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan 430074, P.R. China
*
*Corresponding author. Email addresses:[email protected] (F. Fei), [email protected] (Z. Liu), [email protected] (J. Zhang), [email protected] (C. Zheng)
*Corresponding author. Email addresses:[email protected] (F. Fei), [email protected] (Z. Liu), [email protected] (J. Zhang), [email protected] (C. Zheng)
*Corresponding author. Email addresses:[email protected] (F. Fei), [email protected] (Z. Liu), [email protected] (J. Zhang), [email protected] (C. Zheng)
*Corresponding author. Email addresses:[email protected] (F. Fei), [email protected] (Z. Liu), [email protected] (J. Zhang), [email protected] (C. Zheng)
Get access

Abstract

For gas flows with moderate and low Knudsen numbers, pair-wise collisions in the Boltzmann equation can be approximated by the Langevin model corresponding to the Fokker-Planck equation. Using this simplified collision model, particle numerical schemes, e.g. the Fokker-Planck model (FPM) method, can simulate low Knudsen number gas flows more efficient than those based on the Boltzmann equation, such as the Direct Simulation Monte Carlo (DSMC) method. However, as analyzed in this paper, the transport properties of the FPM method deviate from the physical values as the time step increases, and this problem affects its computational accuracy and efficiency for the simulation of multi-scale flows. Herewe propose a particle Fokker-Planck algorithm with multiscale temporal discretization (MTD-FPM) to overcome the drawbacks of the original FPM method. In the MTD-FPM method, the molecular motion is tracked following the integration scheme of the Langevin model in analogy to the original FPM method. However, to ensure consistent transport coefficients for arbitrary temporal discretization, a time step dependent friction coefficient has been implemented. Several benchmark problems, including Couette, thermal Couette, Poiseuille, and Sod tube flows, are simulated to validate the proposed MTD-FPM method.

Keywords

Particle simulation methodFokker-Planck modeltemporal discretizationmultiscale flow
Type
Research Article
Information
Communications in Computational Physics , Volume 22 , Issue 2 , August 2017 , pp. 338 - 374
Copyright
Copyright © Global-Science Press 2017 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Footnotes

Communicated by Kun Xu

References

[1] La Torre, F., Moerel, J. L., and Kleijn, C. R., Hybrid simulations of rarefied supersonic gas flows in micro-nozzles, Comput. Fluids, 49 (2011), 312–322.Google Scholar
[2] Titov, E., Gallagher-Rogers, A., Levin, D., and Reed, B., Examination of a collision-limiter direct simulation Monte Carlo method for micropropulsion applications, J. Propul. Power, 24 (2008), 311–321.Google Scholar
[3] Ivanov, M. S. and Gimelshein, S. F., Computational hypersonic rarefied flows, Ann. Rev. Fluid Mech., 30 (1998), 469–505.Google Scholar
[4] Votta, R., Schettino, A., and Bonfiglioli, A., Hypersonic high altitude aerothermodynamics of a space re-entry vehicle, Aerosp. Sci. Technol., 25 (2013), 253–265.CrossRefGoogle Scholar
[5] Bird, G. A., Molecular Gas Dynamics and Direct Simulation of Gas Flows, Clarendon Press, Oxford, 1994.Google Scholar
[6] Hash, D. B. and Hassan, H. A., Assessment of schemes for coupling Monte Carlo and Navier-Stokes solution methods, J. Thermophys. Heat Transf., 10 (1996), 242–249.CrossRefGoogle Scholar
[7] Sun, Q., Boyd, I. D., and Candler, G. V., A hybrid continuum/particle approach for modeling rarefied gas flows, J. Comput. Phys., 194 (2004), 256–277.CrossRefGoogle Scholar
[8] Schwartzentruber, T. E. and Boyd, I. D., A hybrid particle-continuum method applied to shock waves, J. Comput. Phys., 215 (2006), 402–416.Google Scholar
[9] Schwartzentruber, T. E., Scalabrin, L. C., and Boyd, I. D., A modular particle-continuum numerical method for hypersonic non-equilibrium gas flows, J. Comput. Phys., 225 (2007), 1159–1174.CrossRefGoogle Scholar
[10] Wijesinghe, H., and Hadjiconstantinou, N., Discussion of hybrid atomistic-continuum methods for multiscale hydrodynamics, Int. J. Multi Comp. Eng., 2 (2004), 189–202.CrossRefGoogle Scholar
[11] Xu, K. and Huang, J. C., A unified gas-kinetic scheme for continuum and rarefied flows, J. Comput. Phys., 229 (2010), 7747–7764.Google Scholar
[12] Huang, J. C., Xu, K., and Yu, P. B., A unified gas-kinetic scheme for continuum and rarefied flows. III: Microflow simulations, Commun. Comput. Phys., 14 (2013), 1147–1173.Google Scholar
[13] Liu, C., Xu, K., Sun, Q., and Cai, Q., A unified gas-kinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations, J. Comput. Phys., 314 (2016), 305–340.Google Scholar
[14] Jenny, P., Torrilhon, M., and Heinz, S., A solution algorithm for the fluid dynamic equations based on a stochastic model for molecular motion, J. Comput. Phys., 229 (2010), 1077–1098.Google Scholar
[15] Gorji, M. H., Torrilhon, M., and Jenny, P., Fokker-Planck model for computational studies of monatomic rarefied gas flows, J. Fluid Mech., 680 (2011), 574–601.Google Scholar
[16] Gorji, M.H. and Jenny, P., An efficient particle Fokker-Planck algorithm for rarefied gas flows, J. Comput. Phys., 263 (2014), 325–343.Google Scholar
[17] Gorji, M. H. and Jenny, P., A kinetic model for gas mixtures based on a Fokker-Planck Equation, J. Phys. Conf. Ser., 362 (2012), 012042.Google Scholar
[18] Gorji, M. H. and Jenny, P., A Fokker-Planck based kinetic model for diatomic rarefied gas flows, Phys. Fluids, 25 (2013), 062006.Google Scholar
[19] Filbet, F. and Pareschi, L., A numerical method for the accurate solution of the Fokker-Planck-Landau equation in the nonhomogeneous case, J. Comput. Phys., 179 (2002), 1–26.Google Scholar
[20] Lemou, M. and Mieussens, L., Implicit schemes for the Fokker-Planck-Landau equation, SIAM J. Sci. Comput., 27 (2005), 809–830.Google Scholar
[21] Jin, S. and Yan, B., A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation, J. Comput. Phys., 230 (2011), 6420–6437.Google Scholar
[22] Heinz, S., Molecular to fluid dynamics: the consequences of stochastic molecular motion, Phys. Rev. E, 70 (2004), 036308.CrossRefGoogle ScholarPubMed
[23] Gorji, M. H. and Jenny, P., Fokker-Planck-DSMC algorithm for simulations of rarefied gas flows, J. Comput. Phys., 287 (2015), 110–129.Google Scholar
[24] Gorji, M. H., Andric, N., and Jenny, P., Variance reduction for Fokker-Planck based particle Monte Carlo schemes, J. Comput. Phys., 295 (2015), 644–664.Google Scholar
[25] Xu, K., A gas-kinetic BGK scheme for the NavierCStokes equations and its connection with artificial dissipation and Godunov method, J. Comput. Phys., 171 (2001), 289–335.Google Scholar
[26] Fei, F. and Fan, J., A diffusive information preservation method for small Knudsen number flows, J. Comput. Phys., 243 (2013), 179–193.Google Scholar
[27] Fei, F., Fan, J., et al., Statistical simulation of decaying and forced homogeneous isotropic turbulence, in: Fan, J. (Ed.), 29th International Symposium on Rarefied Gas Dynamics, AIP Conference Proceedings 1628, 2014, 1356–1362.Google Scholar
[28] Zeng, D. D., Zhong, F. Q., and Fan, J., Molecular simulation of 3D turbulent channel flow, in: Fan, J. (Ed.), 29th International Symposium on Rarefied Gas Dynamics, AIP Conference Proceedings 1628, 2014, 1349–1355.Google Scholar
[29] Fei, F., Liu, Z.H., and Zheng, C. G., Statistical simulation of molecular diffusion effect on turbulent tetrad dispersion, Int. J. Heat Mass Transfer, 103 (2016), 87–98.Google Scholar
[30] Zhang, J., Zeng, D. D., and Fan, J., Statistical simulation of molecular diffusion effect on turbulent tetrad dispersion, Physica A, 411 (2014), 104–112.Google Scholar
[31] Chandrasekhar, S., Stochastic problems in physics and astronomy, Rev. Mod. Phys., 15 (1943), 1–89.Google Scholar
[32] Zwanzig, R., Nonequilibrium Statistical Mechanics, Oxford University Press, New York, 2001.Google Scholar
[33] Gardiner, C.W., Handbook of Stochastic Methods, Springer-Verlag, 1985.Google Scholar
[34] Shen, C., Rarefied Gas Dynamics: Fundamentals, Simulations and Micro Flows, Springer Press, Berlin, 2005.Google Scholar
[35] Helfand, E., Transport coefficients from dissipation in a canonical ensemble, Phys. Rev., 119 (1960), 1–9.Google Scholar
[36] Onskog, T. and Zhang, J., An accurate treatment of diffuse reflection boundary conditions for a stochastic particle Fokker-Planck algorithm with large time steps, Physica A, 440 (2015), 139–146.Google Scholar
[37] Myong, R. S., A full analytical solution for the force-driven compressible Poiseuille gas flow based on a nonlinear coupled constitutive relation, Phys. Fluid, 23 (2011), 012002.Google Scholar
[38] Sod, G. A., A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, J. Comput. Phys., 27 (1978), 1–31.Google Scholar
[39] Tiwari, S., Klar, A., and Hardt, S., A particle-particle hybrid method for kinetic and continuum equations, J. Comput. Phys., 228 (2009), 7109–7124.Google Scholar
[40] De Zarate, J.M.O. and Sengers, J. V., Hydrodynamic fluctuations in fluids and fluid mixtures, Elsevier, 2006.Google Scholar
[41] Hadjiconstantinou, N. G., Garcia, A. L., Bazant, M. Z., and He, G., Statistical error in particle simulations of hydrodynamic phenomena, J. Comput. Phys., 187 (2003), 274–297.Google Scholar
[42] Mathiaud, J. and Mieussens, L., A Fokker-Planck Model of the Boltzmann Equation with Correct Prandtl Number, J. Stat. Phys., 162 (2016), 397–414.Google Scholar
[43] Singh, S. K. and Ansumali, S., Fokker-Planck model of hydrodynamics, Phys. Rev. E, 91 (2015), 033303.Google Scholar
[44] Kannenberg, K. C. and Boyd, I. D., Three-dimensional Monte Carlo simulations of plume impingement, J. Thermophys. Heat Transf., 13 (1999), 226–235.Google Scholar
[45] Sharipov, F. and Kozak, D. V., Rarefied gas flow through a thin slit into vacuum simulated by the Monte Carlo method over the whole range of the Knudsen number, J. Vac. Sci. Technol., A, Vac. Surf. Films, 27 (2009), 479–484.Google Scholar