Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TKE”

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.

54 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
DOE OSTI2024

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Observed Covariations in Boundary Layer and Cumulus Cloud Layer Processes

In this paper we examine variations in boundary-layer processes spanning the shallow to deep cumulus transition. This is accomplished by differentiating boundary layer properties on the basis of convective outcomes, ranging from shallow to deep, as observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, USA. Doppler lidar, radar, and radiosonde data are combined to determine statistical differences in boundary layer and cloud layer properties using a large sample (236) of days with a range of convective outcomes: shallow, congestus, and deep convection. In these analyses, the radar characterizes diurnal cloud depth, the lidar quantifies updraft and downdraft properties in the subcloud layer, and daily radiosonde data provides the convective inhibition (CIN). Combined, these data are used to test the hypothesis that deep convection occurs when the strength of the boundary layer turbulence (i.e., TKE) exceeds the strength of the energy barrier (i.e., CIN) at the top of the CBL. Results show that days with deep convective clouds have significantly lower vertical velocity variance and weaker updrafts within the subcloud layer. However, CIN values are also found to be significantly lower on deep convective days, allowing for these weaker updrafts to penetrate the energy barrier and reach the level of free convection (LFC). In contrast, shallow convective outcomes occur when the updrafts are strong in an absolute sense, but are weak when compared to the strength of the energy barrier. These findings support the use of the CIN/TKE framework in parameterizing convection in coarse resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Fission Product Yields of 233 U and 239 Pu by Neutron-Induced Fission at Neutron Energies from 0.18-140 MeV [Slides]

Nuclear fission is an important process with applications in astrophysics, nuclear reactors, and stockpile stewardship. Despite having been discovered over 80 years ago, fission is not fully described by a predictive model. This study will show fission product mass yields and total kinetic energies (TKE) across a large range of incident neutron energies for future refinement of fission models. 233U(n,f) data is used to verify features and structures on a previous experiment and compared against fission models. 239Pu(n,f) data will be analyzed to show new data at higher neutron energies using techniques to handle fission in the presence of alphas. 252Cf(sf) data will be used for calibration and benchmarking with a thin-backed target.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Evaluated 238 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗
DOE OSTI2022

Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]

Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.

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

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Total kinetic energy release in the fast neutron induced fission of actinide nuclei

The total kinetic energy release and fission mass distributions for the fast neutron (En = 3–100 MeV) induced fission of 232 Th, 233 U, 235 U, 237 Np, 239 Pu, 240 Pu, and 242 Pu have been measured using the LANSCE facility. The neutron energies were deduced from time-of- flight measurements. The fission fragments were detected using Si PIN diode detectors, giving us the fragment energies. The actinide targets were made by vapor deposition leading to high-quality targets, that were thin and uniform with reduced impurities. Corrections were made to the data for pulse height defect and the fragment energy loss in the target and its backing. The TKE distributions were Gaussian in shape and their mean value as a function of incoming neutron energy could be fitted with second order polynomials. In the case of 233 U and 235 U, our measurements agree with prior work. Our measurements for 232 Th are unique. Our data agree with Viola scaling. The constant position of the heavy mass peak is interpreted as being due to the influence of the N = 88 and Z = 50 shells. The GEF model predictions agree with the data in general as do the CGMF model predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The characteristics of atmospheric boundary layer height over the Arctic Ocean during MOSAiC

The important roles that the atmospheric boundary layer (ABL) plays in the central Arctic climate system have been recognized, but the atmospheric boundary layer height (ABLH), defined as the layer of continuous turbulence adjacent to the surface, has rarely been investigated. Using a year-round radiosonde dataset during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, we improve a Richardson-number-based algorithm that takes cloud effects into consideration and subsequently analyze the characteristics and variability of the ABLH over the Arctic Ocean. The results reveal that the annual cycle is clearly characterized by a distinct peak in May and two respective minima in January and July. This annual variation in the ABLH is primarily controlled by the evolution of the ABL thermal structure. Temperature inversions in the winter and summer are intensified by seasonal radiative cooling and warm-air advection with the surface temperature constrained by melting, respectively, leading to the low ABLH at these times. Meteorological and turbulence variables also play a significant role in ABLH variation, including the near-surface potential temperature gradient, friction velocity, and turbulent kinetic energy (TKE) dissipation rate. In addition, the MOSAiC ABLH is more suppressed than the ABLH during the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment in the summer, which indicates that there is large variability in the Arctic ABL structure during the summer melting season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

COURAGE S7 TBS AIRBORNE SONIC WIND AND TURBULENCE DATA

These data were collected with airborne wind instrumentation booms onboard the TBS at CoURAGE S7 during February 2025. The data include 60 Hz wind speed, component wind speed, wind direction, TKE, TI, and altitude measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

BNF M1 TBS AIRBORNE SONIC WIND AND TURBULENCE DATA

These data were collected with airborne wind instrumentation booms onboard the TBS at BNF M1. The data include 60 Hz wind speed, component wind speed, wind direction, TKE, TI, and altitude measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling

The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.

clusterwakes↗
DOE OSTI2024

CROCUS Urban Fluxes of CO₂, H₂O, and Turbulence at University of Illinois Chicago

As of May 12, 2026 this dataset is currently being versioned to include data up to April 2026. Once the versioning process is complete, new data files will be available for access. The dataset metadata will also be updated to reflect the data availability of the new data being versioned. This dataset was collected at the UIC Plant Research Laboratory in Chicago, Illinois, as part of the Community Research on Climate and Urban Science (CROCUS) Urban Integrated Field Laboratory (UIFL) project, led by Argonne National Laboratory. The site provides continuous atmospheric flux measurements, focusing on CO₂, H₂O, and heat and momentum transport in an urban setting. The data is processed at 30 minutes interval using the Eddy Covariance method and includes quality control and diagnostic data generated by EddyPro software. The data is stored in the netCDF files following CF conventions. The UIC Plant Research Laboratory is located near major highways and urban infrastructure, including buildings and parking areas. The surrounding landscape consists of a mix of turf, plants, trees, and impervious surfaces such as concrete and asphalt, making it ideal for studying urban at for studies on urban sustainability, air quality, and the effects of urbanization on atmospheric processes on urban climate dynamics, air quality, and surface-atmosphere exchanges within the city of Chicago. This dataset is funded by the U.S. Department of Energy’s Office of Science, Biological and Environmental Research (BER) program.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Fitting $\overline{\nu}$ for minor Pu isotopes

After successful fitting of prompt $\overline{\nu}$ for 235 U, 238 U, and 239 Pu(n,f) using CGMF, we will move on to the minor plutonium isotopes. Minor isotopes pose a greater challenge both because there is less data available and we do not, by default, have parametrizations in CGMF already. To mitigate these challenges— and provide consistency within CGMF—we will take a stepped approach to the optimization. Continuing from the 239 Pu work, we will then fit 241 Pu(n,f) $\overline{\nu}$, where there is also a number of experimental measurements, keeping consistency between the parameters that are included in the calculation for 239 Pu and 241 Pu. With these parametrizations settled, we can consistently optimize 240 Pu(n,f) $\overline{\nu}$. Following that step, we will move on to 242 Pu(n,f) (and increasing neutron number) and 238 Pu(n,f) (and decreasing neutron number). By including the fissioning systems in this manner, we should be able to minimize the unknown parameters in CGMF. We will possibly also be able to develop systematics for the CGMF input parameters along the Pu isotopic chain. This work can serve as a guide to broadening the reactions available in CGMF. In this short report, we first give an example of how we have updated CGMF to include 240 Pu(n,f) and 242 Pu(n,f), keeping consistency with the current 239 Pu and 241 Pu calculations, but without rigorous optimization (Sec. 2). Then, we will shown in Section 3 what experimental data exist for the various isotopes to provide some insight into why 241 Pu and 239 Pu are used as anchor points.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Wind Plant Flow Physics and Power Performance in Complex Environments: Cooperative Research and Development (Final Report)

Cornell University will partner with NLR on the topic of wind farm wake effects to improve understanding of interactions between complex atmospheric flows, terrain, and wind turbine wakes and plant efficiency. Wind plant flow simulation tools will also be validated. The work performed will help improve wind farm modeling by analyzing data, applying models, designing and performing experiments to acquire additional wind farm data, and develop better models.

17 WIND ENERGY↗
DOE OSTI2025

Turbulent Parameters by airborne measurements over BNF in March 2025

The original data were collected during the AAF Engineering Flights (AEF2025) in the vicinity of the ARM Bankhead National Forest (BNF) Atmospheric Observatory (https://www.arm.gov/capabilities/observatories/bnf ) in northwestern Alabama in March 2025. The ARM Aerial Facility ArcticShark uncrewed aerial system (UAS, https://www.arm.gov/capabilities/observatories/aaf/uas) was based at the public-use airport of Posey Field, Alabama (FAA LID: 1M4, 34.28027778° N, 87.60055556° W, 283m MSL) from March 10 through March 24, 2025. The ArcticShark UAS performed nine flights, including eight research flights over the AMF3 (BNF Main Site) and Supplemental Facilities to measure atmospheric state, turbulence, surface IR temperature and imagery, and aerosol number concentration and size distribution. The current data set presents a collection of turbulent parameters in the atmospheric boundary layer or lower free troposphere based on airborne measurement throughout the field campaign. The primary instruments used to create the current data set were the Aircraft Integrated Meteorological Measurement System (AIMMS-30) and the fine-wire thermocouple probe.

Atmosphere↗
DOE OSTI2025

CHELAX-BNF: Turbulent Parameters by airborne measurements

The original data were collected on board the ARM Aerial Facility ArcticShark uncrewed aerial system (UAS; https://www.arm.gov/capabilities/observatories/aaf/uas ) during the “Characterizing HEterogeneous Land-Atmosphere eXchanges at BNF” field campaign (CHEAX-BNF; https://arm.gov/research/campaigns/aaf2025CHELAX-BNF ). The ARM Aerial Facility ArcticShark UAS was based at the public-use airport of Posey Field, AL (FAA LID: 1M4, 34.28027778° N, 87.60055556° W, 283m MSL) from May 28 through June 23, 2025. The ArcticShark UAS performed 5 flights, including 4 research flights over the BNF Main Site (ARM Mobile Facility 3, https://arm.gov/capabilities/observatories/amf ) and Supplemental Facilities to measure atmospheric state, turbulence, surface IR temperature and imagery, aerosol number concentration, and aerosol size distribution. The current data set presents a collection of turbulent parameters in the atmospheric boundary layer or lower free troposphere based on airborne measurement throughout the field campaign. The primary instruments used to create the current data set were the Aircraft Integrated Meteorological Measurement System (AIMMS-30) and the fine-wire thermocouple probe.

Aircraft Integrated Meteorological Measurement Sys↗