Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetically constrained models”

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.

178 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Pickup Ion Phase Space Distributions at Titan in a Three Dimensional Exosphere

The composition and structure of neutral exospheres imbedded in moving plasmas can be determined by measurements of the velocity distributions of their pickup ion progeny. In turn, the velocity distributions are dependent on the spatial structure of the neutral source gases. Since Titan's neutral exosphere extends into the Saturn's magnetosphere (or solar wind) and well above its ionopause, it serves as a good place to analyze such characteristics. They are analyzed using pickup ion measurements made by the Cassini Plasma Spectrometer (CAPS) at Titan [e.g., Hartle et al., 2006] and an ion kinetic model. An early version of the model [Hartle and Sittler, 2007] is an expression describing the phase space density of pickup ions, which is derived from the Vlasov equation with an ion source that explicitly accounts for the velocity and spatial variation of the exosphere source gases. The current version used here includes exosphere source gases in three dimensions. A fundamental parameter of the phase space densities is the ratio of the gyroradius to the neutral scale height alpha, = r(sub g)/H. Titan's exosphere structure yields pickup ions whose phase space distributions are beam-like when alpha >> 1 and fluid-like when alpha << 1. Downstream from the source peak, the light pickup ions, with alpha << 1, are easily observed because their phase space densities are almost uniform over the orbit phases. On the other hand, the phase space distributions of the heavier ions, with alpha >> 1, peak over narrow velocity and spatial ranges. This beam-like nature makes it considerably more difficult to observe heavy ions because their downstream positions and viewing directions are narrowly constrained. Examples of these extremes will be discussed.

Hartle, Richard↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

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

Characterization of micro-sandwich structures via direct ink writing epoxy based cores

Sandwich structured (SS) composites demonstrate considerable flexural stiffness and high strength-to-weight ratios and can be tailored as functional materials. Historically they have been constrained to specific material types and geometry due to limitations in manufacturing methods. However, employing additive manufacturing (AM), specifically direct ink writing (DIW), can provide an alternative method for making SS composites with complex and controllable micro and mesostructures with multifunctionality targeted at desired mechanical, thermal, and electrical properties. DIW, an extrusion-based AM technique, uses a viscous and thixotropic ink with desired components that, once printed, is cured to obtain the final complex net shape parts. In this paper, a novel hybrid AM technique is employed to manufacture SS composite materials containing bisphenol A-based epoxy core and carbon fiber reinforced polymer (CFRP) face sheets that are fabricated via DIW and vacuum infusion process (VIP), respectfully. We demonstrate that the fabrication of these SS composites can be tailored from a thermosetting material, from which additives and/or various lattice structures can be manufactured to achieve enhanced and desirable mechanical integrity with functional properties. Surface topology and mechanical testing techniques are used to characterize the fabricated hybrid SS composites to study and assess mechanical stability. A rheo-kinetic cure model was developed for the core material to allow for additive manufacturing process requirements while ensuring complete cross-linking for the thermoset-based core material. Because of the ability to obtain relatively small core-thickness and controlled architecture, this method now allows for fabricating layered micro-sandwich structures for realizing further light-weighting in relevant applications.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
DOE OSTI2021

Simulations of Multi-Mode Combustion Regimes Realizable In a Gasoline Direct Injection Engine

Lean and dilute gasoline compression ignition (GCI) operation in spark ignition (SI) engines are an attractive strategy to attain high fuel efficiency and low NOx levels. However, this combustion mode is often limited to low-load engine conditions due to the challenges associated with autoignition controllability. In order to overcome this constrain, multi-mode (MM) operating strategies, consisting of advanced compression ignition (ACI) at low load and conventional SI at high load, have been proposed. In this three-dimensional computational fluid dynamics study, the concept of multi-mode combustion using two RON98 gasoline fuel blends (Co-Optima Alkylate and E30) in a gasoline direct injection (GDI) engine were explored. To this end, a new reduced mechanism for simulating the kinetics of E30 fuel blend is introduced in this study. To cover the varying engine load demands for multi-mode engines, primary combustion dynamics observed in ACI and SI combustion modes was characterized and validated against experimental measurements. In order to implement part-load conditions, a strategy of mode transition between SI and ACI combustion (i.e., mixed-mode combustion) was then explored numerically by creating a virtual test condition. The results obtained from the mixed-mode simulations highlight an important feature that deflagrative flame propagation regime coexists with ignition-assisted end-gas autoignition. This study also identifies a role of turbulent flow property adjacent to premixed flame front in characterizing the mixed-mode combustion. The employed hybrid combustion model was verified to perform simulations aiming at suitable range of multi-mode engine operations.

Kim, Sayop↗
DOE OSTI2025

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2022

Reversible Ligand Protonation in Noninnocent Constrained-Geometry-Like Group 4 Complexes

The proligands C 5 H 5 CMe 2 CHPhOx R (R = Me 2 , CHMe 2 , and CMe 3 ) react with M(NMe 2 ) 4 (M = Ti, Zr, and Hf) to give monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 , doubly deprotonated (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 , or a mixture of both. The observed products depend on reaction conditions, the oxazoline substituent, and the metal center, with 4-t-butyl-oxazoline or Ti giving only the monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 . Amine elimination from (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 to give (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 is reversible for 4,4-dimethyl- and 4-isopropyl-oxazoline-based ligands in Zr or Hf complexes. Temperature-dependent kinetic studies of the equilibration of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 , HNMe 2 and (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 provide the experimental thermodynamic parameters ΔS° = -17.4 ± 2.6 cal·mol –1 K –1 and ΔH° = -6.8 ± 0.8 kcal·mol –1 . An Eyring plot of the rate constants, determined from the system as it approaches equilibrium, gives the activation entropy and activation enthalpy for the addition of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 and HNMe 2 as -36 ± 4 cal·mol –1 K –1 and 8.9 ± 1.3 kcal·mol –1 , respectively; the elimination of HNMe 2 from (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 is characterized by ΔS ‡ = -19 ± 5 cal·mol –1 K –1 and ΔH ‡ = 15.7 ± 1.5 kcal·mol –1 . DFT computational models indicate a single-step, nonlinear transfer of the H between the benzylic position of the noninnocent, oxazoline-coordinated ligand and NMe 2 . Computations also confirm the negative activation entropy and the trends in the barriers support the experimental results. Together, these studies indicate the importance of steric effects from the oxazoline ligand, metal center, ancillary ligands, and leaving group on the shuttling of the proton between HNMe 2 and the noninnocent ligand. Finally, these effects suggest that coordination of oxazoline to the metal center is a key part of the benzylic deprotonation and noninnocent behavior of the cyclopentadienyl-oxazoline ligand.

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

Precision Kinematic Sunyaev--Zel'dovich Measurements Across Halo Mass and Redshift with DESI DR2 and ACT DR6: Part I. Luminous Red Galaxies

We present the most precise measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect around luminous red galaxies to date, detecting the signal at $18σ$ significance in both harmonic and configuration space. Our analysis cross-correlates 2.4 million spectroscopic LRGs from the Dark Energy Spectroscopic Instrument (DESI) DR2 sample with Data Release 6 (DR6) of the Atacama Cosmology Telescope (ACT). We develop a novel harmonic-space cross-correlation approach using momentum-weighted kSZ templates, yielding nearly uncorrelated bandpowers within a framework consistent with other large-scale structure analyses. By incorporating the LRG halo occupation distribution (HOD) and its uncertainty, we convert measured galaxy gas profiles into halo gas profiles and provide generalized Navarro-Frenk-White (GNFW) fitting profiles, providing empirical targets for tuning feedback efficiency in hydrodynamical simulations and for baryonic modeling in large-scale structure analyses. We find strong evidence that gas profiles do not trace dark matter, providing direct evidence for gas redistribution beyond gravitational collapse. Comparing to hydrodynamical simulations, our measurements favor feedback efficiencies exceeding those in the Battaglia profile, suggesting more efficient gas ejection in group-scale halos than previously predicted. Splitting by redshift, we detect the kSZ signal at SNR $\approx 5$--$10$ in each of four bins and find amplitude evolution consistent with the expected decline in mean halo mass at fixed comoving number density. Splitting by stellar mass, we study the scaling of kSZ amplitude with galaxy properties. Together with BGS and ELG measurements in Paper II, these results span $0.1 \lesssim z \lesssim 1.6$ across three galaxy populations, demonstrating the potential of spectroscopic kSZ to map circumgalactic gas and constrain baryonic feedback.

Qu, F. J. [KIPAC, Menlo Park; Stanford U., HEPL; S↗
NASA NTRS2011

The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectrum

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg(exp 2) with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500 < l < 10,000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. At l = 3000, about half the power at 148 GHz comes from primary CMB after masking bright radio sources. The power from thermal and kinetic SZ is estimated to be Beta(sub 3000) is identical to 6.8 +/- 2.9 mu K (exp 2), where Beta (sub l) is identical to l(l + 1) C(sub l)/2pi. The IR Poisson power at 148 GHz is Bewta(sub 3000) 7.8 +/- 0.7 muK(exp 2) (C(sub l) = 5.5 +/- 0.5 nK(exp 2)), and a clustered IR component is required with Beta (sub 3000) = 4.6 +/- 0.9 muK(exp 2), assuming an analytic model for its power spectrum shape. At 218 GHz only about 15% of the power, approximately 27 mu K(exp 2), is CMB anisotropy at l = 3000. The remaining 85% is attributed to IR sources (approximately 50% Poisson and 35% clustered), with spectral index alpha = 3.69 +/- 0.14 for flux scaling as S(nu) varies as nu(sup alpha). We estimate primary cosmological parameters from the less contaminated 148 GHz spectrum, marginalizing over SZ and source power. The ACDM cosmological model is a good fit to the data (chi square/dof = 29/46), and ACDM parameters estimated from ACT+Wilkinson Microwave Anisotropy Probe (WMAP) are consistent with the seven-year WMAP limits, with scale invariant n(sub s) = 1 excluded at 99.7% confidence level (CL) (3 sigma). A model with no CMB lensing is disfavored at 2.8 sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6 sigma detection of primordial helium, with Y(sub p) = 0.313 +/- 0.044, and a 4 sigma detection of relativistic species, assumed to be neutrinos, with N(sub eff) = 5.3 +/- 1.3 (4.6 +/- 0.8 with BAO+H(sub 0) data). From the CMB alone the running of the spectral index is constrained to be d(sub s) / d ln k = -0,034 +/- 0,018, the limit on the tensor-to-scalar ratio is r < 0,25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension G(sub mu) < 1.6 x 10(exp -7) (95% CL),

Dunkley, J.↗
DOE OSTI2020

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

The Fast Evolution of SN 2010bh Associated with XRF 100316D

The first observational evidence of a connection between supernovae (SNe) and gamma-ray bursts (GRBs) was found about a decade ago. Since then, only half a dozen spectroscopically confirmed associations have been discovered and XRF 1003160 associated with the type-Ic SN 20lObh is among the latest. Aims. We constrain the progenitor radius, the host-galaxy extinction, and the physical parameters of the explosion of XRF l00316D/SN 20lObh at z 0.059. We study the SN brightness and colours in the context of GRB-SNe. Methods. We began observations with the Gamma-Ray burst Optical and Near-infrared Detector (GROND) 12 hours after the GRB trigger and continued until 80 days after the burst. GROND provided excellent photometric data in six filter bands covering a wavelength range from approximately 350 to 1800 nm, significantly expanding the pre-existing data set for this event. Combining GROND and Swift data, the early broad-band spectral energy distribution (SED) is modelled with a blackbody and afterglow component attenuated by dust and gas absorption. The temperature and radius evolution of the thermal component are analysed and combined with earlier measurements available from the literature. Templates of SN 1998bw are fitted to the SN itself to directly compare the lightcurve properties. Finally, a two-component parametrized model is fitted to the quasi-bolometric light curve. which delivers physical parameters of the explosion. Results. The best-fit models to the broad-band SEDs imply moderate reddening along the line of sight through the host galaxy (A(sub v.host = 1.2 +/- 0.1 mag). Furthermore, the parameters of the blackbody component reveal a cooling envelope at an apparent initial radius of 7 x 10(exp 11) cm, which is compatible with a dense wind surrounding a Wolf-Rayet star. A multicolor comparison shows that SN 2010bh is 60 - 70% as bright as SN 1998bw. Reaching maximum brightness at 8 - 9 days after the burst in the blue bands, SN 20lObh proves to be the most rapidly evolving GRB-SNe to date. Modelling of the quasi-bolometric light curve yields M(sub Ni) = 0.21 +/- 0.03 solar M and M(sub ej) = 2.6 +/- 0.2 solar M, typical of values within the GRB-SN population. The kinetic energy is E(sub k) = (2.4 +/- 0.7) x 10(exp 52) erg, which is making this SN the second most energetic GRB-SN after SN 1998bw. Conclusions. This supernova has one of the earliest peaks ever recorded and thereafter fades more rapidly than other GRB-SNe, hypernovae, or typical type-Ic SNe. This could be explained by a thin envelope expanding at very high velocities, which is therefore unable to retain the gamma-rays that would prolong the duration of the SN event.

Olivares E., F.↗
DOE OSTI2020

Nanostructures for Electrical Energy Storage (NEES) (2020 Final Technical Report)

Nanostructures for Electrical Energy Storage (NEES, www.efrc.umd.edu) was an Energy Frontier Research Center supported by the DOE Office of Science, Basic Energy Sciences, from 8/1/2009 to 7/31/2020. Led by the University of Maryland, NEES enjoyed extensive collaborations with its funded partners, including two DOE Laboratories and six universities. The NEES vision has been to reveal a set of scientific insights and design principles that can underpin a next-generation electrical energy storage approach, building on advances in nanoscale science and technology to achieve simultaneous high power and high energy over extended charge/discharge cycling. The vision is motivated by the recognition that scaling into the nano regime opens the door to new physical phenomena and that the tools enlisted in nanoscale research provide major new opportunities for the synthesis not only of materials at molecular scale but for structures at nano scale and above. NEES has translated this vision into its research program based on two observations. First, while the behavior of ions and electrons in electrolytes and in electrode materials is crucial to electrical energy storage (or more appropriately electrochemical energy storage), it is the transport of ion and electron charge between different structural components of a storage device that ultimately determine its performance. With it well recognized that the choice of electrode materials typically constrain ion transport kinetics as well as maximum ion concentration, the search for better electrode materials has been a primary driver of battery research. At the same time the synthesis of electrodes is typically based on aggregation of particles with varying size, shape, and orientation in the electrode. Together with the presence of additional materials to impart electrical conductivity and cohesion to the composite electrode, change in electrode materials is necessarily accompanied by structural changes at the nano/micro scale that are difficult to categorize and manage. From the beginning, NEES’ vision has been to create and study simpler, highly controlled spatial arrangements of known materials as battery components (electrodes, current collectors, and electrolyte) and to understand how design and structure above the molecular scale determines the energy storage performance available from known materials. Second, advances in nanoscience dramatically expanded the portfolio of synthesis methods, structural motifs, and new phenomena available for research. Some of these gave rapid access to new building blocks at the deep nanoscale (e.g., carbon nanotubes grown by self-assembly, nanoscale arrays formed by electrochemical self-alignment, monolayer films controlled by self-limiting reaction). Such advances served as the enabler for the NEES vision to be pursued experimentally through study of 3D structures created and controlled at the nano, micro, and meso scales. Here, we use meso as in the BES MESO Report, implying not only intermediate or varying length scales, but very much the way behavior is influenced by other factors including aggregation of nanocomponents at different densities and spatial configurations, statistical variations in the aggregates, hierarchical architectures in which they can be assembled, or local 3D configurations that result from the architectures. Over its life cycle, NEES has pursued two overarching goals: (1) to understand the scientific fundamentals of electrochemical storage from the nanoscale to the mesoscale; and (2) to create and learn from innovative, controlled, heterogeneous nanostructures, where such nanostructures can enable the first goal and serve as models for future paradigms in energy storage. Specific goals have included: Synthesize heterogeneous nanostructures comprised of multiple materials arranged in controlled fashion and characterize their behavior; Demonstrate and elucidate design principles for achieving simultaneous high power and high energy; Develop materials processes which enable precision control of thin layers and 3D structures; Investigate the impact of artificial interphases on electrode stability during ion insertion/deinsertion; Create dense arrays of nanostructures to understand how the architecture of these assemblies, along with nanostructure design, influences energy storage behavior at the mesoscale; Identify and understand the consequences of nanoconfinement and local inhomogeneities in 3D mesoscale arrays; Develop and apply computational models to stimulate, guide and interpret experiments.

25 ENERGY STORAGE↗
DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗