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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Conjugate gradient determination of optimal plane changes for a class of three-impulse transfers between noncoplanar circular orbits

A particular type of three-impulse transfer between two circular orbits is analyzed. The possibility of three plane changes is recognized, and the problem is to optimally distribute these plane changes to minimize the sum of the individual impulses. Numerical difficulties and their solution are discussed. Numerical results obtained from a conjugate gradient technique are presented for both the case where the individual plane changes are unconstrained and for the case where they are constrained. Possibly not unexpectedly, multiple minima are found. The techniques presented could be extended to the finite burn case, but primarily the contents are addressed to preliminary mission design and vehicle sizing.

Burrows, R. R.↗
NASA NTRS2022

Analyzing the Efficacy of Flexible Execution, Replanning, and Plan Optimization for a Planetary Lander

Plan execution in unknown environments poses a number of challenges: uncertainty in domain modeling, stochasticity at execution time, and the presence of exogenous events. These challenges motivate an integrated approach to planning and execution that is able to respond intelligently to variation. We examine this problem in the context of the Europa Lander mission concept, and evaluate a planning and execution framework that responds to feedback and task failure using two techniques: flexible execution and replanning with plan optimization. We develop a theoretical framework to estimate gains from these techniques, and we compare these predictions to empirical results generated in simulation. These results indicate that an integrated approach to planning and execution leveraging flexible execution, replanning, and utility maximization shows significant promise for future tightly-constrained space missions that must address significant uncertainty.

Chien, Steve↗
DOE OSTI2021

Are suspensory ligaments important for middle ear reconstruction?

As the resolution of 3D printing techniques improves, the possibility of individualized, 3-ossicle constructions adds a new dimension to middle ear prostheses. In order to optimize these designs, it is essential to understand how the ossicles and ligaments work together to transmit sound, and thus how ligaments should be replicated in a middle ear reconstruction. The middle ear ligaments are thought to play a significant role in maintaining the position of the ossicles and constraining axis of rotation. Paradoxically, investigations of the role of ligaments to date have shown very little impact on middle ear sound transmission. We explored the role of the two attachments in the gerbil middle ear analogous to human ligaments, the posterior incudal ligament and the anterior mallear process, severing both attachments and measuring change in hearing sensitivity. The impact of severing the attachments on the position of the ossicular chain was visualized using synchrotron microtomography imaging of the middle ear. In contrast to previous studies, a threshold change on the order of 20 dB across a wide range of frequencies was found when both ligaments were severed. Concomitantly, a shift in position of the ossicles was observed from the x-ray imaging and 3D renderings of the ossicular chain. These findings contrast with previous studies, demonstrating that these ligaments play a significant role in the transmission of sound through the middle ear. It appears that both mallear and incudal ligaments must be severed in order to impair sound transmission. The results of this study have significance for middle ear reconstructive surgery and the design of 3D-printed three-ossicle biocompatible prostheses.

3D printing↗
NASA NTRS1977

Linear tracking systems with applications to aircraft control system design

A class of optimal linear time invariant tracking systems, both in continuous time and discrete time, of which the number of inputs (which are restricted to be step functions) is equal to the number of system outputs, is studied. Along with derivation of equations and design procedures, two discretization schemes are presented, constraining either the control or its time derivative, to be a constant over each sampling period. Descriptions are given for the linearized model of the F-8C aircraft longitudinal dynamics, and the C* handling qualities criterion, which then serve as an illustration of the applications of these linear tracking designs. A suboptimal reduced state design is also presented. Numerical results are given for both the continuous time and discrete time designs.

Lee, W. H.↗
NASA NTRS1991

Theoretical characterization of the potential energy surface for the reversible reaction H + O2 yields HO2(asterisk) yields OH + O. III - Computed points to define a global potential energy surface

Computed energies and geometries are reported which, combined with previously published calculations, permit a global representation of the potential energy surface for the reaction H + O2 yields HO2(asterisk) yields OH + O. These new calculations characterize the potential energy surface (PES) for all H atom angles of approach to O2 and for the region of the inner repulsive wall. The region of the T-shaped H-O2 exchange saddle point is connected with the constrained energy minimum (CEM) path, and a new collinear H-O2 exchange saddle point is characterized which lies only 9 kcal/mol above the H + O2 asymptote. A vibrational analysis which utilizes local cubic and quartic polynomial representations of the PES along the CEM path has been carried out. Optimal geometries, energies, and harmonic frequencies are reported along with anharmonic analyses for the O2 and OH asymptotes and for the HO2 minimum region of the PES.

Walch, Stephen P.↗
NASA NTRS2007

Estimating the Inertia Matrix of a Spacecraft

A paper presents a method of utilizing some flight data, aboard a spacecraft that includes reaction wheels for attitude control, to estimate the inertia matrix of the spacecraft. The required data are digitized samples of (1) the spacecraft attitude in an inertial reference frame as measured, for example, by use of a star tracker and (2) speeds of rotation of the reaction wheels, the moments of inertia of which are deemed to be known. Starting from the classical equations for conservation of angular momentum of a rigid body, the inertia-matrix-estimation problem is formulated as a constrained least-squares minimization problem with explicit bounds on the inertia matrix incorporated as linear matrix inequalities. The explicit bounds reflect physical bounds on the inertia matrix and reduce the volume of data that must be processed to obtain a solution. The resulting minimization problem is a semidefinite optimization problem that can be solved efficiently, with guaranteed convergence to the global optimum, by use of readily available algorithms. In a test case involving a model attitude platform rotating on an air bearing, it is shown that, relative to a prior method, the present method produces better estimates from few data.

Acikmese, Behcet↗
NASA NTRS1975

Optimal orbit transfer for out-of-ecliptic mission using a solar electric spacecraft

This paper analyzes the performance capability of a Titan IIIE/Centaur/SEP (15 kw) for out-of-ecliptic missions. Two classes of orbital transfer modes are considered; circular-to-circular and circular-to-elliptic. The results presented compare the achievable heliographic inclination using these different transfer modes and also show the degradations in performance due to constrained launch azimuth and thruster life. It is observed that the latter transfer mode could deliver a spacecraft to a much higher inclination than the former, but the final orbit is very eccentric and could come very close to the sun. In view of the existence of many trajectories by which this mission can be performed, it is suggested that a trajectory other than IAU circular out-of-ecliptic type be used to explore more scientifically rewarding missions.

Yen, C.-W. L.↗
NASA NTRS2000

Reliability Constrained Priority Load Shedding for Aerospace Power System Automation

The need for improving load shedding on board the space station is one of the goals of aerospace power system automation. To accelerate the optimum load-shedding functions, several constraints must be involved. These constraints include congestion margin determined by weighted probability contingency, component/system reliability index, generation rescheduling. The impact of different faults and indices for computing reliability were defined before optimization. The optimum load schedule is done based on priority, value and location of loads. An optimization strategy capable of handling discrete decision making, such as Everett optimization, is proposed. We extended Everett method to handle expected congestion margin and reliability index as constraints. To make it effective for real time load dispatch process, a rule-based scheme is presented in the optimization method. It assists in selecting which feeder load to be shed, the location of the load, the value, priority of the load and cost benefit analysis of the load profile is included in the scheme. The scheme is tested using a benchmark NASA system consisting of generators, loads and network.

Momoh, James A.↗
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times (tE > 120 days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~ 85% by selecting events with both tE > 120 days and a microlensing parallax of πE < 0.08. Unfortunately, such a selection criterion cannot be applied during the event as πE requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of O(100 - 1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Application of multivariable search techniques to structural design optimization

Multivariable optimization techniques are applied to a particular class of minimum weight structural design problems: the design of an axially loaded, pressurized, stiffened cylinder. Minimum weight designs are obtained by a variety of search algorithms: first- and second-order, elemental perturbation, and randomized techniques. An exterior penalty function approach to constrained minimization is employed. Some comparisons are made with solutions obtained by an interior penalty function procedure. In general, it would appear that an interior penalty function approach may not be as well suited to the class of design problems considered as the exterior penalty function approach. It is also shown that a combination of search algorithms will tend to arrive at an extremal design in a more reliable manner than a single algorithm. The effect of incorporating realistic geometrical constraints on stiffener cross-sections is investigated. A limited comparison is made between minimum weight cylinders designed on the basis of a linear stability analysis and cylinders designed on the basis of empirical buckling data. Finally, a technique for locating more than one extremal is demonstrated.

Jones, R. T.↗
NASA NTRS2019

Airfoil Selection for Mars Rotor Applications

This short white paper is aimed at providing an overview of design considerations for airfoil choices for rotor applications in the Martian atmosphere, at very low chord-based Reynolds number flows, around Rec = O (103 -104). The low Reynolds number typical of rotorcraft operation in the Martian atmosphere reduces the rotor lifting force and efficiency, which is only partially compensated for by a lower gravity on Mars. Additionally, the low temperature and largely CO2 based atmosphere of Mars compound the overall aerodynamic problem by resulting in a lower speed of sound, further constraining rotor operation in the Martian atmosphere by limiting the maximum rotor tip speed possible so as not to exceed an acceptable tip Mach number. In light of the expected reduced rotor efficiency, evaluating airfoils for compressible, low-Reynolds number Mars rotor applications is key. Prior research on airfoil optimization and performance evaluation at low Reynolds numbers, especially in the compressible regime, is scarce and deserves further understanding. Specifically, the proposed goal stemming from this overview is to develop airfoils tailored to the unique demands of the second generation of Mars rotorcraft, i.e. the Mars Science Helicopter (MSH). This research focuses on the airfoil performance at low Reynolds numbers and hopes to add to the work performed by, amongst others, Kroo et al., Kunz and Kroo, Oyama and Fujii, Anyoji et al.

Koning, Witold J. F.↗
NASA NTRSDate not supplied

Generating a 4D Global CH(4) Product by Assimilating TROPOMI column CH(4) in NASA’s GEOS GCM

Examination of temporal and spatial CH4 variability is crucial for better understanding the human and natural processes driving climate change and ultimately designing mitigation strategies. Here we present an analysis framework that uses NASA’s GEOS General Circulation Model (GCM) to construct a high-resolution, time varying picture of atmospheric CH4 consistent with measurements from a variety of platforms, both in situ and remotely sensed. The resulting time varying atmospheric CH4 product can (i) support interpretation of high-resolution point source detection approaches, (ii) provide reanalysis fields for CH4 and other greenhouse gases, and (iii) supply boundary conditions for regional models. Our approach starts with a set of CH4 emissions from various inventories that have been adjusted to match the global annual growth rate over recent decades. These emissions are transported by the GEOS GCM, which in turn is constrained by meteorology from NASA’s Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) product. The simulated atmospheric field is compared with CH4 measurements, such as those from the TROPOspheric Monitoring Instrument (TROPOMI), and adjustments calculated following a Bayesian protocol. The accuracy of the resultant optimal atmospheric CH4 field can be demonstrated by its improved agreement (compared to a direct simulation of the CH4 inventories) with a host of independent CH4 measurements, such as those from the Total Carbon Column Observation Network (TCCON) and in situ observations from surface and airborne platforms.

Nikolay V. Balashov↗
DOE OSTI2025

Modeling distributed energy resource aggregations in security constrained unit commitment and economic dispatch

The Federal Energy Regulatory Commission (FERC) recently issued Order 2222, which requires all wholesale electricity markets in the US to allow distributed energy resources (DERs) to participate in the market as aggregated resources. These DER aggregations may be composed of many individual resources that are offered and dispatched by the market as a single entity. We present here a model of a distributed energy resource aggregator (DERA) that is scheduled by a market operator’s security constrained unit commitment (SCUC) and security constrained economic dispatch (SCED). The DERA model includes constraints for battery energy storage systems (BESSs), demand response resources (DRRs), and a simple distributed energy resource (DER). This paper describes a model for each resource type and presents two methods for the DERA to generate market offer curves: a profit-maximizing optimization to compute cost curves and a direct cost algorithm to determine dispatch costs for each resource and combine into cost curves. Once all participating DERAs are scheduled in SCUC/SCED, the model is then modified to dispatch individual DERs to maximize profit or minimize schedule deviation of the DERAs. A simulation of a representative day illustrates the DERA offers, the scheduled generation, and the DERA dispatch. Findings show the potential for unavoidable schedule deviations due to internal DER constraints and due to economic incentives to deviate from the SCUC/SCED schedules. This highlights the importance of DERA offer construction on market efficiency and system reliability. Novel aspects of our approach include: (1) We consider the asymmetry of price incentives impacting DERAs from the wholesale market compared to those impacting consumers from the retail market, as imposed by current regulations and laws. (2) We model aggregate consumer response through statistically parameterizable utility functions rather than a potentially impractical approach of modeling each individual consumer. (3) We show how to use the DERA operational dispatch model to create offers into the wholesale electricity market. (4) We show how DERAs may fail to meet their scheduled dispatch because the market offer format may not permit them to fully express their operational features such as intertemporal costs and constraints to the market.

aggregations↗
DOE OSTI2021

DC-Ripple-Energy Adaptive-Minimization (DREAM) Modulation Scheme for a High Power Density Inverter

The DC bus capacitor is one of the major power-density and reliability hurdles of electric drive systems. It is hard to shrink because it is constrained by the DC bus RMS ripple current, which is only load dependent. A dual-inverter based segmented drive can reduce the ripple current by ~50% compared to a non-segmented case. This paper analyzes the origin of this ripple current and points out the path for minimization. An optimal DC-ripple-energy adaptive-minimization (DREAM) modulation method is proposed to further reduce the ripple current. It is observed in experimental results that the proposed method can achieve additional 38% reduction over the traditional segmented drive system.

Xue, Lincoln↗
NASA NTRS2017

Linear Subpixel Learning Algorithm for Land Cover Classification from WELD using High Performance Computing

In this work, we use a Fully Constrained Least Squares Subpixel Learning Algorithm to unmix global WELD (Web Enabled Landsat Data) to obtain fractions or abundances of substrate (S), vegetation (V) and dark objects (D) classes. Because of the sheer nature of data and compute needs, we leveraged the NASA Earth Exchange (NEX) high performance computing architecture to optimize and scale our algorithm for large-scale processing. Subsequently, the S-V-D abundance maps were characterized into 4 classes namely, forest, farmland, water and urban areas (with NPP-VIIRS-national polar orbiting partnership visible infrared imaging radiometer suite nighttime lights data) over California, USA using Random Forest classifier. Validation of these land cover maps with NLCD (National Land Cover Database) 2011 products and NAFD (North American Forest Dynamics) static forest cover maps showed that an overall classification accuracy of over 91 percent was achieved, which is a 6 percent improvement in unmixing based classification relative to per-pixel-based classification. As such, abundance maps continue to offer an useful alternative to high-spatial resolution data derived classification maps for forest inventory analysis, multi-class mapping for eco-climatic models and applications, fast multi-temporal trend analysis and for societal and policy-relevant applications needed at the watershed scale.

Subpixel↗
DOE OSTI2026

Measurement of the muon neutrino charged-current mesonless cross section in the NOvA near detector

NOvA is a long-baseline accelerator neutrino experiment at Fermilab. Its physics goals include precision neutrino oscillation measurements, neutrino interaction cross-section measurements and beyond Standard Model explorations. We present a measurement of muon neutrino charged-current cross section with zero mesons in the final state at the NOvA near detector. This measurement is performed as a function of the kinematics of the final state muon. Our chosen interaction channel is especially sensitive to quasielastic and meson exchange current interactions and it provides handles for constraining the cross section systematic uncertainties in oscillation analyses in present and future experiments. For particle identification, we use a convolutional neural network (CNN) trained on individual particles simulated in the NOvA Near detector. This allows us to select the desired signal while reducing the potential bias from neutrino interaction modeling. We study strategies for constraining the remaining charged-pion background via Michel electron information in a template fitting approach. The main experimental result is a two-dimensional differential cross section as a function of final-state muon kinetic energy and polar angle. The parameters of this measurement, including binning and unfolding, were optimized to reduce the expected systematic uncertainty in the total cross section. The final result shows good agreement with the main GENIE-based simulation framework that was independently fine-tuned in NOvA. We finally propose improvements and subsequent steps that build on this analysis and further dissect the final states of neutrino interactions. This work has been supported by US DOE grant DE-SC0015684.

Sánchez Falero, Sebastián Jesús [Iowa State U.]↗
DOE OSTI2020

Communication-Constrained Expansion Planning for Resilient Distribution Systems

Distributed generation and remotely controlled switches have emerged as important technologies to improve the resiliency of distribution grids against extreme weather-related disturbances. Therefore it becomes important to study how best to place them on the grid in order to meet a resiliency criteria, while minimizing costs and capturing their dependencies on the associated communication systems that sustain their distributed operations. This paper introduces the Optimal Resilient Design Problem for Distribution and Communication Systems (ORDPDC) to address this need. The ORDPDC is formulated as a two-stage stochastic mixed-integer program that captures the physical laws of distribution systems, the communication connectivity of the smart grid components, and a set of scenarios that specifies which components are affected by potential disasters. The paper proposes an exact branch-and-price algorithm for the ORDPDC that features a strong lower bound and a variety of acceleration schemes to address degeneracy. The ORDPDC model and branch-and-price algorithm were evaluated on a variety of test cases with varying disaster intensities and network topologies. The results demonstrate the significant impact of the network topologies on the expansion plans and costs, as well as the computational benefits of the proposed approach.

97 MATHEMATICS AND COMPUTING↗