Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained optimization”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimization of catholyte for halide-based all-solid-state batteries

Halide solid electrolytes gain significant attention due to their high ionic conductivity, low processing temperature, dry air compatibility, and high-voltage stability. However, low cathode active material (CAM) loading in the composite cathode constrains the realization of high energy density for halide-based all-solid-state batteries. In this study, three halide materials, raw Li 3 YBrCl 5 (LYBC-R, <30μm), milled LYBC (LYBC-M, <5 μm) and freeze-dried Li 3 InCl 6 (LIC, <500 nm), were used as catholytes, combined with LYBC-M as the electrolyte and LiIn alloy as the anode. The CAM:catholyte ratio was investigated as well as stack pressure and operating temperature. Our study demonstrates that particle size of the catholyte plays an important role only for high CAM loading or high C-rate cycling. At moderate CAM loading (65 and 70wt% LiNi 0.83 Mn 0.06 Co 0.11 O 2 ) and 0.1 C-rate, all the three catholytes perform well, providing initial discharge capacities >177 mAh/g. At high CAM loading (85wt%) and 0.1 C-rate, a cathode with the nano-scale LIC catholyte provides discharge capacity of 175 mAh/g, while the larger particle size catholytes suffer significantly reduced capacity. Both LYBC and LIC catholytes provided capacity retention >80% after 200 cycles at 0.5C. These results imply that cathode particle size is critically important for performance at high CAM loading. Furthermore, both electrolyte and cathode were tape cast to scale up size and prepare realistic layer thicknesses. A small amount of binder was used in both layers, to balance the electrochemical performance and mechanical properties. Further, the discharge capacity of a tape cell was 152mA h/g at 0.1C with a capacity retention of 81.8% after 20 cycles at 0.5C. The results demonstrate the excellent performance of LYBC as an electrolyte, and provide guidance for halide-based cathode design.

25 ENERGY STORAGE↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scenario optimization for the tokamak ramp-down phase in RAPTOR: Part B. safe termination of DEMO plasmas

An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Deriving cloud droplet number concentration from surface-based remote sensors with an emphasis on lidar measurements

Abstract. Given the importance of constraining cloud droplet number concentrations (Nd) in low-level clouds, we explore two methods for retrieving Nd from surface-based remote sensing that emphasize the information content in lidar measurements. Because Nd is the zeroth moment of the droplet size distribution (DSD), and all remote sensing approaches respond to DSD moments that are at least 2 orders of magnitude greater than the zeroth moment, deriving Nd from remote sensing measurements has significant uncertainty. At minimum, such algorithms require the extrapolation of information from two other measurements that respond to different moments of the DSD. Lidar, for instance, is sensitive to the second moment (cross-sectional area) of the DSD, while other measures from microwave sensors respond to higher-order moments. We develop methods using a simple lidar forward model that demonstrates that the depth to the maximum in lidar-attenuated backscatter (Rmax⁡) is strongly sensitive to Nd when some measure of the liquid water content vertical profile is given or assumed. Knowledge of Rmax⁡ to within 5 m can constrain Nd to within several tens of percent. However, operational lidar networks provide vertical resolutions of > 15 m, making a direct calculation of Nd from Rmax⁡ very uncertain. Therefore, we develop a Bayesian optimal estimation algorithm that brings additional information to the inversion such as lidar-derived extinction and radar reflectivity near the cloud top. This statistical approach provides reasonable characterizations of Nd and effective radius (re) to within approximately a factor of 2 and 30 %, respectively. By comparing surface-derived cloud properties with MODIS satellite and aircraft data collected during the MARCUS and CAPRICORN II campaigns, we demonstrate the utility of the methodology.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

An expert system for choosing the best combination of options in a general-purpose program for automated design synthesis

An expert system was developed to aid a user of the Automated Design Synthesis (ADS) general-purpose optimization computer program in selecting the best combination of strategy, optimizer, and one-dimensional search options for solving a problem. There are approximately 100 such combinations available in ADS. The knowledge base contains over 200 rules, and is divided into three categories: constrained problems, unconstrained problems, and constrained problems treated as unconstrained problems. The inference engine is written in LISP and is available on DEC-VAX and IBM PC/XT computers.

Rogers, J. L., Jr.↗
NASA NTRS1985

Thin film module electrical configuration versus electrical performance

The as made and degraded states of thin film silicon (TFS) based modules have been modelled in terms of series resistance losses. The origins of these losses lie in interface and bulk regions of the devices. When modules degrade under light exposure, increases occur in both the interface and bulk components of the loss based on series resistance. Actual module performance can thus be simulated by use of only one unknown parameter, shunt losses. Use of the simulation to optimize module design indicates that the current design of 25 cells per linear foot is near optimum. Degradation performance suggests a shift to approx. 35 cells to effect maximum output for applications not constrained to 12 volts. Earlier studies of energy based performance and tandem structures should be updated to include stability factors, not only the initial loss factor tested here, but also appropriate annealing factors.

Morel, D. L.↗
NASA NTRS1986

Motion of a drag-free experiment on the Space Station

This paper describes the analysis and results of a study concerning the motion of a particle in a 'drag-free experiment' (zero-gravity) performed on the Space Station. Two cases were considered: (1) motion of a particle with an at rest initial condition and (2) motion of a particle with an initial velocity. In particular, this study examined the possibility of optimal solutions in terms of displacement from the initial location relative to the station center of gravity point. These solutions were of special interest since the experiment is to be performed in a constrained volume. The results of this study indicate that the total displacement may be made smaller if the particle has an initial velocity than if the particle starts from at rest initial conditions.

German, D. J.↗
NASA NTRS2010

Australian Soil Moisture Field Experiments in Support of Soil Moisture Satellite Observations

Large-scale field campaigns provide the critical fink between our understanding retrieval algorithms developed at the point scale, and algorithms suitable for satellite applications at vastly larger pixel scales. Retrievals of land parameters must deal with the substantial sub-pixel heterogeneity that is present in most regions. This is particularly the case for soil moisture remote sensing, because of the long microwave wavelengths (L-band) that are optimal. Yet, airborne L-band imagers have generally been large, heavy, and required heavy-lift aircraft resources that are expensive and difficult to schedule. Indeed, US soil moisture campaigns, have been constrained by these factors, and European campaigns have used non-imagers due to instrument and aircraft size constraints. Despite these factors, these campaigns established that large-scale soil moisture remote sensing was possible, laying the groundwork for satellite missions. Starting in 2005, a series of airborne field campaigns have been conducted in Australia: to improve our understanding of soil moisture remote sensing at large scales over heterogeneous areas. These field data have been used to test and refine retrieval algorithms for soil moisture satellite missions, and most recently with the launch of the European Space Agency's Soil Moisture Ocean Salinity (SMOS) mission, to provide validation measurements over a multi-pixel area. The campaigns to date have included a preparatory campaign in 2005, two National Airborne Field Experiments (NAFE), (2005 and 2006), two campaigns to the Simpson Desert (2008 and 2009), and one Australian Airborne Cal/val Experiment for SMOS (AACES), just concluded in the austral spring of 2010. The primary airborne sensor for each campaign has been the Polarimetric L-band Microwave Radiometer (PLMR), a 6-beam pushbroom imager that is small enough to be compatible with light aircraft, greatly facilitating the execution of the series of campaigns, and a key to their success. An L-band imaging radar is being added to the complement to provide simultaneous active-passive L-band observations, for algorithm development activities in support of NASA's upcoming Soil Moisture Active Passive (.S"M) mission. This paper will describe the campaigns, their objectives, their datasets, and some of the unique advantages of working with small/light sensors and aircraft. We will also review the main scientific findings, including improvements to the SMOS retrieval algorithm enabled by NAFE observations and the evaluation of the Simpson Desert as a calibration target for L-band satellite missions. Plans for upcoming campaigns will also be discussed.

Kim, Edward↗
NASA NTRS2012

Measuring the Redshift Dependence of The Cosmic Microwave Background Monopole Temperature With Planck Data

We study the capability of Planck data to constrain deviations of the cosmic microwave background (CMB) blackbody temperature from adiabatic evolution using the thermal Sunyaev-Zeldovich anisotropy induced by clusters of galaxies. We consider two types of data sets depending on how the cosmological signal is removed: using a CMB template or using the 217 GHz map. We apply two different statistical estimators, based on the ratio of temperature anisotropies at two different frequencies and on a fit to the spectral variation of the cluster signal with frequency. The ratio method is biased if CMB residuals with amplitude approximately 1 microK or larger are present in the data, while residuals are not so critical for the fit method. To test for systematics, we construct a template from clusters drawn from a hydro-simulation included in the pre-launch Planck Sky Model. We demonstrate that, using a proprietary catalog of X-ray-selected clusters with measured redshifts, electron densities, and X-ray temperatures, we can constrain deviations of adiabatic evolution, measured by the parameter α in the redshift scaling T (z) = T0(1 + z)(sup 1−alpha), with an accuracy of sigma(sub alpha) = 0.011 in the most optimal case and with sigma alpha = 0.018 for a less optimal case. These results represent a factor of 2-3 improvement over similar measurements carried out using quasar spectral lines and a factor 6-20 with respect to earlier results using smaller cluster samples.

cosmic background radiation↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turn-key constrained parameter space exploration for particle accelerators using Bayesian active learning

Abstract Particle accelerators are invaluable discovery engines in the chemical, biological and physical sciences. Characterization of the accelerated beam response to accelerator input parameters is often the first step when conducting accelerator-based experiments. Currently used techniques for characterization, such as grid-like parameter sampling scans, become impractical when extended to higher dimensional input spaces, when complicated measurement constraints are present, or prior information known about the beam response is scarce. Here in this work, we describe an adaptation of the popular Bayesian optimization algorithm, which enables a turn-key exploration of input parameter spaces. Our algorithm replaces the need for parameter scans while minimizing prior information needed about the measurement’s behavior and associated measurement constraints. We experimentally demonstrate that our algorithm autonomously conducts an adaptive, multi-parameter exploration of input parameter space, potentially orders of magnitude faster than conventional grid-like parameter scans, while making highly constrained, single-shot beam phase-space measurements and accounts for costs associated with changing input parameters. In addition to applications in accelerator-based scientific experiments, this algorithm addresses challenges shared by many scientific disciplines, and is thus applicable to autonomously conducting experiments over a broad range of research topics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Intracycle RPM control for vertical axis wind turbines

The wind energy market is currently dominated by horizontal axis wind turbines (HAWTs); however, vertical axis wind turbines (VAWTs) are emerging as a design alternative, especially for deep-water offshore siting due to their low center of gravity, ease of access to drivetrain components, and overall simplicity. Due to the absence of a pitch mechanism in large-scale Darrieus VAWTs, stall control has often been used to manage power and loads. Introducing a pitching mechanism in H-type VAWTs has been studied, but this diminishes the mechanical simplicity advantage, and the use of a pitching mechanism in a large-scale Darrieus-type VAWT is not practical. This work examines an innovative, alternative method to control the rotor dynamics of a large-scale 5 MW VAWT to maximize power while constraining loads without introducing any new or complex mechanical elements. This control strategy is termed intracycle revolution per minute (RPM) control, where the rotational speed of the turbine is allowed to vary in an optimal fashion with the azimuthal location of blades as opposed to typical constant RPM operation. An optimization framework is formulated for an open-loop optimal control problem and solved to maximize power subject to constraints on aerodynamic design loads. Results are presented to demonstrate the benefits and the performance limits of intracycle RPM control for large-scale 5 MW Darrieus VAWTs, namely, (1) power production (quantified in terms of AEP) that can be increased subject to baseline load limits and (2) opportunities to significantly increase AEP or decrease loads via intracycle RPM control that are examined for both two-bladed and three-bladed VAWTs.

17 WIND ENERGY↗
DOE OSTI2026

Arctic Deployment of a Fully Integrated Self-Powered Drifting Buoy Harvesting Wave Energy via a Triboelectric Nanogenerator

The Arctic Ocean remains one of the most poorly sampled regions on Earth, where improved in situ environmental monitoring is vital for advancing oceanographic and atmospheric studies. However, data collection efforts are constrained by the short operational lifespans and high costs of conventional systems. Drifting buoys powered by pendulum-driven wave energy harvesters offer a cost-effective alternative, yet earlier designs have neither been optimized for real-world wave conditions nor validated in the Arctic. In this study, we develop a self-powered drifting buoy that integrates a pendulum-driven triboelectric nanogenerator (TENG) system with a mechanical motion rectifier, a high-gear-ratio transmission, and power management circuits. Through coupled buoy–pendulum dynamic simulations and laboratory testing using a motion simulator, we identify an optimal pendulum mass of 1.6 kg (12.7% of total buoy weight) that maximizes energy output while maintaining buoy stability. Laboratory experiments achieved average power outputs of 12.7 mW under Arctic-like wave and temperature conditions. The system was successfully deployed in the Bering Sea, where it generated 11 J of energy in 3.1 m waves, marking the first Arctic deployment of a TENG-based drifting buoy for sea surface temperature monitoring. This work establishes a cost-effective framework for designing self-powered Arctic monitoring platforms and advances the feasibility of long-term environmental observations in real Arctic waters.

marine enerby↗
DOE OSTI2026

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗
DOE OSTI2022

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters [Slides]

Sensitivity coefficients of prompt neutron decay constant (α) and effective delayed neutron fraction (β eff ) to Pu-239 nuclear data were calculated using a newly available tool (ACEtk). Uncertainty calculations showed the following trends: (1) The prompt neutron decay constant can be used to reduce nuclear data-induced uncertainty in Pu-239(n,f) and (2) the effective delayed neutron fraction can be used to effectively reduce nuclear data-induced uncertainty in Pu-239 total fission $\tilde{v}$. Sensitivity/uncertainty analysis can be used to determine optimal experiment/detector setup to identify compensating errors in nuclear data. Upcoming work includes investigating other nuclide-reaction pairs and response sensitivities, such as Δ$\rho$, neutron leakage spectra, and reaction rate measurements to constrain nuclear data of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Numerical method for design of minimum-drag supersonic wing camber with constraints on pitching moment and surface deformation

A numerical method, based on linearized theory, for designing minimum-drag supersonic wing camber surfaces of arbitrary planform for a given lift, with options for constraining the pitching moment and/or the surface deformation at the trailing edge of the root chord and for selecting any desired combination of eight specified wing-loading distributions to be employed in the optimization procedure is presented. Two examples are given to illustrate applications of the method. The results indicate that relatively small drag penalties are incurred in designing wings to be self-trimming and to have a reasonable camber surface.

Sorrells, R. B.↗
DOE OSTI2022

A Classification of G -invariant Shallow Neural Networks

When trying to fit a deep neural network (DNN) to a G-invariant target function with G a group, it only makes sense to constrain the DNN to be G-invariant as well. However, there can be many different ways to do this, thus raising the problem of “G-invariant neural architecture design”: What is the optimal Ginvariant architecture for a given problem? Before we can consider the optimization problem itself, we must understand the search space, the architectures in it, and how they relate to one another. In this paper, we take a first step towards this goal; we prove a theorem that gives a classification of all G-invariant single-hidden-layer or “shallow” neural network (G-SNN) architectures with ReLU activation for any finite orthogonal group G, and we prove a second theorem that characterizes the inclusion maps or “network morphisms” between the architectures that can be leveraged during neural architecture search (NAS). The proof is based on a correspondence of every G-SNN to a signed permutation representation of G acting on the hidden neurons; the classification is equivalently given in terms of the first cohomology classes of G, thus admitting a topological interpretation. The G-SNN architectures corresponding to nontrivial cohomology classes have, to our knowledge, never been explicitly identified in the literature previously. Using a code implementation, we enumerate the G-SNN architectures for some example groups G and visualize their structure. Lastly, we prove that architectures corresponding to inequivalent cohomology classes coincide in function space only when their weight matrices are zero, and we discuss the implications of this for NAS.

Agrawal, Devanshu↗
DOE OSTI2024

Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction

Abstract Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB polarization B-modes, require exquisite control of instrumental systematics and Galactic foreground contamination. Blind minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds and systematics across diverse sky models without suffering from foreground mismodelling errors. Still, residual foreground contamination from NILC may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios.By adding constraints to NILC to deproject statistical moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has been demonstrated to further enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full automated optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of moments for deprojection, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method.Validated on B-mode simulations of thePICOspace mission concept with four challenging foreground models, ocMILC exhibits lower Galactic foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Astronomy & Astrophysics↗