Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

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 775 records · Page 43

Adaptive state estimation for control of flexible structures

This paper proposes a new approach of obtaining adaptive state estimation of a system in the presence of unknown system disturbances and measurement noise. In the beginning, a non-optimal Kalman filter with arbitrary initial guess for the process and measurement noises is implemented. At the same time, an adaptive transversal predictor (ATP) based on the recursive least-squares (RLS) algorithm is used to yield optimal one- to p- step-ahead output predictions using the previous input/output data. Referring to these optimal predictions the Kalman filter gain is updated and the performance of the state estimation is thus improved. If forgetting factor is implemented in the recursive least-squares algorithm, this method is also capable of dealing with the situation when the noise statistics are slowly time-varying. This feature makes this new approach especially suitable for the control of flexible structures. A numerical example demonstrates the feasibility of this real time adaptive state estimation method.

Chen, Chung-Wen↗

Robotic Recon for Human Exploration

Robotic reconnaissance has the potential to significantly improve scientific and technical return from lunar surface exploration. In particular, robotic recon may increase crew productivity and reduce operational risk for exploration. However, additional research, development and field-testing is needed to mature robot and ground control systems, refine operational protocols, and specify detailed requirements. When the new lunar surface campaign begins around 2020, and before permanent outposts are established, humans will initially be on the Moon less than 10% of the time. During the 90% of time between crew visits, robots will be available to perform surface operations under ground control. Understanding how robotic systems can best address surface science needs, therefore, becomes a central issue Prior to surface missions, lunar orbiters (LRO, Kaguya, Chandrayyan-1, etc.) will map the Moon. These orbital missions will provide numerous types of maps: visible photography, topographic, mineralogical and geochemical distributions, etc. However, remote sensing data will not be of sufficient resolution, lighting, nor view angle, to fully optimize pre-human exploration planning, e.g., crew traverses for field geology and geophysics. Thus, it is important to acquire supplemental and complementary surface data. Robotic recon can obtain such data, using robot-mounted instruments to scout the surface and subsurface at resolutions and at viewpoints not achievable from orbit. This data can then be used to select locations for detailed field activity and prioritize targets to improve crew productivity. Surface data can also help identify and assess terrain hazards, and evaluate alternate routes to reduce operational risk. Robotic recon could be done months in advance, or be part of a continuing planning process during human missions.

Deans, Matthew↗

Travels, optimizations, and ordering

With all the necessary build pathing in place from the preceding steps, the final task is to organize the build pathing via directionality, ordering, and connectivity. Directionality controls closed-loop pathing by deciding whether construction occurs clockwise or counterclockwise. Ordering determines whether construction occurs outside-in or inside-out, by manipulating what path types will be printed first. Connectivity links all the pathing together, which is achieved through the insertion of travels. Connectivity has numerous approaches, but only the most common will be addressed in this chapter. In addition to these manipulations, other considerations, such as spiral pathing and unique cases for infill, will be discussed.

Borish, Michael↗

A Time-Regularized Multiple Gravity-Assist Low-Thrust Bounded-Impulse Model for Trajectory Optimization

The multiple gravity assist low-thrust (MGALT) trajectory model combines the medium-fidelity Sims-Flanagan bounded-impulse transcription with a patched-conics flyby model and is an important tool for preliminary trajectory design. While this model features fast state propagation via Kepler's equation and provides a pleasingly accurate estimation of the total mass budget for the eventual flight-suitable integrated trajectory it does suffer from one major drawback, namely its temporal spacing of the control nodes. We introduce a variant of the MGALT transcription that utilizes the generalized anomaly from the universal formulation of Kepler's equation as a decision variable in addition to the trajectory phase propagation time. This results in two improvements over the traditional model. The first is that the maneuver locations are equally spaced in generalized anomaly about the orbit rather than time. The second is that the Kepler propagator now has the generalized anomaly as its independent variable instead of time and thus becomes an iteration-free propagation method. The new algorithm is outlined, including the impact that this has on the computation of Jacobian entries for numerical optimization, and a motivating application problem is presented that illustrates the improvements that this model has over the traditional MGALT transcription.

Optimization↗

Parametric and sensitivity analysis of a PCM-integrated wall for optimal thermal load modulation in lightweight buildings

Load modulation in buildings is becoming increasingly important due to growing disparity in energy demand during peak and off-peak hours. Integrating phase change material (PCM) in building envelopes and using a controlled precooling strategy can provide substantial thermal load modulation; however, it may greatly increase the total energy use. Previous studies have employed PCM in building envelopes primarily for energy savings and, to some extent, peak load shedding and shifting. However, the load modulation capacity of a PCM-integrated envelope has not been well explored in the literature. In this study, we perform an extensive parametric and sensitivity analysis on PCM-integrated lightweight building walls and examine the combinatory effects of various PCM parameters on thermal load modulation and wall-related heat gains in buildings. Using numerical simulations, we investigate eight PCM parameters: PCM location in the wall, transition temperature, thickness, latent heat, transition range, density, specific heat, and thermal conductivity. Here, we evaluate their impact and relative importance to achieve maximum load modulation in buildings without compromising occupants’ thermal comfort or total energy use. The results show that the optimized PCM proposed in this study can completely invert the transient heat gain profile of the wall, providing up to 70% reduction of wall-related heat gain during peak hours without a major increase in the cumulative heat gain.

42 ENGINEERING↗

Investigations into Balancing Peak-to-Average Power Ratio and Mean Power Extraction for a Two-Body Point-Absorber Wave Energy Converter

The power harnessed by wave energy converters (WECs) in oceans is highly variable and, therefore, has a high peak-to-average power (PTAP) ratio. To minimize the cost of a WEC power take off (PTO) system, it is desirable to reduce the PTAP ratio while maximizing the mean power extracted by WECs. The important issue of how PTAP ratio reduction measures (such as adding an inertia element) can affect the mean power extracted in a reference model has not been thoroughly addressed in the literature. To investigate this correlation, this study focuses on the integration of the U.S. Department of Energy’s Reference Model 3, a two-body point absorber, with a slider-crank WEC for linear-to-rotational conversion. In the first phase of this study, a full-scale numerical model was developed that predicts how PTO system parameters, along with an advanced control algorithm, can potentially affect the proposed WEC’s PTAP ratio as well as the mean power extracted. In the second phase, an appropriate scaled-down model was developed, and extracted power results were successfully validated against the full-scale model. Finally, numerical and hardware-in-the-loop (HIL) simulations based on the scaled-down model were designed and conducted to optimize or make trade-offs between the operational performance and PTAP ratio. The initial results with numerical and HIL simulations reveal that gear ratio, crank radius, and generator parameters substantially impact the PTAP ratio and mean power extracted.

16 TIDAL AND WAVE POWER↗

Contributions of the Transonic Dynamics Tunnel to the Testing of Active Control of Aeroelastic Response

By the 1960s, researchers began to investigate the feasibility of using active controls technology (ACT) for increasing the capabilities of military and commercial aircraft. Since then many researchers, too numerous to mention, have investigated and demonstrated the usefulness of ACT for favorably modifying the aeroelastic response characteristics of flight vehicles. As a result, ACT entered the limelight as a viable tool for answering some very difficult design questions and had the potential for obtaining structural weight reductions optimizing maneuvering performance, and satisfying the multimission requirements being imposed on future military and commercial aircraft designs. Over the past 40 years, the NASA Langley Research Center (LaRC) has played a major role in developing ACT in part by its participation in many wind-tunnel programs conducted in the Transonic Dynamics Tunnel (TDT). These programs were conducted for the purposes of: (1) establishing concept feasibility; (2) demonstrating proof of concept; and (3) providing data for validating new modeling, analysis, and design methods. This paper provides an overview of the ACT investigations conducted in the TDT. For each program discussed herein, the objectives of the effort, the testing techniques, the test results, any, signIficant findings, and the lessons learned with respect to ACT testing are presented.

Perry, Boyd, III↗

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗

Efficient, direct compilation of SU(N) operations into SNAP & Displacement gates

We present a function which connects the parameter of a previously published short sequence of selective number-dependent arbitrary phase (SNAP) and displacement gates acting on a qudit encoded into the Fock states of a superconducting cavity, $V_k(\alpha)=D(\alpha)R_\pi(k)D(-2\alpha)R_\pi(k)D(\alpha)$ to the angle of the Givens rotation $G(\theta)$ on levels $|k\rangle,|k+1\rangle$ that sequence approximates, namely $\alpha=\Phi(\theta) = \frac{\theta}{4\sqrt{k+1}}$. Previous publications left the determination of an appropriate $\alpha$ to numerical optimization at compile time. The map $\Phi$ gives us the ability to compile directly any $d$-dimensional unitary into a sequence of SNAP and displacement gates in $O(d^3)$ complex floating point operations with low constant prefactor, avoiding the need for numerical optimization. Numerical studies demonstrate that the infidelity of the generated gate sequence $V_k$ per Givens rotation $G$ scales as approximately $O(\theta^6)$. We find numerically that the error on compiled circuits can be made arbitrarily small by breaking each rotation into $m$$\theta/m$ rotations, with the full $d\times d$ unitary infidelity scaling as approximately $O(m^{-4})$. This represents a significant reduction in the computational effort to compile qudit unitaries either to SNAP and displacement gates or to generate them via direct low-level pulse optimization via optimal control.

Job, Joshua↗

Numerical analysis of flow through scramjet engine inlets

A set of computer programs has been developed to analyze flow through supersonic combustion ramjet (scramjet) inlets. These programs solve either the two- or three-dimensional Euler/Navier-Stokes equations in full conservation form by MacCormack's explicit or explicit-implicit method. An algebraic two-layer eddy viscosity model is used for turbulent flow calculations. The programs are operational on Control Data CYBER-200 series vector-processing computer system and have been optimized to take maximum advantage of the vector processing capability of the system. Since their development, the programs have been extensively verified and used to analyze a number of very complex inlet configurations. In this paper, results are presented from two-dimensional, quasi-three-dimensional, and three-dimensional analyses of the inlet flow field to illustrate the use of the programs.

Kumar, Ajay↗

An explicit solution to the optimal LQG problem for flexible structures with collocated rate sensors

We present a class of compensators in explicit form (not requiring numerical computer calculations) for stabilizing flexible structures with collocated rate sensors. They are based on the explicit solution, valid for both Continuum and FEM Models, of the LQG problem for minimizing mean square rate. They are robust with respect to system stability (will not destabilize modes even with mismatch of parameters), can be instrumented in state space form suitable for digital controllers, and can be specified directly from the structure modes and mode 'signature' (displacement vectors at sensor locations). Some simulation results are presented for the NASA LaRC Phase-Zero Evolutionary Model - a modal Trust model with 86 modes - showing damping ratios attainable as a function of compensator design parameters and complexity.

Balakrishnan, A. V.↗

Evaluation of the Energy, Hygrothermal, and Thermal Capacity Performance of Cross-Laminated Timber

Cross-laminated timber (CLT) construction is gaining momentum in the US because it offers multiple advantages over traditional construction methods. Benefits that have received the most attention focus on constructability, the environment, and protection (e.g., blast resistance), although CLT construction is likely to offer other benefits, as well. Still, these have not been studied at length because such evaluations are costly, requiring long-term assessments in an actual building and specialized technical knowledge. Among the possible benefits, CLT construction likely provides a higher-performing building envelope. Using CLT panels to enclose a building means fewer joints in the opaque envelope than what is required in traditional stick-framed construction. Fewer joints mean fewer locations where the air- and water-resistive barrier (WRB) could be compromised; thus, a CLT building enclosure may require less maintenance and have a longer lifespan than a traditionally built structure because of fewer air and water leaks. In addition, CLT’s thermal mass moderates indoor temperatures, allowing the heating, ventilation, and air conditioning (HVAC) system to operate more efficiently during peak hours, reducing operational energy consumption throughout the lifetime of the CLT building (Salonvaara et al., 2022). Furthermore, more stable indoor temperatures can increase occupant comfort. The CLT’s thermal mass can also reduce energy costs by adjusting to utility time-of-use pricing without affecting occupant comfort. The ability of CLT buildings to bridge periods without HVAC operation prepares them for future grid interaction and provides a certain level of resilience against power outages. Researchers have attempted to quantify these benefits; however, their work is based on simplified simulations with numerous assumptions. To correctly understand the benefits, an actual building must be monitored. Therefore, information needs to be gathered on indoor and outdoor temperatures, HVAC energy consumption, thermostat setpoints, temperatures, and thermal transport in CLT components to comprehend how these parameters are affected by the CLT’s thermal mass. These data are needed to reduce the number of assumptions and calibrate simulation models to optimize HVAC controls to minimize overall energy consumption, reduce energy use and higher fees during peak demand, and maintain occupant comfort. Additionally, the calibrated simulation model allows the optimization exercise to be repeated in various US climates. Potential benefits can be tailored to buildings in various locations, and decisions can be made on where CLT construction could be most advantageous. Furthermore, monitoring and simulation results are needed to evaluate the durability of the CLT structures in different climates. This project’s researchers gathered information to help understand and quantify the benefits of CLT buildings concerning operational energy, moderated indoor temperatures, and comfort; the dynamic operation to provide grid services; and resilience in times of power outage. Through the corroboration of simulation models with real-world measurements, this study paves the way for extrapolating findings to other climatic zones and building typologies, thereby broadening the understanding of CLT’s multifaceted benefits and reinforcing its position as a material of choice in sustainable construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗

TPSAS-NF1676L-13606-DND

For over two decades, researchers at NASA Langley Research Center have been combining principles of engineering and psychology to explore the cognitive experiences of system operators. Conceptualizing the application of this interdisciplinary combination in the context of operators? extracurricular activities has resulted in the creation of numerous mental training spin-off technologies. Methods derived from feedback control concepts have been employed to assess the brains responses when managing automation and principles of conditioning have been embodied in devices designed to reinforce effective modes of psychophysiological activity. These methods and devices will be discussed and how they are being applied to ameliorate attention and stress disorders in clinical practice and impediments to optimal performance in sports will be demonstrated.

Alan T Pope↗

Examination of a Methane/Diesel RCCI Engine Using Pele: Preprint

Multi-fuel, advanced injection strategies have become increasingly promising as a strategy to mitigate the emissions generated from internal combustion engines. By carefully controlling the combustion phasing in-cylinder, these new multi-pulse, multi-fuel injection strategies are able to burn in the low-temperature combustion regime where both NOx and soot are not readily produced, reducing the need for extensive exhaust gas recirculation systems. In this study, we examine a reactivity-controlled compression ignition (RCCI) strategy that uses an early pre-filled methane-air mixture with low turbulence background as the low-reactivity fuel and a direct injection of four discrete dodecane jets as a surrogate for the high-reactivity diesel fuel. We use the Pele software suite, a highly optimized, exascale-ready, adaptive mesh refinement codebase to perform high-resolution numerical simulations of a scaled down, single cylinder from the RCCI engine. Here, we resolve the ignition kernels down to micrometer scales and present several statistical quantities evaluating the development of the flow and detailing the onset of ignition and subsequent flame development. Particular attention is paid to the conditions surrounding the onset of the first ignition kernels and discussing what led to the development of those conditions.

CFD↗

Fluid Distribution for In-space Cryogenic Propulsion

The ultimate goal of this task is to enable the use of a single supply of cryogenic propellants for three distinct spacecraft propulsion missions: main propulsion, orbital maneuvering, and attitude control. A fluid distribution system is sought which allows large propellant flows during the first two missions while still allowing control of small propellant flows during attitude control. Existing research has identified the probable benefits of a combined thermal management/power/fluid distribution system based on the Solar Integrated Thermal Management and Power (SITMAP) cycle. Both a numerical model and an experimental model are constructed in order to predict the performance of such an integrated thermal management/propulsion system. This research task provides a numerical model and an experimental apparatus which will simulate an integrated thermal/power/fluid management system based on the SITMAP cycle, and assess its feasibility for various space missions. Various modifications are done to the cycle, such as the addition of a regeneration process that allows heat to be transferred into the working fluid prior to the solar collector, thereby reducing the collector size and weight. Fabri choking analysis was also accounted for. Finally the cycle is to be optimized for various space missions based on a mass based figure of merit, namely the System Mass Ratio (SMR). -. 1 he theoretical and experimental results from these models are be used to develop a design code (JETSIT code) which is able to provide design parameters for such a system, over a range of cooling loads, power generation, and attitude control thrust levels. The performance gains and mass savings will be compared to those of existing spacecraft systems.

Lear, William↗

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗