Kinetic modeling of multiphase fow based on simplified Enskog equation  

作　　者：Yu-Dong Zhang Ai-Guo Xu Jing-Jiang Qiu Hong-Tao Wei Zung Hang Wei 

机构地区：[1]School of Mechanics and Safety Engineering,Zhengzhou University,Zhengzhou 450001,China [2]Laboratory of Computational Physics,Institute of Applied Physics and Computational Mathematics,P.O.Bor 8009-26,Beijing 100088,China [3]Center for Applied Physics and Technology,MOE Key Center for High Energy Density Physics Simulations,College of Engineering,Peking Universily,Beifing 100871,China

出　　处：《Frontiers of physics》2020年第6期133-145,共13页物理学前沿（英文版）

基　　金：Many thanks to the anonymous referees for thelr helpful comments and suggestions.This work was supported by the China Postdoctoral Science Foundation Grant No.2019M662521,the National Natural Science Foundation of China under Grant No.11772064,CAEP Foundation(Grant No.CX201903);the Opening Project of State Key Laboratory of Explosion Science and Technology(Beiing Institute of Technology)under Grant No.KFJJ19-01 M.

摘　　要：A new kinetic model for multiphase flow was presented under the framework of the discrete Boltzmann method(DBM).Significantly different from the previous DBM,a bottom-p approach was adopted in this model.The effects of molecular size and repulsion potential were described by the Enskog collision model;the attraction potential was obtained through the mean field approximation method.The molecular interactions,which result in the non-ideal equation of state and surface tension,were directly introduced as an external force term.Several typical benchmark problems,including Couette flow,two-phase coexistence curve,the Laplace law,phase separation,and the collision of two droplets,were simulated to verify the model.Especially,for two types of droplet collisions,the strengths of two non-equilibrium ffcts,D*2 and D*3,defined through the second and third order non-conserved kinetic moments of(f-feq),are comparatively investigated,where f(feq)is the(equilibrium)distribution function.It is interesting to find that during the collision process,D*2is always significantly larger than D*3,D*2 can be used to identify the different stages of the collision process and to distinguish different types of collisions.The modeling method can be directly extended to a higher.order model for the case where the non-equilibrium pffect is strong,and the linear constitutive law of viscous stress is no longer valid.

关 键 词：multiphase flow discrete Boltzmann method Enskog equation non-equilibrium charac-teristics 

分 类 号：O35[理学—流体力学]