Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “trajectory 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,693 records · Page 94

HIPED: Machine learning framework for spherical tokamak pedestal prediction and optimization

We introduce a Machine Learning framework, HIPED (HeIght and width Predictor for Edge Dynamics), for predicting and optimizing pedestal and core performance in spherical tokamak plasmas. Trained on pedestal and core datasets from the third MAST-U campaign, HIPED provides accurate estimates of pedestal height and width. The results reveal notable differences compared with conventional aspect-ratio studies; for instance, a simple power-law relation between pedestal width and height has very low accuracy. Instead, additional parameters such as normalized plasma pressure, elongation, and Greenwald fraction significantly improve accuracy. HIPED can also be trained only on `control room parameters' to inform experimentalists of which controllable parameters to adjust for improving core-integrated performance. The framework further includes a multi-objective optimization scheme that helps guide experimental planning and optimization. We find Pareto-optimal discharges with respect to various features, including distance from edge-localized modes and normalized plasma pressure, track their parameter trajectories over time, and identify the control room parameters required for these Pareto-optimal discharges. This provides a framework for systematically optimizing core and edge performance according to different experimental priorities.

Parisi, Jason F. [Princeton Plasma Physics Laborat↗

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↗

Multi-Objective Hybrid Optimal Control for Interplanetary Mission Planning

Preliminary design of low-thrust interplanetary missions is a highly complex process. The mission designer must choose discrete parameters such as the number of flybys, the bodies at which those flybys are performed, and in some cases the final destination. Because low-thrust trajectory design is tightly coupled with systems design, power and propulsion characteristics must be chosen as well. In addition, a time-history of control variables must be chosen which defines the trajectory. There are often may thousands, if not millions, of possible trajectories to be evaluated. The customer who commissions a trajectory design is not usually interested in a point solution, but rather the exploration of the trade space of trajectories between several different objective functions. This can be a very expensive process in terms of the number of human analyst hours required. An automated approach is therefore very desirable. This work presents such an approach by posing the mission design problem as a multi-objective hybrid optimal control problem. The method is demonstrated on hypothetical mission to the main asteroid belt and to Deimos.

Low Thrust↗

Exploring Transfers Between Earth-Moon Halo Orbits via Multi-Objective Reinforcement Learning

Multi-Reward Proximal Policy Optimization, a multi-objective deep reinforcement learning algorithm, is used to examine the design space of low-thrust trajectories for a SmallSat transferring between two libration point orbits in the Earth-Moon system. Using Multi-Reward Proximal Policy Optimization, multiple policies are simultaneously and efficiently trained on three distinct trajectory design scenarios. Each policy is trained to create a unique control scheme based on the trajectory design scenario and assigned reward function: a unique combination of weights scaling competing objectives that guide the spacecraft to the target mission orbit, incentivize faster flight times, and penalize propellant mass usage. Then, the policies are evaluated on the same set of perturbed initial conditions in each scenario to generate the propellant mass usages, flight times, and state discontinuities from a reference trajectory for each control scheme. This solution space of low-thrust trajectories for a SmallSat is used to examine the multi-objective trade space for the trajectory design scenario. By autonomously constructing the solution space, insights into the required propellant mass, flight time, and transfer geometry are rapidly achieved.

Christopher J Sullivan↗

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation that is based on a coupling between inertial navigation and flush air data sensing pressure measurements. In this approach, the full navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to directly estimate atmospheric winds and density using a nonlinear weighted least-squares algorithm. The approach uses a high fidelity model of atmosphere stored in table-look-up form, along with simplified models of that are propagated along the trajectory within the algorithm to provide prior estimates and covariances to aid the air data state solution. Thus, the method is essentially a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere and winds are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the discrete-time observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content to the system. The algorithm is then applied to the design of the pressure measurement system for the Mars 2020 mission. The pressure port layout is optimized to maximize the observability of atmospheric states along the trajectory. Linear covariance analysis is performed to assess estimator performance for a given pressure measurement uncertainty. The results indicate that the new tightly-coupled estimator can produce enhanced estimates of atmospheric states when compared with existing algorithms.

Karlgaard, Christopher D.↗

Computation of optimal feedback strategies for interception in a horizontal plane

The problem of minimum-time interception of a target moving in a horizontal plane is studied. The target may fly a trajectory that is known at the start of the interception or try to evade the intercepting aircraft. The interceptor's optimal motion is described independently of the target in terms of an extremal trajectory map. The latter is used to develop guidelines for suboptimal approximations to the extremals. A method of constructing the feedback solution by drawing isochrones (constant minimum-time loci) is developed. Sections of the feedback solution for the interception of a straight-flying target are presented, and the construction of an isochrone for the interception of an actively evading target is demonstrated.

Rajan, N.↗

Path Planning: Differential Dynamic Programming and Model Predictive Path Integral Control on VTOL Aircraft

This paper explores two optimal control approaches, widely used in robotics, to establish their viability as real-time trajectory planners for vehicle configurations envisioned for the emerging aviation sector of Urban Air Mobility (UAM). Differential Dynamic Programming (DDP) enables planning over highly nonlinear dynamics using second-order approximations along a nominal trajectory, and displays quadratic convergence to a local solution. Model Predictive Path Integral (MPPI) is a stochastic sampling-based algorithm that can optimize for general cost criteria, including potentially highly nonlinear formulations, and supports parallel computation through the use of modern GPU hardware. In this work, DDP and MPPI were implemented using model predictive control (MPC), and the results indicate they are able to successfully transition the aircraft over different flight envelopes and generate trajectories unique to UAM vehicles.

Differential Dynamic Programming↗

Influence of pilot-fuel mixing on the spatio-temporal progression of two-stage autoignition of diesel-sprays in low-reactivity ambient fuel-air mixture

The spatial and temporal locations of autoignition for direct-injection compression-ignition engines depend on fuel chemistry, temperature, pressure, and mixing trajectories in the fuel jets. Dual-fuel systems can provide insight into both fuel-chemistry and physical effects by varying fuel reactivities and engine operating conditions. In this context, the spatial and temporal progression of two-stage autoignition of a diesel-fuel surrogate, n-heptane, in a lean-premixed charge of synthetic natural-gas (NG) and air is imaged in an optically accessible heavy-duty diesel engine. The lean-premixed charge of NG is prepared by fumigation upstream of the engine intake manifold. Optical diagnostics include high-speed (15kfps) cool-flame chemiluminescence-imaging as an indicator of low-temperature heat-release (LTHR) and OH* chemiluminescence-imaging as an indicator high-temperature heat-release (HTHR). NG prolongs the ignition delay of the pilot fuel and increases the combustion duration. Zero-dimensional chemical-kinetics simulations provide further understanding by replicating a Lagrangian perspective for mixtures evolving along streamlines originating either at the fuel nozzle or in the ambient gas, for which the pilot-fuel concentration is either decreasing or increasing, respectively. The zero-dimensional simulations predict that LTHR initiates most likely on the air streamlines before transitioning to HTHR, either on fuel-streamlines or on air-streamlines in regions of near-constant φ. Due to the relatively short pilot-fuel injection-durations, the transient increase in entrainment near the end of injection (entrainment wave) is important for quickly creating auto-ignitable mixtures. To achieve desired combustion characteristics, e.g., multiple ignition-kernels and favorable combustion phasing and location (e.g., for reducing wall heat-transfer or optimizing charge stratification), adjusting injection parameters could tailor mixing trajectories to offset changes in fuel ignition chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HP-FLEX MPC v0.1.0

HP-FLEX MPC is control software developed by Lawrence Berkeley National Laboratory with support from the California Energy Commission (CEC) through EPIC-19-301. HP-FLEX aims to provide load flexibility for heat pumps (HPs) in response to dynamic grid signals (including Time-of-Use, Dynamic Pricing, and Critical Peak Pricing) while maintaining thermostat temperatures within user-specified bounds. The software includes a system-identification module, which models the dynamics of the building envelope with thermostat data, and a control module based on a model predictive controller (MPC) to make optimal decisions. HP-FLEX receives forecasts of outdoor air temperature, solar irradiation, and internal gain (if available), as well as trajectories of energy price, temperature lower and upper bounds over a prediction horizon. It then optimizes heating and cooling capacities to minimize energy cost and peak power (with a user-defined weight on peak power) over the prediction horizon, while maintaining room air temperature within the temperature constraints, and outputs the optimal thermostat setpoints.

Kim, Donghun↗

Performance Impacts to the NASA Artemis II Trajectory Correction Burn Placement

As NASA embarks to return humans to the lunar vicinity with the upcoming Artemis II mission, the selected free return trajectory taking the crew to the Moon in the Orion spacecraft is impacted by the execution of small trajectory correction burns to ensure the spacecraft stays on course for a successful return to Earth. The placement of these nominally zero translational maneuvers must account for the crew schedule, navigation tracking constraints, spacecraft venting, thermal and communication requirements, and a host of other programmatic factors. Understanding the influence the placement of these periodic burn corrections have on the integrated GN\&C performance can provide valuable insight to mission controllers and trajectory planning processes to untangle the complex trade space considered for both baseline and contingency scenarios. This sensitivity information can also be utilized to facilitate the optimized placement of these burns. This paper utilizes several techniques to systematically generate the performance impacts to the NASA Artemis II trajectory correction burn placement and demonstrate how to derive optimized locations that make the system robust to crew activity, maneuver execution errors, navigation uncertainty, orbit insertion errors, disturbance accelerations, and other system limitations.

GN&C↗

Wassersplines for Neural Vector Field-Controlled Animation

Much of computer-generated animation is created by manipulating meshes with rigs. While this approach works well for animating articulated objects like animals, it has limited flexibility for animating less structured free-form objects. Here we introduce Wassersplines, a novel trajectory inference method for animating unstructured densities based on recent advances in continuous normalizing flows and optimal transport. The key idea is to train a neurally-parameterized velocity field that represents the motion between keyframes. Trajectories are then computed by advecting keyframes through the velocity field. We solve an additional Wasserstein barycenter interpolation problem to guarantee strict adherence to keyframes. Our tool can stylize trajectories through a variety of PDE-based regularizers to create different visual effects. We demonstrate our tool on various keyframe interpolation problems to produce temporally-coherent animations without meshing or rigging.

97 MATHEMATICS AND COMPUTING↗

GEANT4 Simulations on Faraday Cup Design for PIP-II Laser Wire Scanner System

The PIP-II accelerator upgrade at Fermilab represents a groundbreaking leap forward in high-energy physics research. This ambitious initiative involves enhancing Fermilab's accelerator complex by replacing the current linear accelerator with a warm front end (WFE) capable of accelerating H- beams up to 2.1 MeV. Subsequently, a superconducting linac further accelerates these beams up to 800 MeV. To precisely measure the transverse beam profile, a combination of traditional wire scanners at the WFE section and Laser wire scanners along the superconducting linac are planned for implementation. This investigation centers on refining the Faraday cup design for the PIP-II Laser wire scanners by utilizing GEANT4, a Monte Carlo simulation toolkit. Leveraging this method enables a comprehensive analysis of particle trajectories, energy deposition, secondary electron emission, backscattering, etc., facilitating optimization through adjustments to cup geometries, materials, and placement to maximize its efficiency.

43 PARTICLE ACCELERATORS↗

The reverse Goddard problem.

Fuel-time optimal retrothrust control for vertical and gravity turn ballistic trajectories of soft landing nonlifting bodies based on Pontryagin principle

Bejczy, A. K.↗

On-line optimization of aircraft altitude and flight path angle dynamics

The optimization of aircraft altitude and flight path angle dynamics is addressed in a form suitable for on-line computation and control. The approach here is a direct extension of the work reported by Calise (1979), where singular perturbation methods were used to optimize position, energy and heading dynamics; it thus represents an optimal control solution that models all of the primary trajectory related dynamics. It is pointed out that the resulting algorithm can be regarded as a nonlinear feedback control law. The minimum time intercept of a fixed terminal point is used in setting the framework in which the analytical results are developed. The main theoretical result is that the dynamics, while not completely separable, can be approximated by singular perturbation methods when the control model includes relative position dynamics.

Calise, A. J.↗