Application Note: Coupled Simulation with the Xyce™ General External Interface
The Xyce General External Coupling API is documented.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “computer interfaces”
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.
Cycle accurate and cycle reproducible memory for an FPGA based hardware accelerator
A method, system and computer program product are disclosed for using a Field Programmable Gate Array (FPGA) to simulate operations of a device under test (DUT). The DUT includes a device memory having a number of input ports, and the FPGA is associated with a target memory having a second number of input ports, the second number being less than the first number. In one embodiment, a given set of inputs is applied to the device memory at a frequency Fd and in a defined cycle of time, and the given set of inputs is applied to the target memory at a frequency Ft. Ft is greater than Fd and cycle accuracy is maintained between the device memory and the target memory. In an embodiment, a cycle accurate model of the DUT memory is created by separating the DUT memory interface protocol from the target memory storage array.
97 MATHEMATICS AND COMPUTING↗
Discontinuous Galerkin based weighted condition number relaxation for dynamic interface mesh adaptation
Abstract not provided
97 MATHEMATICS AND COMPUTING↗
Ionic density distributions near the charged colloids: Spherical electric double layers
We have studied the structure of the spherical electric double layers on charged colloids by a density functional perturbation theory, which is based both on the modified fundamental-measure theory for the hard spheres and on the one-particle direct correlation functional (DCF) for the electronic residual contribution. The contribution of one-particle DCF has been approximated as the functional integration of the second-order correlation function of the ionic fluids in a bulk phase. The calculated result is in very good agreement with the computer simulations for the ionic density distributions and the zeta potentials over a wide range of macroion sizes and electrolyte concentrations, and compares with the results of Yu et al. [J. Chem. Phys. 120, 7223 (2004)] and modified Poisson-Boltzmann approximation [L. B. Bhuiyan and C. W. Outhwaite, Condens. Matter Phys. 8, 287 (2005)]. The present theory is able to provide interesting insights about the charge inversion phenomena occurring at the interface.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Kokkos Ecosystem [Brief]
In 2016/2017, the field of High-Performance Computing (HPC) entered a new era driven by fundamental physics challenges to produce ever more energy and cost-efficient processors. Since the convergence on the Message-Passing Interface (MPI) standard in the mid-1990s, application developers enjoyed a seemingly static view of the underlying machine — that of a distributed collection of homogeneous nodes executing in collaboration. However, after almost two decades of dominance, the sole use of MPI to derive parallelism acted as a limiter to improved future performance. While MPI is widely expected to continue to function as the basic mechanism for communication between compute nodes for the immediate future, additional parallelism is required on the computing node itself if high performance and efficiency goals are to be realized. When reviewing the architectures of the top HPC systems today, the change in paradigm is clear: the compute nodes of the leading machines in the world are either powered by many-core chips with a few dozen cores each, or use heterogeneous designs, where traditional CPUs marshal work to massively parallel compute accelerators which has as many as 200,000 processing threads in flight simultaneously. Complicating matters further for application developers, each processor vendor has its own preferred way of writing code for their architecture.The Kokkos EcoSystem was released by Sandia in 2017 to address this new era in HPC system design by providing a vendor independent performance portable programming system for scientific, engineering, and mathematical software applications written in the C++ programming language. Using Kokkos, application developers can be more productive because they will not have to create and maintain separate versions of their software for each architecture, nor will they have to be experts in each architecture's peculiar requirements. Instead, they will have a single method of programming for the diverse set of modern HPC architectures. While Kokkos started in 2011 as a programming model only, it soon became clear that complex applications needed more. It is also critical to have a portable mathematical functions and developers need tools to debug their applications, gain insight into the performance characteristics of their codes and tune algorithm performance parameters through automated processes. The Kokkos EcoSystem addresses those needs through its three main components: the Kokkos Core programming model, the Kokkos Kernels math library, and the Kokkos Tools project.
97 MATHEMATICS AND COMPUTING↗
Optimizing Communication in 2D Grid-Based MPI Applications at Exascale
The new reality of exascale computing faces many challenges in achieving optimal performance on large numbers of nodes. A key challenge is the efficient utilization of the message-passing interface (MPI), a critical component for process communication. This paper explores communication optimization strategies to harness the GPU-accelerated architectures of these supercomputers. We focus on MPI applications where processors form a two-dimensional process grid, a common arrangement in applications involving dense matrix operations. This configuration offers a unique opportunity to implement innovative strategies to improve performance and maintain effective load distribution. We study two applications— Dist-FW (Apsp:all-pair-shortest-path) and HPL-MxP (LU factorization with Mixed precision)—on two accelerated systems: Summit (IBM Power and NVIDIA V100) and Frontier (AMD EPYC and MI250X). These supercomputers are operated by the Oak Ridge Leadership Computing Facility (OLCF) and are currently ranked #1 and #5 on the Top500 list. We show how to scale up both applications to exascale levels and tackle the MPI challenges related to implementation, synchronization, and performance. We also compare the performance of several communication strategies at an unprecedented scale. Accurately predicting application performance becomes crucial for cost reduction as the computational scale grows. To address this, we suggest a hyperbolic model as a better alternative to the traditional one-sided asymptotic model for predicting future application performance at such large scales.
Lu, Hao↗
Structural and electrochemical evolution of alloy interfacial layers in anode-free solid-state batteries
“Anode-free” solid-state batteries feature high energy density since there is no anode active material. Although the beneficial effects of interfacial layers at the anode-solid electrolyte interface have been demonstrated, the mechanisms through which they influence lithium plating/stripping are unclear. Here, in this study, we reveal the evolution of 100-nm silver and gold interfacial layers during lithium plating/stripping using electrochemical methods, electron microscopy, X-ray microcomputed tomography, and modeling. The alloy layers improve Coulombic efficiency and resistance to short circuiting, even though the alloys form solute regions or particulates that detach from the current collector during plating. In situ electrochemical impedance spectroscopy shows that the alloys return to the interface and mitigate contact loss at the end of stripping, avoiding a critical vulnerability of anode-free cells. Contact retention is driven by even Li thickness that promotes spatially uniform stripping, as well as local alloy delithiation in response to current concentrations that homogenizes current and diminishes voiding.
25 ENERGY STORAGE↗
Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies
Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.
36 MATERIALS SCIENCE↗
IJS procedure for RELAP5 to TRACE input model conversion using SNAP
The TRAC/RELAP Advanced Computational Engine (TRACE) advanced, best-estimate reactor systems code developed by the U.S. Nuclear Regulatory Commission comes with a graphical user interface called Symbolic Nuclear Analysis Package (SNAP). Much of efforts have been done in the past to develop the RELAP5 input decks. The purpose of this study is to demonstrate the Institut 'Josef Stefan' (IJS) conversion procedure from RELAP5 to TRACE input model of BETHSY facility. The IJS conversion procedure consists of eleven steps and is based on the use of SNAP. For calculations of the selected BETHSY 6.2TC test the RELAP5/MOD3.3 Patch 4 and TRACE V5.0 Patch 1 were used. The selected BETHSY 6.2TC test was 15.24 cm equivalent diameter horizontal cold leg break in the reference pressurized water reactor without high pressure and low pressure safety injection. The application of the IJS procedure for conversion of BETHSY input model showed that it is important to perform the steps in proper sequence. The overall calculated results obtained with TRACE using the converted RELAP5 model were close to experimental data and comparable to RELAP5/MOD3.3 calculations. Therefore it can be concluded, that proposed IJS conversion procedure was successfully demonstrated on the BETHSY integral test facility input model. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
The Open High Throughput Computing Content Delivery Network
LHC experiments make extensive use of web proxy caches, especially for software distribution via the CernVM File System and for conditions data via the Frontier Distributed Database Caching system. Since many jobs read the same data, cache hit rates are high and hence most of the traffic flows efficiently over Local Area Networks. However, it is not always possible to have local web caches, particularly for opportunistic cases whereexperiments have little control over site services. The Open High Throughput Computing (HTC) Content Delivery Network (CDN), openhtc.io, aims to address this by using web proxy caches from a commercial CDN provider. Cloudflare provides a simple interface for registering DNS aliases of any web server and does reverse proxy web caching on those aliases. The openhtc.io domain is hosted on Cloudflare's free tier CDN which has no bandwidth limit and makes use of data centers throughout the world, so the average performance for clients is much improved compared to reading from CERN or a Tier 1. The load on WLCG servers is also significantly reduced. WLCG Web Proxy Auto Discovery is used to select local web caches when they are available and otherwise select openhtc.io caching. This paper describes the Open HTC CDN in detail and provides initial results from its use for LHC@Home and USCMS opportunistic computing.
97 MATHEMATICS AND COMPUTING↗
Fast Transient Simulations for Multi-Segment Transmission Lines with a Graphical Model
This paper studies a computationally efficient algebraic graph theory engine for simulating time-domain one-dimensional waves in a multi-segment transmission line, such as for reflectometry applications. Efficient simulation of time-domain signals in multi-segment transmission lines is challenging because the number of propagation paths (and therefore the number of operations) increases exponentially with each new interface. We address this challenge through the use of a frequency-domain, algebraic graphical model of wave propagation, which is then converted to the time domain via the Fourier transform. We use this model to achieve an exact, stable, and computationally efficient (O(NQ), where N is the number of segments and Q is the bandwidth) approach for studying one-dimensional wave propagation. Our approach requires the reflection and transmission coefficients for each interface and each segment's complex propagation constant. Furthermore, we compare our simulation results with known analytical solutions.
42 ENGINEERING↗
SUNDIALS Multiphysics+MPIManyVector Performance Testing
In this report we document performance test results on a SUNDIALS-based multiphysics demonstration application. We aim to assess the large-scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically: 1. SUNDIALS’ new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources. 2. ARKode’s new multirate integration module, MRIStep, allowing high-order accurate calculations that subcycle "fast" processes within "slow" ones. 3. SUNDIALS’ new flexible linear solver interfaces, that allow streamlined specification of problem-specific linear solvers. 4. SUNDIALS’ new N_Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI_Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS’ generic SUNNonlinearSolver interface and accelerator-enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).
97 MATHEMATICS AND COMPUTING↗
Asperity generation and its relationship to seismicity on a planar fault: A laboratory simulation
Earthquake faults, and all frictional surfaces, establish contact through asperities. A detailed knowledge of how asperities form will enable a better understanding of the manner in which they communicate during foreshock failure sequences that are observed, leading to the larger main shock. We present results of experiments where a pressure sensitive film was used to map, size and measure the magnitudes of the normal stresses at asperities along a seismogenic section of a laboratory simulated fault. We measured seismicity acoustically and foreshocks were found to be the result of localized asperity failure during the nucleation phase of gross fault rupture. Since surface roughness plays an important role in how asperities are formed, two Hurst exponents were measured to characterize a highly worn interface using roughness profiles: (i) long wavelength estimates (H ~ 0.45) and (ii) short wavelength estimates (H ~ 0.8-1.2). The short wavelength roughness estimates were computed at the scale of single asperity junction points. Macroscopically, the number of asperities and real contact area increased with additional application of normal force while the mean normal stress remained constant. The ratio of real to nominal contact area was low - ranging from 0.02 < A r /A 0 < 0.05-predicting that the asperities should be elastically independent of each other. Results from the pressure sensitive film showed that asperities were closely spaced and could not be treated as mechanically independent. Larger asperities carried both higher levels of average normal stress and higher levels of normal stress heterogeneity than smaller ones. Using linear stability theorem, the critical slip distance on foreshocking asperities was estimated to be d 0 ~ 0.65-3 μm. The critical slip distance d 0 was ~1.8-11.5 per cent of the premonitory slip needed to initiate gross fault rupture of the interface (20-40 μm) and the overall slip necessary to initiate gross fault rupture was on the order of the average asperity diameter (52 μm). In conclusion, foreshocks may be due to a change in the critical slip distance, at localized sections of the fault, caused by the two distinct roughness profiles measured at short and long length scales.
58 GEOSCIENCES↗
Systems and methods for performing wireless financial transactions
A secure computing module (SCM) is configured for connection with a host device. The SCM includes a processor for performing secure processing operations, a host interface for coupling the processor to the host device, and a memory connected to the processor wherein the processor logically isolates at least some of the memory from access by the host device. The SCM also includes a proximate-field wireless communicator connected to the processor to communicate with another SCM associated with another host device. The SCM generates a secure digital signature for a financial transaction package and communicates the package and the signature to the other SCM using the proximate-field wireless communicator. Financial transactions are performed from person to person using the secure digital signature of each person's SCM and possibly message encryption. The digital signatures and transaction details are communicated to appropriate financial organizations to authenticate the transaction parties and complete the transaction.
97 MATHEMATICS AND COMPUTING↗
Institute for Sustained Performance, Energy, and Resilience (SuPER)
The University of Tennessee (UTK) and University of Texas at El Paso (UTEP) partnership supported the three main thrusts of the SUPER project---performance, energy, and resilience. The UTK-UTEP effort thus helped advance the main goal of SUPER, which was to ensure that DOE's computational scientists can successfully exploit the emerging generation of high performance computing (HPC) systems. This goal is being met by providing application scientists with strategies and tools to productively maximize performance, conserve energy, and attain resilience. The primary vehicle through which UTK provided performance measurement support to SUPER and the larger HPC community is the Performance Application Programming Interface (PAPI). PAPI is an ongoing project that provides a consistent interface and methodology for collecting hardware performance information from various hardware and software components, including most major CPUs, GPUs and accelerators, interconnects, I/O systems, and power interfaces, as well as virtual cloud environments. The PAPI software is widely used for performance modeling of scientific and engineering applications---for example, the HOMME (High Order Methods Modeling Environment) climate code, and the GAMESS and NWChem computational chemistry codes---on DOE supercomputers. PAPI is widely deployed as middleware for use by higher-level profiling, tracing, and sampling tools (e.g., CrayPat, HPCToolkit, Scalasca, Score-P, TAU, Vampir, PerfExpert), making it the de facto standard for hardware counter analysis. PAPI has established itself as fundamental software infrastructure in every application domain (spanning academia, government, and industry), where improving performance can be mission critical. Ultimately, as more application scientists migrate their applications to HPC platforms, they will benefit from the extended capabilities this grant brought to PAPI to analyze and optimize performance in these environments, whether they use PAPI directly, or via third-party performance tools. Capabilities added to PAPI through this grant include support for new architectures such as the lastest GPU and Xeon Phi accelerators, and advanced power measurement and management features. Another important topic for the UTK team was providing support for a rich ecosystem of different fault management strategies in the context of parallel computing. Our long term efforts have been oriented toward proposing flexible strategies and providing building boxes that application developers can use to build the most efficient fault management technique for their application. These efforts span across the entire software spectrum, from theoretical models of existing strategies to easily assess their performance, to algorithmic modifications to take advantage of specific mathematical properties for data redundancy and to extensions to widely used programming paradigms to empower the application developers to deal with all types of faults. We have also continued our tight collaborations with users to help them adopt these technologies to ensure their application always deliver meaningful scientific data. Large supercomputer systems are becoming more and more power and energy constrained, and future systems and applications running on them will need to be optimized to run under power caps and/or minimize energy consumption. The UTEP team contributed to the SUPER energy thrust by developing power modeling methodologies and investigating power management strategies. Scalability modeling results showed that some applications can scale better with respect to an increasing power budget than with respect to only the number of processors. Power management, in particular shifting power to processors on the critical path of an application execution, can reduce perturbation due to system noise and other sources of runtime variability, which are growing problems on large-scale power-constrained computer systems.
42 ENGINEERING↗
Theory of the spatial transfer of interface-nucleated changes of dynamical constraints and its consequences in glass-forming films
We formulate a new theory for how caging constraints in glass-forming liquids at a surface or interface are modified and then spatially transferred, in a layer-by-layer bootstrapped manner, into the film interior in the context of the dynamic free energy concept of the Nonlinear Langevin Equation (NLE) theory approach. The dynamic free energy at any mean location (cage center) involves contributions from two adjacent layers where confining forces are not the same. At the most fundamental level of the theory, the caging component of the dynamic free energy varies essentially exponentially with distance from the interface, saturating deep enough into the film with a correlation length of modest size and weak sensitivity to the thermodynamic state. This imparts a roughly exponential spatial variation of all the key features of the dynamic free energy required to compute gradients of dynamical quantities including the localization length, jump distance, cage barrier, collective elastic barrier, and alpha relaxation time. The spatial gradients are entirely of dynamical, not structural or thermodynamic, origin. The theory is implemented for the hard sphere fluid and diverse interfaces which can be a vapor, a rough pinned particle solid, a vibrating (softened) pinned particle solid, or a smooth hard wall. Their basic description at the level of the spatially heterogeneous dynamic free energy is identical, with the crucial difference arising from the first layer where dynamical constraints can be weakened, softened, or hardly changed depending on the specific interface. Numerical calculations establish the spatial dependence and fluid volume fraction sensitivity of the key dynamical property gradients for five different model interfaces. A comparison of the theoretical predictions for the dynamic localization length and glassy modulus with simulations and experiments for systems with a vapor interface reveals good agreement. Finally, the present advance sets the stage for using the Elastically Collective NLE theory to make quantitative predictions for the alpha relaxation time gradient, decoupling phenomena, T g gradient, and many film-averaged properties of both model and experimental (colloids, molecules, and polymers) systems with diverse interfaces and chemical makeup.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
THEMIS: A Parameter Estimation Framework for the Event Horizon Telescope
The Event Horizon Telescope (EHT) provides the unprecedented ability to directly resolve the structure and dynamics of black hole emission regions on scales smaller than their horizons. This has the potential to critically probe the mechanisms by which black holes accrete and launch outflows, and the structure of supermassive black hole spacetimes. However, accessing this information is a formidable analysis challenge for two reasons. First, the EHT natively produces a variety of data types that encode information about the image structure in nontrivial ways; these are subject to a variety of systematic effects associated with very long baseline interferometry and are supplemented by a wide variety of auxiliary data on the primary EHT targets from decades of other observations. Second, models of the emission regions and their interaction with the black hole are complex, highly uncertain, and computationally expensive to construct. As a result, the scientific utilization of EHT observations requires a flexible, extensible, and powerful analysis framework. We present such a framework, Themis, which defines a set of interfaces between models, data, and sampling algorithms that facilitates future development. We describe the design and currently existing components of Themis, how Themis has been validated thus far, and present additional analyses made possible by Themis that illustrate its capabilities. Importantly, we demonstrate that Themis is able to reproduce prior EHT analyses, extend these, and do so in a computationally efficient manner that can efficiently exploit modern high-performance computing facilities. Themis has already been used extensively in the scientific analysis and interpretation of the first EHT observations of M87.
79 ASTRONOMY AND ASTROPHYSICS↗
The LLNL nuclear data infrastructure for the GNDS data format
The next generation of nuclear data infrastructure tools at the Livermore National Laboratory (LLNL) consists of pipeline of codes that read and process nuclear data from evaluated files saved in the new GNDS (Generalised Nuclear Data Structure) nuclear data format. The processing code FUDGE (For Updating Data and Generating Evaluations) is at the front-end of this pipeline as it reads and process the evaluated data for use in downstream transport codes. FUDGE is Python based with C and C++ extensions for computationally intensive tasks. As is the case for the evaluated data, the processed output is also saved in the GNDS format and the GIDI+ API is provided as the interface between the processed data and the transport codes. GIDI+ is a C++ based suite of codes and it includes GIDI (General Interaction Data Interface), a library for reading and writing GNDS data, and MCGIDI which is the cross section lookup, and reaction and product distribution sampling interface between Monte Carlo transport codes and the GNDS data. GIDI provides methods for easy access to the multi-group processed GNDS data and this is demonstrated through its implementation in ARDRA, the LLNL deterministic transport code. The evaluation and sampling methods in MCGIDI are available as both CPU and GPU methods which facilitates the use of MCGIDI in both traditional CPU-based as well as the next generation mixed model computational architectures. This is demonstrated through the GIDI+ implementation in MERCURY, the LLNL Monte Carlo transport code. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗