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

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