幂律流体圆柱绕流的格子波尔兹曼模拟简

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明：插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.     

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.      

Comparative study of forced convection of a power-law fluid in a channel with a built-in square cylinder

A detailed comparison between the lattice Boltzmann method and the finite element method is presented for an incompressible steady laminar flow and heat transfer of a power-law fluid past a square cylinder between two parallel plates. Computations are performed for three different blockage ratios (ratios of the square side length to the channel width) and different values of the power-law index n covering both pseudo-plastic fluids (n < 1) and dilatant fluids (n > 1). The methodology is validated against the exact solution. The local and averaged Nusselt numbers are also presented. The results show that the relatively simple lattice Boltzmann method is a good alternative to the finite element method for analyzing non-Newtonian fluids.     

Electroviscous effect on non-Newtonian fluid flow in microchannels

Understanding non-Newtonian flow in microchannels is of both fundamental and practical significance for various microfluidic devices. A numerical study of non-Newtonian flow in microchannels combined with electroviscous effect has been conducted. The electric potential in the electroviscous force term is calculated by solving a lattice Boltzmann equation. And another lattice Boltzmann equation without derivations of the velocity when calculating the shear is employed to obtain flow field. The simulation of commonly used power-law non-Newtonian flow shows that the electroviscous effect on the flow depends significantly on the fluid rheological behavior. For the shear thinning fluid of the power-law exponent n < 1, the fluid viscosity near the wall is smaller and the electroviscous effect plays a more important role. And its effect on the flow increases as the ratio of the Debye length to the channel height increases and the exponent n decreases. While the shear thickening fluid of n > 1 is less affected by the electroviscous force, it can be neglected in practical applications.     

Multiple-relaxation-time lattice Boltzmann model for generalized Newtonian fluid flows

The generalized Newtonian fluid, as an important kind of non-Newtonian fluids, has been widely used in both science and engineering. In this paper, we present a multiple-relaxation-time lattice Boltzmann model for generalized Newtonian fluid, and validate the model through a detailed comparison with analytical solutions and some published results. The accuracy and stability of the present model are also studied, and compared with those of the popular single-relaxation-time lattice Boltzmann model. Finally, the limit and potential of the multiple-relaxation-time lattice Boltzmann model are briefly discussed.     

Non-Newtonian flow in microporous structures under the electroviscous effect

Single phase non-Newtonian microporous flow combined with the electroviscous effect is investigated in the pore-scale under conditions of various rheological properties and electrokinetic parameters. The lattice Boltzmann method is employed to solve both the electric potential field and flow velocity field. The simulation of commonly used power-law non-Newtonian flow shows that the electroviscous effect on the flow depends on both the fluid rheological behavior and pore surface area ratio significantly. For the shear thinning fluid with power-law exponent n < 1, the fluid viscosity near the wall is smaller and the electrovicous effect plays a more important role compared to the Newtonian fluid and shear thickening fluid. The high pore surface area ratio in the porous structure increases the electroviscous force and the induced flow resistance becomes important even to the flow of Newtonian and shear thickening fluids.     

Bingham fluid simulation with the incompressible lattice Boltzmann model

The Bingham fluid flow is numerically studied using the lattice Boltzmann method by incorporating the Papanastasiou exponential modification approach. The He–Luo incompressible lattice Boltzmann model is employed to avoid numerical instability usually encountered in non-Newtonian fluid simulations due to a strong non-linear relationship between the shear rate tensor and the rate-of-strain tensor. First, the value of the regularization parameter in Bingham fluid mimicking is analyzed and a method to determine the value is proposed. Then, the model is validated by pressure-driven planar channel flow and planar sudden expansion flow. The velocity profiles for the pressure-driven planar channel flow are in good agreement with analytical solutions. The calculated reattachment lengths for a 2:1 planar sudden expansion flow also agree well with the available data. Finally, the Bingham flow over a cavity is studied, and the streamlines and yielded/unyielded regions are discussed.     

Pressure-driven and electroosmotic non-Newtonian flows through microporous media via lattice Boltzmann method

The lattice Boltzmann method is developed to simulate the pressure-driven flow and electroosmotic flow of non-Newtonian fluids in porous media based on the representative elementary volume scale. The flow through porous media was simulated by including the porosity into the equilibrium distribution function and adding a non-Newtonian force term to the evolution equation. The non-Newtonian behavior is considered based on the Herschel–Bulkley model. The velocity results for pressure-driven non-Newtonian flow agree well with the analytical solutions. For the electroosmotic flow, the influences of porosity, solid particle diameter, power law exponent, yield stress and electric parameters are investigated. The results demonstrate that the present lattice Boltzmann model is capable of modeling non-Newtonian flow through porous media.     

An approach to control the spurious currents in a multiphase lattice Boltzmann method and to improve the implementation of initial condition

Multiphase lattice Boltzmann methods are known to generate spurious or parasitic currents at the fluid–fluid interfaces. This nonphysical phenomenon has to be avoided, or at least controlled, in order to achieve reliable solutions. In this article, a method to control these fictitious velocities via lattice refinement is proposed, which is based on interface thickness control for which both the spurious currents and the physical fluid–fluid interface thickness vanishes as the spatial resolution increases. It has been found that a proper interface thickness adjustment is required as the lattice refinement is applied, or an increase in spurious currents, instead of a reduction, can occur. By combining the new method with an appropriate multiphase flow initialization, the overall stability for high density O(1000) and viscosity O(100) ratios is greatly improved. Although this research has been conducted with a Rothman and Keller type lattice Boltzmann model, it is believed that other types of multiphase lattice Boltzmann models could benefit from the basic ideas underlying this research. Copyright © 2015 John Wiley & Sons, Ltd.     

Pore-Scale Simulation of Shear Thinning Fluid Flow Using Lattice Boltzmann Method

The present work attempts to identify the roles of flow and geometric variables on the scaling factor which is a necessary parameter for modeling the apparent viscosity of non-Newtonian fluid in porous media. While idealizing the porous media microstructure as arrays of circular and square cylinders, the present study uses multi-relaxation time lattice Boltzmann method to conduct pore-scale simulation of shear thinning non-Newtonian fluid flow. Variation in the size and inclusion ratio of the solid cylinders generates wide range of porous media with varying porosity and permeability. The present study also used stochastic reconstruction technique to generate realistic, random porous microstructures. For each case, pore-scale fluid flow simulation enables the calculation of equivalent viscosity based on the computed shear rate within the pores. It is observed that the scaling factor has strong dependence on porosity, permeability, tortuosity and the percolation threshold, while approaching the maximum value at the percolation threshold porosity. The present investigation quantifies and proposes meaningful correlations between the scaling factor and the macroscopic properties of the porous media.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

A variational approach to non-steady non-newtonian flow in a circular pipe

The transient response of a non-Newtonian power-law fluid to several assumed forms of pressure pulse in a circular tube is analysed by the semi-direct variational method of Kanntovorich. Velocity profiles are shown for several power-law indices, and by comparing the results for the Newtonian case with the exact solution given by Szymanski, it is observed that the results are good to 5%. More accurate solutions have been found for the case involving Newtonian fluid flow. New results are reported concerning the effect of a triangular pressure pulse on the development and transient response of the flow field of a non-Newtonian fluid.     

Self-similar solutions of non-Newtonian fluid boundary layer in MHD

The self-similar solutions of the boundary layer for a non-Newtonian fluid in MHD were considered in [1, 2] for a power-law velocity distribution along the outer edge of the layer and constant electrical conductivity through the entire flow. However, the MHD flows of many conducting media, which are solutions or molten metals, cannot be described by the MHD equations for non-Newtonian fluids.The self-similar solutions of the boundary layer for a non-Newtonian fluid without account for interaction with the electromagnetic field were studied in [3].In the following we present the self-similar solutions for the boundary layer of pseudoplastic and dilatant fluids with account for the interaction with an electromagnetic field for the case of a power-law velocity distribution along the outer edge of the layer, when the conductivity of the fluid is constant throughout the flow and the magnetic Reynolds number is small.Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 2, No. 6, pp. 77–82, 1967The author wishes to thank S. V. Fal'kovich for his interest in this study.     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

Lattice Boltzmann simulation of non-Newtonian flows past confined cylinders

A second-order lattice Boltzmann algorithm is used for Power-Law non-Newtonian flow simulation. The shear dependent behavior of the fluid is implemented through calculating the shear locally from the lattice distribution functions. A step by step verification procedure is taken to ensure the accuracy and the physical correctness of the numerical simulation. The flow past a series of tandem arrangement of two cylinders is computed in a confined domain. The effects of Reynolds number, the Power-Law index, and the distance between two cylinders on both the flow field and the drag coefficients of the cylinders are examined in detail.     

A modified lattice boltzmann method for herschel-bulkley fluids

As one kind of the most important non-Newtonian fluids, Herschel-Bulkley fluids have been applied in practical engineering widely. In order to improve the stability and accuracy in the simulation of Herschel-Bulkley fluids with lattice Boltzmann method (LBM), a modified lattice Boltzmann method was proposed. With Poiseuille flow as the example, the shear-thinning fluids and the shear-thickening fluids were used, respectively, to introduce the method in detail. The comparison between the velocity distributions and the analytical solutions demonstrated the feasibility of the modified method. Also the effect of the power law index, lattice nodes, and the region of the shear rate on the relative error was discussed. Then, the method was applied into the simulation of the cement paste flow in the 3D printing extruder. The streamline figure was obtained and then conducted the flow simulation by the modified method and multiple-relaxation-time lattice Boltzmann method (MRT-LBM), respectively; the comparison further proved the modified method was feasible for Herschel-Bulkley fluids.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

基于格子Boltzmann两相模型的粘弹流体挤出胀大分析

挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象,采用颜色模型模拟出口处粘弹流体和空气的两相流动,通过重新标色获得两种流体的界面,并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较,发现格子Boltzmann两相模型结果与解析解和实验结果相吻合,相比于单相模型,收敛速度更快,解的稳定性更高．研究了流道尺寸对胀大率的影响,并对挤出胀大的内在机理进行了分析．     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

Validation of the Local Thermal Equilibrium Assumption in Forced Convection of Non-Newtonian Fluids through Porous Channels

In this paper, we assess the validity of the local thermal equilibrium assumption in the non-Newtonian forced convection flow through channels filled with porous media. For this purpose, the problem is solved numerically using local thermal non-equilibrium and non-Darcian models. Numerical solutions obtained over broad ranges of representative dimensionless parameters are utilized to map conditions at which the local thermal equilibrium assumption can or cannot be employed. The circumstances of a higher modified Peclet number, a lower modified Biot number, a lower fluid-to-solid thermal conductivity ratio, a lower power-law fluid index, and a lower microscopic and macroscopic frictional flow resistance coefficients, are identified as unfavorable circumstances for the local thermal equilibrium (LTE) condition to hold. Quantitative LTE validity maps that reflect the proportional effect of each parameter as related to others are presented.