期刊
COMPUTERS & FLUIDS
卷 233, 期 -, 页码 -出版社
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.compfluid.2021.105224
关键词
Kinetic equations; Boltzmann model; Discontinuous Galerkin schemes on; unstructured meshes; Implicit-Explicit Runge Kutta; Implicit-Explicit Linear Multistep; High order of accuracy in space and time
资金
- Italian Ministry of Instruction, University and Research (MIUR) [2017KKJP4X]
In this work, deterministic approximate solutions of the Boltzmann model were successfully computed by combining different technologies for discretization and suitable time integration techniques, showing good precision with lower computational cost than traditional methods.
Deterministic solutions of the Boltzmann equation represent a real challenge due to the enormous computational effort which is required to produce such simulations and often stochastic methods such as Direct Simulation Monte Carlo (DSMC) are used instead due to their lower computational cost. In this work, we show that combining different technologies for the discretization of the velocity space and of the physical space coupled with suitable time integration techniques, it is possible to compute very precise deterministic approximate solutions of the Boltzmann model in different regimes, from extremely rarefied to dense fluids, with CFL conditions only driven by the hyperbolic transport term. To that aim, we develop modal Discontinuous Galerkin (DG) Implicit-Explicit Runge Kutta schemes (DG-IMEX-RK) and Implicit-Explicit Linear Multistep Methods based on Backward-Finite-Differences (DG-IMEX-BDF) for solving the Boltzmann model on multidimensional unstructured meshes. The solution of the Boltzmann collision operator is obtained through fast spectral methods, while the transport term in the governing equations is discretized relying on an explicit shock-capturing DG method on polygonal tessellations in the physical space. A novel class of WENO-type limiters, based on a shifting of the moments of inertia for each zone of the mesh, is used to control spurious oscillations of the DG solution across discontinuities. The use of Linear Multistep Methods (LMM) allows the Boltzmann solutions to be consistent not only with the compressible Euler limit but also with the Navier-Stokes asymptotic regime. In addition, as numerically proven, they also permit to strongly reduce the computational effort compared to Runge-Kutta approaches while maintaining the same or even larger accuracy. The performances of these different time discretization techniques are measured comparing both precision and efficiency. At the same time, comparisons against simpler relaxation type kinetic models such as the BGK model are proposed. The order of convergence is numerically measured for different regimes and found to agree with the theoretical findings. The new methods are validated considering two-dimensional benchmark test cases typically used in the fluid dynamics community. A prototype engineering problem consisting of a supersonic flow around a NACA 0012 airfoil with space-time-dependent boundary conditions is also presented for which the pressure coefficients are measured.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据