Adaptive wave-particle decomposition in UGKWP method for high-speed flow simulations

Advances in Aerodynamics

Research

Yufeng Wei¹  Kun Xu²  Junzhe Cao³  Xing Ji⁴ 

[1] Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong, China;Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong, China;Department of Mechanical and Aerospace Engineering, Hong Kong University of Science and Technology, Hong Kong, China;HKUST Shenzhen Research Institute, 518057, Shenzhen, China;School of Aeronautics, Northwestern Polytechnical University, 710072, Xi’an, China;Shaanxi Key Laboratory of Environment and Control for Flight Vehicle, Xi’an Jiaotong University, 710049, Xi’an, China;

关键词: Adaptive wave-particle decomposition;    Multiscale modeling;    Acceleration method;    Non-equilibrium transport;   

DOI  :  10.1186/s42774-023-00156-y

 received in 2023-04-03, accepted in 2023-07-23,  发布年份 2023

来源: Springer

【 摘 要 】

With wave-particle decomposition, a unified gas-kinetic wave-particle (UGKWP) method has been developed for multiscale flow simulations. With the variation of the cell Knudsen number, the UGKWP method captures the transport process in all flow regimes without the kinetic solver’s constraint on the numerical mesh size and time step being determined by the kinetic particle mean free path and particle collision time. In the current UGKWP method, the cell Knudsen number, which is defined as the ratio of particle collision time to numerical time step, is used to distribute the components in the wave-particle decomposition. The adaptation of particles in the UGKWP method is mainly for the capturing of the non-equilibrium transport. In this aspect, the cell Knudsen number alone is not enough to identify the non-equilibrium state. For example, in the equilibrium flow regime with a Maxwellian distribution function, even at a large cell Knudsen number, the flow evolution can be still modelled by the Navier-Stokes solver. More specifically, in the near space environment both the hypersonic flow around a space vehicle and the plume flow from a satellite nozzle will encounter a far field rarefied equilibrium flow in a large computational domain. In the background dilute equilibrium region, the large particle collision time and a uniform small numerical time step can result in a large local cell Knudsen number and make the UGKWP method track a huge number of particles for the far field background flow in the original approach. But, in this region the analytical wave representation can be legitimately used in the UGKWP method to capture the nearly equilibrium flow evolution. Therefore, to further improve the efficiency of the UGKWP method for multiscale flow simulations, an adaptive UGKWP (AUGKWP) method is developed with the introduction of an additional local flow variable gradient-dependent Knudsen number. As a result, the wave-particle decomposition in the UGKWP method is determined by both the cell and gradient Knudsen numbers, and the use of particles in the UGKWP method is solely to capture the non-equilibrium flow transport. The current AUGKWP method becomes much more efficient than the previous one with the cell Knudsen number only in the determination of wave-particle composition. Many numerical tests, including Sod shock tube, normal shock structure, hypersonic flow around cylinder, flow around reentry capsule, and an unsteady nozzle plume flow, have been conducted to validate the accuracy and efficiency of the AUGKWP method. Compared with the original UGKWP method, the AUGKWP method achieves the same accuracy, but has advantages in memory reduction and computational efficiency in the simulation for flows with the co-existing of multiple regimes.

【 授权许可】

CC BY   
© Chinese Aerodynamics Research Society 2023

【 图 表 】

40517_2023_273_Article_IEq2.gif

Fig. 2

42004_2023_1025_Article_IEq7.gif

Fig. 9

Fig. 9

Fig. 4

Fig. 3

41512_2023_158_Article_IEq26.gif

Fig. 2

Fig. 2

Fig. 1

41512_2023_158_Article_IEq20.gif

Fig. 8

Fig. 2

Fig. 6

41512_2023_158_Article_IEq9.gif

Fig. 8

12951_2017_323_Article_IEq1.gif

41512_2023_158_Article_IEq1.gif

Fig. 5

Fig. 1

Fig. 7

Fig. 4

Fig. 3

Fig. 6

Fig. 1

12951_2015_155_Article_IEq4.gif

12936_2017_2014_Article_IEq78.gif

Fig. 4

Fig. 2

Fig. 9

12867_2016_60_Article_IEq2.gif

Fig. 2

Fig. 3

Fig. 5

Fig. 2

Fig. 8

Fig. 1

Fig. 4

Fig. 4

Fig. 1

12951_2015_155_Article_IEq1.gif

12936_2015_894_Article_IEq86.gif

Fig. 7

Fig. 2

Fig. 2

Fig. 2

12936_2017_1885_Article_IEq1.gif

Fig. 3

Fig. 6

Fig. 2

Fig. 4

Fig. 3

12888_2023_5299_Article_IEq6.gif

12888_2023_5299_Article_IEq5.gif

12888_2023_5299_Article_IEq4.gif

Fig. 6

12888_2023_5299_Article_IEq1.gif

Fig. 1

12936_2017_2051_Article_IEq71.gif

40708_2023_205_Article_IEq12.gif

Fig. 1

Fig. 3

12944_2017_533_Article_IEq2.gif

Fig. 2

Fig. 6

Fig. 1

Fig. 1

12864_2023_9693_Article_IEq1.gif

Fig. 3

Fig. 3

Fig. 5

Fig. 4

Fig. 4

Fig. 2

Fig. 3

Fig. 2

Fig. 4

12936_2017_2014_Article_IEq61.gif

Fig. 3

Fig. 3

12936_2015_966_Article_IEq12.gif

Fig. 4

12944_2023_1932_Article_IEq4.gif

Fig. 1

Fig. 5

12944_2023_1932_Article_IEq1.gif

Fig. 3

Fig. 1

Fig. 5

Fig. 4

Fig. 2

Fig. 5

Fig. 3

Fig. 2

Fig. 1

Fig. 2

12864_2017_4133_Article_IEq35.gif

Fig. 3

Fig. 3

Fig. 2

Fig. 1

Fig. 2

Fig. 5

12951_2016_171_Article_IEq6.gif

12951_2016_225_Article_IEq2.gif

Fig. 3

Fig. 6

Fig. 1

Fig. 6

12951_2015_155_Article_IEq88.gif

Fig. 1

Fig. 4

12936_2023_4742_Article_IEq70.gif

Fig. 6

Fig. 4

12951_2016_171_Article_IEq3.gif

Fig. 10

13731_2023_319_Article_IEq4.gif

13731_2023_319_Article_IEq3.gif

Fig. 4

12951_2015_155_Article_IEq86.gif

Fig. 4

Fig. 9

Fig. 6

Fig. 1

Fig. 2

Fig. 1

12951_2016_171_Article_IEq1.gif

12936_2023_4742_Article_IEq36.gif

Fig. 5

Fig. 5

Fig. 3

Fig. 2

Fig. 3

Fig. 5

12951_2017_252_Article_IEq1.gif

Fig. 4

Fig. 2

12936_2023_4742_Article_IEq20.gif

Fig. 1

12936_2023_4742_Article_IEq18.gif

Fig. 2

Fig. 5

Fig. 1

Fig. 4

Fig. 3

12888_2023_5283_Article_IEq1.gif

Fig. 1

Fig. 2

Fig. 1

Fig. 4

Fig. 2

Fig. 7

Fig. 5

Fig. 2

Fig. 1

12951_2016_246_Article_IEq15.gif

12951_2016_246_Article_IEq13.gif

Fig. 1

Fig. 9

12951_2017_255_Article_IEq52.gif

12951_2017_255_Article_IEq51.gif

Fig. 2

Fig. 1

Fig. 1

12951_2016_223_Article_IEq1.gif

12936_2017_1904_Article_IEq1.gif