Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constraining”

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 55 records · Page 3

Constraining Black Carbon Aging in Global Models to Reflect Timescales for Internal Mixing

The radiative effects of black carbon depend critically on its atmospheric lifetime, which is controlled by the rate at which freshly emitted combustion particles become internally mixed with other aerosol components. Global aerosol models strive to represent this process, but the timescale for aerosol mixing is not easily constrained using observations. In this study, we apply a timescale parameterization derived from particle‐resolved simulations to quantify, in a global aerosol model, the timescale for internal mixing. We show that, while highly variable, the average timescale for internal mixing is approximately 3 hr, which is much shorter than the 24‐hr aging timescale traditionally applied in bulk aerosol models. We then use the mixing timescale to constrain the aging criterion in the Modal Aerosol Module. Our analysis reveals that, to best reflect timescales for internal mixing, modal models should assume that particles transition from the hydrophobic (fresh) to the hydrophilic (aged) class once they accumulate a coating thickness equal to four monolayers of sulfuric acid, as opposed to the model's current aging criterion of eight monolayers. We show that, in remote regions like the Arctic and Antarctic, predictions of black carbon loading and its seasonal variation are particularly sensitive to the model representation of aging. By constraining aging in global models to reflect mixing timescales simulated by the particle‐resolved model, we eliminate one of the free parameters governing black carbon's long‐range transport and spatiotemporal distribution.

54 ENVIRONMENTAL SCIENCES↗

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗

Constraining scale dependent growth with redshift surveys

Ongoing and future redshift surveys have the capability to measure the growth rate of large scale structure at the percent level over a broad range of redshifts, tightly constraining cosmological parameters. Beyond general relativity, however, the growth rate in the linear density perturbation regime can be not only redshift dependent but scale dependent, revealing important clues to modified gravity. We demonstrate that a fully model independent approach of binning the gravitational strength G eff (k,z) matches scalar-tensor results for the growth rate fσ 8 (k,z) to 0.02%–0.27% rms accuracy. For data of the quality of the Dark Energy Spectroscopic Instrument (DESI) we find the bin values can be constrained to 1.4%–28%. We also explore the general scalar-tensor form, constraining the amplitude and past and future scalaron mass/shape parameters. Perhaps most interesting is the strong complementarity of low redshift peculiar velocity data with DESI-like redshift space distortion measurements, enabling improvements up to a factor 6–7 on 2D joint confidence contour areas. Finally, we quantify some issues with gravity parametrizations that do not include all the key physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining explosive nucleosynthesis by indirect reaction methods at storage rings using unstable beams in batch mode

Nuclear reaction studies on unstable isotopes can strongly help in improving our understanding of nucleosynthesis in stars. Indirect approaches to determining astrophysical reaction rates are increasingly common-place and undergoing continuous refinement. Of particular interest is the use of such indirect techniques at storage rings, which, among other aspects, allow to recycle rare unstable beams. Here, we propose to investigate the reaction rates of astrophysical interest using indirect methods (surrogate, Trojan horse, etc.) in reverse kinematics at the IMP-CAS storage ring. Long lived radioactive ion beams, produced remotely, can be accelerated, and made to interact with light targets. The proposed reactions are 85 Kr(p, p’γ), 85 Kr(d, pγ), constraining the neutron flux in an s-process branching point, 79 Se(p, p’γ), 79 Se(d, pγ), constraining the temperature in s-process nucleosyntheses, and 59 Fe(d, pγ), constraining core collapse supernovae.

Angelis, G. de [National Inst. of Nuclear Physics ↗

Limits to forests-based mitigation in integrated assessment modelling: global potentials and impacts under constraining factors

Forests-based measures such as afforestation/reforestation (A/R) and reducing deforestation (RDF) are considered promising options to mitigate climate change, yet their mitigation potentials are limited by economic and biophysical factors that are largely uncertain. The range of mitigation potential estimates from integrated assessment models raises concerns about the capacity of land systems to provide realistic, cost-effective and permanent land-based mitigation. We use the Global Change Analysis Model to quantify the economic mitigation potential of forests-based measures by simulating a climate policy including a tax on greenhouse gas emissions from agriculture, forestry, and other land uses. In addition, we assess how constraining unused arable land (UAL) availability, forestland expansion rates, and global bioenergy demand may influence the forests-based mitigation potential by simulating scenarios with alternative combinations of constraints. Results show that the average forests-based mitigation potential in 2020–2050 increases from 738 MtCO 2 .yr -1 through a forestland increase of 86 Mha in the fully constrained scenario to 1394 MtCO 2 .yr -1 through a forestland increase of 146 Mha when all constraints are relaxed. Regional potentials in terms of A/R and RDF differ strongly between scenarios: unconstrained forest expansion rates mostly increase A/R potentials in northern regions (e.g., +120 MtCO 2 .yr -1 in North America); while unconstrained UAL conversion and low bioenergy demand mostly increase RDF potentials in tropical regions (e.g., +76 and +68 MtCO 2 .yr -1 in Southeast Asia, respectively). This study shows that forests-based mitigation is limited by many factors that constrain the rates of land use change across regions. These factors, often overlooked in modelling exercises, should be carefully addressed for understanding the role of forests in global climate mitigation and defining pledges towards the Paris Agreement.

54 ENVIRONMENTAL SCIENCES↗

Constraining the Synthesis of the Lightest 𝑝 Nucleus 74 Se

We provide the first experimental cross section of the 73 As(p, γ) 74 Se reaction to constrain one of the main destruction mechanisms of the p nucleus 74 Se in explosive stellar environments. The measurement was done using a radioactive 73 As beam at incident energies of 3.1 and 3.7 MeV/nucleon. Along with the total cross section measurement, statistical properties of the 74 Se compound nucleus were extracted, constraining the reaction cross section across the entire Gamow window of the γ process. The impact of the experimentally constrained reaction rate on 74 Se production in Type II supernovae is investigated through Monte Carlo one-zone network simulations. The results indicate that the overproduction of 74 Se by Type II supernova models cannot be resolved by nuclear physics and point towards the description of the astrophysical conditions of the γ process.

59 ≤ A ≤ 89↗

Constraining the galaxy-halo connection of infrared-selected unWISE galaxies with galaxy clustering and galaxy-CMB lensing power spectra

We present the first detailed analysis of the connection between galaxies and their dark matter halos for the unWISE galaxy catalog—a full-sky, infrared-selected sample built from WISE data, containing over 500 million galaxies. Using unWISE galaxy-galaxy autocorrelation and Planck CMB lensing-galaxy cross-correlation measurements down to 10 arcmin angular scales, we constrain the halo occupation distribution (HOD), a model describing how central and satellite galaxies are distributed within dark matter halos, for three unWISE galaxy samples at mean redshifts z ¯ ≈ 0.6 , 1.1, and 1.5, assuming a fixed cosmology at the best-fit Planck Λ CDM values. We constrain the characteristic minimum halo mass to host a central galaxy, M min HOD = 1.8 3 - 1.63 + 0.41 × 10 12 M ⊙ / h , 5.22 - 4.80 + 0.34 × 10 12 M ⊙ / h , 6.60 - 1.11 + 0.30 × 10 13 M ⊙ / h and the mass scale at which one satellite galaxy per halo is found, M 1 ' = 1.1 3 - 0.70 + 0.32 × 10 13 M ⊙ / h , 1.18 - 1.11 + 0.30 × 10 13 M ⊙ / h , 1.23 - 1.17 + 0.14 × 10 14 M ⊙ / h for the unWISE samples at z ¯ ≈ 0.6 , 1.1, and 1.5, respectively. We find that all three samples are dominated by central galaxies, rather than satellites. Using our constrained HOD models, we infer the effective linear galaxy bias for each unWISE sample, and find that it does not evolve as steeply with redshift as found in previous perturbation-theory-based analyses of these galaxies. We discuss possible sources of systematic uncertainty in our results, the most significant of which is the uncertainty on the galaxy redshift distribution. Our HOD constraints provide a detailed, quantitative understanding of how the unWISE galaxies populate the underlying dark matter halo distribution. These constraints will have a direct impact on future studies employing the unWISE galaxies as a cosmological and astrophysical probe, including measurements of ionized gas thermodynamics and dark matter profiles via Sunyaev-Zel’dovich and lensing cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutrino-induced charged-current coherent pion production for constraining the muon neutrino flux at DUNE

We study neutrino-induced charged-current coherent pion production ( ν μ CC - Coh π ) as a tool for constraining the neutrino flux at the Deep Underground Neutrino Experiment (DUNE). The neutrino energy and flavor in the process can be directly reconstructed from the outgoing particles, making it especially useful to specifically constrain the muon neutrino component of the total flux. The cross section of this process can be obtained using the Adler relation with the π -Ar elastic scattering cross section taken either from external data or, as we explore, from a simultaneous measurement in the DUNE near detector. We develop a procedure that leverages ν μ CC - Coh π events to fit for the neutrino flux while simultaneously accounting for relevant effects in the cross section. We project that this method has the statistical power to constrain the uncertainty on the normalization of the flux at its peak to a few percent. This study demonstrates the potential utility of a ν μ CC - Coh π flux constraint, though further work will be needed to determine the range of validity and precision of the Adler relation upon which it relies, as well as to measure the π -Ar elastic scattering cross section to the requisite precision. We discuss the experimental and phenomenological developments necessary to unlock the ν μ CC - Coh π process as a “standard candle” for neutrino experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Failure Probability Constrained AC Optimal Power Flow

Despite cascading failures being the central cause of blackouts in power transmission systems, existing operational and planning decisions are made largely by ignoring their underlying cascade potential. This paper posits a reliability-aware AC Optimal Power Flow formulation that seeks to design a dispatch point which has a low operator-specified likelihood of triggering a cascade starting from any single component outage. By exploiting a recently developed analytical model of the probability of component failure, our Failure Probability-constrained ACOPF (FP-ACOPF) utilizes the system's expected first failure time as a smoothly tunable and interpretable signature of cascade risk. Here, we use techniques from bilevel optimization and numerical linear algebra to efficiently formulate and solve the FP-ACOPF using off-the-shelf solvers. Extensive simulations on the IEEE 118-bus case show that, when compared to the unconstrained and N-1 security-constrained ACOPF, our probability-constrained dispatch points can significantly lower the probabilities of long severe cascades and of large demand losses, while incurring only minor increases in total generation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

CHMMPP: A c++ library for constrained Hidden Markov Models

SAND2024-13027O The CHMMPP: A c++ Library for Constrained Hidden Markov Models (HMM) software supports the analysis of multivariate time series data to detect patterns using HMM. Many applications involve the detection and characterization of hidden or latent states in a complex system using observable states and variables. This software supports inference of latent states integrating both an HMM and application-specific constraints that reflect known relationships in hidden states. The CHMMPP software supports application-specific and generic methods for constrained inference. This includes a framework for customized Viterbi methods, constrained inference of hidden states with A* and integer programming methods, and various constraint-informed methods for learning HMM model parameters. CHMMPP focuses on supporting generic methods that enable the agile expression of complex sets of constraints that naturally arise in many real-world applications.

Hart, William↗

arco (Assembled Resource-Constrained Optimization) [SWR-26-030]

Arco (Assembled Resource-Constrained Optimization) is a memory-smart optimization DSL and solver for LP and MIP problems on constrained hardware. The software is an optimization framework built around a KDL-based domain-specific language and a CLI compiler/solver. You write optimization models in .kdl files, and the arco CLI compiles, validates, inspects, and solves them. Language bindings (Python today, more planned) provide programmatic access to the same engine. Built for harder optimization problems on constrained resources, Arco is intentional about every allocation, careful with stack and heap behavior, and relentless about minimizing memory usage so more systems can run real workloads. Arco is built primarily for internal use within our organization. You are welcome to try it, but we make no guarantees about API stability or robustness at this stage

Sanchez Perez, Pedro Andres [National Laboratory o↗

Extratropical Cloud Feedback Constrained by Cloud Sources and Sinks in Cyclones

Constraining cloud feedback in global climate models (GCMs) using observations is important for establishing accurate predictions of future climate. Uncertainty in shortwave cloud feedback (SW FB ) dominates uncertainty in total cloud feedback. Recent studies show a shift toward more positive extratropical SW FB in the latest generations of GCMs leading to the emergence of very high equilibrium climate sensitivity (ECS). In this study, we use precipitation efficiency and albedo susceptibility to constrain liquid water path (LWP) response to warming and SW FB in the Southern Ocean (SO; 50°–80°S). We analyze precipitation in extratropical cyclones (ECs) to learn about extratropical condensed water sink processes, combined with observations of clouds and moisture convergence, and use the analysis to better understand and constrain SW FB . We utilize a perturbed parameter ensemble (PPE) hosted in the Community Atmosphere Model, version 6 (CAM6), to provide a constraint on SW FB based on observations from Clouds and the Earth’s Radiant Energy System (CERES) and Multisensor Advanced Climatology of LWP (MAC-LWP). We apply Gaussian process regression to emulate the model response to all parameters perturbed in the PPE. Confronting the emulator output with observations provides a new estimated response of Earth to global warming. Furthermore, our new estimates of SO LWP reduce the PPE range by 66%–72%, which results in a shortwave cloud radiative effect estimated range that is 27%–34% less than the PPE range. Observations suggest a more positive SO SW FB than the Community Earth System Model, version 2 (CESM2), and consequently do not reject the high climate sensitivity GCMs emerging from the Coupled Model Intercomparison Project phase 6 (CMIP6).

Atmosphere↗

Hyper-Differential Sensitivity Analysis of Uncertain Parameters in PDE-Constrained Optimization

Many problems in engineering and sciences require the solution of large scale optimization constrained by partial differential equations (PDEs). Though PDE-constrained optimization is itself challenging, most applications pose additional complexity, namely, uncertain parameters in the PDEs. Uncertainty quantification (UQ) is necessary to characterize, prioritize, and study the influence of these uncertain parameters. Sensitivity analysis, a classical tool in UQ, is frequently used to study the sensitivity of a model to uncertain parameters. In this article, we introduce "hyper-differential sensitivity analysis" which considers the sensitivity of the solution of a PDE-constrained optimization problem to uncertain parameters. Our approach is a goal-oriented analysis which may be viewed as a tool to complement other UQ methods in the service of decision making and robust design. We formally define hyper-differential sensitivity indices and highlight their relationship to the existing optimization and sensitivity analysis literatures. Assuming the presence of low rank structure in the parameter space, computational efficiency is achieved by leveraging a generalized singular value decomposition in conjunction with a randomized solver which converts the computational bottleneck of the algorithm into an embarrassingly parallel loop. Two multi-physics examples, consisting of nonlinear steady state control and transient linear inversion, demonstrate efficient identification of the uncertain parameters which have the greatest influence on the optimal solution.

97 MATHEMATICS AND COMPUTING↗

Neutrino Induced Charged Current Coherent Pion Production for Constraining the Muon Neutrino Flux at DUNE

We study neutrino induced charge current coherent pion production (νμ CC-Coh π) as a tool for constraining the neutrino flux at the Deep Underground Neutrino Experiment (DUNE). The neutrino energy and flavor in the process can be directly reconstructed from the outgoing particles, making it especially useful to specifically constrain the muon neutrino component of the total flux. The cross section of this process can be obtained using the Adler relation with the π-Ar elastic scattering cross section, taken either from external data or, as we explore, from a simultaneous measurement in the DUNE near detector. We develop a procedure that leverages νμ CC-Coh π events to fit for the neutrino flux while simultaneously accounting for relevant effects in the cross section. We project that this method has the statistical power to constrain the uncertainty on the normalization of the flux at its peak to a few percent. This study demonstrates the potential utility of aνμ CC-Coh π flux constraint, though further work will be needed to determine the range of validity and precision of the Adler relation upon which it relies, as well as to measure the π-Ar elastic scattering cross section to the requisite precision. We discuss the experimental and phenomenological developments necessary to unlock the νμ CC-Coh π process as a "standard candle'' for neutrino experiments.

Putnam, Gray [Fermilab]↗

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗

Process-based and observation-constrained SOA simulations in China: the role of semivolatile and intermediate-volatility organic compounds and OH levels

Organic aerosol (OA) is a major component of tropospheric submicron aerosol that contributes to air pollution and causes adverse effects on human health. Chemical transport models have difficulties in reproducing the variability in OA concentrations in polluted areas, hindering understanding of the OA budget and sources. Herein, we apply both process-based and observation-constrained schemes to simulate OA in GEOS-Chem. Comprehensive data sets of surface OA, OA components, secondary organic aerosol (SOA) precursors, and oxidants were used for model–observation comparisons. The base models generally underestimate the SOA concentrations in China. In the revised schemes, updates were made on the emissions, volatility distributions, and SOA yields of semivolatile and intermediate-volatility organic compounds (SVOCs and IVOCs) and additional nitrous acid sources. With all the model improvements, both the process-based and observation-constrained SOA schemes can reproduce the observed mass concentrations of SOA and show spatial and seasonal consistency with each other. Our best model simulations suggest that anthropogenic SVOCs and IVOCs are the dominant source of SOA, with a contribution of over 50% in most of China, which should be considered for pollution mitigation in the future. The residential sector may be the predominant source of SVOCs and IVOCs in winter, despite large uncertainty remaining in the emissions of IVOCs from the residential sector in northern China. The industry sector is also an important source of IVOCs, especially in summer. More SVOC and IVOC measurements are needed to constrain their emissions. Besides, the results highlight the sensitivity of SOA to hydroxyl radical (OH) levels in winter in polluted environments. The addition of nitrous acid sources can lead to over 30% greater SOA mass concentrations in winter in northern China. It is important to have good OH simulations in air quality models.

54 ENVIRONMENTAL SCIENCES↗

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗

Interpreting Primal-Dual Algorithms for Constrained Multiagent Reinforcement Learning

Constrained multiagent reinforcement learning (C-MARL) is gaining importance as MARL algorithms find new applications in real-world systems ranging from energy systems to drone swarms. Most C-MARL algorithms use a primal-dual approach to enforce constraints through a penalty function added to the reward. In this paper, we study the structural effects of this penalty term on the MARL problem. First, we show that the standard practice of using the constraint function as the penalty leads to a weak notion of safety. However, by making simple modifications to the penalty term, we can enforce meaningful probabilistic (chance and conditional value at risk) constraints. Second, we quantify the effect of the penalty term on the value function, uncovering an improved value estimation procedure. We use these insights to propose a constrained multiagent advantage actor critic (C-MAA2C) algorithm. Simulations in a simple constrained multiagent environment affirm that our reinterpretation of the primal-dual method in terms of probabilistic constraints is effective, and that our proposed value estimate accelerates convergence to a safe joint policy.

chance constraints↗