A fictitious domain method with distributed Lagrange multipliers on adaptive quad/octrees for the direct numerical simulation of particle-laden flows

JOURNAL OF COMPUTATIONAL PHYSICS

卷:430

Article

Selcuk, Can²  Ghigo, Arthur R.^2,3,4  Popinet, Stephane¹  Wachs, Anthony^2,3 

[1] Sorbonne Univ, Inst Jean Rond Alembert, CNRS, F-75005 Paris, France

[2] Univ British Columbia, Dept Math, 1984 Math Rd, Vancouver, BC V6T 1Z2, Canada

[3] Univ British Columbia, Dept Chem & Biol Engn, 2360 East Mall, Vancouver, BC V6T 1Z3, Canada

[4] Univ British Columbia, PIMS CNRS, 4176-2207 Main Mall, Vancouver, BC V6T 1Z4, Canada

关键词: Fictitious domain;    Adaptive grid;    Quadtree/octree;    Particle-laden flow;    Parallel computing;   

DOI  :  10.1016/j.jcp.2020.109954

来源: Elsevier

【 摘 要 】

In this article, we extend our Distributed Lagrange Multiplier/Fictitious Domain method previously implemented on simple regular Cartesian grids to quadtree/octree adaptive grids. The objective is to improve both the accuracy and efficiency of our DLM/FD particle resolved simulation method by extending its computing capabilities through dynamic local mesh refinement. The main features of our numerical method, such as a first-order operator splitting time algorithm and a second-order reconstruction of the velocity field close to the boundary of the immersed rigid bodies (of arbitrary shape), are unchanged. We implemented our adaptive DLM/FD algorithm within Basilisk, a parallel platform to solve partial differential equations on dynamic quadtree/octree grids. The quadtree/octree structure of the grid and specific design rules of Basilisk impose a special treatment of some of the operations performed on the grid in the DLM/FD-Uzawa algorithm. The new computational method is then tested and validated on a set of flow configurations including the challenging problem of accurately computing lubrication interaction forces without resorting to using any ad hoc correction. Finally, we illustrate the potential of our code to compute complex particle-laden flow configurations that were not attainable in the past with a DLM/FD algorithm implemented on a simple regular Cartesian grid. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free