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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optimization with Neural Network Feasibility Surrogates: Formulations and Application to Security-Constrained Optimal Power Flow

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Constraining AGN Torus Sizes with Optical and Mid-infrared Ensemble Structure Functions

We propose a new method to constrain the size of the dusty torus in broad-line active galactic nuclei (AGNs) using optical and mid-infrared (MIR) ensemble structure functions (SFs). Because of the geometric dilution of the torus, the MIR response to optical continuum variations has suppressed variability with respect to the optical that depends on the geometry (e.g., size, orientation, opening angle) of the torus. More extended tori have steeper MIR SFs with respect to the optical SFs. We demonstrate the feasibility of this SF approach using simulated AGN light curves and a geometric torus model. While it is difficult to use SFs to constrain the orientation and opening angle due to the insensitivity of the SF on these parameters, the size of the torus can be well determined. Applying this method to the ensemble SFs measured for 587 SDSS quasars, we measure a torus R–L relation of $\mathrm{log}\,{R}_{\mathrm{eff}}(\mathrm{pc})={0.51}_{-0.04}^{+0.04}\times \mathrm{log}({{L}}_{\mathrm{bol}}/{10}^{46}\,\mathrm{erg}\ {{\rm{s}}}^{-1})-{0.38}_{-0.01}^{+0.01}$ in the WISE W1 band and sizes ~1.4 times larger in the W2 band, which are in good agreement with dust reverberation mapping measurements. Compared with the reverberation mapping technique, the SF method is much less demanding in data quality and can be applied to any optical+MIR light curves for which a lag measurement may not be possible, as long as the variability process and torus structure are stationary. While this SF method does not extract all information contained in the light curves (i.e., the transfer function), it provides an intuitive interpretation for the observed trends of AGN MIR SFs compared with optical SFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Aircraft in-situ measurements from SOCRATES constrain the anthropogenic perturbations of cloud droplet number

Aerosol-cloud interactions (ACI) in warm clouds alter reflected shortwave radiation by influencing cloud microphysical and macrophysical properties. The variable of state controlling ACI is the cloud droplet number concentration (N d ). Here, we examine the perturbations in N d due to anthropogenic aerosols (ΔN d, PD-PI ) using a perturbed parameter ensemble (PPE) hosted in the sixth Community Atmosphere Model (CAM6). Surrogate models are created for the CAM6 PPE outputs and are used to generate 1 million model variants of CAM6 by sampling 45 sources of parameter uncertainty. The range of uncertain physical parameters related to ACI are constrained with observations of aerosol and cloud properties from SOCRATES. The likely range of uncertain parameters and the associated range of ΔN d, PD-PI are more strongly constrained with observations of N d relative to observations of cloud condensation nuclei. We conduct sensitivity tests of how constraints on ΔN d, PD-PI are affected by systematic uncertainties in observations and our limitations in our surrogate models created for CAM6 PPE outputs. Based on this, we provide guidance on the impact of reducing systematic uncertainty in airborne microphysical observations and in surrogate models.

Song, Ci [Univ. of Wyoming, Laramie, WY (United St↗
DOE OSTI2020

Seawater pH reconstruction using boron isotopes in multiple planktonic foraminifera species with different depth habitats and their potential to constrain pH and <i>p</i>CO<sub>2</sub> gradients

Boron isotope systematics of planktonic foraminifera from core-top sediments and culture experiments have been studied to investigate the sensitivity ofδ11B of calcite tests to seawater pH. However, our knowledge of the relationship between δ11B and pH remains incomplete for many taxa. Thus, to expand the potential scope of application of this proxy, we report δ11B data for seven different species of planktonic foraminifera from sediment core tops. We utilize a method for the measurement of small samples of foraminifera and calculate the δ11B-calcite sensitivity to pH for Globigerinoides ruber, Trilobus sacculifer (sacc or without sacc), Orbulina universa, Pulleniatina obliquiloculata, Neogloboquadrina dutertrei, Globorotalia menardii, and Globorotalia tumida, including for unstudied core tops and species. These taxa have diverse ecological preferences and are from sites that span a range of oceanographic regimes, including some that are in regions of air–sea equilibrium and others that are out of equilibrium with the atmosphere. The sensitivity of δ11Bcarbonate to δ11Bborate (e.g., Δδ11Bcarbonate/Δδ11Bborate) in core tops is consistent with previous studies for T. sacculifer and G. ruber and close to unity for N. dutertrei, O. universa, and combined deep-dwelling species. Deep-dwelling species closely follow the core-top calibration fo rO. universa, which is attributed to respiration-driven microenvironments likely caused by light limitation and/or symbiont–host interactions. Our data support the premise that utilizing boron isotope measurements of multiple species within a sediment core can be utilized to constrain vertical profiles of pH andpCO2 at sites spanning different oceanic regimes, thereby constraining changes in vertical pH gradients and yielding insights into the past behavior of the oceanic carbon pumps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The collapse of the Cordilleran–Laurentide ice saddle and early opening of the Mackenzie Valley, Northwest Territories, Canada, constrained by 10 Be exposure dating

Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ~15.8 ka (17.1–14.6 ka), ~1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ~14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ~11.2 m between 16 and 13 ka.

58 GEOSCIENCES↗
DOE OSTI2025

3D Continuous Forcing Dataset from 3D Constrained Variational Analysis at SGP

The continuous 3D large-scale forcing (VARANAL3D) data set derived from 3D constrained variational analysis (3DCVA) extends the conventional constrained variational analysis method by incorporating multiple sub-columns within the analysis domain. This advancement introduces spatial variability into the large-scale forcing fields, thereby enriching the data set’s applicability. The VARANAL3D data set spans from 2004 to 2018 and covers a region of 5˚×4.5˚ domain around the ARM SGP site. The analysis domain is divided into 10×9 sub-columns with 0.5˚ resolution. The 3D large-scale forcing data provides necessary variables to drive and evaluate single-column models (SCM), cloud-resolving models (CRM) ,and large-eddy simulations (LES), as well as information for testing model sensitivity to spatial variability of the large-scale forcing data, facilitating more rigorous testing and refinement of physical processes in SCM/CRM/LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Availability Estimates of Low-Cost, Dispatch-Constrained Electricity (LDE)

This dataset accompanies the report "Potential Availability of and Supply Curves for Low-Cost, Dispatch-Constrained Electricity", which provides initial estimates of the quantity and availability of the potential low-cost, dispatch-constrained electricity (LDE) resource in the United States under scenarios with high variable renewable energy (VRE) penetrations. The analysis also estimates supply curves that can be used in subsequent analysis of the opportunity to use the LDE. The estimated LDE resource could be used for various applications that value low-cost electricity and can operate at reduced capacity factors, including electrolytic hydrogen production and carbon capture. The analysis accompanying this dataset was performed in 2018-2019 to support an analysis of "The Technical and Economic Potential of the H2@Scale Concept within the United States" ( https://www.nlr.gov/docs/fy21osti/77610.pdf ). This dataset includes the quantity of estimated LDE resource for each hour of the year in 2050, aggregated to both the national level and the ReEDS balancing area level. As described in the accompanying report, the ReEDS model was used to project the evolution of the power sector over time. Links to the mapping of ReEDS balancing areas to states and ReEDS balancing area shape files are provided. This dataset also includes the report figure data.

08 HYDROGEN↗
DOE OSTI2021

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Constraining the Charge-Sign and Rigidity-Dependence of Solar Modulation

Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from AMS-02 and BESS Polar II collected between 2007 and 2012, during which the heliospheric magnetic field was in a state of negative polarity. Our model uses measurements of the heliospheric magnetic field and the tilt angle of the heliospheric current sheet to accurately predict the effects of solar modulation as a function of time, charge, and rigidity. By incorporating data from a period of negative polarity, we have been able to robustly observe and constrain the charge-dependent effects of solar modulation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Physics-constrained Deep Recurrent Neural Models of Building Thermal Dynamics

We develop physics-constrained and control-oriented predictive deep learning models for the thermal dynamics of a real-world commercial office building. The proposed method is based on the systematic encoding of physics-based prior knowledge into a structured recurrent neural architecture. Specifically, our model mimics the structure of the building thermal dynamics model and leverages penalty methods to model inequality constraints. Additionally, we use constrained matrix parameterization based on the Perron-Frobenius theorem to bound the eigenvalues of the learned network weights. We interpret the stable eigenvalues as dissipativeness of the learned building thermal model. We demonstrate the effectiveness of the proposed approach on a dataset obtained from an office building with $20$ thermal zones.

Building Energy, structured neural networks↗
DOE OSTI2023

Interpreting Primal-Dual Algorithms for Constrained Multiagent Reinforcement Learning: Preprint

We study multiagent reinforcement learning (MARL) with constraints. This setting is gaining importance as MARL algorithms find new applications in real-world systems ranging from power grids to drone swarms. Most constrained MARL (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 the primal-dual approach on the constraints and value function. First, we show that using the constraint evaluation as the penalty leads to a weak notion of safety, but by making simple modifications to the penalty function, we can enforce meaningful probabilistic safety constraints. Second, we show that the penalty term changes the value function in a way that is easy to model, and demonstrate the consequences of not doing so. We conclude with simulations in a simple constrained multiagent environment to back up the theoretical results.

data-driven control↗
DOE OSTI2025

Compressing Vision Transformers in Geospatial Transfer Learning with Manifold-Constrained Optimization

Deploying geospatial foundation models on resource-constrained edge devices demands compact architectures that maintain high downstream performance. However, their large parameter counts and the accuracy loss often induced by compression limit practical adoption.In this work, we leverage manifold-constrained optimization framework DLRT to compress large vision transformer–based geospatial foundation models during transfer learning. By enforcing structured low-dimensional parameterizations aligned with downstream objectives, this approach achieves strong compression while preserving task-specific accuracy. We show that the method outperforms of-the-shelf low-rank methods as LoRA. Experiments on diverse geospatial benchmarks confirm substantial parameter reduction with minimal accuracy loss, enabling high-performing, on-device geospatial models.

Snyder, Thomas [Yale University]↗
DOE OSTI2020

Moment preserving constrained resampling with applications to particle-in-cell methods

The Moment Preserving Constrained Resampling (MPCR) algorithm for particle resampling is introduced and applied to particle-in-cell (PIC) methods to increase simulation accuracy, reduce compute cost, and/or avoid numerical instabilities. The general algorithm partitions the system space into smaller subsets and resamples the distribution within each subset. Further, the algorithm is designed to conserve any number of particle and grid moments with a high degree of accuracy (i.e. machine accuracy). The effectiveness of MPCR is demonstrated with several numerical tests, including a use-case study in gyrokinetic fusion plasma simulations. Finally, the computational cost of MPCR is negligible compared to the cost of particle evolution in PIC methods, and the tests demonstrate that periodic particle resampling yields a significant improvement in the accuracy and stability of the results.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

PyCDFT: A Python package for constrained density functional theory

In this paper, we present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.

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

aphBO-2GP-3B: a budgeted asynchronous parallel multi-acquisition functions for constrained Bayesian optimization on high-performing computing architecture

High-fidelity complex engineering simulations are often predictive, but also computationally expensive and often require substantial computational efforts. The mitigation of computational burden is usually enabled through parallelism in high-performance cluster (HPC) architecture. Optimization problems associated with these applications is a challenging problem due to the high computational cost of the high-fidelity simulations. In this paper, an asynchronous parallel constrained Bayesian optimization method is proposed to efficiently solve the computationally expensive simulation-based optimization problems on the HPC platform, with a budgeted computational resource, where the maximum number of simulations is a constant. The advantage of this method are three-fold. Firstly, the efficiency of the Bayesian optimization is improved, where multiple input locations are evaluated parallel in an asynchronous manner to accelerate the optimization convergence with respect to physical runtime. This efficiency feature is further improved so that when each of the inputs is finished, another input is queried without waiting for the whole batch to complete. Second, the proposed method can handle both known and unknown constraints. Third, the proposed method samples several acquisition functions based on their rewards using a modified GP-Hedge scheme. The proposed framework is termed aphBO-2GP-3B, which means asynchronous parallel hedge Bayesian optimization with two Gaussian processes and three batches. The numerical performance of the proposed framework aphBO-2GP-3B is comprehensively benchmarked using 16 numerical examples, compared against other 6 parallel Bayesian optimization variants and 1 parallel Monte Carlo as a baseline, and demonstrated using two real-world high-fidelity expensive industrial applications. The first engineering application is based on finite element analysis (FEA) and the second one is based on computational fluid dynamics (CFD) simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

PDE-constrained high-order mesh optimization

Here, we present a novel framework for PDE-constrained r-adaptivity of high-order meshes. The proposed method formulates mesh movement as an optimization problem, with an objective function defined as a convex combination of a mesh quality metric and a measure of the accuracy of the PDE solution obtained via finite element discretization. The proposed formulation achieves optimized, well-defined high-order meshes by integrating mesh quality control, PDE solution accuracy, and robust gradient regularization. We adopt the Target-Matrix Optimization Paradigm to control geometric properties across the mesh, independent of the PDE of interest. To incorporate the accuracy of the PDE solution, we introduce error measures that control the finite element discretization error. The implicit dependence of these error measures on the mesh nodal positions is accurately captured by adjoint sensitivity analysis. Additionally, a convolution-based gradient regularization strategy is used to ensure stable and effective adaptation of high-order meshes. We demonstrate that the proposed framework can improve mesh quality and reduce the error by up to 10 times for the solution of Poisson and linear elasto-static problems. The approach is general with respect to the dimensionality, the order of the mesh, the types of mesh elements, and can be applied to any PDE that admits well-defined adjoint operators.

Computer science↗
DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗
DOE OSTI2021

Gaussian process regression constrained by boundary value problems

We develop a framework for Gaussian processes regression constrained by boundary value problems. The framework may be applied to infer the solution of a well-posed boundary value problem with a known second-order differential operator and boundary conditions, but for which only scattered observations of the source term are available. Scattered observations of the solution may also be used in the regression. The framework combines co-kriging with the linear transformation of a Gaussian process together with the use of kernels given by spectral expansions in eigenfunctions of the boundary value problem. Furthermore, it benefits from a reduced-rank property of covariance matrices. We demonstrate that the resulting framework yields more accurate and stable solution inference as compared to physics-informed Gaussian process regression without boundary condition constraints.

42 ENGINEERING↗