Search NASA⌕ Search

SEARCH · Search NASA

Results for “RANS 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Transient Efficiency, Flexibility, and Reliability Optimization of Coal-Fired Power Plants - Final Program Review

This program developed an advanced model-based monitoring and model-predictive control algorithms for a coal fired power plant (CFPP), and deployed these algorithms in a real-time platform to demonstrate performance benefits for transient flexibility and plant operation efficiency. More specifically, the objectives were successfully achieved through a combination of (i) developing a high-fidelity transient plant model in Apros, which was used as a high-fidelity plant simulation between $100-50\% TMCR$ where TMCR denotes the turbine maximum continuous rating, i.e., baseload, (ii) developing a very fast physics-based reduced-order model (ROM) of the plant, which ran more than $100\times$ faster than real-time, enabling its use as real-time embedded model for model-based estimation (MBE) and model predictive control (MPC) (iii) implementing a real-time MBE based on ROM using a robust unscented Kalman filter (UKF) to continuously tune the ROM to match the measurements from high-fidelity Apros plant model despite significant plant-model mismatch, and thus, obtain a Digital Twin of the plant (iv) designing and implementing a real-time MPC with dual objectives of transient plant load tracking with high ramp rates and minimizing coal consumption, i.e., improving plant efficiency in the baseload-partload operation range of $100-50\% TMCR$. Each key element above was developed and tested individually, and has been reported in corresponding Topical Reports. Finally, all the individual elements were integrated in an overall closed-loop system, that was successfully tested in desktop Simulink test harness simulations with ROM or high-fidelity model as the plant. Thereafter, the Simulink implementation was used to auto-generate C-code and deploy as real-time Docker microservice containers in Linux, and validate that they can run in real-time in the hardware-in-the loop (HIL) setup and produce the same results as in Simulink. The results of the integrated simulation tests in Simulink as well as the real-time HIL deployment are documented in this final report, showing good load tracking for load ramps at $3-4\%/min$ ramp rates, and achieving up to $5.5\%$ reduction in coal relative to baseline operation at $50\% TMCR$ load. The desktop and HIL simulations show successful performance of the overall model based estimation and control solution and achieve the key objectives of the program for flexible, efficient and reliable operation of sub-critical coal fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
DOE OSTI2023

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗
DOE OSTI2022

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗
DOE OSTI2022

Updates to Nek5000: RANS wall functions, NRC support, and documentation

As a key part of the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program (NEAMS), the next generation of tools for nuclear reactor design are being developed. A primary metric for success of this program is the adoption of those tools by industry stakeholders. To help meet this metric for Nek5000, development efforts this year have focused on the expansion of the available Reynolds-averaged Navier-Stokes (RANS) turbulence models to include wall function models, continued support of the U.S. Nuclear Regulatory Commission in simulating a hydrogen mitigation benchmark experiment, and expansion of the Nek5000 documentation to enhance the code’s usability. The implementation of wall functions is intended to significantly decrease time to-solution for a wide range of problems and it has been identified as an important feature by stakeholders across the nuclear community. The initial implementation of standard and pressure corrected wall functions has been tested in Nek5000 on the classical problems of channel flow and a backward facing step. It has also been tested on a 2D molten salt fast reactor core and a T-junction. The recommended formulation of choosing y + = 30 as a boundary value combined with a Neumann-Neumann formulation of the boundary condition for k and τ produces the most consistent results compared to a wall resolved approach. The collaboration with the NRC is a continuing exercise that has been ongoing for multiple years. This year we have concluded investigations into appropriate inlet conditions, showing that fully developed turbulent conditions are adequate for the full PANDA domain. We have also concluded simulations for an unobstructed jet and are now refocusing to the original obstructed jet case. The expansion of the code documentation has also been an ongoing effort over the past few years to address feedback from the various stakeholders. The documentation now more fully covers many of the available features and closes identified SQA gaps. Finally we report on three training sessions that were offered over the past year in an effort to expand the Nek5000 user base.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

WUS324: Multiscale Full Waveform Inversion Approaching Convergence Improves Waveform Fits While Imaging Seismic Structure of the Western United States

Abstract We report a new model of radially anisotropic crustal and upper mantle structure of the western United States (WUS324) obtained from full waveform inversion of earthquake data. We ran three multiscale inversion stages beyond model WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) allowing them to approach convergence to fit a larger data set to a shorter minimum period of 16 s. WUS324 is based on 324 total iterations from its starting model, significantly more (16 times) than previous studies. Waveform misfit reductions are 66%–70% for both the inversion data and an independent validation data set providing confidence in the predictive power of the model. WUS324 provides much better fits and reveals shear wavespeed, v S , structure of this large region with more detail than previous waveform tomography models. We show representative images demonstrating the resolution of diverse seismic structure across this highly heterogeneous region including oceanic lithosphere, subducting slabs and continental magmatism.

58 GEOSCIENCES↗
DOE OSTI2023

A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model

Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

StaNdaRT: a repository of standardised test models and outputs for supernova radiative transfer

We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardised test models is simulated by currently used RT codes. We ran a total of ten codes on a set of four benchmark ejecta models of Type Ia SNe. We consider two sub-Chandrasekhar-mass (M tot = 1.0 M ⊙ ) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several Python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe the test models, radiative-transfer codes, and output formats in detail, and provide access to the repository. We present example results of several key diagnostic features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns: Modeling Archive

This Modeling Archive is in support of a TES-SFA publication “Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns” (Craig et al., 2021). We ran and evaluated a multi-assumption soil organic carbon (SOC) model to investigate whether alternative assumptions regarding constraints on soil microbial biomass could lead to soil carbon saturation patterns. We developed this model in the Multi-Assumption Architecture and Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_Craig2021; Walker et al. 2018). Using MAAT, we embedded three alternative hypotheses in a microbially explicit three-pool SOC model: 1) the efficiency of mineral-associated SOC formation decreases as mineral-associated SOC approaches a maximum value (“Mineral saturation”), 2) the microbial biomass turnover rate increases with increasing microbial biomass (“Density-dependent turnover”), and 3) community carbon use efficiency decreases as microbial biomass increases toward an upper limit (“Density-dependent growth”). We ran a factorial combination of these hypotheses resulting in eight models for three different classes of model (linear decay, Michaelis-Menten decay, or reverse Michaelis-Menten decay), resulting in 24 models, 12 of which are presented or discussed in the related publication. Models were parameterized using values from previous studies with similar models (Wang et al. 2013, Wieder et al. 2014, Li et al. 2014, Georgiou et al. 2017, Hassink and Whitmore 1997) and ran to an approximate steady state (200 years) at six (6) different C input rates corresponding to 0.5, 1, 2, 4, 7, and 10 times the default input value. Further model details are available in the related publication. This archive contains output from three MAAT simulations, and scripts to run these simulations and process and plot the data. Simulations are labeled “lin”, “MM_highKm”, and “RMM_highKm” reflecting factorial runs for linear, Michealis-Menten, and reverse Michaelis-Menten models, respectively. This archive contains: • 3 R scripts prepended with “init_MAAT_” to initialize runs (1 for each simulation), • 1 csv file containing years over which to run simulations (“met_year.csv”), • 1 bash script (.bs) to run MAAT, • 6 XML files that are output from MAAT describing a run (2 for each simulation), • 3 model output csv files prepended by “out_” (1 for each simulation), and • 1 analysis R script for reproducing figures 3 and 4 in Craig et al. 2021. See included user guide (Craig_2021_modeling_archive_20210315.pdf) for file organization details.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Floods due to Atmospheric Rivers along the U.S. West Coast: The Role of Antecedent Soil Moisture in a Warming Climate

Precipitation extremes are projected to become more frequent along the U.S. West Coast due to increased atmospheric river (AR) activity, but the frequency of less intense precipitation events may decrease. Antecedent soil moisture (ASM) conditions can have a large impact on flood responses, especially if prestorm precipitation decreases. Taken together with increased antecedent evaporative demand due to warming, this would result in reduced soil moisture at the onset of extreme precipitation events. In this work, we examine the impact of ASM on AR-related floods in a warming climate in three basins that form a transect along the U.S. Pacific Coast: the Chehalis River basin in Washington, the Russian River basin in Northern California, and the Santa Margarita River basin in Southern California. Additionally, we ran the Distributed Hydrology Soil Vegetation Model (DHSVM) over the three river basins using forcings downscaled from 10 global climate models (GCMs). We examined the dynamic role of ASM by comparing the changes in the largest 50, 100, and 150 extreme events in two periods, 1951–2000 and 2050–99. In the Chehalis basin, the projected fraction of AR-related extreme discharge events slightly decreases. In the Russian basin, this fraction increases, however, and more substantially so in the Santa Margarita basin. This is due to increases in AR-related extreme precipitation events, as well as the fact that the relationship of extreme precipitation to extreme discharge is strengthened by projected increases in year-to-year volatility of annual precipitation in California, which increases the likelihood of concurrent occurrence of large storms and wet ASM conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2023

The myth of URANS

Since the 1990s, RANS practitioners have observed spontaneous unsteadiness in RANS simulations. Some have suggested deliberately using this as a method of resolving large turbulent structures. However, to date, no one has produced a theoretical justification for this unsteady RANS (URANS) approach. Here, we extend the dynamical system fixed point analysis to create a theoretical model for URANS dynamics. The results are compared to URANS simulations for homogeneous isotropic decaying turbulence. The model shows that URANS can predict incorrect decay rates and that the solution tends towards steady RANS over time. Similar analysis for forced turbulence shows a fixed modelled energy of about 30% of total energy, regardless of the model parameters. The same analysis can be used to show how hybrid type models can begin to address these issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Illuminating the pathways to carbon liberation: a systems approach to characterizing the consequential unknowns of carbon transformation and loss from thawing permafrost peatlands (Final Report)

The IsoGenie3 Project delivered new systems-level insights into carbon cycling in thawing permafrost landscapes, with an emphasis on methane and carbon dioxide emissions. From >200 samples from the site collected over a decade, co-analyzed for geochemistry and microbiology, the team recovered ~1,500 assembled microbial genomes and ~1,900 viral population genomes, revealing appreciable genetic novelty - from a new highly abundant bacterial phylum, to novel methane consumers and their activities, to rampant viral novelty. IsoGenie3 linked these organisms to carbon compound transformations (which define the cycling of organic matter in soils, and the loss of the greenhouse gases carbon dioxide and methane), and saw that the microbes at each stage of permafrost thaw had different genetic potential to degrade categories of compounds, expressed that genetic potential differently, and actually transformed carbon compounds into greenhouse gases in different ways. IsoGenie 3 identified that some of the thaw-stage differences were due to plant-microbiome relationships; the plant species across the thaw gradient contributed different carbon compounds into the soil, and hosted distinct microbiota (differing among parts of plants as well as species). Lastly, microbes in the saturated post-thaw conditions appeared likely to contribute to the mobilization and toxification of mercury released during thaw. In parallel with ongoing field sampling and analysis, hypotheses arising from field observations were tested via lab incubation experiments. When communities are taken out of their native habitats, they behave differently, and the team first rigorously quantified the magnitude of this effect on microbiome composition and functional capacity, organic matter composition, and gas production; overall the main system processes were maintained in the lab incubations under the conditions tested. Further, the microbial data could inform geochemical reaction network models of those processes. Then, the team ran experiments with additions of compounds, varying temperature, and “live” vs. “dead” peat (the latter having been gamma irradiated, with a few additional variants to control for methodological artifacts). From these, we (a) determined the importance of plant-derived soluble phenolic compounds in bogs’ extraordinary recalcitrance of organic matter, and carbon gas emissions skewed to carbon dioxide; (b) proposed an abiotic ‘tanning’ mechanism, which could contribute to Sphagnum’s inhibitory effect on anaerobic decomposition through alteration of N availability. IsoGenie3 illuminated longer-term and landscape-scale interactions of permafrost thaw and carbon cycling, advancing knowledge of the drivers of methane dynamics not only across in the permafrost-associated peatland (where hydrology and plant communities dictate microbiomes) but also their interconnected lakes (where sediment carbon quality and resident microbiota are determined by position within lake, and lake features). By leveraging observations of site methane dynamics extending well before this project, the team was able to construct a 44-year portrait of the interplay of permafrost thaw, hydrology, vegetation dynamics, and carbon gas emissions, and the doubling of the fully-thawed fens over this time. From the detailed study of this focal site, IsoGenie3 also aimed to improve model representation of these kinds of sites and processes. To improve predictions of methane transformations, we incorporated acetate and isotope dynamics into the ‘DNDC’ biogeochemistry model. In addition, recovered genomes were grouped into ‘functional groups’, i.e. the genomes that perform a specific function of interest, then used to parameterize maximum growth rate and optimum growth temperature (via signatures in their sequence composition) for the BioCrunch model. The BioCrunch model was then in turn used to test the impact of increasing functional resolution of the microbes, on the carbon gas emissions. Lastly for modeling, the ecosys model was parameterized from the microbial and other data, and used to evaluate drivers of e.g. change in methane emissions. Finally, this project also led to the development of a range of new methods and tools, a new metric of organic matter decomposability, as well as a graph-database solution to multidisciplinary data storage and querying. This project’s ongoing analyses at our focal site also contributed to broader advancements in understanding elements of genetic plasticity and methane metabolism, climate change microbiology and community assembly, global peatland geochemistry and Arctic lakes’ roles in climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Film Cooling Performance Predictions For Air And Supercritical CO 2

The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling geometries. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing a supercritical fluid. The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. A computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state RANS equations were solved along with the SST kω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness will be presented at different operating conditions. Recommendations for future studies will be given wherever applicable.

42 ENGINEERING↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal adaptation of soil microbial growth traits in response to chronic warming

ABSTRACT Adaptation of soil microbes due to warming from climate change has been observed, but it remains unknown what microbial growth traits are adaptive to warming. We studied bacterial isolates from the Harvard Forest Long-Term Ecological Research site, where field soils have been experimentally heated to 5°C above ambient temperature with unheated controls for 30 years. We hypothesized that Alphaproteobacteria from warmed plots have (i) less temperature-sensitive growth rates; (ii) higher optimum growth temperatures; and (iii) higher maximum growth temperatures compared to isolates from control plots. We made high-throughput measurements of bacterial growth in liquid cultures over time and across temperatures from 22°C to 37°C in 2–3°C increments. We estimated growth rates by fitting Gompertz models to the growth data. Temperature sensitivity of growth rate, optimum growth temperature, and maximum growth temperature were estimated by the Ratkowsky 1983 model and a modified Macromolecular Rate Theory (MMRT) model. To determine evidence of adaptation, we ran phylogenetic generalized least squares tests on isolates from warmed and control soils. Our results showed evidence of adaptation of higher optimum growth temperature of bacterial isolates from heated soils. However, we observed no evidence of adaptation of temperature sensitivity of growth and maximum growth temperature. Our project begins to capture the shape of the temperature response curves, but illustrates that the relationship between growth and temperature is complex and cannot be limited to a single point in the biokinetic range. IMPORTANCE Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗
DOE OSTI2025

Hydroboost

HydroBoost is the most realistic revenue optimization tool for the hybridization of hydropower and battery energy storage systems to date. The innovative representation of how operators actually schedule hydropower in practice results in more realistic predictions of revenue and operations. Unlike other optimization tools, HydroBoost generates forecast energy prices with uncertainty to use in the optimization. This allows HydroBoost to give users a range of potential revenue with an upper bound using the perfect foresight pricing and a lower bound using a naive persistence forecast model. Additional forecast can be generated and used in the optimization, such as additive models, random forest, and neural networks to give further insight into potential revenue. HydroBoost has been designed to be applicable for both run-of-river and reservoir storage sites. The primary focus is on the day-ahead market and requires year-long data with an hour time-step. All time-series input and constraints are contained in an Excel worksheet for convince. The user will run the forecasting generation first with a Python script to give the optimization model the necessary requirements. Next the optimization is ran using Julia and results are generated and stored into a directory as csv files. HydroBoost includes an additional module to generate figures based on the results of the optimization simulation. The results help analyze the results and users to draw insights into how the hydro and battery systems are operated and the revenue each is producing. Additionally, the difference between the perfect foresight model and models that include forecast can easily be inspected.

Phillips, TylerB. [Idaho National Laboratory (INL)↗