Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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,405 records · Page 78

mystic : software for autonomous discovery and design under uncertainty

Throughout the diverse range of science and engineering applications, there is a growing desire to develop computational methods that can reliably predict the behavior of complex systems. Specifically, there is a strategic need for tools that can robustly forecast the behavior of complex physical systems, where data may be high-dimensional, noisy, or sparse, and models of the system may be time-dependent or include uncertainty. We use mystic to build tools that leverage statistical learning, physics-informed learning, and active learning in the efficient generation of reliably predictive surrogates for complex physical systems. mystic is a robust, proven, open-source optimization and uncertainty quantification toolkit with over a decade of use in the design and optimization of neutron instrumentation, solar-powered drones, and gasguns, and in iterative tuning of models for Raman spectroscopy and elastoplastic materials strength. Recent developments have focused on automated learning of statistically robust surrogates under uncertainty, with applications in materials in extreme environments, nanostructures, materials simulations and strength models, and the failure of shielding under particle radiation. In 2020, McKerns demonstrated active learning of optimally robust surrogates with respect to new simulated data for molecular dynamics simulations of materials mixing in warm dense matter, and is currently applying active learning to the automated steering of particle accelerator beams and the optimal design and control of quantum optical sensor instrumentation.

42 ENGINEERING↗

Mitigation of Motor Stalling and FIDVR via Energy Storage Systems with Temporal Logic Specifications

The fault-induced delayed voltage recovery (FIDVR) phenomenon has been very common from the distribution system through the transmission system. It causes a delay on recovering significantly depressed local voltage after the fault is cleared, and it can also lead to more widespread cascading system failures. Mitigating this event with current control approaches is challenging and becoming a crucial issue. Here, a model predictive control-based strategy employing signal temporal logic specifications is proposed to help mitigate FIDVR. To this end, it investigates and extends a dynamic performance model allowing analytic insights into the system-wide impact of motor stalling and FIDVR. The proposed controller provides richer descriptions of voltage specifications addressing both magnitude and time simultaneously. We consider different control specifications with reactive power support from energy storage systems to prevent the voltage during/after the fault from dropping too low, and reduce the delay time of voltage recovery. The simulation results conducted with the IEEE 57 bus test network validate the proposed method and demonstrate the effectiveness of the mitigation strategy on FIDVR.

25 ENERGY STORAGE↗

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Prediction of pilot-aircraft stability boundaries and performance contours

Control-theoretic pilot models can provide important new insights regarding the stability and performance characteristics of the pilot-aircraft system. Optimal-control pilot models can be formed for a wide range of flight conditions, suggesting that the human pilot can maintain stability if he adapts his control strategy to the aircraft's changing dynamics. Of particular concern is the effect of sub-optimal pilot adaptation as an aircraft transitions from low to high angle-of-attack during rapid maneuvering, as the changes in aircraft stability and control response can be extreme. This paper examines the effects of optimal and sub-optimal effort during a typical 'high-g' maneuver, and it introduces the concept of minimum-control effort (MCE) adaptation. Limited experimental results tend to support the MCE adaptation concept.

Stengel, R. F.↗

Automatic Generation of Algorithms for the Statistical Analysis of Planetary Nebulae Images

Analyzing data sets collected in experiments or by observations is a Core scientific activity. Typically, experimentd and observational data are &aught with uncertainty, and the analysis is based on a statistical model of the conjectured underlying processes, The large data volumes collected by modern instruments make computer support indispensible for this. Consequently, scientists spend significant amounts of their time with the development and refinement of the data analysis programs. AutoBayes [GF+02, FS03] is a fully automatic synthesis system for generating statistical data analysis programs. Externally, it looks like a compiler: it takes an abstract problem specification and translates it into executable code. Its input is a concise description of a data analysis problem in the form of a statistical model as shown in Figure 1; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. Internally, however, it is quite different: AutoBayes derives a customized algorithm implementing the given model using a schema-based process, and then further refines and optimizes the algorithm into code. A schema is a parameterized code template with associated semantic constraints which define and restrict the template s applicability. The schema parameters are instantiated in a problem-specific way during synthesis as AutoBayes checks the constraints against the original model or, recursively, against emerging sub-problems. AutoBayes schema library contains problem decomposition operators (which are justified by theorems in a formal logic in the domain of Bayesian networks) as well as machine learning algorithms (e.g., EM, k-Means) and nu- meric optimization methods (e.g., Nelder-Mead simplex, conjugate gradient). AutoBayes augments this schema-based approach by symbolic computation to derive closed-form solutions whenever possible. This is a major advantage over other statistical data analysis systems which use numerical approximations even in cases where closed-form solutions exist. AutoBayes is implemented in Prolog and comprises approximately 75.000 lines of code. In this paper, we take one typical scientific data analysis problem-analyzing planetary nebulae images taken by the Hubble Space Telescope-and show how AutoBayes can be used to automate the implementation of the necessary anal- ysis programs. We initially follow the analysis described by Knuth and Hajian [KHO2] and use AutoBayes to derive code for the published models. We show the details of the code derivation process, including the symbolic computations and automatic integration of library procedures, and compare the results of the automatically generated and manually implemented code. We then go beyond the original analysis and use AutoBayes to derive code for a simple image segmentation procedure based on a mixture model which can be used to automate a manual preproceesing step. Finally, we combine the original approach with the simple segmentation which yields a more detailed analysis. This also demonstrates that AutoBayes makes it easy to combine different aspects of data analysis.

Fischer, Bernd↗

The evaluation of several agility metrics for fighter aircraft using optimal trajectory analysis

Several functional agility metrics, including the combat cycle time metric, dynamic speed turn plots, and relative energy state metric, are used to compare turning performance for generic F-18, X-29, and X-31-type aircraft models. These three-degree-of-freedom models have characteristics similar to the real aircraft. The performance comparisons are made using data from optimal test trajectories to reduce sensitivities to different pilot input techniques and to reduce the effects of control system limiters. The turn performance for all three aircraft is calculated for simulated minimum time 180 deg heading captures from simulation data. Comparisons of the three aircraft give more insight into turn performance than would be available from traditional measures of performance. Using the optimal test technique yields significant performance improvements as measured by the metrics. These performance improvements were found without significant increases in turn radius.

Ryan, George W., III↗

Optimal inputs for system identification

Identification criteria are presented for linear dynamic systems with and without process noise. With process noise, the state equations are replaced by the Kalman filter equations. If the identification performance index is expanded in a Taylor's series with respect to the parameters to be identified, then maximizing the weighting factor of the quadratic term with respect to the inputs will insure that an identification algorithm will converge more rapidly and to a more accurate result than with nonoptimal inputs. The expectation of this weighting factor is the Fisher information matrix, and its inverse is a lower bound for the covariance of the parameters. Direct and indirect methods of calculating the information matrix are presented for systems with and without process noise.

Reid, D. B.↗

Design of helicopter rotor blades for optimum dynamic characteristics

The mass and stiffness distributions for helicopter rotor blades are tailored in such a way to give a predetermined placement of blade natural frequencies. The optimal design is pursued with respect of minimum weight, sufficient inertia, and reasonable dynamic characteristics. Finite element techniques are used as a tool. Rotor types include hingeless, articulated, and teetering.

Peters, D. A.↗

An Automatic Medium to High Fidelity Low-Thrust Global Trajectory Toolchain; EMTG-GMAT

Solving the global optimization, low-thrust, multiple-flyby interplanetary trajectory problem with high-fidelity dynamical models requires an unreasonable amount of computational resources. A better approach, and one that is demonstrated in this paper, is a multi-step process whereby the solution of the aforementioned problem is solved at a lower-fidelity and this solution is used as an initial guess for a higher-fidelity solver. The framework presented in this work uses two tools developed by NASA Goddard Space Flight Center: the Evolutionary Mission Trajectory Generator (EMTG) and the General Mission Analysis Tool (GMAT). EMTG is a medium to medium-high fidelity low-thrust interplanetary global optimization solver, which now has the capability to automatically generate GMAT script files for seeding a high-fidelity solution using GMAT's local optimization capabilities. A discussion of the dynamical models as well as thruster and power modeling for both EMTG and GMAT are given in this paper. Current capabilities are demonstrated with examples that highlight the toolchains ability to efficiently solve the difficult low-thrust global optimization problem with little human intervention.

low thrust↗

Efficient sensitivity analysis and optimization of a helicopter rotor

Aeroelastic optimization of a system essentially consists of the determination of the optimum values of design variables which minimize the objective function and satisfy certain aeroelastic and geometric constraints. The process of aeroelastic optimization analysis is illustrated. To carry out aeroelastic optimization effectively, one needs a reliable analysis procedure to determine steady response and stability of a rotor system in forward flight. The rotor dynamic analysis used in the present study developed inhouse at the University of Maryland is based on finite elements in space and time. The analysis consists of two major phases: vehicle trim and rotor steady response (coupled trim analysis), and aeroelastic stability of the blade. For a reduction of helicopter vibration, the optimization process requires the sensitivity derivatives of the objective function and aeroelastic stability constraints. For this, the derivatives of steady response, hub loads and blade stability roots are calculated using a direct analytical approach. An automated optimization procedure is developed by coupling the rotor dynamic analysis, design sensitivity analysis and constrained optimization code CONMIN.

Lim, Joon W.↗

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗

Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river

Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗

Chamber Flux and Porewater Concentration of CH4, CO2 and N2O, 2018, Columbia River bank at the Hanford site, WA, USA

This is the corresponding observations data for the paper: Jorge A Villa; Garrett J Smith; Yang Ju; Lupita Renteria; Jordan C Angle; Evan Arntzen; Samuel F Harding; Huiying Ren; Xingyuan Chen; Audrey H Sawyer; Emily B Graham; James C Stegen; Kelly C Wrighton; Gil Bohrer (2020) Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river. Science of the Total Environment.Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow, intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during the low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗

Parameter identification and sensitivity analysis for a robotic manipulator arm

The development of a nonlinear dynamic model for large oscillations of a robotic manipulator arm about a single joint is described. Optimization routines are formulated and implemented for the identification of electrical and physical parameters from dynamic data taken from an industrial robot arm. Special attention is given to difficulties caused by the large sensitivity of the model with respect to unknown parameters. Performance of the parameter identification algorithm is improved by choosing a control input that allows actuator emf to be included in an electro-mechanical model of the manipulator system.

Brewer, D. W.↗

Utilization of Existing Pipelines in Hydrogen Transport: Literature Review Report

This report critically reviews the flow behavior of hydrogen-natural gas (H 2 -NG) mixtures in pipelines and examines the critical factors of hydrogen integration into existing natural gas infrastructure. It addresses the choking behavior characterized by velocity increase and pressure drop, as well as the effects of flow restrictions and pressure losses during hydrogen transport. Computational and analytical models are used to investigate these effects, and their effects on thermodynamic properties and system performance are evaluated. The study also reviews the energy efficiency and flow dynamics of hydrogen and methane-hydrogen mixtures and optimizes the hydrogen flow rate. In addition, the effects of these mixtures on the flow characteristics are discussed in detail, with special emphasis on the compressibility factor (z factor) and fluid properties based on equations of state for hydrogen-natural gas mixtures. The study also analyzes the mixture ratios and highlights the thermophysical properties, flow dynamics, and hydrogen-blended natural gas application potential. These investigations assess flow stability, material interactions, and operational feasibility of transporting hydrogen mixtures through natural gas pipelines, which contribute to developing sustainable and efficient energy systems.

08 HYDROGEN↗

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rotorcraft Optimization Tools: Incorporating Rotorcraft Design Codes into Multi-Disciplinary Design, Analysis, and Optimization

One of the goals of NASA's Revolutionary Vertical Lift Technology Project (RVLT) is to provide validated tools for multidisciplinary design, analysis and optimization (MDAO) of vertical lift vehicles. As part of this effort, the software package, RotorCraft Optimization Tools (RCOTOOLS), is being developed to facilitate incorporating key rotorcraft conceptual design codes into optimizations using the OpenMDAO multi-disciplinary optimization framework written in Python. RCOTOOLS, also written in Python, currently supports the incorporation of the NASA Design and Analysis of RotorCraft (NDARC) vehicle sizing tool and the Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics II (CAMRAD II) analysis tool into OpenMDAO-driven optimizations. Both of these tools use detailed, file-based inputs and outputs, so RCOTOOLS provides software wrappers to update input files with new design variable values, execute these codes and then extract specific response variable values from the file outputs. These wrappers are designed to be flexible and easy to use. RCOTOOLS also provides several utilities to aid in optimization model development, including Graphical User Interface (GUI) tools for browsing input and output files in order to identify text strings that are used to identify specific variables as optimization input and response variables. This paper provides an overview of RCOTOOLS and its use

Analysi↗