Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
DOE OSTI2020

Evaluation of semi-empirical soot models for nonpremixed flames with increased stoichiometric mixture fraction and strain

Due to the size and complexity of industrial-scale systems and computing limitations, semi-empirical soot models are often employed in CFD simulations rather than computationally-expensive detailed models. Many of these models were developed for specific applications with unique characteristic timescales and/or validated only under fuel-air combustion conditions. Hence, their use in different contexts, such as oxy-combustion, could lead to inaccurate predictions. In this study, twelve semi-empirical models (1-step or 2-step) are evaluated on their ability to respond to changes in stoichiometric mixture fraction (Z st ) and strain in a series of ethylene counterflow flames which span the experimental sooting-to-non-sooting (yellow to blue) transition. A unique approach of plotting soot formation rate in normalized, local equivalence ratio space is introduced to aid in the analysis. Results show that no existing model is able to predict a blue flame when Z st is increased beyond the experimentally-measured sooting limit. Many models give the counter-intuitive result of increasing peak and/or integrated soot volume fraction (svf) as Z st increased, contrary to experimental observations. Two models are able to successfully predict a blue flame upon increase of strain to the sooting limit value. For two-step models, the strong dependence of growth rate on surface area results in high sensitivity of svf to changing flame boundary conditions, but with mixed accuracy. There is still a significant need for a robust semi-empirical model which can accurately predict soot fraction in systems with oxygen-enrichment or variable strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗
DOE OSTI2025

Toward integrating high-fidelity CFD approaches in the thermal-hydraulic analysis of turbulent dry cask systems

Nuclear power plants have been supplying resilient and reliable electricity for decades, contributing to energy independence of the U.S.. However, nuclear waste management remains one of the most significant challenges in the industry. The safety of dry cask storage systems relies heavily on their thermal-hydraulic performance. Computational Fluid Dynamics (CFD) simulations are often used to demonstrate this performance and ensure that the system design meets safety standards. This study presents reduced numerical models for various types of dry cask systems. These numerical models can produce efficient and fast results based on the employed modeling strategies. Additionally, the study uses a novel approach to high-fidelity simulations to evaluate modeling assumptions in dry cask modeling. Large Eddy Simulations (LES) are used for this purpose, particularly in regions where fluid velocity is relatively high and the turbulence characteristics become important. Furthermore, the results of these high-fidelity simulations will enhance the interpretation of outcomes produced from a lower-fidelity CFD model.

CFD↗
DOE OSTI2023

Classroom aerosol dispersion modeling: experimental assessment of a low-cost flow simulation tool

The purpose of this study was to assess the utility of a low-cost flow simulation tool for an indoor air modeling application by comparing its outputs with the results of a physical experiment, as well as those from a more advanced computational fluid dynamics (CFD) software package. In this study, five aerosol dispersion tests were performed in two different classrooms by releasing a CO 2 tracer gas from six student locations. Resultant steady-state concentrations were monitored at 13 locations around the periphery of the room. Subsequently, the experiments were modeled using both a low-cost tool (SolidWorks Flow Simulation) and a more sophisticated tool (STAR-CCM+). Models were evaluated based on their ability to predict the experimentally measured concentrations at the 13 monitoring locations by calculating four performance parameters commonly used in the evaluation of dispersion models: fractional mean bias (FB), normalized mean-square error (NMSE), fraction of predicted value within a factor of two (FAC2), and normalized absolute difference (NAD). The more sophisticated model performed better in 15 of the 20 possible cases (five tests at four parameters each), with parameters meeting acceptance criteria in 19 of 20 cases. However, the lower-cost tool was only slightly worse, with parameters meeting acceptance criteria in 18 of 20 cases, and it performed better than the other tool in 3 of 20 cases. Because it provides useful results at a fraction of the monetary and training cost and is already widely accessible to many institutions, such a tool may be worthwhile for many indoor aerosol dispersion applications, especially for students or researchers just beginning CFD modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiscale CFD simulation of biomass fast pyrolysis with a machine learning derived intra-particle model and detailed pyrolysis kinetics

Coupling particle and reactor scale models is as essential as reactor fluid dynamics and particle motion for accurate Computational Fluid Dynamic (CFD) simulations of biomass fast pyrolysis reactors due to intraparticle heat transfer and chemical reactions controlling conversion time and product distributions. Direct online coupling of a particle model with a reactor model is computationally expensive, while offline coupling is case-dependent. In this research, solutions from a series of particle pyrolysis simulations were regressed with Artificial Neural Network (ANN). Furthermoer, this machine learning-derived model predicted the same temperature and conversion profiles compared with particle resolved simulation while the isothermal approach overpredicted the temperature by 130 K and underpredicted the conversion time by 30 s. The ANN model was then integrated into CFD simulations of fluidized bed biomass fast pyrolysis with varied feedstocks via coupling PyTorch and MFiX. The averaged error of simulation predicted bio-oil yields with four feedstocks is 6.4%. This multi-scale approach provides an efficient tool for the coupled particle and reactor scale simulations of biomass pyrolysis.

09 BIOMASS FUELS↗
DOE OSTI2024

Modeling Dense Particle Flow in Multistage and Obstructed Flow Receivers Using High Fidelity Simulations

Particles are a leading contender for next-generation, concentrating solar power technologies, and the design of the particle receiver is critical to minimize the levelized cost of electricity. Falling particle receivers (FPRs) are a viable receiver concept, but many new designs feature complex particle obstructions that include dense discrete phase flows. This creates additional challenges for modeling as particle-to-particle interactions (i.e., collisions) and particle drag become more complex. To improve upon existing modeling strategies, a CFD-DEM simulation capability was created by coupling two independent codes: Sierra/Fuego and LAMMPS. A suitable receiver model was then defined using a traditional continuum-based model for the air and a granular model for the particle curtain. A sensitivity study was executed using this model to determine the relevance of different granular model inputs on important quantities of interest in obstructed flow FPRs: the particle velocity and curtain opacity. The study showed that the granular model inputs had little effect on the particle velocity magnitude and curtain opacity after an obstruction.

Mills, Brantley↗
DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗
DOE OSTI2022

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗
DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗
DOE OSTI2023

Modeling the impact of the fuel injection strategy on the combustion and performance characteristics of a heavy-duty GCI engine

Gasoline compression ignition (GCI) is a promising strategy to achieve high thermal efficiency and low emissions with limited modifications to the conventional diesel engine hardware. It is a partially premixed concept, which derives its superiority from higher volatility and longer ignition delay of gasoline-like fuels combined with higher compression ratio typical of diesel engines. Here, the present study investigates the combustion process in a GCI engine operating with different injection strategies using computational fluid dynamics (CFD). Simulations are carried out on a single cylinder of a multi cylinder heavy-duty compression ignition engine, which operates at a compression ratio of 17:1 and an engine speed of 1038 rev/min. Two different injection strategies viz., late injection (LI), and early pilot injection (EP) are investigated to understand their impact on combustion and performance of the engine. Renormalized group (RNG) k-ε model is used to describe in-cylinder turbulence and KH RT model is used to simulate the fuel spray breakup. The developed CFD methodology is validated against relevant experimental data under a wide range of operating conditions for each injection strategy. The developed CFD methodology was found to capture the engine combustion behavior quite well. Based on the validated CFD model, the differences in the progress of combustion event for the two injection strategies is highlighted. It was found that a larger pilot fuel mass fraction results in a steeper rise in the initial heat release rate which in turn influences the transition to mixing controlled combustion. In line with the experimental data, the study showed that the late pilot injection strategy with three injection pulses, results in higher performance compared to the other conditions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal Analysis Evaluations Using the Dry Cask Simulator

Different thermal analysis models were developed to simulate the dry cask simulator (DCS). The DCS is an experiment designed to simulate dry storage of a single boiling water reactor fuel assembly under a variety of heat loads and internal pressures. The DCS was set up and tested in both a vertical and horizontal configuration to determine cladding temperatures in vertical and horizontal dry cask storage systems. The models included a detailed STAR-CCM+ model with the fuel assembly geometry explicitly modeled, a porous STAR-CCM+ model with the fuel assembly geometry modeled as porous media region with calculated effective properties, and a COBRA-SFS model. COBRA-SFS is a thermal-hydraulic code developed for steady-state and transient analysis of multi assembly spent-fuel storage and transportation systems. STAR-CCM+ is a commercial computation fluid dynamics (CFD) code. Both a detailed and porous STARCCM+ model were developed to look at the effective thermal conductivity (keff) approach to modeling a fuel assembly. A keff fuel model is typically model in CFD thermal analyses due to its significantly lower computational costs. The models were run for a combination of low and high canister pressures (100 kPa and 800 kPa) and low and high internal heat loads (0.5 kW and 5 kW). Results from all three models were compared against experimental data taken from the DCS for the peak cladding temperature (PCT) and inlet air mass flow.

Suffield, Sarah R.↗
DOE OSTI2025

Benchmark study of the DTU OWC chamber with both two-way and one-way absorption

This paper reports on a benchmark study based on small-scale (1:50) measurements of a single, oscillating water column chamber mounted sideways in a long flume. The geometry of the OWC chamber is extracted from a barge-like, attenuator-type floating concept “KNSwing” with 40 chambers targeted for deployment in the Danish part of the North Sea. In addition to traditional two-way energy extraction we also consider one-way energy extraction with passive venting and compare chamber response, pressures and total absorbed energy between the two methods. A blind study was established for the numerical modeling, with participants applying several implementations of weakly nonlinear potential flow theory and commercial Navier–Stokes solvers (CFD). Both compressible and incompressible models were used for the air phase. Potential flow calculations predict more energy absorption near the chamber resonance for one-way absorption than for two-way absorption, but the opposite is found from the experimental measurements. This outcome is mainly attributed to energy losses in the experimental passive valve system, but this conclusion must be confirmed by better experimental measurements. Modeling the one-way valve in CFD proved to be very challenging and only one team was able to provide results which were generally closer to the experiments. The study illustrates the challenges associated with both numerical and experimental analysis of OWC chambers. Air compressibility effects were not found to be important at this scale, even with the large volume of additional air used for the one-way case.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

CFD and systems engineering to minimize membrane-based carbon capture costs

In this study, we propose a superstructure for membrane-based carbon capture, addressing these challenges through the formulation of an optimization problem. The case study considers a flue gas with a flowrate of 1000 mol/s and a 30 % CO2 (molar basis), which is representative of steel plants. Reduced models, derived from rigorous CFD simulations of membrane modules, are utilized to enhance the efficiency of the optimization process. The reduced model includes specific input variables such as inlet flow rate, retentate and permeate pressures, and inlet CO2 concentration; while the CO2 recovery and purity in the permeate side are considered as output variables. In this way, we exploit the information from the CFD simulations in the superstructure, aiding in the selection of the optimal configuration for the multi-stage membrane process. The results underscore the efficiency of a multi-stage membrane design, showing capture costs around 45 $/t-CO2.

Pedrozo, Hector A.↗
DOE OSTI2025

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Hotel Room Computational Fluid Dynamics to Investigate Airborne Pathogen Dispersal Patterns

A hotel room unit consisting of a bedroom and bathroom was modelled using computational fluid dynamics (CFD) to investigate airborne pathogen dispersal patterns. The full-scale model includes a ‘typical’ hotel room configuration, furniture, and vents. The air sources and sinks include a bathroom vent, a heating, ventilation, and cooling (HVAC) unit located in the bedroom, and a ½” gap at the bottom of the entry door. In addition, the entry door and window can be opened or closed, as desired. Three key configuration simulations were conducted: 1) both the bathroom vent and HVAC were on, 2) only the HVAC was on, and 3) only the bathroom vent was on. If the HVAC air is from a fresh, clean source, or passes through a high-efficiency filter/UV device, then the first configuration is the safest, as contaminated air is highly reduced. The second configuration is also safe, but does not benefit from the outsourcing of potentially-infected air, such as contaminated air flowing through an ineffective filter. The third configuration should be avoided, as the bathroom vent causes air to flow from the hallway, which can be of dubious origin. The CFD simulations also showed that recirculation and swirling regions tend to accumulate the largest concentrations of heavier airborne particles, pathogens, dust, etc. These regions are associated with the largest turbulence kinetic energy (TKE) , and tend to occur in areas with flow recirculation and corners. Therefore, TKE presents a reasonable metric to guide the strategic location of pathogen mitigation devices. The simulations show complex flow patterns with distinct upper and lower flow regions, swirling flow, and significant levels of turbulent mixing. These simulations provide intriguing insights that can be applied to help mitigate pathogen aerosol dispersal, generate building design guidelines, as well as provide insights for the strategic placement of mitigation devices, such as ultraviolet (UV) light, supplemental fans, and filters.

54 ENVIRONMENTAL SCIENCES↗