SEARCH · Search NASA
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.
A computer program has been developed to extract three dimensional information of dislocations, interfaces, and precipitates/ cavities from two transmission electron microcopy (TEM) images.
Interface stabilization is critical to the development of working all‐solid‐state batteries. Rigid cathode/solid electrolyte interfaces often disintegrate due to anisotropic volume change of cathode‐active materials, resulting in irreversible capacity loss. Herein, we demonstrate that Li 1.211 Mo 0.467 Cr 0.3 O 2 (LMCO), a pioneering cation‐disordered rocksalt oxide (DRX) cathode that has intrinsically small volume change upon lithium intercalation, can be integrated with a thiophosphate‐based solid electrolyte for all‐solid‐state batteries. Interface stability of the all‐solid LMCO cell was investigated by electrochemical impedance spectroscopy, X‐ray micro‐computed tomography, and electron microscopy. Since LMCO was initially synthesized as a layered phase exhibiting a large volume change, interface disintegration can be observable in the early cycles. As layered LMCO phase‐transformed into DRX LMCO in subsequent cycles, reintegration of the interfaces occurs within a pressurized cell as a result of its zero‐stain behavior. Consequently, the DRX LMCO cathode maintains interface integrity, and thus electrical wiring, over an extended number of cycles, leading to improved capacity retention with small internal cell resistance.
Improvements in scientific instrumentation allow imaging at mesoscopic to atomic length scales, many spectroscopic modes, and now—with the rise of multimodal acquisition systems and the associated processing capability—the era of multidimensional, informationally dense data sets has arrived. Technical issues in these combinatorial scientific fields are exacerbated by computational challenges best summarized as a necessity for drastic improvement in the capability to transfer, store, and analyze large volumes of data. The Bellerophon Environment for Analysis of Materials (BEAM) platform provides material scientists the capability to directly leverage the integrated computational and analytical power of High Performance Computing (HPC) to perform scalable data analysis and simulation via an intuitive, cross-platform client user interface. This framework delivers authenticated, “push-button” execution of complex user workflows that deploy data analysis algorithms and computational simulations utilizing the converged compute-and-data infrastructure at Oak Ridge National Laboratory's (ORNL) Compute and Data Environment for Science (CADES) and HPC environments like Titan at the Oak Ridge Leadership Computing Facility (OLCF). In this work we address the underlying HPC needs for characterization in the material science community, elaborate how BEAM's design and infrastructure tackle those needs, and present a small sub-set of user cases where scientists utilized BEAM across a broad range of analytical techniques and analysis modes.
Many computationally unintensive models for material process-structure-property connections exist in the literature (e.g. journal articles and textbooks). However, to our knowledge there is no existing compilation of these models in one package. Many industrially and research relevant workflows involve combinations of many of these models and if the models have inconsistent interfaces, building these workflows can be challenging. Ramen provides a collection of computationally unintensive models for calculating material process-structure-property relationships from the literature, as well as related tools (e.g. for plotting results). These models use a consistent interface structure and all leverage material data files and data structures from the Mist library (also developed at ORNL). Ramen provides a structured package for combining many standard models in one place, with a standardized interface. This approach makes combining models in long/complex chains easy and straightforward.
The subsurface flow and reactive transport code PFLOTRAN can accommodate fracture–matrix interaction for multicomponent reactive transport using a Dual Continuum formulation. In this work we present two examples to assess this approach. Through these examples, we also illustrate how this Dual Continuum model can be parameterised and what kind of grid spacing is needed to accurately capture reactive transport in the secondary continuum. The results of the calculations show that when geochemical reactions in the matrix (or the secondary continuum) are included, a very fine grid spacing close to the fracture–matrix interface is required to properly capture geochemical gradients at the interface. We also show that one can reduce the computational cost by using a variable grid spacing, that limits the total number in the secondary continuum, and yet obtain accurate results. Furthermore, this feature makes the Dual Continuum approach suitable for large-scale reactive transport modelling in fractured crystalline rock.
Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.
CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.
To test the feasibility and safety of injecting a high-contrast hydrogel marker at the head of the pancreas (HOP) and duodenum interface and assesses the marker visibility on cone beam computed tomography (CBCT) to localize this important boundary during image guided radiation therapy in a porcine model.
Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.
Nontoxic, biocompatible nanodiamonds (ND) have recently been implemented in rational, systematic design of optimal therapeutic use in nanomedicines. However, hydrophilicity of the ND surface strongly influences structure and dynamics of biomolecules that restrict in situ applications of ND. Therefore, fundamental understanding of the impact of hydrophilic ND surface on biomolecules at the molecular level is essential. For tRNA, we observe an enhancement of dynamical behavior in the presence of ND contrary to generally observed slow motion at strongly interacting interfaces. We took advantage of neutron scattering experiments and computer simulations to demonstrate this atypical faster dynamics of tRNA on ND surface. The strong attractive interactions between ND, tRNA, and water give rise to unlike dynamical behavior and structural changes of tRNA in front of ND compared to without ND. As a result, our new findings may provide new design principles for safer, improved drug delivery platforms.
A seamless chip-to-world photonic interface enables broad advancements in optical ranging, display, communication, computation and quantum information science. The ideal solution enables two-dimensional scanning of a diffraction-limited beam from anywhere on a photonic integrated circuit to a large number of resolvable spots. Current beam-scanning technologies are limited by a fundamental trade-off: photonic-integrated-circuits with diffractive optics offer scalability but have poor mode quality, whereas inertially limited micromechanical scanners provide high-quality beams but lack scalable integration. Here we report a photonic ski-jump—a nanoscale waveguide monolithically integrated on a piezoelectric cantilever—to overcome these limitations. It passively curls ~90° out-of-plane within a less-than-0.1 mm 2 footprint, emits a submicrometre, broadband diffraction-limited beam, and exhibits kilohertz-rate mechanical resonances with quality factors of over 10,000. Fabricated in a volume complementary metal–oxide–semiconductor (CMOS) foundry, our device enables scalable two-dimensional beam scanning. Driven on-resonance at CMOS-level voltages, it achieves a footprint-adjusted spot rate of 68.6 mega spots s –1 mm–², exceeding state-of-the-art micro-electro-mechanical systems mirrors by more than 50-fold, which is sufficient for one million pixels at 100 Hz from an approximately 1.5 mm diameter footprint. We demonstrate full-colour image and video projection, and single-photon initialization and readout from silicon vacancy centres in diamond. Finally, by demonstrating uniformity across a 64 ski-jump array, we establish a pathway to achieving greater than one gigaspot resolution at kilohertz rates within a sub-5-cm-diameter footprint, creating a seamless optical pipeline between integrated photonic processors and the free-space world.
Remarkable success has been achieved in the discovery of ceramic alkali superionic conductors as electrolytes in solid-state batteries; however, obtaining a stable interface between these electrolytes and electrodes is difficult. Only limited studies on the compatibility between electrodes and solid electrolytes have been reported, partially because of the need for expensive instrumentation and special cell designs. Without simple yet powerful tools, these compatibility issues cannot be systematically investigated, thus hindering the generalization of design rules for the integration of solid-state battery components. Herein, we present a methodology that combines density functional theory calculations and simple experimental techniques such as X-ray diffraction, simultaneous differential scanning calorimetry and thermal gravimetric analysis, and electrochemistry to efficiently screen the compatibility of numerous electrode/electrolyte pairs. We systemically distinguish between the electrochemical stability of the solid-state conductor, which is relevant wherever the electrolyte contacts an electron pathway, and the electrochemical stability of the electrode/electrolyte interfaces. For the solid electrolyte, we are able to computationally derive an absolute thermodynamic stability voltage window, which is small for Na 3 PS 4 and Na 3 PSe 4 , and a larger voltage window which can be kinetically stabilized. The experimental stability, when measured with reliable techniques, falls between these thermodynamic and kinetic limits. Employing a Na solid-state system as an example, we demonstrate the efficiency of our method by finding the most stable system (NaCrO 2 |Na 3 PS 4 |Na–Sn) within a selected chemical space (more than 20 different combinations of electrodes and electrolytes). Important selection criteria for the cathode, electrolyte, and anode in solid-state batteries are also derived from this study. The current method not only provides an essential guide for integrating all-solid-state battery components but can also significantly accelerate the expansion of the electrolyte/electrode compatibility data.
The Production Operations Management System (POMS) software allows production teams and analysis groups across multiple Fermilab experiments to launch, modify and monitor large scale campaigns of related Monte Carlo or data processing jobs, and currently manages the majority of production computing of experiments at Fermilab.POMS provides a web service interface that enables automated job submission on distributed resources according to customers’ requests, as well as subsequent monitoring of those submissions as well as recovery of failed submissions, debugging and record keeping.POMS interfaces with existing HEP data access, processing, movement and monitoring tools at Fermilab, including Jobsub (for job submission), SAM (for data management), and Landscape (for monitoring), and in combination with them handles the creation of intermediate datasets of files for multistage campaigns.POMS visualizes these campaigns with an interactive GUI interface. An important feature of POMS is a one-to-one connection between the GUI campaign visualizer and a text based representation of the complete campaign configuration. An extensible library of template campaign configurations is available. The templates are user modifiable and map cleanly to and from the GUI description of the campaign.This paper will discuss the current usage, history, and some technical details of POMS, as well as planned future extensions, particularly interfacing with Rucio as well as SAM for file movement.
We present the simulation framework CRPropa version 3 designed for efficient development of astrophysical predictions for ultra-high energy particles. Users can assemble modules of the most relevant propagation effects in galactic and extragalactic space, include their own physics modules with new features, and receive on output primary and secondary cosmic messengers including nuclei, neutrinos and photons. In extension to the propagation physics contained in a previous CRPropa version, the new version facilitates high-performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques, an improved photonuclear interaction calculation, and propagation in time dependent environments to take into account cosmic evolution effects in anisotropy studies and variable sources. First applications using highlighted features are presented as well.
When pH and alkalinity increase, calcite frequently precipitates and hence modifies the petrophysical properties of porous media. The complex conductivity method can be used to directly monitor calcite precipitation in porous media because it is sensitive to the evolution of the mineralogy, pore structure and its connectivity. Here, we have developed a mechanistic grain polarization model considering the electrochemical polarization of the Stern and diffuse layers surrounding calcite particles. Our complex conductivity model depends on the surface charge density of the Stern layer and on the electrical potential at the onset of the diffuse layer, which are computed using a basic Stern model of the calcite/water interface. The complex conductivity measurements of Wu et al. on a column packed with glass beads where calcite precipitation occurs are reproduced by our surface complexation and complex conductivity models. The evolution of the size and shape of calcite particles during the calcite precipitation experiment is estimated by our complex conductivity model. At the early stage of the calcite precipitation experiment, modelled particles sizes increase and calcite particles flatten with time because calcite crystals nucleate at the surface of glass beads and grow into larger calcite grains. Finally, at the later stage of the calcite precipitation experiment, modelled sizes and cementation exponents of calcite particles decrease with time because large calcite grains aggregate over multiple glass beads and only small calcite crystals polarize.
Quantum computing vendors are beginning to open up application programming interfaces for direct pulse-level quantum control. With this, programmers can begin to describe quantum kernels of execution via sequences of arbitrary pulse shapes. This opens new avenues of research and development with regards to smart quantum compilation routines that enable direct translation of higher-level digital assembly representations to these native pulse instructions. In this work, we present an extension to the XACC system-level quantum-classical software framework that directly enables this compilation lowering phase via user-specified quantum optimal control techniques. This extension enables the translation of digital quantum circuit representations to equivalent pulse sequences that are optimal with respect to the backend system dynamics. Our work is modular and extensible, enabling third party optimal control techniques and strategies in both C++ and Python. We demonstrate this extension with familiar gradient-based methods like gradient ascent pulse engineering (GRAPE), gradient optimization of analytic controls (GOAT), and Krotov's method. Our work serves as a foundational component of future quantum-classical compiler designs that lower high-level programmatic representations to low-level machine instructions.
MPI library that supports serializing (and check-pointing) its MPI Stages is a communication library that allows nodes and programs in a computer cluster to communicate via the Message-Passing Interface (MPI). The internal data objects of the library are serializable, i.e., the user can convert the MPI state or objects, such as groups or communicators, into byte streams and save them into a checkpoint file. Later users can load the checkpoint file and de-serialize previous objects into the library. Serialization/de-serialization is useful in MPI applications to recover from node failures.