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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Method for Fabricating Composite Structures Using Continuous Press Forming

A method for fabricating composite structures at a low-cost. moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a net-shape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electron-beam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates. resulting in lower cost and high structural performance.

Farley, Gary L.↗
NASA NTRS2000

Method for Fabricating Composite Structures Using Pultrusion Processing

A method for fabricating composite structures at a low-cost, moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a net-shape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electron-beam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates, resulting in lower cost and high structural performance.

Farley, Gary L.↗
NASA NTRS2000

Method for Fabricating Composite Structures Using Pultrusion Processing

A method for fabricating composite structures at a low-cost, moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a netshape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electronbeam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates, resulting in lower cost and high structural performance.

Farley, Gary L.↗
NASA NTRS2007

Fabrication of Polyimide-Matrix/Carbon and Boron-Fiber Tape

The term HYCARB denotes a hybrid composite of polyimide matrices reinforced with carbon and boron fibers. HYCARB and an improved process for fabricating dry HYCARB tapes have been invented in a continuing effort to develop lightweight, strong composite materials for aerospace vehicles. Like other composite tapes in this line of development, HYCARB tapes are intended to be used to build up laminated structures having possibly complex shapes by means of automated tow placement (ATP) - a process in which a computer-controlled multiaxis machine lays down prepreg tape or tows. The special significance of the present process for making dry HYCARB for ATP is that it contributes to the reduction of the overall cost of manufacturing boron-reinforced composite-material structures while making it possible to realize increased compression strengths. The present process for making HYCARB tapes incorporates a "wet to dry" process developed previously at Langley Research Center. In the "wet to dry" process, a flattened bundle of carbon fiber tows, pulled along a continuous production line between pairs of rollers, is impregnated with a solution of a poly(amide acid) in N-methyl-2-pyrrolidinone (NMP), then most of the NMP is removed by evaporation in hot air. In the present case, the polyamide acid is, more specifically, that of LaRC. IAX (or equivalent) thermoplastic polyimide, and the fibers are, more specifically, Manganite IM7 (or equivalent) polyacrylonitrile- based carbon filaments that have a diameter of 5.2 m and are supplied in 12,000-filament tows. The present process stands in contrast to a prior process in which HYCARB tape was made by pressing boron fibers into the face of a wet carbon-fiber/ poly(amide acid) prepreg tape . that is, a prepreg tape from which the NMP solvent had not been removed. In the present process, one or more layer(s) of side-by-side boron fibers are pressed between dry prepreg tapes that have been prepared by the aforementioned gwet to dry h process. The multilayer tape is then heated to imidize the matrix material and remove most of the remaining solvent, and is pressed to consolidate the multiple layers into a dense tape. For tests, specimens of HYCARB tapes and laminated composite panels made from HYCARB tape were prepared as follows: HYCARB tapes were fabricated as described above. Each panel was made by laying down ten layers of tape, containing, variously, one, two, or three boron-fiber plies and the remainder carbon- fiber-only plies (see figure). Each panel was made by laying down ten layers of tape. Each panel was then cured by heating to a temperature of 225 C for 15 minutes, then pressing at 200 psi (A1.4 MPa) while heating to 371 C, holding at 371 C for 1 hour, then continuing to hold pressure during cooling. Control specimens that were otherwise identical except that they did not contain boron fibers also were prepared. In room-temperature flexural tests, the HYCARB specimens performed comparably to the control specimens; in room-temperature, open-hole compression tests, the HYCARB specimens performed slightly better, by amounts that increased with boron content.

Belvin, Harry L.↗
DOE OSTI2020

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)↗
NASA NTRS2011

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗
DOE OSTI2020

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
DOE OSTI2021

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↗
NASA NTRS1972

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗
NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗
DOE OSTI2024

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↗
DOE OSTI2020

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↗
DOE OSTI2021

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↗
DOE OSTI2022

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↗
NASA NTRSDate not supplied

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗
NASA NTRSDate not supplied

Efficient Preconditioning of a High-Order Solver for Multiple Physics

This work addresses preconditioning approaches for an implicit high-order solver frame-work applied to multiple physics. The solver is based on a space-time spectral element method and matrix-free Newton-Krylov solver developed at NASA over the recent years. Within this context, most preconditioning methods are impractical, as the computational time and memory requirements scale poorly with increasing polynomial orders. To improve computational efficiency, we first describe a novel entity-based Block Jacobi preconditioner for the continuous-Galerkin solution of the linear-elasticity and linear-shell equations. Second, we introduce a multigrid algorithm to further reduce time-to-solution on stiff cases arising from continuous-and discontinuous-Galerkin discretizations. Results obtained on relevant single-physics reference solutions, demonstrate the feasibility of the methods, paving the way for high-order solutions of fully coupled multi-physics problems.

STMD↗