Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Asymptotic Behavior of Memristive Circuits

The interest in memristors has risen due to their possible application both as memory units and as computational devices in combination with CMOS. This is in part due to their nonlinear dynamics, and a strong dependence on the circuit topology. We provide evidence that also purely memristive circuits can be employed for computational purposes. In the present paper we show that a polynomial Lyapunov function in the memory parameters exists for the case of DC controlled memristors. Such a Lyapunov function can be asymptotically approximated with binary variables, and mapped to quadratic combinatorial optimization problems. This also shows a direct parallel between memristive circuits and the Hopfield-Little model. In the case of Erdos-Renyi random circuits, we show numerically that the distribution of the matrix elements of the projectors can be roughly approximated with a Gaussian distribution, and that it scales with the inverse square root of the number of elements. This provides an approximated but direct connection with the physics of disordered system and, in particular, of mean field spin glasses. Using this and the fact that the interaction is controlled by a projector operator on the loop space of the circuit. We estimate the number of stationary points of the approximate Lyapunov function and provide a scaling formula as an upper bound in terms of the circuit topology only.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

MODELING THE AMBIENT CONDITION EFFECTS OF AN AIR-COOLED NATURAL CIRCULATION SYSTEM

The Reactor Cavity Cooling System (RCCS) is a passive safety concept under consideration for the overall safety strategy of advanced reactors such as the High Temperature Gas-Cooled Reactor (HTGR). One such variant, air-cooled RCCS, uses natural convection to drive the flow of air from outside the reactor building to remove decay heat during normal operation and accident scenarios. The Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne National Laboratory (“Argonne”) is a half-scale model of the primary features of one conceptual air-cooled RCCS design. The facility was constructed to carry out highly instrumented experiments to study the performance of the RCCS concept for reactor decay heat removal that relies on natural convection cooling. Parallel modeling and simulation efforts were performed to support the design, operation, and analysis of the natural convection system. Throughout the testing program, strong influences of ambient conditions were observed in the experimental data when baseline tests were repeated under the same test procedures. Thus, significant analysis efforts were devoted to gaining a better understanding of these influences and the subsequent response of the NSTF to ambient conditions. It was determined that air humidity had negligible impacts on NSTF system performance and therefore did not warrant consideration in the models. However, temperature differences between the building exterior and interior air, along with the outside wind speed, were shown to be dominant factors. Combining the stack and wind effects together, an empirical model was developed based on theoretical considerations and using experimental data to correlate zero-power system flow rates with ambient meteorological conditions. Some coefficients in the model were obtained based on best fitting the experimental data. The predictive capability of the empirical model was demonstrated by applying it to the new set of experimental data. The empirical model was also implemented in the computational models of the NSTF using both RELAP5-3D and STARCCM+ codes. Accounting for the effects of ambient conditions, simulations from both codes predicted the natural circulation flow rates very well.

Hu, Rui↗
DOE OSTI2024

Thermodynamic Limits of Redox-Based Thermochemical Processes (REDOTHERM)

Solar thermochemical fuel production is a potential pathway for the production of sustain liquid drop-in fuels, which can help decarbonize the aviation and maritime sectors. In an attempt to analyze the commercial viability of this technology, several studies have been conducted, including system and technoeconomic analysis (TEA) modeling. However, most studies to date simply assume a given redox reactor efficiency, which is significantly higher than demonstrated values to date. While it is widely recognized that utilizing a counter-current flow (CF) configuration could increase the redox reactor efficiency, an over-simplification in the thermodynamic modeling may lead to unphysical results which has been included in multiple publications. The fact that the solar redox reactor is the least developed component in the process chain makes it hard to identify technology gaps and evaluate pathways to deployment at scale using this approach. In this work, a thermodynamic model for a moving oxide system has been developed, in a general form that allows to analyze the system for different redox-active materials, under a wide range of operating conditions, for both parallel and countercurrent flows. The model capabilites are demonstrated, and the model's code will be shared as an open-source on GitHub in the next few months.

chemical looping↗
DOE OSTI2022

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE [Slides]

Depending on the system, energy, and/or nuclide, altering the number of latent generations can significantly impact generated sensitivity coefficients. In general, as the number of latent generations increases, the accuracy of the generated sensitivity value compared with the DP value increases and there is a significant increase in the uncertainty with generated sensitivity values. IFP generated sensitivities appear to be more stable than those with CLUTCH. Complex systems, even with IFP and increased latent generations can still produce poor results (i.e., MPC- 32). This reiterates the importance of performing DPs to confirm sensitivities. IFP-Shift and CLUTCH can take advantage of parallel computing abilities. Work continues in this area as additional parameters and calculational methods are examined to provide analysts insights to successfully generating sensitivity coefficients for confirmatory analyses and validation efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Scalable semitransparent prototype organic photovoltaic module with minimal resistance loss

Semi-transparent organic photovoltaic (ST-OPV) cells are regarded as an attractive solution to building-integrated solar energy harvesting. Both the power conversion efficiency (PCE) and average photopic transmission (APT) of ST-OPVs have shown substantial increases in recent years. However, less attention has been paid to the area scaling of ST-OPV cells. In this work, we investigate the scalability of ST-OPV cells from 4 mm 2 to 1 cm 2 , as well as 9 cm 2 active area prototype module. By integrating both top and bottom metal grids onto the transparent electrodes, the series resistance loss of 1 cm 2 ST-OPV cell is significantly reduced, and is comparable to that of the grid-free 4 mm 2 cell. Nine 1 cm 2 cells are then connected in a series-parallel circuit to realize a prototype ST-OPV module. A 100% fabrication yield with only 5% PCE deviation among discrete cells is achieved. The semitransparent module shows PCE = 7.2 ± 0.1% under simulated AM 1.5G illumination at 1 sun intensity, which exhibits no connection resistance loss compared to that of the individual cells. The ST-OPV module exhibits an APT = 38.1 ± 1.1%, which enables a light utilization efficiency, LUE = 2.74 ± 0.09%. Here, the method demonstrates a promising way for ST-OPV modules to scale without compromising performance.

14 SOLAR ENERGY↗
DOE OSTI2022

Developing an ELM Ecosystem Dynamics Model on GPU with OpenACC

Porting a complex scientific code, such as the E3SM land model (ELM), onto a new computing architecture is challenging. The paper presents design strategies and technical approaches to develop an ELM ecosystem dynamics model with compiler directives (OpenACC) on NVIDIA GPUs. The code has been refactored with advanced OpenACC features (such as deepcopy and routine directives) to reduce memory consumption and to increase the levels of parallelism through parallel loop reconstruction and new data structures. As a result, the optimized parallel implementation achieved more than a 140-time speedup (50 ms vs 7600 ms), compared to a naive implementation that uses OpenACC routine directive and parallelizes the code across existing loops on a single NVIDIA V100. On a fully loaded computing node with 44 CPUs and 6 GPUs, the code achieved over a 3.0-times speedup, compared to the original code on the CPU. Furthermore, the memory footprint of the optimized parallel implementation is 300 MB, which is around 15% of the 2.15 GB of memory consumed by a naive implementation. This study is the first effort to develop the ELM component on GPUs efficiently to support ultra-high-resolution land simulations at continental scales.

Schwartz, Peter↗
DOE OSTI2023

Multipoint Observations of Dispersive Scale Alfvénic Field‐Line Resonances Associated With Substorm Auroral Beads

We present a case study of the field-aligned current (FAC) systems that transpire within the high-altitude auroral acceleration region of an “auroral bead” initiated double oval substorm observed on 23 February 2001 by the Cluster fleet. Conjunctive Cluster measurements and auroral images from IMAGE reveal that auroral bead current system formation and evolution is a multi-scale, injection-mediated process. The FACs at large scales vary on substorm evolution time scales (~minutes) in response to the injection and evolution of hotter denser magnetospheric plasma. Embedded within the large-scale FACs are intense short-scale (≲ few 10s of km) currents comprising dispersive scale Alfvén wave (DAW) fluctuations. The DAWs are a complex mixture of ingoing and reflected components that regularly interfere to form a broad spectrum of kinetic (dispersive) scale Alfvénic field-line resonances (KFLRs). The Alfvénic currents appear as a nested series of upward and downward FAC densities with amplitudes reaching a few 100 nA/m 2 . Energized field-aligned or counterstreaming electrons near keV energies and below are observed with parallel skews that vary in concert with variations in the DAW current sense. Positive correlations between DAW electric field energy densities and the energies of energized H + , He + , and O + outflow are observed, indicative of ion energization within the DAW fields. Due to their L-shell location (L~5.8–7.0) and associations with injections, the KFLRs are interpreted as the high-altitude auroral zone analog of KFLRs observed in the equatorial inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Towards Highly Scalable Ab Initio Molecular Dynamics (AIMD) Simulations on the Intel Knights Landing Manycore Processor

The Ab Initio Molecular Dynamics (AIMD) method allows scientists to treat the dynamics of molecular and condensed phase systems while retaining a first-principles-based description of their interactions. This extremely important method has tremendous computational requirements, because the electronic Schr¨odinger equation, approximated using Kohn-Sham Density Functional Theory (DFT), is solved at every time step. With the advent of manycore architectures, application developers have a significant amount of processing power within each compute node that can only be exploited through massive parallelism. A compute intensive application such as AIMD forms a good candidate to leverage this processing power. In this paper, we focus on adding thread level parallelism to the plane wave DFT methodology implemented in NWChem. Through a careful optimization of tall-skinny matrix products, which are at the heart of the Lagrange multiplier and nonlocal pseudopotential kernels, as well as 3D FFTs, our OpenMP implementation delivers excellent strong scaling on the latest Intel Knights Landing (KNL) processor. We assess the efficiency of our Lagrange multiplier kernels by building a Roofline model of the platform, and verify that our implementation is close to the roofline for various problem sizes. Finally, we present strong scaling results on the complete AIMD simulation for a 64 water molecules test case, that scales up to all 68 cores of the Knights Landing processor.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Roll-to-Roll SOFC Manufacturing for Low Cost Energy Generation

Oak Ridge National Laboratory (ORNL) and Redox Power Systems, LLC. collaborated on this project to develop an integrated, continuous roll-to-toll (R2R) manufacturing processes to produce solid oxide fuel cells (SOFCs). SOFCs are ceramic-based electrochemical devices that efficiently produce electricity from hydrogen-containing fuels, including natural gas. Standard processes for manufacturing SOFCs include tape casting of the anode and electrolyte, followed by stacking of the layers, and lamination using a uniaxial press or isostatic press. After high temperature firing, the cathode and electrical contact are screen printed and fired once again to complete the cell. This effort aimed to further develop R2R tape fabrication processes for green cell fabrication with the DOE Battery Manufacturing Facility (BMF) at ORNL as next generation improvements to current methods, focusing on R2R calendering and slot die coating for lower cost, high volume approaches for multilayer green cell fabrication. Two paths were investigated in parallel in this project. The first path was to optimize R2R processing using traditional tape cast materials from the Redox production chain on the BMF’s heated calender for compressing and laminating multilayer materials. The second path was to use BMF’s pilot-scale slot-die coater for R2R manufacturing as an alternative to tape casting of some or all of the layers in the green cell. Processes have been developed and and process conditions have been identified. Additionally, basic process economics were estimated for comparison to Redox current manufacturing approaches for the applicable process modules. The SBV program results will enable Redox to advance lower cost R2R manufacturing processes to produce millions of green cells per year on a single production line. As a result of this work, one invention disclosure has been submitted. In addition, one paper is under preparation for submission to a peer-reviewed journal.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A Platform Technology for High-throughput Atomically Precise Manufacturing: Mechatronics at the Atomic Scale

The main objective of this project is to invent the necessary enabling technologies for high throughput atomically precise manufacturing (APM). APM is an emerging technology that refers to any manufacturing capability that enables fabrication of atomically precise structures, components, and devices under programmable control. APM will require positional assembly at the atomic and/or molecular scale, as well as at the nano and microscales using hierarchical assembly to create products ranging from nanoscale and quantum devices to macroscale systems and materials. This project builds on the expectation that commercial viability of APM will depend on a high level of parallelism to achieve the required throughput, a capability that currently does not exist. This project is the first attempt to address this key technological bottleneck.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Electron Acceleration at Pulsar Wind Termination Shocks

We study the acceleration of electrons and positrons at an electromagnetically modified, ultrarelativistic shock in the context of pulsar wind nebulae. We simulate the outflow produced by an obliquely rotating pulsar in proximity of its termination shock with a two-fluid code that uses a magnetic shear wave to mimic the properties of the wind. We integrate electron trajectories in the test-particle limit in the resulting background electromagnetic fields to analyze the injection mechanism. We find that the shock-precursor structure energizes and reflects a sizable fraction of particles, which becomes available for further acceleration. We investigate the subsequent first-order Fermi process sustained by small-scale magnetic fluctuations with a Monte Carlo code. We find that the acceleration proceeds in two distinct regimes: when the gyroradius r{sub g} exceeds the wavelength of the shear λ, the process is remarkably similar to first-order Fermi acceleration at relativistic, parallel shocks. This regime corresponds to a low-density wind that allows the propagation of superluminal waves. When r{sub g}<λ, which corresponds to the scenario of driven reconnection, the spectrum is softer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modeling UV Radiation Feedback from Massive Stars. I. Implementation of Adaptive Ray-tracing Method and Tests [Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Algorithm]

Here, we present an implementation of an adaptive ray-tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses nonblocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems, including the propagation of radiation in vacuum and the expansions of various types of H II regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim {10}^{3}$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is 12% when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with those from the M 1 closure relation. Although the ART and M 1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Myth of Fungible FTE: A Quantitative Assessment of Matrixed Resource Allocation

Matrix organizations allow scientific facilities to share specialized personnel across projects, operations, maintenance, and strategic initiatives. Nominal staffing allocations, however, may not capture the schedule consequences of fragmented individual commitments, limited access to specialist groups, and intermittent availability of key decision makers. We developed a stochas- tic, daily-time-step simulation of a hypothetical medium-sized accelerator-facility project com- prising sequential phases and parallel tasks. Each task requires role-specific work measured in FTE-days. Ordinary personnel may be unavailable because they contribute concurrently to other institutional activities, while designated key roles have independently specified daily un- availability probabilities. An organization-wide priority factor scales the number of people from each functional group who can effectively contribute to the project. It is interpreted as a composite proxy for project access and workforce fragmentation across competing commit- ments. We examined project completion time as a function of this factor and Project Lead unavailability using 100 Monte Carlo runs per condition. Increasing priority factor from 0.1 to 1.0 reduced median completion time from 1708.5 days (interquartile range 1681.5–1735.25) to 390 days (interquartile range 379–399). At priority factor = 0.1, increasing Project Lead unavailability from 0.5 to 0.9 increased median completion time from 1713.5 days (interquartile range 1691–1733.25) to 4,417 days (interquartile range 4271.75–4550.5). The model quantifies the commonly expected sensitivity of project schedules to fragmented resource commitments and limited coordination availability. Within this model, the results also indicate a possible threshold regime in which small increases in workforce availability yield only modest sched- ule improvements until sufficient capacity becomes accessible, after which project performance improves sharply. With further validation and calibration, this quantitative framework could support resource-allocation decisions during initial project planning and subsequent schedule rebaselining.

Bai, Mei [SLAC National Accelerator Laboratory (SL↗
DOE OSTI2016

Parallel electric fields are inefficient drivers of energetic electrons in magnetic reconnection

We present two-dimensional kinetic simulations, with a broad range of initial guide fields, which isolate the role of parallel electric fields (E- ∥ ) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, E- ∥ drives essentially all of the electron energy gains, yet fails to generate an energetic component. Here, we suggest that this is due to the weak energy scaling of particle acceleration from E- ∥ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. Lastly, this result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗
DOE OSTI2017

AMPS: An Augmented Matrix Formulation for Principal Submatrix Updates with Application to Power Grids

We present an augmented matrix approach to update the solution to a linear system of equations when the coefficient matrix is modified by a few elements within a principal submatrix. This problem arises in the dynamic security analysis of a power grid, where operators need to perform $N-x$ contingency analysis, i.e., determine the state of the system when up to $x$ links from $N$ fail. Our algorithms augment the coefficient matrix to account for the changes in it, and then compute the solution to the augmented system without refactoring the modified matrix. We provide two algorithms, a direct method, and a hybrid direct-iterative method for solving the augmented system. We also exploit the sparsity of the matrices and vectors to accelerate the overall computation. Our algorithms are compared on three power grids with PARDISO, a parallel direct solver, and CHOLMOD, a direct solver with the ability to modify the Cholesky factors of the coefficient matrix. We show that our augmented algorithms outperform PARDISO (by two orders of magnitude), and CHOLMOD (by a factor of up to 5). Further, our algorithms scale better than CHOLMOD as the number of elements updated increases. The solutions are computed with high accuracy. Our algorithms are capable of computing $N-x$ contingency analysis on a $778K$ bus grid, updating a solution with $x=20$ elements in $1.6 \times 10^{-2}$ seconds on an Intel Xeon processor.

Yeung, Yu-Hong↗
DOE OSTI2016

FY 2016 Fusion Energy Sciences (FES) Theory and Simulation Performance Target

Annual target: Predicting the magnitude and scaling of the divertor heat load width in magnetically confined burning plasmas is a high priority for the fusion program and ITER. One of the key unresolved physics issues is what sets the heat flux width at the entrance to the divertor region. Perform massively parallel simulations using 3D edge kinetic and fluid codes to determine the parameter dependence of the heat load width at the divertor entrance and compute the divertor plate heat flux applicable to moderate particle recycling conditions. Comparisons will be made with data from DIII-D, NSTX-U, and C-Mod.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗