American Institute of Mathematical Sciences

Advanced
June  2013, 6(2): 291-315. doi: 10.3934/krm.2013.6.291

The moment guided Monte Carlo method for the Boltzmann equation

Giacomo Dimarco 1,

1. 

Institut de Mathématiques de Toulouse, UMR 5219, Université Paul Sabatier, 118, route de Narbonne 31062 TOULOUSE Cedex, France

Received  July 2012 Revised  October 2012 Published  February 2013

In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. Here, at the contrary to the previous work in which we considered the simplified BGK operator, we deal with the full Boltzmann operator. Moreover, we introduce an hybrid setting which permits to entirely remove the resolution of the kinetic equation in the limit of infinite number of collisions and to consider only the solution of the compressible Euler equation. This modification additionally reduce the statistical error with respect to our previous work and permits to perform simulations of non equilibrium gases using only a few number of particles. We show at the end of the paper several numerical tests which prove the efficiency and the low level of numerical noise of the method.
Keywords: Monte Carlo methods, hybrid methods, variance reduction, Boltzmann equation, moments closure, fluid equations..
Mathematics Subject Classification: Primary: 76P05, 65C20, 65C0.
Citation: Giacomo Dimarco. The moment guided Monte Carlo method for the Boltzmann equation. Kinetic & Related Models, 2013, 6 (2) : 291-315. doi: 10.3934/krm.2013.6.291
References:
[1]

H. Babovsky, On a simulation scheme for the Boltzmann equation, Math. Methods Appl. Sci., 8 (1986), 223-233. doi: 10.1002/mma.1670080114.  Google Scholar

[2]

G. A. Bird, "Molecular Gas Dynamics and Direct Simulation of Gas Flows," Oxford Engineering Science Series, 42, Oxford Science Publications, The Clarendon Press, Oxford University Press, New York, 1995.  Google Scholar

[3]

C. K. Birsdall and A. B. Langdon, "Plasma Physics Via Computer Simulation," Institute of Physics (IOP), Series in Plasma Physics, 2004. Google Scholar

[4]

J.-F. Bourgat, P. Le Tallec and M. D. Tidriri, Coupling Boltzmann and Navier-Stokes equations by friction, J. Comput. Phys., 127 (1996), 227-245. doi: 10.1006/jcph.1996.0172.  Google Scholar

[5]

J. Burt and I. Boyd, A hybrid particle approach for continuum and rarefied flow simulation, J. Comput. Phys., 228 (2009), 460-475. Google Scholar

[6]

R. E. Caflisch, Monte Carlo and Quasi-Monte Carlo Methods, in "Acta Numerica, 1998," Acta Numer., 7, Cambridge Univ. Press, Cambridge, (1998), 1-49. doi: 10.1017/S0962492900002804.  Google Scholar

[7]

R. E. Caflisch and L. Pareschi, Towards an hybrid method for rarefied gas dynamics, in "Transport in Transition Regimes" (Minneapolis, MN, 2000), IMA Vol. App. Math., 135, Springer, New York, (2004), 57-73. doi: 10.1007/978-1-4613-0017-5_3.  Google Scholar

[8]

C. Cercignani, "The Boltzmann Equation and its Applications," Applied Mathematical Sciences, 67, Springer-Verlag, New York, 1988. doi: 10.1007/978-1-4612-1039-9.  Google Scholar

[9]

A. Crestetto, N. Crouseilles and M. Lemou, Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles, to appear in KRM, 2012. Google Scholar

[10]

N. Crouseilles, P. Degond and M. Lemou, A hybrid kinetic-fluid model for solving the gas-dynamics Boltzmann BGK equation, J. Comput. Phys., 199 (2004), 776-808. doi: 10.1016/j.jcp.2004.03.007.  Google Scholar

[11]

P. Degond, G. Dimarco and L. Pareschi, The moment-guided Monte Carlo method, Int. J. Num. Meth. Fluids, 67 (2011), 189-213. doi: 10.1002/fld.2345.  Google Scholar

[12]

P. Degond, J.-G. Liu and L. Mieussens, Macroscopic fluid models with localized kinetic upscaling effects, Multiscale Model. Simul., 5 (2006), 940-979. doi: 10.1137/060651574.  Google Scholar

[13]

P. Degond and G. Dimarco, Fluid simulations with localized Boltzmann upscaling by direct Monte Carlo, J. Comp. Phys., 231 (2012), 2414-2437. doi: 10.1016/j.jcp.2011.11.030.  Google Scholar

[14]

P. Degond, G. Dimarco and L. Mieussens, A multiscale kinetic-fluid solver with dynamic localization of kinetic effects, J. Comp. Phys., 229 (2010), 4907-4933. doi: 10.1016/j.jcp.2010.03.009.  Google Scholar

[15]

G. Dimarco and L. Pareschi, Hybrid multiscale methods. I. Hyperbolic relaxation problems, Comm. Math. Sci., 1 (2006), 155-177.  Google Scholar

[16]

G. Dimarco and L. Pareschi, Fluid solver independent hybrid methods for multiscale kinetic equations, SIAM J. Sci. Comput., 32 (2010), 603-634. doi: 10.1137/080730585.  Google Scholar

[17]

G. Dimarco and L. Pareschi, Exponential Runge-Kutta methods for stiff kinetic equations, SIAM J. Num. Anal., 49 (2011), 2057-2077. doi: 10.1137/100811052.  Google Scholar

[18]

F. Filbet and S. Jin, A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources, J. Comp. Phys., 229 (2010), 7625-7648. doi: 10.1016/j.jcp.2010.06.017.  Google Scholar

[19]

F. Filbet and S. Jin, An asymptotic preserving scheme for the ES-BGK model of the Boltzmann equation, J. Sci. Comput., 46 (2011), 204-224. doi: 10.1007/s10915-010-9394-x.  Google Scholar

[20]

W. E and B. Engquist, The heterogeneous multiscale methods, Comm. Math. Sci., 1 (2003), 87-132.  Google Scholar

[21]

D. B. Hash and H. A. Hassan, Assessment of schemes for coupling Monte Carlo and Navier-Stokes solution methods, J. Thermophys. Heat Transf., 10 (1996), 242-249. Google Scholar

[22]

T. Homolle and N. Hadjiconstantinou, A low-variance deviational simulation Monte Carlo for the Boltzmann equation, J. Comp. Phys., 226 (2007), 2341-2358. doi: 10.1016/j.jcp.2007.07.006.  Google Scholar

[23]

T. Homolle and N. Hadjiconstantinou, Low-variance deviational simulation Monte Carlo, Phys. Fluids, 19 (2007), 041701. Google Scholar

[24]

S. Jin, Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations, SIAM J. Sci. Comput., 21 (1999), 441-454. doi: 10.1137/S1064827598334599.  Google Scholar

[25]

S. Jin and Z. P. Xin, Relaxation schemes for systems of conservation laws in arbitrary space dimensions, Comm. Pure Appl. Math., 48 (1995), 235-276. doi: 10.1002/cpa.3160480303.  Google Scholar

[26]

P. Le Tallec and F. Mallinger, Coupling Boltzmann and Navier-Stokes by half fluxes, J. Comput. Phys., 136 (1997), 51-67. doi: 10.1006/jcph.1997.5729.  Google Scholar

[27]

R. J. LeVeque, "Numerical Methods for Conservation Laws," Lectures in Mathematics, Birkhäuser Verlag, Basel, 1990. Google Scholar

[28]

S. Liu, "Monte Carlo Strategies in Scientific Computing," Springer, 2004. Google Scholar

[29]

K. Nanbu, Direct simulation scheme derived from the Boltzmann equation, J. Phys. Soc. Japan, 49 (1980), 2042-2049. Google Scholar

[30]

D. I. Pullin, Direct simulation methods for compressible inviscid ideal gas flow, J. Comput. Phys., 34 (1980), 231-244. Google Scholar

[31]

S. Tiwari, Coupling of the Boltzmann and Euler equations with automatic domain decomposition, J. Comput. Phys., 144 (1998), 710-726. doi: 10.1006/jcph.1998.6011.  Google Scholar

[32]

S. Tiwari, A. Klar and S. Hardt, A particle-particle hybrid method for kinetic and continuum equations, J. Comput. Phys., 228 (2009), 7109-7124. doi: 10.1016/j.jcp.2009.06.019.  Google Scholar

show all references

References:
[1]

H. Babovsky, On a simulation scheme for the Boltzmann equation, Math. Methods Appl. Sci., 8 (1986), 223-233. doi: 10.1002/mma.1670080114.  Google Scholar

[2]

G. A. Bird, "Molecular Gas Dynamics and Direct Simulation of Gas Flows," Oxford Engineering Science Series, 42, Oxford Science Publications, The Clarendon Press, Oxford University Press, New York, 1995.  Google Scholar

[3]

C. K. Birsdall and A. B. Langdon, "Plasma Physics Via Computer Simulation," Institute of Physics (IOP), Series in Plasma Physics, 2004. Google Scholar

[4]

J.-F. Bourgat, P. Le Tallec and M. D. Tidriri, Coupling Boltzmann and Navier-Stokes equations by friction, J. Comput. Phys., 127 (1996), 227-245. doi: 10.1006/jcph.1996.0172.  Google Scholar

[5]

J. Burt and I. Boyd, A hybrid particle approach for continuum and rarefied flow simulation, J. Comput. Phys., 228 (2009), 460-475. Google Scholar

[6]

R. E. Caflisch, Monte Carlo and Quasi-Monte Carlo Methods, in "Acta Numerica, 1998," Acta Numer., 7, Cambridge Univ. Press, Cambridge, (1998), 1-49. doi: 10.1017/S0962492900002804.  Google Scholar

[7]

R. E. Caflisch and L. Pareschi, Towards an hybrid method for rarefied gas dynamics, in "Transport in Transition Regimes" (Minneapolis, MN, 2000), IMA Vol. App. Math., 135, Springer, New York, (2004), 57-73. doi: 10.1007/978-1-4613-0017-5_3.  Google Scholar

[8]

C. Cercignani, "The Boltzmann Equation and its Applications," Applied Mathematical Sciences, 67, Springer-Verlag, New York, 1988. doi: 10.1007/978-1-4612-1039-9.  Google Scholar

[9]

A. Crestetto, N. Crouseilles and M. Lemou, Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles, to appear in KRM, 2012. Google Scholar

[10]

N. Crouseilles, P. Degond and M. Lemou, A hybrid kinetic-fluid model for solving the gas-dynamics Boltzmann BGK equation, J. Comput. Phys., 199 (2004), 776-808. doi: 10.1016/j.jcp.2004.03.007.  Google Scholar

[11]

P. Degond, G. Dimarco and L. Pareschi, The moment-guided Monte Carlo method, Int. J. Num. Meth. Fluids, 67 (2011), 189-213. doi: 10.1002/fld.2345.  Google Scholar

[12]

P. Degond, J.-G. Liu and L. Mieussens, Macroscopic fluid models with localized kinetic upscaling effects, Multiscale Model. Simul., 5 (2006), 940-979. doi: 10.1137/060651574.  Google Scholar

[13]

P. Degond and G. Dimarco, Fluid simulations with localized Boltzmann upscaling by direct Monte Carlo, J. Comp. Phys., 231 (2012), 2414-2437. doi: 10.1016/j.jcp.2011.11.030.  Google Scholar

[14]

P. Degond, G. Dimarco and L. Mieussens, A multiscale kinetic-fluid solver with dynamic localization of kinetic effects, J. Comp. Phys., 229 (2010), 4907-4933. doi: 10.1016/j.jcp.2010.03.009.  Google Scholar

[15]

G. Dimarco and L. Pareschi, Hybrid multiscale methods. I. Hyperbolic relaxation problems, Comm. Math. Sci., 1 (2006), 155-177.  Google Scholar

[16]

G. Dimarco and L. Pareschi, Fluid solver independent hybrid methods for multiscale kinetic equations, SIAM J. Sci. Comput., 32 (2010), 603-634. doi: 10.1137/080730585.  Google Scholar

[17]

G. Dimarco and L. Pareschi, Exponential Runge-Kutta methods for stiff kinetic equations, SIAM J. Num. Anal., 49 (2011), 2057-2077. doi: 10.1137/100811052.  Google Scholar

[18]

F. Filbet and S. Jin, A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources, J. Comp. Phys., 229 (2010), 7625-7648. doi: 10.1016/j.jcp.2010.06.017.  Google Scholar

[19]

F. Filbet and S. Jin, An asymptotic preserving scheme for the ES-BGK model of the Boltzmann equation, J. Sci. Comput., 46 (2011), 204-224. doi: 10.1007/s10915-010-9394-x.  Google Scholar

[20]

W. E and B. Engquist, The heterogeneous multiscale methods, Comm. Math. Sci., 1 (2003), 87-132.  Google Scholar

[21]

D. B. Hash and H. A. Hassan, Assessment of schemes for coupling Monte Carlo and Navier-Stokes solution methods, J. Thermophys. Heat Transf., 10 (1996), 242-249. Google Scholar

[22]

T. Homolle and N. Hadjiconstantinou, A low-variance deviational simulation Monte Carlo for the Boltzmann equation, J. Comp. Phys., 226 (2007), 2341-2358. doi: 10.1016/j.jcp.2007.07.006.  Google Scholar

[23]

T. Homolle and N. Hadjiconstantinou, Low-variance deviational simulation Monte Carlo, Phys. Fluids, 19 (2007), 041701. Google Scholar

[24]

S. Jin, Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations, SIAM J. Sci. Comput., 21 (1999), 441-454. doi: 10.1137/S1064827598334599.  Google Scholar

[25]

S. Jin and Z. P. Xin, Relaxation schemes for systems of conservation laws in arbitrary space dimensions, Comm. Pure Appl. Math., 48 (1995), 235-276. doi: 10.1002/cpa.3160480303.  Google Scholar

[26]

P. Le Tallec and F. Mallinger, Coupling Boltzmann and Navier-Stokes by half fluxes, J. Comput. Phys., 136 (1997), 51-67. doi: 10.1006/jcph.1997.5729.  Google Scholar

[27]

R. J. LeVeque, "Numerical Methods for Conservation Laws," Lectures in Mathematics, Birkhäuser Verlag, Basel, 1990. Google Scholar

[28]

S. Liu, "Monte Carlo Strategies in Scientific Computing," Springer, 2004. Google Scholar

[29]

K. Nanbu, Direct simulation scheme derived from the Boltzmann equation, J. Phys. Soc. Japan, 49 (1980), 2042-2049. Google Scholar

[30]

D. I. Pullin, Direct simulation methods for compressible inviscid ideal gas flow, J. Comput. Phys., 34 (1980), 231-244. Google Scholar

[31]

S. Tiwari, Coupling of the Boltzmann and Euler equations with automatic domain decomposition, J. Comput. Phys., 144 (1998), 710-726. doi: 10.1006/jcph.1998.6011.  Google Scholar

[32]

S. Tiwari, A. Klar and S. Hardt, A particle-particle hybrid method for kinetic and continuum equations, J. Comput. Phys., 228 (2009), 7109-7124. doi: 10.1016/j.jcp.2009.06.019.  Google Scholar

[1]

Zhiyan Ding, Qin Li. Constrained Ensemble Langevin Monte Carlo. Foundations of Data Science, 2021  doi: 10.3934/fods.2021034
[2]

Thomas Chen, Ryan Denlinger, Nataša Pavlović. Moments and regularity for a Boltzmann equation via Wigner transform. Discrete & Continuous Dynamical Systems, 2019, 39 (9) : 4979-5015. doi: 10.3934/dcds.2019204
[3]

Milana Pavić-Čolić, Maja Tasković. Propagation of stretched exponential moments for the Kac equation and Boltzmann equation with Maxwell molecules. Kinetic & Related Models, 2018, 11 (3) : 597-613. doi: 10.3934/krm.2018025
[4]

M. Sumon Hossain, M. Monir Uddin. Iterative methods for solving large sparse Lyapunov equations and application to model reduction of index 1 differential-algebraic-equations. Numerical Algebra, Control & Optimization, 2019, 9 (2) : 173-186. doi: 10.3934/naco.2019013
[5]

Guillaume Bal, Ian Langmore, Youssef Marzouk. Bayesian inverse problems with Monte Carlo forward models. Inverse Problems & Imaging, 2013, 7 (1) : 81-105. doi: 10.3934/ipi.2013.7.81
[6]

Ajay Jasra, Kody J. H. Law, Yaxian Xu. Markov chain simulation for multilevel Monte Carlo. Foundations of Data Science, 2021, 3 (1) : 27-47. doi: 10.3934/fods.2021004
[7]

Theodore Papamarkou, Alexey Lindo, Eric B. Ford. Geometric adaptive Monte Carlo in random environment. Foundations of Data Science, 2021, 3 (2) : 201-224. doi: 10.3934/fods.2021014
[8]

Yves Bourgault, Damien Broizat, Pierre-Emmanuel Jabin. Convergence rate for the method of moments with linear closure relations. Kinetic & Related Models, 2015, 8 (1) : 1-27. doi: 10.3934/krm.2015.8.1
[9]

Wenxiong Chen, Shijie Qi. Direct methods on fractional equations. Discrete & Continuous Dynamical Systems, 2019, 39 (3) : 1269-1310. doi: 10.3934/dcds.2019055
[10]

Jérôme Droniou. Remarks on discretizations of convection terms in Hybrid mimetic mixed methods. Networks & Heterogeneous Media, 2010, 5 (3) : 545-563. doi: 10.3934/nhm.2010.5.545
[11]

Michael B. Giles, Kristian Debrabant, Andreas Rössler. Analysis of multilevel Monte Carlo path simulation using the Milstein discretisation. Discrete & Continuous Dynamical Systems - B, 2019, 24 (8) : 3881-3903. doi: 10.3934/dcdsb.2018335
[12]

Jiakou Wang, Margaret J. Slattery, Meghan Henty Hoskins, Shile Liang, Cheng Dong, Qiang Du. Monte carlo simulation of heterotypic cell aggregation in nonlinear shear flow. Mathematical Biosciences & Engineering, 2006, 3 (4) : 683-696. doi: 10.3934/mbe.2006.3.683
[13]

B. S. Goh, W. J. Leong, Z. Siri. Partial Newton methods for a system of equations. Numerical Algebra, Control & Optimization, 2013, 3 (3) : 463-469. doi: 10.3934/naco.2013.3.463
[14]

Tengteng Yu, Xin-Wei Liu, Yu-Hong Dai, Jie Sun. Variable metric proximal stochastic variance reduced gradient methods for nonconvex nonsmooth optimization. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021084
[15]

Samira Shahsavari, Saeed Ketabchi. The proximal methods for solving absolute value equation. Numerical Algebra, Control & Optimization, 2021, 11 (3) : 449-460. doi: 10.3934/naco.2020037
[16]

Andrea Malchiodi. Topological methods for an elliptic equation with exponential nonlinearities. Discrete & Continuous Dynamical Systems, 2008, 21 (1) : 277-294. doi: 10.3934/dcds.2008.21.277
[17]

Noboru Okazawa, Toshiyuki Suzuki, Tomomi Yokota. Energy methods for abstract nonlinear Schrödinger equations. Evolution Equations & Control Theory, 2012, 1 (2) : 337-354. doi: 10.3934/eect.2012.1.337
[18]

Lijian Jiang, Yalchin Efendiev, Victor Ginting. Multiscale methods for parabolic equations with continuum spatial scales. Discrete & Continuous Dynamical Systems - B, 2007, 8 (4) : 833-859. doi: 10.3934/dcdsb.2007.8.833
[19]

Hui Liang, Hermann Brunner. Collocation methods for differential equations with piecewise linear delays. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1839-1857. doi: 10.3934/cpaa.2012.11.1839
[20]

Zalman Balanov, Carlos García-Azpeitia, Wieslaw Krawcewicz. On variational and topological methods in nonlinear difference equations. Communications on Pure & Applied Analysis, 2018, 17 (6) : 2813-2844. doi: 10.3934/cpaa.2018133

2020 Impact Factor: 1.432

Tools

Metrics

  • PDF downloads (136)
  • HTML views (0)
  • Cited by (8)

Other articles
by authors

[Back to Top]

Copyright © 2022 American Institute of Mathematical Sciences | Privacy Policy