Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic feature selection”

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 397 records · Page 22

From Pseudostreamer Jets to Coronal Mass Ejections:: Observations of the Breakout Continuum

The magnetic breakout model, in which reconnection in the corona leads to destabilization of a filament channel, explains numerous features of eruptive solar events, from small-scale jets to global-scale coronal mass ejections (CMEs). The underlying multipolar topology, pre-eruption activities, and sequence of magnetic-reconnection onsets (first breakout, then flare) of many observed fast CMEs/eruptive flares are fully consistent with the model. Recently, we have demonstrated that most observed coronal-hole jets in fan/spine topologies also are induced by breakout reconnection at the null point above a filament channel (with or without a filament).For these two types of eruptions occurring in similar topologies, the key question is, why do some events generate jets while others form CMEs? We focused on the initiation of eruptions in large bright points/small active regions that were located in coronal holes and clearly exhibited null-point (fan/spine) topologies: such configurations are referred to as pseudostreamers. We analyzed and compared SDO/AIA,SOHO/LASCO, and RHESSI observations of three events. Our analysis of the events revealed two new observable signatures of breakout reconnection prior to the explosive jet/CME outflows and flare onset: coronal dimming and the opening-up of field lines above the breakout current sheet. Most key properties were similar among the selected erupting structures, thereby eliminating region size, photospheric field strength, magnetic configuration, and pre-eruptive evolution as discriminating factors between jets and CMEs. We consider the factors that contribute to the different types of dynamic behavior, and conclude that the determining factor is the ratio of the magnetic free energy associated with the filament channel compared to the energy associated with the overlying flux inside and outside the pseudostreamer dome. Our results imply that a pseudostreamer eruption must be fast in order to produce an interplanetary CME.

Sun↗

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

Compiling Quantum Circuits for Dynamically Field-Programmable Neutral Atoms Array Processors

Dynamically field-programmable qubit arrays (DPQA) have recently emerged as a promising platform for quantum information processing. In DPQA, atomic qubits are selectively loaded into arrays of optical traps that can be reconfigured during the computation itself. Leveraging qubit transport and parallel, entangling quantum operations, different pairs of qubits, even those initially far away, can be entangled at different stages of the quantum program execution. Such reconfigurability and non-local connectivity present new challenges for compilation, especially in the layout synthesis step which places and routes the qubits and schedules the gates. In this paper, we consider a DPQA architecture that contains multiple arrays and supports 2D array movements, representing cutting-edge experimental platforms. Within this architecture, we discretize the state space and formulate layout synthesis as a satisfiability modulo theories problem, which can be solved by existing solvers optimally in terms of circuit depth. For a set of benchmark circuits generated by random graphs with complex connectivities, our compiler OLSQ-DPQA reduces the number of two-qubit entangling gates on small problem instances by 1.7x compared to optimal compilation results on a fixed planar architecture. To further improve scalability and practicality of the method, we introduce a greedy heuristic inspired by the iterative peeling approach in classical integrated circuit routing. Using a hybrid approach that combined the greedy and optimal methods, we demonstrate that our DPQA-based compiled circuits feature reduced scaling overhead compared to a grid fixed architecture, resulting in 5.1X less two-qubit gates for 90 qubit quantum circuits. These methods enable programmable, complex quantum circuits with neutral atom quantum computers, as well as informing both future compilers and future hardware choices.

Physics↗

Metal–Organic Frameworks as a Tunable Platform to Deconvolute Stereoelectronic Effects on the Catalytic Activity of Thioanisole Oxidation

The local environment of a metal active site plays an important role in affecting catalytic activity and selectivity. In re-cent studies, tailoring the behavior of a molybdenum-based active site via modulation of the first coordination sphere has led to improved thioanisole oxidation performance, but disentangling electronic effects from steric influences that arise from these modifications is nontrivial, especially in heterogeneous systems. To this end, the tunability of metal–organic frameworks (MOFs) makes them promising scaffolds for controlling the coordination sphere of a heterogene-ous, catalytically-active metal site while offering additional, attractive features such as crystallinity and high porosity. Herein, we report a variety of MOF-supported Mo species, which were investigated for catalytic thioanisole oxidation to methyl phenyl sulfoxide and/or methyl phenyl sulfone using tert-butyl hydroperoxide (tBHP) as the oxidant. In particular, MOFs of contrasting node architectures were targeted, presenting a unique opportunity to investigate ste-reoelectronic control of Mo active sites in a systematic manner. A Zr 6 -based MOF, NU-1000, was employed along with its sulfated analog Zr 6 -based NU-1000-SO 4 to anchor a dioxymolybdenum species, which enabled examination of sup-port-mediated active site polarizability on catalytic performance. In addition, a MOF containing a mixed metal node, Mo-MFU-4l, was used to probe the stereoelectronic impact of an N-donor ligand environment on the catalytic activity of the transmetalated Mo center. Furthermore, characterization techniques including single crystal X-ray diffraction, were concomi-tantly used with reaction kinetic profiles to better comprehend the dynamics of different Mo active sites, thus correlat-ing structural change with activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mars Ascent Vehicle Solid Propulsion Configuration

As part of a Mars Sample Return (MSR) campaign, two Mars Ascent Vehicle (MAV) configurations have been designed in parallel. Each ascent vehicle configuration has a different propulsion system, which ultimately leads to two unique vehicle designs. As part of a Preliminary Architecture Assessment (PAA), these vehicle designs were developed to the same level of maturity in order to inform the selection of one of the vehicles as the point of departure design for the campaign. The selection will be made in November 2019. Although the initial MSR architecture called for a hybrid-based propulsion MAV featuring solid wax fuel with liquid oxidizer, a configuration using more traditional solid propulsion was developed as an additional risk mitigation option. Though lacking in the single stage to orbit (SSTO) and throttle flexibility of a hybrid configuration, a solid configuration vehicle allows a simpler design with significantly longer flight heritage and higher Technology Readiness Level (TRL). This paper describes the design of the solid propulsion configuration. An additional paper will be published describing the design of the hybrid propulsion configuration. The solid propulsion configuration MAV was developed in 2019 by NASA Marshall Space Flight Center (MSFC) in association with NASA Jet Propulsion Laboratory (JPL). It features two stages with a modified STAR-17 motor for the second stage and a traditional electromechanical actuator Thrust Vector Controller (TVC). The vehicle was designed to deliver approximately 0.47kg of Martian geological samples to a circular orbit at Mars of 343km at a 25° inclination. Although solid motor designs in general are at a relatively high TRL, the integrated vehicle subsystems that operate in conjunction with these propulsion elements do not typically operate in a Martian environment, which in this application can get as cold as -40°C. The PAA advanced the maturity of these subsystems by performing detailed design and analysis on the vehicle with respect to structures and mechanisms, Guidance/Navigation/Control (GNC) systems, avionics, Reaction Control System (RCS), TVC, thermal environments, and advanced Computational Fluid Dynamics (CFD). This paper will summarize the results of these studies.

Yaghoubi, Darius↗

Mars Ascent Vehicle Solid Propulsion Configuration

As part of a Mars Sample Return (MSR) campaign, two Mars Ascent Vehicle (MAV) configurations have been designed in parallel. Each ascent vehicle configuration has a different propulsion system, which ultimately leads to two unique vehicle designs. As part of a Preliminary Architecture Assessment (PAA), these vehicle designs were developed to the same level of maturity in order to inform the selection of one of the vehicles as the point of departure design for the campaign. The selection will be made in November 2019. Although the initial MSR architecture called for a hybrid-based propulsion MAV featuring solid wax fuel with liquid oxidizer, a configuration using more traditional solid propulsion was developed as an additional risk mitigation option. Though lacking in the single stage to orbit (SSTO) and throttle flexibility of a hybrid configuration, a solid configuration vehicle allows a simpler design with significantly longer flight heritage and higher Technology Readiness Level (TRL). This paper describes the design of the solid propulsion configuration. An additional paper will be published describing the design of the hybrid propulsion configuration. The solid propulsion configuration MAV was developed in 2019 by NASA Marshall Space Flight Center (MSFC) in association with NASA Jet Propulsion Laboratory (JPL). It features two stages with a modified STAR-17 motor for the second stage and a traditional electromechanical actuator Thrust Vector Controller (TVC). The vehicle was designed to deliver approximately 0.47kg of Martian geological samples to a circular orbit at Mars of 343km at a 25° inclination. Although solid motor designs in general are at a relatively high TRL, the integrated vehicle subsystems that operate in conjunction with these propulsion elements do not typically operate in a Martian environment, which in this application can get as cold as -40°C. The PAA advanced the maturity of these subsystems by performing detailed design and analysis on the vehicle with respect to structures and mechanisms, Guidance/Navigation/Control (GNC) systems, avionics, Reaction Control System (RCS), TVC, thermal environments, and advanced Computational Fluid Dynamics (CFD). This paper will summarize the results of these studies.

Yaghoubi, Darius↗

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗

Seismoball: A Small Europa Orbiter Drop-Off Probe for Early Exploration of the Europan Surface

Recent magnetometry data received from Galileo indicate that the most likely explanation for the magnetic signature there is indeed a global conducting layer below the surface. This conducting layer is well matched by a salty, mineral rich strata beneath the Europan ice crust or a salt water ocean. Galileo imaging results show a variety of terrain types thought to contain young material; for example, lineaments, chaotic terrain, and eruption features. Additionally, Galileo images have shown indications of areas of up-welling where subsurface material periodically gets pushed to the surface due to the forces of fracturing, butting, and refreezing of the ice sheet. While Europa Orbiter will provide close-flyby high resolution images, as well as magnetometry, spectroscopy and other remote sensing data of the surface, it will not be able to provide essential engineering data like surface hardness and surface ice structure needed to support eventual landed missions. Additionally, ice chemical composition at microscopic scales can only be studied in detail through in situ instrumentation. Seismoball is a small probe designed to be injected into a surface intersect orbit around Europa. Using small reverse thrusters, the probe will be capable of nulling the high horizontal injection velocity as it approaches the 2 km surface injection altitude, thus allowing it to fall to the surface at an impact velocity of < 100m/sec (much less than the DS-2 impact velocities). The external breakaway thruster structure and crushable exterior shell absorb the impact energy while allowing the science instrument suite to remain intact. JPL has already started analyzing the entry dynamics and designing/building a small, low mass probe which will withstand the impact g-forces and fit as a 'carry-on' on board the Europa Orbiter. The probe will carry a suite of 5-6 micro-instruments for imaging the surface (both microscopic and far-field), surface and shallow subsurface ice temperatures, surface hardness, crustal dynamics and periodicity, and compositional chemistry. If selected, this flight development activity will provide a unique science opportunity and adjunct to the primary Orbiter science mission. The final flight system will be designed to accommodate orbiter mass, volume, and power interface constraints, as well as entry dynamics, g-load mitigation, and arbitrary landing orientation.

Tamppari, L.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Use of Design of Experiments and Rule-Based Inference in Determining Neural Network Architectures for Loss of Control Detection

In this work, we describe methods for selecting the neural network architectures and input spaces to implement belief state inference on generic commercial transport aircraft. First, we highlight a case study on the planning, execution, and analysis of a set of experiments to determine the configurations of a conditional variational autoencoder (CVAE). We present a structured method that can be used in a number of aerospace applications, to optimize the structure and training parameters of the CVAE for belief state inference, using Design of Experiments (DOE) statistical methodologies. The motivation for this specific DOE was to identify the appropriate hyperparameters for measuring the CVAE reconstruction probability and latent space, such that the measurements can be used to infer qualitative state changes for the aircraft. We demonstrate that this process yields information about a trained neural network’s utility for this specific application, along with a quantifiable range of certainty. We execute 84 experiments using loss-of-control flight maneuver data from a NASA T-2 aircraft, demonstrating that this empirical process allows us to construct cheap and simple models with specific attributes amenable to belief state inference in aerospace applications. While theoretically, we could create a single CVAE with an input space the size of all measurable flight variables and environmental dynamics, it becomes intractable to use such a neural network in an in-situ intelligent multi-agent system. Using the recommendations from our case study, we introduce a technical approach for feasibly describing the belief space by (1) identifying significant statistical relationships among flight variables using rule induction, (2) using a set of rules that cover all features to define the input space of multiple CVAEs, and (3) forming a belief space based on the joint probability density of their collective latent spaces. This results in a series of relatively small matrix multiplications that can be performed in real time, as opposed to large matrix computations in a single CVAE. We demonstrate the application of this approach on the T-2 flight loss-of control experiments, using the architecture and hyperparameter recommendations from the case study. We compare the utilities of an individual CVAE trained on all flight variables and multiple CVAEs defined on subsets of flight variables for detecting qualitative changes in flight. We demonstrate that the use of multiple CVAEs with smaller input spaces permits the CVAE to capture more granular relationships in the latent space, permitting better state space characterization and loss-of-control detection.

Design of experiments↗

Control, Fault Management, and Grid Support Functionality of an MV AC-DC Solid State Transformer based EV Extreme Fast Charging Station

Electric vehicles (EVs) have become increasingly popular in recent times while revolutionizing the consumer and commercial transportation market. The development of charging infrastructure has become one of the priorities for increasing the adoption of EVs. Extreme fast charging (XFC) technology can reduce the so-called ’range anxiety’ of consumers as they significantly reduce the charging time. With the advent of wide band-gap (WBG) power devices and improvement in power electronic converters, medium voltage (MV) solid state transformer (SST) based XFC system has the potential to replace the traditional XFC stations because of the lower footprint, ease of installation, enhanced control feature, and better system efficiency. The control system design is one of the critical aspects of the SST development process. Careful consideration and detailed analysis are required to find out suitable control method for the SST based on its topology among different centralized and decentralized control architectures. Also, the control parameters selection and potential improvement to the transient response of the controller ought to be investigated. Another major concern of the SST is different types of internal fault which reduces the overall reliability of the XFC system. As a result, designing a robust protection system is essential. Among different fault modes, open circuit switch faults have received significant attention as an active research area because of their likelihood and severe effects on converters. Therefore, the power stages used in the XFC system require functional and accurate open circuit switch fault management methods. An equally significant aspect of this SST based XFC is its compatibility in a microgrid where there is no synchronous generator present. When the grid is not available, the XFC SSTs can provide grid forming capability and continue supplying the critical loads in islanded mode. The transition between grid connected and islanded mode, especially the grid resynchronization process has to be carefully performed for the safety of the microgrid components. The challenges posed by the aforementioned issues have inspired the work done in this dissertation. Here, a 13.2 kV, 1 MVA, AC/DC SST for the XFC system is examined and a comparative analysis is conducted to select the control architecture based on feasibility of implementation and performance. A detailed control parameter design process is demonstrated considering the sensor dynamics and delay. The selected decentralized control method is augmented by introducing a novel sensor-less load current feedforward method to provide better voltage regulation at the DC bus during a change of load. Next, in the fault management section, a hierarchical failure mode effect analysis (FMEA) is proposed to enable a systematic design of the internal fault protection of the XFC SST as there are limited examples in the literature regarding the analysis of the safety and design of the protection of a power electronic converter system. Novel open circuit switch fault management methods for the converters in the system are presented. Finally, XFC SST based MV microgrid operations in grid connected mode and islanded mode are explored. A secondary control method for grid resynchronization is presented and a design process of control parameters is shown to ensure the stability of the secondary voltage and frequency regulation.

30 DIRECT ENERGY CONVERSION↗

A molecular line study of NGC 1333/IRAS 4

Molecular line surveys and fully sampled spectral line maps at 1.3 and 0.87 mm are used to examine the physical and chemical characteristics of the extreme Class I sources IRAS 4A and 4B in the L1450/NGC 1333 molecular cloud complex. A very well collimated, jetlike molecular flow emanates from IRAS 4A, with a dynamical age of a few thousand years. Symmetric, clumpy structure along the outflow lobes suggests that there is considerable variability in the mass-loss rate or wind velocity even at this young age. Molecular emission lines toward IRAS 4A and 4B are observed to be weak in the velocity range corresponding to quiescent material surrounding the young stellar objects (YSOs). Depletion factors of 10-20 are observed for all molecules, including CO, even for very conservative mass estimates from the measured millimeter and submillimeter dust continuum. However, abundance scaled with respect to CO are similar to other dark molecular cloud cores. Such depletions could be mimicked by high dust and optical depths or increased grain emissivities at the observing frequencies of 230 and 435 GHz, but the millimeter and submillimeter spectral energy distributions suggest that this is unlikely over the single-dish size scales of 5000-10,000 AU. Dense, outflowing gas is found to be kinematically, but not spatially, distinct from the quiescent material on these size scales. If CO is used as a chemical standard for the high-velocity gas, we find substantial enhancements in abundances of several molecules in outflowing material, most notably CS, SIO, and CH3OH. The SiO emission is kinematically well displaced from the bulk cloud velocity and likely arises from directly shocked material. As is the case for CO, however, the outflow features from more volatile species are centered near the cloud velocity and are often characterized by quite low rotational temperatures. We suggest that grain-grain collisions induced by velocity shear zones surrounding the outflow axes transiently desorb the grain mantles, resulting in large abundance enhancements of selected species. Similar results have recently been obtained in several other low-mass YSOs, where outflowing gas is often both kinematically and spatailly distinct, and are illustrative of the ability of accretion and outflow processes to simultaneously modify the composition of the gas and dust surrounding young stars.

Blake, Geoffrey A.↗

FY21 Progress Report: SRNL Analysis of ICCWR LCM and WAMS data for Corrosion and Cracking

The development of algorithms for machine learning and data analysis for the 3013 Surveillance Program is a collaborative effort by the Savannah River National Laboratory (SRNL) and the University of South Carolina (USC). For corrosion detection, Laser Confocal Microscope (LCM) or Wide Area 3D Measurement System (WAMS) data is extracted from large binary files, with software written to convert the data to physical attributes (e.g., height, color and grayscale values; all as functions of a location in a plane projection). A user-friendly Matlab Graphical User Interface (GUI) that reads data from either LCM or WAMS files was developed to integrate input data with software developed for processing and evaluation. The GUI can selectively download binary data, interrogate data attributes, label data, flag significant features, execute Machine Learning (ML) algorithms, output parameters for trained ML algorithms, report ML model accuracy with respect to labeled data, and generate graphical representations for various analyses. Features can be called out by user-specified thresholds, manual labeling or machine learning algorithms when they have been completed. The ability to rapidly label data is important because of the volume of data required for training machine learning algorithms. The GUI has the flexibility to allow addition of improved ML algorithms, methods for data visualization, and statistical computations. Statistical analyses via the GUI include areas of pits within a defined range of pit depths, correlations between Red-Green-Blue (RGB) or grayscale intensity and relative surface height, covariances between values associated with features, and feature histograms. The development of supervised machine learning algorithms, however, has been hindered by a lack of training data. The machine learning algorithms for crack identification are being refined but require improvements to the true positive rate for crack detection. This shortcoming is an artifact of the limited training data currently available, perhaps more so than the structure of the neural networks. At present, the best results are had from a consensus over an ensemble of randomly generated Deep Neural Network (DNN) or Convolutional Neural Network (CNN) algorithms. Although the consensus accuracy method has yielded optimum true positive and true negative rates in excess of 80%, additional validation testing is necessary. In addition to the suite of LCM data that was initially used, and which represents the majority of the work presented in this report, WAMS image data was also reviewed at a preliminary level. The review included a comparison between image resolution and dynamic range for each method. WAMS (ZON file) image data was found to have a pixel pitch of 3.69μm compared to 1 μm for the LCM (vk4 file) data, which implies a lower resolution for the WAMS images. Conversely, the ratio of dynamic range of the WAMS data to the LCM data was approximately 41:20 for height data, suggesting that information from WAMS should more accurately determine the depth of pits. At present, the significance of the greater dynamic range of the WAMS data relative to the LCM data has not yet been evaluated.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A simulation survey of galaxy interactions

Many carefully selected samples of interacting galaxies have been observed extensively in attempts to clarify whether interaction produces activity in galaxies. Because the sample members represent a wide range of encounter parameters and times, one can then study whether there are correlations between observable encounter features and, for example, Seyfert activity. On the other hand, in theoretical studies, simulations typically deal with either time-consuming detailed modelling of single galaxy pairs or tracing a few model encounters over time. The authors extend the observational survey approach by combining it with a simulation survey. The authors are conducting a survey of model encounters, covering the most important encounter parameters over a wide range. Some parameters, such as companion structure and initial velocity, are demonstratably less important and can be ignored in a first pass. The parameter range must be richly enough sampled so that the authors can evaluate the uniqueness of the observable morphology and velocity structure of the resulting simulated pairs to diagnose unobservable companion orbit parameters. They are using a self-gravitating polar n-body code run on the Cray X-MP at the Alabama Supercomputer Network. For each simulation, the authors have stellar and gas distributions predicted over, typically, a billion years, along with information on gas motions within the disk and any material captured by the companion or lost to the system. Features of disturbed spiral galaxies are sensitive enough to time and encounter parameters so that a match of the simulation survey results to observations can be applied as starting points to infer unobservable orbital or system parameters in actual sample members. This should enable them to examine whether interesting observed properties (Seyfert activity, nuclear star-formation rate) are functions of unobservable dynamical properties which characterize each encounter. Any correlations (or lack of some expected ones) will provide strong clues as to how or whether these phenomena are related to interactions. Aside from its use with such observed samples, this survey should greatly speed determination of initial orbital parameters for more detailed subsequent simulations of individual systems.

Byrd, Gene G.↗

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗

Illuminating the Photophysics of Oxygen Atom Transfer in Rare Earth Metal-Organic Complexes

Industrial separations of rare earth metals are reliant on selectivity based on ionic radii, which are time, energy, and resource consuming.1 The unique and discontinuous manifold of energy states associated with the f-electron configurations of the rare earth ions2 offer an opportunity to overcome some of these limitations by incorporating light into the separation process. Here, we use transient absorption spectroscopy to understand the differences in mechanism of light-initiated oxygen atom transfer (OAT) in yttrium (Y) and dysprosium (Dy) rare earth complexes. These metal-organic complexes are functionalized with hexfluoroacetyl-acetonate (hfac) and 4-methylmorpholine-N-oxide (NMMO). After light absorption in the presence of triphenyl phosphine (TPP), the NMMO in the Y-NMMO complex is replaced by the TPPO OAT reaction product. Because the Dy f-electron manifold is accessible to the ligand excited states, while there are no accessible states associated with f-electrons for Y, there are significant differences in reactivity. Transient absorption spectroscopy (TAS) was used to monitor changes after 340 nm photoexcitation of Dy/Y-NMMO or the product complex Dy/Y-TPPO. TAS is a pump-probe technique where a white light probe pulse is delayed relative to a visible pump pulse, which photoexcites a sample. Differences between the ground and excited states are monitored over time. With TAS, we access time scales from 100 fs-400 us to determine the photophysics of the complexes of interest. Photoexcitation of Y-NMMO (Fig. 1A) creates a broad photoinduced absorption (PIA) that evolves into a sharp positive feature over 4.2 ps centered at 450 nm. A secondary feature centered at 550 nm grows in with a 200 ps time constant and remains through the 5 ns time window of our experiment. Control experiments for Y-TPPO reveal that the feature at 450 nm is due to internal dynamics of the hfac ligand, while the 550 nm shoulder is due to the presence of the NMMO. Comparing kinetics at 550 nm for Y-NMMO and Dy-NMMO (Fig. 1B) shows that while Y exhibits a distinct growth in the NMMO feature, the Dy decays almost completely within 5 ns. Control experiments with Dy-TPPO prove that the ligand excited states centered on the hfac ligand are deactivated via energy transfer to the Dy manifold within 200 ps. This deactivation leads to significant differences in OAT reactivity between Y and Dy.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Mechanisms of Retention and Transport of Fission Products in Virgin and Irradiated Nuclear Graphite

To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗