-
On h-refined meshless solution to Navier-Stokes problem in porous media: comparing meshless Lattice Boltzman Method with ACM RBF-FD approach
Authors:
Dawid Strzelczyk,
Miha Rot,
Gregor Kosec,
Maciej Matyka
Abstract:
In this paper, two mesh-free CFD solvers for pore-scale fluid flow through porous media are considered, namely the Lattice Boltzmann Method with the two relaxation time collision term and the direct Navier-Stokes solver under the artificial compressibility limit. The porous media is built with a regular arrangement of spherical grains with variable radii, which allows control of the porosity. Both…
▽ More
In this paper, two mesh-free CFD solvers for pore-scale fluid flow through porous media are considered, namely the Lattice Boltzmann Method with the two relaxation time collision term and the direct Navier-Stokes solver under the artificial compressibility limit. The porous media is built with a regular arrangement of spherical grains with variable radii, which allows control of the porosity. Both solvers use the same h-refined meshless spatial discretization to adequately capture the underlying geometry and the same Radial Basis Function (RBF) method to approximate the involved fields and partial differential operators. First, the results are compared with the data from the literature in terms of drag coefficient and permeability at different porosities achieving excellent agreement with the reported results. Next, the simulations are extended beyond the porosity range reported in the literature using proposed h-refined CFD solvers. The results are supported by convergence and timing analyses and discussions on meshless parameters such as stencil size and refinement settings.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Exploiting decoupled discretization in meshless lattice Boltzmann method for high-Reynolds number flows
Authors:
Dawid Strzelczyk,
Maciej Matyka
Abstract:
One of the most severe limitations of the Lattice Boltzmann Method in the context of simulating inertial flows is the coupling of the discretization of space to the velocity discretization. It causes the need to increase the size of computational meshes whenever one wants to increase the Reynolds number in the system for a fixed velocity and viscosity. This work aims at adopting the recently propo…
▽ More
One of the most severe limitations of the Lattice Boltzmann Method in the context of simulating inertial flows is the coupling of the discretization of space to the velocity discretization. It causes the need to increase the size of computational meshes whenever one wants to increase the Reynolds number in the system for a fixed velocity and viscosity. This work aims at adopting the recently proposed meshless LBM formulation to the problem of high Reynolds number flows by using its internal property of decoupled space and velocity discretizations. In meshless formulation, one can change the reference length of the problem by scaling the distances between points that discretize the domain while the streaming distance remains fixed. By doing that, we increase the Reynolds number at no cost: additional nodes need not be inserted in the discretization. We measure the accuracy and efficiency of this approach in the flow around a circular obstacle at $Re=84$--$168$, the lid-driven cavity flow at $Re=100$--$5000$, and the flow through a porous sample at $Re=0.8$--$1500$. Additionally, we apply the meshless streaming step to the recently proposed fixed relaxation time $τ=1$ LBM to extend its applicability to model inertial flows.
△ Less
Submitted 5 July, 2024; v1 submitted 22 December, 2023;
originally announced December 2023.
-
Study of the convergence of the Meshless Lattice Boltzmann Method in Taylor-Green and annular channel flows
Authors:
Dawid Strzelczyk,
Maciej Matyka
Abstract:
The Meshless Lattice Boltzmann Method (MLBM) is a numerical tool that relieves the standard Lattice Boltzmann Method (LBM) from regular lattices and, at the same time, decouples space and velocity discretizations. In this study, we investigate the numerical convergence of MLBM in two benchmark tests: the Taylor-Green vortex and annular (bent) channel flow. We compare our MLBM results to LBM and to…
▽ More
The Meshless Lattice Boltzmann Method (MLBM) is a numerical tool that relieves the standard Lattice Boltzmann Method (LBM) from regular lattices and, at the same time, decouples space and velocity discretizations. In this study, we investigate the numerical convergence of MLBM in two benchmark tests: the Taylor-Green vortex and annular (bent) channel flow. We compare our MLBM results to LBM and to the analytical solution of the Navier-Stokes equation. We investigate the method's convergence in terms of the discretization parameter, the interpolation order, and the LBM streaming distance refinement. We observe that MLBM outperforms LBM in terms of the error value for the same number of nodes discretizing the domain. We find that LBM errors at a given streaming distance $δx$ and timestep length $δt$ are the asymptotic lower bounds of MLBM errors with the same streaming distance and timestep length. Finally, we suggest an expression for the MLBM error that consists of the LBM error and other terms related to the semi-Lagrangian nature of the discussed method itself.
△ Less
Submitted 31 October, 2023; v1 submitted 2 June, 2023;
originally announced June 2023.
-
Deep learning for diffusion in porous media
Authors:
Krzysztof M. Graczyk,
Dawid Strzelczyk,
Maciej Matyka
Abstract:
We adopt convolutional neural networks (CNN) to predict the basic properties of the porous media. Two different media types are considered: one mimics the sand packings, and the other mimics the systems derived from the extracellular space of biological tissues. The Lattice Boltzmann Method is used to obtain the labeled data necessary for performing supervised learning. We distinguish two tasks. I…
▽ More
We adopt convolutional neural networks (CNN) to predict the basic properties of the porous media. Two different media types are considered: one mimics the sand packings, and the other mimics the systems derived from the extracellular space of biological tissues. The Lattice Boltzmann Method is used to obtain the labeled data necessary for performing supervised learning. We distinguish two tasks. In the first, networks based on the analysis of the system's geometry predict porosity and effective diffusion coefficient. In the second, networks reconstruct the concentration map. In the first task, we propose two types of CNN models: the C-Net and the encoder part of the U-Net. Both networks are modified by adding a self-normalization module [Graczyk \textit{et al.}, Sci Rep 12, 10583 (2022)]. The models predict with reasonable accuracy but only within the data type, they are trained on. For instance, the model trained on sand packings-like samples overshoots or undershoots for biological-like samples. In the second task, we propose the usage of the U-Net architecture. It accurately reconstructs the concentration fields. In contrast to the first task, the network trained on one data type works well for the other. For instance, the model trained on sand packings-like samples works perfectly on biological-like samples. Eventually, for both types of the data, we fit exponents in the Archie's law to find tortuosity that is used to describe the dependence of the effective diffusion on porosity.
△ Less
Submitted 6 June, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.