Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The airborne transmission of infection between flats in high-rise residential buildings: Particle simulation</title><author><name sortKey="Gao, N P" sort="Gao, N P" uniqKey="Gao N" first="N. P." last="Gao">N. P. Gao</name><affiliation><nlm:aff id="aff1">Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong</nlm:aff></affiliation></author><author><name sortKey="Niu, J L" sort="Niu, J L" uniqKey="Niu J" first="J. L." last="Niu">J. L. Niu</name><affiliation><nlm:aff id="aff1">Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong</nlm:aff></affiliation></author><author><name sortKey="Perino, M" sort="Perino, M" uniqKey="Perino M" first="M." last="Perino">M. Perino</name><affiliation><nlm:aff id="aff2">DENER, Politecnico di Torino, Corso Duca degli Abruzzi, 24, I-10129 Torino, Italy</nlm:aff></affiliation></author><author><name sortKey="Heiselberg, P" sort="Heiselberg, P" uniqKey="Heiselberg P" first="P." last="Heiselberg">P. Heiselberg</name><affiliation><nlm:aff id="aff3">Hybrid Ventilation Centre, Aalborg University, Sohngaardsholmsvej 57, DK-9000 Aalborg, Denmark</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmc">7116967</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC7116967</idno><idno type="RBID">PMC:7116967</idno><idno type="doi">10.1016/j.buildenv.2008.03.016</idno><idno type="pmid">NONE</idno><date when="2008">2008</date><idno type="wicri:Area/Pmc/Corpus">000E64</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000E64</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The airborne transmission of infection between flats in high-rise residential buildings: Particle simulation</title><author><name sortKey="Gao, N P" sort="Gao, N P" uniqKey="Gao N" first="N. P." last="Gao">N. P. Gao</name><affiliation><nlm:aff id="aff1">Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong</nlm:aff></affiliation></author><author><name sortKey="Niu, J L" sort="Niu, J L" uniqKey="Niu J" first="J. L." last="Niu">J. L. Niu</name><affiliation><nlm:aff id="aff1">Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong</nlm:aff></affiliation></author><author><name sortKey="Perino, M" sort="Perino, M" uniqKey="Perino M" first="M." last="Perino">M. Perino</name><affiliation><nlm:aff id="aff2">DENER, Politecnico di Torino, Corso Duca degli Abruzzi, 24, I-10129 Torino, Italy</nlm:aff></affiliation></author><author><name sortKey="Heiselberg, P" sort="Heiselberg, P" uniqKey="Heiselberg P" first="P." last="Heiselberg">P. Heiselberg</name><affiliation><nlm:aff id="aff3">Hybrid Ventilation Centre, Aalborg University, Sohngaardsholmsvej 57, DK-9000 Aalborg, Denmark</nlm:aff></affiliation></author></analytic><series><title level="j">Building and Environment</title><idno type="ISSN">0360-1323</idno><idno type="eISSN">1873-684X</idno><imprint><date when="2008">2008</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Several case clusters occurred in high-rise residential buildings in Hong Kong in the 2003 SARS (the severe acute respiratory syndrome) epidemic, which motivated a series of engineering investigations into the possible airborne transport routes. It is suspected that, driven by buoyancy force, the polluted air that exits the window of the lower floor may re-enter the immediate upper floor through the window on the same side. This cascade effect has been quantified and reported in a previous paper, and it is found that, by tracer gas concentration analysis, the room in the adjacent upstairs may contain up to 7% of the air directly from the downstairs room. In this study, after validation against the experimental data from literatures, Eulerian and Lagrangian approaches are both adopted to numerically investigate the dispersion of expiratory aerosols between two vertically adjacent flats. It is found that the particle concentration in the upper floor is two to three orders of magnitude lower than in the source floor. 1.0 μm particles disperse like gaseous pollutants. For coarse particles larger than 20.0 μm, strong deposition on solid surfaces and gravitational settling effect greatly limit their upward transport.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Yu, I T S" uniqKey="Yu I">I.T.S. Yu</name></author><author><name sortKey="Li, Y G" uniqKey="Li Y">Y.G. Li</name></author><author><name sortKey="Wong, T W" uniqKey="Wong T">T.W. Wong</name></author><author><name sortKey="Tam, W" uniqKey="Tam W">W. Tam</name></author><author><name sortKey="Chan, A" uniqKey="Chan A">A. Chan</name></author><author><name sortKey="Lee, J H W" uniqKey="Lee J">J.H.W. Lee</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, Y" uniqKey="Li Y">Y. Li</name></author><author><name sortKey="Huang, X" uniqKey="Huang X">X. Huang</name></author><author><name sortKey="Yu, I T S" uniqKey="Yu I">I.T.S. Yu</name></author><author><name sortKey="Wong, T W" uniqKey="Wong T">T.W. Wong</name></author><author><name sortKey="Qian, H" uniqKey="Qian H">H. Qian</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, Y" uniqKey="Li Y">Y. Li</name></author><author><name sortKey="Duan, S" uniqKey="Duan S">S. Duan</name></author><author><name sortKey="Yu, I T S" uniqKey="Yu I">I.T.S. Yu</name></author><author><name sortKey="Wong, T W" uniqKey="Wong T">T.W. Wong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Niu, J L" uniqKey="Niu J">J.L. Niu</name></author><author><name sortKey="Tung, C W" uniqKey="Tung C">C.W. Tung</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Brankston, G" uniqKey="Brankston G">G. Brankston</name></author><author><name sortKey="Gitterman, L" uniqKey="Gitterman L">L. Gitterman</name></author><author><name sortKey="Hirji, Z" uniqKey="Hirji Z">Z. Hirji</name></author><author><name sortKey="Lemieux, C" uniqKey="Lemieux C">C. Lemieux</name></author><author><name sortKey="Gardam, M" uniqKey="Gardam M">M. Gardam</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tellier, R" uniqKey="Tellier R">R. Tellier</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Yakhot, V" uniqKey="Yakhot V">V. Yakhot</name></author><author><name sortKey="Orszag, S" uniqKey="Orszag S">S. Orszag</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Elghobashi, E" uniqKey="Elghobashi E">E. Elghobashi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hinds, W C" uniqKey="Hinds W">W.C. Hinds</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Abadie, M" uniqKey="Abadie M">M. Abadie</name></author><author><name sortKey="Limam, K" uniqKey="Limam K">K. Limam</name></author><author><name sortKey="Allard, F" uniqKey="Allard F">F. Allard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Morawska, L" uniqKey="Morawska L">L. Morawska</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nicas, M" uniqKey="Nicas M">M. Nicas</name></author><author><name sortKey="Nazaroff, W W" uniqKey="Nazaroff W">W.W. Nazaroff</name></author><author><name sortKey="Hubbard, A" uniqKey="Hubbard A">A. Hubbard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Loth, E" uniqKey="Loth E">E. Loth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z. Zhang</name></author><author><name sortKey="Chen, Q" uniqKey="Chen Q">Q. Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, B" uniqKey="Zhao B">B. Zhao</name></author><author><name sortKey="Yang, C" uniqKey="Yang C">C. Yang</name></author><author><name sortKey="Yang, X" uniqKey="Yang X">X. Yang</name></author><author><name sortKey="Liu, S" uniqKey="Liu S">S. Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chao, C Y H" uniqKey="Chao C">C.Y.H. Chao</name></author><author><name sortKey="Wan, M P" uniqKey="Wan M">M.P. Wan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gao, Z" uniqKey="Gao Z">Z. Gao</name></author><author><name sortKey="Mashayek, F" uniqKey="Mashayek F">F. Mashayek</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tian, L" uniqKey="Tian L">L. Tian</name></author><author><name sortKey="Ahmadi, G" uniqKey="Ahmadi G">G. Ahmadi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z. Zhang</name></author><author><name sortKey="Chen, Q" uniqKey="Chen Q">Q. Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Q" uniqKey="Wang Q">Q. Wang</name></author><author><name sortKey="Squires, K D" uniqKey="Squires K">K.D. Squires</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Y" uniqKey="Wang Y">Y. Wang</name></author><author><name sortKey="James, P W" uniqKey="James P">P.W. James</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hinze, J O" uniqKey="Hinze J">J.O. Hinze</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lai, A C K" uniqKey="Lai A">A.C.K. Lai</name></author><author><name sortKey="Nazaroff, W W" uniqKey="Nazaroff W">W.W. Nazaroff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lu, W Z" uniqKey="Lu W">W.Z. Lu</name></author><author><name sortKey="Howarth, A T" uniqKey="Howarth A">A.T. Howarth</name></author><author><name sortKey="Adam, N" uniqKey="Adam N">N. Adam</name></author><author><name sortKey="Riffat, S B" uniqKey="Riffat S">S.B. Riffat</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, F Z" uniqKey="Chen F">F.Z. Chen</name></author><author><name sortKey="Yu, S C M" uniqKey="Yu S">S.C.M. Yu</name></author><author><name sortKey="Lai, A C K" uniqKey="Lai A">A.C.K. Lai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y. Liu</name></author><author><name sortKey="Matida, E A" uniqKey="Matida E">E.A. Matida</name></author><author><name sortKey="Gu, J" uniqKey="Gu J">J. Gu</name></author><author><name sortKey="Johnson, M R" uniqKey="Johnson M">M.R. Johnson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Matida, E A" uniqKey="Matida E">E.A. Matida</name></author><author><name sortKey="Finlay, W H" uniqKey="Finlay W">W.H. Finlay</name></author><author><name sortKey="Lange, C F" uniqKey="Lange C">C.F. Lange</name></author><author><name sortKey="Grgic, B" uniqKey="Grgic B">B. Grgic</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lai, A C K" uniqKey="Lai A">A.C.K. Lai</name></author><author><name sortKey="Chen, F" uniqKey="Chen F">F. Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Murakami, S" uniqKey="Murakami S">S. Murakami</name></author><author><name sortKey="Kato, S" uniqKey="Kato S">S. Kato</name></author><author><name sortKey="Nagano, S" uniqKey="Nagano S">S. Nagano</name></author><author><name sortKey="Tanaka, Y" uniqKey="Tanaka Y">Y. Tanaka</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, B" uniqKey="Zhao B">B. Zhao</name></author><author><name sortKey="Guan, P" uniqKey="Guan P">P. Guan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gao, N P" uniqKey="Gao N">N.P. Gao</name></author><author><name sortKey="Niu, J L" uniqKey="Niu J">J.L. Niu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Uijttewaal, W S J" uniqKey="Uijttewaal W">W.S.J. Uijttewaal</name></author><author><name sortKey="Oliemans, R V A" uniqKey="Oliemans R">R.V.A. Oliemans</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Beghein, C" uniqKey="Beghein C">C. Béghein</name></author><author><name sortKey="Jiang, Y" uniqKey="Jiang Y">Y. Jiang</name></author><author><name sortKey="Chen, Q" uniqKey="Chen Q">Q. Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, X X" uniqKey="Li X">X.X. Li</name></author><author><name sortKey="Liu, C H" uniqKey="Liu C">C.H. Liu</name></author><author><name sortKey="Leung, D Y C" uniqKey="Leung D">D.Y.C. Leung</name></author><author><name sortKey="Lam, K M" uniqKey="Lam K">K.M. Lam</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duguid, J P" uniqKey="Duguid J">J.P. Duguid</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Build Environ</journal-id><journal-id journal-id-type="iso-abbrev">Build Environ</journal-id><journal-title-group><journal-title>Building and Environment</journal-title></journal-title-group><issn pub-type="ppub">0360-1323</issn><issn pub-type="epub">1873-684X</issn><publisher><publisher-name>Elsevier Ltd.</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmc">7116967</article-id><article-id pub-id-type="publisher-id">S0360-1323(08)00053-X</article-id><article-id pub-id-type="doi">10.1016/j.buildenv.2008.03.016</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>The airborne transmission of infection between flats in high-rise residential buildings: Particle simulation</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Gao</surname><given-names>N.P.</given-names></name><xref rid="aff1" ref-type="aff">a</xref></contrib><contrib contrib-type="author"><name><surname>Niu</surname><given-names>J.L.</given-names></name><email>Bejlniu@polyu.edu.hk</email><xref rid="aff1" ref-type="aff">a</xref><xref rid="cor1" ref-type="corresp">⁎</xref></contrib><contrib contrib-type="author"><name><surname>Perino</surname><given-names>M.</given-names></name><xref rid="aff2" ref-type="aff">b</xref></contrib><contrib contrib-type="author"><name><surname>Heiselberg</surname><given-names>P.</given-names></name><xref rid="aff3" ref-type="aff">c</xref></contrib></contrib-group><aff id="aff1"><label>a</label>Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong</aff><aff id="aff2"><label>b</label>DENER, Politecnico di Torino, Corso Duca degli Abruzzi, 24, I-10129 Torino, Italy</aff><aff id="aff3"><label>c</label>Hybrid Ventilation Centre, Aalborg University, Sohngaardsholmsvej 57, DK-9000 Aalborg, Denmark</aff><author-notes><corresp id="cor1"><label>⁎</label>Corresponding author. Tel.: +852 2766 7781; fax: +852 2774 6146. <email>Bejlniu@polyu.edu.hk</email></corresp></author-notes><pub-date pub-type="pmc-release"><day>6</day><month>4</month><year>2008</year></pub-date><pmc-comment> PMC Release delay is 0 months and 0 days and was based on .</pmc-comment>
<pub-date pub-type="ppub"><month>2</month><year>2009</year></pub-date><pub-date pub-type="epub"><day>6</day><month>4</month><year>2008</year></pub-date><volume>44</volume><issue>2</issue><fpage>402</fpage><lpage>410</lpage><history><date date-type="received"><day>11</day><month>1</month><year>2008</year></date><date date-type="rev-recd"><day>26</day><month>3</month><year>2008</year></date><date date-type="accepted"><day>27</day><month>3</month><year>2008</year></date></history><permissions><copyright-statement>Copyright © 2008 Elsevier Ltd. All rights reserved.</copyright-statement><copyright-year>2008</copyright-year><copyright-holder>Elsevier Ltd</copyright-holder><license><license-p>Since January 2020 Elsevier has created a COVID-19 resource centre with free information in English and Mandarin on the novel coronavirus COVID-19. The COVID-19 resource centre is hosted on Elsevier Connect, the company's public news and information website. Elsevier hereby grants permission to make all its COVID-19-related research that is available on the COVID-19 resource centre - including this research content - immediately available in PubMed Central and other publicly funded repositories, such as the WHO COVID database with rights for unrestricted research re-use and analyses in any form or by any means with acknowledgement of the original source. These permissions are granted for free by Elsevier for as long as the COVID-19 resource centre remains active.</license-p></license></permissions><abstract><p>Several case clusters occurred in high-rise residential buildings in Hong Kong in the 2003 SARS (the severe acute respiratory syndrome) epidemic, which motivated a series of engineering investigations into the possible airborne transport routes. It is suspected that, driven by buoyancy force, the polluted air that exits the window of the lower floor may re-enter the immediate upper floor through the window on the same side. This cascade effect has been quantified and reported in a previous paper, and it is found that, by tracer gas concentration analysis, the room in the adjacent upstairs may contain up to 7% of the air directly from the downstairs room. In this study, after validation against the experimental data from literatures, Eulerian and Lagrangian approaches are both adopted to numerically investigate the dispersion of expiratory aerosols between two vertically adjacent flats. It is found that the particle concentration in the upper floor is two to three orders of magnitude lower than in the source floor. 1.0 μm particles disperse like gaseous pollutants. For coarse particles larger than 20.0 μm, strong deposition on solid surfaces and gravitational settling effect greatly limit their upward transport.</p></abstract><kwd-group><title>Keywords</title><kwd>Particles</kwd><kwd>Lagrangian method</kwd><kwd>Eulerian method</kwd><kwd>Transport</kwd><kwd>Deposition</kwd></kwd-group></article-meta></front><body><sec><title>Nomenclatures</title><def-list><def-item><term><italic>C</italic></term><def><p>particle concentration (g/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>C</italic><sub>0</sub></term><def><p>initial indoor particle concentration (g/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>C</italic><sub>C</sub></term><def><p>cunningham correction factor</p></def></def-item><def-item><term><italic>C</italic><sub>D</sub></term><def><p>drag coefficient</p></def></def-item><def-item><term><italic>C<sub>j</sub></italic></term><def><p>particle concentration in the <italic>j</italic>th cell (g/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>C<sub>t</sub></italic></term><def><p>indoor particle concentration at time <italic>t</italic> (g/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>d</italic><sub>p</sub></term><def><p>particle diameter (μm)</p></def></def-item><def-item><term><italic>D</italic><sub>p</sub></term><def><p>Brownian diffusivity (m<sup>2</sup>/s)</p></def></def-item><def-item><term><italic>d<sub>t</sub></italic>(<italic>i</italic>, <italic>j</italic>)</term><def><p>the <italic>i</italic>th particle residence time in the <italic>j</italic>th cell (s)</p></def></def-item><def-item><term><inline-formula><mml:math id="M1" altimg="si1.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>correction factor for turbulent kinetic energy</p></def></def-item><def-item><term><inline-formula><mml:math id="M2" altimg="si2.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">addi</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>other additional force per unit mass (N/kg)</p></def></def-item><def-item><term><inline-formula><mml:math id="M3" altimg="si3.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">drag</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>drag force per unit mass (N/kg)</p></def></def-item><def-item><term><inline-formula><mml:math id="M4" altimg="si4.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">grav</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>gravity force per unit mass (N/kg)</p></def></def-item><def-item><term><inline-formula><mml:math id="M5" altimg="si5.gif" overflow="scroll"><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:math></inline-formula></term><def><p>gravitational acceleration (m/s<sup>2</sup>)</p></def></def-item><def-item><term><italic>k</italic></term><def><p>turbulent kinetic energy (m<sup>2</sup>/s<sup>2</sup>)</p></def></def-item><def-item><term><italic>L</italic><sub>e</sub></term><def><p>eddy length scale (mm)</p></def></def-item><def-item><term><italic>M</italic></term><def><p>flow rate of each trajectory (g/s)</p></def></def-item><def-item><term><italic>S</italic><sub>C</sub></term><def><p>particle source term (g/sm<sup>3</sup>)</p></def></def-item><def-item><term><italic>S<sub>ϕ</sub></italic></term><def><p>pollutant source term (g/sm<sup>3</sup>)</p></def></def-item><def-item><term><italic>t</italic></term><def><p>time (<italic>t</italic>)</p></def></def-item><def-item><term><italic>t</italic><sub>cross</sub></term><def><p>eddy crossing time (s)</p></def></def-item><def-item><term><italic>t</italic><sub>eddy</sub></term><def><p>eddy lifetime (s)</p></def></def-item><def-item><term><italic>t</italic><sub>L</sub></term><def><p>Lagrangian time scale (s)</p></def></def-item><def-item><term><inline-formula><mml:math id="M6" altimg="si6.gif" overflow="scroll"><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:math></inline-formula></term><def><p>velocity vector (m/s)</p></def></def-item><def-item><term><inline-formula><mml:math id="M7" altimg="si7.gif" overflow="scroll"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula></term><def><p>fluctuation part of instantaneous velocity in <italic>i</italic> direction (m/s)</p></def></def-item><def-item><term><inline-formula><mml:math id="M8" altimg="si8.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>particle velocity (m/s)</p></def></def-item><def-item><term><inline-formula><mml:math id="M9" altimg="si9.gif" overflow="scroll"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula></term><def><p>fluctuation part of instantaneous velocity normal to the wall (m/s)</p></def></def-item><def-item><term><italic>u<sub>τ</sub></italic></term><def><p>wall shear velocity (m/s)</p></def></def-item><def-item><term><italic>V<sub>j</sub></italic></term><def><p>the volume of the <italic>j</italic>th cell (m<sup>3</sup>)</p></def></def-item><def-item><term><inline-formula><mml:math id="M10" altimg="si10.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>gravitational settling velocity calculated from Stokes's law (m/s)</p></def></def-item><def-item><term><italic>V</italic><sub>met</sub></term><def><p>meteorological wind speed (m/s)</p></def></def-item><def-item><term><italic>V<sub>y</sub></italic></term><def><p>wind speed at the height of <italic>y</italic> (m/s)</p></def></def-item><def-item><term><italic>y</italic></term><def><p>height above the ground in the calculation of wind speed (m)</p></def></def-item><def-item><term><italic>y</italic><sub>cell</sub></term><def><p>the distance between the wall and the first cell center (m)</p></def></def-item><def-item><term><italic>y</italic><sup>+</sup></term><def><p>non-dimensional wall distance (<inline-formula><mml:math id="M11" altimg="si11.gif" overflow="scroll"><mml:mrow><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mi>ρ</mml:mi><mml:msub><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>τ</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">cell</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:mi>μ</mml:mi></mml:mrow></mml:math></inline-formula>)</p></def></def-item></def-list><def-list><title>Greek letters</title><def-item><term><italic>α</italic></term><def><p>air change rate (h<sup>−1</sup>)</p></def></def-item><def-item><term><italic>ε</italic></term><def><p>viscous dissipation rate (m<sup>2</sup>/s<sup>3</sup>)</p></def></def-item><def-item><term><italic>ε</italic><sub>p</sub></term><def><p>particle eddy diffusivity (m<sup>2</sup>/s)</p></def></def-item><def-item><term><italic>κ</italic></term><def><p>particle loss-rate coefficient caused by deposition (h<sup>−1</sup>)</p></def></def-item><def-item><term><italic>λ</italic></term><def><p>molecular mean free length (μm)</p></def></def-item><def-item><term><italic>μ</italic></term><def><p>molecular viscosity of the air (g/ms)</p></def></def-item><def-item><term><italic>μ</italic><sub>eff</sub></term><def><p>turbulent effective viscosity (g/ms)</p></def></def-item><def-item><term><inline-formula><mml:math id="M12" altimg="si12.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>ν</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></term><def><p>turbulent viscosity (m<sup>2</sup>/s)</p></def></def-item><def-item><term><italic>ς<sub>i</sub></italic></term><def><p>a random number with Gaussian distribution</p></def></def-item><def-item><term><italic>ρ</italic></term><def><p>air density (kg/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>ρ</italic><sub>p</sub></term><def><p>particle material density (kg/m<sup>3</sup>)</p></def></def-item><def-item><term><italic>σ</italic><sub>c</sub></term><def><p>non-dimensional number 1.0</p></def></def-item><def-item><term><italic>τ</italic><sub>p</sub></term><def><p>particle response time (s)</p></def></def-item></def-list></sec><sec><label>1</label><title>Introduction</title><p>The SARS outbreak in 2003 stimulated a series of engineering investigations into the airborne infectious diseases transmission mechanisms in the built environment. A few studies attempted to combine epidemiologic investigations with airflow analysis, but the airflow analysis remained to be preliminary <xref rid="bib1" ref-type="bibr">[1]</xref>, <xref rid="bib2" ref-type="bibr">[2]</xref>. In the investigation of the outbreak in a high-rise residential estate in Hong Kong, the CFD simulation did indicate that the transmission pattern between adjacent high-rise blocks qualitatively agrees with the wind direction. In the related study aiming to reveal the vertical transmission pattern in one building, multi-zone modeling technique was employed, which had inadequacies in modeling airflow through large window openings <xref rid="bib3" ref-type="bibr">[3]</xref>. Niu and Tung <xref rid="bib4" ref-type="bibr">[4]</xref> used on-site tracer gas measurement to experimentally investigate the re-entry possibility of exhaust air from the lower floor into the immediate upper floor for high-rise residential buildings with single-sided natural ventilation conditions (<xref rid="fig1" ref-type="fig">Fig. 1</xref>). It was revealed that the upstairs room may contain 7% of the exhaust air from the lower floor, which can help explain the presence of SARS-CoV RNA found on the window sill deposits in the upstairs of an index patient. Gao et al.'s <xref rid="bib5" ref-type="bibr">[5]</xref> CFD study on the dispersion of one tracer gas agrees with the on-site measurement. With the help of the Wells–Riley infection risk model, Gao et al. <xref rid="bib5" ref-type="bibr">[5]</xref> estimated the vertical cross-household infection risk in the upper floor using the data of tuberculosis. It was found that a significant infection risk up to 6% for the upstairs residents can be resulted in, when merely assuming an 8 h continuous stay at home.<fig id="fig1"><label>Fig. 1</label><caption><p>Schematic view of pollutant transport by the cascade effect in natural ventilation.</p></caption><graphic xlink:href="gr1"></graphic></fig></p><p>On the other hand, recent systematic reviews <xref rid="bib6" ref-type="bibr">[6]</xref>, <xref rid="bib7" ref-type="bibr">[7]</xref> concluded that person to person respiratory virus transmission could occur in three possible modes: airborne transmission, droplet transmission, and direct contact with secretions (or fomites). Currently, by droplets transmission it is meant short-distance (usually less than 1 m) transmission via large droplets (⩾5 μm diameter) generated during coughing, sneezing, or talking, and by airborne transmission it is meant long-distance transmission via the dissemination of virus-laden droplet nuclei (particles less than 5 μm that result from the evaporation of large droplets). But whether the three transmission modes are mutually exclusive and which mode is the most significant transmission route for a particular disease are debatable. In general, it is believed that respiratory droplets movements in the room air play a key role for airborne transmitted diseases. Arguably, airborne transmission does not necessarily cause long range infections since a wide range of other confounding factors are also playing important roles, and these include the shedding rate, viability and infectiousness of the pathogen and the wind, temperature and humidity conditions. In general, aerosol dynamics differs from gaseous pollutants due to the effects of gravity, inertia and deposition at solid surfaces. Therefore, as a companion paper of Gao et al. <xref rid="bib5" ref-type="bibr">[5]</xref>, which studied the airflow and gaseous pollutants between vertically adjacent floors, this study simulates the transport possibilities of droplet residuals by using both Lagrangian and Eulerian methods. The aim is to explore whether the virus-laden aerosols could be carried upward by inter-flat airflows between floors in high-rise residential buildings when the gravity and surface deposition effects are considered.</p></sec><sec><label>2</label><title>Numerical approaches</title><sec><label>2.1</label><title>Airflow field</title><p>Prior to modeling particle movements, the airflow field is simulated by using a commercial code, Fluent <xref rid="bib8" ref-type="bibr">[8]</xref>, which solves the mass, momentum, and energy equations in a finite-volume procedure with a staggered grid system. The convection terms are discretized by second order upwind scheme and the diffusion term by central difference also with second-order accuracy. The variables at the near-wall cells and the corresponding quantities on the wall are bridged by the standard logarithmic law wall function. The turbulent effect is modeled by the RNG k–ε model <xref rid="bib9" ref-type="bibr">[9]</xref>. Although large eddy simulation (LES) has been adopted to investigate wind field and indoor–outdoor airflows due to its capacity of handling the unsteadiness and intermittency of the flow as well as providing detailed information on the turbulence structure, the RANS models are extensively used in industrial and engineering practice because they generally use about one order of magnitude less CPU time than LES while keeping an acceptable accuracy.</p></sec><sec><label>2.2</label><title>Assumptions for particle modeling</title><p>When particle loadings, denoted by volume fraction, are higher than 10<sup>−6</sup>, particles significantly affect the turbulent flow and two-way coupling between the phases has to be taken into account. If particle loadings are even higher than 10<sup>−3</sup>, particle–particle interaction becomes important <xref rid="bib10" ref-type="bibr">[10]</xref>. Since indoor particulate levels normally vary from several to hundreds of μg/m<sup>3</sup> (about 10<sup>5</sup>–10<sup>7</sup> particles/m<sup>3</sup>) the two-phase flow can be deemed as dilute systems and one-way coupling is rational. A single-phase turbulence model can be used to represent the mixture.</p><p>In both numerical approaches, all particles are assumed to be spherical. Particle coagulation, reflection at walls, re-suspension, and phase-change such as evaporation, which may actually occur in indoor environments, are ignored. According to Hinds <xref rid="bib11" ref-type="bibr">[11]</xref>, coagulation rate is proportional to the square of the number concentration. Supposing a diameter of 1.0 μm and a concentration of 10<sup>7</sup> particles per m<sup>3</sup>, 1% decrement of particle number concentration needs about 500 days. This time scale is much larger than the indoor airflow time scale in that the indoor air can be renewed several times per hour. Particles’ bounce back from a solid surface is influenced by particle size, surface material, orientation, initial momentum and even electrostatic forces. Fine particles have a strong probability to stick on the solid surfaces, especially on those made of synthetic fibers <xref rid="bib12" ref-type="bibr">[12]</xref>. Morawska <xref rid="bib13" ref-type="bibr">[13]</xref> demonstrated that pure water droplets may completely evaporate within 0.003 and 0.3 s for 1 and 10 μm droplets respectively, at a condition of 80% relative humidity (RH). These periods are much shorter than their settling time, the indoor airflow time scale, and even shorter than the sneezing and coughing processes (around 0.5 s). Nicas et al. <xref rid="bib14" ref-type="bibr">[14]</xref> suggested that evaporation from respiratory particles with an initial diameter less than 20 μm be treated as an instantaneous process. The droplets produced from sneezing and coughing could have a significant amount of pathogens-carrying residues. The equilibrium diameter of the desiccated particles is generally estimated to be around one-half of its initial diameter <xref rid="bib14" ref-type="bibr">[14]</xref>. Small pathogens-carrying particles (diameter less than 10 μm) can reach the alveolar region with a different efficiency, posing an infection risk for a susceptible person <xref rid="bib11" ref-type="bibr">[11]</xref>. Based on this information, in this study we simulated the movements of particles with diameter of 1, 10, and 20 μm.</p></sec><sec><label>2.3</label><title>Lagrangian method—discrete random walk model</title><p>Mathematical treatments of indoor particle movements include: the Euler–Lagrange method and the Euler–Euler method. Each method has its own pros and cons <xref rid="bib15" ref-type="bibr">[15]</xref>, <xref rid="bib16" ref-type="bibr">[16]</xref>.</p><p>The Lagrangian approach considers the fluid phase as a continuum and predicts the trajectory of each discrete phase particle by integrating the force balance on the particle. The particle momentum equation can be expressed as:<disp-formula id="eq1"><label>(1)</label><mml:math id="M13" altimg="si13.gif" overflow="scroll"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi></mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">drag</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">grav</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">addi</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mtext>,</mml:mtext></mml:mrow></mml:math></disp-formula><inline-formula><mml:math id="M14" altimg="si14.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">drag</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M15" altimg="si15.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">grav</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> can be written as <inline-formula><mml:math id="M16" altimg="si16.gif" overflow="scroll"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>18</mml:mn><mml:mi>μ</mml:mi><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">C</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M17" altimg="si17.gif" overflow="scroll"><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi>ρ</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, respectively, where the Cunningham correction to Stokes’ drag law is computed from <inline-formula><mml:math id="M18" altimg="si18.gif" overflow="scroll"><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>2</mml:mn><mml:mi>λ</mml:mi><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>1.257</mml:mn><mml:mo>+</mml:mo><mml:mn>0.4</mml:mn><mml:mspace width=".16em"></mml:mspace><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="normal">e</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>1.1</mml:mn><mml:msub><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:mn>2</mml:mn><mml:mi>λ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>. The third term in Eq. <xref rid="eq1" ref-type="disp-formula">(1)</xref> represents the additional forces exerted on the particles. In the current study, Brownian force, and Saffman's Lift force are taken into account in that they may be significant for fine particles. Thermophoretic force is included as well to reflect the non-isothermal condition. As the density ratio between particles and air is in the order of O(10<sup>3</sup>), other forces, such as Basset force, pressure gradient force, virtual mass force, are neglected because their fraction in <inline-formula><mml:math id="M19" altimg="si19.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">addi</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is infinitesimal <xref rid="bib17" ref-type="bibr">[17]</xref>. Detailed descriptions of these forces are available in the Fluent manual <xref rid="bib8" ref-type="bibr">[8]</xref>. The effect of turbulence on particle dispersion is modeled by using a stochastic method, i.e. the discrete random walk (DRW) model to simulate the fluctuating velocity. Stochastic realizations of turbulence have been successfully applied in many works <xref rid="bib18" ref-type="bibr">[18]</xref>, <xref rid="bib19" ref-type="bibr">[19]</xref>. The DRW model assumes turbulence is isotropic and the random velocity fluctuations follow a Gaussian probability distribution:<disp-formula id="eq2"><label>(2)</label><mml:math id="M20" altimg="si20.gif" overflow="scroll"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>ς</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">¯</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:msqrt><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>ς</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>The fluctuating velocities remain constant over a time interval that is the minimum between two time scales: the characteristic lifetime of the eddy and the time required by the particle to move across the eddy. This time interval is also called interaction time between the particle and eddies. The eddy life time is defined as:<disp-formula id="eq3"><label>(3)</label><mml:math id="M21" altimg="si21.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">eddy</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">L</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.30</mml:mn><mml:mfrac><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>ε</mml:mi></mml:mrow></mml:mfrac><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>The crossing time is calculated by<disp-formula id="eq4"><label>(4)</label><mml:math id="M22" altimg="si22.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">cross</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mspace width=".16em"></mml:mspace><mml:mi>ln</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">e</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mtext>,</mml:mtext></mml:mrow></mml:math></disp-formula><disp-formula id="eq5"><label>(5)</label><mml:math id="M23" altimg="si23.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">C</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn>18</mml:mn><mml:mi>μ</mml:mi></mml:mrow></mml:mfrac><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>In a new interaction time, the particle is at a new location, and the calculation of the particle velocity and trajectory is continued using the new local instantaneous air velocities.</p><p>The ordinary differential Eq. <xref rid="eq1" ref-type="disp-formula">(1)</xref> is integrated with the fifth-order Runge–Kutta method. The particle motion is advanced with a time step which is one order of magnitude smaller than the crossing time. For the treatment of particle deposition at solid surfaces in Lagrangian method, particles reaching a wall surface are assumed to be trapped, meaning the calculation of trajectories is terminated.</p><p>However, in DRW model isotropic turbulence assumption is not applicable at near-wall regions. In the boundary layer, velocity fluctuations in the direction normal to the wall are substantially smaller than in the spanwise and streamwise directions. Thus, isotropic decomposition of kinetic energy over-predicts the deposition rate of small particles <xref rid="bib20" ref-type="bibr">[20]</xref>, <xref rid="bib21" ref-type="bibr">[21]</xref>. In this scenario, Reynolds stress model or LES are much more preferable <xref rid="bib22" ref-type="bibr">[22]</xref>. If the flow field is simulated in the framework of two-equation models, to settle this problem Wang and James <xref rid="bib23" ref-type="bibr">[23]</xref> developed the following equations to modify the turbulent kinetic energy for the first cells close to the wall:<disp-formula id="eq6"><label>(6)</label><mml:math id="M24" altimg="si24.gif" overflow="scroll"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>ς</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt><mml:mo>×</mml:mo><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mtext>,</mml:mtext></mml:mrow></mml:math></disp-formula><disp-formula id="eq7"><label>(7)</label><mml:math id="M25" altimg="si25.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mi>exp</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mo>-</mml:mo><mml:mn>0.02</mml:mn><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula>With the help of Eqs. <xref rid="eq6" ref-type="disp-formula">(6)</xref>, <xref rid="eq7" ref-type="disp-formula">(7)</xref>, the near-wall kinetic energy is artificially reduced, prior to tracking the particles’ trajectories. Its improvement on the prediction of deposition rate will be examined in the validation section.</p></sec><sec><label>2.4</label><title>Eulerian method—the drift-flux model</title><p>In the Eulerian approach, both the fluid phase and the particulate phase are treated as interpenetrating continua. In order to save computational load, a simplified Eulerian method, i.e. the drift-flux method, is used. In this model, a drift velocity consisting of gravitational setting and diffusion is taken into account between particles and bulk air. The governing equation for particle transport in turbulent flow is written as:<disp-formula id="eq8"><label>(8)</label><mml:math id="M26" altimg="si26.gif" overflow="scroll"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mo mathvariant="normal">∂</mml:mo></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>ρ</mml:mi><mml:mi>C</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mo mathvariant="normal">∂</mml:mo></mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mo>∇</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>ρ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mi>C</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mo>∇</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mi>ρ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mi>ε</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>∇</mml:mo><mml:mi>C</mml:mi><mml:mo stretchy="false">]</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">C</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>The gravitational settling velocity of particles (<inline-formula><mml:math id="M27" altimg="si27.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>) is calculated by Stokes equation <xref rid="bib11" ref-type="bibr">[11]</xref>. In view of that in the present study, the Reynolds number for particles is usually less than 1.0, <inline-formula><mml:math id="M28" altimg="si28.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo>⇒</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> equals <inline-formula><mml:math id="M29" altimg="si29.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">c</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mo>¯</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo stretchy="false">/</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>18</mml:mn><mml:mi>μ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>. For particles with very small relaxation time, particle eddy diffusivity <inline-formula><mml:math id="M30" altimg="si30.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>ε</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> approximates the carried fluid turbulent diffusivity <inline-formula><mml:math id="M31" altimg="si31.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>ν</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula><xref rid="bib24" ref-type="bibr">[24]</xref>. Generally speaking, for particles larger than 1.0 μm, Bownian diffusion is much weaker than the turbulent diffusion and gravitational settling effect is significant. Therefore, in the implementation of Eq. <xref rid="eq8" ref-type="disp-formula">(8)</xref>, we replace <inline-formula><mml:math id="M32" altimg="si32.gif" overflow="scroll"><mml:mrow><mml:mi>ρ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mi>ε</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula> by <italic>μ</italic><sub>eff</sub>. Eq. <xref rid="eq8" ref-type="disp-formula">(8)</xref> is discretized directly into an algebraic equation by the finite volume method, which is unlike the other numerical treatment where the settling term <inline-formula><mml:math id="M33" altimg="si33.gif" overflow="scroll"><mml:mrow><mml:mi>ρ</mml:mi><mml:msub><mml:mrow><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo>¯</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mi>C</mml:mi></mml:mrow></mml:math></inline-formula> is moved into the source term. For the boundary condition at walls, a semi-empirical deposition model by Lai and Nazaroff <xref rid="bib25" ref-type="bibr">[25]</xref> is adopted, which is able to link the local airflow condition with the deposition process, by expressing the deposition velocity as a function of particle density, size, and friction velocity.</p></sec></sec><sec><label>3</label><title>Validations</title><p>Lu et al. <xref rid="bib26" ref-type="bibr">[26]</xref> measured particles decay in a two-room chamber with a size of width×depth×height=5 m×3 m×2.4 m (<xref rid="fig2" ref-type="fig">Fig. 2</xref>). The chamber is divided into two rooms by a partition in the middle. There is a large opening with height×width=0.95 m×0.70 m between the rooms. The clean air is supplied into one room through the inlet (1.0 m×0.5 m) and exhausted from the other room through the outlet (1.0 m×0.5 m) at the low level. The air change rate equals 9.216 h<sup>−1</sup>. In the experiments, the oil–smoke particles with a material density of 865.0 kg/m<sup>3</sup> are equally divided into five size groups (1, 2, 3, 4, and 5 μm). Once the released particles are uniformly distributed in room 1 without initial motion, dampers at both inlet and outlet diffusers are opened and the supply fan is switched on. The particle mass concentrations at the centers of the two rooms are measured to represent the average concentrations in each room.<fig id="fig2"><label>Fig. 2</label><caption><p>Geometry of the two rooms connected by an internal door.</p></caption><graphic xlink:href="gr2"></graphic></fig></p><p>To numerically simulate this case, each zone is divided into 84,000 cells and 8000 sample particles (1600 particles in each size group) are tracked. The particle mass carried by each sample particle is the ratio of the total particle mass injected into room 1 in the experiment to the samples (8000 particles). By counting the suspended particles, the numerical mass concentration in each room can be determined. The temporal evolutions of particle concentrations from measurements, the DRW model, the DRW model plus near-wall modification, and the drift-flux model are illustrated in <xref rid="fig3" ref-type="fig">Fig. 3</xref>. For comparison, the decay process of tracer gas (CO<sub>2</sub>) is simulated as well. It is observed that the simulation results from tracer gas and the drift-flux model are in agreement with the experiments although the mean concentration in room 2 is slightly under-predicted. Particle loss in rooms is caused by the combined effect of ventilation and deposition. Using a well-mixed model, the mass concentration in room 1 can be calculated by the following equation:<disp-formula id="eq9"><label>(9)</label><mml:math id="M34" altimg="si34.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mspace width=".16em"></mml:mspace><mml:mi>exp</mml:mi><mml:mo stretchy="false">[</mml:mo><mml:mo>-</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mi>κ</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">]</mml:mo><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula><fig id="fig3"><label>Fig. 3</label><caption><p>Temporal evolutions of particle concentrations in rooms 1 and 2.</p></caption><graphic xlink:href="gr3"></graphic></fig></p><p>Because the deposition rate of particles ranging from 1 to 5 μm is much lower than the current air change rate, the particle mass concentration decrement should be similar with the decay of CO<sub>2</sub>. The DRW model greatly under-predicts the concentrations in both rooms by over-predicting the deposition rates. Near-wall modification can improve the prediction, but a large discrepancy with measured data still appears. We believe that this amplified deposition is mainly caused by the anisotropic turbulence in the boundary layer, rather than by the “all trap” assumption at walls.</p><p>Chen et al. <xref rid="bib27" ref-type="bibr">[27]</xref> experimentally investigated the concentration distributions of mono-dispersed 10 μm particles with a material density of 1400 kg/m<sup>3</sup> in a simple model room using a phase Doppler anemometry system (<xref rid="fig4" ref-type="fig">Fig. 4</xref>). Particles are mixed into the supply air with a stable particle flow rate by a solid particle disperser. Particle concentration is normalized by the inlet concentration, which is considered to be uniform. Air inlet and outlet (0.04 m×0.04 m) are symmetrical with the central plane <italic>y</italic>=0.2 m. The supply velocity is 0.225 m/s.<fig id="fig4"><label>Fig. 4</label><caption><p>Geometry of the scaled chamber.</p></caption><graphic xlink:href="gr4"></graphic></fig></p><p><xref rid="fig5" ref-type="fig">Fig. 5</xref>compares the simulated and measured concentration profiles at the center plane. Because the DRW model does not directly provide the particle concentration, the scheme of particle source in-cell (PSI-C) is adopted to correlate the concentration with the trajectories:<disp-formula id="eq10"><label>(10)</label><mml:math id="M35" altimg="si35.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>M</mml:mi><mml:msubsup><mml:mrow><mml:mo>∑</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mi mathvariant="italic">dt</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula><fig id="fig5"><label>Fig. 5</label><caption><p>Comparison of the measured and simulated particle concentrations at the center plane.</p></caption><graphic xlink:href="gr5"></graphic></fig></p><p>Statistically stable concentration requires sufficient particle trajectories. These required trajectories also depend on the number of computational cells. In the present computational domain of 16,000 cells, 40,000 particle tracks are found to be stable. Compared with 20,000 particles, the difference of the proportion of particles with each kind of fate (deposited onto a certain wall or exhausted) to the total released numbers is less than 0.5%. <xref rid="fig5" ref-type="fig">Fig. 5</xref> indicates again that the drift-flux model performs better than the DRW model, and that the DRW model under-predicts the concentration, especially in the region close to the floor. <xref rid="fig6" ref-type="fig">Fig. 6</xref>shows the concentrations in the plane <italic>y</italic>=0.2 m after near-wall kinetic energy modification. In the right part under the ceiling, an unreasonably high concentration appears. It implies that the modification does not improve the simulation in this case. Inherently, the DRW model will show a tendency for fine particles to concentrate in low-turbulence regions of the flow, and it may give non-physical results in strongly non-homogeneous diffusion-dominated flows.<fig id="fig6"><label>Fig. 6</label><caption><p>Normalized particle concentration distributions from Lagrangian simulation with near-wall modification of kinetic energy.</p></caption><graphic xlink:href="gr6"></graphic></fig></p><p>Based on our experiences and a small-scale literature review, the performances of Lagrangian and Eulerian approaches are summarized in <xref rid="tbl1" ref-type="table">Table 1</xref>. In the following simulations, we use the drift-flux model to evaluate the concentration distributions in the adjacent flats and only use the DRW model to visualize the movement of individual particles.<table-wrap position="float" id="tbl1"><label>Table 1</label><caption><p>Comparison of Lagrangian and Eulerian approaches in modeling particle dispersion and deposition</p></caption><table frame="hsides" rules="groups"><thead><tr><th rowspan="2">Method number</th><th colspan="2">Primary phase<hr></hr></th><th colspan="2">Second phase<hr></hr></th><th colspan="2">Performance<hr></hr></th><th rowspan="2">References</th></tr><tr><th>Turbulence model</th><th>Near wall treatment</th><th>Particle modeling</th><th>Boundary condition at walls</th><th>Dispersion</th><th>Deposition</th></tr></thead><tbody><tr><td>1</td><td>k–ε models</td><td>Wall functions</td><td>DRW model</td><td>Particles trapped at walls</td><td>Acceptable</td><td>Greatly over-predicted</td><td><xref rid="bib20" ref-type="bibr">[20]</xref>, <xref rid="bib28" ref-type="bibr">[28]</xref></td></tr><tr><td>2</td><td>k–ε models</td><td>Wall functions and kinetic energy modification</td><td>DRW model</td><td>Particles trapped at walls</td><td><xref rid="tblfn1a" ref-type="table-fn">a</xref></td><td>Over-predicted</td><td><xref rid="bib20" ref-type="bibr">[20]</xref>, <xref rid="bib29" ref-type="bibr">[29]</xref>, <xref rid="bib30" ref-type="bibr">[30]</xref></td></tr><tr><td>3</td><td>k–ε models</td><td>Wall functions</td><td>Drift-flux model</td><td>Particle concentration at wall=0 or flux=0</td><td>Acceptable</td><td><xref rid="tblfn1b" ref-type="table-fn">b</xref></td><td><xref rid="bib31" ref-type="bibr">[31]</xref>, <xref rid="bib32" ref-type="bibr">[32]</xref></td></tr><tr><td>4</td><td>k–ε models</td><td>Wall functions</td><td>Drift-flux model</td><td>Empirical deposition models</td><td>Acceptable</td><td>Depending on deposition models</td><td><xref rid="bib27" ref-type="bibr">[27]</xref>, <xref rid="bib33" ref-type="bibr">[33]</xref></td></tr><tr><td>5</td><td>Reynolds stress model</td><td>Wall functions</td><td>DRW model</td><td>Particles trapped at walls</td><td>Acceptable</td><td>Slightly over-predicted</td><td><xref rid="bib20" ref-type="bibr">[20]</xref></td></tr><tr><td>6</td><td>Large eddy simulation</td><td>Wall functions</td><td>Lagrangian method</td><td>Particles trapped at walls</td><td>Good</td><td>Good</td><td><xref rid="bib28" ref-type="bibr">[28]</xref>, <xref rid="bib34" ref-type="bibr">[34]</xref>, <xref rid="bib35" ref-type="bibr">[35]</xref></td></tr></tbody></table><table-wrap-foot><fn id="tblfn1a"><label>a</label><p>No information.</p></fn></table-wrap-foot><table-wrap-foot><fn id="tblfn1b"><label>b</label><p>Particle deposition is under-predicted if the concentration gradient is assumed to be zero, and over-predicted if the concentration is assumed to be zero.</p></fn></table-wrap-foot></table-wrap></p></sec><sec><label>4</label><title>Case set-up</title><p>A four-storey building is placed in a large computational domain (<xref rid="fig7" ref-type="fig">Fig. 7</xref>). Windows are opened on the windward side on the second and third floor. The room dimensions are height (<italic>Y</italic>)×length (<italic>X</italic>)×width (<italic>Z</italic>)=2.7 m×3.1 m×2.4 m and the window height (<italic>Y</italic>)×width (<italic>Z</italic>)=1.2 m×0.75 m. The bottom of the window is 0.8 m above the room floor. The domain boundary at <italic>x</italic>–<italic>y</italic> plane is defined as symmetry. It means the building and its surroundings are extended in the negative and positive <italic>z</italic>-direction. This treatment is based on the fact that in many high residential blocks in Hong Kong, the shape of the whole building is like a vertical slab.<fig id="fig7"><label>Fig. 7</label><caption><p>Description of the building model and the computational domain.</p></caption><graphic xlink:href="gr7"></graphic></fig></p><p>Heat is released only from the internal walls on the second and third floors. The atmospheric air temperature is 20 °C and the indoor wall surface temperature is 25 °C. In the simulation of a windless case, a uniform velocity profile (0.1 m/s) is set at the domain inlet. However, in a normal urban environment, even in a windless day in meteorology, the wind speeds are usually higher than 0.1 m/s. Therefore, a small meteorological wind speed, <italic>V</italic><sub>met</sub> (wind speed usually taken at 10 m above the ground level), at 1.0 and 4.0 m/s is considered. The wind profile in an urban environment is calculated by<disp-formula id="eq11"><label>(11)</label><mml:math id="M36" altimg="si36.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.35</mml:mn><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">met</mml:mi></mml:mrow></mml:mrow></mml:mrow></mml:msub><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mn>0.25</mml:mn></mml:mrow></mml:mrow></mml:msup><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>The turbulence on the inlet boundary is characterized by turbulence intensity and length scale, which are 8% and 1 m (for a wind speed of 1.0 m/s, the kinetic energy is 9.6×10<sup>−3</sup> m<sup>2</sup>/s<sup>2</sup>, and dissipation rate 1.55×10<sup>−4</sup> m<sup>2</sup>/s<sup>3</sup>). Given the low-speed conditions, 8% is acceptable since a general estimation of it in the wind field is 10% <xref rid="bib36" ref-type="bibr">[36]</xref>. The length scale is taken to be around 20% of the height of the two flats. Particles of 1.0, 10.0, and 20.0 μm, without an initial speed, are generated at a rate of 8 mg/s in the middle of the second floor at the height of 1.6 m. The material density of the particles is 1000 kg/m<sup>3</sup>.</p></sec><sec><label>5</label><title>Results and discussions</title><p>The airflow fields indoors and outdoors, the air change rates (ACHs) of the lower and upper floor, and the ratio of the gaseous pollutant concentrations between the floors have been discussed in an earlier publication <xref rid="bib5" ref-type="bibr">[5]</xref>. In the pure-buoyancy driven case, the ACHs through the windows are highest for both floors. A gentle approaching wind of 1.0 m/s normal to the window has very little influence on the ACH. However, 4.0 m/s wind could lower the ACH by 30–50% by suppressing the buoyancy force. In the windless situation, 7.5% of the air leaving from the window of the lower floor can re-enter into the upper floor. 1.0 m/s wind deflates the plume emitted from the lower window into the upper window, resulting in a higher re-entry ratio (10.9%). But 4.0 m/s wind confines the convection of air and thus reduces the ratio to 3.5%.</p><p>The particle concentrations are calculated based on the airflow conditions in Gao et al. <xref rid="bib5" ref-type="bibr">[5]</xref>. <xref rid="fig8" ref-type="fig">Fig. 8</xref>compares the particle concentrations and shows the concentration ratios of the upper floor to the lower floor. It is observed that the concentrations and their ratio for 1.0 μm particles are close to the values of CO<sub>2</sub>. The results of 10.0 and 20.0 μm particles differ from the 1.0 μm particle remarkably. Different from gaseous pollutants, two factors constrain the cascade effect of coarse particles: the particle’ gravity which counteracts with the upward buoyancy force and the deposition at solid surfaces. The latter one is a self-cleaning mechanism of particles which contributes to the decreased suspended concentrations in the lower floor and the reduced outlet concentrations. For an approximate estimation of particle’ settling effect, the ratio of two time scales is indicative: the time required by the particle to settle down across the room height and the time needed to refresh the total volume of room air once <xref rid="bib31" ref-type="bibr">[31]</xref>. For the current air change rate of 7.2–10.1 h<sup>−1</sup>, the ratio of the falling time of 2.7 m (room height) for particles larger than 15 μm to the air change time is comparable.<fig id="fig8"><label>Fig. 8</label><caption><p>Mean mass fractions (kg/kg) of particles in the lower and upper floors and their ratios (in the legend, “Lag ratio” is the ratio of numbers of particles deposited onto the inner surfaces of the upper floor to the total 8000 particles produced in the middle of the lower floor).</p></caption><graphic xlink:href="gr8"></graphic></fig></p><p>As wind speed increases from 0.1 to 4.0 m/s, the windward airflow field simulated by the k–ε model becomes more diffusion-dominated. The ratio of particle concentration in the upper floor to that in the lower floor decreases. In case of 4.0 m/s, the ratio is 0.3% for 20.0 μm particles. It implies that in the design of ventilation and air distribution systems, the indoor air speed level is essential to determine the suspension of particles in the room air. Air at higher speeds could entrain larger particles, and even makes deposited particles re-suspended. Although, according to the deposition model by Lai and Nazaroff <xref rid="bib25" ref-type="bibr">[25]</xref>, high friction velocity in the boundary layer will also increase the deposition velocity, this effect is negligible for particles larger than 2.5 μm because deposition onto the upward facing surfaces caused by gravitational force constitutes most of the overall deposition.</p><p>The ratios predicted from the DRW model are much less than those predicted from the drift-flux model. In spite of the different definitions of the “ratio” in the two approaches, excluding the factor that some particles may enter into the upper floor and then move out of the same floor through the window, the particulate mass in the upper floor predicted by the DRW model is still lower than that by the drift-flux model. <xref rid="fig8" ref-type="fig">Fig. 8</xref> shows that the discrepancy becomes large for fine particles, and small for coarse particles. For a close view of the results from the DRW model, <xref rid="tbl2" ref-type="table">Table 2</xref>gives the numbers of deposited particles in both flats. Even for 1.0 μm particles, more than 70% of the particles have been predicted to be absorbed by the inner surfaces of the lower flat. This result is obviously irrational. Therefore in current settings, the application of the DRW model may falsely give us an overoptimistic result with regard to the upward transport of viruses by natural ventilation airflows. In the best case when wind velocity is 4.0 m/s, three out of 8000 particles of 20.0 μm can get into the upper floor. In terms of infection risk, this may still be indicative, given that one infected person may sneeze several times and one sneeze produces more than one million particles <xref rid="bib37" ref-type="bibr">[37]</xref>.<table-wrap position="float" id="tbl2"><label>Table 2</label><caption><p>Particles deposition onto the inner surfaces of the lower and upper floors in the simulations using DRW model (totally 8000 particles are released)</p></caption><table frame="hsides" rules="groups"><thead><tr><th></th><th>Cases</th><th><italic>d</italic><sub>p</sub>=1.0 μm</th><th><italic>d</italic><sub>p</sub>=10.0 μm</th><th><italic>d</italic><sub>p</sub>=20.0 μm</th></tr></thead><tbody><tr><td rowspan="2">0.1 m/s</td><td>Lower floor</td><td align="char">5828</td><td align="char">6223</td><td align="char">7252</td></tr><tr><td>Upper floor</td><td align="char">226</td><td align="char">187</td><td align="char">76
</td></tr><tr><td rowspan="2">1.0 m/s</td><td>Lower floor</td><td align="char">5437</td><td align="char">5876</td><td align="char">7247</td></tr><tr><td>Upper floor</td><td align="char">110</td><td align="char">95</td><td align="char">77
</td></tr><tr><td rowspan="2">4.0 m/s</td><td>Lower floor</td><td align="char">7232</td><td align="char">7426</td><td align="char">7819</td></tr><tr><td>Upper floor</td><td align="char">22</td><td align="char">14</td><td align="char">3</td></tr></tbody></table></table-wrap></p><p>In principle, the Lagrangian method has an advantage to track the trajectory of each particle. <xref rid="fig9" ref-type="fig">Fig. 9</xref>shows the locations of 500 representative particles at various times after their generations. In less than 3 min, particles could enter the space of the upper floor. This may well explain the epidemiologic findings that the SARS virus was found within the deposits on the window-sill above the room whose host was infected during the SARS outbreak in Hong Kong.<fig id="fig9"><label>Fig. 9</label><caption><p>Spatial distributions of 1.0 μm particles from Lagrangian simulations in which 500 particles are generated at <italic>t</italic>=0 s in the center of the lower floor: (a) <italic>t</italic>=30 s, (b) <italic>t</italic>=70 s, (c) <italic>t</italic>=100 s, (d) <italic>t</italic>=130 s, (e) <italic>t</italic>=160 s, and (f) <italic>t</italic>=200 s.</p></caption><graphic xlink:href="gr9"></graphic></fig></p><p>In this study, the particle dispersion characteristics between two adjacent flats is studied. According to Duguid <xref rid="bib37" ref-type="bibr">[37]</xref>, in a cough more than 4900 particles are generated, and in a sneeze more than one million. Particles smaller than 24 μm make up about 90% of the total aerosol numbers, while particles greater than 24 μm contain more than 90% of the total aerosol mass. It remains unclear how many pathogens of some airborne diseases such as SARS are carried in each of the particle size bin. In any event, fine pathogen-containing particles can pose a higher infection risk via airborne inhalation of airborne residues and coarse pathogen-containing particles increase the risk via direct contact after rapid deposition.</p><p>Finally, the effect of particle source location in the lower floor and initial particle velocities are not considered here. Occupants may move frequently in their rooms. If the infected one coughs or sneezes at the window, a more serious situation than in this simulation may be envisaged. Sneezed or coughed particles usually have very high initial momentums. This makes the initial direction assumption an important factor in analyzing the movements of particles. Careful analysis of these considerations will be of the interest in the future studies.</p></sec><sec><label>6</label><title>Conclusions</title><p>In this paper, the cascade effect of particles with diameters of 1.0–20.0 μm is numerically investigated using both Eulerian and Lagrangian approaches. The following conclusions could be drawn:<list list-type="simple"><list-item><label>•</label><p>In the Lagrangian frame, tracking particles based on the time-averaged airflow field simulated by k–ε models and wall functions over-predicts the deposition rate. This over-prediction is more substantial for fine particles. Modification of the near-wall kinetic energy may alleviate this problem. However, sometimes the modification could lead to unreasonable spatial distributions of particles.</p></list-item><list-item><label>•</label><p>Both simulation approaches revealed that the cascade effect exists for particulate pollutants. The particle concentration in the upper floor is two to three orders of magnitude lower than that in the lower floor, depending on the particle sizes. 1.0 μm particles disperse like gases. Strong deposition at solid surfaces and gravitational settling of particles larger than 20.0 μm greatly limit the upward transport of them.</p></list-item><list-item><label>•</label><p>Current simulation deals with single-sided natural ventilation where buoyancy is the major driving force. The effect of wind is more complicated. A high wind velocity could suppress the buoyancy force, and limit the dispersion of particles between flats.</p></list-item></list></p><p>In this study, the airflow field is obtained by steady-state simulation using the k–ε model. But it should be noted that the buoyancy-driven naturally ventilated flows are inherently unsteady. The effect of the instantaneously fluctuating velocities on the air exchange between indoor and outdoor spaces and on the movements of particles is unable to be captured by the Reynolds-averaged Navier–Stokes modeling. Here LES can potentially provide more accurate predictions, which will be explored in our future studies.</p></sec></body><back><ref-list><title>References</title><ref id="bib1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname><given-names>I.T.S.</given-names></name><name><surname>Li</surname><given-names>Y.G.</given-names></name><name><surname>Wong</surname><given-names>T.W.</given-names></name><name><surname>Tam</surname><given-names>W.</given-names></name><name><surname>Chan</surname><given-names>A.</given-names></name><name><surname>Lee</surname><given-names>J.H.W.</given-names></name></person-group><article-title>Evidence of airborne transmission of the severe respiratory syndrome virus</article-title><source>New England Journal of Medicine</source><volume>350</volume><issue>17</issue><year>2004</year><fpage>1731</fpage><lpage>1739</lpage><pub-id pub-id-type="pmid">15102999</pub-id></element-citation></ref><ref id="bib2"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>Y.</given-names></name><name><surname>Huang</surname><given-names>X.</given-names></name><name><surname>Yu</surname><given-names>I.T.S.</given-names></name><name><surname>Wong</surname><given-names>T.W.</given-names></name><name><surname>Qian</surname><given-names>H.</given-names></name></person-group><article-title>Role of air distribution in SARS transmission during the largest nosocomical outbreak in Hong Kong</article-title><source>Indoor Air</source><volume>15</volume><issue>2</issue><year>2005</year><fpage>83</fpage><lpage>95</lpage><pub-id pub-id-type="pmid">15737151</pub-id></element-citation></ref><ref id="bib3"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>Y.</given-names></name><name><surname>Duan</surname><given-names>S.</given-names></name><name><surname>Yu</surname><given-names>I.T.S.</given-names></name><name><surname>Wong</surname><given-names>T.W.</given-names></name></person-group><article-title>Multi-zone modeling of probable SARS virus transmission by airflow between flats in Block E, Amoy Gardens</article-title><source>Indoor Air</source><volume>15</volume><issue>2</issue><year>2005</year><fpage>96</fpage><lpage>111</lpage><pub-id pub-id-type="pmid">15737152</pub-id></element-citation></ref><ref id="bib4"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Niu</surname><given-names>J.L.</given-names></name><name><surname>Tung</surname><given-names>C.W.</given-names></name></person-group><article-title>On-site quantification of re-entry ratio of ventilation exhausts in multi-family residential buildings and implications</article-title><source>Indoor Air</source><volume>18</volume><issue>1</issue><year>2008</year><fpage>12</fpage><lpage>26</lpage><pub-id pub-id-type="pmid">18093125</pub-id></element-citation></ref><ref id="bib5"><label>5</label><mixed-citation publication-type="other">Gao NP, Niu JL, Perino M, Heiselberg P. The airborne transmission of infection between flats in high-rise residential buildings-tracer gas simulation. Building and Environment 2007. <ext-link ext-link-type="doi" xlink:href="10.1016/j.buildenv.2007.10.023">doi:10.1016/j.buildenv.2007.10.023</ext-link>.</mixed-citation></ref><ref id="bib6"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brankston</surname><given-names>G.</given-names></name><name><surname>Gitterman</surname><given-names>L.</given-names></name><name><surname>Hirji</surname><given-names>Z.</given-names></name><name><surname>Lemieux</surname><given-names>C.</given-names></name><name><surname>Gardam</surname><given-names>M.</given-names></name></person-group><article-title>Transmission of influenza A in human beings</article-title><source>Lancet Infectious Diseases</source><volume>7</volume><issue>4</issue><year>2007</year><fpage>257</fpage><lpage>265</lpage><pub-id pub-id-type="pmid">17376383</pub-id></element-citation></ref><ref id="bib7"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tellier</surname><given-names>R.</given-names></name></person-group><article-title>Review of aerosol transmission of influenza A virus</article-title><source>Emerging Infectious Diseases</source><volume>12</volume><issue>11</issue><year>2006</year><fpage>1657</fpage><lpage>1662</lpage><pub-id pub-id-type="pmid">17283614</pub-id></element-citation></ref><ref id="bib8"><label>8</label><mixed-citation publication-type="other">Fluent. Fluent 6.2 user's guide. Lebanon NH: Fluent Inc.; 2005.</mixed-citation></ref><ref id="bib9"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yakhot</surname><given-names>V.</given-names></name><name><surname>Orszag</surname><given-names>S.</given-names></name></person-group><article-title>Renormalization group analysis of turbulence: I. Basic theory</article-title><source>Journal of Scientific Computing</source><volume>1</volume><issue>1</issue><year>1986</year><fpage>1</fpage><lpage>51</lpage></element-citation></ref><ref id="bib10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elghobashi</surname><given-names>E.</given-names></name></person-group><article-title>On predicting particle-laden turbulent flows</article-title><source>Applied Scientific Research</source><volume>52</volume><issue>4</issue><year>1994</year><fpage>309</fpage><lpage>329</lpage></element-citation></ref><ref id="bib11"><label>11</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hinds</surname><given-names>W.C.</given-names></name></person-group><chapter-title>Aerosol technology: properties, behavior, and measurement of airborne particles</chapter-title><year>1999</year><publisher-name>Wiley, Inc</publisher-name></element-citation></ref><ref id="bib12"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abadie</surname><given-names>M.</given-names></name><name><surname>Limam</surname><given-names>K.</given-names></name><name><surname>Allard</surname><given-names>F.</given-names></name></person-group><article-title>Indoor particle pollution: effect of wall textures on particle deposition</article-title><source>Building and Environment</source><volume>36</volume><issue>7</issue><year>2001</year><fpage>821</fpage><lpage>827</lpage></element-citation></ref><ref id="bib13"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Morawska</surname><given-names>L.</given-names></name></person-group><article-title>Droplet fate in indoor environments, or can we prevent the spread of infection?</article-title><source>Indoor Air</source><volume>16</volume><issue>5</issue><year>2006</year><fpage>335</fpage><lpage>347</lpage><pub-id pub-id-type="pmid">16948710</pub-id></element-citation></ref><ref id="bib14"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nicas</surname><given-names>M.</given-names></name><name><surname>Nazaroff</surname><given-names>W.W.</given-names></name><name><surname>Hubbard</surname><given-names>A.</given-names></name></person-group><article-title>Toward understanding the risk of secondary airborne infection: emission of respirable pathogens</article-title><source>Journal of Occupational and Environmental Hygiene</source><volume>2</volume><issue>3</issue><year>2005</year><fpage>143</fpage><lpage>154</lpage><pub-id pub-id-type="pmid">15764538</pub-id></element-citation></ref><ref id="bib15"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Loth</surname><given-names>E.</given-names></name></person-group><article-title>Numerical approaches for motion of dispersed particles, droplets and bubbles</article-title><source>Progress in Energy and Combustion Science</source><volume>26</volume><issue>3</issue><year>2000</year><fpage>161</fpage><lpage>223</lpage></element-citation></ref><ref id="bib16"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>Z.</given-names></name><name><surname>Chen</surname><given-names>Q.</given-names></name></person-group><article-title>Comparison of the Eulerian and Lagrangian methods for predicting particle transport in enclosed spaces</article-title><source>Atmospheric Environment</source><volume>41</volume><issue>25</issue><year>2007</year><fpage>5236</fpage><lpage>5248</lpage></element-citation></ref><ref id="bib17"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname><given-names>B.</given-names></name><name><surname>Yang</surname><given-names>C.</given-names></name><name><surname>Yang</surname><given-names>X.</given-names></name><name><surname>Liu</surname><given-names>S.</given-names></name></person-group><article-title>Particle dispersion and deposition in ventilated rooms: testing and evaluation of different Eulerian and Lagrangian models</article-title><source>Building and Environment</source><volume>43</volume><issue>4</issue><year>2008</year><fpage>388</fpage><lpage>397</lpage></element-citation></ref><ref id="bib18"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chao</surname><given-names>C.Y.H.</given-names></name><name><surname>Wan</surname><given-names>M.P.</given-names></name></person-group><article-title>A study of the dispersion of expiratory aerosols in unidirectional downward and ceiling-return type airflows using a multiphase approach</article-title><source>Indoor Air</source><volume>16</volume><issue>4</issue><year>2006</year><fpage>296</fpage><lpage>312</lpage><pub-id pub-id-type="pmid">16842610</pub-id></element-citation></ref><ref id="bib19"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname><given-names>Z.</given-names></name><name><surname>Mashayek</surname><given-names>F.</given-names></name></person-group><article-title>Stochastic modeling of evaporating droplets poly-dispersed in turbulent flows</article-title><source>International Journal of Heat and Mass Transfer</source><volume>47</volume><issue>19–20</issue><year>2004</year><fpage>4339</fpage><lpage>4348</lpage></element-citation></ref><ref id="bib20"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tian</surname><given-names>L.</given-names></name><name><surname>Ahmadi</surname><given-names>G.</given-names></name></person-group><article-title>Particle deposition in turbulent duct flows-comparison of different model predictions</article-title><source>Journal of Aerosol Science</source><volume>38</volume><issue>4</issue><year>2007</year><fpage>377</fpage><lpage>397</lpage></element-citation></ref><ref id="bib21"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>Z.</given-names></name><name><surname>Chen</surname><given-names>Q.</given-names></name></person-group><article-title>Experimental measurement and numerical simulation of particle transport and deposition in the ventilated rooms</article-title><source>Atmospheric Environment</source><volume>40</volume><issue>18</issue><year>2006</year><fpage>3396</fpage><lpage>3408</lpage></element-citation></ref><ref id="bib22"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Q.</given-names></name><name><surname>Squires</surname><given-names>K.D.</given-names></name></person-group><article-title>Large eddy simulation of particle deposition in a vertical turbulent channel flow</article-title><source>International Journal of Multiphase Flow</source><volume>22</volume><issue>4</issue><year>1996</year><fpage>667</fpage><lpage>683</lpage></element-citation></ref><ref id="bib23"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Y.</given-names></name><name><surname>James</surname><given-names>P.W.</given-names></name></person-group><article-title>On the effect of anisotropy on the turbulent dispersion and deposition of small particles</article-title><source>International Journal of Multiphase Flow</source><volume>25</volume><issue>3</issue><year>1999</year><fpage>551</fpage><lpage>558</lpage></element-citation></ref><ref id="bib24"><label>24</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hinze</surname><given-names>J.O.</given-names></name></person-group><chapter-title>Turbulence</chapter-title><year>1975</year><publisher-name>McGraw-Hill</publisher-name><publisher-loc>New York</publisher-loc></element-citation></ref><ref id="bib25"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lai</surname><given-names>A.C.K.</given-names></name><name><surname>Nazaroff</surname><given-names>W.W.</given-names></name></person-group><article-title>Modeling indoor particle deposition from turbulent flow onto smooth surfaces</article-title><source>Journal of Aerosol Science</source><volume>31</volume><issue>4</issue><year>2000</year><fpage>463</fpage><lpage>476</lpage></element-citation></ref><ref id="bib26"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname><given-names>W.Z.</given-names></name><name><surname>Howarth</surname><given-names>A.T.</given-names></name><name><surname>Adam</surname><given-names>N.</given-names></name><name><surname>Riffat</surname><given-names>S.B.</given-names></name></person-group><article-title>Modelling and measurement of airflow and aerosol particle distribution in a ventilated two-zone chamber</article-title><source>Building and Environment</source><volume>31</volume><issue>5</issue><year>1996</year><fpage>417</fpage><lpage>423</lpage></element-citation></ref><ref id="bib27"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>F.Z.</given-names></name><name><surname>Yu</surname><given-names>S.C.M.</given-names></name><name><surname>Lai</surname><given-names>A.C.K.</given-names></name></person-group><article-title>Modeling particle distribution and deposition in indoor environments with a new drift-flux model</article-title><source>Atmospheric Environment</source><volume>40</volume><issue>2</issue><year>2006</year><fpage>357</fpage><lpage>367</lpage></element-citation></ref><ref id="bib28"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>Y.</given-names></name><name><surname>Matida</surname><given-names>E.A.</given-names></name><name><surname>Gu</surname><given-names>J.</given-names></name><name><surname>Johnson</surname><given-names>M.R.</given-names></name></person-group><article-title>Numerical simulation of aerosol deposition in a 3-D human nasal cavity using RANS, RANS/EIM, and LES</article-title><source>Journal of Aerosol Science</source><volume>38</volume><issue>7</issue><year>2007</year><fpage>683</fpage><lpage>700</lpage></element-citation></ref><ref id="bib29"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matida</surname><given-names>E.A.</given-names></name><name><surname>Finlay</surname><given-names>W.H.</given-names></name><name><surname>Lange</surname><given-names>C.F.</given-names></name><name><surname>Grgic</surname><given-names>B.</given-names></name></person-group><article-title>Improved numerical simulation of aerosol deposition in an idealized mouth-throat</article-title><source>Journal of Aerosol Science</source><volume>35</volume><issue>7</issue><year>2004</year><fpage>1</fpage><lpage>19</lpage></element-citation></ref><ref id="bib30"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lai</surname><given-names>A.C.K.</given-names></name><name><surname>Chen</surname><given-names>F.</given-names></name></person-group><article-title>Modeling of particle deposition and distribution in a chamber with a two-equation Reynolds-averaged Navier–Stokes model</article-title><source>Journal of Aerosol Science</source><volume>37</volume><issue>12</issue><year>2006</year><fpage>1770</fpage><lpage>1780</lpage></element-citation></ref><ref id="bib31"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Murakami</surname><given-names>S.</given-names></name><name><surname>Kato</surname><given-names>S.</given-names></name><name><surname>Nagano</surname><given-names>S.</given-names></name><name><surname>Tanaka</surname><given-names>Y.</given-names></name></person-group><article-title>Diffusion characteristics of airborne particles with gravitational settling in a convection-dominant indoor flow field</article-title><source>ASHRAE Transactions</source><volume>98</volume><issue>1</issue><year>1992</year><fpage>82</fpage><lpage>97</lpage></element-citation></ref><ref id="bib32"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname><given-names>B.</given-names></name><name><surname>Guan</surname><given-names>P.</given-names></name></person-group><article-title>Modeling particle dispersion in personalized ventilated room</article-title><source>Building and Environment</source><volume>42</volume><issue>3</issue><year>2007</year><fpage>1099</fpage><lpage>1109</lpage></element-citation></ref><ref id="bib33"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname><given-names>N.P.</given-names></name><name><surname>Niu</surname><given-names>J.L.</given-names></name></person-group><article-title>Modeling particle dispersion and deposition in indoor environments</article-title><source>Atmospheric Environment</source><volume>41</volume><issue>18</issue><year>2007</year><fpage>3862</fpage><lpage>3876</lpage></element-citation></ref><ref id="bib34"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Uijttewaal</surname><given-names>W.S.J.</given-names></name><name><surname>Oliemans</surname><given-names>R.V.A.</given-names></name></person-group><article-title>Particle dispersion and deposition in direct numerical and large eddy simulations of vertical flows</article-title><source>Physics of Fluids</source><volume>8</volume><issue>10</issue><year>1996</year><fpage>2590</fpage><lpage>2604</lpage></element-citation></ref><ref id="bib35"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Béghein</surname><given-names>C.</given-names></name><name><surname>Jiang</surname><given-names>Y.</given-names></name><name><surname>Chen</surname><given-names>Q.</given-names></name></person-group><article-title>Using large eddy simulation to study particle motions in a room</article-title><source>Indoor Air</source><volume>15</volume><issue>4</issue><year>2005</year><fpage>281</fpage><lpage>290</lpage><pub-id pub-id-type="pmid">15982275</pub-id></element-citation></ref><ref id="bib36"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>X.X.</given-names></name><name><surname>Liu</surname><given-names>C.H.</given-names></name><name><surname>Leung</surname><given-names>D.Y.C.</given-names></name><name><surname>Lam</surname><given-names>K.M.</given-names></name></person-group><article-title>Recent progress in CFD modelling of wind field and pollutant transport in street canyons</article-title><source>Atmospheric Environment</source><volume>40</volume><issue>29</issue><year>2006</year><fpage>5640</fpage><lpage>5658</lpage></element-citation></ref><ref id="bib37"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duguid</surname><given-names>J.P.</given-names></name></person-group><article-title>The size and duration of air-carriage of respiratory droplets and droplet-nuclei</article-title><source>Journal of Hygiene</source><volume>44</volume><issue>6</issue><year>1946</year><fpage>471</fpage><lpage>479</lpage><pub-id pub-id-type="pmid">20475760</pub-id></element-citation></ref></ref-list><ack><title>Acknowledgments</title><p>The project was funded by the Central Earmarked Research Grant from Research Grant Committee of Hong Kong SAR government under the project no. PolyU 5125/04E and the Research Fund for Infectious Disease Control (RFIDC) under project no. 03040642.</p></ack></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/SrasV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000E649 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000E649 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= SrasV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |