0
Research Papers

Lattice Boltzmann Simulation of Airflow and Heat Transfer in a Model Ward of a Hospital

[+] Author and Article Information
Md. Farhad Hasan, Taasnim Ahmed Himika

Department of Mathematics & Physics,
North South University,
Dhaka 1229, Bangladesh

Md. Mamun Molla

Department of Mathematics & Physics,
North South University,
Dhaka 1229, Bangladesh
e-mail: mamun.molla@northsouth.edu

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF THERMAL SCIENCE AND ENGINEERING APPLICATIONS. Manuscript received January 30, 2016; final manuscript received July 14, 2016; published online November 2, 2016. Assoc. Editor: Giulio Lorenzini.

J. Thermal Sci. Eng. Appl 9(1), 011011 (Nov 02, 2016) (8 pages) Paper No: TSEA-16-1028; doi: 10.1115/1.4034817 History: Received January 30, 2016; Revised July 14, 2016

In this research, a very popular alternative computational technique, the lattice Boltzmann method (LBM), has been used to simulate the indoor airflow and heat transfer in a model hospital ward. Different Reynolds numbers have been used to study the airflow pattern. Boundary conditions for velocity and temperature have also been discussed in detail. Several tests have been conducted for code validation. LBM is demonstrated through simulation in forced convection inside hospital ward with six beds for two different situations: ward without partition and ward with partition. Changes in average rate of heat transfer in terms of average Nusselt numbers have also been recorded for those situations. Average Nusselt numbers were found to differ for different cases. In terms of airflow, it has been found that, for various Reynolds numbers, airflow changes its pattern and leads to few recirculations for relatively higher Reynolds number but remains steady for low Reynolds number. It was observed that partition narrowed the channel for airflow and once the air overcame this barrier, it gets free space and recirculation appears more. For higher Reynolds number, the average rate of heat transfer increases and patients near the recirculation zone release maximum heat and will feel more comfortable.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Hasslacher, B. , Pomeau, Y. , and Frisch, U. , 1986, “ Lattice-Gas Automata for the Navier–Stokes Equation,” Phys. Rev. Lett., 56(14), p. 1505. [CrossRef] [PubMed]
Peng, Y. , Shu, C. , and Chew, Y. T. , 2003, “ Simplified Thermal Lattice Boltzmann Model for Incompressible Thermal Flows,” Phys. Rev. E, 68(2), p. 026701. [CrossRef]
Szucki, M. , and Suchy, J. S. , 2012, “ A Method for Taking Into Account Local Viscosity Changes in Single Relaxation Time the Lattice Boltzmann Model,” Metall. Foundry Eng., 38(1), pp. 33–42. [CrossRef]
Lee, T. S. , Huang, H. , and Shu, C. , 2005, “ An Axisymmetric Incompressible Lattice BGK Model for Simulation of the Pulsatile Flow in a Circular Pipe,” Int. J. Numer. Meth. Fluids, 49(1), pp. 99–116. [CrossRef]
Sidik, N. A. C. , 2008, “ The Development of Thermal Lattice Boltzmann Models in Incompressible Limit,” Malays. J. Fundam. Appl. Sci., 3(2), pp. 193–202.
Jeong, H. K. , Yoon, H. S. , Ha, M. Y. , and Tsutahara, M. , 2010, “ An Immersed Boundary-Thermal Lattice Boltzmann Method Using an Equilibrium Internal Energy Density Approach for the Simulation of Flows With Heat Transfer,” J. Comp. Phys., 229(7), pp. 2526–2543. [CrossRef]
Mondal, B. , and Li, X. , 2010, “ Effect of Volumetric Radiation on Natural Convection in a Square Cavity Using Lattice Boltzmann Method With Non-Uniform Lattices,” Int. J. Heat Mass Transfer, 53(21), pp. 4935–4948. [CrossRef]
Ma, Y. , Bhattacharya, A. , Kuksenok, O. , Perchak, D. , and Balazs, A. C. , 2012, “ Modeling the Transport of Nanoparticle-Filled Binary Fluids Through Micropores,” Langmuir, 28(31), pp. 11410–11421. [CrossRef] [PubMed]
Chen, S. , and Doolen, G. D. , 1998, “ Lattice Boltzmann Method for Fluid Flows,” Ann. Rev. Fluid Mech., 30(1), pp. 329–364. [CrossRef]
Breyiannis, G. , and Valougeorgis, D. , 2004, “ Lattice Kinetic Simulations in Three-Dimensional Magnetohydrodynamics,” Phys. Rev. E, 69(6), p. 065702. [CrossRef]
Halliday, I. , Hammond, L. A. , Care, C. M. , Good, K. , and Stevens, A. , 2001, “ Lattice Boltzmann Equation Hydrodynamics,” Phys. Rev. E, 64(1), p. 011208. [CrossRef]
Azwadi, C. S. N. , and Tanahashi, T. , 2008, “ Development of 2-D and 3-D Double-Population Thermal Lattice Boltzmann Models,” Matematika, 24, pp. 53–66.
Crouse, B. , Krafczyk, M. , Kühner, S. , Rank, E. , and Van Treeck, C. , 2002, “ Indoor Air Flow Analysis Based on Lattice Boltzmann Methods,” Energy Build., 34(9), pp. 941–949. [CrossRef]
Zhang, S. J. , and Lin, C. X. , 2010, “ Application of Lattice Boltzmann Method in Indoor Airflow Simulation,” HVAC R Res., 16(6), pp. 825–841. [CrossRef]
Liu, J. , Wang, H. , and Wen, W. , 2009, “ Numerical Simulation on a Horizontal Airflow for Airborne Particles Control in Hospital Operating Room,” Build. Environ., 44(11), pp. 2284–2289. [CrossRef]
Chen, Q. , Lee, K. , Mazumdar, S. , Poussou, S. , Wang, L. , Wang, M. , and Zhang, Z. , 2010, “ Ventilation Performance Prediction for Buildings: Model Assessment,” Build. Environ., 45(2), pp. 295–303. [CrossRef]
Chen, Q. , 2009, “ Ventilation Performance Prediction for Buildings: A Method Overview and Recent Applications,” Build. Environ., 44(4), pp. 848–858. [CrossRef]
Wang, L. L. , and Chen, Q. , 2008, “ Evaluation of Some Assumptions Used in Multizone Airflow Network Models,” Build. Environ., 43(10), pp. 1671–1677. [CrossRef]
Zhang, Z. , Chen, X. , Mazumdar, S. , Zhang, T. , and Chen, Q. , 2009, “ Experimental and Numerical Investigation of Airflow and Contaminant Transport in an Airliner Cabin Mockup,” Build. Environ., 44(1), pp. 85–94. [CrossRef]
Guo, Z. , Shi, B. , and Zheng, C. , 2002, “ A Coupled Lattice BGK Model for the Boussinesq Equations,” Int. J. Numer. Methods Fluids, 39(4), pp. 325–342. [CrossRef]
Chen, S. , Chen, H. , Martnez, D. , and Matthaeus, W. , 1991, “ Lattice Boltzmann Model for Simulation of Magnetohydrodynamics,” Phys. Rev. Lett., 67(27), p. 3776. [CrossRef] [PubMed]
Qian, Y. , dHumieres, D. , and Lallemand, P. , 1992, “ Recovery of Navier–Stokes Equations Using a Lattice-Gas Boltzmann Method,” Europhys. Lett, 17(6), p. 479. [CrossRef]
Chen, H. , Chen, S. , and Matthaeus, W. H. , 1992, “ Recovery of the Navier–Stokes Equations Using a Lattice-Gas Boltzmann Method,” Phys. Rev. A, 45(8), p. R5339. [CrossRef] [PubMed]
Himika, T. A. , Hasan, M. F. , and Molla, M. M. , 2016, “ Lattice Boltzmann Simulation of Airflow and Mixed Convection in a General Ward of Hospital,” J. Comput. Eng., 2016, p. 5405939. [CrossRef]
Zou, Q. , and He, X. , 1997, “ On Pressure and Velocity Boundary Conditions for the Lattice Boltzmann BGK Model,” Phys. Fluids, 9(6), pp. 1591–1598. [CrossRef]
Ghia, U. K. N. G. , Ghia, K. N. , and Shin, C. T. , 1982, “ High-Resolutions for Incompressible Flow Using the Navier–Stokes Equations and a Multigrid Method,” J. Comp. Phys., 48(3), pp. 387–411. [CrossRef]
Erturk, E. , 2008, “ Numerical Solutions of 2-D Steady Incompressible Flow Over a Backward-Facing Step—Part I: High Reynolds Number Solutions,” Comput. Fluids, 37(6), pp. 633–655. [CrossRef]
Armaly, B. F. , Durst, F. , Pereira, J. C. F. , and Schönung, B. , 1983, “ Experimental and Theoretical Investigation of Backward-Facing Step Flow,” J. Fluid Mech., 127, pp. 473–496. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Schematic diagram of hospital ward: (a) without any block or partition, (b) with blocks, and (c) with both blocks and partition

Grahic Jump Location
Fig. 2

Discrete velocity vectors for the D2Q9 model for 2D LBM

Grahic Jump Location
Fig. 3

Comparison with the benchmark results of lid-driven cavity flow by Ghia et al [26] in terms of the (a) U velocity and (b) V velocity while Re = 400

Grahic Jump Location
Fig. 4

Comparison with the benchmark results of flow through a backward facing step by Erturk [27] in terms of the streamlines while Re = 400

Grahic Jump Location
Fig. 5

Present results for the velocity profiles u of a backward facing step at the same location recorded by Armaly et al. [28] and Erturk [27]: (a) x/h=0.0, (b) x/h=2.55, (c) x/h=3.06, (d) x/h=3.57, (e) x/h=4.18, (f) x/h=4.80, (g) x/h=5.41, (h) x/h=6.12, (i) x/h=7.14, (j) x/h=7.76, (k) x/h=8.52, (l) x/h=9.18, (m) x/h=9.74, (n) x/h=11.07, (o) x/h=11.84, and (p) x/h=13.57 while Re = 389

Grahic Jump Location
Fig. 6

Velocity profiles u at the different location of a backward facing step by Armaly et al. [28] (top) and Erturk [27] (bottom) while Re = 389

Grahic Jump Location
Fig. 7

Streamlines appended on u velocity contour for (a) Re = 100, (b) Re = 250, and (c) Re = 350

Grahic Jump Location
Fig. 8

Streamlines appended on u velocity contour for (a) Re = 100, (b) Re = 250, and (c) Re = 350

Grahic Jump Location
Fig. 9

Isotherms for (a) Re = 100, (b) Re = 250, and (c) Re = 350

Grahic Jump Location
Fig. 10

u velocity profiles (without partition case) at different position of x: (a) x = 0.5, (b) x = 1, (c) x = 1.5, (d) x = 2, (e) x = 2.5, (f) x = 3, (g) x = 3.5, (h) x = 4, (i) x = 5, (j) x = 5.5, (k) x = 6, (l) x = 6.5, (m) x = 7, and (n) x = 7.5

Grahic Jump Location
Fig. 11

Velocity distribution (without partition case) at different Re: (a) u and (b) v at mid y

Grahic Jump Location
Fig. 12

Streamlines appended on u velocity contour for (a) Re = 100, (b) Re = 250, and (c) Re = 350

Grahic Jump Location
Fig. 13

Isotherms for (a) Re = 100, (b) Re = 250, and (c) Re = 350

Grahic Jump Location
Fig. 14

u velocity profiles (with partition case) at different position of x: (a) x = 0.5, (b) x = 1, (c) x = 1.5, (d) x = 2, (e) x = 2.5, (f) x = 3, (g) x = 3.5, (h) x = 4, (i) x = 5, (j) x = 5.5, (k) x = 6, (l) x = 6.5, (m) x = 7, and (n) x = 7.5

Grahic Jump Location
Fig. 15

Velocity distribution (with partition case) for the different Re: (a) u velocity and (b) vvelocity at mid y

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In