一种新的模拟渗流运动的数值方法

根据格子Boltzmann方法及相关理论，建立了一个新的模拟渗流运动的数值模型，所得模型没有在边界上采取相应平均措施，同时还避免了一些非物理副产品的出现-实例计算数值结果与精确解符合较好，证明模型可靠- 关键词： 渗流 格子Boltzmann方法 数值模型     

Experiment and Lattice Boltzmann numerical study on nanofluids flow in a micromodel as porous medium

Al₂O₃ nanofluids flow has been studied in etched glass micromodel which is idealization of porous media by using a pseudo 2D Lattice Boltzmann Method (LBM). The predictions were compared with experimental results. Pressure drop / flow rate relations have been measured for pure water and Al₂O₃ nanofluids. Because the size of Al₂O₃ nanoparticles is tiny enough to permit through the pore throats of the micromodel, blockage does not occur and the permeability is independent of the nanofluid volume fraction. Therefore, the nanofluid behaves as a single phase fluid, and a single phase LBM is able to simulate the results of this experiment. Although the flow in micromodels is 3D, we showed that 2D LBM can be used provided an effective viscous drag force, representing the effect of the third dimension, is considered. Good qualitative and quantitative agreement is seen between the numerical and experimental results.     

多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀)过程, 当雷诺数越大时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩散控制, 此时反应主要发生在上游及大孔隙区域.     

A Lattice Boltzmann Kinetic Model for Microflow and Heat Transfer

In this paper, we propose a lattice Boltzmann BGK model for simulation of micro flows with heat transfer based on kinetic theory and the thermal lattice Boltzmann method (He et al., J. Comp. Phys. 146:282, 1998). The relaxation times are redefined in terms of the Knudsen number and a diffuse scattering boundary condition (DSBC) is adopted to consider the velocity slip and temperature jump at wall boundaries. To check validity and potential of the present model in modelling the micro flows, two two-dimensional micro flows including thermal Couette flow and thermal developing channel flow are simulated and numerical results obtained compare well with previous studies of the direct simulation Monte Carlo (DSMC), molecular dynamics (MD) approaches and the Maxwell theoretical analysis     

Heuristic optimality criterion algorithm for shape design of fluid flow

This paper presents a heuristic optimality criterion algorithm for shape design of fluid flow. In this algorithm, the lattice Boltzmann method (LBM) is utilized to calculate the flow field of a fluid domain which is divided into elemental cells. A heuristic optimality criterion is applied for cells at the solid–fluid interface, i.e. the dynamic pressure for fluid cells and the viscous stress on their neighboring solid cells. An automatic program is processed step by step to exchange the positions of solid and fluid cells identified by the optimality criterion, with the objective of decreasing the flow resistance at the constraint of constant fluid volume. To illustrate the procedure of this algorithm for shape design of fluid flow, two simple examples are presented: one with fluid flowing through a right angle elbow and the other through a converging T-junction. Numerical results show that this algorithm can successfully reduce the total pressure drop of the system, demonstrating its potential applications in engineering optimal design.     

用格子玻尔兹曼方法研究流动-反应耦合的非线性渗流问题

根据格子玻尔兹曼计算技术以及相应渗流理论，对多孔介质内流动-反应(矿物介质的溶解等)耦合这一非线性渗流问题进行了数值研究，计算结果与解析解基本符合.数字图像重构技术反映的结果表明流体流动和反应之间可以发生强烈的耦合和反耦合作用，同时可以形成条带结构这一自组织现象，与实验和其他理论分析结果符合也很好. 关键词： 非线性渗流 耦合反应 数值模型     

Lattice Boltzmann method for the generalized Kuramoto-Sivashinsky equation

In this paper, a lattice Boltzmann model with an amending function is proposed for the generalized Kuramoto-Sivashinsky equation that has the form u_t+uu_x+αu_xx+βu_xxx+γu_xxxx=0. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. It is found that the numerical results agree well with the analytical solutions.     

Three-dimensional lattice Boltzmann method for simulating blood flow in aortic arch

The three-dimensional （3D） lattice Boltzmann models, 3DQ15, 3DQ19 and 3DQ27, under different wall boundary conditions and lattice resolutions have been investigated by simulating Poiseuille flow in a circular cylinder for a wide range of Reynolds numbers. The 3DQ19 model with improved Fillippova and Hanel （FH） curved boundary condition represents a good compromise between computational efficiency and reliability. Blood flow in an aortic arch is then simulated as a typical haemodynamic application. Axial and secondary fluid velocity and effective wall shear stress profiles in a 180° bend are obtained, and the results also demonstrate that the lattice Boltzmann method is suitable for simulating the flow in 3D large-curved vessels.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

Lattice Boltzmann Model for Free Surface Flow for Modeling Foaming

We present a 2D- and 3D-lattice Boltzmann model for the treatment of free surface flows including gas diffusion. Interface advection and related boundary conditions are based on the idea of the lattice Boltzmann equation. The fluid dynamic boundary conditions are approximated by using the mass and momentum fluxes across the interface, which do not require explicit calculation of gradients. A similar procedure is applied to fulfill the diffusion boundary condition. Simple verification tests demonstrate the correctness of the algorithms. 2D- and 3D-foam evolution examples demonstrate the potential of the method.     

An improved immersed boundary-lattice Boltzmann method for simulating three-dimensional incompressible flows

The recently proposed boundary condition-enforced immersed boundary-lattice Boltzmann method (IB-LBM) [14] is improved in this work to simulate three-dimensional incompressible viscous flows. In the conventional IB-LBM, the restoring force is pre-calculated, and the non-slip boundary condition is not enforced as compared to body-fitted solvers. As a result, there is a flow penetration to the solid boundary. This drawback was removed by the new version of IB-LBM [14], in which the restoring force is considered as unknown and is determined in such a way that the non-slip boundary condition is enforced. Since Eulerian points are also defined inside the solid boundary, the computational domain is usually regular and the Cartesian mesh is used. On the other hand, to well capture the boundary layer and in the meantime, to save the computational effort, we often use non-uniform mesh in IB-LBM applications. In our previous two-dimensional simulations [14], the Taylor series expansion and least squares-based lattice Boltzmann method (TLLBM) was used on the non-uniform Cartesian mesh to get the flow field. The final expression of TLLBM is an algebraic formulation with some weighting coefficients. These coefficients could be computed in advance and stored for the following computations. However, this way may become impractical for 3D cases as the memory requirement often exceeds the machine capacity. The other way is to calculate the coefficients at every time step. As a result, extra time is consumed significantly. To overcome this drawback, in this study, we propose a more efficient approach to solve lattice Boltzmann equation on the non-uniform Cartesian mesh. As compared to TLLBM, the proposed approach needs much less computational time and virtual storage. Its good accuracy and efficiency are well demonstrated by its application to simulate the 3D lid-driven cubic cavity flow. To valid the combination of proposed approach with the new version of IBM [14] for 3D flows with curved boundaries, the flows over a sphere and torus are simulated. The obtained numerical results compare very well with available data in the literature.     

超声波对多孔介质中两相流动的影响

近年来，超声技术已被应用于采油工程中，在油井解堵，水井增注等方面发挥了重要的作用。     

Direct pore-level modeling of incompressible fluid flow in porous media

We present a dynamic particle-based model for direct pore-level modeling of incompressible viscous fluid flow in disordered porous media. The model is capable of simulating flow directly in three-dimensional high-resolution micro-CT images of rock samples. It is based on moving particle semi-implicit (MPS) method. We modify this technique in order to improve its stability for flow in porous media problems. Using the micro-CT image of a rock sample, the entire medium, i.e., solid and fluid, is discretized into particles. The incompressible Navier–Stokes equations are then solved for each particle using the MPS summations. The model handles highly irregular fluid–solid boundaries effectively. An algorithm to split and merge fluid particles is also introduced. To handle the computational load, we present a parallel version of the model that runs on distributed memory computer clusters. The accuracy of the model is validated against the analytical, numerical, and experimental data available in the literature. The validated model is then used to simulate both unsteady- and steady-state flow of an incompressible fluid directly in a representative elementary volume (REV) size micro-CT image of a naturally-occurring sandstone with 3.398 μm resolution. We analyze the quality and consistency of the predicted flow behavior and calculate absolute permeability using the steady-state flow rate.     

Numerical analysis of fluid flow through a cylinder array using a lattice Boltzmann model

In this paper we present a detailed computational study of an incompressible Newtonian fluid flow across a periodic array of two-dimensional cylinders which is a simplest non-trivial representation of a porous media. A two-dimensional Lattice Boltzmann Method is used to solve the governing Navier－Stokes equation taking into account of viscous dissipation effects and influence of nonlinear fluid drag. Both the flow fields and the Darcy－Forchheimer drag coefficient as a function of the solid volume fraction are calculated for a wide range of flow Reynolds numbers. The predictions were compared with the results from conventional numerical and empirical models for verification. Apart from confirming that inertial effects can cause a significant deviation from Darcy's law for large velocities the results also show that the characteristics of the vorticity field vary considerably as the Reynolds number increases, which will have major implications to the transport of passive particulate substances within the pores and their removal rate.     

Computational modeling of MR flow imaging by the lattice Boltzmann method and Bloch equation

In this work, a computational model of magnetic resonance (MR) flow imaging is proposed. The first model component provides fluid dynamics maps by applying the lattice Boltzmann method. The second one uses the flow maps and couples MR imaging (MRI) modeling with a new magnetization transport algorithm based on the Eulerian coordinate approach. MRI modeling is based on the discrete time solution of the Bloch equation by analytical local magnetization transformations (exponential scaling and rotations).     

A cumulant analysis for non-Gaussian displacement distributions in Newtonian and non-Newtonian flows through porous media

We use displacement encoding pulsed field gradient (PFG) nuclear magnetic resonance to measure Fourier components S(q) of flow displacement distributions P(zeta) with mean displacement (zeta) for Newtonian and non-Newtonian flows through rocks and bead packs. Displacement distributions are non-Gaussian; hence, there are finite terms above second order in the cumulant expansion of ln(S(q)). We describe an algorithm for an optimal self-consistent cumulant analysis of data, which can be used to obtain the first three (central) moments of a non-Gaussian P(zeta), with error bars. The analysis is applied to Newtonian and non-Newtonian flows in rocks and beads. Flow with shear-thinning xanthan solution produces a 15.6+/-2.3% enhancement of the variance sigma(2) of displacement distributions when compared to flow experiments with water.     

Study of force-dependent and time-dependent transition of secondary flow in a rotating straight channel by the lattice Boltzmann method

A numerical study using the lattice Boltzmann method has been carried out for flow through a rotating straight channel with a rectangular cross section. With different forces applied, the secondary flow exhibits two-cell states, four-cell states or six-cell states at a range of low rotational Reynolds number, however, within which only the two-cell states have been commonly reported. In addition to the force-dependent flow transition, a time-dependent flow transition of the secondary flow among two-cell states, four-cell states and six-cell states is also discovered during flow development. These newly found flow transitions and their regulations by force application have been analyzed. Based on numbers of case studies, it is found that a dimensionless number, the ratio of the driving pressure gradient to the centrifugal force, regulates such flow transitions. This study not only releases new phenomena of flow transition, but also indicates new applications in flow control, particle separation and heat transfer.     

Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media

The multiscale finite-volume (MSFV) method was originally developed for the solution of heterogeneous elliptic problems with reduced computational cost. Recently, some extensions of this method for parabolic problems have been proposed. These extensions proved effective for many cases, however, they are neither general nor completely satisfactory. For instance, they are not suitable for correctly capturing the transient behavior described by the parabolic pressure equation. In this paper, we present a general multiscale finite-volume method for parabolic problems arising, for example, from compressible multiphase flow in porous media. Opposed to previous methods, here, the basis and correction functions are solutions of full parabolic governing equations in localized domains. At the same time, to enhance the computational efficiency of the scheme, the basis functions are kept pressure independent and do not have to be recalculated as pressure evolves. This general approach requires no additional assumptions and its good efficiency and high accuracy is demonstrated for various challenging test cases. Finally, to improve the quality of the results and also to extend the scheme for highly anisotropic heterogeneous problems, it is combined with the iterative MSFV (i-MSFV) method for parabolic problems. As one iterates, the i-MSFV solutions of compressible multiphase problems (parabolic problems) converge to the corresponding fine-scale reference solutions in the same way as demonstrated recently for incompressible cases (elliptic problems). Therefore, the proposed MSFV method can also be regarded as an efficient linear solver for parabolic problems and studies of its efficiency are presented for many test cases.     

耦合不可压流场输运方程的格子Boltzmann方法研究

基于格子Boltzmann方法,提出了求解耦合不可压缩流场输运方程的一种改进数值方法. 该方法使用格子Boltzmann方法求解流场方程,并根据流场格子模型的密度分布函数构建了输运方程的二阶离散格式. 通过二维平板通道流场输运系统验证了该方法的有效性.数值结果表明,该方法可以有效地减少计算过程中出现的非物理耗散, 并克服了传统模型所需巨大存储量的缺点.     

An interaction potential based lattice Boltzmann method with adaptive mesh refinement (AMR) for two-phase flow simulation

The lattice Boltzmann method (LBM) for two-phase flow simulation is often hindered by insufficient resolution at the interface. As a result, the LBM simulation of bubbles in bubbling flows is commonly limited to spherical or slightly deformed bubble shapes. In this study, the adaptive mesh refinement method for the LBM is developed to overcome such a problem. The approach for this new method is based on the improved interaction potential model, which is able to maintain grid-independent fluid properties in the two-fluid phases and at the interface. The LBM–AMR algorithm is described, especially concerning the LBM operation on a non-uniform mesh and the improved interaction potential model. Numerical simulations have been performed to validate the method in both single phase and multiphase flows. The 2D and 3D simulations of the buoyant rise of bubbles are conducted under various conditions. The agreement between the simulated bubble shape and velocity with experiments illustrates the capability of the LBM–AMR approach in predicting bubble dynamics even under the large bubble deformation conditions. Further, the LBM–AMR technique is capable of simulating a complex topology change of the interface. Integration of LBM with AMR can significantly improve the accuracy and reduce computation cost. The method developed in this study may appreciably enhance the capability of LBM in the simulation of complex multiphase flows under realistic conditions.