An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Free convection of nanofluid filled enclosure using lattice Boltzmann method (LBM)

The lattice Boltzmann method （LBM） is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles： TiO2, Ag, Cu, and A1203. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical meth- ods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ =2.5 when Ra = 106 while at A = 1.5 for other Rayleigh numbers.

Numerical simulation of dynamic process for liquid film spreading by lattice Boltzmann method and its experimental verification

Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.

NEW STUDYING OF LATTICE BOLTZMANN METHOD FOR TWO-PHASE DRIVEN IN POROUS MEDIA

By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore, the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.

Study on Free Surface and Channel Flow Induced by Low Temperature Plasma via Lattice Boltzmann Method

Active boundary layer flow control and boundary layer manipulation in the channel flow that was based on low temperature plasma were studied by means of a lattice Boltzmann method. Two plasma actuators were placed in a row to obtain the influence rule of their separation distance on the velocity profile at three locations and maximum velocity in the flow field. Two plasma actuators were placed symmetrically inside a channel to examine the effect of channel height and voltage on the velocity profile and flow rate. It was found that the channel height controls the distribution of flow velocity, which affected the flow rate and its direction. Increasing plasma voltage had a negative effect on the flow rate due to the generation of a larger and stronger flow vortex.

Rotational Flow in a Narrow Annular Gap Based on Lattice Boltzmann Method

For comprehensive characteristics of flow in a gas bearing,lattice Boltzmann method(LBM)is applied for study of the two-dimensional flow between two eccentric cylinders with the inner one rotating at a high speed.The flow pattern and circumferential pressure distribution are discussed based on critical issues such as eccentricity ranging from 0.2 to 0.9,clearance ratio varying from 0.005 to 0.01 and rotating speed in the range of 3×104—1.8×105 r/min.The analysis and discussion on the circumferential pressure distribution affirmed the quasilinear relation between the extremum pressure and rotating speed.Furthermore,a high eccentricity and small clearance ratio contributes most to the fluctuation of the circumferential pressure distribution.The flow pattern inside the channel exhibits separation vortex under a large eccentricity.The conclusions drawn in this work give rise to prediction of the flow pattern in the gas bearing which is beneficial for evaluating the performance of as well as instructing the design and development.

Data-driven optimization study of the multi-relaxation-time lattice Boltzmann method for solid-liquid phase change

Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical performance during simulation;however,the value method of the relaxation parameters needs to be specified.Therefore,in this study,a random forest(RF)model is used to discriminate the importance of different relaxation parameters to the convergence,and a support vector machine(SVM)is used to explore the decision boundary of the convergent samples in each dimensional model.The results show that the convergence of the samples is consistent with the sign of the decision number,and two types of the numerical deviations appear,i.e.,the phase mushy zone and the non-physical heat transfer.The relaxation parameters chosen on the decision boundary can further suppress the numerical bias and improve numerical accuracy.

Airfoil design optimization based on lattice Boltzmann method and adjoint approach

We present a new aerodynamic design method based on the lattice Boltzmann method （LBM） and the adjoint approach. The flow field and the adjoint equation are numerically simulated by the GILBM （generalized form of interpolation supplemented LBM） on non-uniform meshes. The first-order approximation for the equilibrium dis- tribution function on the boundary is proposed to diminish the singularity of boundary conditions. Further, a new treatment of the solid boundary in the LBM is described par- ticularly for the airfoil optimization design problem. For a given objective function, the adjoint equation and its boundary conditions are derived analytically. The feasibility and accuracy of the new approach have been perfectly validated by the design optimization of NACA0012 airfoil.

High velocity flow simulation using lattice Boltzmann method with no-free-parameter dissipation scheme

A finite difference lattice Boltzmann method of second-order accuracy in time is developed based on non-oscillatory scheme with no-free-parameter dissipation （NND） difference scheme in this paper. The NND lattice Boltzmann method is used to simulate high-speed flows by constructing a new equilibrium distribution function of the lattice Boltzmann method. Compared with a variation of lattice Boltzmann method developed by Qu, et al., the present method can capture shock waves and handle oscillations of high velocity flows accurately in larger time steps and in shorter computing time. Numerical results indicate the correctness and capability of simulating shock wave interactions of the NND lattice Boltzmann method.

双重网格的Lattice Boltzmann方法

针对均匀网格的LatticeBoltzmann方法(LBM)计算效率较低的不足,引入有限差分方法(FDM)中双重网格算法的基本思想,结合LBM在连续流场上变量也连续分布的特点,推导出双重网格的LBM.这种算法充分利用粗、细两种网格各自不同的优点,加快了收敛速度,提高了流场的计算效率.通过模拟二维Poiseulle流动和后台阶流动两个算例,证明双重网格算法的计算时间明显地少于均匀网格算法的计算时间,计算结果正确可信.

Lattice Boltzmann方法不同边界条件的混用

在Lattice Boltzmann方法数值模拟中,比较重要的一个环节就是边界条件的确定.通过对不同边界条件混用情况的研究,发现不同的边界条件确定方法可以相互混用,说明了Lattice Bolt-zmann方法对于边界的处理是比较灵活的.

结合亚格子模型的Lattice-Boltzmann算法

介绍一种亚格子模型与LB(Lattic-Boltzmann)方法相结合的算法,并用其对亚临界雷诺数情况下的圆柱绕流问题进行数值模拟,得到流场的流线图以及圆柱表面的压力系数。通过与实验值相比较,发现数值结果与实验值是比较吻合的,证明该算法可以用于计算雷诺数较大的流体问题,同时具有LB方法计算效率高的优点。

用Lattice Boltzmann方法模拟1／4圆腔内的定常层流运动

分析格子Boltzmann方法中二阶精度的曲线边界处理方法．应用格子Boltzmann方法及其边界处理方法模拟1／4圆腔内的定常层流运动，引入流线图和等涡线图分析流场随如数的变化．并且发现当Re数在10～100区间内变化时，随着Re数的增大，顺时针旋转流场的涡心位置偏离X轴的角度逐渐减小，而逆时针旋转流场的涡心位置偏离X轴的角度却越来越大．

一种高效建模数值算法及其应用——Lattice Boltzmann方法

近年来,格子Boltzmann方法(LBM)已发展为一种模拟流体和物理问题的新颖的、有前景的数值方法,在许多领域的各种数值问题求解上取得很大的成功。文章介绍了一种模拟复杂流动的高效建模数值算法Lattice Boltzmann方法和它的基本原理及其应用。并通过两个实例数值模拟计算,说明Lattice Boltzmann方法正确、有效,并展示了广阔的应用前景,为今后更深入的研究和广泛应用打下基础。

基于Lattice Boltzmann方法的页岩张性裂缝渗流特征研究

为深入了解张性裂缝结构中流体渗流特征,采用物理实验和数值模拟相结合的方法,系统研究流体在页岩张性裂缝中的渗流特征。将Barnett页岩岩心劈裂成张性裂缝,并采集裂缝面数据,将其作为物理模型,同时采用LBM-BKG方法 D3Q19模型研究流体在张性裂缝中的流动。研究尺度依次为整体、截面、截线、离散单元体,考察各尺度下渗流场特征参数分布,并分析原因。研究表明:整体上,沿渗流横向和纵向上压力均呈线性分布,横向各截线上压力呈波动式线性分布,纵向截线上压力分布呈非线性;缝面端部局部区域会产生回流区,回流区与粗糙性及粗糙性变化程度有关;近平均速度的速度等值面间基本呈平行分布,速度较大时,等值面基本对应于端面相对平缓的区域;沿同一流线,渗流速度略有差别;整体渗流为稳态,不代表渗流场各处速度相等,离散的渗流单元体的渗流速度各不相同;三维粗糙性造成了渗流速度的三维分布,但主要以渗流方向和垂向上渗流速度为主。该项研究对深入理解页岩等裂缝发育储层的裂缝系统渗流具有重要意义。

基于神威加速计算架构的LBM多级并行计算

格子玻尔兹曼方法(lattice Boltzmann method,LBM)是一种基于分子运动理论计算流体力学(computational fluid dynamics,CFD)的方法,提高LBM的并行计算能力是高性能计算领域的一项重要的研究内容.本文基于SW26010Pro处理器,通过区域分解、数据重构、双缓冲、向量化等优化方法,实现了LBM的多级并行.基于以上优化方案,测试了5600万网格规模,实现结果显示,相比于MPI进行级并行,碰撞过程的平均加速倍数达到61.737、迁移过程的平均加速倍数达到17.3,同时对方腔流案例做了强扩展测试,网格规模为1200×1200×1200,以6.2万计算核心为基准,百万核心的并行效率超过60.5%.

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.

Boundary scheme for lattice Boltzmann modeling of micro-scale gas flow in organic-rich pores considering surface diffusion

We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.

Lattice Boltzmann simulation of MHD natural convection in a cavity with porous media and sinusoidal temperature distribution

The lattice Boltzmann method(LBM) is used to simulate the effect of magnetic field on the natural convection in a porous cavity. The sidewalls of the cavity are heated sinusoidally with a phase derivation, whereas the top and bottom walls are thermally insulated. Numerical simulation is performed, and the effects of the pertinent parameters, e.g., the Hartmann number, the porosity, the Darcy number, and the phase deviation, on the fluid flow and heat transfer are investigated. The results show that the heat transfer is affected by the temperature distribution on the sidewalls clearly. When the Hartmann number is 0, the maximum average Nusselt number is obtained at the phase deviation 90?. Moreover, the heat transfer enhances when the Darcy number and porosity increase, while decreases when the Hartman number increases.

Comparative study of two lattice Boltzmann multiphase models for simulating wetting phenomena: implementing static contact angles based on the geometric formulation

Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such systems is important to both scientific research and practical applications. It is challenging due to the complexity of the phenomena and difficulties in choosing an appropriate numerical method. To provide some detailed guidelines for selecting a suitable multiphase lattice Boltzmann model, two kinds of lattice Boltzmann multiphase models, the modified S-C model and the H-C-Z model, are used in this paper to investigate the static contact angle on solid surfaces with different wettability combined with the geometric formulation(Ding, H. and Spelt, P.D. M. Wetting condition in diffuse interface simulations of contact line motion. Physical Review E, 75(4), 046708(2007)). The specific characteristics and computational performance of these two lattice Boltzmann method(LBM) multiphase models are analyzed including relationship between surface tension and the control parameters, the achievable range of the static contact angle, the maximum magnitude of the spurious currents(MMSC), and most importantly, the convergence rate of the two models on simulating the static contact angle. The results show that a wide range of static contact angles from wetting to non-wetting can be realized for both models. MMSC mainly depends on the surface tension. With the numerical parameters used in this work, the maximum magnitudes of the spurious currents of the two models are on the same order of magnitude. MMSC of the S-C model is universally larger than that of the H-C-Z model. The convergence rate of the S-C model is much faster than that of the H-C-Z model. The major foci in this work are the frequently-omitted important details in simulating wetting phenomena. Thus, the major findings in this work can provide suggestions for simulating wetting phenomena with LBM multiphase models along with the geometric formulation.