All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.

Related Documents

Share

Description

Lattice Boltzmann Simulation of Non-newtonian Fluid Flow in a Lid Driven Cavit

Transcript

International Journal of Mechanical Engineering and Technology (IJMET), ISSN 0976 – 6340(Print), ISSN 0976 – 6359(Online), Volume 5, Issue 8, August (2014), pp. 20-33 © IAEME
20
LATTICE BOLTZMANN SIMULATION OF NON-NEWTONIAN FLUID FLOW IN A LID DRIVEN CAVITY
M. Y. Gokhale
1
, Ignatius Fernandes
2 1
Department of Mathematics, Maharashtra Institute of Technology, Pune, India
2
Department of Mathematics, Rosary College, Navelim Goa, India
ABSTRACT
Lattice Boltzmann Method (LBM) is used to simulate the lid driven cavity flow to explore the mechanism of non-Newtonian fluid flow. The power law model is used to represent the class of non-Newtonian fluids (shear-thinning and shear-thickening fluids) by considering a range of 0.8 to 1.6. Investigation is carried out to study the influence of power law index and Reynolds number on the variation of velocity profiles and streamlines plots. Velocity profiles and the streamline patterns for various values of power law index at Reynolds numbers ranging 100 to 3200 are presented. Half way bounce back boundary conditions are employed in the numerical method. The LBM code is validated against the results taken from the published sources for flow in lid driven cavity and the results show fine agreement with established theory and the rheological behavior of the fluids.
Keywords:
Lid Driven Cavity, Non-Newtonian Fluids, Power Law, Lattice Boltzmann Method.
1. INTRODUCTION
Non-Newtonian fluid flow is an important subject in various natural and engineering processes which include applications in packed beds, petroleum engineering and purification processes. A wide range of research is available for Newtonian and non-Newtonian fluid flow in these areas. In general, the hydrodynamics of non-Newtonian fluids is much complex compared to that of their Newtonian counterpart because of the complex rheological properties. An important factor in understanding the mechanism of non-Newtonian fluids is to identify a local profile of non-Newtonian properties corresponding to the shear rate. Power law model is generally used to represent a class of non-Newtonian fluids which are inelastic and exhibit time independent shear stress. Though analytical solutions for the flow of power law fluids through simple geometry is available, computational approach becomes unavoidable in most of the situations particularly if the
INTERNATIONAL JOURNAL OF MECHANICAL ENGINEERING AND TECHNOLOGY (IJMET)
ISSN 0976 – 6340 (Print) ISSN 0976 – 6359 (Online) Volume 5, Issue 8, August (2014), pp. 20-33 © IAEME: www.iaeme.com/IJMET.asp Journal Impact Factor (2014): 7.5377 (Calculated by GISI) www.jifactor.com
IJMET
© I A E M E
International Journal of Mechanical Engineering and Technology (IJMET), ISSN 0976 – 6340(Print), ISSN 0976 – 6359(Online), Volume 5, Issue 8, August (2014), pp. 20-33 © IAEME
21
flow field is not one dimensional. Over the past several years, various computational methods have been applied to simulate the power law fluid flows in different geometries. Among different geometries, the lid driven cavity flow is considered to be one of the benchmark fluid problems in computational fluid dynamics. With recent advances in mathematical modeling and computer technology, lattice Boltzmann method (LBM) has been developed as an alternative approach to common numerical methods which are based on discretization of macroscopic continuum equations. LBM is effective for investigating the local non-Newtonian properties since important information to non-Newtonian fluids can be locally estimated. The fundamental idea of this method is to construct simplified kinetic models that include the essential physics of mesoscopic processes so that the macroscopic averaged properties obey the desired macroscopic properties [2]. In recent years, the problem of lid driven cavity flow has been widely used to understand the behavior of non-Newtonian fluid flow using power law model. Patil et al [5] applied the LBM for simulation of lid-driven flow in a two-dimensional, rectangular, deep cavity. They studied the location and strength of the primary vortex, the corner-eddy dynamics and showed the existence of corner eddies at the bottom, which come together to form a second primary-eddy as the cavity aspect-ratio is increased above a critical value. Bhaumik et al. [6] investigated lid-driven swirling flow in a confined cylindrical cavity using LBM by studying steady, 3-dimensional flow with respect to height-to-radius ratios and Reynolds numbers using the multiple-relaxation-time method. Nemati et al [7] applied LBM to investigate the mixed convection flows utilizing nanofluids in a lid-driven cavity. They investigated a water-based nanofluid containing Cu, CuO or Al
2
O
3
nanoparticles and the effects of Reynolds number and solid volume fraction for different nanofluids on hydrodynamic and thermal characteristics. Mendu et al. [9] used LBM to simulate non-Newtonian power law fluid flows in a double sided lid driven cavity. They investigated two different cases-parallel wall motion and anti-parallel wall motion of two sided lid driven cavity and studied the influence of power law index
)(
n
and Reynolds number
(Re)
on the variation of velocity and center of vortex location of fluid with the help of velocity profiles and streamline plots. In another paper, Mendu et al. [10] applied LBM to simulate two dimensional fluid flows in a square cavity driven by a periodically oscillating lid. Yang et al. [25] investigated the flow pattern in a two-dimensional lid-driven semi-circular cavity based on multiple relaxation time lattice Boltzmann method (MRT LBM) for Reynolds number ranging from 5000 to 50000. They showed that, as Reynolds number
increases, the flow in the cavity undergoes a complex transition. Erturk [8] discussed, in detail, the 2-D driven cavity flow problem by investigating the incompressible flow in a 2-D driven cavity in terms of physical, mathematical and numerical aspects, together with a survey on experimental and numerical studies. The paper also presented very fine grid steady solutions of the driven cavity flow at very high Reynolds numbers. The application of LBM to non-Newtonian fluid flow has been aggressively intensified in last few decades. Though, the problem of fluid flow in lid driven cavity has been studied rigorously, most of these studies have been confined to either laminar fluid flow or Newtonian fluids. The present paper uses LBM to simulate the lid driven cavity flow to explore the mechanism of non-Newtonian fluid flow which is laminar as well as under transition for a wide range of shear-thinning and shear-thickening fluids. The power law model is used to represent the class of non-Newtonian fluids (shear-thinning and shear-thickening). The influence of power law index
)(
n
and Reynolds number
(Re)
on the variation of velocity and center of vortex location of fluid with the help of velocity profiles and streamline plots is studied.
International Journal of Mechanical Engineering and Technology (IJMET), ISSN 0976 – 6340(Print), ISSN 0976 – 6359(Online), Volume 5, Issue 8, August (2014), pp. 20-33 © IAEME
22
Fig.1:
Geometry of the lid driven cavity.
2. BACKGROUND AND PROBLEM FORMULATION 2.1. Governing equations
In continuum domain, fluid flow is governed by Navier-Stokes (NS) equations along with the continuity equation. For incompressible, two-dimensional flow, the conservative form of the NS equations and the continuity equation can be written in Cartesian system as [1].
∂∂∂∂+
∂∂∂∂+∂∂−=∂∂+∂∂+∂∂
yu y xu x x p yvu xuut u
µ µ ρ ρ
(1)
∂∂∂∂+
∂∂∂∂+∂∂−=∂∂+∂∂+∂∂
yv y xv x x p yvv xuvt v
µ µ ρ ρ
(2)
0)()(
=∂∂+∂∂
yv xu
ρ ρ
(3) The two-dimensional lid driven square cavity with the top wall moving from left to right with a uniform velocity
1.0
0
==
uU
is considered, as shown in Fig 1. The left, right and bottom walls are kept stationary i.e. velocities at all other nodes are set to zero. We consider fluid to be non-Newtonian represented by power law model.
2.2. Physical boundary conditions
Top moving lid :
0
),(
u H xu
=
and
0),(
=
H xv
(4a) Bottom stationary side :
0)0,(
=
xu
and
0)0,(
=
xv
(4b) Left stationary side:
0),0(
=
yu
and
0),0(
=
yv
(4c) Right stationary side:
0),(
=
y Lu
and
0),(
=
y Lv
(4d)
International Journal of Mechanical Engineering and Technology (IJMET), ISSN 0976 – 6340(Print), ISSN 0976 – 6359(Online), Volume 5, Issue 8, August (2014), pp. 20-33 © IAEME
23
3. NUMERICAL METHOD AND FORMULATION 3.1. Lattice Boltzmann Method
In the present study, we cover incompressible fluid flows and a nine-velocity model on a two-dimensional lattice (D2Q9). The lattice Boltzmann method can be used to model hydrodynamic or mass transport phenomena by describing the particle distribution function
),(
t x f
i
giving the probability that a fluid particle with velocity
i
e
enters the lattice site
x
at a time
t
[1]. The subscript
i
represents the number of lattice links and
0
=
i
corresponds to the particle at rest residing at the center. The evolution of the particle distribution function on the lattice resulting from the collision processes and the particle propagation is governed by the discrete Boltzmann equation [11, 12, 13, 14].
),(),(),(
t xt x f dt t dt e x f
iiii
Ω=−++
8,....,1,0
=
i
(5) where
dt
is the time step and
i
Ω
is the collision operator which accounts for the change in the distribution function due to the collisions. The Bhatnagar-Gross-Krook (BGK) model [19] is used for the collision operator
[ ]
),(),(
1),(
t x f t x f t x
eqiii
−−=Ω
τ
8,.....,1,0
=
i
(6) where
τ
is the relaxation time and is related to the kinematic viscosity
υ
by
−=
21
2
τ υ
dt c
s
(7) Here
s
c
is the sound speed expressed by
( )
3 / 3 /
cdt dxc
s
==
(
c
is the particle speed and
dx
is the lattice spacing).
),(
t x f
eqi
, in equation (6), is the corresponding equilibrium distribution function for D2Q9 given by
( )
−
++=
),().,((
21),(.(
21),(.(
11,),(
2242
t xut xu
ct xue
ct xue
ct xwt x f
sisisieqi
ρ
(8a) where
),(
t xu
is the velocity and
i
w
is the weight coefficient with values
====
8,7,6,536 / 1
4,3,2,19 / 1
09 / 4
iiiw
i
(8b) Local particle density
),(
t x
ρ
and local particle momentum
u
ρ
are given by
∑
=
=
80
),(),(
ii
t x f t x
ρ
and
∑
=
=
80
),(),(
iii
t x f et xu
ρ
. (9)

We Need Your Support

Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks