A Direct Effective Viscosity Approach for Modeling and Simulating Bingham Fluids with the Cumulant Lattice Boltzmann Method

Modeling of fluids with complex rheology in the lattice Boltzmann method (LBM) is typically realized through the introduction of an effective viscosity. For fluids with a yield stress behavior, such as so-called Bingham fluids, the effective viscosity has a singularity for low shear rates and may become negative. This is typically avoided by regularization such as Papanastasiou’s method. Here we argue that the effective viscosity model can be re-interpreted as a generalized equilibrium in which no violation of the stability constraint is observed. We implement a Bingham fluid model in a three-dimensional cumulant lattice Boltzmann framework and compare the direct analytic effective viscosity/generalized equilibrium method to the iterative approach first introduced by Vikhansky which avoids the singularity in viscosity that can arise in the analytic method. We find that both methods obtain similar results at coarse resolutions. However, at higher resolutions the accuracy of the regularized method levels off while the accuracy of the direct method continuously improves. We find that the accuracy of the proposed direct method is not limited by the singularity in viscosity indicating that a regularization is not strictly necessary.

Analysis of Bingham fluid radial flow in smooth fractures

Solutions for radial flow of a Bingham fluid are analyzed in this paper.It aims to eliminate confusions in the literature concerning the plug flow region in different solutions for analysis and design of grouting in rock fractures.The analyses based on the force balance equation reveal that the plug flow region in Bingham radial flow is independent of the fracture radius,and is not a growth function adapted from the solution of one-dimensional(1D)slit flow according to‘similarity’.Based on the shear stress distribution,we analytically proposed that a non-uniform plug flow region cannot exist.The Bingham fluid(grout)penetration and flowrate evolution as functions of grouting time are given using the correct expression for the plug flow region.The radius-independent plug flow region and the presented flowrate evolution equation are also verified numerically.For radial flow,the relative penetration length is equal to the relative width of plug flow region,which is the same as that for 1D channel flow.Discrepancies in analytical solutions for grout penetration and flowrate evolution were also illustrated.The clarification of the plug flow region and evaluation of discrepancies in analytical solutions presented in this work could simplify modeling and design of grouting in rock engineering applications.

Discussion on “Analysis of Bingham fluid radial flow in smooth fractures”[J Rock Mech Geotech Eng 12(2020) 1112e1118]

Recently,Zou et al.(2020 a)published a theoretical analysis on the radial flow of a Bingham fluid,where they argued that the classical analysis by Dai and Bird(1981)violates the mass conservation.The present discussion aims to clarify this conflict between those two studies.It is noted that Zou et al.(2020 a)presumed the gap-wise mass flux is negligible in the mass conservation equation,while Dai and Bird(1981)did not require so in their model,and this is found to be the origin of the conflict.In fact,Dai and Bird(1981)’s model is shown to not violate the mass conservation.Therefore,those two models should be viewed as separate models derived from different perspectives.Details of the major difference between the two models are discussed.

Heat transfer and Helmholtz-Smoluchowski velocity in Bingham fluid flow

A mathematical study is developed for the electro-osmotic flow of a nonNewtonian fluid in a wavy microchannel in which a Bingham viscoplastic fluid model is considered.For electric potential distributions,a Poisson-Boltzmann equation is employed in the presence of an electrical double layer(EDL).The analytical solutions of dimensionless boundary value problems are obtained with the Debye-Huckel theory,the lubrication theory,and the long wavelength approximations.The effects of the Debyelength parameter,the plug flow width,the Helmholtz-Smoluchowski velocity,and the Joule heating on the normalized temperature,the velocity,the pressure gradient,the volumetric flow rate,and the Nusselt number for heat transfer are evaluated in detail using graphs.The analysis provides important findings regarding heat transfer in electroosmotic flows through a wavy microchannel.

Motion Law of a Sphere in Bingham Fluid with Thixotropy

The thixotropy properties and the motion law of a sphere in the Bingham fluid have been studied. Through observation of the settling motion of a single sphere in the Bingham fluid on the X-ray screen, it has been discovered that the mud in estuaries and along sea bay, and the hyperconcentrated flow all behave as the Bingham fl fluid with thixotropy properties as the large sediment concentration. Through derivation, the theoretical relationship between the yield stress and non-settling maximum sphere supported by the stress for the Bingham fluid has been developed, the equations for calculating the increasing yield stress and the non-settling maximum sphere diameter with the duration at rest of the slurry have been obtained. In consideration of the effect of thixotropy on fluid motion, the Navier-Stokes equation group for the Bingham thixotropy fluid has been developed. Through further study of the flow boundary condition of settling motion of ii single sphere in the Bingham thixotropy fluid, and the solving of the Navier-Stokes equation group, under the small Reynolds number, the theoretical equation of the drag force of the Bingham thixotropy fluid flowing around a sphere has been deduced. The theoretical relationship between drag coefficient and Reynolds number has been derived. By use of the experimental data of rheological test of various slurries measured with viscometer and those of single sphere motion observed on the X-ray screeen, the above equations have been verified. The equations are in good agreement with the experimental data for various slurries.

Effects of particle fractions on the Bingham yield stress and viscosity of fine-coarse particle suspensions

The rheological behaviors of highly concentrated fine particle suspensions(clay-silt-water mixtures)and coarse particle suspensions(coarse particles within a fine particle suspension)were investigated in this study.Experimental results demonstrated that the Bingham Fluid Model with two rheological parameters,Bingham yield stressand viscosity,wellcharacterized the rheological behavior of fine particle suspensions at shear rates between 4 and 20 s^(-1).The inclusion of coarse particles within a fine particle suspension induced an enhancement to the rheological parameters.The rheological parameters of a coarse particle suspension not only depend on its total particle fraction but also on its relative fine/coarse particle fractions.Empirical equations of these two parameters were proposed,quantitatively related to both fine and coarse particle fractions.Results indicated that the Bingham yield stress and viscosity are much more(an order larger)sensitive to changes in fine particle fraction than to changes in coarse particle fraction.

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.

Reply to Discussion on “Analysis of Bingham fluid radial flow in smooth fractures”

Recently,Hoang et al.(2021)discussed our paper Zou et al.(2020).In our paper,we made a statement that Dai and Bird(1981)’s solution for two-dimensional(2 D)radial Bingham fluid flow between parallel plates violates mass balance.Hoang et al.pointed out that Dai and Bird(1981)’s solution does not violate the mass balance because Dai and Bird(1981)’s solution and our analysis are based on different assumptions,i.e.with consideration of the vertical velocity component in the continuity equation or not,which leads to two different approximation models.In this sense,the mass balance of Dai and Bird(1981)’s solution should not be checked using our solution as a reference.In this reply,we add remarks on the two approximation models and their implication for rock grouting analysis.The discussion by Hoang et al.and this reply are helpful to thoroughly eliminate the existing confusion regarding the two solutions in the rock grouting research community.

A NUMERICAL STUDY OF BINGHAM TURBULENT FLOW IN SUDDEN-EXPANSION STRAIGHT CIRCULAR PIPE

The control equations for the turbulent flow of Bingham fluid are established according to Bingham fluid constitution equation. Pressure field and velocity field are correlted by pressure-correction equation. The numerical computations are performed on Bingham fluid turbulent flow in sudden-expansion straight circular pipe, and the flow mechanisms are discussed.

The Analysis of Stability of Bingham Fluid Flowing Down an Inclined Plane

In this paper, the stability problem of Bingham fluids flowing down an inclinedplane is studied with respect to two dimensional disturbances. The crilical Reynolodsnumber is given in ihe case of long waves, and the effect of yield stress on stability isanalysed.

THE ANALYSIS OF STABILITY OF BINGHAM FLUID FLOWIGN DOWN AN INCLINED PLANE

In this paper the stability problem of Bingham flowing down an inclinedplane is studied with respect to two dimensional disturbances, The critical Reynolodsnumber is given in the. case. of long wayes and the effect of yield stress on stability isanalysed.

AN IMPLICIT SOLUTION OF BI-PENALTY APPROXIMATION WITH ORTHOGONALITY PROJECTION FOR THE NUMERICAL SIMULATION OF BINGHAM FLUID FLOW

An implicit algorithm of Bi-penalty approximation with orthogonality projection for the numerical simulation of Bingham fluid flow problems is proposed in this paper. A Newton fluid flow with two kinds of artificial viscosity subjected to the inequality constraint is introduced to approximate the Bingham fluid flow. This approach can effectively simulate the Bingham fluid flow with floating rigid cores or fixing rigid cores.

DESIGN METHOD OF MAGNETORHEOLOGICAL FLUID SHOCK ABSORBER FOR CAR SUSPENSION

The Bingham constitutive model, which is previously used in depiction of magnetorheological （MR） fluids rheological behaviors for design devices, exhibits discontinuous characteristics in representation of pre-yield behaviors and post-yield behaviors. A Biviscous constitutive model is presented to depict rheological behaviors of MR fluids and design automotive shock absorber. Quasi-static flow equations of MR fluids in annular channels are set theoretically up based on Navier-Stokes equations and several rational simplifications are made. And both flow boundary conditions and flow compatibilities conditions are established. Meantime, analytical velocity profiles of MR fluids though annular channels are obtained via solution of the quasi-static flow equations using Biviscous constitutive model. The prediction methodology of damping force offered by MR fluid shock absorber is formulated and damping performances are predicated in order to determine design parameters. MR fluid shock absorber for Mazda 323 car suspension is designed and fabricated in Chongqing University, China. Measurements from sinusoidal displacement cycle by Shanchuan Shock Absorber Ltd. of China North Industry Corporation reveal that the analytical methodology and design theory are reasonable.

A Review of Studies on the Interaction Between Waves and Muddy Bottom

This paper presents a brief review of the results on the interaction between waves and muddy bottom obtained during the last decade including the results obtained by the author at the Coastal Engineering Laboratory of Tianjin University.

Simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

Fresh cement mortar is a type of workable paste, which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering. In this paper, Papanastasiou＇s model for Bingham fluids is solved by using the multiple- relaxation-time lattice Boltzmann model （MRT-LB）. Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou＇s model provides a good approximation of realistic Bingham plastics for values of m 〉 108. For lower values of m, Papanastasiou＇s model is valid for fluids between Bingham and Newtonian fluids. The MRT-LB model is validated by two benchmark problems： 2D steady Poiseuille flows and lid-driven cavity flows. Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability. We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle. Besides the rich flow structures obtained in this work, the dynamics fhi d force on the round particle is calculated. Results show that both the Reynolds number Re and the Bingham number Bn affect the drag coefficients Co, and a drag coefficient with Re and Bn being taken into account is proposed. The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed. Finally, the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields. These results help simulation of fresh concrete flowing in porous media.

Bingham Flows in Periodic Domains of Infinite Length

The Bingham fluid model has been successfully used in modeling a large class of non-Newtonian fluids. In this paper, the authors extend to the case of Bingham fluids the results previously obtained by Chipot and Mardare, who studied the asymptotics of the Stokes flow in a cylindrical domain that becomes unbounded in one direction, and prove the convergence of the solution to the Bingham problem in a finite periodic domain, to the solution of the Bingham problem in the infinite periodic domain, as the length of the finite domain goes to infinity. As a consequence of this convergence, the existence of a solution to a Bingham problem in the infinite periodic domain is obtained, and the uniqueness of the velocity field for this problem is also shown. Finally, they show that the error in approximating the velocity field in the infinite domain with the velocity in a periodic domain of length 2l has a polynomial decay in , unlike in the Stokes case （see [Chipot, M. and Mardare, S., Asymptotic behaviour of the Stokes problem in cylinders becoming unbounded in one direction, Journal de Mathgmatiques Pures et Appliqudes, 90（2）, 2008, 133-159]） where it has an exponential decay. This is in itself an important result for the numerical simulations of non-Newtonian flows in long tubes.

Effects of heat transfer on the peristaltic MHD flow of a Bingham fluid through a porous medium in a channel

In this paper, we study the effects of heat transfer on the peristaltic magneto- hydrodynamic （MHD） flow of a Bingham fluid through a porous medium in a channel. Long wavelength approximation （that is, the wavelength of the peristaltic wave is large in comparison with the radius of the channel） and low Reynolds number are used to linearize the governing equations. The velocity field for the model of interest is solved by Adomian decomposition method. The expressions for pressure rise, flow rate and frictional force are obtained. The effect of magnetic field, Darcy number, yield stress, amplitude ratio and the temperature on the axial pressure gradient, pumping charac- teristics and frictional force are discussed through graphs.

Asymptotic Analysis of a Bingham Fluid in a Thin Domain with Fourier and Tresca Boundary Conditions

In this paper we prove first the existence and uniqueness results for the weak solution,to the stationary equations for Bingham fluid in a three dimensional bounded domain with Fourier and Tresca boundary condition;then we study the asymptotic analysis when one dimension of the fluid domain tend to zero.The strong convergence of the velocity is proved,a specific Reynolds limit equation and the limit of Tresca free boundary conditions are obtained.