The present study investigates aerosol transport and surface deposition in a realistic classroom environment using computational fluid-particle dynamics simulations. Effects of particle size, aerosol source location, glass barriers, and windows are explored. While aerosol transport in air exhibits some stochasticity, it is found that a significant fraction (24%–50%) of particles smaller than 15 µ m exit the system within 15 min through the air conditioning system. Particles larger than 20 µ m almost entirely deposit on the ground, desks, and nearby surfaces in the room. Source location strongly influences the trajectory and deposition distribution of the exhaled aerosol particles and affects the effectiveness of mitigation measures such as glass barriers. Glass barriers are found to reduce the aerosol transmission of 1 µ m particles from the source individual to others separated by at least 2.4 m by ∼92%. By opening windows, the particle exit fraction can be increased by ∼38% compared to the case with closed windows and reduces aerosol deposition on people in the room. On average, ∼69% of 1 µ m particles exit the system when the windows are open.

No abstract

Identifying economically viable intervention measures to reduce COVID-19 transmission on aircraft is of critical importance especially as new SARS-CoV2 variants emerge. Computational fluid-particle dynamic simulations are employed to investigate aerosol transmission and intervention measures on a Boeing 737 cabin zone. The present study compares aerosol transmission in three models: (a) a model at full passenger capacity (60 passengers), (b) a model at reduced capacity (40 passengers), and (c) a model at full capacity with sneeze guards/shields between passengers. Lagrangian simulations are used to model aerosol transport using particle sizes in the 1–50 μ m range, which spans aerosols emitted during breathing, speech, and coughing. Sneeze shields placed between passengers redirect the local air flow and transfer part of the lateral momentum of the air to longitudinal momentum. This mechanism is exploited to direct more particles to the back of the seats in front of the index patient (aerosol source) and reduce lateral transfer of aerosol particles to other passengers. It is demonstrated that using sneeze shields on full capacity flights can reduce aerosol transmission to levels below that of reduced capacity flights without sneeze shields.

Background Human snores are caused by vibrating anatomical structures in the upper airway. The glottis is a highly variable structure and a critical organ regulating inhaled flows. However, the effects of the glottis motion on airflow and breathing sound are not well understood, while static glottises have been implemented in most previous in silico studies. The objective of this study is to develop a computational acoustic model of human airways with a dynamic glottis and quantify the effects of glottis motion and tidal breathing on airflow and sound generation. Methods Large eddy simulation and FW-H models were adopted to compute airflows and respiratory sounds in an image-based mouth-lung model. User-defined functions were developed that governed the glottis kinematics. Varying breathing scenarios (static vs. dynamic glottis; constant vs. sinusoidal inhalations) were simulated to understand the effects of glottis motion and inhalation pattern on sound generation. Pressure distributions were measured in airway casts with different glottal openings for model validation purpose. Results Significant flow fluctuations were predicted in the upper airways at peak inhalation rates or during glottal constriction. The inhalation speed through the glottis was the predominating factor in the sound generation while the transient effects were less important. For all frequencies considered (20–2500 Hz), the static glottis substantially underestimated the intensity of the generated sounds, which was most pronounced in the range of 100–500 Hz. Adopting an equivalent steady flow rather than a tidal breathing further underestimated the sound intensity. An increase of 25 dB in average was observed for the life condition (sine-dynamic) compared to the idealized condition (constant-rigid) for the broadband frequencies, with the largest increase of approximately 40 dB at the frequency around 250 Hz. Conclusion Results show that a severely narrowing glottis during inhalation, as well as flow fluctuations in the downstream trachea, can generate audible sound levels.

No abstract

Despite extensive efforts in studying radioactive aerosols, including the transmission of radionuclides in different chemical matrices throughout the body, the internal organ-specific radiation dose due to inhaled radioactive aerosols has largely relied on experimental deposition data and simplified human phantoms. Computational fluid-particle dynamics (CFPD) has proven to be a reliable tool in characterizing aerosol transport in the upper airways, while Monte Carlo based radiation codes allow accurate simulation of radiation transport. The objective of this study is to numerically assess the radiation dosimetry due to particles decaying in the respiratory tract from environmental radioactive exposures by coupling CFPD with Monte Carlo N-Particle code, version 6 (MCNP6). A physiologically realistic mouth-lung model extending to the bifurcation generation G9 was used to simulate airflow and particle transport within the respiratory tract. Polydisperse aerosols with different distributions were considered, and deposition distribution of the inhaled aerosols on the internal airway walls was quantified. The deposition mapping of radioactive aerosols was then registered to the respiratory tract of an image-based whole-body adult male model (VIP-Man) to simulate radiation transport and energy deposition. computer codes were developed for geometry visualization, spatial normalization, and source card definition in MCNP6. Spatial distributions of internal radiation dosimetry were compared for different radionuclides (131 i, 134,137 cs, 90 Sr-90 Y, 103 Ru and 239,240 Pu) in terms of the radiation fluence, energy deposition density, and dose per decay. Radioactive aerosols arise from various sources such as nuclear accidents, natural decay processes, and the decommissioning of nuclear reactors 1-3. Highly radioactive micro-particles were released to the surrounding environment in the Chernobyl and Fukushima Daiichi accidents 1,4. The aerosol particles released from accidents and natural decay processes may remain suspended in the air for extended periods, be incorporated into soil particles that can reenter the air, or be inhaled by humans or animals. When aerosol particles are inhaled, a fraction of the inhaled particles deposit in the respiratory tract, while the rest is exhaled out 5-9. The deposition fraction is dependent on many parameters such as the geometry of the respiratory airways, the particle size of the inhaled aerosol, and the breathing condition 10,11. Inhaled radioactive aerosols can expose internal organs to radiation for extended periods and may induce a spectrum of functional and morphological changes, such as genetic mutations and carcinogenesis 12,13. Their severity can be related to the absorbed radiation dose in the different specific tissues, which is further related to the number of radioactive particles inhaled, how long they stay, and the type of radiation they emit. Unlike fat and muscles being the main internal defense to external exposures, there is nothing to protect cells and tissues of sus...

Simulations are performed using the reverse non-equilibrium molecular dynamics (rNEMD) method and the Stillinger-Weber (SW) potential to determine the input parameters for achieving ±1% convergence of the calculated thermal conductivity of silicon. These parameters are then used to investigate the effects of the interatomic potentials of SW, Tersoff II, Environment Dependent Interatomic Potential (EDIP), Second Nearest Neighbor, Modified Embedded-Atom Method (MEAM), and Highly Optimized Empirical Potential MEAM on determining the bulk thermal conductivity as a function of temperature (400–1000 K). At temperatures > 400 K, data collection and swap periods of 15 ns and 150 fs, system size ≥6 × 6 UC2 and system lengths ≥192 UC are adequate for ±1% convergence with all potentials, regardless of the time step size (0.1–0.5 fs). This is also true at 400 K, except for the SW potential, which requires a data collection period ≥30 ns. The calculated bulk thermal conductivities using the rNEMD method and the EDIP potential are close to, but lower than experimental values. The 10% difference at 400 K increases gradually to 20% at 1000 K.

No abstract