Boosting the convergence of low-variance DSMC by GSIS

Advances in Aerodynamics

Research

Qi Li¹  Liyan Luo¹  Lei Wu¹ 

[1] Department of Mechanics and Aerospace Engineering, Southern University of Science and Technology, 518055, Shenzhen, China;

关键词: Rarefied gas dynamics;    Direct simulation Monte Carlo;    General synthetic iterative scheme;    Multiscale simulation;   

DOI  :  10.1186/s42774-023-00138-0

 received in 2022-10-24, accepted in 2022-12-28,  发布年份 2022

来源: Springer

【 摘 要 】

The low-variance direct simulation Monte Carlo (LVDSMC) is a powerful method to simulate low-speed rarefied gas flows. However, in the near-continuum flow regime, due to limitations on the time step and spatial cell size, it takes plenty of time to find the steady-state solution. Here we remove these deficiencies by coupling the LVDSMC with the general synthetic iterative scheme (GSIS) which permits the simulation at the hydrodynamic scale rather than the much smaller kinetic scale. As a proof of concept, we propose the stochastic-deterministic coupling method based on the Bhatnagar-Gross-Krook kinetic model. First, macroscopic synthetic equations are derived exactly from the kinetic equation, which not only contain the Navier-Stokes-Fourier constitutive relation, but also encompass the higher-order terms describing the rarefaction effects. Then, the high-order terms are extracted from LVDSMC and fed into synthetic equations to predict the macroscopic properties which are closer to the steady-state solution than LVDSMC. Finally, the state of simulation particles in LVDSMC is updated to reflect the change of macroscopic properties. As a result, the convergence to steady state is greatly accelerated, and the restrictions on cell size and the time step are removed. We conduct the Fourier stability analysis and simulate several canonical rarefied gas flows to demonstrate the advantages of LVDSMC-GSIS: when the Knudsen number is lower than 0.1, it can use the grid size about 10 times larger than that in traditional DSMC, and it can reduce the computational cost by two orders of magnitude in the flow regime.

【 授权许可】

CC BY   
© The Author(s) 2023

【 图 表 】

42774_2023_138_Article_IEq297.gif

42774_2023_138_Article_IEq296.gif

42774_2023_138_Article_IEq295.gif

42774_2023_138_Article_IEq294.gif

42774_2023_138_Article_IEq293.gif

42774_2023_138_Article_IEq292.gif

42774_2023_138_Article_IEq291.gif

42774_2023_138_Article_IEq290.gif

42774_2023_138_Article_IEq289.gif

42774_2023_138_Article_IEq288.gif

42774_2023_138_Article_IEq287.gif

42774_2023_138_Article_IEq286.gif

42774_2023_138_Article_IEq285.gif

42774_2023_138_Article_IEq284.gif

42774_2023_138_Article_IEq283.gif

42774_2023_138_Article_IEq282.gif

42774_2023_138_Article_IEq281.gif

42774_2023_138_Article_IEq280.gif

42774_2023_138_Article_IEq279.gif

42774_2023_138_Article_IEq278.gif

42774_2023_138_Article_IEq277.gif

42774_2023_138_Article_IEq276.gif

42774_2023_138_Article_IEq275.gif

42774_2023_138_Article_IEq274.gif

42774_2023_138_Article_IEq273.gif

42774_2023_138_Article_IEq272.gif

42774_2023_138_Article_IEq271.gif

42774_2023_138_Article_IEq270.gif

42774_2023_138_Article_IEq269.gif

42774_2023_138_Article_IEq268.gif

42774_2023_138_Article_IEq267.gif

42774_2023_138_Article_IEq266.gif

42774_2023_138_Article_IEq265.gif

42774_2023_138_Article_IEq264.gif

42774_2023_138_Article_IEq263.gif

42774_2023_138_Article_IEq262.gif

42774_2023_138_Article_IEq261.gif

42774_2023_138_Article_IEq260.gif

42774_2023_138_Article_IEq259.gif

42774_2023_138_Article_IEq258.gif

42774_2023_138_Article_IEq257.gif

42774_2023_138_Article_IEq256.gif

42774_2023_138_Article_IEq255.gif

13753_2022_461_Article_IEq3.gif

42774_2023_138_Article_IEq253.gif

42774_2023_138_Article_IEq252.gif

42774_2023_138_Article_IEq251.gif

42774_2023_138_Article_IEq250.gif

42774_2023_138_Article_IEq249.gif

42774_2023_138_Article_IEq248.gif

42774_2023_138_Article_IEq247.gif

42774_2023_138_Article_IEq246.gif

42774_2023_138_Article_IEq245.gif

42774_2023_138_Article_IEq244.gif

42774_2023_138_Article_IEq243.gif

42774_2023_138_Article_IEq242.gif

42774_2023_138_Article_IEq241.gif

42774_2023_138_Article_IEq240.gif

42774_2023_138_Article_IEq239.gif

42774_2023_138_Article_IEq238.gif

42774_2023_138_Article_IEq237.gif

42774_2023_138_Article_IEq236.gif

13753_2022_461_Article_IEq1.gif

42774_2023_138_Article_IEq234.gif

42774_2023_138_Article_IEq233.gif

42774_2023_138_Article_IEq232.gif

42774_2023_138_Article_IEq231.gif

42774_2023_138_Article_IEq230.gif

42774_2023_138_Article_IEq229.gif

Fig. 6

42774_2023_138_Article_IEq227.gif

42774_2023_138_Article_IEq226.gif

Fig. 5

42774_2023_138_Article_IEq224.gif

42774_2023_138_Article_IEq223.gif

42774_2023_138_Article_IEq222.gif

Fig. 4

42774_2023_138_Article_IEq220.gif

42774_2023_138_Article_IEq219.gif

42774_2023_138_Article_IEq218.gif

42774_2023_138_Article_IEq217.gif

42774_2023_138_Article_IEq216.gif

42774_2023_138_Article_IEq215.gif

42774_2023_138_Article_IEq214.gif

Fig. 2

Fig. 2

42774_2023_138_Article_IEq211.gif

42774_2023_138_Article_IEq210.gif

Fig. 1

42774_2023_138_Article_IEq208.gif

42774_2023_138_Article_IEq207.gif

Fig. 1

42774_2023_138_Article_IEq205.gif

42774_2023_138_Article_IEq204.gif

42774_2023_138_Article_IEq203.gif

42774_2023_138_Article_IEq202.gif

42774_2023_138_Article_IEq201.gif

42774_2023_138_Article_IEq200.gif

42774_2023_138_Article_IEq199.gif

42774_2023_138_Article_IEq198.gif

42774_2023_138_Article_IEq197.gif

42774_2023_138_Article_IEq196.gif

42774_2023_138_Article_IEq195.gif

42774_2023_138_Article_IEq194.gif

42774_2023_138_Article_IEq193.gif

42774_2023_138_Article_IEq192.gif

42774_2023_138_Article_IEq191.gif

42774_2023_138_Article_IEq190.gif

42774_2023_138_Article_IEq189.gif

42774_2023_138_Article_IEq188.gif

13753_2023_465_Article_IEq5.gif

13753_2023_465_Article_IEq4.gif

13753_2023_465_Article_IEq3.gif

13753_2023_465_Article_IEq2.gif

13753_2023_465_Article_IEq1.gif

Fig. 15

42774_2023_138_Article_IEq179.gif

Fig. 14

Fig. 2

42774_2023_138_Article_IEq176.gif

42774_2023_138_Article_IEq175.gif

Fig. 1

42774_2023_138_Article_IEq173.gif

42774_2023_138_Article_IEq172.gif

42774_2023_138_Article_IEq171.gif

42774_2023_138_Article_IEq170.gif

Fig. 13

42774_2023_138_Article_IEq168.gif

42774_2023_138_Article_IEq167.gif

Fig. 12

42774_2023_138_Article_IEq165.gif

42774_2023_138_Article_IEq164.gif

42774_2023_138_Article_IEq163.gif

42774_2023_138_Article_IEq162.gif

42774_2023_138_Article_IEq161.gif

42774_2023_138_Article_IEq160.gif

42774_2023_138_Article_IEq159.gif

42774_2023_138_Article_IEq158.gif

42774_2023_138_Article_IEq157.gif

42774_2023_138_Article_IEq156.gif

42774_2023_138_Article_IEq155.gif

42774_2023_138_Article_IEq154.gif

42774_2023_138_Article_IEq153.gif

42774_2023_138_Article_IEq152.gif

42774_2023_138_Article_IEq151.gif

42774_2023_138_Article_IEq150.gif

42774_2023_138_Article_IEq149.gif

42774_2023_138_Article_IEq148.gif

42774_2023_138_Article_IEq147.gif

42774_2023_138_Article_IEq146.gif

42774_2023_138_Article_IEq145.gif

42774_2023_138_Article_IEq144.gif

42774_2023_138_Article_IEq143.gif

42774_2023_138_Article_IEq142.gif

42774_2023_138_Article_IEq141.gif

42774_2023_138_Article_IEq140.gif

42774_2023_138_Article_IEq139.gif

42774_2023_138_Article_IEq138.gif

42774_2023_138_Article_IEq137.gif

42774_2023_138_Article_IEq136.gif

42774_2023_138_Article_IEq135.gif

Fig. 11

42774_2023_138_Article_IEq131.gif

Fig. 2

Fig. 1

42774_2023_138_Article_IEq128.gif

42774_2023_138_Article_IEq127.gif

42774_2023_138_Article_IEq126.gif

Fig. 10

42774_2023_138_Article_IEq124.gif

42774_2023_138_Article_IEq123.gif

42774_2023_138_Article_IEq122.gif

42774_2023_138_Article_IEq121.gif

42774_2023_138_Article_IEq120.gif

42774_2023_138_Article_IEq119.gif

Fig. 8

42774_2023_138_Article_IEq117.gif

42774_2023_138_Article_IEq116.gif

42774_2023_138_Article_IEq115.gif

42774_2023_138_Article_IEq114.gif

42774_2023_138_Article_IEq113.gif

42774_2023_138_Article_IEq112.gif

42774_2023_138_Article_IEq111.gif

42774_2023_138_Article_IEq110.gif

Fig. 7

42774_2023_138_Article_IEq108.gif

42774_2023_138_Article_IEq107.gif

42774_2023_138_Article_IEq106.gif

42774_2023_138_Article_IEq105.gif

42774_2023_138_Article_IEq104.gif

42774_2023_138_Article_IEq103.gif

42774_2023_138_Article_IEq102.gif

42774_2023_138_Article_IEq101.gif

42774_2023_138_Article_IEq100.gif

Fig. 6

42774_2023_138_Article_IEq98.gif

42774_2023_138_Article_IEq97.gif

Fig. 5

42774_2023_138_Article_IEq95.gif

42774_2023_138_Article_IEq94.gif

42774_2023_138_Article_IEq93.gif

42774_2023_138_Article_IEq92.gif

42774_2023_138_Article_IEq91.gif

42774_2023_138_Article_IEq90.gif

42774_2023_138_Article_IEq89.gif

42774_2023_138_Article_IEq88.gif

42774_2023_138_Article_IEq87.gif

42774_2023_138_Article_IEq86.gif

Fig. 1

42774_2023_138_Article_IEq84.gif

42774_2023_138_Article_IEq83.gif

42774_2023_138_Article_IEq82.gif

42774_2023_138_Article_IEq81.gif

42774_2023_138_Article_IEq80.gif

42774_2023_138_Article_IEq79.gif

42774_2023_138_Article_IEq78.gif

42774_2023_138_Article_IEq77.gif

Fig. 3

42774_2023_138_Article_IEq75.gif

42774_2023_138_Article_IEq74.gif

42774_2023_138_Article_IEq73.gif

Fig. 2

42774_2023_138_Article_IEq71.gif

42774_2023_138_Article_IEq70.gif

42774_2023_138_Article_IEq69.gif

42774_2023_138_Article_IEq68.gif

42774_2023_138_Article_IEq67.gif

42774_2023_138_Article_IEq66.gif

Fig. 1

42774_2023_138_Article_IEq64.gif

42774_2023_138_Article_IEq63.gif

42774_2023_138_Article_IEq62.gif

42774_2023_138_Article_IEq61.gif

42774_2023_138_Article_IEq60.gif

42774_2023_138_Article_IEq59.gif

40100_2022_240_Article_IEq37.gif

40100_2022_240_Article_IEq36.gif

40100_2022_240_Article_IEq35.gif

40100_2022_240_Article_IEq34.gif

40100_2022_240_Article_IEq33.gif

40100_2022_240_Article_IEq32.gif

40100_2022_240_Article_IEq31.gif

40100_2022_240_Article_IEq30.gif

40100_2022_240_Article_IEq29.gif

40100_2022_240_Article_IEq28.gif

40100_2022_240_Article_IEq27.gif

40100_2022_240_Article_IEq25.gif

40100_2022_240_Article_IEq24.gif

42774_2023_138_Article_IEq45.gif

40100_2022_240_Article_IEq23.gif

40100_2022_240_Article_IEq22.gif

40100_2022_240_Article_IEq21.gif

40100_2022_240_Article_IEq20.gif

40100_2022_240_Article_IEq19.gif

40100_2022_240_Article_IEq18.gif

40100_2022_240_Article_IEq17.gif

40100_2022_240_Article_IEq16.gif

40100_2022_240_Article_IEq15.gif

42774_2023_138_Article_IEq35.gif

40100_2022_240_Article_IEq14.gif

40100_2022_240_Article_IEq13.gif

40100_2022_240_Article_IEq12.gif

40100_2022_240_Article_IEq11.gif

40100_2022_240_Article_IEq10.gif

40100_2022_240_Article_IEq9.gif

40100_2022_240_Article_IEq8.gif

40100_2022_240_Article_IEq7.gif

40100_2022_240_Article_IEq6.gif

40100_2022_240_Article_IEq5.gif

40100_2022_240_Article_IEq4.gif

40100_2022_240_Article_IEq3.gif

42774_2023_138_Article_IEq22.gif

40100_2022_240_Article_IEq2.gif

42774_2023_138_Article_IEq20.gif

42774_2023_138_Article_IEq19.gif

42774_2023_138_Article_IEq18.gif

42774_2023_138_Article_IEq17.gif

40100_2022_240_Article_IEq1.gif

42774_2023_138_Article_IEq15.gif

42774_2023_138_Article_IEq14.gif

42774_2023_138_Article_IEq13.gif

42774_2023_138_Article_IEq12.gif

42774_2023_138_Article_IEq11.gif

42774_2023_138_Article_IEq10.gif

42774_2023_138_Article_IEq9.gif

42774_2023_138_Article_IEq8.gif

42774_2023_138_Article_IEq7.gif

42774_2023_138_Article_IEq6.gif

Fig. 2

42774_2023_138_Article_IEq4.gif

42774_2023_138_Article_IEq2.gif

42774_2023_138_Article_IEq1.gif

Fig. 1

Fig. 10

Fig. 9

Fig. 8

Fig. 7

Fig. 1

Fig. 5

Fig. 4

Fig. 3

Fig. 2

Fig. 1