Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convex programming”

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 181 records · Page 10

Global optimization methods for engineering design

The problem is to find a global minimum for the Problem P. Necessary and sufficient conditions are available for local optimality. However, global solution can be assured only under the assumption of convexity of the problem. If the constraint set S is compact and the cost function is continuous on it, existence of a global minimum is guaranteed. However, in view of the fact that no global optimality conditions are available, a global solution can be found only by an exhaustive search to satisfy Inequality. The exhaustive search can be organized in such a way that the entire design space need not be searched for the solution. This way the computational burden is reduced somewhat. It is concluded that zooming algorithm for global optimizations appears to be a good alternative to stochastic methods. More testing is needed; a general, robust, and efficient local minimizer is required. IDESIGN was used in all numerical calculations which is based on a sequential quadratic programming algorithm, and since feasible set keeps on shrinking, a good algorithm to find an initial feasible point is required. Such algorithms need to be developed and evaluated.

Arora, Jasbir S.↗

NASA/USRA advanced space design program: The laser powered interorbital vehicle

A preliminary design is presented for a low-thrust Laser Powered Interorbital Vehicle (LPIV) intended for cargo transportation between an earth space station and a lunar base. The LPIV receives its power from two iodide laser stations, one orbiting the earth and the other located on the surface of the moon. The selected mission utilizes a spiral trajectory, characteristic of a low-thrust spacecraft, requiring 8 days for a lunar rendezvous and an additional 9 days for return. The ship's configuration consists primarily of an optical train, two hydrogen plasma engines, a 37.1 m box beam truss, a payload module, and fuel tanks. The total mass of the vehicle fully loaded is 63300 kg. A single plasma, regeneratively cooled engine design is incorporated into the two 500 N engines. These are connected to the spacecraft by turntables which allow the vehicle to thrust tangentially to the flight path. Proper collection and transmission of the laser beam to the thrust chambers is provided through the optical train. This system consists of the 23 m diameter primary mirror, a convex parabolic secondary mirror, a beam splitter and two concave parabolic tertiary mirrors. The payload bay is capable of carrying 18000 kg of cargo. The module is located opposite the primary mirror on the main truss. Fuel tanks carrying a maximum of 35000 kg of liquid hydrogen are fastened to tracks which allow the tanks to be moved perpendicular to the main truss. This capability is required to prevent the center of mass from moving out of the thrust vector line. The laser beam is located and tracked by means of an acquisition, pointing and tracking system which can be locked onto the space-based laser station. Correct orientation of the spacecraft with the laser beam is maintained by control moment gyros and reaction control rockets. Additionally an aerobrake configuration was designed to provide the option of using the atmospheric drag in place of propulsion for a return trajectory.

Source record↗

Volumetric Assessment of UPRITE Exercises From Marker-Based Motion Capture

BACKGROUND Lack of volumetric data on full-body movement of exercises presents a challenge to ensuring the fit of crew member’s full range of motion on the International Space Station (ISS). The Upright Proprioception Retention via In-flight Training and Evaluation (UPRITE) is a sensorimotor countermeasure device designed for maintaining crew members’ proprioception in a microgravity environment. A footplate—attached to a static base—rotates in two degrees of freedom (pitch and roll) up to a 20 deg angle. An initial volumetric assessment assuming an upright standing posture produced a cone-like shape with a narrow bottom and wide top. Such general volumetric assessments risk creating an overly conservative volume estimate, taking up more space than is necessary on the already limited interior space of the ISS, and neglecting necessary volume due to oversimplifying assumptions. Rather, higher-fidelity volumetric assessments offer more comprehensive insights in an environment where every area counts. The main objective of this work is to provide the spatial parameters of exercises on the UPRITE such that it is placed on the ISS according to its volumetric demands or that usage is adjusted to fit the available space. METHODS In 2023, a data collection was performed originally to inform loads and dynamics of system use and was recently leveraged for volumetric assessment. Three human subjects representing different body types (~63-76 inches in stature) performed a variety of board manipulations using UPRITE with body weight offload. The test collected the 3D positional data of a modified full-body Plug-in Gait marker set [1] via a 16-camera OptiTrack MoCap system. After processing – filling marker gaps and trimming data – in OptiTrack Motive, the recorded marker location data, which included device markers, was exported to a readable trajectory file. To accurately represent the full volume defining landmarks, additional markers were digitally added to an unscaled Modified Full Body Model [2]. The model was then scaled according to its subject parameters upon which an inverse kinematics analysis was performed. A custom plugin yielded model marker location data files. Volumetric analyses were performed on the recorded trajectory and model trajectory files using a custom Python-built tool that extracted the marker location data and plotted it in a 3D space. Concerned with only the maximum volume of the motion, a 3D convex hull analysis was applied to the plot, extracting the vertices or external points of the eventual 3D CAD output, dubbed aptly as a “volume shell”. This overall approach was based on guidance in a NASA-STD-3001 Technical Brief [3]. RESULTS AND DISCUSSION Batch volumetric assessment on the exercises for each subject was performed, producing high-fidelity volume shells in minimal time. Preliminary results highlighted the value in higher-fidelity volumes based on collected data when possible. For example, revolving a single posture in the cone assessment would not have sufficiently captured a single leg stance; rather, it would need to involve swinging the leg both forward and back. Additional observations and the maximal dimensions of the volumes, including those based on scaled data for ISS anthropometric requirements, will be presented at the Human Research Program Investigator’s Workshop. CONCLUSIONS While this work’s primary objective was for the UPRITE-to-ISS integration, the tool built to conduct this analysis has wide applications for future exercise systems as an informational tool for optimal device placement. The tool and its findings also have implications for exercise device design and spacecraft interior considerations on Gateway, the Lunar Pressurized Rover, and beyond. REFERENCES [1] Bell, C. A., et al. (2023) Recent Improvements and Verification of a Full Body Model in OpenSim. NASA Human Research Program Investigator’s Workshop. https://ntrs.nasa.gov/citations/20230001080 [2] Lostroscio, K., et al (2023) The Digital Astronaut Simulation. AHFE International Conference on Human Factors in Design, Engineering, and Computing for All. [3] Exercise Overview. (2023) NASA-STD-3001 Technical Brief. https://www.nasa.gov/wp-content/uploads/2023/12/ochmo-tb-031-exercise-overview.pdf?emrc=9d454c?emrc=9d454c

L D Quinto↗

A Risk-Averse Approach for Distribution Grid Expansion Planning

Recent episodes of natural disasters have challenged the resilience of power grids. Adequate distribution grid planning that properly captures the risk aversion of the utility system planner is a key factor to increase the flexibility of distribution networks to circumvent these events. In this paper, we propose a methodology to determine the optimal portfolio of investments in lines and storage devices in order to minimize a convex combination between expected value and CVaR of operational costs, including energy not served, while taking into account the multistage nature of the energy storage management within this context. While the expected value of energy not served has been traditionally employed to tackle routine failures, we also minimize the CVaR of energy not served to address high-impact, low-probability (HILP) events. We illustrate the performance of the proposed methodology with a 54-Bus system test case.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The pattern of secondary root formation in curving roots of Arabidopsis thaliana (L.) Heynh

A gravitational stimulus was used to induce the curvature of the main root of Arabidopsis thaliana. The number of secondary roots increased on the convex side and decreased on the concave side of any curved main root axes in comparison with straight roots used as the control. The same phenomenon was observed with the curved main roots of plants grown on a clinostat and of mutant plants exhibiting random root orientation. The data suggest that the pattern of lateral root formation is associated with curvature but is independent of the environmental stimuli used to induce curvature.

Non-NASA Center↗

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗

The effect of direct and reflected radiation on the temperatures of space structures

As an orbiting structure passes from the bright side to the dark side of the earth, substantial differences in surface heating may take place due to direct sunlight and albedo. Most structures are built of convex curved elements to spread out the reflected energy. A change in the shortwave absorptivity and emissivity maintains the temperature of electronic components within a narrow range. Temperature extremes may also be reduced by the use of thermal capacity of the system or by shielding the surface from direct sunlight. For an instrument deployed from a Shuttle, a temperature must be established for it in the cargo bay before deployment. The effect of obstructed views by surfaces such as the Shuttle arm complicates solar flux computations because of the object's mobility. For typical orbits, the earth's albedo and the direct long wavelength radiation may be as much as 30 percent of the direct solar flux. Irradiation due to direct and reflected solar, long wave radiation and reflected radiation from earth, and from Shuttle surfaces is largely dependent on the specular reflection from the radiators as the Shuttle rolls and its surfaces block the sun. Through a judicious choice of surface properties, thermal capacitance and orbit attitude programming, the temperature response of a system can be controlled.

Emery, A. F.↗

Computer-based video digitizer analysis of surface extension in maize roots: kinetics of growth rate changes during gravitropism

We used a video digitizer system to measure surface extension and curvature in gravistimulated primary roots of maize (Zea mays L.). Downward curvature began about 25 +/- 7 min after gravistimulation and resulted from a combination of enhanced growth along the upper surface and reduced growth along the lower surface relative to growth in vertically oriented controls. The roots curved at a rate of 1.4 +/- 0.5 degrees min-1 but the pattern of curvature varied somewhat. In about 35% of the samples the roots curved steadily downward and the rate of curvature slowed as the root neared 90 degrees. A final angle of about 90 degrees was reached 110 +/- 35 min after the start of gravistimulation. In about 65% of the samples there was a period of backward curvature (partial reversal of curvature) during the response. In some cases (about 15% of those showing a period of reverse bending) this period of backward curvature occurred before the root reached 90 degrees. Following transient backward curvature, downward curvature resumed and the root approached a final angle of about 90 degrees. In about 65% of the roots showing a period of reverse curvature, the roots curved steadily past the vertical, reaching maximum curvature about 205 +/- 65 min after gravistimulation. The direction of curvature then reversed back toward the vertical. After one or two oscillations about the vertical the roots obtained a vertical orientation and the distribution of growth within the root tip became the same as that prior to gravistimulation. The period of transient backward curvature coincided with and was evidently caused by enhancement of growth along the concave and inhibition of growth along the convex side of the curve, a pattern opposite to that prevailing in the earlier stages of downward curvature. There were periods during the gravitropic response when the normally unimodal growth-rate distribution within the elongation zone became bimodal with two peaks of rapid elongation separated by a region of reduced elongation rate. This occurred at different times on the convex and concave sides of the graviresponding root. During the period of steady downward curvature the elongation zone along the convex side extended farther toward the tip than in the vertical control. During the period of reduced rate of curvature, the zone of elongation extended farther toward the tip along the concave side of the root. The data show that the gravitropic response pattern varies with time and involves changes in localized elongation rates as well as changes in the length and position of the elongation zone. Models of root gravitropic curvature based on simple unimodal inhibition of growth along the lower side cannot account for these complex growth patterns.

NASA Discipline Plant Biology↗

Touch and gravitropic set-point angle interact to modulate gravitropic growth in roots

Plant roots must sense and respond to a variety of environmental stimuli as they grow through the soil. Touch and gravity represent two of the mechanical signals that roots must integrate to elicit the appropriate root growth patterns and root system architecture. Obstacles such as rocks will impede the general downwardly directed gravitropic growth of the root system and so these soil features must be sensed and this information processed for an appropriate alteration in gravitropic growth to allow the root to avoid the obstruction. We show that primary and lateral roots of Arabidopsis do appear to sense and respond to mechanical barriers placed in their path of growth in a qualitatively similar fashion. Both types of roots exhibited a differential growth response upon contacting the obstacle that directed the main axis of elongation parallel to the barrier. This growth habit was maintained until the obstacle was circumvented, at which point normal gravitropic growth was resumed. Thus, the gravitational set-point angle of the primary and lateral roots prior to encountering the barrier were 95 degrees and 136 degrees respectively and after growing off the end of the obstacle identical set-point angles were reinstated. However, whilst tracking across the barrier, quantitative differences in response were observed between these two classes of roots. The root tip of the primary root maintained an angle of 136 degrees to the horizontal as it traversed the barrier whereas the lateral roots adopted an angle of 154 degrees. Thus, this root tip angle appeared dependent on the gravitropic set-point angle of the root type with the difference in tracking angle quantitatively reflecting differences in initial set-point angle. Concave and convex barriers were also used to analyze the response of the root to tracking along a continuously varying surface. The roots maintained the a fairly fixed angle to gravity on the curved surface implying a constant resetting of this tip angle/tracking response as the curve of the surface changed. We propose that the interaction of touch and gravity sensing/response systems combine to strictly control the tropic growth of the root. Such signal integration is likely a critical part of growth control in the stimulus-rich environment of the soil. c2003 COSPAR. Published by Elsevier Ltd. All rights reserved.

NASA Discipline Plant Biology↗

Physiological asymmetry in etiolated pea epicotyls: relation to patterns of auxin distribution and phototropic behavior

Etiolated pea seedlings require transformation of Pr phytochrome to Pfr before they display optimal phototropic response to unilateral blue light. This study investigates the possible role of auxin transport in explaining these phenomena. Labeled [2-14C]IAA applied to the intact terminal buds of dark-grown and red light-treated pea seedlings was measured 210 min later on the shaded and illuminated sides of the epicotyl as a function of direction and duration of irradiation with blue light. Totally darkened epicotyls show an asymmetry in distribution of radioactivity in the upper growth zone of the epicotyl, in favor of the side under the concave part of the apical hook. Red light, which greatly potentiates curvature toward subsequent unilateral blue light, lowers this asymmetry. Blue light directed to the epicotyl of red-pretreated plants in a plane parallel to the hook and from the side bearing the convex portion of the hook induces positive phototropic curvature as well as a surplus of radioactivity on the illuminated side of the upper epicotyl and on the shaded side of the lower growth zone of the epicotyl. Light directed to the side bearing the concave part of the hook also causes an accumulation of counts in the upper part of the lighted side but produces neither curvature of the epicotyl nor accumulation of counts in the lower shaded side. Because of this built-in physiological asymmetry in the growth zone just below the apical hook, it is difficult to explain the effects of red and blue light on curvature in terms of patterns of auxin distribution alone.

NASA Program Space Biology↗

Hardware in the Loop Performance of Terrestrial Powered Descent Dual Quaternion Guidance With A Custom First-Order Solver

The Safe & Precise Landing Integrated Capabilities Evolution (SPLICE) program continues to push the development of advanced descent and landing technologies. This paper will focus on the improvements made to the SPLICE Dual Quaternion Guidance (DQG) algorithm, which is a critical software component for generating approach phase and hazard avoidance trajectories. The trajectory performance of SPLICE DQG will be evaluated in preparation for a terrestrial flight test that incorporates a hazard avoidance maneuver. Additionally, the implementation of a customized first order Proportional-Integral Projected Gradient solver will be reviewed, along with preliminary execution results on the SPLICE Descent and Landing Computer (DLC).

SPLICE↗

CometBoards Users Manual Release 1.0

Several nonlinear mathematical programming algorithms for structural design applications are available at present. These include the sequence of unconstrained minimizations technique, the method of feasible directions, and the sequential quadratic programming technique. The optimality criteria technique and the fully utilized design concept are two other structural design methods. A project was undertaken to bring all these design methods under a common computer environment so that a designer can select any one of these tools that may be suitable for his/her application. To facilitate selection of a design algorithm, to validate and check out the computer code, and to ascertain the relative merits of the design tools, modest finite element structural analysis programs based on the concept of stiffness and integrated force methods have been coupled to each design method. The code that contains both these design and analysis tools, by reading input information from analysis and design data files, can cast the design of a structure as a minimum-weight optimization problem. The code can then solve it with a user-specified optimization technique and a user-specified analysis method. This design code is called CometBoards, which is an acronym for Comparative Evaluation Test Bed of Optimization and Analysis Routines for the Design of Structures. This manual describes for the user a step-by-step procedure for setting up the input data files and executing CometBoards to solve a structural design problem. The manual includes the organization of CometBoards; instructions for preparing input data files; the procedure for submitting a problem; illustrative examples; and several demonstration problems. A set of 29 structural design problems have been solved by using all the optimization methods available in CometBoards. A summary of the optimum results obtained for these problems is appended to this users manual. CometBoards, at present, is available for Posix-based Cray and Convex computers, Iris and Sun workstations, and the VM/CMS system.

Guptill, James D.↗

Polyhedral Relaxations for Optimal Pump Scheduling of Potable Water Distribution Networks

The classic pump scheduling or optimal water flow (OWF) problem for water distribution networks (WDNs) minimizes the cost of power consumption for a given WDN over a fixed time horizon. In its exact form, the OWF is a computationally challenging mixed-integer nonlinear program (MINLP). It is complicated by nonlinear equality constraints that model network physics, discrete variables that model operational controls, and intertemporal constraints that model changes to storage devices. To address the computational challenges of the OWF, this paper develops tight polyhedral relaxations of the original MINLP, derives novel valid inequalities (or cuts) using duality theory, and implements novel optimization-based bound tightening and cut generation procedures. The efficacy of each new method is rigorously evaluated by measuring empirical improvements in OWF primal and dual bounds over 45 literature instances. The evaluation suggests that our relaxation improvements, model strengthening techniques, and a thoughtfully selected polyhedral relaxation partitioning scheme can substantially improve OWF primal and dual bounds, especially when compared with similar relaxation-based techniques that do not leverage these new methods.

bound tightening↗

An Incremental Gradient Method for Optimization Problems With Variational Inequality Constraints

We consider minimizing a sum of agent-specific nondifferentiable merely convex functions over the solution set of a variational inequality (VI) problem in that each agent is associated with a local monotone mapping. This problem finds an application in computation of the best equilibrium in nonlinear complementarity problems arising in transportation networks. We develop an iteratively regularized incremental gradient method where at each iteration, agents communicate over a directed cycle graph to update their solution iterates using their local information about the objective and the mapping. The proposed method is single-timescale in the sense that it does not involve any excessive hard-to-project computation per iteration. We derive nonasymptotic agent-wise convergence rates for the suboptimality of the global objective function and infeasibility of the VI constraints measured by a suitably defined dual gap function. Finally, the proposed method appears to be the first fully iterative scheme equipped with iteration complexity that can address distributed optimization problems with VI constraints over cycle graphs.

convergence↗

Architecture-Preserving Provable Repair of Deep Neural Networks

Deep neural networks (DNNs) are becoming increasingly important components of software, and are considered the state-of-the-art solution for a number of problems, such as image recognition. However, DNNs are far from infallible, and incorrect behavior of DNNs can have disastrous real-world consequences. This paper addresses the problem of architecture-preserving V-polytope provable repair of DNNs. A V-polytope defines a convex bounded polytope using its vertex representation. V-polytope provable repair guarantees that the repaired DNN satisfies the given specification on the infinite set of points in the given V-polytope. An architecture-preserving repair only modifies the parameters of the DNN, without modifying its architecture. The repair has the flexibility to modify multiple layers of the DNN, and runs in polynomial time. It supports DNNs with activation functions that have some linear pieces, as well as fully-connected, convolutional, pooling and residual layers. To the best our knowledge, this is the first provable repair approach that has all of these features. We implement our approach in a tool called APRNN. Using MNIST, ImageNet, and ACAS Xu DNNs, we show that it has better efficiency, scalability, and generalization compared to PRDNN and REASSURE, prior provable repair methods that are not architecture preserving.

97 MATHEMATICS AND COMPUTING↗

Grid Optimization Competition on Synthetic and Industrial Power Systems

This paper summarizes a grid optimization (GO) competition effort in the United States to find the best solution strategies for up to interconnect-scale power system networks with around 32,000 buses. The optimization problem is a mixedinteger, non-convex non-linear problem, (MINLP) and includes discrete variables such as unit commitment and line switching, control settings (transformer taps and phase shifters with impedance correction tables), and bus shunts. The case study includes six actual industry grids as well as 16 realistic synthetic grids created by three different dataset teams. The winners are selected and ranked based on scoring criteria, which consider the solution quality (such as objective functions) within time limits. Nine winner teams are selected from 26 competitor teams. The results achieved by different teams are described and the performance of different algorithms on synthetic grids and actual industry grids are compared and analyzed.

mixed-integer non-linear programming↗

Automated Resonance Fitting for Nuclear Data Evaluation

Global and national efforts to deliver high-quality nuclear data to users have a wide-ranging impact, affecting applications in national security, reactor operations, basic science, medicine, and more. Cross section evaluation is a major part of this effort, combining theory and experimentation to produce recommended values and uncertainties for reaction probabilities. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant manual effort and months of time from expert scientists. In this article, non-convex non-linear optimization methods are combined with concepts of inferential statistics to infer a resonance model from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and the potential for bias from prior assumptions, while enhancing reproducibility and documentation, thereby addressing well-known challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-Frequency Testing of Composite Fan Vanes With Erosion-Resistant Coating Conducted

The mechanical integrity of hard, erosion-resistant coatings were tested using the Structural Dynamics Laboratory at the NASA Glenn Research Center. Under the guidance of Structural Mechanics and Dynamics Branch personnel, fixturing and test procedures were developed at Glenn to simulate engine vibratory conditions on coated polymer-matrix- composite bypass vanes using a slip table in the Structural Dynamics Laboratory. Results from the high-frequency mechanical bench testing, along with concurrent erosion testing of coupons and vanes, provided sufficient confidence to engine-endurance test similarly coated vane segments. The knowledge gained from this program will be applied to the development of oxidation- and erosion-resistant coatings for polymer matrix composite blades and vanes in future advanced turbine engines. Fan bypass vanes from the AE3007 (Rolls Royce America, Indianapolis, IN) gas turbine engine were coated by Engelhard (Windsor, CT) with compliant bond coatings and hard ceramic coatings. The coatings were developed collaboratively by Glenn and Allison Advanced Development Corporation (AADC)/Rolls Royce America through research sponsored by the High-Temperature Engine Materials Technology Project (HITEMP) and the Higher Operating Temperature Propulsion Components (HOTPC) project. High-cycle fatigue was performed through high-frequency vibratory testing on a shaker table. Vane resonant frequency modes were surveyed from 50 to 3000 Hz at input loads from 1g to 55g on both uncoated production vanes and vanes with the erosion-resistant coating. Vanes were instrumented with both lightweight accelerometers and strain gauges to establish resonance, mode shape, and strain amplitudes. Two high-frequency dwell conditions were chosen to excite two strain levels: one approaching the vane's maximum allowable design strain and another near the expected maximum strain during engine operation. Six specimens were tested per dwell condition. Pretest and posttest inspections were performed optically at up to 60 magnification and using a fluorescent-dye penetrant. Accumulation of 10 million cycles at a strain amplitude of two to three times that expected in the engine (approximately 670 Hz and 20g) led to the development of multiple cracks in the coating that were only detectable using fluorescent-dye penetrant inspection. Cracks were prevalent on the trailing edge and on the convex side of the midsection. No cracking or spalling was evident using standard optical inspection at up to 60 magnification. Further inspection may reveal whether these fine cracks penetrated the coating or were strictly on the surface. The dwell condition that simulated actual engine conditions produced no obvious surface flaws even after up to 80 million cycles had been accumulated at strain amplitudes produced at approximately 1500 Hz and 45g.

Bowman, Cheryl L.↗