Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal sizing method”

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 937 records · Page 52

Television monitor field shifter and an opto-electronic method for obtaining a stereo image of optimal depth resolution and reduced depth distortion on a single screen

A method and apparatus for obtaining a stereo image with reduced depth distortion and optimum depth resolution. Static and dynamic depth distortion and depth resolution trade-off is provided. Cameras obtaining the images for a stereo view are converged at a convergence point behind the object to be presented in the image, and (a) the collection-surface-to-object distance, (b) the camera separation distance, and (c) the focal lengths of zoom lenses for the camera are all increased. Doubling the distances cuts the static depth distortion in half while maintaining image size and depth resolution. Dynamic depth distortion is minimized by panning a stereo view-collecting camera system about a circle which passes through the convergence point and the cameras' first nodal points. Horizontal field shifting of the television fields on a television monitor brings both the monitor and the stereo views within the viewer's limit of binocular fusion.

Daniel B Diner↗

A Globally Convergent Augmented Lagrangian Pattern Search Algorithm for Optimization with General Constraints and Simple Bounds

We give a pattern search adaptation of an augmented Lagrangian method due to Conn, Gould, and Toint. The algorithm proceeds by successive bound constrained minimization of an augmented Lagrangian. In the pattern search adaptation we solve this subproblem approximately using a bound constrained pattern search method. The stopping criterion proposed by Conn, Gould, and Toint for the solution of this subproblem requires explicit knowledge of derivatives. Such information is presumed absent in pattern search methods; however, we show how we can replace this with a stopping criterion based on the pattern size in a way that preserves the convergence properties of the original algorithm. In this way we proceed by successive, inexact, bound constrained minimization without knowing exactly how inexact the minimization is. So far as we know, this is the first provably convergent direct search method for general nonlinear programming.

Lewis, Robert Michael↗

Multicompartment Polymeric Nanoreactors for Non-Orthogonal Cascade Catalysis

Spatial confinement of multiple catalysts presents an effective strategy for performing sequential or tandem chemical transformations in a one-pot reaction. These methods may be used to catalyze numerous reactions in conditions that are otherwise incompatible between catalyst and solvent, different catalysts, or reagents. Appropriate site isolation or support structure design would lead to significant advantages in atom economy, purification, and costs and the development of the interface between a catalyst and its confined microenvironment is paramount for realizing the next-generation of nanoreactors. Polymer scaffolds can create tailor-made microenvironments resulting in catalyst compartmentalization. Through the optimization of a number of variables such as size, solubility, functionality, and morphology of the nanoreactor, catalyst activity and selectivity can be tuned. In this Feature Article, we introduce design principles and early strategies for polymer supports for catalyst site-isolation and discuss current strategies towards multi-compartment polymer nanoreactors for non-orthogonal cascade catalysis. We conclude by outlining future design trends in this burgeoning field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Carbon Ferrochrome Effects on Microstructure and Mechanical Properties of Powder Metallurgy Titanium Alloys

In this paper, the powder metallurgy method was employed to prepare titanium alloys with high-carbon ferrochrome (HCFeCr). With the HCFeCr addition, the size of the β grains decreased and the lamellae of the α phase became thinner. The TiC strengthening phase precipitated at the grain boundaries, resulting in the inhibition of the β grain growth and hardness enhancement. In addition, Fe and Cr as β-stabilizing elements existed in the enriched β-Ti phase, which increased the β phase amount. The strength improvement was attributed to grain size optimization and strengthening phase formation. With the 9 wt.% HCFeCr addition, the corresponding Vickers hardness increased to 480 HV, which is 60% higher than the Ti6Al4V alloy. The tensile and yield strength was 1228 and 1140 MPa, respectively. The addition of HCFeCr can effectively enhance the strength and hardness while reducing the cost. This is a promising additive to obtain low-cost Ti alloys with high mechanical properties.

36 MATERIALS SCIENCE↗

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗

A comprehensive method for preliminary design optimization of axial gas turbine stages. II - Code verification

The present effort represents an extension of previous work wherein a calculation model for performing rapid pitchline optimization of axial gas turbine geometry, including blade profiles, is developed. The model requires no specification of geometric constraints. Output includes aerodynamic performance (adiabatic efficiency), hub-tip flow-path geometry, blade chords, and estimates of blade shape. Presented herein is a verification of the aerodynamic performance portion of the model, whereby detailed turbine test-rig data, including rig geometry, is input to the model to determine whether tested performance can be predicted. An array of seven (7) NASA single-stage axial gas turbine configurations is investigated, ranging in size from 0.6 kg/s to 63.8 kg/s mass flow and in specific work output from 153 J/g to 558 J/g at design (hot) conditions; stage loading factor ranges from 1.15 to 4.66.

Jenkins, R. M.↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Advanced manufacturing of 3D custom boron-carbide collimators designed for complex environments for neutron scattering

Scattered-beam collimation is a very useful method to reduce unwanted backgrounds and to boost the desired sample signal instead. This approach is of particular interest for samples contained within a complex environment that gives rise to much unwanted parasitic scatter. As neutron scattering instrument and techniques advances, small samples are becoming of more and more interest, which necessitates optimized collimation. Here, in this work, we describe a concept for the design and fabrication of advanced scattered-beam collimation 3D printed from B 4 C specifically tailored for samples contained within a complex environment. This concept is demonstrated through the use of a diamond anvil cell for high pressure experimentation, a technique that very typically requires small samples. The collimators here are designed through a modeling procedure via Monte Carlo neutron ray tracing that encompasses the entire system: the instrument, the complex environment and the collimator. Since the first approach of simply scaling up of the print-size was not successful, a novel concept of a multi-part alternate-blade collimator was developed. This approach addresses printing constraints but gives greater flexibility in design. Its performance is computationally compared against an unprintable progressively tighter blade collimator to assess the effect of alternating blades. No strong difference was observed. Its performance was validated through experimentation at the Spallation Neutron Source. The results emphasize the critical importance of ultra-high precision alignment while showing good overall agreement between simulation and experiment and underscore the feasibility of the method and its real-world application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Structural Analysis Using Computer Based Methods

The stiffness of a flex hose that will be used in the umbilical arms of the Space Launch Systems mobile launcher needed to be determined in order to properly qualify ground umbilical plate behavior during vehicle separation post T-0. This data is also necessary to properly size and design the motors used to retract the umbilical arms. Therefore an experiment was created to determine the stiffness of the hose. Before the test apparatus for the experiment could be built, the structure had to be analyzed to ensure it would not fail under given loading conditions. The design model was imported into the analysis software and optimized to decrease runtime while still providing accurate rest~lts and allow for seamless meshing. Areas exceeding the allowable stresses in the structure were located and modified before submitting the design for fabrication. In addition, a mock up of a deep space habitat and the support frame was designed and needed to be analyzed for structural integrity under different loading conditions. The load cases were provided by the customer and were applied to the structure after optimizing the geometry. Once again, weak points in the structure were located and recommended design changes were made to the customer and the process was repeated until the load conditions were met without exceeding the allowable stresses. After the stresses met the required factors of safety the designs were released for fabrication.

Dietz, Matthew R.↗

Power Electronics for a Miniaturized Arcjet

A 0.3 kW Power Processing Unit (PPU) was designed, tested on resistive loads, and then integrated with a miniaturized arcjet. The main goal of the design was to minimize size and mass while maintaining reasonable efficiency. In order to obtain the desired reductions in mass, simple topologies and control methods were considered. The PPU design incorporates a 50 kHz, current-mode-control, pulse-width-modulated (PWM), push-pull topology. An input voltage of 28 +/- 4V was chosen for compatibility with typical unregulated low voltage busses anticipated for smallsats. An efficiency of 0.90 under nominal operating conditions was obtained. The component mass of the PPU was 0.475 kg and could be improved by optimization of the output filter design. The estimated mass for a flight PPU based on this design is less than a kilogram.

Pinero, Luis R.↗

Optimization of a Solver for Computational Materials and Structures Problems on NVIDIA Volta and AMD Instinct GPUs

The Scalable Implementation of Finite Elements by NASA (ScIFEN) is a software package developed to solve complex computational materials and structures problems using the finite element method (FEM). In this paper, we describe optimization techniques to speed up the linear solver computation that occurs within the ScIFEN application. We consider GPUs from two different vendors, NVIDIA and AMD as our target platforms for optimization and highlight differences in performance and optimization techniques. The NVIDIA GPU Volta V100 is used in the Summit system deployed at Oak Ridge National Laboratory, and the new exascale system, Frontier, will be using AMD Radeon Instinct GPU. We evaluated the performance of various optimization techniques on test matrices, ranging in size from100K to 4M, that are representative of ScIFEN applications. The linear solver computation is memory-bound on both GPUs. Our experiments show that on the NVIDIA GPU we obtained up to79%of the theoretical peak bandwidth, while the AMD GPU achieved 59%. Overall, the NVIDIA V100 GPU outperforms the AMD MI 25 GPU1. We observed an overall speedup of up to37X on an NVIDIA V100 compared to an Intel Skylake 12-coremachine. The solver for a 4M degree of freedom system took under 2.5 seconds.

Mohammad Zubair↗

Green Design of a Cellulosic Butanol Supply Chain Network: A Case Study of Sorghum Stem Bio-butanol in Missouri

Cellulosic butanol is a very promising renewable fuel to consider for the future transportation market. However, the seasonal availability of the raw materials, high maintenance cost, and high logistical cost of the biomass energy supply chain are the main factors impeding the commercialization and large-scale-production of this energy source. Furthermore, research focusing on an environmental or green supply chain network design of cellulosic butanol has been insufficient. Here, this study focused on designing a green supply chain network for cellulosic butanol. A life cycle analysis was integrated into a multi-objective linear programming model to optimize the cellulosic butanol supply chain network. With the objectives of maximizing the economic profits and minimizing the greenhouse gas emissions, the proposed model can optimize the location and size of a bio-butanol production plant. The mathematical model was applied to a case study in the state of Missouri, and solved the tradeoff between the feedstock and market availabilities of sorghum stem bio-butanol.

09 BIOMASS FUELS↗

Printable Fiber Reinforced Cement Composites – Feasibility Study

Additive manufacturing is enabling the manufacturability of structures with previously unattainable complexity or functionality, and there is growing interest in additive manufacturing of “printed” concrete structures. The focus of this Phase 1 Technical Collaboration (TC) project was to evaluate feasibility of printing hybrid cement composite structures reinforced with textile carbon fibers (tCF). This project leverages other (non-IACMI) projects on cement formulation and printing process development, as well as on the production process for tCF. This project’s primary focus was to explore cement composite mix design with textile carbon fibers to be manufactured by MonteFibre (TC partner) and evaluate suitable fiber-matrix interface or sizing for cement composites working with Michelman (TC Partner). This project supports IACMI’s goal of reducing the cost and embodied energy of carbon fiber composites. Cost is one of the fundamental challenges to carbon fiber reinforced cement composites. Cement is an extremely inexpensive material (approximately $\$$0.05/lb). Adding 1 wt% of conventional carbon fiber to cement quadruples its cost. Therefore, the need to use low-cost carbon fiber and ensure that the additional cost of the carbon fiber has a greater cost benefit to the final product. This was the first preliminary evaluation to integrate tCF reinforcement in cement composites, and such potential tCF utilization should significantly reduce materials cost. Cement composite production is energy and emissions intensive, thus by strengthening it less material will be required. Hence, the embodied energy and production time of the resulting structures will be reduced. Additionally, integrating these new materials into additive processes can enable selective use of the material in high load or stress areas. It is noteworthy that past work in this field of fiber reinforced cement composites did not consider the optimization of fiber-matrix interface using suitable sizing. Carbon fiber reinforcement offers potential added benefits of thermal conductivity (which affects cure rate) and flow behavior that could provide opportunities for site specific utilization of carbon fiber on hybrid cement structures (e.g. use the fiber reinforcement on outer surfaces to enhance strength and modulus and then infiltrating the internal structures with conventional concrete). MonteFibre was the industry lead and planned on supplying the tCF for this project. However, during the short Phase-1 duration of this project, MonteFibre was unable to produce tCF for this project due to manufacturing plant being off-line throughout the course of the project. The project team decided to pursue an alternate option which involved demonstrating printable concrete with steel fibers by the ORNL lead, Dr. Brian Post. The University of Tennessee collaboration team focused on evaluating the suitable chemical sizing for carbon fibers working with Michelman and also developed methods for material characterization of cement-based composites to evaluate the material response for compression, shear, flexure, and tension. The two milestones for University of Tennessee, Knoxville were realized related to identification of one sizing suitable for carbon fiber reinforced cement composite and developing data associated with mechanical behavior of unreinforced (neat) and carbon fiber reinforced cement composites. ORNL could not complete the task of carbon fiber reinforced printed cement composites due to the reasons mentioned earlier, but was able to replace tCF with steel fibers to demonstrate the feasibility of printing with fiber reinforced cement composites. The Project team reviewed possible sizing chemistry available in collaboration with Michelman for use on carbon fiber reinforcement in cement composites and concrete applications. Our initial goal was to identify a sizing most promising for formulation with textile carbon fibers (tCF) to deliver excellent mechanical properties in composite material state. Since tCF was not available for this project as originally envisioned, the team continued this task to identify a suitable sizing for carbon fiber applications by applying such sizing to lower cost carbon fibers currently available commercially from Zoltek called Panex fibers. At a future time this can be optimized for textile carbon fibers from Montefibre. The bulk of previous work on carbon fiber reinforced cement has neglected the importance of fiber-matrix adhesion on mechanical properties of the cement composite and identifying this missing link was an important accomplishment for future research. Tensile behavior of fiber reinforced concrete is important to evaluate in order to realize the dream of concrete products that do not need reinforcing steel. Important sample preparation and testing procedures were addressed in this study and it was concluded that substantial improvements in tensile behavior, without compromising compressive strength, and improved ductility can result from the use of carbon fiber reinforcement.

36 MATERIALS SCIENCE↗

Simulating self-powered neutron detector responses to infer burnup-induced power distribution perturbations in next-generation light water reactors

Understanding how 3D power distribution will be monitored throughout reactor core volumetric space in next-generation nuclear power reactors is crucial to the design, deployment, and licensing of these reactors. Although numerous techniques exist for 3D power distribution monitoring based on the response of both in situ and ex situ sensors currently implemented or proposed for use in the US reactor fleet, crucial details about these techniques are often unclear. The publicly available documentation does not include information such as how well these techniques are characterized and optimized in their implementations and the levels of uncertainty in the inferred 3D power distribution. The work described herein investigated a recently developed 3D power distribution inferencing method as applied to two next-generation reactor simulations: (1) the NuScale small modular reactor design and (2) the Westinghouse AP1000 design, both of which contain in-core strings of vanadium self-powered neutron detectors (SPNDs). This investigation considered a range of SPND string sensor densities, as well as a range of 3D power distribution axial segment sizes. In this work, SPND response simulation is informed by neutron flux calculations in representative homogenized cores. For the different sensor densities and power distribution axial segment sizes in these simulations, the average solution error, solver iterations, and run time were tracked to parameterize the sensor-core configuration.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

Machine Learning Solutions for a Stable Grid Recovery

Grid operating security studies are typically employed to establish operating boundaries, ensuring secure and stable operation for a range of operation under NERC guidelines. However, if these boundaries are severely violated, existing system security margins will be largely unknown, as would be a secure incremental dispatch path to higher security margins while continuing to serve load. As an alternative to the use of complex optimizations over dynamic conditions, this work employs the use of machine learning to identify a sequence of secure state transitions which place the grid in a higher degree of operating security with greater static and dynamic stability margins. Several reinforcement learning solution methods were developed using deep learning neural networks, including Deep Q-learning, Mu-Zero, and the continuous algorithms Proximal Reinforcement Learning, and Advantage Actor Critic Learning. The work is demonstrated on a power grid with three control dimensions but can be scaled in size and dimensionality, which is the subject of ongoing research.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗