Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective turbulence”

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 199 records · Page 11

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

When does turbulence spreading matter?

Few, if any, of the many papers on turbulence spreading address the key question of how turbulence spreading actually affects the profile structure. In this work, we are using a reduced model to answer that question. Turbulence spreading is most relevant near regions where the profiles support a strong intensity gradient ∇I. One such case is at the edge of an L mode ischarge, near a source of turbulence [i.e., either a localized source of edge turbulence or an influx of turbulence from the scrape-off layer (SOL)]. Another is in “No Man’s Land” (NML), which connects the pedestal to the stiff core in H mode. In the case of L mode, without an edge intensity source, the turbulence intensity profile is nearly flat and spreading has a weak effect. An edge localized source increases the edge ∇I, which then drives inward spreading. Invasion of turbulence from the SOL to the edge softens the edge pressure gradient. In H mode, the strong shear suppression of pedestal turbulence necessarily forces a sharp ∇I in NML. This sharp∇I drives a significant flux of turbulence from the core to the pedestal, where it is ultimately dissipated by shearing. Counter-intuitively, the results indicate that spreading actually increases the pedestal height and width and hence the energy content in H mode. This suggests that models of the pedestal structure should include NML turbulence spreading effects. The relation of avalanches to spreading is studied. Spreading weakly affects the avalanche distribution, but the spatiotemporal correlation of intensity increases with spreading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of radial pressure corrugations and profile shearing on turbulence in Fusion plasmas

Microturbulence can produce stationary fine-scale radial corrugations on the plasma density and temperature gradients in magnetic confinement fusion devices. We show that these structures play a significant role in regulating turbulent transport. We focus on the pedestal, studying electron-temperature-gradient (ETG) mode destabilisation and saturation in the presence of radial corrugations on the electron temperature gradient that could result from microtearing turbulence. A linear dispersion relation is derived for a shearless slab case, which indicates that in the presence of a sinusoidal background corrugation, each ETG mode splits into three distinct eigenvalues, with one being the original, one being more unstable and one being less unstable. However, despite the presence of more unstable linear modes, nonlinear gyrokinetic simulations of ETG with corrugated background electron temperature show a reduction of fluxes. Our investigation reveals a radial variation of the phase velocity of the modes that is proportional to the diamagnetic drift velocity and the local pressure gradient. The associated profile shearing breaks the turbulent eddies apart, reducing the transport level. This profile shearing resulting from fine-scale pressure corrugations could be a ubiquitous turbulence saturation mechanism not just in Fusion plasmas, but in Astrophysics and other areas.

FOS: Physical sciences↗

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗

Dynamic Wind Loading on CSP Collectors Caused by Turbulent Wind Fluctuations: Insights from a 2-Year Field Campaign

Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.

admittance functions↗

Simulating Turbulence-aided Neutrino-driven Core-collapse Supernova Explosions in One Dimension

The core-collapse supernova (CCSN) mechanism is fundamentally 3D, with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain among the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixinglength theory approach, introducing a few free parameters to the model that are fit to the results of 3D simulations. For sufficiently large values of the mixing-length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 M ⊙ . We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational and Experimental Investigation of Discrete “W” Turbulators for sCO2 Turbine Internal Cooling

Direct-fired supercritical CO2 (sCO2 ) cycles are of particular interest for power generation due to the potential for higher thermal efficiencies with inherent CO2 capture. Direct fired sCO2 turbines must withstand high temperatures and the density of the working fluid is much higher than conventional Brayton cycles. Based on conventional Brayton cycle applications, internal serpentine cooling channels with ribbed turbulator strips are a common approach to achieve the required cooling, while minimizing coolant utilization. Most of the previous studies of coolant channel heat transfer enhancements have focused on serpentine channels using air as the cooling media at Reynold’s numbers typical of Brayton cycles. This paper investigates the Nusselt number augmentation of coolant channels with no features, ribbed turbulators, and discrete “W” turbulators with surface roughness effects included and CO2 as the coolant. The numerical results are compared to experimental results using CO2 as the cooling media at Reynold’s numbers ranging from 100,000 to 300,000. The results show the Nusselt number augmentation plateaus with increasing Reynolds number, in contrast to the decay observed in prior publications with air as the coolant media. Additionally, the CFD results indicate the discrete “W” geometry provides a 10% increase in Nusselt number augmentation relative to the angled ribbed turbulators.

computational fluid dynamics (CFD)↗

Overall Cooling Effectiveness of Internally Cooled Additively Manufactured Blades

Power generation gas turbine blades must utilize highly effective cooling schemes to withstand moderate heat loads for long operating times. Prior work indicates that gas turbines in the 5-10 MW range primarily utilize internal cooling to prevent further increase of aerodynamic losses that result from film cooling, but few studies focus solely on internal cooling. In this study, three additively manufactured internal cooling designs were tested in the 50 percent span section of the National Experimental Turbine (NExT) blade: a baseline cooling design with pin-fin arrays and rib turbulators, a double wall design with impingement jets and serpentine channels, and an incremental impingement design. The designs were tested in a high-speed linear cascade at exit Reynolds numbers ranging from 700,000 to 1,400,000 at an exit Mach of 0.8. In addition, freestream turbulence and build angle effects on cooling performance were studied.

gas turbine blades↗

Turbulent entrainment in finite-length wind farms

In this article, we present an entrainment-based model for predicting the flow and power output of finite-length wind farms. The model is an extension of the three-layer approach of Luzzatto-Fegiz & Caulfield ( Phys. Rev. Fluids , vol. 3, 2018, 093802) for wind farms of infinite length, and assumes dependence of key flow quantities, such as the wind farm bulk velocity, on the streamwise distance from the farm entrance. To assist our analysis and validate the proposed model, we undertake a series of large-eddy simulations with different turbine spacing arrangements and layouts. Comparisons are also made with the top-down model with entrance effects of Meneveau ( J. Turbul. , vol. 13, 2012, N7) and data from the literature. The finite-length entrainment model is shown to be capable of capturing the power drop between contiguous rows of turbines as well as describing the advection and turbulent transport of kinetic energy in both the entrance and fully developed regions. The fully developed regime is approximated only deep in the wind farm, after approximately 15 rows of turbines. Our data suggest that for the cases considered in this study, the empirical coefficients that can be used to describe turbulent entrainment and transfers above the wind farm exhibit little dependence on the farm layout and may be considered constant for modelling purposes. However, the flow field within the wind farm layer can be strongly modulated by the turbine density (spacing) as well as the array layout, and to that extent it can be argued that they are both primary factors determining the wind farm power output.

17 WIND ENERGY↗

Partnership Center for High-fidelity Boundary Plasma Simulation

The University of Colorado Boulder contributed to the DOE Partnership Center for High‑fidelity Boundary Plasma Simulation (award DE‑SC0018271) by advancing gyrokinetic modeling of core–edge plasma interactions in tokamaks. Using the GEM and XGC codes, the project investigated turbulence spreading from the plasma edge to the core, finding that while edge‑driven trapped electron mode turbulence enhances core ion heat flux in DIII‑D L‑mode plasmas, it does not fully explain experimentally observed transport shortfalls, indicating a need for flux‑driven models. The team also demonstrated that edge micro‑instabilities are highly sensitive to small variations in magnetic equilibrium, particularly near the pedestal and edge, highlighting the importance of accurately characterizing equilibrium uncertainties. Significant progress was made in spatially coupling the core GEM and edge XGC simulations, with successful validation against standalone XGC results in realistic tokamak geometries, laying the groundwork for future extensions to kinetic electrons and electromagnetic effects.

42 ENGINEERING↗

Enhancements in Cloud Condensation Nuclei Activity From Turbulent Fluctuations in Supersaturation

Abstract The effect of aerosols on the properties of clouds is a large source of uncertainty in predictions of weather and climate. These aerosol‐cloud interactions depend critically on the ability of aerosol particles to form cloud droplets. A challenge in modeling aerosol‐cloud interactions is the representation of interactions between turbulence and cloud microphysics. Turbulent mixing leads to small‐scale fluctuations in water vapor and temperature that are unresolved in large‐scale atmospheric models. To quantify the impact of turbulent fluctuations on cloud condensation nuclei (CCN) activation, we used a high‐resolution Large Eddy Simulation of a convective cloud chamber to drive particle‐based cloud microphysics simulations. We show small‐scale fluctuations strongly impact CCN activity. Once activated, the relatively long timescales of evaporation compared to fluctuations causes droplets to persist in subsaturated regions, which further increases droplet concentrations.

54 ENVIRONMENTAL SCIENCES↗

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗

Measuring the Effects of Wake Steering With Nacelle-Mounted Scanning Lidar

In large wind plants, wakes from upwind turbines affect downwind turbines by reducing wind speeds and increasing turbulence. Wake steering may mitigate this effect by deflecting the wakes of upwind turbines away from downwind units. Here, we characterize the impact of wakes from yawed turbines at a commercial-scale wind plant under varying atmospheric and turbine operating conditions. Six months of nacelle-based lidar measurements were collected as part of a field campaign in 2019-2020 in the northern US Great Plains to test the effectiveness of wake steering. We separate these lidar scans by atmospheric stability and turbine operating condition to summarize how yawed wind turbine wakes vary with these input parameters in the atmospheric boundary layer. We summarize the impact of wake steering on various wake characteristics including velocity deficit, wake width, and wake center as retrieved from these lidar data. Yawed wakes have significantly different centerline characteristics compared to unyawed wakes, with large regions of meander in the mid wake region. Yawed wakes are also deflected farther than unyawed wakes in less turbulent conditions with velocity deficits persisting further downstream in stable atmospheric conditions. Overall, yawed wind turbine wakes are larger and wider in lower wind speed environments, and yawed wakes are deflected farther in less turbulent conditions, suggesting that wake steering is most effective in stable atmospheric stratification.

17 WIND ENERGY↗

Evaluation of Tritium Behavior in Forest Vegetation for the Purpose of Determining Appropriate Non-Zero Deposition Velocities for Tritium Oxide

Modeling and field measurements have been conducted between 2019 and 2021 to assess the practicality of using a non-zero deposition velocity for safety basis estimates of tritium oxide fate and transport modeling. A model was developed which used a complex deposition algorithm designed to assess how tritium oxide would mix based on the turbulent motions and wind speed effects that the forest canopy has on the atmosphere. The model was driven based on measurements of wind and turbulence taken from the Aiken AmeriFlux Tower which measures these properties at five levels located within and just above the forest. The model was then validated against a series of field experiments which were designed to test the model predictions and estimate the deposition velocity occurring over the forest environment at the Savannah River Site. The field releases used deuterium oxide as a surrogate for tritium oxide and was released as a fine mist which rapidly evaporated, creating a gaseous tracer in the atmosphere. Using air samplers, the elevations in deuterium concentration in the air relative to background measurements was assessed and then modeled. Generally, the numerical model tended to underpredict the amount of deuterium being mixed from above the canopy to the forest floor, indicating that the predictions it provides are still conservative relative to what was measured during the field experiments. Across a suite of modeling runs, the 95 th and 99 th percentile deposition velocities were estimated to be 1.2 and 0.7 cm s -1 , respectively. Estimated deposition velocities in the 2021 field experiments, which specifically assessed a release above the forest canopy and its mixing to the surface, predicted deposition velocities ranging from 1.75 to 6.61 cm s -1 . While these field releases do not cover all possible meteorological conditions, it seems appropriate to use a non-zero deposition velocity when performing safety-basis modeling of tritium oxide. The recommendation presented in this report is to use 1.0 cm s -1 . This is between the 95 th and 99 th percentile value estimated from the modeling study, suggesting it should be appropriate for the majority of release scenarios given the model’s apparent conservatism relative to field measurements.

54 ENVIRONMENTAL SCIENCES↗