Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Content Model”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2020

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↗
DOE OSTI2021

Topographical shifts in fine root lifespan in a mixed, mesic temperate forest

Root lifespan, often is estimated in landscape- and ecosystem-level carbon models using linear approximations. In water manipulation experiments, fine root lifespan can vary with soil water content. Soil water content is generally structured by complex topography, which is largely unaccounted for in landscape- and ecosystem-scale carbon models. Topography governs the range of soil water content experienced by roots which may impact their lifespan. We hypothesized that root lifespan varied nonlinearly across a temperate, mesic, forested catchment due to differences in soil water content associated with topographic position. We expected regions of the landscape that were too wet or too dry would have soils that were not optimal for roots and thus result in shorter root lifespans. Specifically, we hypothesized that root lifespan would be longest in areas that consistently had soil water content in the middle of the soil water content spectrum, while in soils at either very low or very high soil water content, root lifespan would be relatively short. We tested this hypothesis by collecting and analyzing two years of minirhizotron and soil moisture data in plots widely distributed in the Shale Hills catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Pennsylvania. We found that fine root lifespans were longer in traditionally wetter topographic regions, but detected no short term (biweekly) effect of soil moisture on root lifespan. Additionally, depth in soil, soil series, slope face orientation, and season of birth strongly affected root lifespans across the catchment. In contrast, lifespan was unaffected by root diameter or mycorrhizal association. Failure to account for these variables could result in erroneous estimates of fine root lifespan and, consequentially, carbon flux in temperate forested regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Finite Element Modeling of Spent Nuclear Fuel in Dry Storage under Seismic Excitation

Spent nuclear fuel (SNF) is currently stored at independent spent fuel storage installations (ISFSIs) at nuclear power plant sites across the US. The US Department of Energy (DOE) Spent Fuel and Waste Science and Technology (SFWST) program is currently working to understand the response of SNF during a seismic event. This investigation is part of a long-term effort by DOE to quantify the external mechanical loads on SNF. The SFWST program has recently completed campaigns evaluating SNF behavior during normal conditions of transport and in 30 cm package drop scenarios; both campaigns included a combination of testing and modeling. Likewise, the investigation into seismic response of SNF will include full-size shake table tests of a vertical storage overpack containing instrumented surrogate fuel assemblies, as well as pre- and post-test finite element modeling. This work focuses on the development of finite element models simulating the planned test configuration. The seismic response is calculated at two scales: at the package level, and at the fuel level. The package level model evaluates the behavior of the storage overpack and its contents in response to the ground motions of a representative earthquake. The motion of the package is input to the fuel level model, which consists of a detailed finite element model of a 17x17 or 16x16 pressurized water reactor (PWR) fuel assembly or of a 10x10 boiling water reactor (BWR) fuel assembly. The stress and strain on individual components such as fuel rods and guide tubes are extracted from this model. This work details the SNF response to various input ground motions, considering representative earthquakes in the Central and Eastern US under various site conditions. Also discussed is the development of a “reduced scope” finite element model of the SNF, which consists of a single fuel rod while retaining the dynamic response of the entire fuel assembly. This reduced scope model addresses the significant computational expense of modeling the SNF under long-duration input ground motions and facilitates simulating additional parameter variations of the model beyond the as-tested configuration. The modeling results provide important insight that will help to inform future testing and modeling efforts by the SFWST program to close the knowledge gap related to external loads on SNF.

Kadooka, Kevin↗
DOE OSTI2022

Persistent, extensive channelized drainage modeled beneath Thwaites Glacier, West Antarctica

Abstract. Subglacial hydrology is a leading control on basal friction and the dynamics of glaciers and ice sheets. At low discharge, subglacial water flows through high-pressure, sheet-like systems that lead to low effective pressures. However, at high discharge, subglacial water melts the overlying ice into localized channels that efficiently remove water from the bed, thereby increasing effective pressure and basal friction. Recent observations suggest channelized subglacial flow exists beneath Thwaites Glacier, yet it remains unclear if stable channelization is feasible in West Antarctica, where surface melting is nonexistent and water at the bed is limited. Here, we use the MPAS-Albany Land Ice model to run a suite of over 130 subglacial hydrology simulations of Thwaites Glacier across a wide range of physical parameter choices to assess the likelihood of channelization. We then narrow our range of viable simulations by comparing modeled water thicknesses to previously observed radar specularity content, which indicates flat, spatially extensive water bodies at the bed. In all of our data-compatible simulations, stable channels reliably form within 100–200 km of the grounding line and reach individual discharge rates of 35–110 m3 s−1 at the ice–ocean boundary. While only one to two channels typically form across the 200 km width of the glacier in our simulations, their high efficiency drains water across the entire lateral extent of the glacier. We posit the large catchment size of Thwaites Glacier, its funnel-like geometry, and high basal melt rates together accumulate enough water to form stable channels. No simulations resembled observed specularity content when channelization was disabled. Our results suggest channelized subglacial hydrology has two consequences for Thwaites Glacier dynamics: (i) amplifying submarine melting of the terminus and ice shelf while (ii) simultaneously raising effective pressure within 100 km of the grounding line and increasing basal friction. The distribution of effective pressure implied from our modeling differs from parameterizations typically used in large-scale ice sheet models, suggesting the development of more process-based parameterizations may be necessary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗
DOE OSTI2021

Spatial and temporal variability of tritium in precipitation within South Africa and it's bearing on hydrological studies

Tritium, the radioactive isotope of hydrogen, has been used to understand groundwater recharge processes for decades. The current variation of tritium in the atmosphere is largely attributed to stratospheric production and fall out rates as well as global circulation phenomena controlling the hydrological cycle. Global controls on the variability in atmospheric tritium activity are poorly suited to explain local variation and tritium activities in precipitation are often assumed to be uniform over both local and regional catchments and watersheds. This assumption can result in both over and under estimation of modern recharge within an aquifer when using tritium as the recharge proxy. In order to minimize the inherent prediction residuals associated with tritium based recharge investigations, the variability of tritium in precipitation was modelled from 127 spatial precipitation samples taken over a two year period, combined with a 76 precipitation sample group-set taken over a one year period in a single location. Precipitation events were traced backward in time, from the point of collection, using HYSPLIT modelling to ascertain the origins of moisture content as well as the altitudes of moisture origin reached along the particle track. Tritium activities, collected over a one year period in Paarl, range from 0.45 to 4.16 TU and have a mean of 1.59 TU. Spatial storm events in the Western Cape in 2017 and 2018 had a range from 0 to 2.2 and 0.37 to 3.27 TU, respectively, with mean activities of 1.18 (n = 34) and 1.25 TU (n = 32). Both storm events had similar tritium variability (σ = 0.5 n = 35 and 0.48 n = 32). Regional precipitation events had the largest range of tritium activities (0.55–12.2 TU). Although not all tritium activities can be explained by interrogating the water mass origin, this study suggests that approximately 90% of events can be completely or partially attributed to the origin of the water mass. The variability of tritium, both spatially and temporally, was higher than expected, confirming that when uniform tritium inputs are used, the groundwater system would provide inaccurate modern recharge estimates. Higher spatial resolution of tritium variation in precipitation for a particular region will improve our ability to relate tritium activities in groundwater to local precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measurement of f1(1285) production in pp collisions at $\sqrt{s}$ = 13 TeV

This study presents the first measurement of the f 1 (1285) resonance using the ALICE detector in inelastic proton–proton collisions at a center-of-mass energy of 13 TeV. The resonance is reconstructed at midrapidity (|y| < 0.5) through the hadronic decay channel f 1 (1285)→ K$^0_S$K ± π ± . Key measurements include the determination of its mass, transverse-momentum integrated yield, and average transverse momentum. Additionally, the ratio of the transverse-momentum integrated yield of f 1 (1285) to pion is compared with calculations from the canonical statistical hadronization model. The model calculation, assuming a zero total strangeness content for f 1 (1285), reproduces the data within 1σ deviation, shedding light on the quark composition of f 1 (1285).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

WACCM6 Projections of Polar Mesospheric Cloud Abundance Over the 21st Century

Polar mesospheric clouds (PMC), or noctilucent clouds, can be observed over high latitudes with the naked eye from the ground or from space near the summer solstice. PMC are considered a direct and sensitive indicator of climate change and have been reported to appear more frequently in recent decades. How PMC will change in the future under the influence of natural variability and anthropogenic forcing is uncertain. In this study, we utilize model output from the Whole Atmosphere Community Climate Model under several shared socioeconomic pathway (SSP) scenarios and input the water vapor, temperature, and pressure information into a 0-d PMC model to project the trend and variation of PMC over the 21st century, and their relationship to future changes of temperature, water vapor, and the solar cycle. The 0-d model calculations indicate that PMC ice water content (IWC) will increase and PMC will extend to lower latitudes under high SSP scenarios. Under these scenarios, more mesospheric water vapor leads to an increased IWC of PMC over the polar region, and colder mesopause temperature leads to more PMC over the mid-latitudes. There is a significant anti-correlation between the solar cycle and PMC IWC over the 21st century, but the anti-correlation is not always significant on the decadal scale. Finally, methane oxidation in the stratosphere and water vapor entering from the troposphere are both responsible for future changes in mesospheric water vapor and thus PMC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantifying uncertainty in physics-based predictions of rare-isotope production cross sections via Bayesian-inspired model averaging across nuclear mass tables

Accurate prediction of fragmentation cross sections is essential for rare-isotope beam production, planning new-isotope searches, and designing experiments to study the most exotic regions of the nuclear chart. However, existing reaction models and phenomenological cross-section parametrizations often exhibit significant deviations over broad regions of mass and charge. In this work, a Bayesian-inspired model-averaging framework is developed to combine abrasion-ablation (AA) calculations based on multiple nuclear mass tables into a single statistically weighted estimate. For the calibrated systems, the model weights are assigned empirically according to the relative quality of fit to measured cross sections, thereby reducing systematic model bias while preserving the underlying physics content of the AA description. The weights are constrained using proton-rich fragmentation data for the 78 Kr and 124 Xe projectiles. The resulting parameter trends are then propagated to the 92 Mo and 144 Sm systems through a controlled scaling procedure. In the present implementation, the excitation-energy prescription is fixed, while the averaging is performed across nuclear-mass inputs; the framework provides both weighted cross sections and associated uncertainty estimates. Applied to proton-rich fragmentation, the present approach provides a practical basis for interpolation and limited extrapolation in regions relevant to rare-isotope production. The resulting predictions are used to assess the production of very proton-rich nuclei, and candidate new isotopes are discussed.

Bayesian methods↗
DOE OSTI2023

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Demonstration of the On-the-Fly Shielding Analysis Method: Spent Fuel and Waste Disposition

This report documents work performed supporting the US Department of Energy (DOE) Office of Nuclear Energy (NE) Spent Fuel and Waste Disposition (SFWD) Integrated Waste Management activities under work breakdown structure element 1.08.02.04.01, “Data and Tools Development, Validation, and Maintenance.” In particular, this report fulfills milestone M3SF-21OR020401016, “Implement on-the-fly dose analysis methodology in UNF-ST&DARDS” within work package SF-21OR02040101, “Commercial SNF Characterization - ORNL.” The Used Nuclear Fuel - Storage, Transportation & Disposal Analysis Resource and Data System (UNFST& DARDS) enables automated dose rate calculations for spent nuclear fuel (SNF) transportation packages and storage casks using a Monte Carlo radiation transport code. The explicit method uses a detailed model of the SNF system and its contents. Therefore, a dose rate calculation is required for each as-loaded transportation package or storage cask because the SNF assemblies within a canister typically have unique irradiation characteristics. An alternate method, referred to as the “on-the-fly” shielding analysis method, has been proposed that requires only a set of Monte Carlo dose rate calculations for each transportation packaging/storage cask design. The results of the Monte Carlo dose rate calculations are independent of the SNF assembly irradiation and decay characteristics. The dose rate values may then be combined with the radiation source strength of the SNF assemblies associated with a particular transportation packaging/storage cask design to determine actual dose rates. This report presents on-the-fly dose rate calculations for a representative SNF storage cask and verification of the on-the-fly dose rate calculation results by comparison with reference dose rate calculations using the explicit Monte Carlo dose rate calculation. The on-the-fly shielding analysis method was implemented in UNF-ST&DARDS. A Python program was developed to process the MAVRIC dose rate results obtained by source particle type, energy group, and fuel geometry region. A Python processor created binary files, which were saved as a special UNF-ST&DARDS library for on-the- fly shielding analyses. UNF-ST&DARDS uses the precalculated on-the-fly binary libraries generated by the Python data processor and directly executes the Python code for on-the-fly dose analysis. This Python code unzips the pre-generated binary files mentioned above, reads the data, and combines them with user-specified sources for dose and uncertainty calculations. The Python programs were verified using Excel calculations and by comparison with the values obtained with the MAVRIC post-processing utilities applied to the 3dmap files. This method can currently be used to determine dose rates for as-loaded HI-STORM FW storage casks. The UNF-ST&DARDS analysis wizard for on-the-fly shielding analysis is described in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Carbon Intensities of Refining Products in Petroleum Refineries with Co-Processed Biofeedstocks

Petroleum refineries increasingly seek to generate fuels with lower carbon intensities (CIs; a measure of life cycle greenhouse gas [GHG] emissions per unit of energy of fuel; well-to-wheel [WTW]) to meet growing demand. Co-processing refers to a process that adds biomass-derived feedstocks to the fossil-based feedstocks of existing petroleum refinery process units. With the use of biofeedstocks, it is expected that co-processed fuels would have lower CIs than their petroleum counterparts without requiring changes in the existing infrastructure for producing, transporting, and using fuels. To quantify the GHG emissions reduction benefits of co-processing, this study uses a linear programming model to simulate petroleum refinery conditions with and without co-processing. The co-processing cases include three biofeedstocks (soy oil, tallow, and used cooking oil or UCO) used as 10 vol.% of the feedstock to a hydrotreater or hydrocracker since these lipid-based feedstocks have favorable properties to be treated in a hydrotreater or hydrocracker. In addition, we considered pyrolysis oil used as 10 vol.% of the feedstock to a fluid catalytic cracking (FCC) unit in the modeled refinery due to its higher oxygen content compared to other lipid-based feedstocks. Life cycle analysis (LCA) using two distinct approaches—process-level energy allocation and a refinery-level marginal approach—has been conducted for each case. The LCA results using process-level allocation show that there are no noticeable changes in emissions or energy use impacts at the facility level. The life cycle GHG emission reductions of co-processing cases are mainly related to the fraction of biogenic carbon embedded in each fuel product. For example, co-processed jet fuels (a mixture of fossil and biogenic fuels) made via hydrotreating or hydrocracking have higher biogenic carbon, which results in jet fuel CI reductions of 3.9%–8.6% compared to the CI of baseline petroleum jet fuels on a WTW basis. However, analysis of co-processed pyrolysis oil in an FCC shows that a higher fraction of biofeedstocks (29%) becomes process emissions (i.e., CO and CO 2 ), mainly due to the oxygenates in pyrolysis oil, and so it generates less renewable fuel than biofeedstocks co-processed via hydrotreating or hydrocracking. Using the refinery-level marginal approach, the changes in energy use and emissions of co-processing cases compared to the petroleum-only baseline case are allocated to the changes in fuel production (assuming renewable fuels). This approach generates the life cycle GHG emission values of co-processed renewable fuels, which are comparable to the CIs of standalone biofuel production pathways. However, as this approach relies on a rough assumption that product yields and emissions from co-processing units on fossil feedstocks remain the same with and without biofeedstock inputs, co-processing cases like FCC pyrolysis oil may generate quite skewed results.

02 PETROLEUM↗
DOE OSTI2020

Quantifying Information without Entropy: Identifying Intermittent Disturbances in Dynamical Systems

A system’s response to disturbances in an internal or external driving signal can be characterized as performing an implicit computation, where the dynamics of the system are a manifestation of its new state holding some memory about those disturbances. Identifying small disturbances in the response signal requires detailed information about the dynamics of the inputs, which can be challenging. This paper presents a new method called the Information Impulse Function (IIF) for detecting and time-localizing small disturbances in system response data. The novelty of IIF is its ability to measure relative information content without using Boltzmann’s equation by modeling signal transmission as a series of dissipative steps. Since a detailed expression of the informational structure in the signal is achieved with IIF, it is ideal for detecting disturbances in the response signal, i.e., the system dynamics. Those findings are based on numerical studies of the topological structure of the dynamics of a nonlinear system due to perturbated driving signals. The IIF is compared to both the Permutation entropy and Shannon entropy to demonstrate its entropy-like relationship with system state and its degree of sensitivity to perturbations in a driving signal.

42 ENGINEERING↗
DOE OSTI2022

Porewater Geochemical Assessment of Seismic Indications for Gas Hydrate Presence and Absence: Mahia Slope, East of New Zealand’s North Island

We compare sediment vertical methane flux off the Mahia Peninsula, on the Hikurangi Margin, east of New Zealand’s North Island, with a combination of geochemical, multichannel seismic and sub-bottom profiler data. Stable carbon isotope data provided an overview of methane contributions to shallow sediment carbon pools. Methane varied considerably in concentration and vertical flux across stations in close proximities. At two Mahia transects, methane profiles correlated well with integrated seismic and TOPAS data for predicting vertical methane migration rates from deep to shallow sediment. However, at our “control site”, where no seismic blanking or indications of vertical gas migration were observed, geochemical data were similar to the two Mahia transect lines. This apparent mismatch between seismic and geochemistry data suggests a potential to underestimate gas hydrate volumes based on standard seismic data interpretations. To accurately assess global gas hydrate deposits, multiple approaches for initial assessment, e.g., seismic data interpretation, heatflow profiling and controlled-source electromagnetics, should be compared to geochemical sediment and porewater profiles. A more thorough data matrix will provide better accuracy in gas hydrate volume for modeling climate change and potential available energy content.

03 NATURAL GAS↗
DOE OSTI2020

Information Content of JWST NIRSpec Transmission Spectra of Warm Neptunes

Warm Neptunes offer a rich opportunity for understanding exo-atmospheric chemistry. With the upcoming James Webb Space Telescope (JWST), there is a need to elucidate the balance between investments in telescope time versus scientific yield. We use the supervised machine-learning method of the random forest to perform an information content (IC) analysis on a 11-parameter model of transmission spectra from the various NIRSpec modes. The three bluest medium-resolution NIRSpec modes (0.7–1.27 μm, 0.97–1.84 μm, 1.66–3.07 μm) are insensitive to the presence of CO. The reddest medium-resolution mode (2.87–5.10 μm) is sensitive to all of the molecules assumed in our model: CO, CO{sub 2}, CH{sub 4}, C{sub 2}H{sub 2}, H{sub 2}O, HCN, and NH{sub 3}. It competes effectively with the three bluest modes on the information encoded on cloud abundance and particle size. It is also competitive with the low-resolution prism mode (0.6–5.3 μm) on the inference of every parameter except for the temperature and ammonia abundance. We recommend astronomers to use the reddest medium-resolution NIRSpec mode for studying the atmospheric chemistry of 800–1200 K warm Neptunes; its corresponding high-resolution counterpart offers diminishing returns. We compare our findings to previous JWST IC analyses that favor the blue orders and suggest that the reliance on chemical equilibrium could lead to biased outcomes if this assumption does not apply. A simple, pressure-independent diagnostic for identifying chemical disequilibrium is proposed based on measuring the abundances of H{sub 2}O, CO, and CO{sub 2}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Reconnection-driven Flares in M87*: Proton–Synchrotron-powered GeV Emission

Magnetic reconnection in current layers that form intermittently in radiatively inefficient accretion flows onto black holes is a promising mechanism for particle acceleration and high-energy emission. It has been recently proposed that such layers, arising during flux eruption events, can power the rapid TeV flares observed from the core of M87. In this scenario, inverse-Compton scattering of soft radiation from the accretion flow by energetic electron–positron pairs produced near the reconnection layer was suggested as the primary emission mechanism. However, detailed calculations show that radiation from pairs alone cannot account for the GeV emission detected by the Fermi observatory. In this work, we combine analytic estimates with 3D radiative particle-in-cell simulations of pair–proton plasmas to show that the GeV emission can be naturally explained by synchrotron radiation from protons accelerated in the current sheet. Although the exact proton content of the layer is uncertain, our model remains robust across a broad range of proton-to-pair number density ratios. While protons are subdominant in number compared to pairs, our simulations demonstrate that they can be accelerated more efficiently, leading to a self-regulated steady state in which protons dominate the energy budget. Ultimately, proton synchrotron emission accounts for approximately 5%–20% of the total dissipation power. The majority is radiated as MeV photons via pair synchrotron emission, with a smaller fraction emitted as TeV photons through inverse-Compton scattering.

Active galactic nuclei↗
DOE OSTI2025

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗