Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence simulation”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulated wind farm wake sensitivity to configuration choices in the Weather Research and Forecasting model version 3.8.1

Wakes from wind farms can extend over 50 km downwind in stably stratified conditions. These wakes can undermine power production at downwind turbines, adversely impacting revenue. As such, wind farm wake impacts must be considered in wind resource assessments, especially in regions of dense wind farm development. The open-source Weather Research and Forecasting (WRF) numerical weather prediction model includes a wind farm parameterization to estimate wind farm wake effects, but model configuration choices can influence the resulting predictions of wind farm wakes. These choices include vertical resolution, horizontal resolution, and whether or not to include the addition of turbulent kinetic energy generated by the rotating wind turbines. Despite the sensitivity to model configuration, no clear guidance currently exists for these options. Here we compare simulated wind farm wakes produced by varying model configurations with meteorological observations near a land-based wind farm in flat terrain over several diurnal cycles. A WRF configuration comprised of horizontal resolutions of 3 km or 1 km paired with a vertical resolution of 10 m provides the most accurate representation of wind farm wake effects, such as the correct surface warming and elevated wind speed deficit. The inclusion of turbine-generated turbulence is also critical to produce accurate surface warming and should not be omitted.

17 WIND ENERGY↗
DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measurement of Wind Loading on Heliostats at the Crescent Dunes Power Plant: An Overview

The cost of solar collectors constitutes almost one third of the total cost of a CSP plant. One of the ongoing challenges in the design of these collectors is wind loading on mirrors, support structures, and drives. A particular challenge is dynamic wind loading, caused by the turbulent wind flow. To date, the design of solar collector structures has relied on wind tunnel experiments and numerical simulations that do not entirely capture the dynamic effects observed at scale. The CSP industry has shown increased interest in validating the idealized assumptions with measurements obtained in operational settings to improve wind load assumptions and increase reliability and cost-efficiency of the collector design. Further, performance models need realistic assumptions about wind loading and its impact on optical performance. In a parabolic trough field campaign, NREL successfully collected a wealth of long-term, high-resolution wind and loads data [2], that are publicly available and that can be used for the above-mentioned purposes. Heliostats are impacted differently by wind than parabolic troughs, due to their different shape, size, and field layout. To study the impact of wind and turbulence on heliostats, we initiated another field campaign in an operational power- tower plant, Crescent Dunes, in Nevada, USA. In this work, we present an overview of the measurements, first results, and potential implications of wind driven loads on Heliostats.

concentrated solar power↗
DOE OSTI2022

Pele: An Exascale-Ready Suite of Combustion Codes

High fidelity simulations of realistic combustion devices are extremely demanding computationally because of the requirements to capture complex fuel chemical decomposition, its intricate interactions with turbulent, often multiphase, flows, and the wide separation of space and time scales between the thin flame and the device boundaries. Software required to carry out such computations tends to be extremely complex, particularly when designed to exploit hardware accelerators, and can be difficult to port and maintain. We present Pele, a performance portable suite of tools for the simulation of combustion systems, including codes to evolve reactive multiphase configurations in the low Mach number and compressible flow regimes, along with a set of inter-compatible post processing and in situ analysis tools. The Pele suite of tools is built on top of the AMReX framework for block-structured adaptive mesh refinement, which provides efficient data structures and algorithms that enable the development of a wide variety of efficient mesh and particle based PDE integration schemes. A hierarchical MPI+X parallelism scheme supports CPU-only and accelerated architectures, where X can be OpenMP, CUDA, and HIP based approaches for intra-node computational work distribution. The algorithms and data structures underlying the Pele simulation and analysis tools are highly scalable and performant across a wide variety of high-performance computing platforms, including DOEs newest exascale-class machines, Frontier and Aurora. The simulation and analysis tools are fully documented and freely distributed as open source via GitHub. We present key algorithmic and software challenges, solution strategies, performance and resulting set of capabilities.

AMReX↗
DOE OSTI2024

Multicycle large-eddy simulations of a direct-injection hydrogen-fueled optical engine

Hydrogen (H 2 ) is a carbon-free chemical energy carrier and one promising solution for achieving effective decarbonization of the transportation sector, particularly for internal combustion engines (ICEs). With a focus on ICEs, and compared to port-fuel injection, direct injection (DI) of gaseous H 2 during the compression stroke offers potential advantages, which include backfire avoidance and reduction of preignition occurrence. In these last two decades, much research, experimental and numerical, has been devoted to understanding H 2 's mixing and combustion processes in ICEs. Computational fluid dynamics modeling efforts commonly rely on unsteady Reynolds-averaged Navier Stokes (URANS) turbulence frameworks, mostly due to their computational affordability. However, many authors have pointed out the opportunity to perform large-eddy simulations (LESs) to investigate the cyclic variability of H 2 engines and assess potential advantages of using LES in place of URANS, especially for lean operation. This study addresses this knowledge gap and presents a computational fluid dynamics (CFD) study of the H 2 DI process in an optical engine operating at relatively low tumble conditions, using multicycle LESs. In conclusion, the manuscript presents a thorough validation of the results against experimental data available from the literature as well as direct comparison with URANS, demonstrating the feasibility of multicycle LESs for CFD modeling of DI H 2 -fueled ICEs.

Direct injection↗
DOE OSTI2026

Leveraging generative artificial intelligence to bridge domain gaps in wind turbine research

A central challenge in wind turbine health monitoring is the scarcity of real-world data due to limited instrumentation, leading researchers to rely on simulation models that often suffer from reduced fidelity. However, even within simulation environments, discrepancies arise because of modeling assumptions, and configuration fidelities, creating domain gaps that limit the transferability of learned representations. Here, to investigate domain translation under controlled conditions, this project explores the use of generative artificial intelligence, specifically cycle-consistent generative adversarial networks (CGANs), to bridge the gap between OpenFAST simulation models representing 1.5 MW and 5 MW wind turbines. A physics-informed CGAN architecture is introduced, where a simplified turbine tower dynamics model is incorporated into the training loss to ensure physically consistent outputs. Quantitative results showed moderate to high agreement in frequency-domain features. Incorporating the physics-informed loss function improved the R 2 values by 30%, reduced the RMSE from 1.39 to 1.1 m/s 2 , and reduced training time by 82%. Furthermore, under increased turbulence intensity (IEC Category A), the RMSE remained stable at approximately 1.1 m/s 2 . While the present study is entirely simulation-based, it establishes a pipeline for evaluating physics-informed generative domain translation, which may serve as a foundation for future simulation-to-reality validation studies.

17 WIND ENERGY↗
DOE OSTI2024

Seasonal variability of wake impacts on US mid-Atlantic offshore wind plant power production

Abstract. The mid-Atlantic will experience rapid wind plant development due to its promising wind resource located near large population centers. Wind turbines and wind plants create wakes, or regions of reduced wind speed, that may negatively affect downwind turbines and plants. We evaluate wake variability and annual energy production with the first yearlong modeling assessment using the Weather Research and Forecasting model, deploying 12 MW turbines across the domain at a density of 3.14 MW km−2, matching the planned density of 3 MW km−2. Using a series of simulations with no wind plants, one wind plant, and complete build-out of lease areas, we calculate wake effects and distinguish the effect of wakes generated internally within one plant from those generated externally between plants. We also provide a first step towards uncertainty quantification by testing the amount of added turbulence kinetic energy (TKE) by 0 % and 100 %. We provide a sensitivity analysis by additionally comparing 25 % and 50 % for a short case study period. The strongest wakes, propagating 55 km, occur in summertime stable stratification, just when New England's grid demand peaks in summer. The seasonal variability of wakes in this offshore region is much stronger than the diurnal variability of wakes. Overall, yearlong simulated wake impacts reduce power output by a range between 38.2 % and 34.1 % (for 0 %–100 % added TKE). Internal wakes cause greater yearlong power losses, from 29.2 % to 25.7 %, compared to external wakes, from 14.7 % to 13.4 %. The overall impact is different from the linear sum of internal wakes and external wakes due to non-linear processes. Additional simulations quantify wake uncertainty by modifying the added amount of turbulent kinetic energy from wind turbines, introducing power output variability of 3.8 %. Finally, we compare annual energy production to New England grid demand and find that the lease areas can supply 58.8 % to 61.2 % of annual load. We note that the results of this assessment are not intended to make nor are they suitable to make commercial judgments about specific wind projects.

17 WIND ENERGY↗
DOE OSTI2025

Characterization of HRRR-simulated rotor layer wind speeds and clouds along the coast of California

Stratocumulus clouds, with their low cloud base and top, affect the atmospheric boundary layer wind and turbulence profile, thereby modulating wind energy resources. GOES satellite data reveal an abundance of stratocumulus clouds in the late spring and summer months off the coast of northern and central California, where there are active plans to deploy floating offshore wind farms at two lease areas (near Morro Bay and Humboldt). Since the fall of 2020, two buoys equipped with multiple instruments, including Doppler lidar, have been deployed for about 1 year in these wind farm lease areas to assess the rotor layer wind conditions in these locations. The objective of this study is to evaluate how well the High-Resolution Rapid Refresh (HRRR) model represents stratocumulus cloud characteristics and turbine-relevant rotor layer winds (surface to 300 m) by comparing HRRR simulations with buoy and satellite observations. We first find that the HRRR model reproduces the seasonal cycle of cloud top height reasonably well in these regions. However, during the warm season – especially at Morro Bay – the HRRR-simulated stratocumulus clouds tend to have lower tops by about 150 m and exhibit weaker diurnal cycles than satellite observations. Our analysis also shows that rotor layer wind speeds and vertical shear are stronger at Humboldt than at Morro Bay, and both are generally stronger under clear-sky conditions. Finally, the HRRR model bias in rotor layer wind speed is small under cloudy conditions but larger and dependent on observed wind speed under clear skies. Specifically, HRRR underestimates wind speeds at Morro Bay and overestimates them at Humboldt under clear-sky conditions.

17 WIND ENERGY↗
DOE OSTI2022

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A two-equation soot-in-flamelet modeling approach applied under Spray A conditions

Soot production (including formation and oxidation) is studied in the transient, high-pressure and turbulent n-dodecane Spray A flames from the Engine Combustion Network (ECN) using computational fluid dynamics (CFD) simulations. A two-equation soot-in-flamelet modeling approach is applied within the framework of the Unsteady Flamelet Progress Variable (UFPV) model and results are validated against experimental data. Equations for soot mass fraction and soot number density derived in the mixture fraction space are solved in the context of detailed flamelet calculations. Source terms for the different steps in the soot chemistry are tabulated and incorporated in the flamelet manifold. For the reference condition, the modeling approach based on the tabulated flamelet manifold reduces the computational cost of a CFD calculation by approximately 40 times compared to a non-tabulated well-mixed (WM) modeling approach. The soot-in-flamelet approach is then extended to study the effect of ambient oxygen concentration, ambient mixture composition and ambient temperature on soot production. Results show that the modeling approach is able to capture the experimental trends for the soot volume fraction (SVF) with good quantitative agreement, especially in the soot ramp-up region.

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

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme

Abstract. An earlier study evaluating dust life cycle in the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a revised numerical process coupling scheme that requires the least amount of code changes, in which the surface emissions are applied before turbulent mixing instead of before dry removal. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments show that the revised process coupling substantially weakens dry removal and strengthens vertical mixing in dust source regions. It also strengthens the large-scale transport from source to non-source regions, strengthens dry removal outside the source regions, and strengthens wet removal and activation globally. In transient simulations of the years 2000–2009 conducted using 1∘ horizontal grid spacing, 72 vertical layers, and unchanged tuning parameters of emission strength, the revised process coupling leads to a 40 % increase in the global total dust burden and an increase of dust lifetime from 1.8 to 2.5 d in terms of 10-year averages. Weakened dry removal and increased mixing ratios are also seen for other aerosol species that have substantial surface emissions, although the changes in mixing ratio are considerably smaller for the submicron species than for dust and sea salt. Numerical experiments confirm that the revised coupling scheme significantly reduces the strong and non-physical sensitivities of model results to vertical resolution in the original EAMv1. This provides a motivation for adopting the revised scheme in EAM as well as for further improvements on the simple revision presented in this paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Turbulent exchanges between near-inertial waves and balanced flows

Wind generated near-inertial waves are ubiquitous in the upper ocean. An improved understanding of near-inertial wave dynamics following their excitation in the ocean and their subsequent interaction with mesoscale geostrophic balanced flows is key to decoding oceanic energy flow pathways. In this regard, multiple oceanic data sets accumulated over the past few decades reveal that the relative strength of near-inertial waves and geostrophic balanced eddy fields is highly variable, both geographically and seasonally. Inspired by these observations, we investigate turbulent interactions and energy exchanges between near-inertial waves and balanced flows using freely evolving numerical simulations of the non-hydrostatic Boussinesq equations. Here, we find accelerated vertical propagation and dissipation of the waves in regimes where balanced and wave fields have comparable strengths. In such regimes we also find that near-inertial waves directly extract energy from balanced flows, with O(10%) being the amount of balanced energy loss. In contrast, we find that near-inertial waves transfer energy to balanced flows in regimes where balance-to-wave energy ratio is small, with the gain in balanced energy being dependent on the relative strength of waves. Furthermore, these regimes are characterized by relatively weaker vertical propagation and dissipation of the near-inertial wave field. One of the key outcomes of this study is the demonstration of the lack of a unique direction for near-inertial wave-balanced flow energy transfers. Depending on the balance-to-wave energy ratio, near-inertial waves can act as an energy sink or energy source for the geostrophic balanced flow.

58 GEOSCIENCES↗
DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Role of the Toroidal Vortex in Cumulus Clouds' Entrainment and Mixing

Abstract Shallow convective clouds play a crucial role in Earth's energy budget, as they modulate the radiative transfer in the atmosphere and participate in the vertical transport of aerosols, energy, and humidity. The parameterizations representing these complex, vital players in weather and climate models are mostly based on a description of steady‐state plumes and are a source of major uncertainty. Recently, several studies have shown that buoyant thermals are inherent in atmospheric convection and contain a toroidal (ring) vortex. This work studies those vortices in growing shallow cumulus (Cu) clouds using high‐resolution (10 m) Large Eddy Simulations that resolve these vortices in much detail. Recent analysis of such data showed that small‐scale turbulent diffusion is unable to explain the large diluted portion of the cloud. Here we advocate for the important role of the Cu toroidal vortex (TV) in cloud dilution and present the complex dynamics and structure of a Cu TV. Nevertheless, since the vortex dominates the cloud's dilution, simplicity emerges when considering the cloud's lateral mass flux profile. The cloud mixing is quantified using direct flux calculations and Eulerian tracers. In addition, Lagrangian tracers are used to identify the origin of the entrained air and its thermodynamic properties. It shows that most of the air entrained by the vortex is not recycled by the vortex, yet is significantly more humid than the environment. We suggest that the development of new models describing thermals, together with their toroidal vortices, might improve cloud parameterizations in weather and climate models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.

fusion plasmas↗
DOE OSTI2024

Numerical study on aerosol sampling in a nuclear facility duct with a 90-degree elbow

Due to the challenging design requirements, elbows are often unavoidable in duct configuration, and these 90-degree bends introduce swirl, velocity variations, recirculation, and secondary flow. These disturbances make it difficult for nuclear facilities to meet particle sampling standards. A series of numerical analyses are conducted to track aerosols in a nuclear facility duct having a 90-degree elbow with the assistance of computational fluid dynamics (CFD). A turbulence model, a continuity, and a momentum, a discrete phase model, and species transport equations are solved simultaneously to track aerosols in the duct. The effect of turbulence models, turbulent dispersion models, droplet drag model, aerosol amount, aerosol spray configuration, guide vanes, and mixers are investigated. Simulation results are analyzed per relevant testing codes such as DOE-HDBK-1169, ASME AG1, ISO 14644-3, ACGIH, and ANSI/HPS N13.1.

Han, Kai [Savannah River Nuclear Solutions (SRNS),↗