Detailed analysis of the lattice Boltzmann method on unstructured grids

Marek Krzysztof Misztal,Anier Hernandez-Garcia,Rastin Matin,Henning Osholm Sørensen,Joachim Mathiesen

DOI: https://doi.org/10.48550/arXiv.1409.2754

2014-09-09

Abstract:The lattice Boltzmann method has become a standard for efficiently solving problems in fluid dynamics. While unstructured grids allow for a more efficient geometrical representation of complex boundaries, the lattice Boltzmann methods is often implemented using regular grids. Here we analyze two implementations of the lattice Boltzmann method on unstructured grids, the standard forward Euler method and the operator splitting method. We derive the evolution of the macroscopic variables by means of the Chapman-Enskog expansion, and we prove that it yields the Navier-Stokes equation and is first order accurate in terms of the temporal discretization and second order in terms of the spatial discretization. Relations between the kinetic viscosity and the integration time step are derived for both the Euler method and the operator splitting method. Finally we suggest an improved version of the bounce-back boundary condition. We test our implementations in both standard benchmark geometries and in the pore network of a real sample of a porous rock.

Fluid Dynamics

What problem does this paper attempt to address?

The problem this paper attempts to address is how to effectively simulate fluid dynamics problems in complex geometries when implementing the lattice Boltzmann method (LBM) on unstructured grids. Specifically, the paper analyzes two time discretization schemes (forward Euler method and operator splitting method) for LBM implementation on unstructured grids, derives the evolution process of macroscopic variables, proves that these methods can obtain the Navier-Stokes equations, and analyzes the relationship between numerical viscosity and integration time step. Additionally, the paper proposes an improved bounce-back boundary condition and tests these implementations through standard benchmark geometries and pore networks of actual porous rock samples.

### Main Issues:
1. **LBM Implementation on Unstructured Grids**:
- Unstructured grids allow for more efficient representation of complex geometric boundaries, whereas traditional LBM typically uses regular grids.
- The paper analyzes the implementation of two time discretization schemes (forward Euler method and operator splitting method) on unstructured grids.

2. **Evolution of Macroscopic Variables**:
- Through Chapman-Enskog expansion, the evolution process of macroscopic variables is derived, proving that these methods can obtain the Navier-Stokes equations.
- The relationship between numerical viscosity and integration time step is analyzed, concluding that the time accuracy of the forward Euler method is first-order, and the spatial accuracy is second-order.

3. **Handling of Boundary Conditions**:
- An improved bounce-back boundary condition is proposed to better handle boundary issues in complex geometries.

4. **Experimental Validation**:
- These implementations are tested through standard benchmark geometries and pore networks of actual porous rock samples, validating the correctness of the theoretical derivations and the effectiveness of the methods.

### Key Contributions:
- **Analysis of Time Discretization Schemes**: Detailed analysis of the performance of the forward Euler method and operator splitting method on unstructured grids, deriving the expression for numerical viscosity.
- **Improved Boundary Conditions**: Proposing an improved bounce-back boundary condition that enhances the handling of complex geometries.
- **Experimental Validation**: Validating the effectiveness and accuracy of the methods through multiple experiments, particularly in the application of flow in porous media.

Overall, this paper aims to improve the efficiency and accuracy of simulating fluid dynamics problems in complex geometries by enhancing the LBM implementation on unstructured grids.