Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constraints”

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

Strong constraints on the b → sγ photon polarisation from B0 → K*0e+e-decays

An angular analysis of the B 0 → K* 0 e + e - decay is performed using a data sample corresponding to an integrated luminosity of 9 fb - 1 of pp collisions collected with the LHCb experiment. The analysis is conducted in the very low dielectron mass squared ( q 2 ) interval between 0 . 0008 and 0 . 257 GeV 2 , where the rate is dominated by the B 0 → K * 0 γ transition with a virtual photon. The fraction of longitudinal polarisation of the K * 0 meson, F L , is measured to be F L = (4 . 4 ± 2 . 6 ± 1 . 4)%, where the first uncertainty is statistical and the second systematic. The \( {A}_{\mathrm{T}}^{\mathrm{Re}} \) observable, which is related to the lepton forward-backward asymmetry, is measured to be \( {A}_{\mathrm{T}}^{\mathrm{Re}} \) = - 0 . 06 ± 0 . 08 ± 0 . 02. The \( {A}_{\mathrm{T}}^{(2)} \) and \( {A}_{\mathrm{T}}^{\mathrm{Im}} \) transverse asymmetries, which are sensitive to the virtual photon polarisation, are found to be \( {A}_{\mathrm{T}}^{(2)} \) = 0 . 11 ± 0 . 10 ± 0 . 02 and \( {A}_{\mathrm{T}}^{\mathrm{Im}} \) = 0 . 02 ± 0 . 10 ± 0 . 01. The results are consistent with Standard Model predictions and provide the world’s best constraint on the b → sγ photon polarisation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Positivity and geometric function theory constraints on pion scattering

This paper presents the fascinating correspondence between the geometric function theory and the scattering amplitudes with O(N) global symmetry. A crucial ingredient to show such correspondence is a fully crossing symmetric dispersion relation in the z-variable, rather than the fixed channel dispersion relation. We have written down fully crossing symmetric dispersion relation for O(N) model in z-variable for three independent combinations of isospin amplitudes. We have presented three independent sum rules or locality constraints for the O(N) model arising from the fully crossing symmetric dispersion relations. We have derived three sets of positivity conditions. We have obtained two-sided bounds on Taylor coefficients of physical Pion amplitudes around the crossing symmetric point (for example, π + π – → π 0 π 0 ) applying the positivity conditions and the Bieberbach-Rogosinski inequalities from geometric function theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hydraulic constraints to stomatal conductance in flooded trees

Stomatal closure is a pervasive response among trees exposed to flooded soil. We tested whether this response is caused by reduced hydraulic conductance in the soil-to-leaf hydraulic continuum (k total ), and particularly by reduced root hydraulic conductance (k root ), which has been widely hypothesized. We tracked stomatal conductance at the leaf level (g s ) and canopy scale (G s ) along with physiological conditions in two temperate tree species, Magnolia grandiflora and Quercus virginiana, that were subjected to flood and control conditions in a greenhouse experiment. Flooding reduced g s , G s , k root and k total . Path analysis showed strong support for direct effects of k total on g s and for flood duration on k total , but not k root on k total . A process-based model that accounted for the k total reduction predicted the timeseries of G s in flood and control treatment trees reasonably well (predicted versus observed G s R 2 = 0.80 and 0.51 for M. grandiflora and Q. virginiana, respectively). However, accounting only for k root reduction in flooded trees was insufficient for predicting observed G s reduction. Together, these results suggest that hydraulic constraints were not limited to roots and highlight the need to account for flooding effects on k total when projecting forest ecosystem function using process-based models.

Plant stress↗
DOE OSTI2025

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton↗
DOE OSTI2023

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

U-Pb zircon age constraints on the earliest eruptions of the Deccan Large Igneous Province, Malwa Plateau, India

Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.

58 GEOSCIENCES↗
DOE OSTI2023

Gallium isotopic constraints for the origin of the Earth-Moon system

Comparatively heavy isotopic compositions of moderately volatile elements (MVE) in lunar rocks have been advocated to reflect the loss of light isotopes during devolatilization processes from the Moon. In this study we present new gallium (Ga) isotope data for lunar highland rocks, with a focus on the Ferroan Anorthosite Suite (FAS). These are commonly thought to be direct crystallization products from the late lunar magma ocean (LMO) and should contain the majority of the Ga inventory of the Moon. As such, FAS rocks are crucial for identifying the processes that drove Ga isotope fractionation as well as for inferring the Ga isotopic composition of the bulk Moon. Our data reveal that FAS samples have a range in δ 71 Ga from -0.27 to 0.22‰ and are generally isotopically light in Ga compared to other lunar lithologies, but straddle values typical of terrestrial rocks. Although Ga is defined as an MVE, these Ga isotope variations do not correspond with concentrations of more volatile elements, indicating that Ga isotope variations in the FAS are not primarily controlled by devolatilization processes. Instead, the Ga isotopic compositions of bulk FAS rocks broadly correlate with the composition of plagioclase, with the calcium content of plagioclase decreasing as Ga becomes isotopically heavier. This suggests that fractionation of Ga isotopes in FAS rocks was caused by the preferential incorporation of isotopically light Ga into plagioclase during the later solidification stages of the LMO. The progressive crystallization and extraction of plagioclase forces the residual melt towards increasingly heavier Ga isotope ratios, corroborating similar conclusions derived from correlations between δ 71 and Eu* in the mare basalt suite. Using Ga isotope partitioning calculations, we demonstrate that an isotope fractionation coefficient between plagioclase and coexisting melt of -0.3 to -0.4‰ could explain the observed range of δ 71 values in FAS, mare basalt suite rocks, and KREEP. These calculations allow for a first order estimate of the Ga isotopic composition of the bulk silicate Moon prior to plagioclase fractionation and suggest it was close to the composition of the bulk silicate Earth. This would imply that the Moon did not lose a substantial fraction of its Ga inventory during accretion, consistent with new constraints from Rb-Sr isotope systematics that indicate the Moon's volatile deficit was primarily inherited from Theia. In conjunction with the overlap in non-mass-dependent isotope ratios, these collective observations could be reconciled if Theia and the proto-Earth formed in similar regions of the inner Solar System that were already volatile-depleted.

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

Cost-optimal evaluation of centralized and distributed microgrid topologies considering voltage constraints

Optimal design of hybrid renewable mini-grids requires both economic and power quality considerations. Existing modeling approaches address these considerations via separate or loosely coupled models. Here, we extend REopt—a techno-economic optimization model developed at the National Renewable Energy Laboratory—to consider both within a single model. REopt formulates the design problem as a mixed-integer linear program that solves for a site's optimal technology mix, sizing, and operation to minimize life cycle cost. REopt has traditionally assumed a single node system. In the work presented here, we expand the REopt platform to consider multiple connected nodes with associated voltage constraints. In order to do this, we model power flow using a fixed-point linear approximation method. Additionally, we then use the model to explore design considerations of mini-grids in Sub-Saharan Africa. Specifically, we evaluate under what combinations of transmission line distance and capacity it is technically viable and economically preferable to build multiple isolated mini-grids versus an interconnected, centralized system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Post-extreme-event restoration using linear topological constraints and DER scheduling to enhance distribution system resilience

In this paper, a post-extreme-event restoration (PEER) algorithm is proposed to improve distribution system resilience. Linear topological constraints are proposed to ensure radial topology after N-k contingencies, possibly in multiple islands. The approach is made comprehensive by considering dispatchable distributed energy resources (DERs), non-dispatchable DERs, and demand responses, as well as on-load tap changers (OLTCs) and shunt capacitors. The goal is to minimize the accumulative expense caused by load reduction payment or penalty, as well as DER operation cost. As a result, the overall system will survive longer with higher resilience during an extreme event. To verify the effectiveness of the PEER algorithm, we proposed a resilience evaluation algorithm using Monte Carlo simulation (MCS) with reduced scenarios. This is based on a probabilistic model for generating random scenarios which consider the uncertainty of line faults and solar irradiance. Combined with the proposed PEER algorithm, this reduced-scenario MCS can evaluate the expected energy not served (EENS) which is an essential index for distribution system resilience. Case studies of the IEEE 33-bus and 123-bus test systems validate the proposed algorithm in reducing EENS.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Optimal CO 2 storage management considering safety constraints in multi-stakeholder multi-site GCS projects: A Markov game perspective

Geological carbon storage (GCS) projects could involve a diverse array of stakeholders or players from public, private, and regulatory sectors, each with different objectives and responsibilities. Given the complexity, scale, and long-term nature of GCS operations, determining whether individual stakeholders can independently optimize their interests — or whether collaborative coalition agreements are needed — remains a central question for effective GCS project planning and management. To access large, high-quality storage resources, future GCS deployment may increasingly occur in geologically connected sites, where shared geological features such as pressure space and reservoir pore capacity can lead to competitive behavior among stakeholders. In this work, we propose a paradigm based on Markov games to quantitatively investigate how different coalition structures affect the goals of stakeholders. We frame this multi-stakeholder multi-site problem as a multi-agent reinforcement learning problem with safety constraints. Our approach enables agents to learn optimal strategies while complying with safety regulations. We present an example where multiple operators are injecting CO 2 into their respective project areas in a geologically connected basin. To address the high computational cost of repeated simulations of high fidelity models, a previously developed surrogate model based on the Embed-to-Control (E2C) framework is employed. Our results demonstrate the effectiveness of the proposed framework in addressing optimal management of CO 2 storage when multiple stakeholders with different objectives and goals are involved.

58 GEOSCIENCES↗
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning with knowledge constraints for process optimization of open-air perovskite solar cell manufacturing

Perovskite photovoltaics (PV) have achieved rapid development in the past decade in terms of power conversion efficiency of small-area lab-scale devices; however, successful commercialization still requires further development of low-cost, scalable, and high-throughput manufacturing techniques. One of the critical challenges of developing a new fabrication technique is the high-dimensional parameter space for optimization, but machine learning (ML) can readily be used to accelerate perovskite PV scaling. Herein, we present an ML-guided framework of sequential learning for manufacturing process optimization. We apply our methodology to the Rapid Spray Plasma Processing (RSPP) technique for perovskite thin films in ambient conditions. With a limited experimental budget of screening 100 process conditions, we demonstrated an efficiency improvement to 18.5% as the best-in-our-lab device fabricated by RSPP, and we also experimentally found 10 unique process conditions to produce the top-performing devices of more than 17% efficiency, which is 5 times higher rate of success than the control experiments with pseudo-random Latin hypercube sampling. Our model is enabled by three innovations: (a) flexible knowledge transfer between experimental processes by incorporating data from prior experimental data as a probabilistic constraint; (b) incorporation of both subjective human observations and ML insights when selecting next experiments; (c) adaptive strategy of locating the region of interest using Bayesian optimization first, and then conducting local exploration for high-efficiency devices. Furthermore, in virtual benchmarking, our framework achieves faster improvements with limited experimental budgets than traditional design-of-experiments methods (e.g., one-variable-at-a-time sampling). This framework shows the capability of incorporating researchers’ domain knowledge into the ML-guided optimization loop; therefore, it has the potential to facilitate the wider adoption of ML in scaling to perovskite PV manufacturing.

14 SOLAR ENERGY↗
DOE OSTI2020

Revisiting Bayesian constraints on the transport coefficients of QCD

Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD, through methodical model-to-data comparisons. In this manuscript, we discuss recent developments in Bayesian analyses of heavy ion collision data. Furthermore, we highlight the essential role of closure tests in validating a Bayesian analysis before comparison with measurements. We discuss the role of the emulator that is used as proxy for the multistage theoretical model. We use an ongoing Bayesian analysis of soft hadron measurements by the JETSCAPE Collaboration as context for the discussion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors $G^c_{E,M}(Q^2)$ in the momentum transfer range $0\leq Q^2 \leq 1.4$ $\rm GeV^2$. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment $\mu^c_M=-0.00127(38)_{\rm stat}(5)_{\rm sys}$, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero $G^c_{E}(Q^2)$ indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the $[c(x)-\bar{c}(x)]$ distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Assessing observational constraints on dark energy

Observational constraints on time-varying dark energy (e.g., quintessence) are commonly presented on a w 0 –w a plot that assumes the equation of state of dark energy strictly satisfies w(z) = w 0 + w a z/(1 + z) as a function of the redshift z. Recent observations favor a sector of the w 0 –w a plane in which w 0 > –1 and w 0 + w a < –1, suggesting that the equation of state underwent a transition from violating the null energy condition (NEC) at large z to obeying it at small z. In this paper, we demonstrate that this impression is misleading by showing that simple quintessence models satisfying the NEC for all z predict an observational preference for the same sector. We also find that quintessence models that best fit observational data can predict a value for the dark energy equation of state at present that is significantly different from the best-fit value of w 0 obtained assuming the parameterization above. In addition, the analysis reveals an approximate degeneracy of the w 0 –w a parameterization that explains the eccentricity and orientation of the likelihood contours presented in recent observational studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Probing the limits of statistical neutron capture for the r process: Experimental constraints on 141 Cs nuclear level densities

The r-process abundance peaks, particularly near mass number A ∼ 130, reflect underlying nuclear structure effects such as closed neutron shells, yet modeling the nucleosynthesis in this region remains hindered by uncertain neutron-capture rates. These rates are especially sensitive to nuclear level densities (NLDs) and γ-ray strength functions of neutron-rich nuclei, where experimental data are scarce. We present the first experimental constraint on the NLD of 141 Cs using the β-Oslo method, extending sensitivity to the neutron-rich regime near the N = 82 closed shell. Our data allow for critical calibration of microscopic NLD models and reveal that 141 Cs lies near the limit of statistical model applicability. Using this experimental input, we evaluate radiative neutron-capture rates across neighboring isotones using both Hauser–Feshbach (HF) and High Fidelity Resonance (HFR) models. Our results show order-of-magnitude rate increases for nuclei along the N = 86 line, signaling a transition to resonance-dominated capture in this region. These findings underscore the importance of constraining NLDs to improve r-process reaction network predictions, particularly in environments where the validity of statistical models breaks down.

Nuclear level density↗
DOE OSTI2025

Observation of W + W – γ production in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector and constraints on anomalous quartic gauge-boson couplings

This Letter reports the observation of W + W – γ triboson production in 140 fb –1 of data collected by the ATLAS detector from proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV at the LHC. Events with an opposite-charge eμ pair, a high transverse-momentum photon, and significant missing transverse momentum are considered. The observed (expected) significance of the signal is 5.9 (6.0) standard deviations. The measured fiducial cross-section, defined for the W + W – γ → e ± μ ∓ $ν\bar{ν}$γ final state is 6.2 ± 0.8 (stat.) ± 0.6 (sys.) fb, in good agreement with the Standard Model prediction of 6.1$^{+1.0}_{–0.7}$ fb. Constraints on the Wilson coefficients of 13 dimension-8 operators describing physics beyond the Standard Model through anomalous quartic gauge-boson couplings are derived using the effective field theory framework.

Anomalous quartic gauge couplings↗
DOE OSTI2022

Multidisciplinary Constraints on the Thermal‐Chemical Boundary Between Earth's Core and Mantle

Abstract Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core‐mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co‐constrain core‐mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports Q CMB ∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long‐lived basal mantle molten layer through much of Earth's history.

58 GEOSCIENCES↗