SEARCH · Engineering Papers

Results for “multiple interface dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2018

A unified framework for polycrystal plasticity with grain boundary evolution

Plastic deformation in polycrystals is governed by the interplay between intra-granular slip and grain boundary-mediated plasticity. However, while the role played by bulk dislocations is relatively well-understood, the contribution of grain boundaries (GBs) has only recently begun to be studied. GB plasticity is known to play a key role along with bulk plasticity under a wide range of conditions, such as dynamic recovery, superplasticity, severe plastic deformation, etc., and developing models capable of simultaneously capturing GB and bulk plasticity has become a topic of high relevance. In this paper we develop a thermodynamically-consistent polycrystal plasticity model capable of simulating a variety of grain boundary-mediated plastic processes in conjunction with bulk dislocation slip. The model starts from the description of a single crystal and creates lattice strain-free polycrystalline configurations by using a specially-designed multiplicative decomposition developed by the authors. This leads to the introduction of a particular class of geometrically necessary dislocations (GND) that define fundamental GB features such as misorientation and inclination. The evolution of the system is based on an energy functional that uses a non-standard function of the GND tensor to account for the grain boundary energy, as well as for the standard elastic energy. Our implementation builds on smooth descriptions of GBs inspired on diffuse-interface models of grain evolution for numerical convenience. Here we demonstrate the generality and potential of the methodology by simulating a wide variety of phenomena such as shear-induced GB sliding, coupled GB motion, curvature-induced grain rotation and shrinkage, and polygonization via dislocation sub-grain formation.

36 MATERIALS SCIENCE↗

NASA NTRS2017

Milli-Watt Radioisotope Power to Enable Small, Long-Term Robotic “Probe” Space Exploration

Milli-watt Radioisotope Power Systems (RPS) based on Radioisotope Heater Units (RHUs) could be an ideal power source for certain spacecraft that cannot use solar power due to large distances from the sun, or other environmental constraints, and where they enable or significantly enhance the ability of a mission to meet its scientific or operational goals. Various modular, compact RHU-based thermoelectric (TE) generator concepts developed or derived from current NASA Small Business Innovation Research (SBIR) projects satisfying this need have been investigated. These modular, compact and low mass power systems could support small, highly-mobile robotic exploration packages, and could be incorporated into different robotic package concepts, spacecraft or satellites. Current modular RPS design concepts with 40mW, 80mW and 120mW power levels use RHUs and Bi2Te3 TE converters. Skutterudites materials could be used in the future if new higher thermal energy output and higher temperature miniature heat sources were developed, for example, using technologies currently in the General Purpose Heat Source (GPHS) used in higher electric power output RTGs. Small (a.k.a., “mice-like”) robotic packages could effectively utilize these RHU-driven power levels to accommodate crawling, climbing, monitoring, taking measurements, and communicating during long-term planetary missions aimed at gathering environmental and geologic data (i.e., over multiple decades). Waste heat from the cold side of the TE converter could also be directed toward the electronics and / or energy storage (e.g. batteries) to keep them within design temperature ranges. In addition to power generation and electronics / battery heating, the RHU / TE configuration could be designed to survive an external 500°C bake out procedure for critical spacecraft sterilization, environmental certification and planetary protection. Analytical studies have been performed to optimize various design configurations for power, mass, volume and robotic mobility. Specific power (mW/kg) and volumetric specific power (mW/cm3) characteristics of various design configurations will be presented and key conceptual design tradeoffs will be discussed. Hot- and cold-side thermal interfaces required to meet power and mass goals and associated design sensitivities will also be discussed. RHU / TE systems must overcome critical design challenges to survive high-g loadings in some robotic applications and we will examine the mass impacts required to satisfy various dynamic loading environments up to 10,000 g’s. Power can be generated for a minimum of 30 years or more using plutonium-238 dioxide heat sources (given that Pu-238 has an 87.7 year half-life) with some reduction in power as the heat source naturally degrades.

Bahrami, Poyan↗

DOE OSTI2021

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Ir(hkl) Surface Electrochemistry in a Nonadsorbing Acidic Medium

The fundamental properties of electrochemical materials depend on the multiple and often complex interactions between electrode surface sites and electrolyte species at the electrochemical interface. Despite Iridium use in electrolyzer systems, much of its surface electrochemistry remains underexplored. This study investigates the surface electrochemistry of Ir(111), Ir(100), and Ir(110) surfaces in acidic media. Using cyclic voltammetry and CO charge displacement experiments, we establish the charge states and adsorbate coverages as a function of the electrode potential, revealing the presence of hydrogen and hydroxyl co-adsorption at low potentials on (111), and almost no coverage of H ad on (110) facet. In situ Shell Isolated Nanoparticle Enhanced Raman Spectroscopy experiments provide direct evidence of the formation of key adsorbate species, such as hydrogen, hydroxyl, and oxygen, but most importantly, their interactions with interfacial water, confirmed by Density Functional Theory calculations. Our findings highlight the role of co-adsorption and interspecies interactions, with microkinetic adsorption voltammetry simulations corroborating the influence of lateral interactions on adsorption dynamics, particularly for Ir(100) where the OHad formation occurs as a sharp adsorption/desorption current. Our results underscores the importance of interfacial water and hydrogen bonding networks in shaping the electrochemical behavior on Ir surfaces, refining our baseline understanding of the Ir surface electrochemistry necessary for the development of advanced Ir-based electrochemical materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Toward Multi-scale Modeling and simulation of conduction in heterogeneous materials

This report summarizes a project in which the authors sought to develop and deploy: (i) experimental techniques to elucidate the complex, multiscale nature of thermal transport in particle-based materials; and (ii) modeling approaches to address current challenges in predicting performance variability of materials (e.g., identifying and characterizing physical- chemical processes and their couplings across multiple length and time scales, modeling information transfer between scales, and statically and dynamically resolving material structure and its evolution during manufacturing and device performance). Experimentally, several capabilities were successfully advanced. As discussed in Chapter 2 a flash diffusivity capability for measuring homogeneous thermal conductivity of pyrotechnic powders (and beyond) was advanced; leading to enhanced characterization of pyrotechnic materials and properties impacting component development. Chapter 4 describes success for the first time, although preliminary, in resolving thermal fields at speeds and spatial scales relevant to energetic components. Chapter 7 summarizes the first ever (as far as the authors know) application of TDTR to actual pyrotechnic materials. This is the first attempt to actually characterize these materials at the interfacial scale. On the modeling side, new capabilities in image processing of experimental microstructures and direct numerical simulation on complicated structures were advanced (see Chapters 3 and 5). In addition, modeling work described in Chapter 8 led to improved prediction of interface thermal conductance from first principles calculations. Toward the second point, for a model system of packed particles, significant headway was made in implementing numerical algorithms and collecting data to justify the approach in terms of highlighting the phenomena at play and pointing the way forward in developing and informing the kind of modeling approach originally envisioned (see Chapter 6). In both cases much more remains to be accomplished.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Segregation competition and complexion coexistence within a polycrystalline grain boundary network

Interfacial segregation can stabilize grain structures and even lead to grain boundary complexion transitions. However, understanding of the complexity of such phenomena in polycrystalline materials is limited, as most studies focus on bicrystal geometries. In this work, we investigate interfacial segregation and subsequent complexion transitions in polycrystalline Cu-Zr alloys using hybrid Monte Carlo/molecular dynamics simulations. No significant change in the grain size or structure is observed upon Zr dopant addition to a pure Cu polycrystal at moderate temperature, where grain boundary segregation is the dominant behavior. Segregation within the boundary network is inhomogeneous, with some boundaries having local concentrations that are an order of magnitude larger than the global value and others having almost no segregation, and changes to physical parameters such as boundary free volume and energy are found to correlate with dopant concentration. Further, another alloy sample is investigated at a higher temperature to probe the occurrence of widespread transitions in interfacial structure, where a significant fraction of the originally ordered boundaries transition to amorphous complexions, demonstrating the coexistence of multiple complexion types, each with their own distribution of boundary chemical composition. Overall, this work highlights that interfacial segregation and complexion structure can be diverse in a polycrystalline network. Here, the findings shown here complement existing computational and experimental studies of individual interfaces and help pave the way for unraveling the complexity of interfacial structure in realistic microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Cross-scale Efficient Tensor Contractions for Coupled Cluster Computations Through Multiple Programming Model Backends

Coupled-cluster methods provide highly accurate models of molecular structure by explicit numerical calculation of tensors representing the correlation between electrons. These calculations are dominated by a sequence of tensor contractions, motivating the development of numerical libraries for such operations. While based on matrix-matrix multiplication, these libraries are specialized to exploit symmetries in the molecular structure and in electronic interactions, and thus reduce the size of the tensor representation and the complexity of contractions. The resulting algorithms are irregular and their parallelization has been previously achieved via the use of dynamic scheduling or specialized data decompositions. We introduce our efforts to extend the Libtensor framework to work in the distributed memory environment in a scalable and energy efficient manner. We achieve up to 240 speedup compared with the best optimized shared memory implementation. We attain scalability to hundreds of thousands of compute cores on three distributed-memory architectures, (Cray XC30&XC40, BlueGene/Q), and on a heterogeneous GPU-CPU system (Cray XK7). As the bottlenecks shift from being compute-bound DGEMM's to communication-bound collectives as the size of the molecular system scales, we adopt two radically different parallelization approaches for handling load-imbalance. Nevertheless, we preserve a uni ed interface to both programming models to maintain the productivity of computational quantum chemists.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1991

Implementation issues for a compact 6 degree of freedom force reflecting handcontroller with cueing of modes

Teleoperated control requires a master human interface device that can provide haptic input and output which reflect the responses of a slave robotic system. The effort reported in this paper addresses the design and prototyping of a six degree-of-freedom (DOF) Cartesian coordinate hand controller for this purpose. The device design recommended is an XYZ stage attached to a three-roll wrist which positions a flight-type handgrip. Six degrees of freedom are transduced and control brushless DC motor servo electronics similar in design to those used in computer controlled robotic manipulators. This general approach supports scaled force, velocity, and position feedback to aid an operator in achieving telepresence. The generality of the device and control system characteristics allow the use of inverse dynamics robotic control methodology to project slave robot system forces and inertias to the operator (in scaled form) and at the same time to reduce the apparent inertia of the robotic handcontroller itself. The current control design, which is not multiple fault tolerant, can be extended to make flight control or space use possible. The proposed handcontroller will have advantages in space-based applications where an operator must control several robot arms in a simultaneous and coordinated fashion. It will also have applications in intravehicular activities (within the Space Station) such as microgravity experiments in metallurgy and biological experiments that require isolation from the astronauts' environment. For ground applications, the handcontroller will be useful in underwater activities where the generality of the proposed handcontroller becomes an asset for operation of many different manipulator types. Also applications will emerge in the Military, Construction, and Maintenance/Manufacturing areas including ordnance handling, mine removal, NBC (Nuclear, Chemical, Biological) operations, control of vehicles, and operating strength and agility enhanced machines. Future avionics applications including advanced helicopter and aircraft control may also become important.

Jacobus, Heidi↗

DOE OSTI2025

Tactical Analysis for Calculating Contextual Risk at Boundaries: Summary of Laboratory Directed Research & Development Effort

The Tactical Analysis for Calculating Contextual Risk at Boundaries (TACCRAB) tool is an innovative digital twin (DT) platform and automated risk algorithm designed to transform operational decision-making in structured screening environments, with an initial focus on Southern Border Land Ports of Entry (POEs). The invention provides integration points for advanced artificial intelligence, predictive modeling, and real-time data analysis to produce a comprehensive risk management tool that enables proactive, data-informed security strategies. The core inventive features of TACCRAB center on its unique risk algorithm, which dynamically calculates contextual risk by synthesizing historical data, near real-time streaming data from the checkpoints themselves, and AI-generated predictions. Unlike traditional risk assessment methods, TACCRAB utilizes a DT to provide comprehensive operational insights, allowing stakeholders to visualize, simulate, and optimize checkpoint configurations with unprecedented speed and contextual awareness. TACCRAB's key innovation lies in its ability to combine multiple complex inputs - including technology detection probabilities, resource availability, screening pathway characteristics, and threat actor behavioral patterns - into a unified risk calculation and update these inputs based on changing operational and environmental conditions. By leveraging a DT that continuously updates and learns from linked data, TACCRAB can suggest adaptive mitigation strategies that minimize risk while maintaining operational efficiency. Particularly novel is the platform's approach to decision support, which goes beyond static risk assessment. The DT provides dynamic metrics such as wait times, resource allocation effectiveness, and potential emerging threat scenarios, enabling users to view sophisticated, relevant what-if simulations and optimize checkpoint operations in near real-time. The system's architecture allows for generalized application across different screening environments, such as secure facilities, ports of entry, and soft targets, making it a versatile tool for security and operational management. The invention distinguishes itself through its comprehensive integration of predictive modeling, AI-driven pattern discovery, and user-friendly interface design. By combining these elements, TACCRAB transforms complex risk data into actionable insights, supporting decision-makers at various organizational levels - from booth agents making split-second screening decisions to checkpoint managers optimizing the day's resource allocation to strategic planners managing long-term investments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS2003

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

DOE OSTI2019

Soft costs and knowledge pathways in the US solar photovoltaic ecosystem

An enduring puzzle in the US solar industry is that despite substantial reductions in the cost of photovoltaic hardware, soft costs (defined here as all costs not directly driven by hardware) have not done the same (US DOE, 2012). In fact, soft costs now account for the majority of the cost of total cost of installed solar systems, despite decades of efforts by firms and policy-makers to reduce these costs. In this study, we address two questions: a) what drives soft costs within the US solar industry, and b) how do firms gain the knowledge required to address these costs? We address this question through an explorative inductive multiple-case study (e.g., Eisenhardt, 1989; Appleyard, 1996). Our sample consists of 15 actors from across various roles within the industry, including residential and commercial installers (8), distribution intermediaries (3), and business software vendors (4). We conducted extensive interviews with executives at each firm (76 interviews in total), and supplemented this information with archival resources such as press releases, news coverage, analyst reports, and internet articles. This data was synthesized into a comprehensive case for each actor, detailing the activities conducted within each firm, the drivers of soft costs within that firm, and knowledge flows (e.g., new technologies, best practices, etc) within and across firms. We then compared across these cases to identify common patterns in cost and knowledge dynamics. Our study contributes a nuanced perspective on the drivers of soft costs within the solar industry as well as the patterns and mechanisms of knowledge flows within the industry. We find for example that the majority of the soft costs occur at the interface between firms (e.g., utilities, suppliers). We also observe substantive variation in costs across different types of actors and even within firm types (e.g., small vs. large installers). We also provide insight into how firms acquire knowledge and when they are able to apply it in order to reduce these costs. Our data reveal heterogeneity in a) the types of knowledge, and b) the pathways of knowledge flow between the actors in the industry. For example, we observe that distributors and suppliers play a key role in the dissemination of new technologies as well as general management knowledge for smaller firms. In contrast, informal personal networks and internal experimentation appear to be more important for larger firms. We also reveal critical distinctions across markets – e.g., differences in organizational culture that facilitate learning and experimentation more readily in residential vs. commercial markets. Overall, these results contribute preliminary insights into the knowledge dynamics of the solar industry, the persistence of soft costs over time, heterogeneity across firm types and regions, and the effective design of policy to support the solar industry.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2024

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗

NASA NTRS2013

Tracker Toolkit

This software can track multiple moving objects within a video stream simultaneously, use visual features to aid in the tracking, and initiate tracks based on object detection in a subregion. A simple programmatic interface allows plugging into larger image chain modeling suites. It extracts unique visual features for aid in tracking and later analysis, and includes sub-functionality for extracting visual features about an object identified within an image frame. Tracker Toolkit utilizes a feature extraction algorithm to tag each object with metadata features about its size, shape, color, and movement. Its functionality is independent of the scale of objects within a scene. The only assumption made on the tracked objects is that they move. There are no constraints on size within the scene, shape, or type of movement. The Tracker Toolkit is also capable of following an arbitrary number of objects in the same scene, identifying and propagating the track of each object from frame to frame. Target objects may be specified for tracking beforehand, or may be dynamically discovered within a tripwire region. Initialization of the Tracker Toolkit algorithm includes two steps: Initializing the data structures for tracked target objects, including targets preselected for tracking; and initializing the tripwire region. If no tripwire region is desired, this step is skipped. The tripwire region is an area within the frames that is always checked for new objects, and all new objects discovered within the region will be tracked until lost (by leaving the frame, stopping, or blending in to the background).

Lewis, Steven J.↗

DOE OSTI2024

In Situ Characterization of Interface Evolution in Argyrodite‐Based All‐Solid‐State Li Batteries

Interfacial stability is one of the critical challenges in all-solid-state Li metal batteries. Multiple processes such as solid electrolyte (SE) decomposition and lithium dendrite growth take place at the solid interfaces during cycling, leading to the overall cell failure. To deconvolute these complex processes, in situ characterization is of paramount importance to elucidate the interfacial evolution on the SE upon Li plating/stripping. Herein, an all-solid-state asymmetric in situ cell is developed that allows the direct visualization of the highly localized Li plating/stripping processes under the optical microscope. Moreover, this cell configuration enables reliable post-mortem chemical and morphological analysis of the intact SE/Li interface. Using combined scanning electron microscopy and energy-dispersive X-ray spectroscopy, the study reveals that the evolution of the Li argyrodite interface is strongly influenced by the current density, particularly in terms of chemical distribution and Li plating morphology. More specifically, the solid interface is LiCl-rich with the formation of Li cubes at low current densities, while high currents result in more uniform elemental distribution and filament morphology. These findings elucidate the dynamic evolution mechanism at solid interfaces and offer valuable guidance for developing stable solid interfaces in all-solid-state Li metal batteries.

Huang, Di↗

DOE OSTI2015

DYNAMICS OF ON-DISK PLUMES AS OBSERVED WITH THE INTERFACE REGION IMAGING SPECTROGRAPH, THE ATMOSPHERIC IMAGING ASSEMBLY, AND THE HELIOSEISMIC AND MAGNETIC IMAGER

We examine the role of small-scale transients in the formation and evolution of solar coronal plumes. We study the dynamics of plume footpoints seen in the vicinity of a coronal hole using the Atmospheric Imaging Assembly (AIA) images, the Helioseismic and Magnetic Imager magnetogram on board the Solar Dynamics Observatory and spectroscopic data from the Interface Region Imaging Spectrograph (IRIS). Quasi-periodic brightenings are observed in the base of the plumes and are associated with magnetic flux changes. With the high spectral and spatial resolution of IRIS, we identify the sources of these oscillations and try to understand what role the transients at the footpoints can play in sustaining the coronal plumes. IRIS “sit-and-stare” observations provide a unique opportunity to study the evolution of footpoints of the plumes. We notice enhanced line width and intensity, and large deviation from the average Doppler shift in the line profiles at specific instances, which indicate the presence of flows at the footpoints of plumes. We propose that outflows (jet-like features) as a result of small-scale reconnections affect the line profiles. These jet-like features may also be responsible for the generation of propagating disturbances (PDs) within the plumes, which are observed to be propagating to larger distances as recorded from multiple AIA channels. These PDs can be explained in terms of slow magnetoacoustic waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Cross-scale efficient tensor contractions for coupled cluster computations through multiple programming model backends

Coupled-cluster methods provide highly accurate models of molecular structure through explicit numerical calculation of tensors representing the correlation between electrons. These calculations are dominated by a sequence of tensor contractions, motivating the development of numerical libraries for such operations. While based on matrix–matrix multiplication, these libraries are specialized to exploit symmetries in the molecular structure and in electronic interactions, and thus reduce the size of the tensor representation and the complexity of contractions. The resulting algorithms are irregular and their parallelization has been previously achieved via the use of dynamic scheduling or specialized data decompositions. We introduce our efforts to extend the Libtensor framework to work in the distributed memory environment in a scalable and energy-efficient manner. We achieve up to 240× speedup compared with the optimized shared memory implementation of Libtensor. We attain scalability to hundreds of thousands of compute cores on three distributed-memory architectures (Cray XC30 and XC40, and IBM Blue Gene/Q), and on a heterogeneous GPU-CPU system (Cray XK7). As the bottlenecks shift from being compute-bound DGEMM's to communication-bound collectives as the size of the molecular system scales, we adopt two radically different parallelization approaches for handling load-imbalance, tasking and bulk synchronous models. Nevertheless, we preserve a unified interface to both programming models to maintain the productivity of computational quantum chemists.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Approaches for Synthesis and Deployment of Controller Models on Automated Vehicles for Car-following in Mixed Autonomy

This paper describes the software design patterns and vehicle interfaces that were employed to transition vehicle controllers from simulation environments to open-road field experiments. The approach relies on a life cycle that utilizes model-based design and code generation, along with agile software development, and both software and hardware-in-the-loop testing, with additional safety margins. Autonomous designs should consider the dynamics of mixed autonomy in traffic to safely operate among humans. The software that provides a vehicle’s behavior intelligence is often developed through simulation, which may have a mismatch between dynamics, or as a result of a reinforcement learning workflow, which may be a black box with challenges to analyze. In each of these cases, it is important to have research interfaces that provide strongly typed data streams accessible to researchers who are not software experts while continuing to satisfy safety and liveness constraints. This paper describes how we design the hardware platform interfaces and software design process for a mixed autonomy traffic experiment with a leader-follower scenario. Controller synthesis for these vehicles requires clearly articulated vehicle interfaces and software design patterns for successful onboard deployment. Testing strategies for such controllers are also described before algorithms are transitioned to full-scale field experiments with safety operators for the vehicles. Testing strategies include software-in-the-loop simulation testing, hardware-in-the-loop simulation, ghost-car testing, and read-only testing in live traffic. With our approach, we were not only able to validate our controller synthesized in scripts and simulation, but also able to scale deployment to multiple vehicles.

Bhadani, Rahul↗

DOE OSTI2017

Mechanisms of interfacial reactivity in near surface and extreme environments

The local water structure surrounding ions in aqueous solutions greatly affects their chemical properties such as reaction rates, ion association, and proton and electron transport. These properties result in the behavior of ions in natural aqueous environments. For example ore transport is facilitated by chloride ion pair formation and the reaction of ions in an interface is strongly dependent on the dehydration of the ion hydration shell. We are developing the use of high-resolution XAFS observations and 1st principles based MD-XAFS analysis (spectra simulated using 1st principle methods with no adjustable parameters, AIMD) to interpret the solution properties of strongly interacting aqueous solutes under arbitrary pressure and temperature conditions. In the 1st principle MD-XAFS method density functional theory (DFT) based MD simulations(Car and Parrinello, 1985) are used to generate a large ensemble of structural snap shots of the hydration region. These are then used to generate scattering intensities. I emphasize three points about this novel approach to analyzing XAFS data. 1st: As illustrated in Figure 1, the level of agreement between the calculated and observed intensities is considerably higher than has been obtained by any XAFS analysis to date (note 2nd shell region, R> 2 Å). 2nd: This result was obtained from a parameter free simulation with no fitting of the interaction potentials to any data. This supports the use of these methods for more difficult environments and more complex solutes (polyions). 3rd: New information about the shell structure (Figure 1) is now available because of this more detailed agreement. We note also that both multiple scattering and second shell features are well represented in the analysis. As far as we know this is the 1st analysis of second shell structure and multiple scattering. Excellent agreement has been obtained for most of the third row metal ions: Ca 2+ , Zn 2+ , Cu 2+ , Ni 2+ , Co 2+ , Mn 2+ , Fe 3+ , Cr 3+ . Calculations on these systems are demanding because of their open electronic shells, and high ionic charge. Principal Investigator: Professor John Weare (University of California, San Diego) The prediction of the interactions of geochemical fluids with minerals, nanoparticles, and colloids under extreme near surface conditions of temperature (T) and pressure (P) is a grand challenge research need in geosciences (U.S. DOE 2007, Basic Research Needs for Geosciences: Facilitating the 21st Energy Systems.). To evaluate the impact of these processes on energy production and management strategies it is necessary to have a high level of understanding of the interaction between complex natural fluids and mineral formations. This program emphasizes 1st principle parameter free simulations of complex chemical processes in solutions, in the mineral phase, and in the interfaces between these phases The development of new computational tools (with emphasis on oxide materials and reaction dynamics) tailored to treat wide range of conditions and time scales experienced in such geochemical applications is has been developed. Because of the sensitivity of the interaction in these systems to electronic structure and local bonding environments, and of the need to describe bond breaking/formation, our simulations are based on interactions calculated at the electronic structure level (ab-initio molecular dynamics, AIMD). The progress in the computational aspects of program may be summarized in terms of the following themes (objectives); Development of efficient parameter free dynamical simulation technology based on 1st principles force and energy calculations especially adapted for geochemical applications (e.g., mineral, interfaces and aqueous solutions) (continuing program); Calculation of the dynamics of water structure of in the surface-water interface of transition metal oxides and oxihydroxides; and Development of improved (beyond DFT+GGA) electronic structure calculations for minerals and the interface region that more accurately calculate electron correlation, spin density, and localization. The focus of the program is also on the iron oxide and oxihydroxide minerals and Fe 2+ (aq)/Fe 3+ (cr) oxidation in the mineral solution interface region. These methods included the development of model Hamilitonian methods that can be solved to near convergence for single site models (DMFT) and many-body perturbation methods (MP2, GW); Development of time decomposition methods to extend time scales of molecular dynamics (MD) simulations and support the use of high complexity electronic structure calculations (MP2, CCSD(T)) of forces for use in dynamical simulations where very high chemical accuracy is required (microsolvated reactions in absorbed surface layers); and The development of a new linear scaling finite element solver for eigenvalue problem that supports solution of quantum problems with unusual potential and boundary values. Application progress of the above new simulation technology to problems of geochemical interests includes; The prediction of metal oxide surface structure and the reduction/oxidation of Fe 3+ (cr)/Fe 2+ (aq) in metal oxide (hematite, goethite)/solution interfaces. Result: water interacts strongly with the 001 Hematite surface; interaction of water with the 100 goethite is weak; The study of ion solvation and the composition of ion hydration shells under extreme conditions (focus on Fe 3+/2+ , Al 3+ and Mg 2+ and their hydroxide speciation). Result: Ion association in water solutions can be calculated from 1st principle methods. Efficient sampling of the free energy requires more development; The continued development of new high resolution analysis of XAFS scattering of disordered systems (particularly Al, Mg) and of XANES calculations for aqueous ions. Result: EXAFS spectra can be calculated to high accuracy with DFT level dynamic simulations; The exploration of electron localization and electron transport in metal oxides (highly correlated materials). Result: Proper description of electron localization requires levels of calculation beyond DFT; and Localization of electrons in DFT type Hamiltonians was studied. Result: For very Dirac high exchange new solutions (New unphysical bifrucations) to the eigenvalue problem are found. The program was highly collaborative involving faculty and students in mathematics, physics and computer science departments as well as coworkers at the Pacific Northwest National Laboratories (PNNL). The students in this program had the opportunity to develop skills in the development of methods, the implementation of method on high performance parallel computers and the application of these methods to problem in geochemical science. Much of the software that was developed was incorporated in the NWchem software package maintained by PNNL.

58 GEOSCIENCES↗