Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributionally robust 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 127 records · Page 7

New nonlinear control algorithms for multiple robot arms

Multiple coordinated robot arms are modeled by considering the arms as closed kinematic chains and as a force-constrained mechanical system working on the same object simultaneously. In both formulations, a novel dynamic control method is discussed. It is based on feedback linearization and simultaneous output decoupling technique. By applying a nonlinear feedback and a nonlinear coordinate transformation, the complicated model of the multiple robot arms in either formulation is converted into a linear and output decoupled system. The linear system control theory and optimal control theory are used to design robust controllers in the task space. The first formulation has the advantage of automatically handling the coordination and load distribution among the robot arms. In the second formulation, it was found that by choosing a general output equation it became possible simultaneously to superimpose the position and velocity error feedback with the force-torque error feedback in the task space.

Tarn, T. J.↗

Ultra-Long Distance BOTDA Sensor System Employing Hybrid Amplification and Advanced Noise Reduction Techniques

Brillouin Optical Time Domain Analysis (BOTDA) sensor systems play a pivotal role in distributed sensing, which enables precise measurements of strain and temperature across extensive fiber lengths. This research offers a comprehensive strategy to extend the sensing range of BOTDA systems beyond tens of kilometers while maintaining high spatial resolutions. Such enhanced sensing is realized through the integration of distributed Raman amplification, inline amplification using erbium-doped fiber amplifiers (EDFA), and advanced noise reduction techniques. Optimizing Raman pump and probe wave profiles ensures the robustness of Brillouin scattering signals, effectively countering attenuation-induced losses. Additionally, enhancing weak Brillouin signals by strategically placing EDFAs at optimal fiber locations keeps sensor sensitivity high. Leveraging inherent redundancy in measured data as a function of frequency and fiber distance, the non-local means (NLM) filter removes noise while preserving essential physical information. In summary, this research has shown a holistic exploration of extending BOTDA’s distance sensing capabilities up to 150 km with spatial resolutions of 8 meters.

Bhatta, Hari↗

Network-Cognizant Time-Coupled Aggregate Flexibility of Distribution Systems Under Uncertainties

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. To exploit this flexibility and to use the capacity potential of aggregate DERs, feasible substation power injection trajectories need to be efficiently characterized. This paper provides an ellipsoidal inner approximation of the set of feasible power injection trajectories at the substation such that for any point in the set, there exists a feasible disaggregation strategy of DERs for any load uncertainty realization. The problem is formulated as one of finding the robust maximum volume ellipsoid inside the flexibility region under uncertainty. Though the problem is NP-hard even in the deterministic case, this paper derives novel approximations of the resulting adaptive robust optimization problem based on optimal second-stage policies. The proposed approach yields less conservative flexibility characterization than existing flexibility region approximation formulations. The efficacy of the proposed method is demonstrated on a realistic distribution feeder.

aggregates↗

Multi-scale planning model for robust urban drought response

Increasingly severe droughts are straining municipal water resources and jeopardizing urban water security, but uncertainty in their duration, frequency, and intensity challenges drought planning and response. We develop the Drought Resilient Interscale Portfolio Planning model (DRIPP) to generate optimal planning responses to urban drought. DRIPP is a generalizable multi-scale framework for optimizing dynamic planning strategies of long-term infrastructure deployment and short-term drought response. It integrates climate and hydrological variability with high-fidelity representations of urban water distribution, available technology options, and demand reduction measures to yield robust and cost-effective water supply portfolios that are location-specific. We apply DRIPP in Santa Barbara, California to assess how least cost water supply portfolios vary under different drought scenarios and identify portfolios that are robust across drought scenarios. In Santa Barbara, we find that drought intensity, not duration or frequency, drives cost increases, reliability risk, and regret of overbuilding infrastructure. Under uncertain drought conditions, a diversified technology portfolio that includes both rapidly deployable, decentralized technologies alongside larger centralized technologies minimizes water supply cost while maintaining high robustness to climate uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗

Microphysical Retrievals from Simultaneous Measurements by Airborne and Ground Radars during OLYMPEX

The OLYMPEX field campaign took place over the Olympic Peninsula of Washington during winter 2015-2016. During the intensive observing period, several aircraft flights obtained multi-frequency airborne radar measurements at X-, Ku-, Ka, and W-band from radars aboard the ER-2 and DC-8 aircraft. In addition, ground radars at S- and X-band performed RHI scans under the aircraft ground tracks. These coincident datasets provide a wealth of complementary information about the hydrometeor particle sizes, shapes, and orientations.In order to synthesize these measurements and test the robustness of scattering models, an optimal estimation retrieval based upon a Hitschfeld-Bordan profiling algorithm has been developed. This algorithm retrieves profiles of the particle size distribution parameters Nw and Dm, and, in the ice phase, relative proportions of aggregate, pristine, and rimed particles, using scattering models with different size-density and size-aspect ratio relationships. Under this formulation, only the pristine particles are horizontally aligned and capable of producing non-zero ZDR and KDP.From the nadir-looking airborne multifrequency radars, we find that aggregates may be readily distinguished from rimed and pristine particles owing to their uniqueness in triple-frequency space. However, rimed and pristine particles occupy a similar region in this space and thus polarimetric measurements greatly enhance their identification in our retrieval framework. With these capabilities we will present analyses of three-dimensional hydrometeor mapping obtained during various OLYMPEX cases. Some features retrieved in these cases include a layer of enhanced aggregation about 2km above the melting layer, hypothesized to be maintained by orographic uplift. This layer is often situated above a layer of denser, aligned particles. Regions of riming and supercooled liquid water beneath generating cells are also identified by our retrieval algorithm. Comparisons to in-situ observations and evaluation of scattering models will be presented.

Munchak, S. Joseph↗

Distributed Inertia Management Under Communication Constraints

This paper presents a framework for distributed inertia management based on consensus algorithm. We propose a methodology to achieve optimal operation of inertia sources such as distributed energy resources (DERs) and synchronous generators in real time. Additionally, we analyze the algorithm under communication constraints and evaluate its robustness under scenarios involving communication time delays and packet losses. The proposed approach is validated via simulations on a 4-node system test feeder, demonstrating its effectiveness and resilience.

Yadav, Ajay [ORNL] (ORCID:000000016111881X)↗

Numerical Modeling & Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (CRADA Final Report)

This project represents a continuation of the collaboration between iProTech and NLR to simulate, optimize and design the iProTech Pitching Inertial Pump (PIP) device. The objectives of this TEAMER project are twofold: 1. Refining the physical characteristics of the existing iProTech PIP WEC-Sim model to enhance the model’s fidelity and include controllable components. Key model enhancements target the inclusion of Coulomb friction, the introduction of a controllable bypass valve, and the replacement of traditional check valves with advanced motorized ones. 2. Exploring traditional and advanced control algorithms. From traditional methods like latching control to cutting-edge reinforcement learning (RL) algorithms, the goal is to ensure the PIP device's adaptability and optimal performance across a range of ocean conditions. NLR is tasked with augmenting the WEC-Sim model and implementing the control algorithms, culminating in performance comparison analyses. iProTech will update their existing 3D models, advise on model improvements, and determine crucial system metrics. WEC-Sim, developed in MATLAB/SIMULINK with Simscape Multibody, is the main piece of software that will be used in this project. Coupled with the MATLAB RL Toolbox, it offers a robust platform for in-depth simulation and optimization of the iProTech PIP device. Building on previous work to explore the PIP design space and optimize its geometry, mass distribution, center of gravity and other key parameters, this project aims to refine iProTech’s existing numerical models and develop effective control algorithms that can seamlessly integrate into their future hardware testing campaigns.

16 TIDAL AND WAVE POWER↗

Monte Carlo Tree Search Methods for the Earth-Observing Satellite Scheduling Problem

This work explores on-board planning for the single spacecraft, multiple ground station Earth-observing satellite scheduling problem through artificial neural network function approximation of state–action value estimates generated by Monte Carlo tree search (MCTS). An extensive hyperparameter search is conducted for MCTS on the basis of performance, safety, and downlink opportunity utilization to determine the best hyperparameter combination for data generation. A hyperparameter search is also conducted on neural network architectures. The learned behavior of each network is explored, and each network architecture’s robustness to orbits and epochs outside of the training distributions is investigated. Furthermore, each algorithm is compared with a genetic algorithm, which serves to provide a baseline for optimality. MCTS is shown to compute near-optimal solutions in comparison to the genetic algorithm. The state–action value networks are shown to match or exceed the performance of MCTS in six orders of magnitude less execution time, showing promise for execution on board spacecraft.

Adam P. Herrmann↗

Optimal Control Allocation for Distributed Electric Propulsion in a Series/Parallel Partial Hybrid Powertrain

The SUbsonic Single Aft eNgine (SUSAN) Electrofan is a NASA concept jet transport aircraft with a 2040 entry-into-service date. It utilizes electrified aircraft propulsion (EAP) to enable propulsive and aerodynamic benefits to reduce fuel usage, emissions, and cost. The powertrain consists of a single thrust producing, boundary layer-ingesting (BLI) turbofan gas turbine engine (GTE) with generators driving a series/parallel partial hybrid EAP system. The architecture includes 16 underwing contrarotating BLI fans, eight on each side, in a mailslot configuration. The 16 fans run on power extracted from the GTE through four 5 MW motor/generators connected to the Low-Pressure Spool, and a single 1 MW motor/generator on the High-Pressure Spool. The distributed fans can be used by the flight control to augment or replace the rudder function. At top of climb, the power extracted from the GTE for the fans is boosted by batteries. The design provides redundancy, and the capacity for boost means that the fans are designed to be able to provide additional thrust when necessary. These features can be leveraged in case of a fan or generator failure. This paper sets up the optimal control problem of setpoint determination for individual fans in the distributed propulsion system, accounting for electrical string efficiencies, saturations, and failures. The solution minimizes power consumption while maintaining thrust and torque on the airframe for maneuvering. Additionally, thrust that would have been lost due to temporary fan speed or power saturation is optimally redistributed to maintain overall desired thrust and torque on the aircraft. The power extraction range constraints derive from the gas turbine engine design and the small amount of variation allowed for the engine to maintain operability. The problem formulation allows the number and location of fan failures for which the thrust and torque can be maintained to be investigated, which has implications for certification. Simulations of a coordinated turn utilizing the distributed electric propulsion for yaw rate control under different failure scenarios demonstrate the robustness of the powertrain design to failures and help define its limitations.

Distributed Electric Propulsion↗

An optimization program based on the method of feasible directions: Theory and users guide

The theory and user instructions for an optimization code based on the method of feasible directions are presented. The code was written for wide distribution and ease of attachment to other simulation software. Although the theory of the method of feasible direction was developed in the 1960's, many considerations are involved in its actual implementation as a computer code. Included in the code are a number of features to improve robustness in optimization. The search direction is obtained by solving a quadratic program using an interior method based on Karmarkar's algorithm. The theory is discussed, focusing on the important and often overlooked role played by the various parameters guiding the iterations within the program. Also discussed is a robust approach for handling infeasible starting points. The code was validated by solving a variety of structural optimization test problems that have known solutions obtained by other optimization codes. It has been observed that this code is accurate and robust: it has solved a variety of problems from different starting points. However, the code is inefficient in that it takes considerable CPU time as compared with certain other available codes. Further work is required to improve its efficiency while retaining its robustness.

Ashok D. Belegundu↗

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗

Robustness-Based Design Optimization Under Data Uncertainty

This paper proposes formulations and algorithms for design optimization under both aleatory (i.e., natural or physical variability) and epistemic uncertainty (i.e., imprecise probabilistic information), from the perspective of system robustness. The proposed formulations deal with epistemic uncertainty arising from both sparse and interval data without any assumption about the probability distributions of the random variables. A decoupled approach is proposed in this paper to un-nest the robustness-based design from the analysis of non-design epistemic variables to achieve computational efficiency. The proposed methods are illustrated for the upper stage design problem of a two-stage-to-orbit (TSTO) vehicle, where the information on the random design inputs are only available as sparse point and/or interval data. As collecting more data reduces uncertainty but increases cost, the effect of sample size on the optimality and robustness of the solution is also studied. A method is developed to determine the optimal sample size for sparse point data that leads to the solutions of the design problem that are least sensitive to variations in the input random variables.

Zaman, Kais↗

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗

CO 2 fertilization of terrestrial photosynthesis inferred from site to global scales

Significance The magnitude of the CO 2 fertilization effect on terrestrial photosynthesis is uncertain because it is not directly observed and is subject to confounding effects of climatic variability. We apply three well-established eco-evolutionary optimality theories of gas exchange and photosynthesis, constraining the main processes of CO 2 fertilization using measurable variables. Using this framework, we provide robust observationally inferred evidence that a strong CO 2 fertilization effect is detectable in globally distributed eddy covariance networks. Applying our method to upscale photosynthesis globally, we find that the magnitude of the CO 2 fertilization effect is comparable to its in situ counterpart but highlight the potential for substantial underestimation of this effect in tropical forests for many reflectance-based satellite photosynthesis products.

59 BASIC BIOLOGICAL SCIENCES↗

Robust bidding strategy for aggregation of distributed prosumers in flexiramp market

Distributed prosumers (DPs) are the grid customers that own energy production/storage assets. Due to the flexibility and fast response of their assets, they can procure ancillary service products (ASP) in the wholesale market. An appealing ASP offered by California ISO in the real-time market (RTM) is flexiramp for which market participants do not submit direct offers, and the compensation is based on their energy opportunity costs. Here in this report, we propose a bidding strategy model for DP aggregator participation in the RTM considering energy and flexiramp. First, we develop a risk-averse optimization to determine the optimal energy and reserve product to trade in day-ahead market while considering proper amounts of flexiramp to trade in the RTM. In the RTM, to obtain optimal amounts of energy and flexiramp, the aggregator must submit hourly multi-level price-quantity energy bids for multiple RTM intervals with 15 min time-steps. On this basis, we propose a robust hourly economic bidding strategy model that determines the optimal energy bids in the RTM. We develop an adjustable robust counterpart of the model to address the RTM energy and flexiramp price uncertainties. The simulation results justify the efficacy of our proposed framework in gaining profits from the wholesale market.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Extreme-scale stochastic optimization and simulation via learning-enhanced decomposition and parallelization (Final Technical Report)

Stochastic optimization and simulation models ubiquitously arise in designing and operating complex service/engineering systems. They can be extreme in scale due to high-dimensional data and decisions, and can also involve decisions made sequentially in response to newly revealed data, both causing significant computational challenge. The objective of this research is to explore a unified framework that integrates machine learning with discrete optimization and risk-averse modeling, to improve the efficiency of decomposition paradigms for stochastic optimization and simulations at extreme scale. The models we consider represent a broad class of complex decision-making problems, where 0-1 or continuous decisions are made before and/or after knowing multiple sources of uncertainties that could be correlated. We will employ machine learning methods to dynamically decide and prioritize computational procedures, including cut generation, branching, and bounding of the optimal objective. Furthermore, the research will shed new lights on the traditional decomposition algorithms for extreme-scale computing. Deliverables of the research include new modeling and computational methods for advancing the state-of-the-art research in optimization and simulation, bringing many relevant risk-averse, data-driven optimization problems in practice within the range of tractability. Examples include distributed computing server scheduling and sensor deployment for monitoring critical infrastructures. Success in this effort will enable progress in solving multiple extreme-scale problems in the complex system design and operations arising from DoE missions in energy, environment, and national security.

24 POWER TRANSMISSION AND DISTRIBUTION↗