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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Comparison of ice-phase microphysical parameterization schemes using numerical simulations of tropical convection

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

Mccumber, Michael↗
NASA NTRS2020

Multibody Simulation of REMORA CubeSat Docking to and Pushing a Spent Rocket Booster

This paper details the multibody simulation of three phases of critical importance to the feasibility of the REMORA CubeSat space debris mitigation concept: the final approach of the CubeSat to a spent rocket booster; the grappling of the spent rocket booster using a robotic arm; and the pushing of the spent rocket booster to divert its course from another on-orbit asset. The extension of a robotic mobility and manipulation modeling toolkit (M3TK) from multibody dynamics simulation of manipulators and ground vehicles to simulation of orbital robotics is outlined. This includes the identification of the appropriate parameters required to concisely and generically describe thruster loads, thruster mixing, spacecraft control, and spacecraft navigation for the purpose of on-orbit robotics simulation. A high-level spacecraft navigator commands maneuvers to target spacecraft states. A PID spacecraft controller takes the target states and calculates desired forces and torques. A thruster mixer solves a quadratic program to determine the optimal thruster firing times for the propulsion system. Pulse width modulated actuation of eight canted cold gas thrusters is used in the simulated approach to a rocket nozzle from a distance of 200 m. A five degree of freedom robotic arm is controlled to position a pair of pincers to grasp the rocket nozzle. Contact dynamics are used to accurately simulate the grasping of the rocket nozzle by the pincers. A similar, but separate simulation is performed to assess the ability of the REMORA CubeSat to push the large spent rocket booster. This diversion maneuver makes use of an additional, larger thruster, and pushes the rocket booster in excess of 400 m. Appropriate motor control gains on the robotic arm are found to be higher during the pushing phase than those which are appropriate during free motion; this increase promotes rigidity of the arm and allows it to properly direct the pushing force. Challenges encountered in time step selection for numerical stability of the simulation are also discussed.

Mukherjee, Rudranarayan↗
DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2020

Shape Optimization of Vortex Generators to Control Mack Mode Amplification

This paper demonstrates the potential to use shape optimization for the design of vortex generators in an axisymmetric boundary layer. This shape optimization increases the amplitude of stationary streaks created by streamwise vorticity, with the goal of reducing the amplification of Mack mode instabilities. The test case under consideration matches a trajectory point during the ascent phase of the HIFiRE-1 flight experiment. Wall-mounted vortex generators are added to the cone and their shape is optimized to control the amplification of Mack mode instabilities that are known to initiate laminar-turbulent transition in this flow. An empirical objective function is developed based on the previous studies of optimal streaks. A constraint on the maximum streak amplitude is added to avoid the initiation of bypass transition. The shape optimization is conducted using SU2, an open-source suite for multiphysics simulation and design. A significant improvement is observed in an integral metric of the streak amplitude, while the maximum streak amplitude is maintained close to the baseline level. The qualitative features of the optimized geometry are discussed.

Vortex generators↗
DOE OSTI2025

Reducing NOx Emissions in Ammonia Combustors

Ammonia continues to attract growing interest as a carbonneutral replacement fuel, motivating numerous research efforts toward understanding fundamental ammonia combustion characteristics. A major challenge for the use of ammonia is the development of combustor technologies for mitigating potentially high NOx emissions from the fuel-bound nitrogen chemical pathways to acceptable levels. Our work focuses on a staged RQL combustor architecture for minimizing the NOx emission levels through burning fuel-rich in the primary stage to formcombustion products containing significant levels of hydrogen in addition to nitrogen and water with minimal NOx formation. The subsequent quench and burnout stages of the combustor must then quickly burn residual hydrogen with flame-temperatures moderated by nitrogen and water forming in the first stage. Chemical Reactor Network modeling was used to understand and identify optimal stoichiometry and residence times in each stage for minimizing NOx emissions and to quantify pressure and temperature effects. Reducing the overall NOx emissions requires relatively long residence times in the primary stage to achieve near equilibrium NO levels due to kinetically controlling processes. For conditions relevant to gas turbines (e.g., 30 atm), our work indicates that NOx emissions below 20 ppm are theoretically achievable in a staged RQL combustor architecture. However, these emission predictions significantly depend on the accuracies of currently available chemical kinetic mechanisms which have not been extensively validated under elevated pressure and temperature conditions relevant to gas turbines.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Cellular fluidics

The natural world provides many examples of multiphase transport and reaction processes that have been optimized by evolution. These phenomena take place at multiple length and time scales and typically include gas–liquid–solid interfaces and capillary phenomena in porous media. Many biological and living systems have evolved to optimize fluidic transport. However, living things are exceptionally complex and very difficult to replicate, and human-made microfluidic devices (which are typically planar and enclosed) are highly limited for multiphase process engineering. In this paper, we introduce the concept of cellular fluidics: a platform of unit-cell-based, three-dimensional structures—enabled by emerging 3D printing methods—for the deterministic control of multiphase flow, transport and reaction processes. We show that flow in these structures can be ‘programmed’ through architected design of cell type, size and relative density. We demonstrate gas–liquid transport processes such as transpiration and absorption, using evaporative cooling and CO 2 capture as examples. We design and demonstrate preferential liquid and gas transport pathways in three-dimensional cellular fluidic devices with capillary-driven and actively pumped liquid flow, and present examples of selective metallization of pre-programmed patterns. Our results show that the design and fabrication of architected cellular materials, coupled with analytical and numerical predictions of steady-state and dynamic behaviour of multiphase interfaces, provide deterministic control of fluidic transport in three dimensions. Cellular fluidics may transform the design space for spatial and temporal control of multiphase transport and reaction processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Parametric Study of Panel PCM–Air Heat Exchanger Designs

Heat exchangers, devices for the transfer of heat between two or more working fluids, are extensively used in cooling applications and heating applications. Heat exchangers in buildings are typically components of space-conditioning systems, as well as of water-heating applications. Heat exchangers are also sometimes used in applications that require storage and release of energy at specific times. Phase change materials (PCMs) enhance these heat-exchange processes, given their ability to melt and solidify at a fixed range of temperatures, absorbing or releasing significant amounts of latent heat. Five different configurations of PCM–air heat exchangers for thermal control in buildings are analyzed in this work. The heat exchangers were fitted with PCM encapsulated in plastic and composite pouches of various shapes, and packaged in stackable panel layers. Three-dimensional computational fluid dynamics (CFD) modeling of coupled incompressible fluid and conjugate heat transfer were performed on the designs. The phase change process was numerically modelled using the apparent heat capacity method. Steady-state CFD simulations provided quantification of pressure drop as a function of air flow velocity. Transient simulation results describe the thermal evolution of PCM in the pouches, helping to determine the best performing configuration with respect to total thermal charging time.

10 SYNTHETIC FUELS↗
DOE OSTI2026

Optimal Operation of Residential High Performance Water Heater for Reduction of Electricity Cost and Peak Demand Through Field Validation

Water heating accounts for about 18% of a typical US home’s energy use. Modern water heaters have enabled control options through APIs, offering customers the opportunity to reduce their energy cost and peak demand by dynamically adjusting settings. A water heater’s capacity to store energy using its storage tank makes it an asset for peak demand reduction and energy cost savings. For this reason, a mixed-integer linear programming model is proposed to minimize the energy cost of a high-performance water heater while also reducing the peak demand of the residential household under a time-of-use utility rate by dynamically changing the water heater’s running mode. Specifically, a multi-objective optimization model is formulated to determine the mode settings of the water heater considering hot water use, time-of-use rate, and peak demand limit of the residential household. The mode settings are associated with different dead bands of water temperature for triggering on/off action of the heat pump and heating element. A 66-gal hybrid electric high performance water heater was used for numerical simulation and practical experiments. The simulation results were well aligned with measurements of practical experiments, validating the soundness of the thermodynamic model. In addition, reductions of energy cost, enabling affordability, and reducing peak demand are demonstrated. The research team also developed a software framework with dashboards to automatically and continuously monitor and manage devices.

Liu, Guodong [ORNL] (ORCID:0000000213498608)↗
DOE OSTI2024

Bridging length scales in hard materials with ultra-small angle X-ray scattering – a critical review

Owing to their exceptional properties, hard materials such as advanced ceramics, metals and composites have enormous economic and societal value, with applications across numerous industries. Understanding their microstructural characteristics is crucial for enhancing their performance, materials development and unleashing their potential for future innovative applications. However, their microstructures are unambiguously hierarchical and typically span several length scales, from sub-ångstrom to micrometres, posing demanding challenges for their characterization, especially for in situ characterization which is critical to understanding the kinetic processes controlling microstructure formation. This review provides a comprehensive description of the rapidly developing technique of ultra-small angle X-ray scattering (USAXS), a nondestructive method for probing the nano-to-micrometre scale features of hard materials. USAXS and its complementary techniques, when developed for and applied to hard materials, offer valuable insights into their porosity, grain size, phase composition and inhomogeneities. We discuss the fundamental principles, instrumentation, advantages, challenges and global status of USAXS for hard materials. Using selected examples, we demonstrate the potential of this technique for unveiling the microstructural characteristics of hard materials and its relevance to advanced materials development and manufacturing process optimization. We also provide our perspective on the opportunities and challenges for the continued development of USAXS, including multimodal characterization, coherent scattering, time-resolved studies, machine learning and autonomous experiments. Our goal is to stimulate further implementation and exploration of USAXS techniques and inspire their broader adoption across various domains of hard materials science, thereby driving the field toward discoveries and further developments.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Multi-Parameter Aerosol Scattering Sensor

This work relates to the development of sensors that measure specific aerosol properties. These properties are in the form of integrated moment distributions, i.e., total surface area, total mass, etc., or mathematical combinations of these moment distributions. Specifically, the innovation involves two fundamental features: a computational tool to design and optimize such sensors and the embodiment of these sensors in actual practice. The measurement of aerosol properties is a problem of general interest. Applications include, but are not limited to, environmental monitoring, assessment of human respiratory health, fire detection, emission characterization and control, and pollutant monitoring. The objectives for sensor development include increased accuracy and/or dynamic range, the inclusion in a single sensor of the ability to measure multiple aerosol properties, and developing an overall physical package that is rugged, compact, and low in power consumption, so as to enable deployment in harsh or confined field applications, and as distributed sensor networks. Existing instruments for this purpose include scattering photometers, direct-reading mass instruments, Beta absorption devices, differential mobility analyzers, and gravitational samplers. The family of sensors reported here is predicated on the interaction of light and matter; specifically, the scattering of light from distributions of aerosol particles. The particular arrangement of the sensor, e.g. the wavelength(s) of incident radiation, the number and location of optical detectors, etc., can be derived so as to optimize the sensor response to aerosol properties of practical interest. A key feature of the design is the potential embodiment as an extremely compact, integrated microsensor package. This is of fundamental importance, as it enables numerous previously inaccessible applications. The embodiment of these sensors is inherently low maintenance and high reliability by design. The novel and unique features include the underlying computational underpinning that allows the optimization for specific applications, and the physical embodiment that affords the construction of a compact, durable, and reliable integrated package. The advantage appears in the form of increased accuracy relative to existing instruments, and the applications enabled by the physical attributes of the resulting configuration

Greenberg, Paul S.↗
NASA NTRS2019

Workload Considerations in Urban Air Mobility

Urban air mobility (UAM) is receiving increased attention in aviation as a system for passenger and cargo-carrying new entrants in urban airspace. In order to develop a safe and efficient system, numerous possible concepts of operation for UAM are being explored throughout industry and research domains, the features and assumptions of which may differ according to near, medium and far term operations. Much of the current research into the development of UAM has dominantly focused on technological and engineering capabilities, such as vehicle development. Although these areas of research are essential to furthering UAM, research into the role of the human operator in UAM is limited. The research described in this paper aims to begin to address this gap by investigating the capabilities and implications of human operators as traffic managers in the UAM system, focusing on near-term UAM operations. A human in the loop air traffic control simulation was used to investigate the effect of UAM traffic density, airspace routes and communication procedures on subjective workload and efficiency-related task performance. Findings indicate that medium and high-density operations were associated with high workload. A reduction in verbal communications through a letter of agreement, and optimized routes, were associated with reduced workload and increased performance efficiency. However, even with these adjustments, reported workload remained high, particularly during the high-density scenario. Future research should focus on the human operator roles and responsibilities, and the amount of involvement, in UAM system management. Particular focus should be directed on the impact of reduced human operator involvement and increased automation, on the safety and efficiency of UAM operations and the integration of UAM with traditional air traffic management.

Edwards, Tamsyn↗
DOE OSTI2025

Laser-initiated p – 11 B fusion reactions in petawatt high-repetition-rate laser facilities

Driving of the nuclear fusion reaction p + 11 B → 3α + 8.7 MeV under laboratory conditions by interaction between high-power laser pulses and matter has become a popular field of research, owing to its numerous potential applications: as an alternative to deuterium–tritium for fusion energy production, astrophysics studies, and alpha-particle generation for medical treatment. One possible scheme for laser-driven p– 11 B reactions is to direct a beam of laser-accelerated protons onto a boron (B) sample (the so-called “pitcher-catcher” scheme). This technique has been successfully implemented on large high-energy lasers, yielding hundreds of joules per shot at low repetition. We present here a complementary approach, exploiting the high repetition rate of the VEGA III petawatt laser at CLPU (Spain), aiming at accumulating results from many interactions at much lower energy, to provide better control of the parameters and the statistics of the measurements. Despite a moderate energy per pulse, our experiment allowed exploration of the laser-driven fusion process with tens (up to hundreds) of laser shots. The experiment provided a clear signature of the reactions involved and of the fusion products, accumulated over many shots, leading to an improved optimization of the diagnostics for experimental campaigns of this type. In this paper, we discuss the effectiveness of laser-driven p– 11 B fusion in the pitcher–catcher scheme, at a high repetition rate, addressing the challenges of this experimental scheme and highlighting its critical aspects. Our proposed methodology allows evaluation of the performance of this scheme for laser-driven alpha particle production and can be adapted to high-repetition-rate laser facilities with higher energy and intensity.

Alpha particles↗
DOE OSTI2022

Resilient Operating Constraints for Power Distribution Systems under Setpoint Attacks

Integration and operation of distributed generation (DG) and energy storage (ES) in power distribution systems are enabled by communication networks and embedded sensor and control devices that increase the vulnerability of the systems to cyber-threats, broadening the attack surface and making adversary actions more unpredictable. This paper proposes a methodology that uses ellipsoidal approximations to quantify the potential damage caused by successful attacks that affect, directly or indirectly, the desired operation setpoints and may drive the power distribution operation to unsafe states by violating the limits of voltage or line flows. More specifically, a new methodology is introduced to find the optimal non-symmetric operating constraints that can be imposed to each DG and ES in order to guarantee that the power distribution system is resilient to any malicious setpoints. The proposed method takes as inputs the system topology, DG and ES capabilities, and load limits to solve a convex optimization problem formulated using linear matrix inequalities (LMIs) and the power flow equations. The proposed solution is agnostic to the attacker's action or load profile and it does not require any assumption about the location or means of the attack. The numerical results on a test distribution feeder with several DG and ES illustrate how the proposed resilient operating constraints guarantee the security of the power distribution system under setpoint attacks.

Giraldo, Jairo↗
DOE OSTI2024

Tuning Surface Adhesion Using Grayscale Electron-beam Lithography

Surface texturing of manufactured products tailors their properties, such as friction, adhesion, biocompatibility, or fluid interactions. However, advancements in this area are largely the result of trial-and-effort testing and generally lack a science-guided framework for determining the surface topography that will optimize performance. The present investigation explores grayscale electron-beam lithography as a means to create multiscale surface patterns to control surface performance. Here, we created and characterized a set of surface textures on a silicon wafer; the textures were superpositions of sine waves of varying wavelengths and amplitudes. First, the multiscale topography of the patterned surface was characterized, using profilometry and atomic force microscopy, to understand its fidelity to the designed-in pattern. The results of this analysis demonstrated how grayscale lithography accurately controlled the lateral size of features but was less precise on the vertical height of the surface, and also introduced inherent roughness below the scale of patterning. Second, a micromechanical tester was used to characterize the adhesion of the surfaces with large-scale polished silicon spheres. The results showed that adhesion could be tailored, with significant contribution from all of the designed-in length scales of topography. The strength of adhesion did not correlate with conventional roughness parameters but could be accurately modeled using simple numerical integration. Taken together, this investigation demonstrates the promise and challenges of grayscale e-beam lithography with multiscale patterns as a method for the tailoring of surface performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Combined design of structures and controllers for optimal maneuverability

Approaches to the combined design of structures and controllers for achieving optimal maneuverability are presented. A maneuverability index which directly reflects the minimum time required to perform a given set of maneuvers is introduced. By designing the flexible appendages, the maneuver time of the spacecraft is minimized under the constraints of structural properties, and post maneuver spillover is kept within a specified bound. The spillover reduction is achieved by making use of an appropriate reduced order model. The distributed parameter design problem is approached using assumed shape functions, and finite element analysis with dynamic reduction. Solution procedures have been investigated. Approximate design methods have been developed to overcome the computational difficulties. Some new constraints on the modal frequencies of the spacecraft are introduced in the original optimization problem to facilitate the solution process. It is shown that the global optimal design may be obtained by tuning the natural frequencies to satisfy specific constraints. Researchers quantify the difference between a lower bound to the solution for maneuver time associated with the original problem and the estimate obtained from the modified problem, for a specified application requirement. Numerical examples are presented to demonstrate the capability of this approach.

Ling, Jer↗
NASA NTRS2009

Advanced Computational and Experimental Techniques for Nacelle Liner Performance Evaluation

The Curved Duct Test Rig (CDTR) has been developed to investigate sound propagation through a duct of size comparable to the aft bypass duct of typical aircraft engines. The axial dimension of the bypass duct is often curved and this geometric characteristic is captured in the CDTR. The semiannular bypass duct is simulated by a rectangular test section in which the height corresponds to the circumferential dimension and the width corresponds to the radial dimension. The liner samples are perforate over honeycomb core and are installed on the side walls of the test section. The top and bottom surfaces of the test section are acoustically rigid to simulate a hard wall bifurcation or pylon. A unique feature of the CDTR is the control system that generates sound incident on the liner test section in specific modes. Uniform air flow, at ambient temperature and flow speed Mach 0.275, is introduced through the duct. Experiments to investigate configuration effects such as curvature along the flow path on the acoustic performance of a sample liner are performed in the CDTR and reported in this paper. Combinations of treated and acoustically rigid side walls are investigated. The scattering of modes of the incident wave, both by the curvature and by the asymmetry of wall treatment, is demonstrated in the experimental results. The effect that mode scattering has on total acoustic effectiveness of the liner treatment is also shown. Comparisons of measured liner attenuation with numerical results predicted by an analytic model based on the parabolic approximation to the convected Helmholtz equation are reported. The spectra of attenuation produced by the analytic model are similar to experimental results for both walls treated, straight and curved flow path, with plane wave and higher order modes incident. The numerical model is used to define the optimized resistance and reactance of a liner that significantly improves liner attenuation in the frequency range 1900-2400 Hz. A liner impedance descriptor is used to determine the liner parameters that achieve the optimum impedance.

Gerhold, Carl H.↗