An empirical drag coefficient model for simulating the ...

13 downloads 0 Views 5MB Size Report
Building Services Engineering Research and Technology. ...... particles instead of outdoor field measurement in the previous aerobiology studies for a small ...
An empirical drag coefficient model for simulating the dispersion and deposition of bioaerosol particles in ventilated environments

YU HO CHING

Ph.D

The Hong Kong Polytechnic University

2016

1

The Hong Kong Polytechnic University Department of Building Services Engineering

An empirical drag coefficient model for simulating the dispersion and deposition of bioaerosol particles in ventilated environments

Yu Ho Ching

A thesis submitted in partial fulfilment of the requirements

for the degree of Doctor of Philosophy

June 2016

2

Certificate of Originality

I hereby declare that this thesis is my own work and that, to the best of my knowledge and belief, it reproduces no material previously published or written, nor material which has been accepted for the award of any other degree or diploma, except where due acknowledgement has been made in the text.

_______________________________ Yu Ho Ching

Department of Building Services Engineering The Hong Kong Polytechnic University Hong Kong, China June 2016 i

Abstract

Bioaerosol particles in indoor air are related to airborne transmission infections and some pandemic outbreaks such as Severe Acute Respiratory Syndrome (SARS) in 2003 and Middle East Respiratory Syndrome (MERS) in 2015. Several environmental control strategies and parameters for a ventilation system have been suggested to prevent infections in building environments. To design an appropriate ventilation system, the infection risks of the proposed ventilation system were evaluated in the thesis in order to achieve effective infection control. Computational fluid dynamics (CFD) simulation is often used to predict the dispersions and depositions of bioaerosol particles to evaluate the infection risks of ventilation systems. However, there are differences between bioaerosol and aerosol particles in terms of shape, diameter, surface texture and elasticity. In this study, the transport mechanism of a bioaerosol particle was investigated to formulate a bioaerosol particle transport model for CFD simulation.

The empirical bioaerosol drag coefficient model was developed in this study to investigate the transport mechanism of bioaerosol particles. A chamber study was used to collate the empirical data from 13 common indoor bioaerosol species with the three common ventilation rates (1.7, 10.3 and 18.8 ACH). By comparing the experimental and numerical data, the empirical drag constants and coefficients were determined for each bioaerosol species.

ii

The model (i.e. Kdrag,bp=dbp2/2) was developed by correlating the drag constants Kdrag,bp with the equivalent bioaerosol diameters dbp in a range between 0.054 and 6.9 µm. Several validations were done for the generalization of the model for various bioaerosol species, ventilation rates, enclosures and literature. The model simplifies the transport mechanism of bioaerosol particles, for example, dispersion and deposition, in terms of the equivalent bioaerosol diameter dbp and drag coefficient Cdrag,bp. This is beneficial in that only a single morphological characteristic (i.e. dbp) is required to predict the movement of any bioaerosol species.

A numerical bioaerosol transport framework has been extended based on the proposed model to simulate the bioaerosol distribution to enhance the applicability of the model and impact on the ventilation system design for infection control in terms of ventilation rate and other design factors. The impacts of the proposed model and framework were demonstrated by simulating three practical scenarios such as healthcare centre, sanitation and office. The over-predictions of the drag force and ventilation performance by the Stokes drag was recognized, especially in environments with a unidirectional airflow pattern. The ventilation strategies for infection control need to be reviewed urgently because of the over-prediction of the carrying power of the airflow by the Stokes drag coefficient model.

In this study, the correlation between the drag constant Kdrag,bp and the equivalent bioaerosol diameter dbp has been investigated. This study provides iii

a useful source of reference for ventilation system engineers to minimize the infection risk of airborne transmission diseases, and to mitigate the risk of outbreaks. However, some improvements are suggested to enhance the reliability of the model. Furthermore, the development of the atomistic drag model (i.e. kinetic theory) may provide a solid theoretical base to support the model.

iv

Publications SCI and EI Journal papers *

Yu. H. C., Mui, K. W. Wong, L. T. (2017). Numerical simulation of bioaerosol

particle exposure assessment in office environment from MVAC systems, The Journal of Computational Multiphase Flows. [Submitted] *

Yu, H. C., Mui, K. W., Wong, L. T., & Chu, H. S. (2016). Ventilation of general

hospital wards for mitigating infection risks of three kinds of viruses including Middle East Respiratory Syndrome Coronavirus (MERS-CoV). Indoor and Built Environment. doi: 10.1177/1420326X16631596 *

Wong, L. T., Yu, H. C., Mui, K. W., & Chan, W. Y. (2015). Drag constants of

common indoor bioaerosols. Indoor and Built Environment, 24(3), 401-413. doi: 10.1177/1420326X13515897 *

Mui, K. W., Wong, L. T., Yu, H. C., Cheung, C. T., & Li, N. (2016). Exhaust

ventilation performance in residential washroom for bioaerosol particle removal after water closet flushing. Building Services Engineering Research and Technology. doi: 10.1177/0143624416660597

Mui, K. W., Wong, L. T., Cheung, C. T. & Yu, H. C. (2016). Cooling energy for indoor environmental quality (IEQ) acceptance in demand-controlled ventilated and adaptive comfort temperature-controlled air-conditioned offices, HKIE transactions, 24(2), 78-87. doi: 10.1080/1023697X.2017.1312561 *

Lai, A. C. K., Wong, L. T., Mui, K. W., Chan, W. Y., & Yu, H. C. (2012). An

experimental study of bioaerosol (1-10µm) deposition in a ventilated chamber. Building and Environment, 56, 118-126. doi: 10.1016/j.buildenv.2012.02.027.

*

Publication related to the thesis v

International and local conference papers *

Yu, H.C., Mui, K.W., Wong, L.T., & Chan, W.Y. Numerical study of 10μm

bioaerosols (Rhizopus) deposition in a forced-ventilated chamber. Paper presented at the Healthy buildings 2012, The 10th International Conference, 8-12 Jul, Brisbane, Queensland. *

Mui, K. W., Yu, H. C., & Wong, L. T. (2014b). Validation of the bioaerosol

deposition model in ventilated chamber (Paper HP0585). Paper presented at the 13th International Conference on Indoor Air Quality and Climate, Indoor Air 2014, 7-12 July, Hong Kong. *

Mui, K. W., Wong, L. T., & Yu, H. C. (2014a). Determine the aerodynamic

properties of Legionella pneumophila for a drag force expression Paper presented at the 40th International Symposium on Water Supply and Drainage for Buildings, CIBW062 Symposium 2014, 8-10 September, São Paulo. Brazil.

Cheung, C. T., Mui, K. W., Wong, L. T. & Yu, H. C. (2014). Mobile application of indoor environmental quality (IEQ) calculator in air-conditioned offices and university classrooms. Joint Symposium 2014, Change in Building Services Engineering for Future, 25-November: 2-1 to 2-5.

Mui, K. W., Wong, L. T., Xiao, F., Cheung, C. T., & Yu, H. C. (2015). Use of sustainable building environmental model (SBEM) in Hong Kong air-conditioned buildings. Paper presented at the 13th Asia Pacific Conference on the Built Environment, 19-20 November, Hong Kong, China, pp. 555-563

Mui, K. W., Wong, L. T., Xiao, F., Cheung, C. T., & Yu, H. C. (2015). Development of sustainable building environmental model (SBEM) in Hong Kong. Paper presented at the ICEEE 2015: 17th International Conference on Energy and Environmental Engineering, London, United Kingdom.

*

Publication related to the thesis vi

Acknowledgements

I am heartily thankful to my chief supervisor Dr. Mui Kwok-wai, and my cosupervisor, Dr. Wong Ling-tim, for their valuable guidance and suggestions throughout my research work, as well as their indispensable support and encouragement for my personal development from the initial to the final stage of my study.

Thanks are also extended to whom have assisted with this research study, and my friends who have been supportive both mentally and technically, which helped me to pass through all hurdles.

Lastly, I would like to express my most sincere appreciation to my family especially to, my mother. This thesis would not be successfully completed with their unconditional support and endless love.

vii

Table of Content Certificate of Originality Abstract Publications Acknowledgements Table of Content List of Figures List of Tables List of Abbreviations List of Symbols

i ii v vii viii xii xix xxi xxiv

Chapter 1 Introduction 1.1 Background of bioaerosol particle simulation in buildings 1.1.1 Impacts of Infectious disease outbreaks 1.1.2 Spread of airborne infectious diseases 1.1.3 Infection in indoor environments 1.1.4 Infection risk assessment for ventilation system design 1.1.5 Prediction of bioaerosol movement for infection risk 1.2 Research Objectives 1.3 Research scope 1.4 Organization of the thesis Chapter 2 Literature review of Bioaerosol transport model 2.1 Bioaerosol transport models for ventilaton system 2.2 Bioaerosols 2.3 Bioareosol infection risk models for ventilation system 2.3.1 Epidemiological model for infection risk 2.3.2 Wells-Riley quantum model for infection risk 2.3.3 Dose-Response model for infection risk 2.4 Prediction of bioaerosol particle movement 2.5 CFD simulation for bioaerosol particle transport models 2.5.1 Drift-flux model (DFM) for bioaerosol transport 2.5.2 Discrete phase model (DPM) for bioaerosol transport 2.6 Drag coefficient for transport model 2.6.1 Analytical model of drag coefficient of a rigid sphere 2.6.2 Morphology and rheological factors of drag coefficients 2.6.3 Gas-kinetic theory for drag coefficient in rarefied gas 2.7 The uncertainties of drag coefficient for bioaerosol viii

1 3 3 5 8 21 21 22 23 29 33 33 36 39 39 43 48 55 59 60 63 69 70 76 95 97

Chapter 3 Development of bioaerosol drag coefficient model 3.1 Introduction of the research methodology 3.1.1 Empirical approach for bioaerosol drag coefficient model 3.1.2 Chamber study for a bioaerosol particle trajectory 3.2 Experimental study for Chamber A 3.2.1 Experimental setup and procedures in Chamber A 3.3 Numerical study for Chamber A 3.3.1 Eulerian framework for airflow field simulation in Chamber A 3.3.2 Lagrangian framework for bioaerosol particles simulation in Chamber A 3.3.3 User defined function (UDF) for the empricial bioaerosol drag constant simulation 3.4 Summary of the development of the chamber study Chapter 4 Empirical bioaerosol drag coefficient model 4.1 Introduction of the empirical drag coefficient model 4.2 Determination of empirical drag constants from experiment and simulation results of Chamber A 4.2.1 Bioaerosol fractional count of the deposition patterns 4.2.2 Experimental fractional counts 4.2.3 Simulated fractional count of numerical Chamber A 4.2.4 Empiricald drag constants from Chamber A study 4.3 Bioaerosol drag coefficient model 4.3.1 Correlation of the bioaerosol drag coefficient and morphological characteristic of the bioaerosol 4.3.2 Error analysis for the bioaerosol drag coefficient model 4.4 Implementation of the empirical bioaerosol drag coefficient model 4.4.1 Drag coefficient and Reynolds number of the bioaerosol 4.4.2 Comparison of the bioaerosol drag coefficient model and Stokes-Cunningham drag coefficient model 4.5 Validations for the bioaerosol drag coefficient model 4.5.1 Validation for ventilation rate (Validation A) 4.5.2 Validation for bioaerosol species (Validation B) 4.5.3 Validation for viral group (Validation C) 4.5.4 Validation for Chamber B (Validation D) 4.5.5 Validation summary 4.6 Summary of the empirical bioaerosol drag coefficient model

ix

98 98 100 106 109 109 123 124 129 133 136 138 138 140 140 141 147 154 156 156 157 165 165 168 170 171 173 175 178 180 184

Chapter 5

185

5.1

185

5.1.1 5.1.2 5.1.3 5.2 5.2.1 5.2.2 5.2.3 5.2.4 5.2.5 5.3 5.3.1 5.3.2 5.3.3 5.3.4 5.3.5 5.3.6 5.3.7 5.4 5.4.1 5.4.2 5.4.3 5.4.4

Development of numerical bioaerosol transport framework with the empirical bioaerosol drag coefficient model Introduction of the numerical bioareosol transport framework Development of the numerical bioaerosol transport framework Validations for the numerical transport framework by practical environment from the open literature The applicability of the proposed numerical bioaerosol transport framework Transport characteristics in the hospital ward for human expiratory activites Background of a hospital ward CFD simulation of a general ward (i.e. Ward A) Bioaerosol removal performance of ventilation systems Results and discussion of hospital ward simulations Impact of the bioaerosol drag coefficient model in Ward A simulation Exhaust ventilation performance of residential washroom for bioaerosol removal after water closet (WC) flush Introduction of a residential washroom application (i.e. Washroom A) Exhaust ventilation design for washrooms Scenario study for design factors Modification of the framework for water closet flushing emission Results and discussion of the exhaust ventilation performance in residential washrooms Capacity and design procedure for local exhaust system Impact of the bioaerosol drag coefficient model in Washroom A simulation Potential continuous bioaerosol emission in office environment from MVAC systems Introduction of an office application (i.e. Office A)Introduction of an office application Development of the numerical bioaerosol transport framework for continuous emission source Bioaerosol exposure level of an office cubicle of Hong Kong Impact of the bioaerosol drag coefficient model in Office A simulation x

188 189 194 196 198 200 206 208 218 220 220 224 226 229 233 241 243 246 246 248 252 263

Summary of the numerical bioaerosol transport framework with the bioaerosol drag coefficient model Chapter 6 Conclusions and recommendations 6.1 Findings to address the research objectives and question 6.2 Significance of the study 6.3 Limitations, recommendations and future works Appendix A Matlab script program for the calculation of simulated fractional count Appendix B Fluent UDF program for simulating the drag constant Appendix C Fluent Scheme program for batch automation Appendix D Fluent UDF program for the numerical transport framework Appendix E Fluent Scheme program for the continuous emission source Reference 5.5

xi

266 268 270 274 275

278 281 282 285 286 288

List of Figures

Figure 1.1

The structure diagram of Chapter 1

Figure 1.2

Airborne infection environments

Figure 1.3

Bioaerosol systems

Figure 1.4

The outline of the thesis

32

Figure 2.1

Structure diagram of bioaerosol transport model for infection risk models

35

Figure 2.2

Size range of bioaerosol (in µm)

36

Figure 2.3

Illustration of bioaerosol particles in droplet nucleus and droplet forms

37

Figure 2.4

The Wells evaporation-falling airborne and droplets

38

Figure 2.5

The cycle of infection transmission model

40

Figure 2.6

Infection cycle of SIR models

41

Figure 2.7

An example of outbreak for airborne infection prediction using SIR model

42

Figure 2.8

Infection cycle of SEIR model

43

Figure 2.9

Infection probability of the intake dose and ID50

49

Figure 2.10

Exposure and dose-response assessments for risk infection framework

52

Figure 2.11

Creeping flow over a sphere

71

Figure 2.12

Drag coefficient with various Reynolds number Rep

77

factors

transmission

xii

with

2

in

indoor

9

ventilation

18

curve

of

Figure 2.13

Drag coefficient with sphericity ratios

81

Figure 2.14

Drag coefficients for bubbles

84

Figure 2.15

Structure of a bacterial cell

87

Figure 2.16

Surface roughness on drag coefficient

88

Figure 2.17

Mechanism of flagella movement. (a) Relating the local viscous force to the local filament velocity relative to the fluid, (b) a planar sinusoidal flagellar wave and (c) a helical flagellar wave

90

Figure 2.18

Drag coefficient with Knudsen numbers Kn

92

Figure 2.19

Regimes of gas flow and bioaerosol species range

94

Figure 3.1

Structure diagram for the investigation of an empirical drag coefficient model for bioaerosols

99

Figure 3.2

The morphological and rheological characteristics for bioaerosol drag constants

101

Figure 3.3

Reference SEM photos of 13 common indoor bioaerosol species

103

Figure 3.4

Determine the projected image area, length and width by ImageJ

105

Figure 3.5

Trajectory of a bioaerosol particle under ventilation rates in the field of drag and gravity forces

106

Figure 3.6

The chamber study to determine the empirical drag coefficients and constants for the 13 bioaerosol species

108

Figure 3.7

Procedures of the experimental chamber for bacterial group

110

Figure 3.8

Double-vial package for microorganism strain. (a) photo of the package, and (b) illustration of the package

111

Figure 3.9

Streaked dish for the isolation of single colony of the microorganism

112

xiii

Figure 3.10

Setup of Chamber A

114

Figure 3.11

Six-jet collision nebulizer for the experiment. (a) Empty nebulizer, and (b) Nebulizer with the Ringer’solution

116

Figure 3.12

Generation of the bioaerosol by nebulizer

116

Figure 3.13

Formation of bioaerosol droplet nuclei by diffusion dryer

117

Figure 3.14

Photo of the aerosolization box for bioaerosols with equivalent bioareosol diameter dbp> 3m

118

Figure 3.15

Photo of the matured aerosolization box

the

119

Figure 3.16

Colonies of P. citrinum in the chamber with VRchamber of 10.3 ACH

121

Figure 3.17

Eulerian-Lagrangian framework for numerical Chamber A

124

Figure 3.18

Geometry of the numerical Chamber A

126

Figure 3.19

Flowchart of bioaerosol particle procedures in DPM

130

Figure 3.20

Flowchart of the automation process for CFD simulation

135

Figure 4.1

Structure diagram of the empirical bioaerosol drag coefficient

139

Figure 4.2

Experimental fractional counts on the chamber floor

146

Figure 4.3

Air velocity distribution in Chamber A. a) VRchamber=1.7 ACH, (b) VRchamber=10.3 ACH and (c) VRchamber=18.8 ACH

148

Figure 4.4

Simulated bioaerosol fractional counts by Stokes drag with experimental fractional counts

153

Figure 4.5

Bioaerosol drag constants Kdrag,bp and absolute errors for the 13 tested bioaerosol species

154

Figure 4.6

Bioaerosol drag constant Kdrag,bp against equivalent bioaerosol diameter dbp

156

xiv

spores

for

Figure 4.7

Simulated bioaerosol fractional counts by bioaerosol drag coefficient model with experimental fractional counts

163

Figure 4.8

Errors between simulated and measured bioaerosol deposition patterns. (a) єfb and (b) єnmse

164

Figure 4.9

Bioaerosol drag constant against the Reynolds number of the equivalent bioaerosol particles

167

Figure 4.10

Bioaerosol drag constants against Cunningham slip effect

169

Figure 4.11

Bioaerosol fractional count for Validation A

171

Figure 4.12

Absolutes errors with other ventilation rate settings

172

Figure 4.13

Referenced micrograph of bioaerosol for Validation B

173

Figure 4.14

Bioaerosol drag constants and absolute errors for Validation B

174

Figure 4.15

Electron micrograph of Escherichia coli Phage Phi x 174 (ATCC 13706-B1)

176

Figure 4.16

Absolute errors of Escherichia coli Phage Phi x174

177

Figure 4.17

Chamber B and agar dishes setup.

178

Figure 4.18

Experiment setup for Chamber B

179

Figure 4.19

Absolute errors of E. coli (ATCC 13706) in Chamber B

180

Figure 4.20

Bioaerosol drag constant Kdrag,bp against equivalent bioaerosol diameter dbp for validations

181

Figure 4.21

Reynolds number and drag coefficient for: (a) Validation A, (b) Validation B, (c) Validation C and, (d) Validation D

182

Figure 5.1

Structure diagram of the development of numerical bioaerosol transport framework for ventilation system design

187

xv

Figure 5.2

Framework framework

for

the

Figure 5.3

The washroom of Validation E

191

Figure 5.4

The washroom of Validation F

193

Figure 5.5

Structure diagram of bioaerosol simulation in a hospital general ward

197

Figure 5.6

A 6-bedded general ward cubicle (dimensions in mm)

200

Figure 5.7

CFD configurations of a 6-bedded general ward cubicle

201

Figure 5.8

Electron micrograph of reference viruses. (a) MERS-CoV (b) SARS-CoV and (c) H1N1 influenza virus

202

Figure 5.9

Time-varying fractional counts of bioaerosol particle

207

Figure 5.10

Simulation results of the cubicle with VRward = 6 ACH. (a) Temperature distribution, (b) Air velocity distribution, (c) Flow pathlines from air supply inlets

208

Figure 5.11

MERS-CoV pathways for 6 source locations with VRward = 6 ACH. (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 4, (e) Man 5 and (f) Man 6

210

Figure 5.12

MERS-CoV pathways for 6 source locations with VRward=6 ACH. (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 4, (e) Man 5 and (f) Man 6

212

Figure 5.13

MERS-CoV removal processes for different source locations and VRward

214

Figure 5.14

SARS-CoV removal processes for different source locations and VRward

215

Figure 5.15

H1N1 influenza virus removal processes for different source locations and VRward

216

Figure 5.16

Average max. elapsed time tmax,air with design standards

217

xvi

proposed

transport

188

Figure 5.17

Average max. elapsed time tmax,air between the bioaerosol drag and the Stokes drag coefficient models

219

Figure 5.18

Structure diagram for an effective exhaust system design for washroom

223

Figure 5.19

Illustration washroom

inside

225

Figure 5.20

Overall design factors for an effective exhaust system

226

Figure 5.21

Ventilation arrangement of Washroom A

228

Figure 5.22

Modified framework for the Washroom A simulation

230

Figure 5.23

Exhausted fractional count FCexh against ventilation rate VRwashroom in Washroom A

235

Figure 5.24

Max. elapsed time tmax,air against VRwashroom in Washroom A

236

Figure 5.25

Bioaerosol particle fractional counts FCdepos, FCexh and FCair against time t.

237

Figure 5.26

Simulated air velocity distribution

238

Figure 5.27

Bioaerosol exhausted fractional count FCexh against distance from WC to exhaust fan lfan

239

Figure 5.28

Bioaerosol removal rate in Washroom A

240

Figure 5.29

Flow diagram to improve washroom exhaust system

242

Figure 5.30

Difference between the bioaerosol and the Stokes drag in the local exhaust ventilation system for Washroom A

243

Figure 5.31

Difference of bioaerosol removal rate between the bioaerosol and the Stokes drag coefficient models

245

Figure 5.32

Structure diagram for a continuous emission source in an office cubicle

247

Figure 5.33

Typical mixing ventilation for an office

248

of

design

xvii

variables

Figure 5.34

Description of the bioaerosol removal process.

251

Figure 5.35

Floorplan of Office A (in m)

253

Figure 5.36

CFD layout of Office A.

254

Figure 5.37

S. aureus concentration for ventilation rates and emission concentrations

257

Figure 5.38

C. cladosporides concentration for ventilation rates and emission concentrations

258

Figure 5.39

Fractional count of S. aureus for ventilation rates and emission concentrations

259

Figure 5.40

Fractional count of C. cladosporides concentration for ventilation rates and emission concentrations

260

Figure 5.41

Fractional count of bioaerosol particles for single-shot emission at 1 ACH of VRoffice. (a) S. aureus and (b) C. cladosporides

261

Figure 5.42

Bioaerosol exposure level with different emission concentrations. (a) S. aureus and (b) C. cladosporides

262

Figure 5.43

Average exposure level of the empirical bioaerosol and the Stokes drag in the office cubicle simulation

264

Figure 5.44

Flow pathlines in the office cubicle simulation

265

Fgure 6.1

Structure diagram for the investigation in the thesis

269

Figure 6.2

Roadmap of future works

276

xviii

List of Tables

Table 1.1

Some airborne pathogens for communicable, non-communicable opportunistic disease groups

the and

6

Table 1.2

Selected measurement methods of droplet size distribution of expiratory activities

12

Table 1.3

Summary of selected experiments of toilet seeding

14

Table 1.4

Selected literature on building type related to bioaerosol infections

17

Table 1.5

Selected literature on ventilation systems related to bioaerosol infections

19

Table 2.1

Microorganism morphologies

68

Table 2.2

Stokes correction factors for Reynolds numbers in flow regions

78

Table 2.3

Stokes correction factors for spheroid shaped particles

80

Table 2.4

Aerodynamic diameters of some bioaerosol species

83

Table 2.5

Cunningham slip correction parameters from selected experiments

93

Table 3.1

Common indoor bioaerosol species

102

Table 4.1

Bioaerosol drag constant for the 13 common indoor bioaerosol

155

Table 4.2

Particle tracking results by the bioaerosol drag coefficient model

158

Table 4.3

Reynolds numbers of bioaerosol particle and relative velocities of the 13 bioaerosol species

166

xix

Table 4.4

Validation cases for the empirical bioaerosol drag coefficient model

170

Table 4.5

Validation A results for VRchamber

172

Table 4.6

Validation B results for other bioaerosol species

175

Table 4.7

Validation C results for Escherichia coli Phage Phi x 174

178

Table 5.1

Validations E and F results for practical environments

192

Table 5.2

Application cases for the numerical simulation framework

195

Table 5.3

Virus information of general ward simulation

203

Table 5.4

CFD simulation and boundary condition settings for the general inpatient ward cubicle

205

Table 5.5

Selected exhaust fan locations for simulations in Washrooms A

227

Table 5.6

Information of the emitted bioaerosols in flushing

229

Table 5.7

Configurations and boundary settings for Washroom A simulation

233

Table 5.8

Information of the bioaerosol species

252

Table 5.9

CFD simulation conditions.

settings

xx

and

boundary

255

List of Abbreviations

ACGIH

American Conference of Governmental Industrial Hygienists

ACH

Air change rate per hour

AHAM

Association of Home Appliance Manufacturers

AHU

Air handling unit

AIIR

Airborne infection isolation room

ASHRAE

American Society of Heating, Refrigerating and Air-Conditioning Engineers

ASTM

American Society for Testing and Materials

ATCC

American Type Culture Collection

ATD

Air terminal device

BBO

Basset-Boussinesq-Oseen

CAV

Constant air volume

CDC

Centers for Disease Control and Prevention

CDF

Cumulative distribution function

CFD

Computational fluid dynamics

CFU

Colony forming unit

CIBSE

Chartered Institution of Building Services Engineers

CO

Carbon monoxide

CO2

Carbon dioxide

COPD

Chronic obstructive pulmonary disease

CoV

Coronavirus

DES

Detached eddy simulation

DFM

Drift-flux model

DNA

Deoxyribonucleic acid

DPM

Discrete phase model

DRW

Discrete random walk

ELISA

Enzyme-linked immunosorbent assay

FB

Fractional bias

FCU

Fan-coil unit

FVM

Finite volume method

GCI

Grid convergence index xxi

HAI

Hospital acquired infection

HCW

Healthcare worker

HEPA

High efficiency particulate air

HFMD

Hand, foot and mouth disease

HKHA

Hong Kong Hospital Authority

HVAC

Heating, ventilating, and air conditioning

IAQ

Indoor air quality

ICU

Intensive care unit

ID50

Infectious dose

IMI

Interferometric Mie imaging

LES

Large eddy simulation

MEA

Malt extract agar

MEB

Malt extract broth

MERS

Middle East Respiratory Syndrome

MERV

Minimum efficiency reporting value

MRSA

Methicillin-resistant Staphylococcus aureus

MV

Mixing ventilation

MVAC

Mechanical ventilation and air conditioning

MVOC

Microbial volatile organic compounds

NCRP

National Council on Radiation Protection, and Measurements

NIH

National Institute of Health

NIOSH

National Institute for Occupational Safety and Health

NIST

National Institute of Standards and Technology

NMSE

Normalized mean standard error

OPC

Optical particle count

PAU

Primary air unit

PCA

Plate count agar

PCR

Polymerase chain reaction

PDE

Partial differential equation

PFU

Plaque forming unit

PISO

Pressure implicit with split operator

PIV

Particle image velocimetry

PV

Personalized ventilation

qPCR

Quantitative polymerase chain reaction xxii

RA

Relative abundances

RANS

Reynolds-averaged Navier-Stokes

RH

Relative Humidity (%)

RNA

Ribonucleic acid

RNG

Renormalization group

RT-PCR

Reverse transcriptase polymerase chain reaction

SARS

Severe Acute Respiratory Syndrome

SD

Standard deviation

SEIR

Susceptible-exposed-infected-recovered

SEIS

Susceptible-exposed-infected-susceptible

SEM

Scanning electron microscopy

SIMLPE

Semi-implicit method for pressure linked equation

SIS

Susceptible-infected-susceptible

SIR

Susceptible-infected-recovered

SIRS

Susceptible-infected-recovered-susceptible

TB

Tuberculosis

TCID

Tissue culture infectious dose

TEM

Transmission electron microscopy

TLV

Threshold level value

TSA

Tryptone soya agar

TSB

Tryptone soya broth

T-RFLP

Terminal restriction fragment length polymorphism

UA

Ultrasonic anemometry

UDF

User defined function

UV

Ultraviolet

UVGI

Ultraviolet germicidal irradiation

VAV

Variable air volume

VBNC

Viable but non-culturable

VR

Ventilation rate

WC

Water closet

WHO

World Health Organization

xxiii

List of Symbols

A

Area (µm2)

Aproj

Projected image area of a particle (µm2)

Bo

Bond number

C

Coefficient

Cdrag

Drag coefficient

Ccorr

Coefficient of correlation

Cdetm

Coefficient of determination

Conc

Concentration (kg m-3, CFU m-3)

d

Diameter (µm)

daero

Aerodynamic diameter (µm)

dinit

Microbe-laden droplet diameter (µm)

dbp

Equivalent bioaerosol diameter (µm)

Expo

Exposure level (CFU min L-1)

F

Force (N)

Fdrag

Drag force (N s kg-1)

Fadd

Additional acceleration force per unit particle mass (N kg-1)

FC

Fractional count

FCair

Fractional count of bioaerosol particles suspended in the air

FCexh

Fractional count of bioaerosol particles removed through exhaust

FCdepos

Fractional count of bioaerosol particles deposited onto surfaces

Fr

Froude number

f

Fraction

fasymp

Asymptotic range of convergence

favg

Average volume fraction of room air

fbp

Volumetric fraction of bioaerosol particles (mL) over air (L)

fmolecule

Number of molecules per unit volume (m-3)

fp

Fraction of particle density to continuous-fluid density

freflection

Reflection fraction of the molecules from the surface of the rigid sphere

fsafe

Safety factor

fslip

Cunningham slip correction factor xxiv

GCI

Grid convergence index

g

Gravitational acceleration (m s-2)

ID50

Human infectious dose (ID50)

KBoltzmann

Boltzmann constant (J K-1)

Kconv

Theoretical order of convergence

Kdrag

Drag constant

Kgas

Ideal gas constant (i.e. 8.314 J mol-1 K-1)

KSL

Lift force constant

Kn

Knudsen number

Keff,therm_condy

Effective values of thermal conductivity (W m-1 K-1)

lfan

Distance from WC to exhaust (m)

lmean

Molecular mean free path (µm)

lp

Radius of a solid sphere (µm)

l1

Length of a bioaerosol particle (µm)

l2

Width of a bioaerosol particle (µm)

Mach

Mach number

Mo

Morton number

m

Mass (kg, g)

N

Number

Nair

Number of bioaerosol particles suspended in the air

Nexh

Number of bioaerosol particles removed through exhaust

Ndepos

Number of bioaerosol particles deposited onto surfaces

Ninit,quanta

Number of initial quanta value (quanta)

Nsrc

Number of bioaerosol particles emitted from source

Nwc

Number of microorganisms in the WC seal

Prob

Probability

Prz

Pressure (pa)

p

p-value

Q

Flow rate (m3 s-1, L s-1, m3 min-1)

R

Rate (s-1)

Rcontact

Contact rate (person day-1)

Rdepos

Rate deposition loss of the bioaerosol particles

Routbreak

Basic reproduction number

Rquanta

Quanta generation rate (quanta min-1) xxv

Rrecover

Recovery rate (person day-1)

Rdepos

Rate deposition loss of the bioaerosol particles (min-1)

Reair

Reynolds number for air

Rebp

Reynolds number for bioaerosol particles

r

Ratio

raspect

Aspect ratio of bioaerosol particles

rrefine

Refinement ratio

Schturbulent

Turbulent Schmidt number

Src

Source term

T

Temperature (°C, K)

TCID50

Tissue culture infectious dose (TCID)

t

Time (s, min, hr)

tinit

Initial time (s)

texpo

Exposure time interval (s)

tmax,air

Maximum elapsed time (s)

V

Volume (m3)

Vroom

Room volume (m3)

VR

Ventilation rate (ACH)

VRchamber

Ventilatoin rate of a chamber (ACH)

VRoffice

Ventilatoin rate of an office (ACH)

VRward

Ventilatoin rate of a ward (ACH)

VRwashroom

Ventilatoin rate of a washroom (ACH)

v

Velocity (m s-1)

vair

Air velocity (m s-1)

vbp

Bioaerosol particle velocity (m s-1)

vsettle

Settling velocity (m s-1)

vrel

Relative velocity (m s-1)

vwc

Air velocity at a height of 0.2 m above the WC seal (m s-1)

viab

Viability function

Web

Weber number

є

Error

єrms

Relative error of computed average mass flow rate

єfb

Fractional bias xxvi

єnmse

Normalized mean standard error

θ

Inclined angle (rad)

ρ

Density (kg m-3)

ρair

Air density (kg m-3)

ρbp

Density of the bioaerosol particles emitted (kg m-3)

σ

Diffusivity (m2 s-1)

σbin

Dispersed phase diffusivity (m2 s-1)

ξ

Surface tension (N m-1)

Γ

Surface concentration of surfactant (mol m-2)

τ

Shear stress (N m-2),

μ

Molecular viscosity(kg m-1 s-1)

μair

Molecular viscosity of air (kg m-1 s-1)

μkin

Kinematic viscosity (kg m-1 s-1)

μeff

Effective viscosity of the air (kg m-1 s-1)

Ψ

deformation rate tensors

𝜓

Stream function

ω

Vorticity



Difference



Gradient of a scalar function

Subscripts 0, 1, 2…

of conditions 0, 1, 2…

AM

of added mass

add

of additional

aero

of aerodynamic

air

of air

alveolar

of alveolar

aspect

of aspect

avg

of average

axis

of axis

BH

of Basset history

Boltzmann

of Boltzmann xxvii

Brown

of Brownian

bin

of bin number

bp

of bioaerosol particle

Cunningham

of Cunningham

chamber

of chamber

coarse

of coarse grid

crit

of critical

data-pair

of data pair

depos

of deposition

dose

of dose

drag

of drag

Epstein

of Epstein

eff

of effective

emass

of equivalent mass

emission

of emission

emob

of mobility diameter

epa

of equivalent projected area

esurf

of equivalent surface area

exh

of exhausted

exp

of experiment

expo

of exposure

evol

of equivalent volume

FK

of Froude-Krylov

Feret

of Feret

fb

of fractional bias

fine

of fine gird

floor

of floor

grav

of gravity

infect

of infection

inhale

of inhale

init

of initial

intake

of intake

interface

of interface

Martin

of Martin xxviii

max

of maximum

mean

of mean free path

min

of minimum

molecule

of molecule

mvac

of mechanical ventilation and air conditioning (MVAC)

new

of new infection

nmse

of normalized mean standard error

Oseen

of an Oseen approximation

ocp

of occupant

office

of office

outbreak

of outbreak

p

of particle

pop

of population

prz

of pressure

proj

of projected

pulmonary

of pulmonary

quanta

of quanta

random

of random

recover

of recover

refine

of refinement

reflection

of reflection

rel

of relative

relax

of relaxation

removal

of removal

resp

of respiratory

respirator

of respirator

room

of room

SL

of Saffman’s lift

Stokes

of Stokes

scale

of scale

settle

of settling

shape

of shape factor

sim

of simulation

site

of site xxix

slip

of Cunningham slip correction

sph

of sphericity

sphere

of sphere

src

of source

stress

of stress

surface

of surface

suscept

of susceptible

term

of terminal

turbulent

of turbulent

vaero

of vacuum aerodynamic

vol

of volume

wall

of wall

ward

of ward

washroom

of washroom

wc

of water closet

wf

of wall friction

χ

of dynamic shape factor



of infinite

xxx

1

Introduction

In this chapter, the background and rationale are presented to highlight the significance of this thesis. The goal is to develop an empirical drag coefficient model and transport framework for bioaerosol (i.e. biological aerosol) in order to understand the fundamental transport mechanism of bioaerosol particles, to determine effective strategies of infection control and to prevent airborne infectious diseases as illustrated in Figure 1.1. The relation of the public health and ventilation system design is presented regarding the outbreaks of infectious diseases, infection control and prevention and the indoor environmental conditions in Section 1.1. The spatial and temporal predictions of the bioaerosol distribution are recognized to be critical factors for infection control by ventilation systems. The scope and objectives of the thesis are identified with respect to the findings in this study. The outline of the thesis is summarized in Section 1.4.

1

Public health

Outbreaks of infectious disease

Communi -cable disease

Noncommunicable disease

Opportunistic disease

Ventilation system design

Bioaerosol infection risk

Ventilation system

Source emission

Human expiratory activities

Sanitary system operations

Ventilation system operations

Building environments

Distribution type

Ventilation rate

Healthcases

Transport -ations

Temperature

Relative humidity

Sanitary facilities

Offices

Airflow pattern

Recirculation ratio

Disinfection

Filtration

factories

Residences

Figure 1.1 The structure diagram of Chapter 1 2

Schools

Infection risk assessment Spatial and temporal distribution of bioaerosol particles

Prediction of bioaerosol particle movement CFD simulation

1.1 Background of bioaerosol particle simulation in buildings 1.1.1 Impacts of Infectious disease outbreaks Public health has been threatened by the outbreaks of Severe Acute Respiratory Syndrome (SARS) in 2003, avian A/H5N1 in 2003, pandemic influenza A/H1N1 in 2009, Middle East Respiratory Syndrome (MERS) in 2012, avian A/H7N9 in 2013, Ebola virus in 2014 and MERS in 2015 (Braden et al., 2013; Tang et al., 2015). These outbreaks spread rapidly and globally to the human population due to globalization (Heymann, 2003; Smith, 2006). For example, the cumulative number of probable SARS cases was 8,098 with 774 reported deaths (WHO, 2003b). According to the World Health Organization (WHO), within weeks, SARS spread from Hong Kong to infect individuals in 37 countries (Smith, 2006). The emergence from time to time of pandemic outbreaks of these infectious diseases remains an ongoing threat to the human population.

For example, in 2003, several hundred cases of SARS were reported within 10 days in Amoy Gardens, one large housing estate in Hong Kong (Lee et al., 2003). Four out of seven tower blocks reported SARS cases on different storeys. The spread of SARS was suggested to be by way of airborne transmission of the SARS virus from the faeces of a single index case with diarrhoea. It was possible the virus–contaminated aerosols leaked out from the sewer pipes when travelling down from the apartment (Cheng et al., 2008). The viral particles might have been transmitted to other tower blocks by the wind and thermal plume, thereafter entering other apartments in the adjacent 3

towers through window openings if the kitchen or bathroom exhaust fans were operating, creating a negative pressure sink within.

In May 2015, the first case of MERS in South Korea was identified. Totally, 150 infection cases were reported including 18 deaths in June 2015 (Hui et al., 2015). Some confirmed cases were related to healthcare facilities and a remarkable proportion were associated with healthcare workers (HCWs) (17%). The outbreak in South Korea grew to be the second-largest MERS outbreak since the Saudi Arabia outbreak.

Both examples of outbreaks caused considerable morbidity and mortality. Mortality is undesirable to the individual and to communities. Morbidity causes large social and economic impact due to productivity, absenteeism and medical treatment costs. The cost of outbreaks is rarely quantified; however, the benefits of infection control are visible in terms of health and economic benefits (Zimlichman et al., 2013).

The potential health hazards of infectious diseases to the public are significant, not only from outbreaks (Burger, 1990; Douwes et al., 2003). For example, 18.8 billion cases of upper respiratory infections and 1.5 million cases of lower respiratory infection were reported in the Global Burden of Disease Study 2013 (Vos et al., 2015). For chronic respiratory disease (38.6 million cases), 26.1 million, 10.6 million and 0.5 million cases of chronic obstructive pulmonary disease (COPD), asthma and pneumoconiosis were recorded. 4

COPD is associated with a high dose of endotoxin, fungal spores or mycotoxins inhalation within several hours (Srikanth et al., 2008). The relation of bioaerosol (i.e. biological aerosol) to asthma has also been confirmed as exposure agents (Arya & Kaushik, 2013). Furthermore, attention to indoor air quality (IAQ) has increasingly had a large effect on public health, impacting health and productivity (Wong et al., 2009).

Inhalation of a bioaerosol through breathing can result in a respiratory infection if the inhaled microorganism cannot be removed by lung clearance, phagocytosis or dissolution and adsorption by a primary defensive mechanism. The microorganisms may stay in the lungs and begin either replicating or reactions. The results of such growth might be respiratory infection, allergic reaction or even toxic reaction (Cox & Wathes, 1995). However, respiratory irritation may be caused by microbial volatile organic compounds (MVOCs) produced from bacteria and fungi. Only infectious diseases caused by the transmission of pathogens to others (i.e. bioaerosol particles), are focused on in this thesis.

1.1.2 Spread of airborne infectious diseases Most outbreaks are caused by infectious diseases, especially airborne diseases, because they can be transmitted person to person over a short distance through the air (Wells & Stone, 1934). Inhalation and surface contact are two major routes of airborne transmission: 1) Respiratory droplets (i.e. bioaerosol particles) produced when an infected person coughs or sneezes are inhaled 5

into a respiratory tract of persons who are nearby; and 2) the bioaerosol particles can also be spread when a person touches a surface or object contaminated with infectious droplets and then touches his or her mouth, nose or eye(s) (Morawska, 2006).

The viability, infectivity, allergenicity, toxicity or pharmacological properties of diseases are related to their pathogen (i.e. bioaerosol) species (Fiegel et al., 2006). Airborne diseases can be categorized as communicable, noncommunicable and opportunistic diseases groups according to their epidemiological aspects as shown in Table 1.1. Table 1.1 Some airborne pathogens for the communicable, noncommunicable and opportunistic disease groups (Kowalski, 2006) Noncommunicable disease

Communicable diseases

Respiratory infection

SARS-CoV, MERS-CoV, Tuberculosis bacilli, Mycobacterium

Nonrespiratory infection

E. coli, Norwalk virus

Legionella, H5N1, Aspergillus, Bacillus anthracis,

Cladosporium,

Opportunistic disease Pseudomonas aeruginosa, Rhizopus spp., Staphylococcus aureus, Serratia marcescens, Streptococcus pyogenes Enterococcus faecalis, Staphylococcus epideris

Airborne communicable diseases are those transmitted between humans. All airborne human respiratory viruses are communicable except Hantavirus (Roy & Milton, 2004). For communicable respiratory diseases, airborne transmission may occur by means other than direct inhalation into a 6

respiratory tract, including direct and indirect contact with fomites left on a surface (i.e. surface contact). Secondary transmission (i.e. vertical transmission) may occur by respiratory pathogens liberated in a sneeze or a cough from an infected person (Nicas et al., 2005). For example, the SARS and MERS outbreaks in 2003 and 2013 were caused by airborne communicable respiratory disease pathogens (Joshi, 2013). This group of diseases rapidly spread as outbreaks without physical contact with infected people.

Airborne communicable non-respiratory diseases, including infections of the eyes, ears and skin, can be transmitted by direct or indirect contact with bioaerosol particles deposited on a surface, such as pathogenic infections of Escherichia coli via surface contact on infantile enteritis and diarrhoea (Neter & Shumway, 1950). Other common diseases, such as norovirus illness and hand, foot and mouth disease (HFMD) are kinds of airborne communicable non-respiratory diseases (Klausner et al., 2015; Liu et al., 2013).

A non-communicable pathogen is not transmissible from person to person, but it can be transmitted to humans from animals (i.e. zoonotic) or the environment (i.e. water or soil). Secondary infections among humans rarely occur (Cox & Wathes, 1995). For instance, Legionnaires’ disease is not transmitted between humans, but Legionella is common in soil and aquatic environments. The aerosolized Legionella particles can travel at least 6 km from a water system or cooling tower in air. Some zoonotic diseases

7

transmitted from infected animals to humans such as the H5N1 influenza virus, which can cause avian influenza or bird flu, are commonly found in birds. Human-to-human transmission of the virus has not yet been observed (Nikitin et al., 2014), but the chance of mutation of the virus is not guaranteed (Kawaoka, 2012).

Some opportunistic fungal diseases can be transmitted by taking advantage of a host (i.e. human) with a weakened immune system (i.e. immunodeficient) (Mims, 2001), but they cause no disease in a healthy host (i.e. human) with a normal immune system. An opportunistic infection, therefore, often happens in healthcare facilities as a nosocomial infection. Infection may be through inhalation, via deposited airborne pathogens on wounds or burns, or as a result of intrusive procedures being contaminated. For example, outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) are often reported due to its resistance to antibiotics such as penicillins and cephalosporins (WHO, 2002).

1.1.3 Infection in indoor environments Airborne infections occur both outdoors and indoors. Indoor airborne transmissions, however, are of higher concern, since people spend over 90 percent of their time indoors (Jenkins et al., 1992; Spengler & Sexton, 1983). Research on bioaerosol infections, which are caused by bioaerosol particles, has focused on indoor environments, especially buildings. The indoor environment in buildings plays a key role in the wide range of bioaerosol 8

infections, not only the longer occupancy time (Burger, 1990). Investigations have intensively studied bioaerosol emission in building environments and ventilation systems as shown in Figure 1.2 (Lai & Nazaroff, 2000).

Human (i.e. host)

Bioaerosol emission

Bioaerosol particles (i.e. agent)

Infection

Ventilation system

Building environment

Indoor air (i.e. enviro -nment)

Figure 1.2 Airborne infection factors in indoor environments

Bioaerosol emission for infection Microorganisms are easily aerosolized in ambient air by air movement caused by wind, coughing or sneezing due to their submicron size range (Kowalski, 2006). Aerosolization or atomization is a process of producing bioaerosol particles (i.e. airborne or droplets form) from breaking up the liquid surface by a high-velocity air jet (Sirignano, 2010). If microorganisms are contained in the liquid, they may become microbe-laden droplets and dry to become a droplet nucleus due to air jets or evaporation. For example, respiratory droplets are aerosolized from respiratory fluid in the respiratory tract by the high-velocity exhaled air jet (i.e. 50~100 m s-1) in the sneezing process (Mui 9

et al., 2009; Tang et al., 2012a). The aerosolization process also occurs in toilet flushing in addition to the expiratory activities of humans. Sometimes, water may not be involved in the process, for example, resuspension of deposited bioaerosol particles from a floor by humans walking or from a bedsheet by bed-making (Hathway et al., 2011; Kubota & Higuchi, 2012; Leung et al., 2013; You & Wan, 2014a). Many aerosolization processes occur in nature, such as waterfalls and wind that cause about 20~30% of bioaerosol particles in atmospheric particles in outdoor environments. Some potential bioaerosol particle emission sources in indoor environments have been identified in relation to infection control and are listed as follows (Morawska, 2006): 1) Human expiratory activities – breathing, speaking, coughing and sneezing; 2) Sanitary system operations – drainage systems, toilet flushing, showering, water taps; and 3) Ventilation system operations – cooling tower, filter, heat exchange Each of these processes generates bioaerosol particles of different characteristics in terms of their species, size, shape, emission rate, viability loss and initial speed.

Humans and their activities are linked to a number of processes resulting in the introduction of droplets containing bioaerosol particles into indoor air, including coughing and sneezing. Besides expiratory activities, bioaerosol particles can be emitted from skin or by walking (Bhangar et al., 2015). Sneezing and coughing are the two most commonly studied activities 10

regarding the expiratory activities in terms of emission species, amounts and velocities. The high frequency of breathing, talking and laughing is also significant in infection research. Most of the expiratory studies are associated with bacterial and viral groups that are communicable respiratory disease groups which were discussed in Section 1.1.2. The formation of expiratory bioaerosol particles inside the respiratory tract is less interesting than the development of bioaerosol particles in the air (Morawska et al., 2009).

High-speed photography was first used to investigate the velocities and sizes of the bioaerosol droplets from breathing and sneezing (Jennison, 1942; Jennison & Edgerton, 1940). By the use of an optical particle count (OPC) and a transmission electron microscope (TEM), about 80~90% of the droplets were found to be less than 1 µm in diameter (Papineni & Rosenthal, 1997). Besides expiratory bioaerosol particle studies, the dynamic airflow generated by expiratory activities is also important for infection control. Cough peak flow rate (L s-1), peak velocity time (s) and cough expiratory volume (L) have been investigated regarding the general steady airflow characteristics caused by coughing (Leiner et al., 1966; Smina et al., 2003). Selected literature on research techniques of expiratory activities is summarized in Table 1.2.

11

Table 1.2 Selected measurement methods of droplet size distribution of expiratory activities Measurement method or technique

Description

High-speed photography

The advantages of high-speed, stroboscopic-light photography for the subjects are intense illumination, by which extremely minute particles are mode visible, and short exposure-time, by which their motion is “stopped”.

Solid impaction

Using a celluloid-surface, glass slider, bond paper and chamber for collecting the respiratory particles, the samples are measured under a light microscope.

>0.1µm

Radioactive spores

Using labelled spores with a radioactive isotope, the distribution of the spores can be traced after the injection.

-

Optical particle counter

By counting the pulses of scattered light reaching the detector, particle number can be determined through a beam of light.

>1µm

Aerodynamic particle sizer spectrometer

Scanning mobility particle sizer Interferometric Mie imaging technique

Schlieren and Shadow-graph imaging method

Particle time-of-flight methods rapidly determine the number-weighted aerodynamic particle size distribution by accelerating the incoming particles through a well-defined flow field in the measurement zone of the analyzer. The transit time of the particle between two well-defined locations in the measurement zone is a monotonic function of aerodynamic diameter. The travelling time is measured by light scattering signals to associate the aerodynamic particle size. The technique determines particle size by scanning the electrical mobility of particles when traversing an electrical field. By the out-of-focus imaging of particles illuminated by a laser light sheet, the farfield scattering can be calculated by the Mie theory from the overlapping region of the interference fringe pattern. The method is based on high-resolution imaging with pulsed backlight illumination. The measurement volume is defined by the focal plane and the depth of field of the imaging system.

12

Measurement size

Reference

>1 µm

Jennison (1942), Jennison and Edgerton (1940) Duguid (1946), Papineni and Rosenthal (1997) Harper and Morton (1953), NCRP (1997) Fennelly et al. (2004), Xie et al. (2009), Lieberman (2006)

0.1~20µm

Morawska et al. (2009), Johnson et al. (2011)

2~1000nm

Yang et al. (2007)

1~100 µm

Chao et al. (2009)

-

Settles et al. (1995), Tang et al. (2012b)

The aerosolization process is also found in relation to vomiting wherein 107 virus particles per mL of vomit fluid have been reported in infected individuals (Barker et al., 2001). For example, a series of infections are suspected to have come from an SARS-infected person who vomited in the Metropol Hotel in the 2003 SARS outbreak in Hong Kong (Lee et al., 2003). In addition, the Norwalk-like virus outbreak in a school was found to have taken place after a student vomited in a classroom (Marks et al., 2003).

The virus content (i.e. 1012 virus particles per gram) was reported in human faeces (Barker et al., 2001). The mechanisms for aerosolization in sanitary facilities have been reported in toilet flushing and sewage transport in building drainage systems (Johnson et al., 2013a; Morawska, 2006). In general, little quantitative research has looked at the mechanism of sewage aerosolization through the above processes in terms of the size of the bioaerosol particles generated, and thus their fate in the air and the potential for spreading in sanitations (Cheng et al., 2010; Morawska, 2006; Wise & Swaffield, 2002).

Aerosolization in sanitary facilities may have a high potential for airborne infection with each toilet flush or use of water. The bioaerosol particles are inhaled or ingested by hand contact with deposited surfaces indirectly (Johnson et al., 2013b). For instance, bacteria samples have been reported in water, air and surfaces of hospital toilets (Newsom, 1972). Table 1.3 shows the selected seeding experiments of the aerosolization process after flushing.

13

Table 1.3 Summary of selected experiments of toilet seeding Toilet Type Cistern-fed, gravity-flow and mains-fed pressure value Wash-down Wash-down and siphonic

Conducted by

Year

Microorganism

Jessen (1955)

1955

Serratia Marcescens

1959

Serratia Marcescens

1966

Escherichia coli

Darlow and Bale (1959) Bound and Atkinson (1966)

Wash-down and double-trap siphonic

Newsom (1972)

1972

Escherichia coli, Salmonella typhimurium, Shigella sonnei, Proteus mirabilis, Serratia marcescens, Klebsiella aerogenes, Pseudomonas aeruginosa and Achromobacter spp.

Siphonic gravity-flow

Gerba et al. (1975)

1975

Escherichia coli and MS2 bacteriophage

Siphonic gravity-flow

Barker and Bloomfield (2000)

2000

Salmonella enteritidis and PT4

Siphonic gravity-flow

Barker and Jones (2005)

2005

Serratia marcescens and MS2 bacteriophage

Siphonic gravity-flow

Best et al. (2012)

2012

Clostridium difficile

In a drainage system, the leakage of bioaerosol particles outside the system should be prevented if the system is maintained properly (Wise & Swaffield, 2002). However, the leakage of bioaerosol particles from a drainage system may cause the outbreak of a disease such as the spread of SARS-CoV in Amoy Garden, which was described in Section 1.1 (Sobsey & NMeschke, 2003; WHO, 2003b). Dry floor traps and bathroom fans were identified as the major contributory causes of this outbreak by numerical simulation studies (Cheng et al., 2013; Gormley et al., 2012; Jack et al., 2006). Some drainage designs 14

have been recommended for the avoidance of bioaerosol particles leakage within the network and sewer pipework; especially, predictions of airflow, transient network pressures and trap seal retention level are suggested to be incorporated in the design stage (Jack, 2006). The infection risk analysis for drainage systems has been investigated further in high-rise residential buildings through failure mode effects analysis (Cheng et al., 2008). The reuse of greywater has been quantified by experiments (Benami et al., 2016).

The aerosolization processes of bioaerosol particles have been observed in other sanitary activities, such as pavement cleaning (Seidl et al., 2015), showering (Dennis et al., 1984), water taps (Carson, 1996), water systems (Kool et al., 1999) and hot spa springs (Armstrong & Haas, 2007b), especially the Legionella group in these systems (Ellis, 1993).

The aerosolization process of bioaerosol particles found in ventilation systems such as cooling towers is suspected to be a source of many Legionella outbreaks (Buse et al., 2012; Fraser et al., 1977), although Legionella is found less in cooling towers nowadays due to legislation. Other microorganism species, however, are still found in ventilation systems such as Micrococcus spp. Staphylococcus spp., and Aspergillus spp. Growths of microorganisms inside ventilation systems have been found in many places such as mixing chambers, cooling coils, air filters, air ducts and heat exchangers (Bluyssen et al., 2003; Chow et al., 2005; Hugenholtz & Fuerst, 1992; Schmidt et al., 2012; Zuraimi, 2010). The high potential for aerosolization of bioaerosol

15

particles is correlated with the high airflow rate in ventilation systems (Zuraimi, 2010). The aerosolized particles are then dispersed through the system, and eventually arrive in occupied rooms through the air distribution system. Such bioaerosol emissions from ventilation systems may create an impact on infection control strategy, especially if the ventilation system is a major control mechanism of infection control. No study is available to assess quantitatively bioaerosol emissions from ventilation systems and other relevant factors, like ventilation rate (VR) and numbers of occupants. The maintenance plan and disinfection devices in the installation of ventilation systems should be included in the effect of this potential bioaerosol emission in ventilation system design.

Building environments for infection For building environments, potential emission source, occupant density and occupant activity pattern differ among building types. To optimize the performance of infection control, a ventilation system design has to be specified based on the characteristics of the building type. Table 1.4 summarizes some key studies on building types. In particular, some building types with high potential risk have been intensively studied. Standards and guidelines have been developed to specify the ventilation system design in the building type such as healthcare and transportation facilities (ASHRAE, 2013c; CDC, 2003, 2004).

16

Table 1.4 Selected literature on building type related to bioaerosol infections

Healthcare facilities

Offices and factories

Transportation facilities

Washrooms, water systems, and drainage systems

Aircraft, trains, buses, vehicles and ships

Others

Isolation rooms, and elderly centres

Sanitary facilities

Area

Workplace

Building types

Educational and residential facilities

Literature Nosocomial infection or hospital acquired infection (HAI) is a special type of infection which is acquired in healthcare facilities such as hospitals, clinics and specialized care centres. The facilities have a very high infection risk since some patients could be emission sources of the bioaerosol particles and have weak immune systems, so they are infected easily by opportunistic diseases. For example, the high infection rates of HCW of SARS (20.5%) and MERS (17%) outbreaks were recorded in Hong Kong and South Korea (Jack, 2015; Lau et al., 2004). Healthcare facilities consist of various types of function areas such as operating suites, AIIR, intensive care units (ICUs), burns wards, general patient wards, emergency rooms, laboratories and offices. These areas have different ventilation and infection control requirements (Tang et al., 2015; Tang et al., 2006). Some areas are identified as high-risk areas such as operating suites, AIIR, ICUs and burns wards due to the high infection risk for both patients and HCWs (Noakes et al., 2015). Some ventilation suggestions, such as high VR (i.e. >12 ACH), pressurization control, local exhaust system, and the use of HEPA filters and UVGI devices are recommended for these areas (ASHRAE, 2013c; CDC, 2007). A warm and humid environment promotes microorganism breeding in sanitary facilities such as toilets, showers and bathrooms (Wells, 1943). Microorganisms have been highly reported in washrooms, especially on seats and under the flushing rim of water closets (WC) (Flores et al., 2011; Irvine & Robertson, 1964; Mendes & Lynch, 1976; Newsom, 1972). These microorganisms can survive up to six hours in washrooms (Newsom, 1972). Potential infection paths of related diseases have been identified via ingestion, splashing during defecation, surface contact and direct inhalation (Barker & Jones, 2005; Burgess, 1963; Darlow & Bale, 1959). Bioaerosol particles in washrooms sourced from toilet flushing through splashing and frothing were identified as the key medium of infection (Burgess, 1963). Formation of bioaerosol particles by flushing has been investigated experimentally (Jessen, 1955). To reduce the infection risk by flushing, some research of the reduction of microorganisms in water seals and bioaerosol emission has been conducted, such as using disinfectants, closing the lid and using siphonic gravity flow WC (Barker & Bloomfield, 2000; Johnson et al., 2013b). However, these studies only focused on source control rather than a source removal strategy by an exhaust ventilation system which is used to remove odour and humidity in the facilities. The existing exhaust ventilation design practices only focus on the required VR for odour dilution and fresh air intake in a washroom (ASHRAE, 2013b; Chung et al., 1997; CIBSE, 2005; Rock & Zhu, 2002; Tung et al., 2009; Tung et al., 2010). Long working hours and high occupant densities in workplace environments should not be neglected for infection research. Because of the daily occupancy and interaction of occupants within the environments, many infections are regularly transmitted inside these building types (Morawska, 2006). Occupational safety is also associated with the studies (Corzine et al., 2003; HSE, 2003). Offices and factories, except healthcare facilities for HCWs, are two common types of workplace environments for airborne infection research. Various respiratory symptoms and health outcomes of biocontaminated office occupants have been reported, often implicated by elevated bioaerosol exposure (Wu et al., 2005). Microorganism growths inside ventilation systems have been reported, such as cooling coils, mixing chambers, humidifiers and heat exchangers (Bluyssen et al., 2003; Chow et al., 2005; Hugenholtz & Fuerst, 1992; Schmidt et al., 2012; Zuraimi, 2010). The existence of these microorganisms inside ventilation systems can endanger office occupants. The high potential of the microorganisms to be aerosolized to form bioaerosol particles is due to the high airflow rates inside the ventilation systems (Zuraimi, 2010). The bioaerosol particles are then dispersed in the system, and eventually into the occupied offices through the air distribution system. Such emissions from ventilation systems may create an impact for infection control in offices, especially since over 90% of commercial buildings in Hong Kong are equipped with ventilation systems (Chow et al., 2005). Two typical bioaerosol sources in factories are found in hazardous environments, such as work in landfills or farms, and working processes that will generate bioaerosols such as the water spray in wastewater treatment plant. For poultry houses, zoonotic diseases (i.e. transmissible from animals to humans), such as avian influenza and swine disease, have to be considered for infection control. In addition, foodborne and waterborne diseases may be aerosolized during the manufacturing process such as in wastewater treatment (Gormley et al., 2014; Sanchez-Monedero et al., 2008) and textile plants (Eduarda & Heederik, 1998; Su et al., 2002). Local exhaust systems are mostly recommended for factories. However, the ventilation system should be specified for each industry in terms of bioaerosol emission species, loading and locations. Transportation systems are a hub connecting people in different places, and they can also be a hub for spreading diseases over a large region or even on a global scale (Xu et al., 2013). The infection control research on transportation systems, therefore, has been intensively studied, especially aircraft cabins (Gupta et al., 2011; Liu et al., 2012; Sze To et al., 2009; Wan et al., 2009; Yin et al., 2012). Aircraft, trains, buses and other cabin environments can pose infection risks similar to those of other indoor environments except the typically smaller total volumes often increase the risk due to proximity and make air distribution and air cleanliness a more critical factor. Several infection control challenges in the cabin environment of transport vehicles have been addressed (Liu & Zeng, 2012) regarding the necessity for ventilation, comfort and health in commercial aircraft, cars, trains and buses. These challenges can be summarized as involving three issues. Due to the small spaces of the vehicles, a cabin is often highly packed with passengers, the high density of the cabin environment makes the air distributions more complex than in a building, the weight of the air ventilation system is very sensitive to the consumption of energy, and the thermo-fluid boundary condition of the occupied zone is much nearer to that in a building due to the small spaces in cabin environments (Liu & Zeng, 2012). Educational and residential buildings are the two common building types to be investigated in addition to the above building types. For educational facilities, universities and kindergartens are the two most investigated facilities. (Yamamoto et al., 2015). Residential buildings, especially high-rise residential buildings, have been studied under naturally ventilated conditions (Gao et al., 2008, 2009; Mao & Gao, 2015). However, a ventilation system design for infection control has still to be noted in other building types regarding potential bioaerosol emission sources and exhaust mechanisms.

17

Ventilation systems for infection Ventilation systems are important for the control of bioaerosol infection (Kowalski, 2006). These systems provide a significant engineering control measure in indoor air as well as providing fresh air and a comfortable IAQ and thermal environment for the occupants as illustrated in Figure 1.3 (Azimi & Stephens, 2013). Fresh air

Filtration Ventilation

Occupant (i.e. susceptible)

Recirculation

Infected person Dispersion Deposition Figure 1.3 Bioaerosol transmission with ventilation systems (Azimi & Stephens, 2013)

Some guidelines specify certain ventilation configurations and parameters for minimizing the infection for various building types such as healthcare facilities and air transportation (ASHRAE, 2013c; CDC, 2003, 2004). These design parameters include the selection of air distribution system types, VRs and the location of air intakes. Table 1.5 summarizes the key research on how these ventilation parameters affect bioaerosol infection.

18

Table 1.5 Selected literature on ventilation systems related to bioaerosol infections System

Area

Air distribution system type

Mixing Displace -ment

Natural Pressuri -zation Local exhaust Personal ized

Ventilation parameters

Ventilation rate Temperature & RH Recircul -ation ratio

Others

Air flow pattern Filtration Disinfec -tion devices

Literature Mixing ventilation (MV) systems are the most used commercial ventilation systems. The design and operation of the systems are simple and cost-effective if well designed due to the uniform distribution of fresh air concentration and temperature by the mixing of the supply and room air. The dilution effect of the mixing ventilation for the widespread TB outbreaks in the 1960s was investigated (Nardell et al., 1991; Riley et al., 1959; Rudnick & Milton, 2003). Some guidelines and standards recommended minimum ventilaton rates to reduce infection risks for mixing systems, especially in healthcare settings (ASHRAE, 2013b, 2013c; CDC, 2003). Displacement ventilation systems supply air at floor level and plume the pollutants to the upper levels of a room. Compared with mixing ventilation, the system is more efficient in removing gaseous pollutants than the mixing systems (ASHRAE, 2013a; Chen & Glicksman, 2003). However, a longer elapsed time and more depositions for bioaerosol particles have been reported due to the buoyancy force not being sufficient to lift the bioaerosols up to exhaust level (Mui et al., 2009). Bioaerosol particles in the breathing zone are critical for infection control instead of diluting the total air volume (Bolashikov & Melikov, 2009; Pantelic et al., 2009b). Natural ventilation is commonly found in residential buildings. The meteorological conditions are critical for designing naturally ventilated buildings rather than the building envelope, window opening ratio, cross-ventilation or stack effect that is driven by the wind or a buoyant force (Liu & Cheung, 2009). Outdoor infiltration of bioaerosol to indoor space is not controllable by natural ventilation (Gao et al., 2009; Mao & Gao, 2015). The Amoy Gardens outbreak of SARS in 2003, for example, was dispersed by natural ventilation from one block to the others via the single-sided window design (Gao et al., 2009). A trickle vent design was suggested to reduce the infection risk for naturally ventilated buildings (Mao & Gao, 2015; Yang & Gao, 2015). Pressurization control is frequently found in some building types such as healthcare facilities and biological research laboratories. Positive pressure reduces infiltration of bioaerosols into an indoor space and is recommended in operating rooms, while negative pressure is suggested in AIIR to avoid bioaerosols leakage (ASHRAE, 2013c). Local exhaust ventilation is the most effective contaminant removal strategy if it is located near the emission source. Adding a hood sometimes improves the removal performance of exhaust regardless of the distance (Goodfellow & Tahti, 2001). Sterilization and autopsy rooms are recommended to use the system to control the infection risk (ASHRAE, 2013c; Dygert & Dang, 2012; Leung et al., 2006). The system is also used in laboratories, kitchens and industrial factories for removing chemical pollutants as well. Personalized ventilation (PV) supplies clean air to occupants individually. The gaseous pollutant concentration in a breathing zone could be reduced by 2 to 50 times compared with mixing ventilation (Cermak & Melikov, 2006; Melikov, 2004). Improvement of occupant protection for bioaerosol infection has also been reported as a supplemental method for infection control. Also, the system has been suggested in hospitals and aircraft cabins (Qian et al., 2006; Wan et al., 2009). However, the applications of this suggestion are limited to theatres, cinemas, lecture halls and aircraft cabins where the chairs are connected to fresh air pipes (Pantelic et al., 2009b). Ventilation with outdoor air is intended to provide fresh air to occupants and remove pollutants emitted from indoor sources. The association between VR and occupant health has been reported (Sundell et al., 2011). The importance of the VR has also been reviewed in regard to bioaerosol infection risk (Carrer et al., 2015; Rudnick & Milton, 2003; Sundell et al., 2011). In practice, the actual percentage of air which is exchanged inside a room as the ventilation effectiveness of that air is not perfectly mixed (Li et al., 2008; Wang et al., 2008). However, the VR can still be a surrogate indicator to understand the overall room ventilation performance regarding IAQ and infection control. Environmental conditions such as air temperature and relative humidity are strongly associated with the indoor bioaerosol concentration level (i.e. bacteria and fungi level) (Chan, 2012; Tang, 2009; Wong et al., 2008). These thermal and humidity conditions may accelerate the germination, growth and survival rate of bioaerosol particles. For example, lipid-enveloped human coronavirus 299 E remains alive with a half-life of 67 hours at a RH of 50% and air temperature of 20°C The long survival time of bioaerosol particles increases the infection risk due to the chance of re-entry to the room or entry to other rooms by the recirculation process. For energy saving, room air is usually recirculated and mixed with a certain amount of fresh air (i.e. fresh and recirculation rate) in an air handling unit (AHU) room that reduces the energy consumption in cooling and dehumidifying the fresh air. However, bioaerosol particles may re-enter the room with the recirculated air. Bioaerosol particles could be looped in the room until the end of their lifetime by the recirculation process, which enhances the overall infection risk. On the other hand, the infectious risk could be further spread to other rooms that connect with the same AHU (Kupferschmidt, 2015). Reducing the fresh air recirculation ratio might encourage viruses spreading in indoor space (Mendell et al., 2002). Some suggestions of airflow patterns should also be of concern in the design of ventilation systems for infection control and prevention (Morawska, 2006; Sundell et al., 2011). In general, a downward airflow pattern was found to be the best for controlling the lateral dispersion of bioaerosol particles, especially for particles with a diameter of less than 45 µm (Chao & Wan, 2006). The spatial location of the emission source (i.e. coughing or sneezing) and the subjects was identified as a key factor in the dispersion and deposition of bioaerosol particles (Sze To et al., 2008). Regarding the difficulties in quantifying the airflow pattern for a ventilation system, CFD simulation has been introduced to simulate the airflow in bioaerosol infection studies (Sze To, 2010). HEPA filters are commonly recommended to be used in a bioaerosol infection control environment, although they are expensive and have a higher energy consumption due to the pressure drop from the filter. The HEPA filters are effective at capturing most bacterial and fungal species (1, then the outbreak or epidemic is propagating. If Routbreak1 Routbreak 3 µm (Douglas, 1975), 𝐼𝐷50 = 223.5 × 𝑇𝐶𝐼𝐷50

[2.14]

The relationship between the intake dose in the alveolar region (i.e. target infection site) and the inhaled dose from the surrounding air via nose was suggested in the estimation of the airborne infection risk of tuberculosis bacilli (Nicas, 1996). 𝑁𝑖𝑛𝑡𝑎𝑘𝑒 = 𝑁𝑖𝑛ℎ𝑎𝑙𝑒 × 𝑓𝑎𝑙𝑣𝑒𝑜𝑙𝑎𝑟 =

𝑁𝑖𝑛𝑓𝑒𝑐𝑡 × 𝑅𝑠𝑟𝑐 × 𝑄𝑝𝑢𝑙𝑚𝑜𝑛𝑎𝑟𝑦 × 𝑡𝑒𝑥𝑝𝑜 𝑄𝑟𝑜𝑜𝑚

[2.15]

× 𝑓𝑎𝑙𝑣𝑒𝑜𝑙𝑎𝑟

Where Rsrc is the airborne tuberculosis bacilli emission source rate (min-1), falveolar is the deposition fraction of infectious particles in the alveolar region. fsite is the fraction of the intake dose that reaches the target infection site. For example, the alveolar region may be the target infection site for respiratory tract infections or airborne diseases mostly. Equation (2.15) predicts the intake dose Nintake of tuberculosis bacilli by multiplying the falveolar and inhale dose Ninhale, where the inhale dose Ninhale could be estimated from Rsrc, Ninfect, Qroom and texpo. Combining Equations (2.12) and (2.15), the infection probability for airborne is derived in Equation (2.16). 𝑃𝑟𝑜𝑏𝑖𝑛𝑓𝑒𝑐𝑡 = 1 − exp (−𝑓𝑑𝑜𝑠𝑒

[2.16] 𝑁𝑖𝑛𝑓𝑒𝑐𝑡 × 𝑅𝑠𝑟𝑐 × 𝑓𝑎𝑙𝑣𝑒𝑜𝑙𝑎𝑟 × 𝑄𝑝𝑢𝑙𝑚𝑜𝑛𝑎𝑟𝑦 × 𝑡𝑒𝑥𝑝𝑜 × ) 𝑄𝑟𝑜𝑜𝑚

50

Equation (2.16) has demonstrated the application potential of the doseresponse model in assessing the infection risk for airborne transmission. The equation is similar to the Wells-Riley model (i.e. Equation (2.10)) with quanta generation rate Rquanta replaced by Rsrc× falveolar. Sometimes, falveolar and fdose will be combined to falveolar to simplify the equation.

In addition, the adequacy of using dose-response models in assessing airborne infection risk for another disease was demonstrated in a Legionnaires’ disease study (Armstrong & Haas, 2007a; Armstrong & Haas, 2007b; Armstrong & Haas, 2008). The study performed animal tests for extrapolating the interspecies of the dose-response curve of Legionnaires’ disease, especially in low dose conditions (Armstrong & Haas, 2007a). The spatial variation of the exposure level at a steady state was estimated by a near field model for hot spring spa outbreaks (Armstrong & Haas, 2007b). Then the risk assessments were predicted by the dose-response curve and the exposure level and validated by comparing the estimated risk from the reported attack rate in the outbreaks (Armstrong & Haas, 2007b; Armstrong & Haas, 2008).

These studies also highlighted the potential of the dose-response model in assessing the infection risk of exposure to bioaerosol particles instead of the particles generated by an infector since the model separates the exposure and dose-response assessments for the infection risk assessment framework as shown in Figure 2.10 (Nicas, 1996; Sze To, 2010).

51

Bioaerosol Emission

Bioaerosol dispersion

Exposure assessment

Bioaerosol concentration

Exposure time

Bioaerosol exposure

Doseresponse assessment

Dosimetry factors

Intake dose Dose-response curve

Probability of infection

Risk Infection assessment

Risk characterization

INFECTION RISK Figure 2.10 Exposure and dose-response assessments for risk infection framework (Nicas, 1996; Sze To, 2010)

Under this framework, various emission mechanisms, such as a multiple of coughs, could be expressed as an infectious source strength Rsrc in terms of 52

cough frequency, bioaerosol particle concentration in respiratory fluid and the amount of bioaerosol particles generated in coughs (Nicas et al., 2005). Furthermore, by this risk assessment, a spatial and temporal distribution of the exposure assessment model for virus has been proposed to evaluate the infection risk of a non-uniform distribution ventilated space in Equation (2.17) (Sze To et al., 2008; Sze To et al., 2007) 𝑁𝑖𝑛ℎ𝑎𝑙𝑒 (𝑥, 𝑡𝑒𝑥𝑝𝑜 ) = 𝐶𝑜𝑛𝑐𝑏𝑝

[2.17] 𝑡𝑒𝑥𝑝𝑜

𝑓𝑏𝑝 (𝑥, 𝑡)𝑣𝑖𝑎𝑏(𝑡)𝑑𝑡

× 𝑄𝑝𝑢𝑙𝑚𝑜𝑛𝑎𝑟𝑦 ∫ 0

Where Ninhale(x, texpo) is the exposure level of the bioaerosol particles at location x during the exposure time interval (Plaque forming unit (PFU)), Concbp is the bioaerosol particle concentration in the respiratory fluid (PFU mL-1), viabbp(t) is the viability (i.e. survival) rate function of a bioaerosol particle (%), fbp(x, t) is the volumetric fraction of bioaerosol particles (mL) over air (L) at the location (mL L-1 of air). fbp(x, t) can be determined by measurements or CFD simulation. Generally, Equation (2.17) will provide more practical exposure estimations, but it is more timeconsuming than obtaining the exposure level based on the well-mixed assumption or other simple models. The equation could be further incorporated with particle size bins to form an infection risk model for polydispersed bioaerosol particles regarding their species. The modification is especially suitable for parametric studies on the effect of environmental 53

control from multiple emission sources, such as ventilation strategy or airflow pattern on the infection risk for non-uniform distribution environments (Sze To et al., 2008). Furthermore, the equation was examined by crossreferencing the most probable quanta generation rate of the Wells-Riley model with maximum likelihood estimation in a Varicella (i.e. chickenpox) outbreak (Chao, 2011).

The dose-response model provides the bioaerosol infection risk assessment for indirect contact (i.e. surface contact) transmission (Atkinson & Wein, 2008; Nicas & Jones, 2009; Sze-To et al., 2014). A relative contribution between the direct and indirect contract transmission of influenza virus exposure routes was estimated by the inhaled (i.e. airborne) and ingested (i.e. finger contact from virus-contaminated surfaces, then finger contact to the eyes, nostrils or lips) doses. A significant contribution was reported by hand contact with facial membranes (i.e. 31% contribution) (Nicas & Jones, 2009). On the other hand, the effect of surface material was examined in aircraft and healthcare environments. The results showed that the indirect contact risk of non-fabric surfaces may be higher than that contacting on fabric surface by a thousand times. In addition, reducing the contact rate of the surfaces is relatively more effective instead of increasing the VR for the infection control.

Furthermore, the dose-response model also provides an alternative infection risk assessment for airborne transmission rather than the Wells-Riley quantum model. Various types of the emission source, such as coughing, hot 54

spring spa, toilet flushing or drainage system, could be assessed by the doseresponse model due to the separation of the exposure and dose-response assessments in the infection risk framework (Nicas, 1996). An emission source could be expressed more flexibly in terms of emission species, amount, rate, size, shape and velocity. The estimation of the input parameters for the emission, therefore, is an important factor for the exposure assessment. The prediction of the spatial and temporal distributions of bioaerosol particles are critical for both infection risk models, especially for low dose conditions.

2.4 Prediction of bioaerosol particle movement For the transport model, the bioaerosol particle diameter and falling time were empirically correlated to estimate the possible travel distance from the infector (Wells et al., 1948; Wells et al., 1942). The gravity force, settling velocity and evaporation time have been associated to understand the dispersion and deposition of bioaerosol particles as the Wells evaporationfalling curve in Figure 2.4. Only particles less than 100 µm was concluded to be significance for airborne transmission in stagnant air, although the critical diameter of the 100 µm is questioned (Xie et al., 2007). The study led the empirical transport model to study the dispersion and deposition of bioaerosol particles instead of outdoor field measurement in the previous aerobiology studies for a small travel distance. (Buller, 1909; Burge, 1995; Marshall, 1904). In addition, the study suggests the physical properties (i.e. diameter) of bioaerosol particles are critical factors on particles spreading and the ventilation system design for infection control. Various empirical models 55

have investigated the particle dynamics from an infected person to a new host. For example, a study of air duct shapes (i.e. a long straight, an ‘S’ shape and a ‘U’-shape) suggested the bends of the duct could isolate the larger bioaerosol particles instead of the smaller particles (Andrewes & Glover, 1941).

The understanding of size distribution of bioaerosol particles related to various release mechanisms and their subsequent transport was still limited. Modeling of droplet transport based on the Gaussian plume model was developed for bioaerosol particles transport for dispersion in atmospheric conditions (Gregory, 1961). The bioaerosol emission, evaporation, dispersion and survival functions were included in the model. The findings provided a 30 m accuracy from the source (Lighthart & Kim, 1989). Advection and advection-diffusion models were proposed to improve the accuracy of the mathematical analytical transport model (Aylor, 1986). However, these above models could not provide accurate results for indoor environments. A mass balance model was used to predict the deposition, resuspension and penetration of bioaerosol particles (i.e. 1 to 5 µm) for a residence (Thatcher & Layton, 1995). A multi-zone mass balance model was further developed to solve the indoor bioaerosol concentration level with the Wells-Riley quantum model for infection risk assessment (Noakes & Sleigh, 2008). However, a simple geometric structure and no spatial variation of a zone limit the applicability of the multi-zone mass balance model.

56

Since the SARS outbreak in 2003, CFD approach has been used and popularized to simulate the indoor bioaerosol particle transports and recognized the possible transmission routes in the SARS outbreak, especially in healthcare facilities and transportations, as discussed as in Section 1.1.3 (Chao & Wan, 2006; Lee et al., 2003; Li et al., 2005). A two-phase flow (i.e. continuum-dispersed) simulation resolves spatial and temporal solutions for airflow and bioaerosol particle dynamics to evaluate the infection risk assessment by the Wells-Riley quantum or dose-response models (Ishii & Hibiki, 2006). This framework is a hybrid approach utilizing a reference description (i.e. Eulerian or Lagrangian) for continuum airflow fields and another reference description (i.e. Eulerian or Lagrangian) for a bioaerosol particle (Subramaniam, 2013; Zhang & Chen, 2007). In the framework, the airflow field of the ventilated space was solved by continuum phase simulation in terms of air velocities, pressures and turbulences. Then bioaerosol particles were predicted in dispersed phase simulation since the droplet nuclei form of bioaerosol particles were supposed to be dried residue. Prediction movement of individual particles is important for understanding the spatial risks associated with an identified airborne pathogen.

For the continuum phase simulation, DNS and LES approaches provide very detailed simulation results in airflow and particle movements. However, both approaches require heavy computing resources and are very time-consuming. The RANS model is commonly used due to reasonable computational demand, especially in RNG k-ε turbulence model (Chen et al., 2013), 57

although other simulation approaches, such as the lattice Boltzmann method, DES model and population balance model, were proposed for predicting the bioaerosol particle movements (Chen et al., 2013; Fu et al., 2015).

The RANS model is popular for bioaerosol particle simulation in ventilated enclosures since the better performance of near-surface turbulence has been demonstrated for indoor surfaces (Lai & Nazaroff, 2000). Standard and RNG k-ε turbulence have been suggested for use in ventilated environments (Chen, 1995). Various air distributions systems for infection control were compared based on the results of CFD simulation such as displacement and personal ventilation systems (Tham & Pantelic, 2011; Wan & Chao, 2007). The CFD simulation was also used in ventilation system design for healthcare facilities and transportations (Tang et al., 2015; Wan et al., 2009). Various air distribution systems and building types could be predicted and recognize the ventilation system design problem for infection control and prevention by CFD simulation (Pantelic & Tham, 2011; Wan & Lin, 2015).

The interaction between occupants and room environment for the bioaerosol particle transport has been investigated by CFD simulation recently. Due to the dynamic airflow generated by humans or objects (i.e. door) movements, the potential bioaerosol emission or transient breakdown of infection control was reported (Han et al., 2014; Tang et al., 2005) The interaction between human and bioaerosol particles movement have been investigated by dynamic mesh technique to simulate the human walking and door opening in an 58

isolation room (Shih et al., 2007). Due to the heavy loading of the method, an interface-tracking technique was suggested by setting up some interfaces between the human and space (Brohus et al., 2006; Mazumdar & Chen, 2007). The results were more accurate than an interface-capturing method since the fixed mesh reduces the computation power (Tezduyar, 2006). By adding sources of momentum to some zones for human walking, a simple steady state model has been proposed to simulate the human movement in an isolation room (Brohus et al., 2006; Hathway & Papakonstantis, 2015). The results provide a relatively simple method to simulate the influence of movements for a transient condition in an acceptable approximate way. In addition, the re-suspension by human movements was also simulated (Kubota & Higuchi, 2012; Leung et al., 2013; You & Wan, 2014b). The effect of the interaction between occupants and indoor air on the bioaerosol particle transport is a key issue in the future development of CFD simulation.

2.5 CFD simulation for bioaerosol particle transport models For bioaerosol particle movement simulation, the drift-flux model (DFM) and discrete phase model (DPM) are two commonly used approaches in dispersed phase CFD simulation (Holmberg & Li, 1998). Since the size of the bioaerosol particles less than 100 µm, it was assumed to be airborne (Wells, 1934), the gas-gas two-phase CFD simulation was used to estimate bioaerosol particle dispersion by a DFM under a Eulerian-Eulerian framework (Hibiki & Ishii, 2003). However, the behaviors of bioaerosol particles in air differ from gas molecules for the same boundary condition. For example, bioaerosol 59

particles maintain higher momentum (i.e. velocity) along an airstream when compared with the rapid momentum decay of gas molecules (Wan, 2006).

In addition, molecular diffusion is an important transport mechanism for the gas, but the diffusion could be neglected for bioaerosol (Nicas et al., 2005; Pantelic et al., 2009a; Wan, 2006). DPM have been proposed to predict bioaerosol particle movement by the force balance of the interactions with the continuum phase. This approach is sometimes regarded as the EulerianLagrangian framework. Compared with the Eulerian-Eulerian approach, the Eulerian-Lagrangian framework provides a more natural description of the actual physical phenomena since each particle is considered individually. Stochastic fluctuations are complemented with statistical turbulent dispersions by discrete random walk (DRW), although the instantaneous turbulence quantities of the dispersed phase cannot be solved by the RANS model (Wan, 2006). In this section, both models are discussed to understand the fundamental mechanisms of bioaerosol particle transport model.

2.5.1 Drift-flux model (DFM) for bioaerosol transport The drift-flux model (DFM) is a mixed-flow model in which the focus is on the relative motion between the phases rather than the motion of the individual phase motion (Hibiki & Ishii, 2003; Zuber & Findlay, 1965). The dynamics between two phases, therefore, can be expressed by the mixture momentum equation. Four field equations (i.e. mixture continuity, mixture 60

momentum, mixture energy and gas continuity) are only involved in DFM instead of six equations in a two-fluid model (Hibiki & Ishii, 2003).The simplicity of the model makes it widely used in two-phase flow applications such as bubble, slug, droplet, annular and fluidized bed simulations (Ishii, 1977). Furthermore, the model provides an inter-species simulation for various bioaerosol particle sizes.

The void propagation function of the DFM avoided the non-continuum distribution of bioaerosol particles where low particle volume fraction was found in CFD simulation of the particle transport for an indoor environment (Holmberg & Li, 1998). Further, the application of DFM was used for SARS-CoV simulation for the nosocomial outbreak in Hong Kong (Li et al., 2005). Modelling of the dispersed phase can be solved by adding an equation of mixture momentum in Equation (2.18). 𝜕(𝜌𝑎𝑖𝑟 𝑓𝑏𝑝,𝑏𝑖𝑛 ) + ∇ ∙ (𝜌𝑎𝑖𝑟 (𝑣𝑏𝑝,𝑏𝑖𝑛 + 𝑣𝑠𝑒𝑡𝑡𝑙𝑒,𝑏𝑖𝑛 )𝑓𝑏𝑝,𝑏𝑖𝑛 ) 𝜕𝑡 𝜇𝑒𝑓𝑓 =∇∙( ∇𝑓 ) + 𝑆𝑟𝑐𝑏𝑖𝑛 𝜎𝑏𝑖𝑛 𝑏𝑖𝑛

[2.18]

Where fbp,bin is the volumetric fraction of bioaerosol particles (mL) over air (L) of the dispersed phase at the binth size bin, vbp,bin is the velocity of bioaerosol particles at the binth size bin, σbin is the dispersed phase diffusivity (m2 s-1) at the binth size bin, Srcbin is the source term at the binth size bin, µ eff is the effective viscosity of the air (kg m-1 s-1), vsettle,bin is the settling velocity at the binth size bin, which could be obtained by Equation (2.19): 61

0.5

𝑣𝑠𝑒𝑡𝑡𝑙𝑒,𝑏𝑖𝑛

4 𝑔𝑑𝑝 𝜌𝑝 − 𝜌𝑎𝑖𝑟 =[ ] 3 𝐶𝑑𝑟𝑎𝑔,𝑏𝑖𝑛 𝜌𝑎𝑖𝑟

[2.19]

Where dbp,bin is the diameter of the particle in the binth size bin (µm) and ρbp is the particle density in the binth size bin (kg m-3), Cdrag,bin is the drag coefficient in the binth size bin, ρair is the air density (kg m-3). The turbulent effect has been added to improve the equation for near surfaces diffusion for the indoor environment by modifying the dispersed phase diffusivity σbin (m2 s-1) (Shimada et al., 1996). Furthermore, the Stokes-Einstein diffusivity was included into the diffusivity σbin for improving the deposition of human airways simulation (Wang & Lai, 2005). The effects were combined into Equation (2.20) (Wan, 2006).

𝜎𝑏𝑖𝑛 =

𝐾𝐵𝑜𝑙𝑡𝑧𝑚𝑎𝑛𝑛 𝑇𝑎𝑖𝑟 𝑓𝑠𝑙𝑖𝑝 𝜇𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡 + 𝜌𝑎𝑖𝑟 𝑆𝑐ℎ𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡 3𝜋𝜇𝑎𝑖𝑟 𝑑𝑏𝑝

[2.20]

Where µ turbulent is the turbulent viscosity of the air (kg m-1 s-1), Schturbulent is the turbulent Schmidt number, KBoltzmann is the Boltzmann constant (i.e. 1.38×10-23 J K-1), Tair is the temperature (K), fslip is the Cunningham slip correction factor. However, the gaseous mixture of the DFM simulation could not reflect the behaviours of bioaerosol particles in the air, such as the higher velocity of particles (i.e. momentum) maintain better along the pathline of airflow than gas molecules. Furthermore, the diffusion mechanism is important for the gas, but diffusion for bioaerosol particles is rather negligible (Pantelic et al., 2009a; Wan, 2006). 62

2.5.2 Discrete phase model (DPM) for bioaerosol transport The discrete phase model (DPM) is a gas-solid two-phase transport model that simulates the continuum (i.e. air) phase in the Eulerian framework and the dispersed (i.e. bioaerosol particles) phase in the Lagrangian framework (Hutchinson et al., 1971). The trajectory of a bioaerosol particle is described by the balance of acting forces on the particle as point-volume descriptions for interphase transfer of momentum in Equation (2.21) (Loth, 2000). 𝐹𝑝 = 𝐹𝑑𝑟𝑎𝑔 + 𝐹𝑔𝑟𝑎𝑣 + 𝐹𝐹𝐾 + 𝐹𝐴𝑀 + 𝐹𝐵𝐻 + 𝐹𝑎𝑑𝑑

[2.21]

Where Fp is the force acted on a particle (N), Fdrag is the drag force (N), Fgrav is the gravity force, FFK is the Froude-Krylov force due to the fluid stress gradients arising from the continuous-phase acceleration (N), FAM is the added mass (N), FBH is the Basset history terms (N) and Fadd is the additional forces acting on the particle (N).

The Maxey and Riley equation was used to simplify the motion of a small spherical shaped particle in unsteady, non-uniform flow at the low Reynolds numbers (i.e. Rep