Search NASASearch

SEARCH · Search NASA

Results for “constrained estimation”

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.

138 records · Page 3

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING

Reversible Ferroelectric Polarization Modulation of Chiral Molecular Ferroelectrics by Circularly Polarized Light

Abstract The optical modulation of ferroelectric polarization constitutes a transformative, non‐contact strategy for the precise manipulation of ferroelectric properties, heralding advancements in optically stimulated ferroelectric devices. Despite its potential, progress in this domain is constrained by material limitations and the intricate nature of the underlying mechanisms. Recent studies have achieved efficient regulation of ferroelectric polarization and thermal conductivity in chiral ferroelectric thin films through the application of left‐ and right‐handed circularly polarized light (LCP and RCP). Differential absorption of circularly polarized light (CPL) induces nonequilibrium carrier dynamics, generating distinctive interfacial electrostatic fields that enable precise control of ultrathin ferroelectric films. For (R)‐BINOL−DIPASi and (S)‐BINOL−DIPASi (C 26 H 26 O 2 Si), polarization changes surpass 23%, exhibiting opposite response under LCP and RCP excitation. In R chiral films, remnant polarization decreases from 1.05 µC cm − 2 under LCP to 0.85 µC cm − 2 under RCP, whereas in S chiral films, polarization increases from 0.85 µC cm − 2 under LCP to 0.98 µC cm − 2 under RCP. This reversible modulation facilitates reliable switching between ON and OFF states, presenting the potential of chiral ferroelectric materials for flexible, high‐speed integrated photonic sensor technologies.

Chemistry

Radiolytically reworked Archean organic matter in a habitable deep ancient high-temperature brine

Abstract Investigations of abiotic and biotic contributions to dissolved organic carbon (DOC) are required to constrain microbial habitability in continental subsurface fluids. Here we investigate a large (101–283 mg C/L) DOC pool in an ancient (>1Ga), high temperature (45–55 °C), low biomass (10 2 −10 4 cells/mL), and deep (3.2 km) brine from an uranium-enriched South African gold mine. Excitation-emission matrices (EEMs), negative electrospray ionization (–ESI) 21 tesla Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR MS), and amino acid analyses suggest the brine DOC is primarily radiolytically oxidized kerogen-rich shales or reefs, methane and ethane, with trace amounts of C 3 –C 6 hydrocarbons and organic sulfides. δ 2 H and δ 13 C of C 1 –C 3 hydrocarbons are consistent with abiotic origins. These findings suggest water-rock processes control redox and C cycling, helping support a meagre, slow biosphere over geologic time. A radiolytic-driven, habitable brine may signal similar settings are good targets in the search for life beyond Earth.

Science & Technology - Other Topics

Physical thickness characterization of the FRIB production targets

The FRIB heavy-ion accelerator, commissioned in 2022, is a leading facility for producing rare isotope beams (RIBs) and exploring nuclei beyond the limits of stability. These RIBs are produced via reactions between stable primary beams and a graphite target. Approximately 20–40% of the primary beam power is deposited in the target, requiring efficient thermal dissipation. Currently, FRIB operates with a primary beam power of up to 20 kW. To enhance thermal dissipation efficiency, a single-slice rotating graphite target with a diameter of approximately 30 cm is employed. The effective target region is a 1 cm-wide outer rim of the graphite disc. To achieve high RIB production rates, the areal thickness variation must be constrained within 2%. This paper presents physical thickness characterizations of FRIB production targets with various nominal thicknesses, measured using a custom-built non-contact thickness measurement apparatus.

Instrumentation for radioactive beams (fragmentati

Oversizing and Part-Load Problems

Oversizing, the common engineering practice of specifying devices with capacity exceeding the actual load requirement, is a widespread practice across virtually all building technologies end-use categories, including HVAC, electrical systems, lighting, appliances, and plug loads. This practice, driven by factors like design uncertainty, institutional pressures, and risk aversion, results in wasted capital investment, control difficulties, and excessive energy consumption due to inefficient part-load operation. Part-load operation, where devices run below maximum capacity, is the dominant operating mode in most energy systems and presents a complex design challenge. Solutions to match output to load fall into three broad categories: constraining the output, adjusting the device’s internal behavior, and linking output to energy storage or other waste-heat reuse applications. The energy implications of part-load are critical, as efficiency often drops sharply as load decreases across a wide diversity of devices. To quantify the extent of this problem, we derive a dimensionless Part-Load Metric (PLM) based on device efficiency and its frequency distribution of operating hours at various output levels. The PLM quantifies the deviation of a device's actual efficiency from its maximum design efficiency. This metric also serves as a measure of "capital inefficiency," enabling engineers to compare the impact of different part-load solutions and providing a unified framework for evaluating performance across various devices and systems.

Meier, Alan

Coupled Chemical and Mechanical Control of Phase Stability in Lanthanide-Substituted BiVO 4

Doping is widely used to enhance the photoelectrochemical performance of BiVO 4 , yet solubility limits and polymorphic stability constrain compositional tuning. Here, in this study, the role of trivalent cation substitution (Ln = La, Nd, Dy, Ho, Y) on pressure-induced phase transformations in Bi 1–x Ln x VO 4 (x ≤ 0.5) is described. Powder X-ray and neutron diffraction reveal that increasing Ln content stabilizes the tetragonal zircon-type polymorph under ambient conditions, while applied pressures of up to ∼5 GPa promote conversion to the monoclinic fergusonite-type polymorph. In-situ neutron diffraction on Bi 0.8 La 0.2 VO 4 shows a reversible monoclinic to tetragonal transition near 2–3 GPa with a bulk modulus of 147 GPa. The extent of conversion depends strongly on dopant identity, concentration, and synthetic route, with mixed-phase solid-state samples converting more efficiently than phase-pure coprecipitated materials. These results demonstrate pressure as a viable pathway to access metastable, doped BiVO 4 compositions beyond conventional solubility limits.

Sypkes, Kathryn I. [University of Sydney, NSW (Aus

Mechanics and Cooling of the CMS Phase-II Tracker Barrel

The CMS upgrade for the High-Luminosity LHC (HL-LHC) era will allow for an incredibly high rate of data collection, to be used in the next generation of precision measurements and searches for beyond the standard model physics. A key element of this upgrade is the full replacement of the silicon strip tracker to enable level-1 triggering on track level primitives, produced by pixel-strip (PS) and strip-strip (2S) modules, tiled together to form the new detector. Mechanically supporting and cooling these modules is not trivial, as material budget in the tracker volume is highly constrained, and as the modules accrue radiation damage they risk becoming entirely non-functional if the silicon sensors are not sufficiently cooled. To address this, each region of the tracker has a bespoke solution for housing and cooling their modules, mostly depending on carbon-fiber and carbon-foam structures, with dual phase CO2 cooling. For the most central part of the tracker, the Flat Barrel with PS modules (Flat TBPS) covers the pseudorapidity range out to $|\eta|=0.4$, and the mechanical structures are manufactured, assembled, and tested entirely by Fermilab and UC Davis. This presentation will describe the intense design and testing requirements for these components, and how they will ensure strong performance of the tracker through the HL-LHC era.

Cosby, Christopher [Fermilab]

Optimization of the FRIB beam dump: a hybrid genetic algorithm and reinforcement learning approach

The operational envelope of high-power-density systems, such as particle accelerators and advanced nuclear energy systems, is critically constrained by the need to manage extreme thermal loads. To address this, we present a novel hybrid optimization framework combining a genetic algorithm (GA) with a soft actor-critic (SAC) deep reinforcement learning agent. This framework was applied to a practical high-heat-flux problem: redesigning the beam dump at the Facility for Rare Isotope Beams (FRIB) for a power upgrade from 20 kW to 50 kW. The resulting design, validated by three-dimensional conjugate heat transfer simulations, suppresses hazardous hot spots and yields a markedly more uniform temperature distribution. This provides a robust operating margin, increasing the average power-handling capability by 72% relative to the current design, demonstrating the framework’s potential to solve complex thermal management challenges in both accelerator technology and advanced nuclear systems.

Accelerator

Purpose, Design, and Analysis of the NSTX-U Centerstack Casing Lateral Support Shims

The National Spherical Torus Experiment Upgrade (NSTX-U) centerstack casing surrounds the central magnet core of NSTXU. It provides the inner vacuum boundary of the plasma chamber and supports the radially inboard plasma facing components lining the plasma chamber. The casing is supported like a cantilever, fixed at the lower flange by a connection to the TF Bundle through a G-10 ring. It is flexibly supported at the top by bellows, and insulated shims intended to react the largest part of the inventory of the disruption loads at the upper support, thereby protecting the upper and lower connections from excessive loads. A large part of the qualification of the casing and tiles requires addressing and quantifying disruption loads. These have been estimated based on extensions of the behavior of NSTX during its run period roughly 15 years ago. Structural qualification of the Upgrade used these estimates of the loading as input to calculations to demonstrate adequacy for a “most probable” disruption. The estimates are included in a requirements document to guide conservative design of all components loaded by disruptions. However, the changes made for the upgrade made the extrapolation of disruption behavior uncertain. Disruptions include axisymmetric loading from axisymmetric currents induced in the vessel shells. In addition, more complex currents called Halo currents are a result of plasma poloidal currents impinging on the casing. The halo currents add another layer of uncertain behavior and net lateral loading to the casing. As a result, NSTX-U is equipped with a number of instruments intended to assess the nature and severity of disruptions. The casing shims are required for lateral support of the casing, but these have also been converted to load cells to provide another measure of the magnitude of symmetric and asymmetric loading. Some of the logic in the development of the requirements document are presented. Electromagnetic analysis of axisymmetric and halo loading is described. Global simulations of the casing transient dynamic response, and local analysis and testing of the shim/load cell are included. The shim load cells and other instrumentation will be an important guide in planning NSTX-U experiments beyond its next run period expected in late 2026 or early 2027.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE

Photosynthetic capacity is reduced by warming but unaffected by elevated CO2 in seedlings of five boreal tree species

Abstract Increasing atmospheric CO2 concentrations fuel global warming, with boreal regions warming at a faster rate than many other areas. Boreal forests are an important component of the global carbon cycle, yet we have little data on photosynthetic responses of boreal trees to elevated CO2 (EC) and warming. We grew seedlings of 5 widespread North American boreal tree species (from Betula, Larix, Picea, and Pinus) under current (410 ppm) or elevated (750 ppm) CO2 and either ambient (+0 °C) or increased (+4 °C or +8 °C) temperature, then measured photosynthetic traits over a range of leaf temperatures. Our results were generally consistent across species: photosynthetic capacity (maximum rates of Rubisco carboxylation, Vcmax, and electron transport, Jmax) was unaffected by EC but decreased under +8 °C warming. Accordingly, net photosynthesis measured at the growth CO2 concentration (Agrowth) was reduced under warming and increased under EC. The thermal optimum for Agrowth (ToptA) increased by ∼1.8 °C with EC but increased with warming in only two species. In contrast, the activation energies and thermal optima for Vcmax and Jmax, which are used to estimate photosynthesis in Earth System Models, were unaffected by growth environment. There were a few interactions between growth, CO2, and warming. These results suggest increased photosynthesis of widespread boreal tree species under EC may be offset by future reductions in photosynthetic capacity related to warming. We also show that the temperature sensitivities of parameters used to estimate global photosynthesis in large-scale models are generally unaffected by simulated climate change in these species.

Plant Sciences

Automation of Solid Oxyde Electrolyzer Cell (SOEC) & Stack Assembly

The work performed under this agreement before the No/Go decision amounted to the following: 1. Preliminary manufacturing requirements defined by business sensitivity, baseline processes & risk assessment; 2. Process & materials development, in-line gauge exploration, and preliminary automation work used to reduce risk & further refine equipment specifications; 3. “Request for Quote” issued to multiple suppliers for major equipment; 4. Application specific equipment, hardware, and fixturing will be more beneficial to fabricate inhouse; 5.Current proposals & estimates meet cycle time, capital spend, & direct labor; space is on target but requires awareness; 6. Team will continue to explore opportunities to reduce risk (dry time, traceability); 7. Go / No-Go, Purchase Orders, Equipment build & commissioning next All other future task were not completed because during the Go /No-Go decision it was confirmed this project was a "No-Go" and did not proceed past BP1; these included the following: 1. Preliminary manufacturing requirements defined by business sensitivity, baseline processes & risk assessment; 2. Process & materials development, in-line gauge exploration, and preliminary automation work used to reduce risk & further refine equipment specifications; 3. “Request for Quote” issued to multiple suppliers for major equipment; 4. Application specific equipment, hardware, and fixturing will be more beneficial to fabricate inhouse; 5. Current proposals & estimates meet cycle time, capital spend, & direct labor; space is on target but requires awareness; 6. Team will continue to explore opportunities to reduce risk (dry time, traceability) In final preparation of the termination of the Automation of Solid Oxyde Electrolyzer Cell (SOEC) & Stack Assembly, Cummins has purchased no equipment with government funds, and there is no government-owned property in Cummins possession related to this project. Minimal labor and travel were completed and paid. This acts as the Final technical report and concludes our business with DOE on this grant.

36 MATERIALS SCIENCE

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu

Eddy Covariance Theory: A Review

Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.

ecosystem fluxes