Extension of lattice Boltzmann flux solver for simulation of 3D viscous compressible flows
Extension of lattice Boltzmann flux solver for simulation of 3D viscous compressible flows
复制标题
用于模拟 3D 粘性可压缩流的晶格玻尔兹曼通量解算器的扩展
DOI：
10.1016/j.camwa.2016.03.027
复制
发表时间：
2016-05
影响因子：
2.9
通讯作者：
J. Wu
中科院分区：
文献类型：
--
作者：
L. M. Yang;C. Shu;J. Wu
Shu 和他的同事（Yang et al., 2012, 2013; Ji et al., 2009; Shu et al., 2014）提出的用于模拟非粘性可压缩流的格子玻尔兹曼通量求解器（LBFS）在这项工作中被扩展到模拟 3D 粘性可压缩流。在求解器中，通过将非自由参数D1Q4模型应用于黎曼问题，通过一维格子玻尔兹曼解的局部重建来评估细胞界面处的无粘通量，而粘性通量则通过传统的光滑函数近似来评估。在现有的LBFS中（Yang et al., 2012, 2013; Ji et al., 2009; Shu et al., 2014），直接使用从相邻点流出的单元界面处的分布函数来计算无粘通量，其中包含用于模拟粘性流的超量数值耗散。在目前的工作中，我们从 Chapman-Enskog 分析（Guo 和 Shu，2013）开始，考虑细胞界面处分布函数的平衡部分和非平衡部分。众所周知，无粘通量可以完全由平衡部分确定，非平衡部分可以视为无粘通量计算的数值耗散。通过引入范围从 0 到 1 的开关函数来控制数值耗散，消除了现有 LBFS 的缺点。在边界层等光滑区域，开关函数取接近于零的值，而在强激波周围，开关函数则趋于1。通过具有复杂几何形状的测试用例，已经证明本求解器可以很好地模拟 3D 粘性可压缩流。
The lattice Boltzmann flux solver (LBFS), which was presented by Shu and his coworkers (Yang et al., 2012, 2013; Ji et al., 2009; Shu et al., 2014) for simulation of inviscid compressible flows, is extended to simulate 3D viscous compressible flows in this work. In the solver, the inviscid flux at the cell interface is evaluated by local reconstruction of one-dimensional lattice Boltzmann solution through the application of non-free parameter D1Q4 model to the Riemann problem, while the viscous flux is evaluated by conventional smooth function approximation. In the existing LBFS (Yang et al., 2012, 2013; Ji et al., 2009; Shu et al., 2014), the distribution functions at the cell interface streamed from neighboring points are directly used to compute the inviscid flux, which contains superabundant numerical dissipation for simulation of viscous flows. In the present work, we start from the Chapman–Enskog analysis (Guo and Shu, 2013) and consider both the equilibrium part and non-equilibrium part of the distribution function at the cell interface. It is well known that the inviscid flux can be fully determined by the equilibrium part and the non-equilibrium part can be viewed as numerical dissipation for the calculation of inviscid flux. The drawback of the existing LBFS is removed by introducing a switch function which ranges from 0 to 1 in order to control the numerical dissipation. In the smooth region such as in boundary layer, the switch function takes a value close to zero, while around the strong shock wave, it tends to one. Through test cases with complex geometry, it has been demonstrated that the present solver can work very well for simulation of 3D viscous compressible flows.
登录
查看更多内容
影响因子：
2.8
作者：
M. Stiebler;J. Tölke;M. Krafczyk
通讯作者：
M. Stiebler;J. Tölke;M. Krafczyk
影响因子：
6.5
作者：
G. D. V. Albada;B. Leer;W. Roberts
通讯作者：
G. D. V. Albada;B. Leer;W. Roberts
影响因子：
1.4
作者：
L. M. Yang;C. Shu;Jie Wu
通讯作者：
L. M. Yang;C. Shu;Jie Wu
影响因子：
2.4
作者：
Kataoka, T;Tsutahara, M
通讯作者：
Tsutahara, M
DOI：
10.1103/physreve.69.056702
发表时间：
2004-05
期刊：
Physical review. E, Statistical, nonlinear, and soft matter physics
影响因子：
--
作者：
T. Kataoka;M. Tsutahara
通讯作者：
T. Kataoka;M. Tsutahara