Search NASA⌕ Search

SEARCH · Search NASA

Results for “RANS”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗
DOE OSTI2021

Bedrock architecture, soil texture, and hyporheic zone characterization combining electrical resistivity and induced polarization imaging

The shallow (<15 m depth) subsurface environment of a short reach along a bedrock stream was investigated with electrical resistivity and induced polarization (IP) to map details of bedrock and soil sediments. The bedrock is mostly comprised of limestone, and is generally resistive compared to the overlying soil. The soil-bedrock interface was determined through a trial and error approach using a sharp boundary feature in the inversion model. The inferred bedrock surface determined from the inversions exhibited undulating patterns with troughs and ridges. A near continuous trough ran alongside the stream within the floodplain, and is suggestive of a paleochannel among other interpretations of this feature. The structure of the electrical resistivity above the bedrock showed small-scale elongated features. The chargeability from the IP method showed larger scale features. High values of chargeability were associated with the sediments in the floodplain, and low values were associated with bedrock, stream, and soil on the elevated banks above the stream. If the chargeability is associated with membrane polarization characteristic of clayey soils, then IP seems to highlight the mere existence of clay, while the resistivity may be more discernable of the relative proportion of clay. Vegetation differences may also explain the chargeability distribution, where parts of the survey with high chargeability had dense pine with no understory making the soils more organically rich.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A second-order distributed memory parallel fast sweeping method for the Eikonal equation

The Eikonal equation is used to calculate wave propagation and distance fields, and due to its complexity requires numerical treatment for its solution. In this work, we present a second-order distributed memory parallel fast sweeping method. The second-order solution switches on a two-point stencil when two upwind points are available, and reverts to first-order otherwise. In all examples, the second-order method improves the solution over the first-order, allowing for significant savings in memory while achieving the same accuracy. Parallelization over distributed memory saw good weak scaling with optimal convergence. The computational time for second-order was approximately 2.5 times slower than first-order, where the largest amount of mesh points ran on 144 cores (512 GB) was ≈20 billion. The savings in memory from the second-order method combined with the distributed memory algorithm result in the ability to solve problems much larger than are possible with the serial first-order method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗
DOE OSTI2025

Modeling gas, hydrates, and slope stability on the U.S. Atlantic margin during Pleistocene glacial cycles

Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.

Hydrate dissociation↗
DOE OSTI2024

A flexible data acquisition system architecture for the Nab experiment

Here, the Nab experiment will measure the electron–neutrino correlation and Fierz interference term in free neutron beta decay to test the Standard Model and probe Beyond the Standard Model physics. Using National Instrument’s PXIe-5171 Reconfigurable Oscilloscope module, we have developed a data acquisition system that is not only capable of meeting Nab’s specifications, but flexible enough to be adapted in situ as the experimental environment dictates. The L1 and L2 trigger logic can be reconfigured to optimize the system for coincidence event detection at runtime through configuration files and LabVIEW controls. This system is capable of identifying L1 triggers at a rate of at least 1 MHz, while reading out a peak signal rate of approximately 2 GB/s. During the commissioning phase of the experiment, the system ran at a sustained readout rate of 400 MB/s of detector signal data originating from roughly 6 kHz L2 triggers, well within the peak performance of the system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Validation of URANS and STRUCT-ε turbulence models for stratified sodium flow

Simulations of a transient stratified sodium experiment are carried out using a classic unsteady RANS model and the second-generation URANS model, STRUCT-ε. Turbulence modeling challenges and their implications to stratified flow prediction are discussed in the context of other sources of error. Input errors are discussed and addressed; discretization error is calculated to be less than 5% of the inlet velocity, for 80% of the domain; and remaining errors in temperature distributions are attributed to the turbulence model. Qualitative flow features from the simulations are presented and discussed. Compared to the experiment, the STRUCT-ε turbulence model provides a more physically accurate prediction of temperature and momentum mixing in key regions of the domain. Quantitative measures such as the L 2 norm of the temperature discrepancy demonstrate the improved performance of the STRUCT-ε approach. In conclusion, the magnitude of the temperature fluctuations is very well-predicted by the STRUCT-ε, while URANS overpredicts them by approximately 50%.

42 ENGINEERING↗
DOE OSTI2026

Comparison of URANS and LES predictions for the open phase of the OECD NEA CSNI fluid structure interaction CFD benchmark

The OECD NEA CSNI WGAMA CFD Task Group ran a benchmark in 2020 and 2021 to assess the predictive capabilities of coupled fluid structure interaction (FSI) CFD analysis methods. This paper presents the predictions made for the open phase of the benchmark using URANS and LES turbulence modelling approaches, and a comparison of the results to the experimental data. The benchmark comprised a channel containing two inline cylinders in cross-flow. The cylinders were fixed at one end, free at the other, and had measured resonant frequencies and damping properties. The URANS modelling used ANSYS Fluent 2-way coupled to ANSYS Mechanical. The LES modelling used Nek5000, 1-way coupled to Diablo. Comparisons with cross-channel velocity profiles are presented, both for the mean flow and its RMS. Comparisons are also made to the frequency spectra for point measurements of fluid velocity and pressure, and for the accelerations of the free end of each cylinder. URANS predicts the average velocity profiles relatively well, and is able to predict the velocity and acceleration spectra at the shedding frequency. However, the frequency content at the 4th harmonic of the shedding frequency is low in the URANS flow fields, and so does not excite accelerations at the resonant frequency of the cylinders. LES makes better predictions of the average profiles, and the velocity spectra agree well at both the shedding frequency and at higher frequencies. In conclusion, the 1-way coupled LES results show good agreement for acceleration spectra.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A review of computational methods for studying oscillating water columns – the Navier-Stokes based equation approach

This review evaluates the state-of-the-practice numerical tools used to predict the performance of Oscillating Water Column (OWC). The OWC is a widely studied Wave Energy Converter that provides a reliable form of renewable form of electricity that can potentially meet global energy needs. However, the fluid-flow phenomena affecting its hydrodynamic performance are not fully understood. While there has been the successful full-scale deployment of OWCs, various computational methods are being explored to optimize this technology. Potential flow theory is commonly used to evaluate the efficiency of OWCs; however, this assumption tends to over-predict the hydrodynamic performance. Recently, numerical studies using a diverse set of commercial, open-source, or in-house Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier Stokes (RANS) and Large-Eddy Simulation codes show a better comparison to available experimental results but are computationally expensive. ANSYS Fluent was found to be the most widely used CFD code applied to the study of the OWCs, with a high degree of accuracy in terms of experimental validation of numerical results.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Effects of boundary hydraulics, dissolved oxygen, and dissolved organic carbon on growth and death dynamics of aerobic microbes in riverbed dune-induced hyporheic zones

Surface and groundwater interact in the hyporheic zone beneath and adjacent to rivers in the presence of a diverse microbial community. Heterotrophic bacteria mediate a range of environmentally important reactions, yet few studies have quantified bacterial growth and death dynamics in the hyporheic zone, and none have systematically analyzed their response to variations in hydraulic or chemical conditions. Here we used MODFLOW and SEAM3D to simulate hydraulics; dissolved oxygen (DO) and dissolved organic carbon (DOC) transport; and aerobic microbial metabolism, growth, and death in hyporheic zones induced by riverbed dunes. We ran simulations both with and without growth/death processes, and varied hydraulic parameters and DO/DOC boundary concentrations. Microbial biomass reached steady state (t = 3 days) in every simulation, at which time there was greater biomass and DOC biodegradation rates in the hyporheic flowcell (300% and 85% higher for the base case, respectively) when accounting for microbial growth dynamics. This occurred as microbial biomass tailored its spatial distribution to the availability of DO and DOC, demonstrating the importance of simulating growth/death processes. Biomass generally increased with hyporheic flow cell area as upwelling groundwater decreased. When varying surface water DO and DOC source concentrations relative to the base case, the greatest effect on biomass occurred when increasing DOC and decreasing DO. We determined minimum DO and DOC steady-state concentrations required for microbial growth, but the minimums were not absolute or related by stoichiometry. Increasing DOC created a smaller area of microbes with higher concentrations relative to the base case. Increasing DO slightly increased the area occupied by microbes while keeping the total biomass nearly constant. Overall, microbial growth and death dynamics depend on DO and DOC availability in the hyporheic zone, which is dependent on DOC/DO boundary concentrations and hyporheic flow paths, and in turn the hydraulic interaction between surface water and groundwater.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparing Adjoint Waveform Tomography Models of California Using Different Starting Models

Abstract Adjoint waveform tomography (AWT) sits at the cutting edge of seismic tomography on local, regional, and global scales. However, the choice in starting model may have a significant impact on the final inversion results. In this paper, we present 3 AWT models of California that are based on different starting models. We chose three models that were inverted at different scales: SPiRaL, a global travel‐time tomography model (Simmons et al., 2021, 10.1093/gji/ggab277 ), CSEM_NA, a regional adjoint tomography model of North America and the North Atlantic (Krischer et al., 2018, 10.1029/2017JB015289 ), and WUS256, a regional adjoint tomography model of the western US (Rodgers et al., 2022, https://doi.org/10.1029/2022JB024549 ). We then inverted three AWT models using the same source and receiver set. We ran each model over three period bands: 30–100 s, 25–100 s, and 20–80 s. Once the iterations were finalized, we used five methods of testing model similarity in both the model and data space. We conclude that the choice of starting model has a minimal impact on long wavelength models if an appropriate multi‐scale inversion approach is used.

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 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 OSTI2026

Fluid‐Driven Cohesive Strengthening: Critical Role of Reaction Kinetics as the Determinant for Frictional Stability

After an earthquake, faults in the Earth's crust lock up and slowly regain strength. This process, called fault healing, plays a key role in determining how often earthquakes happen and how large they can be. While scientists know that healing involves both mechanical friction and chemical changes in rocks, it has been unclear how these processes work together, especially over shorter timescales in laboratory settings. In our study, we ran controlled experiments on the mineral anhydrite under different fluid conditions. We found that the presence and pressure of water speeds up fault healing by triggering a chemical reaction that strengthens the fault over time. This strengthening makes the fault more resistant to sliding, which could delay future earthquakes—but also sets the stage for larger ones when they do occur. We propose a model that combines frictional and chemical processes to explain this water‐activated, time‐dependent healing. Importantly, our results show that even faults typically considered stable could still produce earthquakes if they undergo sufficient chemical healing. This has broad implications for understanding where and when earthquakes might strike, especially in fluid‐rich environments such as subduction zones and geoenergy settings.

Affinito, R [Pennsylvania State Univ., University ↗
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 OSTI2026

When Do Riverine Systems “Feel the Burn”? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export

Wildfires impact terrestrial landscapes and downstream river corridors through shifts in vegetation and soil properties leading to downstream hydrologic and water quality impacts. The magnitude of these impacts depend on a complex and interconnected set of wildfire, landscape, and aquatic processes. Here, we isolate the impact of post-fire hydrologic changes on streamflow, nitrate, and dissolved organic carbon using the Soil and Water Assessment Tool (SWAT) model. We explore how responses differ across burn severity and area burned in two test basins: a humid forested basin and a semi-arid mixed land use basin. We ran 1830 wildfire simulations testing impacts of area burned, burn severity, and post-fire precipitation on streamflow, nitrate, and dissolved organic carbon. Our work suggests that area burned thresholds differ with burn severity and analyte. Additionally, post-fire transport of dissolved organic carbon was sensitive to both area burned and severity, while nitrate was primarily sensitive to area burned. Despite a muted (−9.5 to 5.7 mm yr −1 change) hydrologic response in the semi-arid basin, the model predicted large (7%–288% increase) shifts in dissolved organic carbon, suggesting that post-fire shifts in flow pathways and soil properties are key in its response. The limited shifts in nitrate responses in the simulations highlight that terrestrial post-fire transformations, rather than hydrologic changes, may control the increases in stream nitrate often observed post-fire. As wildfire regimes are shifting, improving understanding of post-fire nutrient export responses is critical to protect freshwater resources and aquatic ecosystems.

Wampler, Katherine A. [Pacific Northwest National ↗
DOE OSTI2022

Detecting damaged buildings using real-time crowdsourced images and transfer learning

After significant earthquakes, we can see images posted on social media platforms by individuals and media agencies owing to the mass usage of smartphones these days. These images can be utilized to provide information about the shaking damage in the earthquake region both to the public and research community, and potentially to guide rescue work. This paper presents an automated way to extract the damaged buildings images after earthquakes from social media platforms such as Twitter and thus identify the particular user posts containing such images. Using transfer learning and ~ 6500 manually labelled images, we trained a deep learning model to recognize images with damaged buildings in the scene. The trained model achieved good performance when tested on newly acquired images of earthquakes at different locations and when ran in near real-time on Twitter feed after the 2020 M7.0 earthquake in Turkey. Furthermore, to better understand how the model makes decisions, we also implemented the Grad-CAM method to visualize the important regions on the images that facilitate the decision.

58 GEOSCIENCES↗
DOE OSTI2021

Engineering block co-polymer anion exchange membrane domains for highly efficient electrode-decoupled redox flow batteries

We report, anion exchange membranes (AEMs) with high permselectivity, chemical stability and mechanical robustness are a key enabling technology for a variety of electrochemical energy conversion and storage systems such as fuel cells, redox flow batteries and water electrolyzers. Herein we examine well understood chemically stable triblock co-polymer AEMs derived by the chloromethylation (CM-) and subsequent trimethylamine (TMA) functionalization of polystyrene-block-poly(ethylene-ran-butylene)-block-polystyrene (SEBS) block co-polymers. We show using small angle X-ray scattering (SAXS) that the introduction of a hydrophobic reinforcement matrix disrupts phase segregation and decreases ionic domain sizes in these AEM separators. The reinforcement matrix simultaneously causes the ultimate tensile strength (UTS) of the reinforced separator to increase by almost 600% compared to the unreinforced form, while maintaining comparable elasticity (ca. 500% elongation at break). Such reinforced r-SEBS30-TMA separators enable electrode-decoupled redox flow batteries (ED-RFBs) with significantly improved lifetimes by decreasing per cycle capacity fade from 0.5% to <0.05%. This 10× improvement in operando permselectivity opens the door to ED-RFBs employing inexpensive elemental actives achieving decades-long usable lifespans.

25 ENERGY STORAGE↗