Search

The use of flux limiters is widespread within the scientific computing community to capture shock discontinuities and are of paramount importance for the temporal integration of high-speed aerodynamics, multiphase flows, and hyperbolic equations in general. Meanwhile, the breakthrough of new computing architectures and the hybridization of supercomputer systems pose a huge portability challenge, particularly for legacy codes, since the computing subroutines that form the algorithms, the so-called kernels, must be adapted to various complex parallel programming paradigms. From this perspective, the development of innovative implementations relying on a minimalist set of kernels simplifies the deployment of scientific computing software on state-of-the-art supercomputers, while it requires the reformulation of algorithms, such as the aforementioned flux limiters. Equipped with basic algebraic topology and graph theory underlying the classical mesh concept, a new flux limiter formulation is presented based on the adoption of algebraic data structures and kernels. As a result, traditional flux limiters are cast into a stream of only two types of computing kernels: sparse matrix-vector multiplication and generalized pointwise binary operators. The newly proposed formulation eases the deployment of such a numerical technique in massively parallel, potentially hybrid, computing systems and is demonstrated for a canonical advection problem.

In our recent work [H. Zhang, F.X. Trias, A. Oliva, D. Yang, Y. Tan, Y. Sheng. PIBM: Particulate immersed boundary method for fluid-particle interaction problems. Powder Technology. 272(2015), 1-13.], a particulate immersed boundary method (PIBM) for simulating fluid-particle multiphase flow was proposed and assessed in both two- and three-dimensional applications. In this study, the PIBM was extended to solve thermal interaction problems between spherical particles and fluid. The Lattice Boltzmann Method (LBM) was adopted to solve the fluid flow and temperature fields, the PIBM was responsible for the non-slip velocity and temperature boundary conditions at the particle surface, and the kinematics and trajectory of the solid particles were evaluated by the Discrete Element Method (DEM). Four case studies were implemented to demonstrate the capability of the current coupling scheme. Firstly, numerical simulation of natural convection in a two-dimensional square cavity with an isothermal concentric annulus was carried out for verification purpose. The current results were found to have good agreements with previous references. Then, sedimentation of two- and three-dimensional isothermal particles in fluid was numerically studied, respectively. The instantaneous temperature distribution in the cavity was captured. The effect of the thermal buoyancy on particle behaviors was discussed. At last, sedimentation of three-dimensional thermosensitive particles in fluid was numerically investigated. Our results revealed that the LBM-PIBM-DEM is a promising scheme for the solution of complex fluid-particle interaction problems with heat transfer.

It is well known that the number of particles should be scaled up to enable industrial scale simulation. The calculations are more computationally intensive when the motion of the surrounding fluid is considered. Besides the advances in computer hardware and numerical algorithms, the coupling scheme also plays an important role on the computational efficiency. In this study, a particle immersed boundary method (PIBM) for simulating the fluid-particle multiphase flow was presented and assessed in both two- and three-dimensional applications. The idea behind PIBM derives from the conventional momentum exchange-based immersed boundary method (IBM) by treating each Lagrangian point as a solid particle. This treatment enables LBM to be coupled with fine particles residing within a particular grid cell. Compared with the conventional IBM, dozens of times speedup in two-dimensional simulation and hundreds of times in three-dimensional simulation can be expected under the same particle and mesh number. Numerical simulations of particle sedimentation in the Newtonian flows were conducted based on a combined lattice Boltzmann method - particle immersed boundary method - discrete element method scheme, showing that the PIBM can capture the feature of particulate flows in fluid and is indeed a promising scheme for the solution of the fluid-particle interaction problems., Comment: 30 Pages, 15 figures