Accuracy of non-Newtonian Lattice Boltzmann simulations
详细信息 查看全文
- 作者：Daniel Conrad ; daniel.conrad@mv.uni-kl.de" class="auth_mail" title="E-mail the corresponding author ; Andreas Schneider ; Martin Bö ; hle
- 关键词：Lattice Boltzmann Method ; Accuracy ; Non-Newtonian
- 刊名：Journal of Computational Physics
- 出版年：2015
- 出版时间：15 November 2015
- 年：2015
- 卷：301
- 期：Complete
- 页码：218-229
- 全文大小：862 K
文摘
This work deals with the accuracy of non-Newtonian Lattice Boltzmann simulations. Previous work for Newtonian fluids indicate that, depending on the numerical value of the dimensionless collision frequency 惟, additional artificial viscosity is introduced, which negatively influences the accuracy. Since the non-Newtonian fluid behavior is incorporated through appropriate modeling of the dimensionless collision frequency, a 惟 dependent error E惟 is introduced and its influence on the overall error is investigated. Here, simulations with the SRT and the MRT model are carried out for power-law fluids in order to numerically investigate the accuracy of non-Newtonian Lattice Boltzmann simulations. A goal of this accuracy analysis is to derive a recommendation for an optimal choice of the time step size and the simulation Mach number, respectively. For the non-Newtonian case, an error estimate for E惟 in the form of a functional is derived on the basis of a series expansion of the Lattice Boltzmann equation. This functional can be solved analytically for the case of the Hagen–Poiseuille channel flow of non-Newtonian fluids. With the help of the error functional, the prediction of the global error minimum of the velocity field is excellent in regions where the E惟 error is the dominant source of error. With an optimal simulation Mach number, the simulation is about one order of magnitude more accurate. Additionally, for both collision models a detailed study of the convergence behavior of the method in the non-Newtonian case is conducted. The results show that the simulation Mach number has a major impact on the convergence rate and second order accuracy is not preserved for every choice of the simulation Mach number.