Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 235 records · Page 13

Nucleation and growth of polar clusters with in-phase tilts into a long-range ferroelectric matrix in a sodium niobate based complex relaxor

In this study, we have investigated the temperature dependence of atomic ordering at multiple length scales in a lead-free sodium niobate-based relaxor, i.e., 0.75 NaNbO 3 -0.25 Ba 0.9⁢ Ca 0.1⁢ TiO 3 (NN-25BCT) via synchrotron x-ray diffraction, Raman spectroscopy, and pair distribution function analysis. High-resolution synchrotron x-ray powder diffraction (SXRD) measurements reveal a ferroelectric phase transition in the relaxor ferroelectric NN-25BCT below the Vogel-Fulcher freezing temperature (𝑇 VF ≈ 270 K). In addition, SXRD analysis demonstrates the competition between in-phase octahedral tilting and ferroelectric order at the long-range scale using mode crystallography. On the other hand, Raman spectroscopic analysis provides evidence of polar ordering for 𝑇 > 𝑇 VF (with tetragonal symmetry) persisting up to the Burns temperature (𝑇 B ). Furthermore, pair distribution function (PDF) analysis reveals the presence of a polar antiferrodistortive tetragonal phase with 𝑃⁢4⁢𝑏𝑚 space group at short ranges throughout the studied temperatures (i.e., 110 K ≤ 𝑇 ≤500 K), irrespective of nonpolar long-range ordering above 𝑇 VF . Therefore, our measurements provide direct evidence for the presence of polar ordering at short ranges and their gradual transformation into long-range polar ordering using an integrated multiscale structural analysis. In conclusion, as a result of a transition from relaxor to a ferroelectric phase in the vicinity of room temperature, NN-25BCT can be exploited for applications in pyroelectric detectors, electrocaloric devices, and multilayered ceramic capacitors.

36 MATERIALS SCIENCE↗

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗

Improved Three-Dimensional Reconstructions in Electron Ptychography through Defocus Series Measurements

A detailed analysis of ptychography for three-dimensional (3D) phase reconstructions of thick specimens is performed. Here, we introduce multi-focus ptychography, which incorporates a 4D-STEM defocus series to enhance the quality of 3D reconstructions along the beam direction through a higher overdetermination ratio. This method is compared with established multi-slice ptychography techniques, such as conventional ptychography, regularized ptychography, and multi-mode ptychography. Additionally, we contrast multi-focus ptychography with an alternative method that uses virtual optical sectioning through a reconstructed scattering matrix (S-matrix), which offers more precise 3D structure information compared to conventional ptychography. Our findings from multiple 3D reconstructions based on simulated and experimental data demonstrate that multi-focus ptychography surpasses other techniques, particularly in accurately reconstructing the surfaces and interface regions of thick specimens.

3D reconstruction↗

Elastic proton scattering off nonzero spin nuclei

In recent years, we constructed a microscopic optical potential (OP) for elastic nucleon-nucleus ( NA ) scattering using modern approaches based on chiral theories for the nucleon-nucleon ( NN ) interaction. The OP was derived at first order of the spectator expansion in Watson multiple scattering theory and its final expression was a folding integral between the NN t matrix and the nuclear density of the target. Two- and three-body forces are consistently included both in the target and in the projectile description. The purpose of this work is to apply our microscopic OP to nuclei characterized by a ground state of spin-parity quantum numbers J π ≠ 0 + . We extended our formalism to include the spin of the target nucleus. The full amplitudes of the NN reaction matrix are retained in the calculations starting from two- and three-body chiral forces. The microscopic OP can be applied in the energy range 100 ≤ E ≤ 350 MeV. Finally, we show a remarkable agreement with experimental data for the available observables and, simultaneously, provide reliable estimates for the theoretical uncertainties. This work paves the way toward a full microscopic approach to inelastic NA scattering, showing that the derivation of optical potentials between states with J π ≠ 0 + is completely under control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X-ray scattering based scanning tomography for imaging and structural characterization of cellulose in plants

X-ray and neutron scattering have long been used for structural characterization of cellulose in plants. Due to averaging over the illuminated sample volume, these measurements traditionally overlooked the compositional and morphological heterogeneity within the sample. Here, a scanning tomographic imaging method is described, using contrast derived from the X-ray scattering intensity, for virtually sectioning the sample to reveal its internal structure at a resolution of a few micrometres. This method provides a means for retrieving the local scattering signal that corresponds to any voxel within the virtual section, enabling characterization of the local structure using traditional data-analysis methods. This is accomplished through tomographic reconstruction of the spatial distribution of a handful of mathematical components identified by non-negative matrix factorization from the large dataset of X-ray scattering intensity. Joint analysis of multiple datasets, to find similarity between voxels by clustering of the decomposed data, could help elucidate systematic differences between samples, such as those expected from genetic modifications, chemical treatments or fungal decay. The spatial distribution of the microfibril angle can also be analyzed, based on the tomographically reconstructed scattering intensity as a function of the azimuthal angle.

36 MATERIALS SCIENCE↗

Analyzing the Effects of Cyberattacks on Distribution System State Estimation

Key components of power systems—such as energy management systems, automatic generation control, and state estimation—are under serious vulnerability from cyberattacks. Cyber threats in electric grids have increased significantly because of the increased interconnectivity of supervisory control and data acquisition systems and public network infrastructure. As the penetration level of distributed energy resources increases, it is imperative to employ system-monitoring techniques such as state estimation for the reliable operation of distribution systems. Recently, multiple methods have been developed that exploit the low rank property of distribution system state matrix and are robust to bad data, such as matrix completion. This paper analyzes the impact of various realistic cyberattack scenarios on matrix completion. Realistic cyberattack scenarios are converted into data corruption models that are used in an extensive simulation of a custom IEEE 123-bus system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

End-to-End Pipeline for Trigger Detection on Hit and Track Graphs

There has been a surge of interest in applying deep learning in particle and nuclear physics to replace labor-intensive offline data analysis with automated online machine learning tasks. This paper details a novel AI-enabled triggering solution for physics experiments in Relativistic Heavy Ion Collider and future Electron-Ion Collider. The triggering system consists of a comprehensive end-to-end pipeline based on Graph Neural Networks that classifies trigger events versus background events, makes online decisions to retain signal data, and enables efficient data acquisition. Here, the triggering system first starts with the coordinates of pixel hits lit up by passing particles in the detector, applies three stages of event processing (hits clustering, track reconstruction, and trigger detection), and labels all processed events with the binary tag of trigger versus background events. By switching among different objective functions, we train the Graph Neural Networks in the pipeline to solve multiple tasks: the edge-level track reconstruction problem, the edge-level track adjacency matrix prediction, and the graph-level trigger detection problem. We propose a novel method to treat the events as track-graphs instead of hit-graphs. This method focuses on intertrack relations and is driven by underlying physics processing. As a result, it attains a solid performance (around 72% accuracy) for trigger detection and outperforms the baseline method using hit-graphs by 2% higher accuracy.

97 MATHEMATICS AND COMPUTING↗

Nanopore Activity Assays for Detection of Biomarker Protease Activity: Design and Testing of Substrates for Both Nanopore Sequencing and PCR-Based Detection Methods

The work performed in this project has demonstrated the ability to construct proteolytic enzyme substrates that are PCR and sequencing-readable reporter molecules. Specifically, the goal was to detect those reporter molecules via PCR and Oxford Nanopore Technologies MinION sequencing methods following exposure to the biomarker protease thrombin. The assay development focused on binding the constructed peptide-oligonucleotide chimera to immobilized streptavidin. The action of thrombin on the peptide portion of the molecule released the oligonucleotide for detection. Detection of protease activity was demonstrated in a concentration-dependent manner using MALDI-MS, RT-PCR and DNA sequencing. Additional steps to remove background release of reporter molecules during the assay was used to improve the difference in detected oligonucleotide reporter following protease activity. Additional steps in assay development will be to (1) test the assay in an appropriate matrix, (2) investigate detection using additional DNA sequencing platforms and (3) demonstrate multiplexed detection of multiple protease markers in a single reaction.

59 BASIC BIOLOGICAL SCIENCES↗

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗

Analyzing the Effects of Cyberattacks on Distribution System State Estimation: Preprint

Key components of power systems—such as energy management systems, automatic generation control, and state estimation—are under serious vulnerability from cyber attacks. Cyber threats in electric grids have increased significantly because of the increased interconnectivity of supervisory control and data acquisition systems and public network infrastructure. As the penetration level of distributed energy resources increases, it is imperative to employ system-monitoring techniques such as state estimation for the reliable operation of distribution systems. Recently, multiple methods have been developed that exploit the low rank property of distribution system state matrix and are robust to bad data, such as matrix completion. This paper analyzes the impact of various realistic cyber attack scenarios on matrix completion. Realistic cyber attack scenarios are converted into data corruption models that are used in an extensive simulation of a custom IEEE 123-bus system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Response Surface Methodology As a New Approach for Finding Optimal MALDI Matrix Spraying Parameters for Mass Spectrometry Imaging

Automated spraying devices have become ubiquitous in laboratories employing matrixassisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in part because they permit control of a number of matrix application parameters that can easily be reproduced for intra- and interlaboratory studies. Determining the optimal parameters for MALDI matrix application, such as temperature, flow rate, spraying velocity, number of spraying cycles, and solvent composition for matrix application, is critical for obtaining high-quality MALDI-MSI data. However, there are no established approaches for optimizing these multiple parameters simultaneously. Instead optimization is performed iteratively (i.e., one parameter at a time), which is time consuming and can lead to overall non-optimal settings. In this report, we demonstrate the use a novel experimental design and the response surface methodology to optimize five parameters of MALDI matrix application using a robotic sprayer. Thirtytwo combinations of MALDI matrix spraying conditions were tested, which allowed us to elucidate relationships between each of the application parameters as determined by MALDI-MS (specifically, using a 15 Tesla Fourier transform ion cyclotron resonance mass spectrometer). As such, we were able to determine the optimal automated spraying parameters that minimized signal de-localization and enabled high MALDI sensitivity. We envision this optimization strategy can be utilized for other MALDI-MSI applications, tissue types, and matrix application approaches.

60 APPLIED LIFE SCIENCES↗

Rapid AI-based Dissection of Ashes using Raman and XRF Spectroscopy (RADAR-X)

Waste-to-energy (WTE) facilities incinerate ~35 million tons of municipal solid waste annually in the United States. The incineration process reduces the mass and the volume of the waste fraction by over 75 and 95%, respectively. The fraction left after incineration remains as ash residues and is referred to as WTE ash, compromising of bottom and fly ash. In the United States, ~10 million tons of WTE ashes are generated annually and predominantly landfilled because of the lack of secondary end-use pathways. This incurs a significant financial burden (landfilling costs) on U.S. WTE facilities and also results in the loss of materials to the landfill. The primary objective of this research is to understand better the elemental and mineralogical composition of WTE ashes from diverse origins and find composition dependent upcycling pathways for diverting these ashes from landfills. This primary objective was addressed through three research tasks –(Task I) An AI-led Multi-Modal Approach for Compositional Analysis, (Task II) Developing a Dissolution-Based Test for Real Time Analysis, and (Task III) Establishing composition-dependent end uses. The chemical composition of WTE ash is dependent on two factors, i.e., the input waste composition and the operational parameters of a WTE facility (combustion conditions). Amongst these two factors, the input waste composition will likely show spatial and temporal variation. As a result, the chemical composition of WTE ash will also fluctuate. To understand the spatial and temporal variation in WTE ash composition, in Task I, we collected 128 ash samples (62 bottom ash and 66 fly ash samples) from 11 WTE facilities located in 11 U.S. states and characterized them via X-ray Fluorescence, powder X-ray Diffraction, and Raman Spectroscopy. The findings from this extensive characterization work indicated that the key elements in WTE fly ashes are Ca, Cl (greater than 10 wt. %), Si, S, K, Zn ( between 1 and 10 wt. %), Mg, Al, P, Ti, Fe, Cu, Br, and Pb (between 0.1 and 1 wt. %). Similarly, the key elements in WTE bottom ash fraction finer than 45μm are Ca (greater than 10 wt. %), Mg, Al, Si, S, Cl, K, Ti, Fe, Zn (between 1 and 10 wt. %), P, V, Cr, Mn, Cu, Br, and Pb (between 0.1 and 1 wt. %). Here, we note that the dominant fraction of WTE bottom ash is the coarse fraction. The coarse WTE bottom ash fraction (rich in silicon) was not characterized in this study because of excessive grinding requirements and their unsuitability as a supplementary cementitious material due to their coarse nature. The elements in WTE bottom ashes are present as calcite, anhydrite, vaterite, hydroxyapatite, quartz, bassanite, gehlenite, akermanite, hydrocalumite, and portlandite. Similarly, the mineralogical species present in WTE fly ashes are calcium chloride hydroxide, halite, calcite, anhydrite, sylvite, hydrocalumite, vaterite, hannebachite, quartz, and bassanite. Temporal variation in ash composition may also result in significant fluctuations in chemical compositions. Therefore, a WTE facility may need to monitor the ash composition (elemental and mineralogical composition) in real time. In Task I, we evaluated the possibility of using a portable X-ray fluorescence (XRF) spectrometer to monitor the elemental composition in real-time. Specifically, we collected XRF data on identical specimens via a portable XRF spectrometer (low-end) and a lab-based benchtop XRF spectrometer (high-end). The collected data was used to train an A.I. algorithm (portable XRF data as an input and benchtop XRF data as an output) to predict accurate elemental composition using portable XRF data. Finally, we developed a 2-minute photobleaching protocol to monitor the mineralogical characteristics of WTE ashes via Raman spectroscopy. Overall, the activities in Task I improved our understanding of ash composition and developed techniques to monitor elemental and mineralogical composition in near real-time. Based on the findings of Task I, we find that WTE ashes exhibit wide variability in mineralogy. For ICP-based elemental analysis, all the mineralogical species in WTE ashes must be brought into solution. This is traditionally accomplished with acid digestion using a combination of multiple acids. However, acid digestion with multiple acids is time-consuming and often fails to ensure complete digestion of the ash matrix. To address this limitation, in Task II, we developed an alkali-fusion-based digestion protocol using a combination of lithium tetraborate, lithium metaborate, and their combinations as possible alkali fluxes for digesting WTE ashes entirely and rapidly. The validity of the developed method was evaluated on two standard ash specimens, i.e., SRM 1633c coal fly ash and BCR-176R incineration fly ash specimen. The findings suggest that the developed protocol can ensure complete digestion of elements such as Al, Ba, Ca, Cr, Cu, Mg, Mn, P, Sr, V, Zn, Be, K, and rare earth elements. The recent changes in the energy market towards renewables and increased metal recycling have resulted in reduced supplies of supplementary cementitious materials (coal fly ash and slag). Therefore, in Task III, we evaluated the possibility of employing WTE ashes as SCMs. As the chemical composition of WTE ashes varies temporally (on an hourly basis), there was also a need to develop tests that can evaluate the suitability of material to act as supplementary cementitious material rapidly, i.e., in a few minutes. Therefore, in Task II, we also developed a rapid test to assess the suitability of a material to act as an SCM in 5 minutes. This represents a significant advance over the state-of-the-art R 3 test, which takes ~144 hours. This test was initially validated on amorphous aluminosilicates, such as calcined clays, and could be extended to evaluate other industrial by-products, such as WTE ashes. In Task III, we evaluated the possibility of employing WTE ashes for two applications, i.e., as an SCM and a lime substitute for clay stabilization. The findings from Task I indicated that WTE ashes are enriched in chlorine and, therefore, cannot be used directly as an SCM due to corrosion-related risks and altered hydration kinetics. Accordingly, we developed an ash treatment protocol to reduce the solubility of chlorine-containing species in WTE ashes. The developed treatment protocol also immobilized lead in certain mineral forms. As a result of the treatment, WTE ashes can be used as SCMs without any corrosion or heavy metal leaching concerns. The second application examined in this study was clay stabilization. WTE ashes are calcium-rich and can be an adequate lime replacement for clay stabilization. Our findings reveal that the sum of the concentrations of Ca(OH) 2 and CaClOH controls the clay stabilization capability of WTE ashes. In summary, in this work, we evaluated the elemental and mineralogical characteristics of U.S. WTE ashes from diverse origins and developed tests to evaluate the chemical characteristics of these ashes in real time through a portable XRF and a benchtop Raman spectrometer. Based on the chemical characteristics of these ashes, we developed an ash treatment process to enable the use of WTE ashes as an SCM and also evaluated the possibility of employing these ashes for clay stabilization. Overall, the findings from this work enables the diversion of WTE ashes from landfills for multiple end-uses, i.e., as an SCM or a lime substitute for clay stabilization.

36 MATERIALS SCIENCE↗

A graph signal processing‐based multiple model Kalman filter ( GSP‐MMKF ) tool for predictive analytics: An air separation unit process application

Abstract The industrial Air Separations Unit (ASU) is a complicated and tightly operated process. The use of dynamic process analytics is also a key element of safe and economic operation of these processes, with increasing focus on predictive analytics to take preemptive actions. With the availability of real‐time data from hundreds of sensors, the data analysis process should also consider the topology of the data, as seen in sensor networks. In this paper, a novel tool is presented that considers the complex connectivity patterns in the sensor network and uses local adaptive disturbance estimations to predict global network‐scale trends. The paper introduces the emerging field of Graph Signal Processing (GSP) and presents a rigorous derivation of the tool starting from the extraction of the sensor‐network (in a graph theoretical sense) from the data. This network, which is in the form of a matrix, is then used to derive a Kalman‐filter type of state‐space model driven by input disturbances. Multiple disturbance models (e.g., step, ramp, periodic) are included to allow the model to have different kinds of disturbance propagation. Each graph node (representing the sensors used) dynamically adapts to the most recent detected disturbance individually. These estimated disturbances are propagated to the global network using the graph. Modifications to ensure stability are also discussed. The fidelity of the tool is tested on certain downtime events and the paper concludes by discussing the advantages of the method and planned future improvements.

Ghosh, Sambit↗

Single Quadrupole Multiple Fragment Ion Monitoring Quantitative Mass Spectrometry

Single quadrupole mass spectrometry (MS) with enhanced in-source multiple fragment ion monitoring was designed to perform high sensitivity quantitative mass analyses. Enhanced in-source fragmentation amplifies fragmentation from traditional soft electrospray ionization producing fragment ions that have been found to be identical to those generated in tandem MS. We have combined enhanced in-source fragmentation data with criteria established by the European Union Commission Directive 2002/657/EC for electron ionization single quadrupole quantitative analysis to perform quantitative analyses. These experiments were performed on multiple types of complex samples that included a mixture of 50 standards, as well as cell and plasma extracts. The dynamic range for these quantitative analyses was comparable to triple quadrupole multiple reaction monitoring (MRM) analyses at up to 5 orders of magnitude with the cell and plasma extracts showing similar matrix effects across both platforms. Amino acid and fatty acid measurements performed from certified NIST 1950 plasma with isotopically labeled standards demonstrated accuracy in the range of 91–110% for the amino acids, 76–129% for the fatty acids, and good precision (coefficient of variation <10%). To enhance specificity, a newly developed correlated ion monitoring algorithm was designed to facilitate these analyses. This algorithm autonomously processes, aligns, filters, and compiles multiple ions within one chromatogram enabling both precursor and in-source fragment ions to be correlated within a single chromatogram, also enabling the detection of coeluting species based on precursor and fragment ion ratios. Single quadrupole instrumentation can provide MRM level quantitative performance by monitoring/correlating precursor and fragment ions facilitating high sensitivity analysis on existing single quadrupole instrumentation that are generally inexpensive, easy to operate, and technically less complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-matrix correlators and localization

Abstract We study generating functions of$$ \frac{1}{4} $$ 1 4 -BPS states in$$ \mathcal{N} $$ N = 4 super Yang-Mills at finiteNby attempting to generalize the Harish-Chandra-Itzykson-Zuber integral to multiple commuting matrices. This allows us to compute the overlaps of two or more generating functions; such calculations arise in the computation of two-point correlators in the free-field limit. We discuss the four-matrix HCIZ integral in the U(2) context and lay out a prescription for finding a more general formula forN> 2. We then discuss its connections with the restricted Schur polynomial operator basis. Our results generalize readily to arbitrary numbers of matrices, opening up the opportunity to study more generic BPS operators.

Physics↗

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗

Extended cycle life implications of fast charging for lithium-ion battery cathode

Enabling extreme fast charging (XFC, =10–15 min charging) requires a comprehensive understanding of its implications. While lithium plating is a key bottleneck for the anode, the full extent of limitations for the cathode are not well-understood, particularly in extended-cycle settings with well-defined battery designs and conditions. This article presents cycle-life implications of XFC on cathodes at multiple length scales, combining electrochemical analyses, degradation modeling, and post-test characterizations. The comprehensive test matrix includes 41 well-defined gr/NMC pouch cells under varied fast-charge rates (1–9C) and state-of-charges cycled up to 1000 times. Cathode issues remain minimal in early cycling, but begin to accelerate in later life, when distinct cracking is found and identified as a fatigue mechanism. The bulk structure of cathodes remains intact, but distinct particle surface reconstruction is observed; however, this shows less pronounced effect on cathode aging than does cracking.

25 ENERGY STORAGE↗

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗