Preview

Lattice Bgk Model for Incompressible Navier–Stokes Equation

Powerful Essays
Open Document
Open Document
6405 Words
Grammar
Grammar
Plagiarism
Plagiarism
Writing
Writing
Score
Score
Lattice Bgk Model for Incompressible Navier–Stokes Equation
Journal of Computational Physics 165, 288–306 (2000) doi:10.1006/jcph.2000.6616, available online at http://www.idealibrary.com on

Lattice BGK Model for Incompressible Navier–Stokes Equation
Zhaoli Guo,∗ Baochang Shi,† and Nengchao Wang†
∗ National Laboratory of Coal Combustion, and Department of Computer Science, Huazhong University of Science and Technology, Wuhan 430074, People’s Republic of China; and †Department of Mathematics, Huazhong University of Science and Technology, Wuhan 430074, People’s Republic of China E-mail: sbchust@public.wuhan.cngb.com Received May 10, 1999; revised February 28, 2000

Most of the existing lattice Boltzmann BGK models (LBGK) can be viewed as compressible schemes to simulate incompressible fluid flows. The compressible effect might lead to some undesirable errors in numerical simulations. In this paper a LBGK model without compressible effect is designed for simulating incompressible flows. The incompressible Navier–Stokes equations are exactly recovered from this incompressible LBGK model. Numerical simulations of the plane Poiseuille flow, the unsteady 2-D shear decaying flow, the driven cavity flow, and the flow around a circular cylinder are performed. The results agree well with the analytic solutions and the results of previous studies. c 2000 Academic Press Key Words: Lattice BGK method; Incompressible Navier–Stokes equation.

1. INTRODUCTION

The Lattice Boltzmann BGK (LBGK) method is a new numerical scheme for simulating viscous compressible flows in the subsonic regime [2]. In recent years, LBGK has achieved great success in simulations of fluid flows and in modeling physics in fluids. Through multiscaling expansion [7], the compressible Navier–Stokes equations can be recovered from the lattice Boltzmann BGK equation on the assumptions that (i) the Mach number is small, and (ii) the density varies slowly. Therefore, theoretically the LBGK model can only be used to simulate compressible flows in the incompressible limit.



References: 1. Reference removed in proofs. 2. S. Chen and G. Doolen, Lattice Boltzmann method for fluid flows, Ann. Rev. Fluid Mech. 30, 329 (1998). 3. S. Chen, D. Martinez, and R. Mei, On boundary conditions in lattice Boltzmann methods, Phys. Fluids. 8, 2527 (1996). 4. M. Coutanceau and R. Bouard, Experimental determination of the main features of the viscous flow in the wake of a circular cylinder in uniform translation. 1. Steady flow, J. Fluid Mech. 79, 231 (1977). 5. W.-N. E, and J. Liu, Essential compact scheme for unsteady viscous incompressible flows, J. Comput. Phys. 126, 122 (1996). 6. O. Filippova and D. H¨ nel, Grid refinement for lattice-BGK models, J. Comput. Phys. 147, 219 (1998). a 7. U. Frisch, D. d’Humi´ res, B. Hasslacher, P. Lallemand, Y. Pomeau, and J.-P. Rivet, Lattice gas hydrodynamics e in two and three dimensions. Complex Syst. 1, 649 (1987). 8. U. Ghia, K. N. Ghia, and C. T. Shin, High-Re solutions for incompressible flow using the Navier–Stokes equations and a multigrid method, J. Comput. Phys. 48, 387 (1982). 9. X. He and L.-S. Luo, Lattice Boltzmann model for the incompressible Navier–Stokes equation, J. Stat. Phys. 88, 927 (1997). 10. X. He and G. Doolen, Lattice Boltzmann method on curvilinear coordinates system: Flow around a circular cylinder, J. Comput. Phys. 134, 306 (1997). 11. F. J. Higuera and S. Succi, Simulating the flow around a circular cylinder with a lattice Boltzmann equation, Europhys. Lett. 8, 517 (1989). 12. S. Hou and Q. Zou, Simulation of cavity flow by the lattice Boltzmann method, J. Comput. Phys. 118, 329 (1995). 306 GUO, SHI, AND WANG 13. Z. Lin, H. Fang, and R. Tao, Improved lattice Boltzmann model for incompressible two-dimensional steady flows, Phys. Rev. E 54, 6323 (1997). 14. D. O. Martinez, W. H. Matthaeus, S. Chen, and D. C. Montgomery, Comparison of spectral method and lattice Boltzmann simulations of two-dimensional hydrodynamics, Phys. Fluids. 6, 1285 (1994). 15. R. Mei and Q. Shyy, On the finite difference-based lattice Boltzmann method in curvilinear coordinates, J. Comput. Phys. 143, 426 (1998). 16. F. Nieuwstadt and H. B. Keller, Viscous flow past circular cylinders, Comput. Fluids. 1, 59 (1973). 17. Y. Qian, D. d’Humi´ res, and P. Lallemand, Lattice BGK models for Navier–Stokes equation, Europhys. Lett. e 17, 479 (1992). 18. R. Schreiber and H. Keller, Driven cavity flow by efficient numerical techniques, J. Comput. Phys. 49, 310 (1983). 19. S. P. Vanka, Block-implicit multigrid solution of Navier–Stokes equations in primitive variables, J. Comput. Phys. 65, 138 (1986). 20. L. Wagner, Pressure in lattice Boltzmann simulations of flow around a cylinder, Phys. Fluids. 6, 3516 (1994). 21. L. Wagner and F. Hayot, Lattice Boltzmann simulations of flow past a cylindrical obstacle, J. Stat. Phys. 81, 63 (1995). 22. Q. Zou, S. Hou, S. Chen, and G. Doolen, An improved incompressible lattice Boltzmann model for timeindependent flows, J. Stat. Phys. 81, 35 (1995).

You May Also Find These Documents Helpful

  • Good Essays

    Nt1310 Unit 3 Lab Report

    • 3999 Words
    • 16 Pages

    Numerical simulations are carried out employing a commercial CFD code (ANSYS CFX v11). The solver is using a full-scale time-marching 3D viscous model. Underlying equations, three dimensional Navier-Stokes equations in their conservation form, are being solved by using a Finite Volume method, where equations are integrated over the finite control volumes. Thereby, the solution domain is subdivided into a finite number of control volumes employing a suitable grid, which defines the control boundaries around a computational node in each control volume center. 6.1.1 Governing equations In fluid dynamics, the fluid flow is governed by the conservation laws for mass, momentum and energy. The basic conservation laws are formulated by using Leibniz-Reynolds transport theorem, which is an integral relation stating that the changes of some intensive property defined over a control volume must be equal to what is lost (or gained) through the boundaries of the volume plus what is created/consumed by sources and sinks inside the control volume. The…

    • 3999 Words
    • 16 Pages
    Good Essays
  • Satisfactory Essays

    Such these results confirm the findings that obtained recently by Kagumba and Al Mesfer where they reported that at churn turbulent flow regime, the gas holdup that was achieved in bubble column in the absence of vertical internals could be extrapolated to the columns with vertical internals when these columns with internals are operated at the same superficial gas velocity but calculated based on the free cross-sectional area for the flow. This is interesting finding due to they achieved this observation in 6-inch bubble columns while in this work the same observation can also be obtained in 18-inch bubble…

    • 262 Words
    • 2 Pages
    Satisfactory Essays
  • Better Essays

    4. Enhanced Single-Phase Turbulent Tube-side Flows and Heat Transfer. Engineering Data Book III. Wolverine Tube Inc. Web. 4 Apr. 2012.…

    • 1284 Words
    • 9 Pages
    Better Essays
  • Powerful Essays

    Syllabus

    • 2190 Words
    • 9 Pages

    Apply microscopic and macroscopic versions of gas laws to solve problems describing the behavior of gases.…

    • 2190 Words
    • 9 Pages
    Powerful Essays
  • Satisfactory Essays

    Kinetic Molecular Theory

    • 1472 Words
    • 6 Pages

    <i>1. State the five assumptions of the Kinetic-Molecular Theory of gases.</i><br><br>a) Gases consist of large numbers of tiny particles. These particles, usually molecules or atoms, typically occupy a volume about 1000 times larger than occupied by the same number of particles in the liquid or solid state. Thus molecules of gases are much further apart than those of liquids or solids.<br><br>Most of the volume occupied by a gas is empty space. This accounts for the lower density of gases compared to liquids and solids, and the fact that gases are easily compressible.<br><br>b) The particles of a gas are in constant motion, moving rapidly in straight lines in all directions, and thus passes kinetic energy. The kinetic energy of particles overcomes…

    • 1472 Words
    • 6 Pages
    Satisfactory Essays
  • Better Essays

    Hodanbosi, C. (1996). Pascal’s Principle and Hydraulics. Fairman, J. G. (Ed.) Retrieved from National Space and Aeronautics Administration, Glenn Research Center Web site: http://www.grc.nasa.gov/WWW/k-12/WindTunnel/Activities/Pascals_principle.html…

    • 1602 Words
    • 7 Pages
    Better Essays
  • Better Essays

    Ammonia Combustion

    • 2506 Words
    • 11 Pages

    catherineduynslaegher@uclouvain.be de Physico-Chimie de la Combustion, Place Louis Pasteur n° 1 de Thermodynamique et Turbomachines, Place du Levant n° 2 Université Catholique de Louvain - 1348 Louvain-la-Neuve - Belgium…

    • 2506 Words
    • 11 Pages
    Better Essays
  • Powerful Essays

    Drag on a Cylinder

    • 3180 Words
    • 13 Pages

    Munson, B. R., Young, D. F., Okiishi, T. H., & Huebsch, W. W. (2009). Fundamentals of fluid mechanics. (6th ed.). Jefferson City: John Wiley & Sons Inc.…

    • 3180 Words
    • 13 Pages
    Powerful Essays
  • Powerful Essays

    Vortex Induce Vibration

    • 8263 Words
    • 34 Pages

    H. M. Blackburn, R.N.G., D C. H. K. Williamson, A complementary numerical and physical investigation of vortex-induced vibration. Journal of Fluids and Structures, 2000. Vol. 15: p. 481-488.…

    • 8263 Words
    • 34 Pages
    Powerful Essays
  • Good Essays

    We here describe the implementation of the vortex lattice method for two-dimensional flows past thin airfoils. The method is even more useful for three-dimensional wings, i.e., for the flow past wings of finite span, but that problem is not considered here. Instead, the reader is referred to standard aerodynamics texts, e.g., [2]. In this numerical procedure to solve the thin-airfoil problem, we place a finite number of discrete vortices along the chord line, with the boundary condition that the induced vertical velocity dyc v= − α, (1) dx be enforced at selected control points to determine the vortex strengths. Equation (1) simply says that the net velocity vector, comprised of components due to the free stream, at angle of attack α to the chord line, plus that induced by the point vortices, is tangent to the camber line whose slope is dyc / dx; the magnitude of the free stream velocity is taken to be unity. Thus, we discretize the chord line into a finite number N of segments, or panels, as illustrated in Fig. 1 (a). On each panel we place a point vortex and a control point, as illustrated in Fig. 1 (b). The most accurate results are obtained by locating the vortex one-quarter of the panel length, and the control point three-quarters of the panel length, aft of the leading edge of the panel. (This strategy can be shown to reproduce the exact results of analytical…

    • 3124 Words
    • 19 Pages
    Good Essays
  • Better Essays

    In fluid dynamics, Bernoulli’s principles states that an increase in the speed of the fluid corresponds to a decrease in pressure of the same fluid. Similarly, the decrease in pressure corresponds to a loss in the potential energy of the fluid. The principle is applicable to various types of fluids, which leads to Bernoulli’s equation. There are different types of Bernoulli’s equation depending on the nature of fluid in use. The principle is valid for most compressible and incompressible fluids at low Mach numbers (less than 0.3). For this principle, fluid particles are only subject to pressure and its own weight.…

    • 1251 Words
    • 6 Pages
    Better Essays
  • Satisfactory Essays

    Kinetic Theory

    • 2037 Words
    • 3 Pages

    Kinetic theory (or the kinetic or kinetic-molecular theory of gases) is the theory that HYPERLINK http//en.wikipedia.org/wiki/Gasgases are made up of a large number of small particles (HYPERLINK http//en.wikipedia.org/wiki/Atomatoms or HYPERLINK http//en.wikipedia.org/wiki/Moleculemolecules), all of which are in constant, HYPERLINK http//en.wikipedia.org/wiki/Randomnessrandom HYPERLINK http//en.wikipedia.org/wiki/Motion_(physics)motion. The rapidly moving particles constantly collide with each other and with the walls of the container. Kinetic theory explains HYPERLINK http//en.wikipedia.org/wiki/Macroscopicmacroscopic properties of gases, such as pressure, temperature, or volume, by considering their molecular composition and motion. Essentially, the theory posits that pressure is due not to static repulsion between molecules, as was HYPERLINK http//en.wikipedia.org/wiki/Isaac_NewtonIsaac Newtons conjecture, but due to HYPERLINK http//en.wikipedia.org/wiki/Collisioncollisions between molecules moving at different velocities. While the particles making up a gas are too small to be visible, the jittering motion of pollen grains or dust particles which can be seen under a microscope, known as HYPERLINK http//en.wikipedia.org/wiki/Brownian_motionBrownian motion, results directly from collisions between the particle and air molecules. This experimental evidence for kinetic theory, pointed out by HYPERLINK http//en.wikipedia.org/wiki/Albert_EinsteinAlbert Einstein in 1905, is generally seen as having confirmed the existence of atoms and molecules. HYPERLINK http//en.wikipedia.org/wiki/Kinetic_theoryhttp//en.wikipedia.org/wiki/Kinetic_theory Postulates The theory for ideal gases makes the following assumptions The gas consists of very small particles, all with non-zero HYPERLINK http//en.wikipedia.org/wiki/Massmass. The number of molecules is large such that statistical treatment can be applied. These molecules are in constant, HYPERLINK…

    • 2037 Words
    • 3 Pages
    Satisfactory Essays
  • Powerful Essays

    We present a novel, incompressible fluid simulation method based on the Lagrangian Smoothed Particle Hydrodynamics (SPH) model. In our method, incompressibility is enforced by using a predictioncorrection scheme to determine the particle pressures. For this, the information about density fluctuations is actively propagated through the fluid and pressure values are updated until the targeted density is satisfied. With this approach, we avoid the computational expenses of solving a pressure Poisson equation, while still being able to use large time steps in the simulation. The achieved results show that our predictive-corrective incompressible SPH (PCISPH) method clearly outperforms the commonly used weakly compressible SPH (WCSPH) model by more than an order of magnitude while the computations are in good agreement with the WCSPH results. CR Categories: I.3.5 [Computer Graphics]: Computational Geometry and Object Modeling—Physically based modeling; I.3.7 [Computer Graphics]: Three-Dimensional Graphics and Realism— Animation. Keywords: fluid simulation, SPH, incompressibility…

    • 4886 Words
    • 20 Pages
    Powerful Essays
  • Powerful Essays

    Lab Report

    • 12052 Words
    • 49 Pages

    11. Welty, J.R., et al, “Fundamental of momentum, heat and mass transfer”, 4th ed. Wiley (2001)…

    • 12052 Words
    • 49 Pages
    Powerful Essays
  • Good Essays

    Revolutionary research on nanofluid has been done by scientists all around the world. The next generation of flow and heat transfer fluids, called Nanofluid, is equipped with appropriate characteristics of nanoparticles. The research on nanofluids has opened a wide spectrum of interest of physicists, engineers and chemists etc. Enhancing thermal conductivity of fluids using nanoparticles is introduced by Choi and Eastman [17]. Later, Kakac and Pramuanjaroenkij [18] analyzed enhancement of convective heat transfer with nano fluids. Ahmad et al. [19] studied Blasius and Sakiadis problems in nanofluids. Yohannes et al. [20] heat and mass transfer in a flow of nanofluids through porous media due to stretching sheet with viscous dissipation and chemical reaction. Khann et al. [21] investigated ferrofluid with three types of ferroparticles Ferro Oxide, Cobalt Ferrite and Mn-Zn ferrite in the presence of viscous dissipation. Mohsen Sheikholeslami et al. [22] studied two phase simulation of nanofluid and heat and mass transfer in the presence of viscous dissipation. Boundary layer flow of a nanofluid with…

    • 806 Words
    • 4 Pages
    Good Essays

Related Topics