Fast and scalable turbulent flow simulation with two-way coupling
Fast and scalable turbulent flow simulation with two-way coupling
复制标题
DOI：
10.1145/3386569.3392400
复制
发表时间：
2020-07
期刊：
影响因子：
--
通讯作者：
Wei Li;Yixin Chen;M. Desbrun;Changxi Zheng;Xiaopei Liu
中科院分区：
文献类型：
--
作者：
Wei Li;Yixin Chen;M. Desbrun;Changxi Zheng;Xiaopei Liu
尽管它们在电影中很有吸引力，但涉及流固耦合的湍流在动画中仍然是一个计算挑战。这种电流限制的根源是数值色散，大多数精确的纳维斯托克斯解算器都会受到这种色散的影响：流体和固体之间的适当耦合通常会以局部寄生速度振荡的形式产生人为色散，最终导致数值不稳定。虽然多年来的不断改进导致了保守和细节保留的流体积分器，但这些求解器的分散性很少被讨论，尽管它对流体-结构相互作用产生了巨大的影响。在本文中，我们介绍了一种新型的低耗散低分散流体求解器，它可以高效和可扩展地模拟双向耦合，即使是在湍流情况下也是如此。与大多数目前的CG方法形成鲜明对比的是，我们从来自统计力学的流动的动力学公式中构造了我们的求解器。与现有的格子Boltzmann解算器不同，我们的方法利用高阶矩松弛作为控制结果格式的耗散和色散的关键。此外，我们将新的流体求解器与浸没边界法相结合，通过时间自适应模拟来方便地处理流固耦合。我们的动力学求解器本质上是高度可并行化的，使其非常适合在单GPU或多GPU计算平台上实施。在合成测试和真实实验中与现有解算器进行了广泛的比较，以突出我们的工作在准确性、可伸缩性和效率方面相对于传统和最近的方法的多重优势。
Despite their cinematic appeal, turbulent flows involving fluid-solid coupling remain a computational challenge in animation. At the root of this current limitation is the numerical dispersion from which most accurate Navier-Stokes solvers suffer: proper coupling between fluid and solid often generates artificial dispersion in the form of local, parasitic trains of velocity oscillations, eventually leading to numerical instability. While successive improvements over the years have led to conservative and detail-preserving fluid integrators, the dispersive nature of these solvers is rarely discussed despite its dramatic impact on fluid-structure interaction. In this paper, we introduce a novel low-dissipation and low-dispersion fluid solver that can simulate two-way coupling in an efficient and scalable manner, even for turbulent flows. In sharp contrast with most current CG approaches, we construct our solver from a kinetic formulation of the flow derived from statistical mechanics. Unlike existing lattice Boltzmann solvers, our approach leverages high-order moment relaxations as a key to controlling both dissipation and dispersion of the resulting scheme. Moreover, we combine our new fluid solver with the immersed boundary method to easily handle fluid-solid coupling through time adaptive simulations. Our kinetic solver is highly parallelizable by nature, making it ideally suited for implementation on single- or multi-GPU computing platforms. Extensive comparisons with existing solvers on synthetic tests and real-life experiments are used to highlight the multiple advantages of our work over traditional and more recent approaches, in terms of accuracy, scalability, and efficiency.