Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse 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 433 records · Page 24

Modified Cascading Generalized Inverse Control Allocation

The current aviation revolution towards electric propulsion aircraft (e.g., electric vertical takeoff-and-landing) brings unique control challenges. These vehicles are typically over-actuated (more effectors than desired control outcomes), may require control strategies for the three phases of flight (hover, transition and cruise), and currently have limited electric power availability. These vehicle challenges bring the need for optimal control allocation to the forefront of research. A leading control allocation algorithm, used in current flight vehicles, is the Cascading Generalized Inverse (CGI). Unfortunately, the Cascading Generalized Inverse algorithm is unable to achieve some desired outcomes, it intermittently provides non-optimal allocations, and it may fail to preserve moment direction near maximal achievable outcomes. In this research, the shortcomings of the Cascading Generalized Inverse algorithm are addressed by augmenting the algorithm with Scalar Difference Quadratic unsaturation identification and location at each iteration. Rigorous theory is shown that the Modified Cascading Generalized Inverse performs better at obtaining optimal allocations for all attainable outcomes. Numerical case studies for over-actuated vehicles demonstrate resolution to the aforementioned deficiencies.

Control Allocation↗

Multiobjective Optimization of Regional and Teleseismic Data to Constrain the Source of the 12 September 2016 Mw 5.4 Earthquake in South Korea

In this paper, we propose a new approach for source inversion with automatic criteria for the station selection process and using a multiobjective optimization scheme combining teleseismic body-wave and regional surface-wave observations. To illustrate the approach, we determine the source parameters of the 12 September 2016 M w 5.4 South Korea earthquake, as well as its foreshock (same day) and its main aftershock on 19 September 2016. Available teleseismic (short period and broadband) and regional (broadband) seismograms from the Incorporated Research Institutions for Seismology International–Federation of Digital Seismograph Networks (IRIS-FDSN) webservice are selected (see Data and Resources). First, all teleseismic data with signal-to-noise ratio (SNR) above 2 in the 0.3–1.5 Hz frequency band are selected and a cepstrum function is blindly applied to help detect reflected wave arrivals inside a P coda. This method allows the identification of records with weak P-coda and clear depth phase (pP, sP) arrivals and can be used to select an optimal set of teleseismic seismograms to be used for the inversion. Second, regional seismograms are selected according to the fit between their observed and theoretical dispersion curves extracted from local 1D velocity model. Depending on the epicentral distance, two frequency bands are proposed for regional data: 15–50 s for epicentral distances smaller than 2° and 30–80 s otherwise. Finally, we explore the space parameters from exhaustive grid searches for the two objective functions Tf and Rf for teleseismic and regional data, respectively. We show that the seismic moment parameters can be exclusively derived from Rf and can be fixed for Tf. However, the depth parameter is not sensitive to Rf and can be simultaneously inverted from both. The selection of the optimal solution satisfying both cost functions uses the optimum Pareto front that allows a decrease in epistemic uncertainty. Lastly, we discuss the advantage of this approach for improving the characterization of moderate earthquake sources in areas with limited regional instrumentation.

58 GEOSCIENCES↗

Exergy-based Sensitivity Analysis of the Generic Hypersonic Vehicle using FUN3D

In this paper, the implementation of an exergy-based objective function and its adjoint gradient into NASA’s FUN3D solver is discussed and verified. In order to verify that the exergy-based functional is properly implemented, it is used to predict the drag of the Generic Hypersonic Vehicle (GHV), which is then compared to more traditional force-based drag predictions. In addition to the functional implementation, FUN3D’s adjoint capability was extended to obtain sensitivities. Results were verified using FUN3D’s native complex step method for di↵erentiation using a generic wing configuration. The complex and adjoint gradients yielded discrete agreement demonstrating correct implementation and that the functional can be used for gradient-based multidisciplinary analysis and optimization. Next, various trade studies are conducted on the GHV to understand the design space of the vehicle. Finally, an inverse design problem is solved to verify the utilized design optimization framework which is ready to be deployed for exergy-based optimizations in future work.

Neal L. Novotny↗

Nonequilibrium steady state transport of collective-qubit system in strong coupling regime

We investigate the steady state photon transport in a nonequilibrium collective-qubit model. By adopting the noninteracting blip approximation, which is applicable in the strong photon–qubit coupling regime, we describe the essential contribution of indirect qubit–qubit interaction to the population distribution, mediated by the photonic baths. The linear relations of both the optimal flux and noise power with the qubits system size are obtained. Moreover, the inversed power-law style for the finite-size scaling of the optimal photon–qubit coupling strength is exhibited, which is proposed to be universal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Approximate solutions of acoustic 3D integral equation and their application to seismic modeling and full-waveform inversion

Over the recent decades, a number of fast approximate solutions of Lippmann–Schwinger equation, which are more accurate than classic Born and Rytov approximations, were proposed in the field of electromagnetic modeling. Those developments could be naturally extended to acoustic and elastic fields; however, until recently, they were almost unknown in seismology. This paper presents several solutions of this kind applied to acoustic modeling for both lossy and lossless media. We evaluated the numerical merits of those methods and provide an estimation of their numerical complexity. In our numerical realization we use the matrix-free implementation of the corresponding integral operator. We study the accuracy of those approximate solutions and demonstrate, that the quasi-analytical approximation is more accurate, than the Born approximation. Further, we apply the quasi-analytical approximation to the solution of the inverse problem. It is demonstrated that, this approach improves the estimation of the data gradient, comparing to the Born approximation. The developed inversion algorithm is based on the conjugate-gradient type optimization. Numerical model study demonstrates that the quasi-analytical solution significantly reduces computation time of the seismic full-waveform inversion. We also show how the quasi-analytical approximation can be extended to the case of elastic wavefield.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recording of pulsed currents by a fibre-optic Faraday effect-based sensor with limited frequency band

We have studied, experimentally and theoretically, the dependence of the amplitude and noise response parameters of a fibre-optic Faraday effect-based current sensor on the photodetector frequency band in recording rectangular current pulses, the duration of which does not exceed the time of light propagation through the sensing fibre. The presence of an optimal band at which the signal-to-noise ratio has a pronounced maximum is demonstrated. It is shown that the optimal bandwidth is close in magnitude to the inverse time of light passage through the fibre. The revealed feature is explained by the difference in dependences of the response amplitude and the white noise level of optical radiation on the frequency band of the photodetector. A technique for amplitude error correction at the optimal band is proposed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Invertible Neural Networks for Aerodynamic Design of Wind Turbine Blades

The state-of-the-practice methods for aerodynamic design of wind turbine blads use linearized blade element momentum theory (BEM) to optimize the twist and chord profiles from a pre-selected set of 2D airfoil shapes. In this work, we apply invertible neural network (INN) tools to enable the rapid inverse aerodynamic design of wind turbine blades including component airfoils. The INN is trained on data obtained through the use of robust automated mesh generation and the HAMSTR computational fluid dynamics solver with advanced turbulence and transition models validated for turbine applications. Our design technique is a significant improvement over the state-of-the-practice linearized blade element momentum (BEM) techniques in capturing 3D nonlinear aerodynamic effects that are critical for optimal design of the rotors. This is made possible by developing sparse, invertible neural networks (INNs) for inverse design and optimization that realize a 100x cost reduction compared to adjoint-based computational fluid dynamics (CFD) approaches, while enabling increased robustness of the final design. We demonstrate the INN tool for design of a section of the NREL 5-MW blade. All generated shapes satisfy the desired aerodynamic characteristics, demonstrating the success of the INN approach for inverse design of wind turbine blades.

aerodynamic design↗

AutoTandemML: Active Learning Enhanced Tandem Neural Networks for Inverse Design Problems

Inverse design in science and engineering involves determining optimal design parameters that achieve desired performance outcomes, a process often hindered by the complexity and high dimensionality of design spaces, leading to significant computational costs. To tackle this challenge, we propose a novel hybrid approach that combines active learning with Tandem Neural Networks to enhance the efficiency and effectiveness of solving inverse design problems. Active learning allows to selectively sample the most informative data points, reducing the required dataset size without compromising accuracy. We investigate this approach using three benchmark problems: airfoil inverse design, photonic surface inverse design, and scalar boundary condition reconstruction in diffusion partial differential equations. We demonstrate that integrating active learning with Tandem Neural Networks outperforms standard approaches across the benchmark suite, achieving better accuracy with fewer training samples.

97 MATHEMATICS AND COMPUTING↗

Multi-Electron Transfer at H-Terminated p-Si Electrolyte Interfaces: Large Photovoltages under Inversion Conditions

Here, photovoltages for hydrogen terminated p-Si(111) in an acetonitrile electrolyte were quantified with methyl viologen [1,1'-(CH 3 ) 2 -4,4'-bipyridinium](PF 6 ) 2 , abbreviated MV 2+ , and [Ru(bpy) 3 ](PF 6 ) 2 , where bpy is 2,2'-bipyridine, that respectively undergo two and three one-electron transfer reductions. The reduction potentials, E°, of the two MV 2+ reductions occurred at energies within the forbidden bandgap, while the three [Ru(bpy) 3 ] 2+ reductions occurred within the continuum of conduction band states. Bandgap illumination resulted in reduction that were more positive than that measured with a degenerately doped n + -Si demonstrative of a photovoltage, Vph, that increased in the order: MV 2+/+ (260 mV) < MV +/0 (400 mV) < Ru 2+/+ (530 mV) ~ Ru +/0 (540 mV) ~ Ru 0/– (550 mV). Pulsed 532 nm excitation generated electron-hole pairs whose dynamics were nearly constant under depletion conditions and increased markedly as the potential was raised or lowered. A long wavelength absorption feature assigned to conduction band electrons provided additional evidence for the presence of an inversion layer. Collectively, the data reveal that the most optimal photovoltage, as well at the longest electron-hole pair lifetime and the highest surface electron concentration, occur when E° lies energetically within the unfilled conduction band states where an inversion layer is present. The bell-shaped dependence for electron-hole pair recombination with the surface potential was predicted by the time-honored SRH model providing a clear indication that this interface provides access to all four bias conditions, i.e., accumulation, flat band, depletion, and inversion. The implications of these findings for photocatalysis applications and solar energy conversion are discussed.

14 SOLAR ENERGY↗

Engine/airframe integration for waverider cruise vehicles

Engine-airframe integrated hypersonic waveriders have been optimized previously for both cruising flight and accelerators using inverse design flowfields to provide desirable aerodynamic performance and inlet properties. With engine installation, cruisers retain the high L/D advantages of the basic waverider form. To provide an understanding of waverider optimizer trends, an analytical development is presented which includes the presentation of a simplified model of a generic waverider shape. The model is shown to have good agreement with computational results, and it is used to demonstrate engineering tradeoffs in L/D, volume, and volume efficiency. Based on the demonstrated performance trends, and considering the desirability of uniform inlet flow, an improved waverider shape for engine-airframe integration is proposed using a hybrid cone-wedge generating flowfield.

Lewis, Mark J.↗

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗

A hybrid data-driven and model-based approach for computationally efficient stochastic unit commitment and economic dispatch under wind and solar uncertainty

Stochastic unit commitment (UC) and economic dispatch (ED) are imperative in dealing with uncertainty in renewable forecast for power system operation and planning such that the overall expected production cost is minimized over the planning horizon. However, accurate calculation of the expected production cost requires assessment of a very large number of different scenarios of uncertain renewable resources, such as solar and wind, which is practically infeasible to simulate in real time. This article proposes a hybrid datadriven and physics-based model-predictive paradigm to efficiently solve for stochastic unit commitment and economic dispatch considering uncertainty in wind and solar power forecasts. Here, the novelty of the approach lies in decoupling the production cost estimation from the unit commitment and economic dispatch optimization problems under uncertainty without compromising on the fidelity of the solutions. A data-driven machine learning model is first developed to predict the mean optimal production cost. A physics-based inverse problem is then solved to get the stochastic UC and ED profiles from the expected cost. The presented approach considers, for the first time, solar uncertainty in UC/ED determination and enables efficient and accurate propagation of wind and solar uncertainty to estimate the statistics of the production cost. The effectiveness of the developed approach is demonstrated systematically on a stylized RTS-GMLC single-node system. The overall framework predicts the expected cost 62.5% more accurately than the existing state-of-the-art, on unforeseen days during the entire year, and yields, for the first time, the associated physically consistent UC and ED profiles. The solutions are also shown to be flexible in providing adequate daily reserves to address any statistical deviations from probabilistic power forecasts. The computational time associated with the presented method is only about 10 s compared to over 24 h needed for a conventional stochastic UC/ED determination under uncertainty on an Intel Core i9 processor with 32 GB of RAM.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Anisotropic colloidal interactions & assembly in AC electric fields

In this work, we match experimental and simulated configurations of anisotropic epoxy colloidal particles in high frequency AC electric fields by identifying analytical potentials for dipole-field and dipole-dipole interactions. We report an inverse Monte Carlo simulation algorithm to determine optimal fits of analytical potentials by matching simulated and experimental distribution functions for non-uniform liquid, liquid crystal, and crystal microstructures in varying amplitude electric fields. Two potentials that include accurate particle volume and dimensions along with a concentration dependent prefactor quantitatively capture experimental observations. At low concentrations, an effective ellipsoidal point dipole potential works well, whereas a novel stretched point dipole potential is found to be suitable at all concentrations, field amplitudes, and degrees of ordering. The simplicity, accuracy, and adjustability of the stretched point dipole potential suggest it can be applied to model field mediated microstructures and assembly of systematically varying anisotropic particle shapes.

42 ENGINEERING↗

Convolutional Neural Network Surrogate Models for the Mechanical Properties of Periodic Structures

The work in this paper describes neural network surrogate models for calculating the effective mechanical properties of a periodic composites. The models achieve good accuracy even when only provided with training data sampling a small portion of the design space. As an example, the surrogate models are applied to solving the inverse design problem of finding structures with optimal mechanical properties. The surrogate models are sufficiently accurate to recover optimal solutions in general agreement with established topology optimization methods. However, improvements will be required to develop robust, efficient neural network-based surrogate models and several directions for future research are highlighted here.

42 ENGINEERING↗

Optimal low-thrust takeoff from an orbit about an oblate planet

Future space missions to the outer planets may depend upon the use of low-thrust propulsion systems. As these planets are decidedly oblate, the question of the effect of that oblateness on a low-thrust trajectory is of some interest. In this paper the problem of optimal energy increase is attacked under the assumption that the coefficients for the second zonal harmonic, and the nondimensional thrust acceleration are the same order of magnitude. By means of a two-variable asymptotic expansion technique, a near optimal control program is generated and the first-order uniformly valid approximation for the corresponding trajectory is obtained. Tangential thrust is shown to be a good near-optimal thrust program even in the presence of oblateness effects. The optimal control program is found to be oscillatory and quite similar to the optimal control for energy increase in an inverse square gravitational field.

Jacobson, R. A.↗

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗

Design of a mixer for the thrust-vectoring system on the high-alpha research vehicle

One of the advanced control concepts being investigated on the High-Alpha Research Vehicle (HARV) is multi-axis thrust vectoring using an experimental thrust-vectoring (TV) system consisting of three hydraulically actuated vanes per engine. A mixer is used to translate the pitch-, roll-, and yaw-TV commands into the appropriate TV-vane commands for distribution to the vane actuators. A computer-aided optimization process was developed to perform the inversion of the thrust-vectoring effectiveness data for use by the mixer in performing this command translation. Using this process a new mixer was designed for the HARV and evaluated in simulation and flight. An important element of the Mixer is the priority logic, which determines priority among the pitch-, roll-, and yaw-TV commands.

Pahle, Joseph W.↗

Intelligent Control Approaches for UAVs

This talk will present overviews of various intelligent control technologies currently being developed and studied at the NASA Ames Research Center as applicable to Unmanned Aerial Vehicles (UAVs), Mars flyers, and to the next generation of flight controllers for manned aircraft. The approaches being examined include: (a) direct adaptive dynamic inverse controller, (b) adaptive critic-based optimal trajectory generator; (c) optimal allocation technique based on linear programming; (4) immunized maneuvering using autopilot building blocks. These approaches can utilize, but do not require, fault detection and isolation information. Piloted and unmanned simulation studies are performed to examine if the intelligent flight control techniques adequately: 1) match flying qualities of modem fly-by-wire flight controllers under nominal conditions; 2) improve performance under failure conditions when sufficient control authority is available; and 3) achieve intelligent maneuvering capabilities for unmanned vehicles. Results obtained so far will be presented and discussed.

KrisnaKumar, Kalmanje↗