Modeling Indoor Air Pollution

Transcription

1 Modeling Indoor Air Pollution

2 This page intentionally left blank

3 Modeling Indoor Air Pollution Darrell W Pepper University of Nevada, Las Vegas, USA David Carrington Los Alamos National Laboratory, USA ICP Imperial College Press

4 Published by Imperial College Press 57 Shelton Street Covent Garden London WC2H 9HE Distributed by World Scientific Publishing Co. Pte. Ltd. 5 Toh Tuck Link, Singapore USA office: 27 Warren Street, Suite , Hackensack, NJ UK office: 57 Shelton Street, Covent Garden, London WC2H 9HE British Library Cataloguing-in-Publication Data A catalogue record for this book is available from the British Library. MODELING INDOOR AIR POLLUTION Copyright 2009 by Imperial College Press All rights reserved. This book, or parts thereof, may not be reproduced in any form or by any means, electronic or mechanical, including photocopying, recording or any information storage and retrieval system now known or to be invented, without written permission from the Publisher. For photocopying of material in this volume, please pay a copying fee through the Copyright Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to photocopy is not required from the publisher. ISBN ISBN X Printed in Singapore.

5 To our families, friends, and colleagues

6 This page intentionally left blank

7 Acknowledgments This book is the result of many years of developing and applying various models to simulate the dispersion of contaminants. Some of the information presented stems from numerous interactions and collaborations with students and colleagues. Much of the material presented in this book can be found in other texts. Some of these texts briefly describe a few of the more simple numerical techniques, but generally don t provide source codes or describe the latest numerical techniques. We have attempted to integrate fundamental information and the basic formulations regarding indoor air quality and ventilation from the more popular textbooks and monographs. We wish to especially acknowledge Dr. Xiuling Wang, who diligently converted many of our old FORTRAN codes into MATLAB files, and also developed the COMSOL example files. Also we thank Ms. Kathryn Nelson who developed the website for the book and indoor air quality computer codes. We are grateful to Mrs. Jeannie Pepper, who typed the manuscript and made sure that we followed proper procedures and format. Her professionalism and patience in typing and proofreading the chapters are greatly appreciated. Finally, we wish to acknowledge the assistance and support of Miss Elizabeth Bennett, our Editor with ICP. This book was written using Microsoft Word and MathType, with conversion to Adobe pdf. The figures were made using TECPLOT, Microsoft Vizio, and PhotoShop. vii

8 This page intentionally left blank

9 Preface Indoor air pollution is becoming a serious problem. This is especially pertinent in situations where sick building syndrome, under design air flow, and contaminant dispersion occur. Recent interest in homeland security and the aftermath of terrorist activities have created a desire among many governmental agencies to more fully understand interior pollutant dispersion and risk assessment. It is the intention of this book to acquaint the reader with enough information to begin using various modeling tools for assessing indoor air pollution. There are many levels of models, ranging in sophistication from simple analytical expressions to elegant, 3-D schemes for solving the Navier Stokes equations for fluid flow and species transport. The level of modeling effort resides ultimately with the user, and the desired level of accuracy. While 3-D numerical schemes based on finite difference, finite volume, or finite element techniques provide elegant solutions, they also require a great deal of understanding, patience, and computational resources. Analytical solutions, while fast and simple, may be orders of magnitude off in comparison to actual values. This book presents these most common of numerical and analytical tools that can be used for modeling indoor air pollution and the types of problems where a particular model is best suited. Chapter 1 presents an overview of indoor air pollution, types of ventilation systems, exposure, and general modeling techniques. Chapter 2 discusses the governing mathematical equations that serve as the basis for modeling air pollutant and flow patterns. In Chapter 3 a general discussion of contaminant sources routinely associated with indoor air quality assessment studies is given, with a presentation on particulates and evaporation of droplets. Assessment criteria are described in Chapter 4, including what to consider in exposure levels as well as economical issues associated with design. Chapter 5 introduces the fundamental analytical tools, along with advection and the classic box model approach, for performing simple model simulation, including their ix

10 x Modeling Indoor Air Pollution limitations. The dynamics of particle motion, including particle drag and flow in inlets and flanges, are given in Chapter 6. Chapter 7 describes the fundamental numerical approaches commonly used in CFD-type simulations, which are based on finite difference, finite volume, and finite element techniques. Additional discussions are given in Chapter 8 on more advanced methods that include boundary element, particle-incell, and meshless methods, which is relatively new. Several modeling examples are included. In Chapter 9, an extensive description of turbulence modeling is presented with consideration to both finite volume and finite element techniques. A time-dependent, two-equation closure model is presented in detail using the finite element method with adaptive meshing; results are shown for example problems. Issues regarding homeland security and the potential threats attributed to terrorist activities are discussed in Chapter 10, including an example scenario. The examples and computer techniques discussed in the book are available on the web. The website is: We have elected to write the majority of the codes in MATLAB. The website lists locations where you can also find FORTRAN and C/C++ versions of some of the example codes. We have found that most engineering graduates today as well as science and engineering students are familiar with MATLAB, and prefer using it as their primary coding environment. We have also used COMSOL to run the example problems. COMSOL, with headquarters in Sweden, is a very versatile multiphysics finite element package used throughout the world; the package permits easy interface and flexibility in setting up problems, along with MATLAB scripting. Many universities and companies are now using COMSOL. In addition, an adaptive finite-element based model that can be used for indoor air pollution is also available from the authors. This model utilizes h-adaptation (mesh refinement) to accurately simulate the dispersion of contaminant within any shaped interior. Darrell W. Pepper David B. Carrington 2008

11 Contents Acknowledgements Preface vii ix 1. Introduction What is Indoor Air Pollution Ventilation Systems Exposure Risks Numerical Modeling of Indoor Air Flow Comments Fluid Flow Fundamentals Conservation Equations Ideal Fluids Conformal mapping Schwarz Christoffel transform Numerical mapping Superposition for stream functions Turbulence Species Transport Comments Contaminant Sources Types of Contaminants Units Materials Typical Operations The Diffusion Equation Diffusion in Air Evaporation of Droplets Resuspension of Particulate Coagulation of Particulate Comments Assessment Criteria Exposure Economics Comments xi

12 xii Modeling Indoor Air Pollution 5. Simple Modeling Techniques Analytical Tools Advection Model Box Model Comments Dynamics of Particles, Gases and Vapors Drag, Shape, and Size of Particles Particle Motion Deposition of particulate with aerodynamic diameters > 1μ by settling Particle motion in electrostatic field Particle motion induced by temperature gradients Thermophoretic motion for gases and particles with diameter less than the molecular mean free path Thermophoretic transport for particles with diameter greater than the molecular mean free path Particle Flow in Inlets and Flanges Comments Numerical Modeling Conventional Techniques Finite Difference Method Explicit Implicit Upwinding Finite Volume Method FDM FVM The Finite Element Method One-dimensional elements Linear element Quadratic and higher order elements Two-dimensional elements Triangular elements Quadrilateral elements Isoparametric elements Three-dimensional elements Quadrature Time dependence Petrov Galerkin method Mesh generation Bandwidth Adaptation

13 Contents xiii Element subdivision Further CFD Examples Model Verification and Validation Comments Numerical Modeling Advanced Techniques Boundary Element Method Lagrangian Particle Technique Particle-in-cell Meshless Method Application of meshless methods Smoothed particle hydrodynamics (SPH) techniques including Kernel Particle Methods (RKPM), and general kernel reproduction methods (GKR) Meshless Petrov Galerkin (MLPG) methods including moving least squares (MLS), point interpolation methods (PIM), and hp-clouds Local radial point interpolation methods (LRPIM) using finite difference representations Radial basis functions (RBFs) Example cases Heat Transfer Heat transfer in a 2-D plate Singular point in a 2-D domain Heat transfer within an irregular domain Natural Convection Molecular Modeling Boundary Conditions for Mass Transport Analysis Comments Turbulence Modeling Brief History of Turbulence Formulation Physical Model Turbulent flow Two-equation turbulence closure models Two-equation k-ε Two-equation k-w Large Eddy Simulation (LES) Direct Numerical Simulation (DNS) Turbulent transport of energy or enthalpy Derivation of enthalpy transport Turbulent energy transport Turbulent transport species Coupled fluid-thermal flow

14 xiv Modeling Indoor Air Pollution 9.3 Numerical Modeling Projection algorithm Finite volume approach Finite element approach Weak forms of the governing equations Matrix equations Time advancement of the explicit/implicit matrix equations Mass lumping General numerical solution Stability and Time Dependent Solution Boundary Conditions Boundary conditions for velocity under decomposition Viscous boundary condition for velocity Boundary conditions for pressure and velocity correction Boundary conditions for turbulent kinetic energy and specific dissipation rate Boundary conditions for thermal and species transport Thermal and species flux calculation in the presence of Dirichlet boundaries Validation of Turbulence Models Comments Homeland Security Issues Introduction Potential Hazards Prevention and protection A Simple Model Other Indoor Air Quality Models CONTAM 2.4 (NIST) I-BEAM (EPA) COMIS-MIAQ (APTG-LBNL) FLOVENT (Flomerics, Inc.) Comments Appendix A Diffusion Coefficients in Gas 303 Appendix B 2-D Office Simulations: COMSOL and ANSWER Software 309 Bibliography 323 Index 341