Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 541 records · Page 30

Application of mathematical optimization procedures to a structural model of a large finite-element wing

An optimization system is described that consists of the finite-element analysis program known as Engineering Analysis Language (EAL), the CONMIN optimization program, and a piecewise linear analysis. The optimization analysis system (OPTEAL) is used to minimize the mass of a transport wing subject to stress constraints for two static load conditions. The optimization procedure and the optimization analysis system (OPTEAL) are described and the optiomized results with fully stressed design results compared.

Walsh, J. L.↗

Optically resonant magneto-electric cubic nanoantennas for ultra-directional light scattering

Cubic dielectric nanoparticles are promising candidates for futuristic low-loss, ultra-compact, nanophotonic applications owing to their larger optical coefficients, greater packing density, and relative ease of fabrication as compared to spherical nanoparticles; besides possessing negligible heating at nanoscale in contrast to their metallic counterparts. Here, we present the first theoretical demonstration of azimuthally symmetric, ultra-directional Kerker's-type scattering of simple dielectric nanocubes in visible and near-infrared regions via simultaneous excitation and interference of optically induced electric- and magnetic-resonances up to quadrupolar modes. Unidirectional forward-scattering by individual nanocubes is observed at the first generalized-Kerker's condition for backward-scattering suppression, having equal electric- and magnetic-dipolar responses. Both directionality and magnitude of these unidirectional-scattering patterns get enhanced where matching electric- and magnetic-quadrupolar responses spectrally overlap. While preserving azimuthal-symmetry and backscattering suppression, a nanocube homodimer provides further directionality improvement for increasing interparticle gap, but with reduced main-lobe magnitude due to emergence of side-scattering lobes from diffraction-grating effect. We thoroughly investigate the influence of interparticle gap on scattering patterns and propose optimal range of gap for minimizing side-scattering lobes. Besides suppressing undesired side-lobes, significant enhancement in scattering magnitude and directionality is attained with increasing number of nanocubes forming a linear chain. Optimal directionality, i.e., the narrowest main-scattering lobe, is found at the wavelength of interfering quadrupolar resonances; whereas the largest main-lobe magnitude is observed at the wavelength satisfying the first Kerker's condition. These unique optical properties of dielectric nanocubes thus can revolutionize their applications at visible and near-infrared regions in the fields of nanoantennas, nanolasers, photovoltaics, and even in biomedicine.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mission Operations Planning with Preferences: An Empirical Study

This paper presents an empirical study of some nonexhaustive approaches to optimizing preferences within the context of constraint-based, mixed-initiative planning for mission operations. This work is motivated by the experience of deploying and operating the MAPGEN (Mixed-initiative Activity Plan GENerator) system for the Mars Exploration Rover Mission. Responsiveness to the user is one of the important requirements for MAPGEN, hence, the additional computation time needed to optimize preferences must be kept within reasonabble bounds. This was the primary motivation for studying non-exhaustive optimization approaches. The specific goals of rhe empirical study are to assess the impact on solution quality of two greedy heuristics used in MAPGEN and to assess the improvement gained by applying a linear programming optimization technique to the final solution.

Bresina, John L.↗

EnergyPlus Interior Radiant Heat Exchange Runtime Performance Improvements

EnergyPlus is the flagship whole-building energy simulation program developed by the US Department of Energy. This paper describes the various approaches to improve the performance of the interior radiant heat exchange algorithm. Vectorization, optimized Basic Linear Algebra Subprograms (BLAS) library, multithreading, and Graphical Processing Unit (GPU) computation are all investigated. The approach with the best performance, while maintaining drop-in compatibility with existing EnergyPlus code, is the optimized BLAS library. However, GPU computing has the potential to reduce simulation time by orders of magnitude with major refactoring to EnergyPlus.

Adams, Mark↗

Parallelized POD-based suboptimal economic model predictive control of a state-constrained Boussinesq approximation

Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. Finally, to solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.

97 MATHEMATICS AND COMPUTING↗

Multiphase Optimal and Non-Singular Power Flow by Successive Linear Approximations: Preprint

We propose an optimal power flow method for multiphase power systems. The method works for both radial and meshed networks, and is able to accommodate different types of load/source connections such as wye, delta, etc. Moreover, the method generates a sequence of non-singular points in the electrical state space and a corresponding sequence of power injections, which can be used in real-time control settings. Under certain conditions, the sequences attain the local minimum. The method has been tested using several typical networks, including the IEEE 37-bus and 123-bus test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Robust Trajectory Optimization for NRHO Rendezvous Using SPICE Kernel Relative Motion

In this paper, robust optimization is performed on trajectory correction maneuvers during the lunar lander return phase of an Artemis mission, treating the trajectory from one hour after low lunar orbit departure to arrival in the vicinity of the lunar Gateway as a relative motion problem. To enable rapid stochastic optimization techniques requiring many candidate trajectories, SPICE kernel relative motion as implemented by the Quadratic Interpolated State Transition (QIST) system is used as the underlying dynamics propagation. The optimization is performed with a genetic optimizer using linear covariance (LinCov) software in a simplified operational context, taking into account the availability of navigation sensors with varying measurement models, ranges, and accuracies. No numerical integration is used, since the relative motion around Gateway is fully characterized with the a priori computation of the QIST coefficients. Maneuver placements are computed to optimize the minimum 3σ delta-v of the trajectory, the position dispersion at a target point, and a convex combination of these two metrics. An order of magnitude runtime improvement is provided over legacy methods with less than 10% error introduced. All QIST results are shown to be in-family with legacy methods. The tradespace for optimal delta-v design is found to range from 77.0 to 93.9 m/s, while the range of optimal dispersion is between 1.4 and 11.7 km.

Relative Motion↗

Design and simulation of closed-loop ground alignment of inertial platforms with sway motion

This paper is concerned with the application of optimal estimation and control concepts to the ground alignment of inertial platforms subjected to random sway motion. A seventh-order Kalman filter has been designed for estimating the platform misalignments and sway motion. Using the linear quadratic optimal regulator techniques, feedback controller gains have been designed to drive the initial misalignments to a small value. A reconfiguration of the filter controller scheme has been proposed and verified by Monte Carlo simulation for improving the speed of alignment without changing the alignment accuracy.

Vathsal, S.↗

Linear Back-Drive Differentials

Linear back-drive differentials have been proposed as alternatives to conventional gear differentials for applications in which there is only limited rotational motion (e.g., oscillation). The finite nature of the rotation makes it possible to optimize a linear back-drive differential in ways that would not be possible for gear differentials or other differentials that are required to be capable of unlimited rotation. As a result, relative to gear differentials, linear back-drive differentials could be more compact and less massive, could contain fewer complex parts, and could be less sensitive to variations in the viscosities of lubricants. Linear back-drive differentials would operate according to established principles of power ball screws and linear-motion drives, but would utilize these principles in an innovative way. One major characteristic of such mechanisms that would be exploited in linear back-drive differentials is the possibility of designing them to drive or back-drive with similar efficiency and energy input: in other words, such a mechanism can be designed so that a rotating screw can drive a nut linearly or the linear motion of the nut can cause the screw to rotate. A linear back-drive differential (see figure) would include two collinear shafts connected to two parts that are intended to engage in limited opposing rotations. The linear back-drive differential would also include a nut that would be free to translate along its axis but not to rotate. The inner surface of the nut would be right-hand threaded at one end and left-hand threaded at the opposite end to engage corresponding right- and left-handed threads on the shafts. A rotation and torque introduced into the system via one shaft would drive the nut in linear motion. The nut, in turn, would back-drive the other shaft, creating a reaction torque. Balls would reduce friction, making it possible for the shaft/nut coupling on each side to operate with 90 percent efficiency.

Waydo, Peter↗

Numerical optimization of three-dimensional coils for NSTX-U

A tool for the calculation of optimal three-dimensional (3D) perturbative magnetic fields in tokamaks has been developed. The IPECOPT code builds upon the stellarator optimization code STELLOPT to allow for optimization of linear ideal magnetohydrodynamic perturbed equilibrium (IPEC). This tool has been applied to NSTX-U equilibria, addressing which fields are the most effective at driving NTV torques. The NTV torque calculation is performed by the PENT code. Optimization of the normal field spectrum shows that fields with n = 1 character can drive a large core torque. It is also shown that fields with n = 3 features are capable of driving edge torque and some core torque. Coil current optimization (using the planned in-vessel and existing RWM coils) on NSTX-U suggest the planned coils set is adequate for core and edge torque control. In conclusion, comparison between error field correction experiments on DIII-D and the optimizer show good agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model reduction of dynamical systems on nonlinear manifolds using deep convolutional autoencoders

Nearly all model-reduction techniques project the governing equations onto a linear subspace of the original state space. Such subspaces are typically computed using methods such as balanced truncation, rational interpolation, the reduced-basis method, and (balanced) proper orthogonal decomposition (POD). Unfortunately, restricting the state to evolve in a linear subspace imposes a fundamental limitation to the accuracy of the resulting reduced-order model (ROM). In particular, linear-subspace ROMs can be expected to produce low-dimensional models with high accuracy only if the problem admits a fast decaying Kolmogorov n -width (e.g., diffusion-dominated problems). Unfortunately, many problems of interest exhibit a slowly decaying Kolmogorov n-width (e.g., advection-dominated problems). To address this, we propose a novel framework for projecting dynamical systems onto nonlinear manifolds using minimum-residual formulations at the time-continuous and time-discrete levels; the former leads to manifold Galerkin projection, while the latter leads to manifold least-squares Petrov–Galerkin (LSPG) projection. We perform analyses that provide insight into the relationship between these proposed approaches and classical linear-subspace reduced-order models; we also derive a posteriori discrete-time error bounds for the proposed approaches. In addition, we propose a computationally practical approach for computing the nonlinear manifold, which is based on convolutional autoencoders from deep learning. Lastly, we demonstrate the ability of the method to significantly outperform even the optimal linear-subspace ROM on benchmark advection-dominated problems, thereby demonstrating the method's ability to overcome the intrinsic n -width limitations of linear subspaces.

97 MATHEMATICS AND COMPUTING↗

Linear model decision trees as surrogates in optimization of engineering applications

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. Importantly, we observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

42 ENGINEERING↗

Dynamic Optimal Dispatch of Energy Systems with Intermittent Renewables and Damage Model

With the increasing penetration of intermittent renewable energy sources into the grid, there is a growing need for process systems-based strategies that integrate dispatchable and variable energy systems for supplying the demand while maintaining grid reliability. The proposed framework corresponds to a dynamic mixed-integer linear programming optimization approach that integrates coal-fired and natural gas-fired power plants, NaS batteries for energy storage, and solar/wind energy to supply the demand. This optimization approach considers an economic goal and constraints to provide power balance while maintaining the overall damage of the natural gas combined cycle (NGCC) power plant drum under a maximum stress as well as avoiding the overheating of the NGCC superheater and reheater. Renewable curtailment levels are also retained at minimum levels. Case studies are analyzed considering different loads and renewable penetration levels. The results show that the demand was met for all cases. Grid flexibility was mostly provided by the NGCC, while the batteries were used sparingly. In addition, considering a CO 2 equivalent analysis, the environmental performance was intrinsically connected to grid flexibility and the level of renewable penetration. Stress analysis results reinforced the necessity for an equipment health-related constraint.

97 MATHEMATICS AND COMPUTING↗

Efficient BiVO 4 /CoFeO x H y photoanodes using controlled annealing and conformal linear-sweep electrocatalyst photodeposition

Although monoclinic bismuth vanadate (BiVO 4 ) is a promising photoanode for solar water splitting, its practical use is hindered by imperfect photocurrent generation/collection, low photovoltage compared to the bandgap, and corrosion side reactions that limit durability. Here, we introduce a controlled-annealing sol–gel process for BiVO 4 thin-film photoanodes along with an optimized linear-sweep-voltammetry photodeposition of CoFeO x H y cocatalysts. The resulting BiVO 4 films annealed at 550 °C exhibited a photocurrent density of 4.1 mA/cm 2 at 1.23 V RHE under 1 sun AM 1.5G solar simulation and a low onset potential of 0.26 V RHE due to high majority carrier conductivity, a crystalline bulk with reduced defects as evidenced by x-ray photoelectron spectroscopy and photoluminescence lifetime analysis, and thus enhanced photocarrier collection. However, significant degradation in performance was found due to interfacial photocorrosion. To protect the surface and speed the oxygen-evolution reaction CoFeO x H y cocatalyst layers were deposited. By varying the number of consecutive sweeps and adjusting the applied bias range, an ultra-thin (~15 nm) CoFeO x H y cocatalyst layer was uniformly grown deposited over 30 cycles on the BiVO 4 surface. The resulting BiVO 4 /CoFeO x H y yielded 4.03 mA/cm 2 at 1.23 V RHE and onset potential of 0.24 V RHE , with stable operation (~15 % loss in photocurrent at 1.23 V RHE relative to ~60 % loss in the uncatalyzed control sample). These conformal CoFeO x H y catalytic layers function simultaneously to selectively collect photoexcited holes from the BiVO 4 , catalyze the water-oxidation reaction, and protect the BiVO 4 from photodegradation.

42 ENGINEERING↗

Automated design and optimization of flexible booster autopilots via linear programming, volume 1

A nonlinear programming technique was developed for the automated design and optimization of autopilots for large flexible launch vehicles. This technique, which resulted in the COEBRA program, uses the iterative application of linear programming. The method deals directly with the three main requirements of booster autopilot design: to provide (1) good response to guidance commands; (2) response to external disturbances (e.g. wind) to minimize structural bending moment loads and trajectory dispersions; and (3) stability with specified tolerances on the vehicle and flight control system parameters. The method is applicable to very high order systems (30th and greater per flight condition). Examples are provided that demonstrate the successful application of the employed algorithm to the design of autopilots for both single and multiple flight conditions.

Hauser, F. D.↗

Minimum weight design of rotorcraft blades with multiple frequency and stress constraints

Minimum weight designs of helicopter rotor blades with constraints on multiple coupled flap-lag natural frequencies are studied. Constraints are imposed on the minimum value of the blade autorotational inertia to ensure sufficient rotary inertia to autorotate in case of engine failure and on stresses to guard against structural failure due to blade centrifugal forces. Design variables include blade taper ratio, dimensions of the box beam located inside the airfoil and magnitudes of nonstructural weights. The program CAMRAD is used for the blade modal analysis; the program CONMIN is used for the optimization. A linear approximation involving Taylor series expansion is used to reduce the analysis effort. The procedure contains a sensitivity analysis consisting of analytical derivatives for objective function and constraints on autorotational inertia and stresses. Central finite difference derivatives are used for frequency constraints. Optimal designs are obtained for both rectangular and tapered blades. Using this method, it is possible to design a rotor blade with reduced weight, when compared to a baseline blade, while satisfying all the imposed design requirements.

Chattopadhyay, Aditi↗