Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Strong, ductile, and hierarchical hetero-lamellar-structured alloys through microstructural inheritance and refinement

The strength−ductility trade-off exists ubiquitously, especially in brittle intermetallic-containing multiple principal element alloys (MPEAs), where the intermetallic phases often induce premature failure leading to severe ductility reduction. Hierarchical heterogeneities represent a promising microstructural solution to achieve simultaneous strength−ductility enhancement. However, it remains fundamentally challenging to tailor hierarchical heterostructures using conventional methods, which often rely on costly and time-consuming processing. Here, we report a multiscale microstructural inheritance and refinement strategy to process “structural hierarchy precursors” in as-cast heterogeneous Al 0.7 CoCrFeNi MPEAs, which lead directly to a hierarchical hetero-lamellar structure (HLS) after simple rolling and annealing. Interestingly, it takes only 10 min of annealing time, two orders of magnitude less than that required to render the state-of-the-art properties during conventional processing of Al 0.7 CoCrFeNi, for us to achieve record-high strength−ductility combinations via the hierarchical HLS design that sequentially stimulates multiple unusual deformation and reinforcement mechanisms. In particular, the HLS-enabled high hetero-deformation-induced (HDI) internal stress triggers profuse <111>-type dislocations on over five independent slip systems in the supposedly brittle intermetallic phase and activates extensive stacking faults (SFs) and nanotwinning in the adjoining soft phase with a rather high SF energy. These unexpected, dynamically reinforcing hetero-deformation mechanisms across multiple length scales facilitate high sustained HDI strain hardening, along with a salient microcrack-mediated extrinsic ductilization effect, suggesting that the proposed microstructural inheritance and refinement strategy provides an efficient, fast, and low-cost approach to overcome the strength−ductility trade-off in a broad range of structural materials.

Science & Technology - Other Topics↗
DOE OSTI2024

A molecular view of peptoid-induced acceleration of calcite growth

The extensive deposits of calcium carbonate (CaCO 3 ) generated by marine organisms constitute the largest and oldest carbon dioxide (CO 2 ) reservoir. These organisms utilize macromolecules like peptides and proteins to facilitate the nucleation and growth of carbonate minerals, serving as an effective method for CO 2 sequestration. However, the precise mechanisms behind this process remain elusive. In this study, we report the use of sequence-defined peptoids, a class of peptidomimetics, to achieve the accelerated calcite step growth kinetics with the molecular level mechanistic understanding. By designing peptoids with hydrophilic and hydrophobic blocks, we systematically investigated the acceleration in step growth rate of calcite crystals using in situ atomic force microscopy (AFM), varying peptoid sequences and concentrations, CaCO 3 supersaturations, and the ratio of Ca 2+ / HCO 3 − . Mechanistic studies using NMR, three-dimensional fast force mapping (3D FFM), and isothermal titration calorimetry (ITC) were conducted to reveal the interactions of peptoids with Ca 2+ and HCO 3 − ions in solution, as well as the effect of peptoids on solvation and energetics of calcite crystal surface. Our results indicate the multiple roles of peptoid in facilitating HCO 3 − deprotonation, Ca 2+ desolvation, and the disruption of interfacial hydration layers of the calcite surface, which collectively contribute to a peptoid-induced acceleration of calcite growth. These findings provide guidelines for future design of sequence-specific biomimetic polymers as crystallization promoters, offering potential applications in environmental remediation (such as CO 2 sequestration), biomedical engineering, and energy storage where fast crystallization is preferred.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2026

L-VISP: LSTM Visualization for Interpretable Symptom Prediction in Patient Cohorts

Symptom modelling in head and neck cancer is challenged by the complexity of heterogeneous patient data, leading to an interest in deep learning approaches. Although Long Short-Term Memory Networks (LSTMs) have shown great results in patient risk prediction, their low interpretability requires data modellers to collaborate with clinical experts to validate the results. We present L-VISP, a human–machine solution that uses visual analytics for LSTM modelling in clinical research. L-VISP uses custom visual encodings to make multiple LSTM variants interpretable, supporting a full range of analysis, from understanding model operations and evaluating performance to interpreting results in a clinical context. We evaluate L-VISP with data modellers and a clinical oncologist and present the takeaways from this multidisciplinary collaboration.

LSTM modeling↗
DOE OSTI2021

Critical and Californium Source-Driven Noise Analysis Subcritical Measurements with an Unreflected Cylindrical Tank of Mixed Uranium-Plutonium Nitrate Solution

The reported experiments were performed by Oak Ridge National Laboratory (ORNL) at the Battelle Northwest Laboratory’s (now Pacific Northwest National Laboratory) critical experiments facility at Hanford, Washington in 1981 and used 16 days of critical facility time, not including 10 days for setup and removal of ORNL equipment. These measurements were to assess the capability of the Cf source-driven noise analysis (CSDNA) method to measure the subcriticality (keff) of mixed U-Pu nitrate solutions. In addition to the CSNDA measurements, measurements were also performed near delayed criticality where CSDNA measurements cannot be performed. This report documents the experiments that were not reported at that time by presenting the ORNL experimental results and any online analysis performed during and shortly after the measurements. The mixed nitrate solution had a U concentration at 188 grams per liter (g/L), a Pu concentration of 280 g/L, free acid normality of 2.80, H ion molarity of 5.9 and a specific gravity of 1.754 g/cm 3 , a 240 Pu isotopic content of 7.981 wt. %, and a 235 U isotopic content of 0.724 wt. %. The stainless-steel tank for the solution had an inside diameter of 35.38 cm, an outside diameter of 35.53 cm, a height of 56.72 cm, and bottom thickness of 0.9525 cm. A Zircaloy pipe with a 3.1496 cm outside diameter, a 2.7788 cm inside diameter, and bottom thickness of 0.635 cm was available for insertion of the Cf source in the center of the fissile solution. The Cf source was also located at the outside surface of the tank (solution height varied from 10 to 53 cm) and in the center of the solution (solution height varied from 10 to 60.7 cm). The CSDNA measurements were not analyzed online to determine the subcritical neutron multiplication factors. At all subcritical states, the break frequency noise analysis data was fitted to obtain the prompt neutron decay constant. The neutron multiplication factors were determined for the two configurations of the measurements near delayed criticality. The subcritical neutron multiplication factors from the CNSDA measurements can be obtained with further analysis. However, the near delayed critical configuration, the prompt neutron decay constants, the count rates, and the measured cross and auto power spectral densities can be calculated directly for benchmarking. Much of data presented in this report are from ORNL notes—not in the ORNL logbooks. For the final benchmark analysis, the data from the Battelle Northwest Laboratory (which operated the critical facility in 1981) critical facility logbook should be consulted and be incorporated where appropriate. The purpose of this report is to document the experimental information for the measurements performed so that at a later date researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Criticality Safety Benchmark Program (ICSBEP) or Nuclear Energy Agency benchmark. The data from these measurements are available from the ORNL Records Management Services Department, and the logbook is available from ICSBEP at Idaho National Laboratory. Preparation of the present report is part of a larger cooperative effort between Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) to document more than 15 undocumented critical and subcritical experiments enumerated in ORNL/TM-2019/18 and performed by ORNL at ORCEF and other USDOE critical experiments facilities using more than 500 operational days of critical facility time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

As the growing availability of computational power slows, there has been an increasing reliance on algorithmic advances. However, faster algorithms alone will not necessarily bridge the gap in allowing computational scientists to study problems at the edge of scientific discovery in the next several decades. Often, it is necessary to simplify or precondition solvers to accelerate the study of large systems of linear equations commonly seen in a number of scientific fields. Preconditioning a problem to increase efficiency is often seen as the best approach; yet, preconditioners which are fast, smart, and efficient do not always exist. Following the progress of [1], we present a new preconditioner for symmetric diagonally dominant (SDD) systems of linear equations. These systems are common in certain PDEs, network science, and supervised learning among others. Based on spectral support graph theory, this new preconditioner builds off of the work of [2], computing and applying a V-cycle chain of approximate inverse matrices. This preconditioner approach is both algebraic in nature as well as hierarchically-constrained depending on the condition number of the system to be solved. Due to its generation of an Approximate Inverse Chain of matrices, we refer to this as the AIC preconditioner. We further accelerate the AIC preconditioner by utilizing precomputations to simplify setup and multiplications in the con-text of an iterative Krylov-subspace solver. While these iterative solvers can greatly reduce solution time, the number of iterations can grow large quickly in the absence of good preconditioners. Initial results for the AIC preconditioner have shown a very large reduction in iteration counts for SDD systems as compared to standard preconditioners such as Incomplete Cholesky (ICC) and Multigrid (MG). We further show significant reduction in iteration counts against the more advanced Combinatorial Multigrid (CMG) preconditioner. We have further developed no-fill sparsification techniques to ensure that the computational cost of applying the AIC preconditioner does not grow prohibitively large as the depth of the V-cycle grows for systems with larger condition numbers. Our numerical results have shown that these sparsifiers maintain the sparsity structure of our system while also displaying significant reductions in iteration counts.1 2

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

As the growing availability of computational power slows, there has been an increasing reliance on algorithmic advances. However, faster algorithms alone will not necessarily bridge the gap in allowing computational scientists to study problems at the edge of scientific discovery in the next several decades. Often, it is necessary to simplify or precondition solvers to accelerate the study of large systems of linear equations commonly seen in a number of scientific fields. Preconditioning a problem to increase efficiency is often seen as the best approach; yet, preconditioners which are fast, smart, and efficient do not always exist. Following the progress of [1], we present a new preconditioner for symmetric diagonally dominant (SDD) systems of linear equations. These systems are common in certain PDEs, network science, and supervised learning among others. Based on spectral support graph theory, this new preconditioner builds off of the work of [2], computing and applying a V-cycle chain of approximate inverse matrices. This preconditioner approach is both algebraic in nature as well as hierarchically-constrained depending on the condition number of the system to be solved. Due to its generation of an Approximate Inverse Chain of matrices, we refer to this as the AIC preconditioner. We further accelerate the AIC preconditioner by utilizing precomputations to simplify setup and multiplications in the con-text of an iterative Krylov-subspace solver. While these iterative solvers can greatly reduce solution time, the number of iterations can grow large quickly in the absence of good preconditioners. Initial results for the AIC preconditioner have shown a very large reduction in iteration counts for SDD systems as compared to standard preconditioners such as Incomplete Cholesky (ICC) and Multigrid (MG). We further show significant reduction in iteration counts against the more advanced Combinatorial Multigrid (CMG) preconditioner. We have further developed no-fill sparsification techniques to ensure that the computational cost of applying the AIC preconditioner does not grow prohibitively large as the depth of the V-cycle grows for systems with larger condition numbers. Our numerical results have shown that these sparsifiers maintain the sparsity structure of our system while also displaying significant reductions in iteration counts.1 2

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Role of Radiation-Induced Non-Equilibrium Plutonium Oxidation States in Solution

Plutonium plays a key role in global actinide research and nuclear fuel cycle technologies, and yet, our fundamental understanding of its inherent radiation-induced chemical behavior is limited. These radiation-induced processes cannot simply be switched off, as they are as fundamentally inherent to plutonium as the impact of relativistic effects on its f-electrons. In less chemically complex actinide systems, such as aqueous solutions of neptunium and americium, both steady-state and transient radiolysis products play a significant role in the redox cycling of their oxidation states, which ultimately impacts their chemistry and transport. However, plutonium's multiple, coexisting, and chemically active oxidation states, which comprise of bare ions and dioxo cations, provide additional redox pathways that complicate the aforementioned radiation-induced competition kinetics. Mechanistically understanding the response of these multiple oxidation states, both equilibrium and non-equilibrium, to intense ionizing radiation fields is essential for predicting the behavior of plutonium under a host of conditions that support technological innovation in global nuclear energy and non-proliferation efforts. Here, advances in radiation-induced plutonium solution chemistry will be presented, including new results from time-resolved electron pulse and steady-state gamma and alpha irradiations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Projected Multi-Agent Consensus Equilibrium (PMACE) With Application to Ptychography

Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this article, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications. We describe PMACE in terms of an equilibrium problem and an equivalent fixed-point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. In conclusion, our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Greenhouse gas emissions reduction strategies that maximize portfolio-wide life cycle cost reduction, resilience, and environmental justice benefits

While strategies to achieve net-zero emissions at an individual site are well understood, new analysis methods are required for organizations seeking to achieve net-zero across multiple facilities, each with concurrent priority goals. At a portfolio level, distinct locations present varied challenges that cannot be addressed through singular solutions, and competing goals can take precedence with the assumption that net-zero emissions strategies deter from energy resilience and cost savings, therefore negatively impacting nearby communities. This study tests these assumptions by analyzing 16 diverse sites (varying in size, climate, and energy use) to identify strategies that reduce emissions and assess the impact these strategies have on life cycle costs, resilience, and communities with environmental justice concerns. Methods were developed to approximate missing information essential to net-zero evaluation. Established methods were augmented to evaluate life cycle costs, resilience, and environmental justice impacts across a set of strategies and accommodate the multi-criteria analyses. Potential benefits from identified strategies were quantified using site characteristics and a set of corresponding metrics. The net-zero analysis found that 11 sites could use on-site strategies to eliminate all but 2% of emissions generated. The remaining emissions can be offset, for instance through sequestration, executed at the portfolio scale. On-site carbon-free energy was found to reduce 51% of emissions across all sites; efficiency reduced 19% of emissions; sequestration 16%; procured carbon-free energy 15%; fuel switching 1.6%; and fleet electrification 1.3%. Building electrification, however, increased emissions by 4.4%. Different strategies also provide cost, resilience, and/or environmental justice benefits—the degree to which varies with individual site conditions. The findings indicate an advantage to considering the strategies as a comprehensive set, which leads to co-benefits, both in the ability to achieve net-zero goals and in advancing other goals. The results present the case for comprehensive advanced planning at the portfolio level to prioritize investments that will balance the minimization of emissions and life cycle cost with the maximization of resilience and environmental justice benefits. The novel methods for evaluation and integration, valuation of benefits, and consideration at the portfolio scale allow organizations to select investments that simultaneously address multiple key priorities.

Net-Zero Emissions↗
DOE OSTI2021

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Use of a Multiple Hydride Donor To Achieve an n-Doped Polymer with High Solvent Resistance

We report the ability to insolubilize doped semiconducting polymer layers can help enable the fabrication of efficient multilayer solution-processed electronic and optoelectronic devices. Here, we present a promising approach to simultaneously n-dope and largely insolubilize conjugated polymer films using tetrakis[{4-(1,3-dimethyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)phenoxy}methyl]methane (tetrakis-O-DMBI-H), which consists of four 2,3-dihydro-1H-benzoimidazole (DMBI-H) n-dopant moieties covalently linked to one another. Doping a thiophene-fused benzodifurandione-based oligo(p-phenylenevinylene)-co-thiophene polymer (TBDOPV-T) with tetrakis-O-DMBI-H results in a highly n-doped film with bulk conductivity of 15 S cm –1 . Optical absorption spectra provide evidence for film retention of ~93% after immersion in o-dichlorobenzene for 5 min. The optical absorption signature of the charge carriers in the n-doped polymer decreases only slightly more than that of the neutral polymer under these conditions, indicating that the exposure to solvent also results in negligible dedoping of the film. Moreover, thermal treatment studies on a tetrakis-O-DMBI-H-doped TBDOPV-T film in contact with another undoped polymer film indicate immobilization of the molecular dopant in TBDOPV-T. This is attributed to the multiple electrostatic interactions between each dopant tetracation and up to four nearby anionic doped polymer segments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Visualizing the Internal Nanocrystallinity of Calcite Due to Nonclassical Crystallization by 3D Coherent X–Ray Diffraction Imaging

The internal crystallinity of calcite is investigated for samples synthesized using two approaches: precipitation from solution and the ammonium carbonate diffusion method. Scanning electron microscopy (SEM) analyses reveal that the calcite products precipitated using both approaches have a well–defined rhombohedron shape, consistent with the euhedral crystal habit of the mineral. The internal structure of these calcite crystals is characterized using Bragg coherent diffraction imaging (BCDI) to determine the 3D electron density and the atomic displacement field. BCDI reconstructions for crystals synthesized using the ammonium carbonate diffusion approach have the expected euhedral shape, with internal strain fields and few internal defects. In contrast, the crystals synthesized by precipitation from solution have very complex external shapes and defective internal structures, presenting null electron density regions and pronounced displacement field distributions. These heterogeneities are interpreted as multiple crystalline domains, created by a nonclassical crystallization mechanism, where smaller nanoparticles coalescence into the final euhedral particles. The combined use of SEM, X–ray diffraction (XRD), and BCDI allows for structurally differentiating calcite crystals grown with different approaches, opening new opportunities to understand how grain boundaries and internal defects alter calcite reactivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CRCNS22 Learning Rules in the Hippocampus and their Mapping to Neuromorphic Systems (Final Technical Report)

Large scale biologically-realistic computational models are key to investigating the interplay between structure and function in nervous systems, thus paving the way to new clinical methods and neuro-inspired computing solutions. This project focuses on the hippocampus, in particular the CA3-CA1 regions, due to their role in associative learning and memory, pattern separation and completion, and spatial navigation. Investigations into the neuronal organization and learning rule(s) of this circuit can shed light into how declarative memories are formed, stored, recalled and forgotten and inform computational, experimental and clinical neuroscience work. Our project aims at developing a novel data-driven methodology supported by a broad heterogeneous base of neuroscience experimental knowledge and inspired from advances in computer science and engineering. Specifically, this work will benchmark existing and new learning rules within a full-scale spiking neural network simulation of the CA3-CA1 region. The model will be based on an open-source repository, called the Hippocampome, which contains neuronal morphologies, firing patterns, synapse probabilities, and most other required parameters for all known neuron types in the rodent hippocampal formation. The model will be first trained in a supervised fashion for associative memory tasks using backpropagation through time traditionally used in computer science, enhanced with a new technique called the surrogate gradient method. This optimization method will be used to obtain a global loss minimization, but it is not biologically inspired as it assumes the use of data not locally available to the synapses. However, we propose its use as a benchmarking tool, to compare the training performance of local biologically plausible and hardware-mappable learning rules at scale. New rules or combinations will be proposed and tested as needed, based on the obtained results. Progress in this area will also drive the development of novel hardware-mappable algorithms for continual lifelong learning and categorization of new events from few presented examples. This project goes beyond the existing state-of-the-art by looking at large scale realistic neuronal circuits as networks trainable via global optimization methods such as surrogate gradient descent. The objective function of the brain that supports learning is largely unknown, but it is likely that it operates through local learning rules. Studying network trajectories around local minima as proposed in this work represents a useful strategy for understanding whether a network is training by using a specific (set of) learning rule(s). Starting from a completely untrained network is a challenging test since it is difficult to determine how the learning rule affects the trajectory of the network. This interdisciplinary project will help understand what rule governs learning in these regions or if multiple learning rules are involved. The work will develop a robust methodology to measure if the network is converging to the target solution, oscillating around it, or diverging away.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Realizing Complete Solid-Solution Reaction in High Sodium Content P2-Type Cathode for High-Performance Sodium-Ion Batteries

P2-type layered oxides suffer from an ordered Na + /vacancy arrangement and P2→O2/OP4 phase transitions, leading them to exhibit multiple voltage plateaus upon Na + extraction/insertion. The deficient sodium in the P2-type cathode easily induces the bad structural stability at deep desodiation states and limited reversible capacity during Na + de/insertion. These drawbacks cause poor rate capability and fast capacity decay in most P2-type layered oxides. In this work, to address these challenges, a novel high sodium content (0.85) and plateau-free P2-type cathode-Na 0.85 Li 0.12 Ni 0.22 Mn 0.66 O 2 (P2-NLNMO) was developed. The complete solid-solution reaction over a wide voltage range ensures both fast Na + mobility (10 -11 to 10 -10 cm 2 s -1 ) and small volume variation (1.7%). Furthermore, the high sodium content P2-NLNMO exhibits a higher reversible capacity of 123.4 mAhg -1 , superior rate capability of 79.3 mAhg -1 at 20 C, and 85.4% capacity retention after 500 cycles at 5 C. The sufficient Na and complete solid-solution reaction are critical to realizing high-performance P2-type cathodes for sodium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Scalable and accurate multi-GPU-based image reconstruction of large-scale ptychography data

Abstract While the advances in synchrotron light sources, together with the development of focusing optics and detectors, allow nanoscale ptychographic imaging of materials and biological specimens, the corresponding experiments can yield terabyte-scale volumes of data that can impose a heavy burden on the computing platform. Although graphics processing units (GPUs) provide high performance for such large-scale ptychography datasets, a single GPU is typically insufficient for analysis and reconstruction. Several works have considered leveraging multiple GPUs to accelerate the ptychographic reconstruction. However, most of these works utilize only the Message Passing Interface to handle the communications between GPUs. This approach poses inefficiency for a hardware configuration that has multiple GPUs in a single node, especially while reconstructing a single large projection, since it provides no optimizations to handle the heterogeneous GPU interconnections containing both low-speed (e.g., PCIe) and high-speed links (e.g., NVLink). In this paper, we provide an optimized intranode multi-GPU implementation that can efficiently solve large-scale ptychographic reconstruction problems. We focus on the maximum likelihood reconstruction problem using a conjugate gradient (CG) method for the solution and propose a novel hybrid parallelization model to address the performance bottlenecks in the CG solver. Accordingly, we have developed a tool, called PtyGer ( Pty chographic G PU(multipl e )-based r econstruction), implementing our hybrid parallelization model design. A comprehensive evaluation verifies that PtyGer can fully preserve the original algorithm’s accuracy while achieving outstanding intranode GPU scalability.

97 MATHEMATICS AND COMPUTING↗