Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 72, 2017 - Issue 6
290
Views
0
CrossRef citations to date
0
Altmetric
Original Articles

An improved MRT-LBM for Herschel–Bulkley fluids with high Reynolds number

, , &
Pages 409-420 | Received 13 Sep 2017, Accepted 14 Nov 2017, Published online: 15 Dec 2017
 

ABSTRACT

We proposed a truncated multiple-relaxation-time lattice Boltzmann method (MRT-LBM) for the Herschel–Bulkley flow with high Reynolds number. The effectiveness of the proposed model was verified by Poiseuille flow. We analyzed the effect of the initial yield stress and the relative errors with diverse lattice nodes in case of different power-law indexes. Finally, the improved MRT-LBM was used to simulate the lid-driven flow of high Reynolds numbers. The streamlines prove not only the effectiveness of MRT-LBM in simulating the flow with high Reynolds number, but also the feasibility of the method in dealing with the yielding behavior of Herschel–Bulkley fluids.

Nomenclature

Symbols=

Meaning

c=

lattice speed

cs=

sound speed

ei=

discrete velocities

f=

distribution function

h=

width of the channel

H=

distance between two plates

K=

viscosity coefficient

L=

pitch of the single screw

m=

related to the increase in the stress

M=

transformation matrix

n=

power-law index

N=

lattice nodes

Ns=

rotation speed

r=

displacement vector

=

main diagonal matrix related to the relaxation process

Sαβ=

strain-rate tensor

t=

Time

u=

velocity vector

v vx vz=

Velocities

W=

depth of the channel

yl yh yτ=

critical values

α1 α2 α3=

constant terms

=

shearing rate

δx=

lattice spacing

δt=

time spacing

θ=

lead angle

μ=

dynamic viscosity

μB=

viscosity coefficient of Bingham fluids

νap ν νmin νmax=

related kinematic viscosity

ρ=

Density

τ=

relaxation time

τa=

τ0 stress

ωi=

weight coefficient

∇P=

pressure gradient

Nomenclature

Symbols=

Meaning

c=

lattice speed

cs=

sound speed

ei=

discrete velocities

f=

distribution function

h=

width of the channel

H=

distance between two plates

K=

viscosity coefficient

L=

pitch of the single screw

m=

related to the increase in the stress

M=

transformation matrix

n=

power-law index

N=

lattice nodes

Ns=

rotation speed

r=

displacement vector

=

main diagonal matrix related to the relaxation process

Sαβ=

strain-rate tensor

t=

Time

u=

velocity vector

v vx vz=

Velocities

W=

depth of the channel

yl yh yτ=

critical values

α1 α2 α3=

constant terms

=

shearing rate

δx=

lattice spacing

δt=

time spacing

θ=

lead angle

μ=

dynamic viscosity

μB=

viscosity coefficient of Bingham fluids

νap ν νmin νmax=

related kinematic viscosity

ρ=

Density

τ=

relaxation time

τa=

τ0 stress

ωi=

weight coefficient

∇P=

pressure gradient

Additional information

Funding

This work was supported by National Natural Science Foundation of China (Grant No. 51635003) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX17_0938).

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.