Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “physical 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multi-level optimization with the koopman operator for data-driven, domain-aware, and dynamic system security

Cyber-Physical Systems (CPSs) like the power grid are critically important but also increasingly vulnerable; ensuring reliable system operation in the face of disruptions is becoming more and more challenging. Multi-Level Optimization (MLO) is a powerful way to model adversarial interactions, which naturally makes it applicable to studying CPS security. However, MLO typically does not address underlying system dynamics, and incorporating nonlinear dynamics is generally infeasible. In this paper, we show how to combine MLO with the Koopman Operator (KO) to remedy this. The KO maps nonlinear dynamics to a lifted space in which those dynamics are linear, thus making it ideal for use with MLO. Moreover, the structure of the KO also provides convenient ways to incorporate domain knowledge into the data-driven process of learning the KO representation of a given system. Here we then demonstrate the use of MLO-KO on a small example problem taken from the power grid domain, discuss the scalability and computational cost of MLO-KO, and identify future research directions for this work.

42 ENGINEERING↗
DOE OSTI2025

Constraining nuclear mass models using 𝑟-process observables with multiobjective optimization

Modeling nuclear masses, particularly for nuclei far from stability, remains a key objective in nuclear physics. One contemporary approach is machine learning (ML), which trains on experimental data, but can suffer large errors when extrapolating toward neutron-rich species. In nature, such masses shape observables for the rapid neutron capture process (𝑟 process), which in principle could inform ML models. Here, we introduce a multiobjective optimization approach using the Pareto front algorithm. We show that this technique, capable of identifying models that generate 𝑟-process abundances aligning with both solar and stellar data, is a promising method to select ML models with reliable extrapolation power.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Heat Exchanger Design for Waste Heat Recovery using Additive Manufacturing

With anticipated increases in global energy use and their associated emissions, there are many efforts to improve the emissions and efficiency of many industrial processes through waste heat recovery (WHR). To maximize improvements, the heat exchangers in WHR power cycles must be designed and optimized for a given heat source. Meanwhile, additive manufacturing (AM) is a single part manufacturing technique that provides the opportunity to develop unique and complex heat exchanger designs. Consequently, this work centers on the modeling, design, optimization, and testing of AM shell and tube heat exchangers for WHR. Specifically, a modified one-dimensional (1+1-D) and three-dimensional (3-D) heat exchanger models were developed and validated, with the 1+1-D model serving as the basis for design optimization. After modifications to the design to allow for AM, heat exchanger geometry optimizations were completed, targeting increased effectiveness across a range of physical parameters. Of the designs considered, nine were additively manufactured using Ti5553, and then tested using an external heat source and working fluid. Experimentally, the nineteen-tube, zero-baffle heat exchanger displayed the best overall performance, with the effects of dimensional accuracy, surface roughness, and baffles in the other designs also indicated. Next, the experimental results were used as inputs for the 1+1-D model. Although the original model was able to predict fluid outlet temperatures within 5% differences, larger variations in the heat exchanger performance parameters, temperature differences, and pressure drop necessitated optimizing the pressure drop and heat transfer calculations within the model. Thus, the resulting improved model was able to predict the heat transfer rates, heat source temperature difference, and effectiveness within 15% differences, while showing improvements in the pressure drops and working fluid temperature differences. The improved model serves as the foundation for future work to improve the design and performance of AM heat exchangers.

Spickler, Bailey Ann↗
DOE OSTI2020

Automated Adversary Emulation for Cyber-Physical Systems via Reinforcement Learning

Adversary emulation is an offensive exercise that provides a comprehensive assessment of a system’s resilience against cyber attacks. However, adversary emulation is typically a manual process, making it costly and hard to deploy in cyber-physical systems (CPS) with complex dynamics, vulnerabilities, and operational uncertainties. In this paper, we develop an automated, domain-aware approach to adversary emulation for CPS. We formulate a Markov Decision Process (MDP) model to determine an optimal attack sequence over a hybrid attack graph with cyber (discrete) and physical (continuous) components and related physical dynamics. We apply model-based and model-free reinforcement learning (RL) methods to solve the discrete-continuous MDP in a tractable fashion. As a baseline, we also develop a greedy attack algorithm and compare it with the RL procedures. We summarize our findings through a numerical study on sensor deception attacks in buildings to compare the performance and solution quality of the proposed algorithms.

Bhattacharya, Arnab↗
DOE OSTI2021

Multiobjective Dynamic Optimization for Optimal Load-Following of Natural Gas Combined Cycle Power Plants under Stress Constraints

With the increasing penetration of intermittent renewable energy sources into the electric grid, there is an associated need for conventional thermal power plants that are designed to operate at base-load conditions to cycle their load rapidly and frequently and operate under low-load conditions. Fast load-following operations lead to plant efficiency loss and the reduction of equipment life. A high-fidelity dynamic model of a natural gas combined cycle power plant, with rigorous equipment level submodels, is developed to capture the plant transient performance and the thermomechanical stress evolution at the high-pressure drum is computed to assess the drum life consumption. Stress evolution is modeled for the location of the edge at the drum/downcomer junction that experiences high circumferential stress amplitude during the fast load-following operation, leading to higher fatigue damage than locations considered under existing design standards. A dynamic optimization problem is solved for optimal load-following operation. Depending on the value of the stress constraint and the desired ramp rate, satisfying the desired stress constraint may be infeasible without relaxing the ramp rate. A multiobjective dynamic optimization problem is solved using a lexicographic approach that minimizes ramp rate relaxation and maximizes efficiency while satisfying stress constraints, avoiding spraying to saturation, and maintaining main steam and reheat steam temperature within bounds. It was observed that the optimal ramp rate can be highly nonlinear as opposed to the industry-standard linear ramp rate. Nonlinear optimal ramp-rate profiles not only help to avoid stress constraints that may be unavoidable by using the linear profile but also result in higher efficiency than the linear profile. Furthermore, the study shows that there is a strong tradeoff between the relaxation in the plant ramp rate and the time-average thermal efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Computationally efficient control allocation

A computationally efficient method for calculating near-optimal solutions to the three-objective, linear control allocation problem is disclosed. The control allocation problem is that of distributing the effort of redundant control effectors to achieve some desired set of objectives. The problem is deemed linear if control effectiveness is affine with respect to the individual control effectors. The optimal solution is that which exploits the collective maximum capability of the effectors within their individual physical limits. Computational efficiency is measured by the number of floating-point operations required for solution. The method presented returned optimal solutions in more than 90% of the cases examined; non-optimal solutions returned by the method were typically much less than 1% different from optimal and the errors tended to become smaller than 0.01% as the number of controls was increased. The magnitude of the errors returned by the present method was much smaller than those that resulted from either pseudo inverse or cascaded generalized inverse solutions. The computational complexity of the method presented varied linearly with increasing numbers of controls; the number of required floating point operations increased from 5.5 i, to seven times faster than did the minimum-norm solution (the pseudoinverse), and at about the same rate as did the cascaded generalized inverse solution. The computational requirements of the method presented were much better than that of previously described facet-searching methods which increase in proportion to the square of the number of controls.

Durham, Wayne↗
DOE OSTI2022

Numerical gate synthesis for quantum heuristics on bosonic quantum processors

There is a recent surge of interest and insights regarding the interplay of quantum optimal control and variational quantum algorithms. We study the framework in the context of qudits which are, for instance, definable as controllable electromagnetic modes of a superconducting cavity system coupled to a transmon. By employing recent quantum optimal control approaches described in (Petersson and Garcia, 2021), we showcase control of single-qudit operations up to eight states, and two-qutrit operations, mapped respectively onto a single mode and two modes of the resonator. We discuss the results of numerical pulse engineering on the closed system for parametrized gates useful to implement Quantum Approximate Optimization Algorithm (QAOA) for qudits. The results show that high fidelity ( > 0.99) is achievable with sufficient computational effort for most cases under study, and extensions to multiple modes and open, noisy systems are possible. The tailored pulses can be stored and used as calibrated primitives for a future compiler in circuit quantum electrodynamics (cQED) systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

Active Mirror Predictive and Requirements Verification Software (AMP-ReVS)

This software is designed to predict large active mirror performance at various stages in the fabrication lifecycle of the mirror. It was developed for 1-meter class powered mirrors for astronomical purposes, but is extensible to other geometries. The package accepts finite element model (FEM) inputs and laboratory measured data for large optical-quality mirrors with active figure control. It computes phenomenological contributions to the surface figure error using several built-in optimization techniques. These phenomena include stresses induced in the mirror by the manufacturing process and the support structure, the test procedure, high spatial frequency errors introduced by the polishing process, and other process-dependent deleterious effects due to light-weighting of the mirror. Then, depending on the maturity of the mirror, it either predicts the best surface figure error that the mirror will attain, or it verifies that the requirements for the error sources have been met once the best surface figure error has been measured. The unique feature of this software is that it ties together physical phenomenology with wavefront sensing and control techniques and various optimization methods including convex optimization, Kalman filtering, and quadratic programming to both generate predictive models and to do requirements verification. This software combines three distinct disciplines: wavefront control, predictive models based on FEM, and requirements verification using measured data in a robust, reusable code that is applicable to any large optics for ground and space telescopes. The software also includes state-of-the-art wavefront control algorithms that allow closed-loop performance to be computed. It allows for quantitative trade studies to be performed for optical systems engineering, including computing the best surface figure error under various testing and operating conditions. After the mirror manufacturing process and testing have been completed, the software package can be used to verify that the underlying requirements have been met.

Basinger, Scott A.↗
DOE OSTI2022

Non-trivial symmetries in quantum landscapes and their resilience to quantum noise

Very little is known about the cost landscape for parametrized Quantum Circuits (PQCs). Nevertheless, PQCs are employed in Quantum Neural Networks and Variational Quantum Algorithms, which may allow for near-term quantum advantage. Such applications require good optimizers to train PQCs. Recent works have focused on quantum-aware optimizers specifically tailored for PQCs. However, ignorance of the cost landscape could hinder progress towards such optimizers. In this work, we analytically prove two results for PQCs: (1) We find an exponentially large symmetry in PQCs, yielding an exponentially large degeneracy of the minima in the cost landscape. Alternatively, this can be cast as an exponential reduction in the volume of relevant hyperparameter space. (2) We study the resilience of the symmetries under noise, and show that while it is conserved under unital noise, non-unital channels can break these symmetries and lift the degeneracy of minima, leading to multiple new local minima. Based on these results, we introduce an optimization method called Symmetry-based Minima Hopping (SYMH), which exploits the underlying symmetries in PQCs. Our numerical simulations show that SYMH improves the overall optimizer performance in the presence of non-unital noise at a level comparable to current hardware. Overall, this work derives large-scale circuit symmetries from local gate transformations, and uses them to construct a noise-aware optimization method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

ARM-IRL: Adaptive Resilience Metric Quantification Using Inverse Reinforcement Learning

The resilience of safety-critical systems is gaining importance due to the rise in cyber and physical threats, especially within critical infrastructure. Traditional static resilience metrics may not capture dynamic system states, leading to inaccurate assessments and ineffective responses to cyber threats. This work aims to develop a data-driven, adaptive method for resilience metric learning. We propose a data-driven approach using inverse reinforcement learning (IRL) to learn a single, adaptive resilience metric. The method infers a reward function from expert control actions. Unlike previous approaches using static weights or fuzzy logic, this work applies adversarial inverse reinforcement learning (AIRL), training a generator and discriminator in parallel to learn the reward structure and derive an optimal policy. The proposed approach is evaluated on multiple scenarios: optimal communication network rerouting, power distribution network reconfiguration, and cyber–physical restoration of critical loads using the IEEE 123-bus system. The adaptive, learned resilience metric enables faster critical load restoration in comparison to conventional RL approaches.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Implementing Scientific Simulation Codes Highly Tailored for Vector Architectures Using Custom Configurable Computing Machines

The motivation for this work comes from an observation that amidst the push for Massively Parallel (MP) solutions to high-end computing problems such as numerical physical simulations, large amounts of legacy code exist that are highly optimized for vector supercomputers. Because re-hosting legacy code often requires a complete re-write of the original code, which can be a very long and expensive effort, this work examines the potential to exploit reconfigurable computing machines in place of a vector supercomputer to implement an essentially unmodified legacy source code. Custom and reconfigurable computing resources could be used to emulate an original application's target platform to the extent required to achieve high performance. To arrive at an architecture that delivers the desired performance subject to limited resources involves solving a multi-variable optimization problem with constraints. Prior research in the area of reconfigurable computing has demonstrated that designing an optimum hardware implementation of a given application under hardware resource constraints is an NP-complete problem. The premise of the approach is that the general issue of applying reconfigurable computing resources to the implementation of an application, maximizing the performance of the computation subject to physical resource constraints, can be made a tractable problem by assuming a computational paradigm, such as vector processing. This research contributes a formulation of the problem and a methodology to design a reconfigurable vector processing implementation of a given application that satisfies a performance metric. A generic, parametric, architectural framework for vector processing implemented in reconfigurable logic is developed as a target for a scheduling/mapping algorithm that maps an input computation to a given instance of the architecture. This algorithm is integrated with an optimization framework to arrive at a specification of the architecture parameters that attempts to minimize execution time, while staying within resource constraints. The flexibility of using a custom reconfigurable implementation is exploited in a unique manner to leverage the lessons learned in vector supercomputer development. The vector processing framework is tailored to the application, with variable parameters that are fixed in traditional vector processing. Benchmark data that demonstrates the functionality and utility of the approach is presented. The benchmark data includes an identified bottleneck in a real case study example vector code, the NASA Langley Terminal Area Simulation System (TASS) application.

Rutishauser, David↗
NASA NTRS2002

Optimal Design of Functionally Graded Metallic Foam Insulations

The focus of our work has been on developing an insight into the physics that govern the optimum design of thermal insulation for use in thermal protection systems of launch vehicle. Of particular interest was to obtain optimality criteria for designing foam insulations that have density (or porosity) distributions through the thickness for optimum thermal performance. We investigate the optimum design of functionally graded thermal insulation for steady state heat transfer through the foam. We showed that the heat transfer in the foam has competing modes, of radiation and conduction. The problem assumed a fixed inside temperature of 400 K and varied the aerodynamic surface heating on the outside surface from 0.2 to 1.0 MW/sq m. The thermal insulation develops a high temperature gradient through the thickness. Investigation of the model developed for heat conduction in foams showed that at high temperatures (as on outside wall) intracellular radiation dominates the heat transfer in the foam. Minimizing radiation requires reducing the pore size, which increases the density of the foam. At low temperatures (as on the inside wall), intracellular conduction (of the metal and air) dominates the heat transfer. Minimizing conduction requires increasing the pore size. This indicated that for every temperature there was an optimum value of density that minimized the heat transfer coefficient. Two optimization studies were performed. One was to minimize the heat transmitted though a fixed thickness insulation by varying density profiles. The second was to obtain the minimum mass insulation for specified thickness. Analytical optimality criteria were derived for the cases considered. The optimality condition for minimum heat transfer required that at each temperature we find the density that minimizes the heat transfer coefficient. Once a relationship between the optimum heat transfer coefficient and the temperature was found, the design problem reduced to the solution of a simple nonlinear differential equation. Preliminary results of this work were presented at the American Society of Composites meeting, and the final version was submitted for publication in the AIAA Journal. In addition to minimizing the transmitted heat, we investigated the optimum design for minimum weight given an acceptable level of heat transmission through the insulation. The optimality criterion developed was different from that obtained for minimizing beat transfer coefficient. For minimum mass design, we had to find for a given temperature the optimum density, which minimized the logarithmic derivative of the insulation thermal conductivity with respect to its density. The logarithmic derivative is defined as the ratio of relative change in the dependent response (thermal conductivity) to the relative change in the independent variable (density). The results have been documented as a conference paper that will be presented at the upcoming AIAA.

Haftka, Raphael T.↗
DOE OSTI2022

Deep Generative Modeling of Periodic Variable Stars Using Physical Parameters

The ability to generate physically plausible ensembles of variable sources is critical to the optimization of time domain survey cadences and the training of classification models on data sets with few to no labels. Traditional data augmentation techniques expand training sets by reenvisioning observed exemplars, seeking to simulate observations of specific training sources under different (exogenous) conditions. Unlike fully theory-driven models, these approaches do not typically allow principled interpolation nor extrapolation. Moreover, the principal drawback of theory-driven models lies in the prohibitive computational cost of simulating source observables from ab initio parameters. In this work, we propose a computationally tractable machine learning approach to generate realistic light curves of periodic variables capable of integrating physical parameters and variability classes as inputs. Our deep generative model, inspired by the transparent latent space generative adversarial networks, uses a variational autoencoder (VAE) architecture with temporal convolutional network layers, trained using the OGLE-III optical light curves and physical characteristics (e.g., effective temperature and absolute magnitude) from Gaia DR2. A test using the temperature–shape relationship of RR Lyrae demonstrates the efficacy of our generative “physics-enhanced latent space VAE” (PELS-VAE) model. Such deep generative models, serving as nonlinear nonparametric emulators, present a novel tool for astronomers to create synthetic time series over arbitrary cadences.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Evolutionary Design of Controlled Structures

Basic physical concepts of structural delay and transmissibility are provided for simple rod and beam structures. Investigations show the sensitivity of these concepts to differing controlled-structures variables, and to rational system modeling effects. An evolutionary controls/structures design method is developed. The basis of the method is an accurate model formulation for dynamic compensator optimization and Genetic Algorithm based updating of sensor/actuator placement and structural attributes. One and three dimensional examples from the literature are used to validate the method. Frequency domain interpretation of these controlled structure systems provide physical insight as to how the objective is optimized and consequently what is important in the objective. Several disturbance rejection type controls-structures systems are optimized for a stellar interferometer spacecraft application. The interferometric designs include closed loop tracking optics. Designs are generated for differing structural aspect ratios, differing disturbance attributes, and differing sensor selections. Physical limitations in achieving performance are given in terms of average system transfer function gains and system phase loss. A spacecraft-like optical interferometry system is investigated experimentally over several different optimized controlled structures configurations. Configurations represent common and not-so-common approaches to mitigating pathlength errors induced by disturbances of two different spectra. Results show that an optimized controlled structure for low frequency broadband disturbances achieves modest performance gains over a mass equivalent regular structure, while an optimized structure for high frequency narrow band disturbances is four times better in terms of root-mean-square pathlength. These results are predictable given the nature of the physical system and the optimization design variables. Fundamental limits on controlled performance are discussed based on the measured and fit average system transfer function gains and system phase loss.

Masters, Brett P.↗
NASA NTRS2005

Physical and hydraulic properties of baked ceramic aggregates used for plant growth medium

Baked ceramic aggregates (fritted clay, arcillite) have been used for plant research both on the ground and in microgravity. Optimal control of water and air within the root zone in any gravity environment depends on physical and hydraulic properties of the aggregate, which were evaluated for 0.25-1-mm and 1-2-mm particle size distributions. The maximum bulk densities obtained by any packing technique were 0.68 and 0.64 g cm-3 for 0.25-1-mm and 1-2-mm particles, respectively. Wettable porosity obtained by infiltration with water was approximately 65%, substantially lower than total porosity of approximately 74%. Aggregate of both particle sizes exhibited a bimodal pore size distribution consisting of inter-aggregate macropores and intra-aggregate micropores, with the transition from macro- to microporosity beginning at volumetric water content of approximately 36% to 39%. For inter-aggregate water contents that support optimal plant growth there is 45% change in water content that occurs over a relatively small matric suction range of 0-20 cm H2O for 0.25-1-mm and 0 to -10 cm H2O for 1-2-mm aggregate. Hysteresis is substantial between draining and wetting aggregate, which results in as much as a approximately 10% to 20% difference in volumetric water content for a given matric potential. Hydraulic conductivity was approximately an order of magnitude higher for 1-2-mm than for 0.25-1-mm aggregate until significant drainage of the inter-aggregate pore space occurred. The large change in water content for a relatively small change in matric potential suggests that significant differences in water retention may be observed in microgravity as compared to earth.

NASA Center JSC↗
DOE OSTI2022

Calibrating the Classical Hardness of the Quantum Approximate Optimization Algorithm

The trading of fidelity for scale enables approximate classical simulators such as matrix product states (MPSs) to run quantum circuits beyond exact methods. A control parameter, the so-called bond dimension $\mathcal{χ}$ for MPSs, governs the allocated computational resources and the output fidelity. Here, we characterize the fidelity for the quantum approximate optimization algorithm by the expectation value of the cost function that it seeks to minimize and find that it follows a scaling law $\mathscr{F}$(ln $\mathcal{χ}$/N), where N is the number of qubits. With ln $\mathcal{χ}$ amounting to the entanglement that a MPS can encode, we show that the relevant variable for investigating the fidelity is the entanglement per qubit. Importantly, our results calibrate the classical computational power required to achieve the desired fidelity and benchmark the performance of quantum hardware in a realistic setup. For instance, we quantify the hardness of performing better classically than a noisy superconducting quantum processor by readily matching its output to the scaling function. Moreover, we relate the global fidelity to that of individual operations and establish its relationship with $\mathcal{χ}$ and N. We sharpen the requirements for noisy quantum computers to outperform classical techniques at running a quantum optimization algorithm in speed, size, and fidelity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2022

22” ADP Fan Rig Liners Design Report

A liner design study was conducted as part of a cooperative effort between five (5) government/industry teams that together seek to demonstrate the technology for designing and manufacturing acoustic liners that are twenty five percent (25%) more efficient than 1992 technology liners. The study emphasized teaming collaboration. The improved liners were designed by Pratt & Whitney and Boeing Airplane Company and built by Rohr Inc, and will be tested in the NASA/P& W 22-inch ADP fan rig at NASA Lewis Research Center's 9' x 15' wind tunnel. Design guidelines and decisions were made collectively during monthly design review telecons and at the formal final design review. The tools that were used to design the new improved liners were not new, but the process that was developed as part of this study that led into the evaluation and selection of the final and best liner designs is new. Until now, the liner design process as practiced by different industry teams varied greatly. Some procedures were based on empirical liner attenuation databases which could not adequately account for engine-to-engine hardwall fan noise spectral differences, while others were entirely theoretical. And regardless of which method one chose, the major difficulty was still the lack of understanding and knowledge of the actual hardwall fan source noise modal structure and farfield SPL spectra. The NASA-led effort to attempt actual measurements of the fan source noise modal structure by means of the rotating microphone array is expected to contribute significantly to the understanding of the nature of the fan tone noise modes, but the application to broadband noise is still a long way off. The new design process included consideration of the measured fan tone modes, but for the majority of the spectra a separate, systematic process was used. It is a common practice in the engine/nacelle industry that acoustic liners for new products are designed before actual measured hardwall engine far-field noise spectra are available. Target noise spectra are normally derived from existing engine noise databases with some adjustments to absolute SPLs and frequencies. This practice has been acceptable as long as the new engine was a derivative of the base engine. However, in the case of the ADP, the transition from a current engine base is too - great, and the adjustments could not adequately account for the quantum changes in SPL spectral differences. Examples were the 1992 single degree of freedom (SDOF) inlet and aft liners (designed by P&W) that would be used as the baseline liners against which the new improved liners' efficien¬cies would -be measured. As will be shown, these baseline liners have been found to be deeper than the desire optimum depths. The new liner design process begins with the selection of the target hard wall fan noise spectra. These spectra were obtained by scaling up (5.91 scale factor) the measured hardwall fan spectral data from the 22” ADP fan rig. Next, the modal energy contents for each 1/3-octave band center frequency of these hardwall fan noise spectra are estimated. (Within each 1/3-octave band center frequency, the model energy distribution approximation is for both tone and broadband noise). For the inlet noise, P&W uses a derivative of the NASA Lewis (Ed Rice's) method of classifying; and grouping propagating modes by their cutoff ratios. The next step is to assign energy level to each group of modes having the same cutoff ratio. In Ed Rice's model, the modes were grouped according to ten (10) cutoff ratio intervals with center values located at 1.026, 1.085, 1.155, 1.24, 1.35, 1.49, 1.69, 2.0, and 4.47 (with equal number of modes in each interval). The modes were assumed to have equal energy. P&W’s model expands the cutoff ratio-mode grouping into two hundred (200) smaller cutoff ratio intervals with center values located f1".'m 1.003 to 11.5 in 199 increasing incremental intervals. Next, P&W's model uses a "2-parameter" normal distribution as a template to assign energy levels to these 200 pre-determined cutoff ratio values (for each 1/3-octave band center frequency). In the past, P&W had conducted an extensive study to determine what "2-parameter" values are appropriate for fullscale inlet liners, and had developed a set of twenty-four (24) "2-pararneter" values (i.e. one for each 1/3-octave frequency band) that when used in Ed Rice's Inlet Attenuation Prediction Method produced predicted liner attenuations that closely matched measured liner attenuations from several P&W's engines. In the absence of actual measured tone and broadband modal data from the 22" ADP rig, the process will use P&W's proprietary set of "2-parameter" values for this liner design study. For the aft noise, Boeing uses a modal energy approximation that the "transport energy" of each propagating mode is equal. This approximation is almost the same as the equal energy per mode approximation, except for modes that are near cutoff. Boeing's model forces these modes to have lower energy levels. Both P&W and Boeing agreed that the "transport energy" approximation should work well for the ADP aft fan noise which appears to be dominated by broadband noise. The design process then proceeds to calculate the optimum liner impedances for each frequency in both the inlet and the aft. These optimum impedances represent the target impedances that the designed liners should have. Liners with impedances matching the optimum impedances at all frequencies are "ideal" liners. These ideal liners are theoretically the best liners. Unfortunately, it has been showed that it is impossible to design and build such ideal liners. The next best liners are ones that have impedances matching the optimum impedances over some frequency range (not all frequencies as for the ideal). This is accomplished by the use of "frequency weightings". Several of P&W's and Boeing's existing computer decks were used for optimizing and matching the designed liner impedances to the target optimum values (with the various frequency weightings specified). The optimization produces liner candidates with predicted liner impedances and descriptions of their physical liner characteristics. These candidate liner impedances are then used to predict their spectral attenuation characteristics which are then used together with the target hardwall fan noise spectra to determine the resulting treated noise spectra and PNLT values. Further optimization around the selected candidate designs yield final designs that are best in PNLT attenuations. Use of the optimization decks allowed a large number of liner candidates to be screened in a relatively short period of time. This design process is systematic and is efficient. The new inlet SDOF liner design obtained from the improved process was predicted to be 34% more effective (per unit area) than the 1992 baseline inlet liner. This inlet liner design is a 112 rayl wovenwiremesh facesheet over a 0.312-inch deep honeycomb core. The new aft SDOF liner design was predicted to be 52% more efficient than the aft baseline liner. The new aft SDOF liner is "segmented" with a shallower liner on the core cowl (inner duct wall) and a deeper liner on the fan cowl (outer duct wall). The shallower liner is a 70.6 rayl woven-wiremesh facesheet over a 0.141-inch deep honeycomb core. The deeper liner is a 68 rayl woven-wiremesh facesheet over a 0.309-inch deep honeycomb core. All liner dimensions are for the 22-inch model-scale ADP fan rig liners. The selected advanced liners are "segmented" double-layer (DDOF) liners for the inlet and aft locations. Also, for the inlet, a bulk liner with ceramic foam for wider broadband noise absorption was also selected. Triple-layer liner designs were not considered since the model scaled liners ( 1/5. 91) were dimensionally too small to be built correctly, and irrin earlier concept study, Boeing found a triple-layer to have only very small benefits over a double-layer. The inlet DDOF liner was predicted to be 83% more effective than the baseline. The inlet DDOF design is a 78 rayl facesheet over a 0.080 top cavity depth, a 68 rayl septum and a 0.227-inch bottom cavity depth. The inlet bulk liner is a 60 rayl facesheet over a 0.33-inch deep honeycomb filled with high temperature (HTP) ceramic foam with a density of 4.8 lb/cu.ft and a flow resistivity of 167 rayl/cm. The inlet bulk liner was predicted to be 83% more effective than the baseline liner. The aft DDOF liners are segmented with a shallower liner on the core cowl and a deeper liner on the fan cowl. The shallow DDOF liner is a 49.8 rayl facesheet over a 0.093-inch top cavity depth, a septum of 88.2 rayls over a 0.181-inch bottom cavity depth. The deep DDOF liner is a 12.9 rayl facesheet over a 0.140-inch top cavity depth, a septum of 53.1 ray ls over a 0.258-inch botom cavity depth. The segmented aft DDOF liners were predicted to be 86% more efficient than the baseline liner.

turbofan acoustic treatment↗
DOE OSTI2023

Fast model-based scenario optimization in NSTX-U enabled by analytic gradient computation

Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗