Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear projection”

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 91 records · Page 5

Lithium-Ion Battery Diagnostics Using Electrochemical Impedance via Machine-Learning

Diagnosing battery states such as health, state-of-charge, or temperature is crucial for ensuring the safety and reliability of electrochemical energy storage systems. While some states, such as temperature, may be measured using cheap sensors, accurate diagnosis of battery health metrics usually requires time-consuming performance measurements, making them infeasible for use in real-world operation. These health metrics can be measured during lab-testing and then estimated on-line using predictive life models or via state observer algorithms such as Kalman filters, but these predictive methods should be supplemented by actual measurement of battery health whenever possible to ensure reliability. Rapid measurement of battery health may be done by various types of fast diagnostic techniques such as electrochemical impedance spectroscopy (EIS), which can be performed in only a few minutes and require only a fraction of the energy and power needed for a full charge and discharge measurement. But there is a substantial challenge for estimating battery health using EIS data, as EIS is sensitive to cell temperature, state-of-charge, current, and resting time in addition to health. Thus, utilizing EIS data to predict battery capacity requires correcting for all these additional variables, a task that is extremely difficult to handle analytically. This talk utilizes machine-learning methods to estimate the effectiveness of battery capacity prediction from EIS data, leveraging a data set of hundreds of EIS measurements recorded at varying temperature and state-of-charge throughout a 500-day aging study of 32 commercial, large-format NMC-Graphite lithium-ion batteries. Using EIS as input to machine-learning models is complicated by the nonlinear response of impedance to battery health, temperature, and state-of-charge, as well as the collinearity between the impedance response at neighboring frequencies, which can easily lead to overfit models. To train robust models, features from EIS data need to be extracted from the data or some subset of critical frequencies selected. Many approaches for extracting and selecting features from EIS data from electrochemical analysis and machine-learning fields were identified for analysis: using the entire raw spectra; selection of one, two, or many frequencies from the entire spectra; selecting interesting points from the EIS measurement using domain knowledge; fitting EIS with an equivalent-circuit model; calculating statistics on the raw impedance values; and reducing the dimensionality of the data using unsupervised linear (principal component analysis) and non-linear (uniform manifold approximation and projection) methods. These approaches were rigorously compared using a machine-learning pipeline approach, training linear, Gaussian process, and random forest regression models and quantifying performance using cross-validation as well as a held-out test set. An artificial neural network model trained on the raw spectra was also tested. Promising pipelines were fine-tuned via Bayesian hyperparameter optimization using cross-validation loss and training with class-specific weights to counter data set imbalance. The most reliable method for utilizing impedance in this work was the selection of two optimal frequencies through an exhaustive search, resulting in about 2% mean absolute error on test data for both Gaussian process and random forest model architectures. Interrogation of a variety of models reveals critical frequencies of 100 Hz and 103 Hz for this data set, though the optimal set of frequencies is not necessarily intuitive, i.e., the best performing models are not simply those that use impedance at frequencies that have the highest correlation to the relative discharge capacity. The best performing model is an ensemble model, which is able to predict battery capacity with 1.9% mean absolute error for unseen cells using impedance recorded at a variety of temperatures and states-of-charge.

battery↗

Another Geometric Interpretation of Cramer’s Rule

We develop a geometric interpretation of Cramer’s rule as a generalization of projection onto orthogonal basis vectors using the rows of the adjugate. This interpretation makes connections between elementary linear algebra concepts like the solution to linear equations, inner products, and projections. Such connections are useful for introducing broader concepts related to Hilbert spaces and geometric algebras like Grassman algebra. Such connections were essential for the author’s mathematical education as an engineer.

linear algebra↗

Uncovering interpretable relationships in high-dimensional scientific data through function preserving projections

Abstract In many fields of science and engineering, we frequently encounter experiments or simulations datasets that describe the behavior of complex systems and uncovering human interpretable patterns between their inputs and outputs via exploratory data analysis is essential for building intuition and facilitating discovery. Often, we resort to 2D embeddings for examining these high-dimensional relationships (e.g. dimensionality reduction). However, most existing embedding methods treat the dimensions as coordinates for samples in a high-dimensional space, which fail to capture the potential functional relationships, and the few methods that do take function into consideration either only focus on linear patterns or produce non-linear embeddings that are hard to interpret. To address these challenges, we proposed function preserving projections (FPP), which construct 2D linear embeddings optimized to reveal interpretable yet potentially non-linear patterns between the domain and the range of a high-dimensional function. The intuition here is that humans are good at understanding potentially non-linear patterns in 2D but unable to interpret non-linear mapping from high-dimensional space to 2D. Therefore, we should restrict the projection to linear but not the pattern we are seeking. Using FPP on real-world datasets, one can obtain fundamentally new insights about high-dimensional relationships in extremely large datasets that could not be processed with existing dimension reduction methods.

97 MATHEMATICS AND COMPUTING↗

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ARPIST: Provably accurate and stable numerical integration over spherical triangles

Numerical integration on spherical triangles, including the computation of their areas, is a core computation in geomathematics. The commonly used techniques sometimes suffer from instabilities and significant loss of accuracy. We describe a new algorithm, called ARPIST, for accurate and stable integration of functions on spherical triangles. ARPIST is based on an easy-to-implement transformation to the spherical triangle from its corresponding linear triangle via radial projection to achieve high accuracy and efficiency. More importantly, ARPIST overcomes potential instabilities in computing the Jacobian of the transformation, even for poorly shaped triangles that may occur at poles in regular longitude-latitude meshes, by avoiding potential catastrophic rounding errors. We compare our proposed technique with L’Huilier’s Theorem for computing the area of spherical triangles, and also compare it with the recently developed LSQST method (Beckmann et al., 2014) and a radial-basis-function-based technique (Reeger and Fornberg, 2016) for integration of smooth functions on spherical triangulations. In conclusion, our results show that ARPIST enables better or comparable accuracy over previous methods while being easier to implement, significantly faster, and more tolerant of poor element shapes.

97 MATHEMATICS AND COMPUTING↗

Projecting Stratocumulus Transitions on the Albedo—Cloud Fraction Relationship Reveals Linearity of Albedo to Droplet Concentrations

Satellite images show solid marine stratocumulus cloud decks (Sc) that break up over the remote oceans. The Sc breakup is initiated by precipitation and is accompanied by a strong reduction in the cloud radiative effect. Aerosol has been shown to delay the Sc breakup by postponing the onset of precipitation, however its climatic effect is uncertain. Here we introduce a new approach that allows us to re-cast currently observed cloud cover and albedo to their counterfactual cleaner world, enabling the first estimate of the radiative effect due to delayed cloud breakup. Using simple radiative approximation, the radiative forcing with respect to pre-industrial times due to delayed Sc breakup is –0.39 W m –2 . The radiative effect changes nearly linearly with aerosol due to the droplet concentration control on the cloud cover, suggesting a potentially accelerated warming if the current trend of reduction in aerosol emissions continues.

54 ENVIRONMENTAL SCIENCES↗

An interactive machine learning platform for analyzing multi-particle coincidence data from cold target recoil ion momentum spectroscopy

We present SCULPT (Supervised Clustering and Uncovering Latent Patterns with Training), a comprehensive software platform for analyzing tabulated high-dimensional multi-particle coincidence data from Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS) experiments. The software addresses critical challenges in modern momentum spectroscopy by integrating advanced machine learning techniques with physics-informed analysis in an interactive web-based environment. SCULPT implements uniform manifold approximation and projection for non-linear dimensionality reduction to reveal correlations in high-dimensional data. We also discuss potential extensions to deep autoencoders for feature learning and genetic programming for automated discovery of physically meaningful observables. A novel adaptive confidence scoring system provides quantitative reliability assessments by evaluating user-selected clustering quality metrics with predefined weights that reflect each metric’s robustness. The platform features configurable molecular profiles for different experimental systems, interactive visualization with selection tools, and comprehensive data filtering capabilities. Utilizing a subset of SCULPT’s capabilities, we analyze photo-double-ionization data measured using the COLTRIMS method for three-body dissociation of the D 2 O molecule, revealing distinct fragmentation channels and their correlations with physics parameters. The software’s modular architecture and web-based implementation make it accessible to the broader atomic and molecular physics community, significantly reducing the time required for complex multi-dimensional analyses. This opens the door to finding and isolating rare events exhibiting non-linear correlations on the fly during experimental measurements, which can help steer exploration and improve the efficiency of experiments.

Artificial neural networks↗

An Upgrade Path for the Fermilab Accelerator Complex

The completion of the PIP-II project and its superconducting linear accelerator will provide up to 1.2 MW of beam power to the LBNF/DUNE facility for neutrino physics. It will also be able to produce high-power beams directly from the linac that can be used for lower-energy particle physics experiments as well, such as directing beam toward the Muon Campus at Fermilab for example. Any further significant upgrade of the beam power to DUNE, however, will be impeded by the limitations of the present Booster synchrotron at the facility. To increase the power to DUNE by a factor of two would require a new accelerator arrangement to feed the Main Injector that does not include the Booster. In what follows, a path toward upgrading the Fermilab accelerator complex to bring the beam power for DUNE to 2.4 MW is presented, using a new rapid-cycling synchrotron plus an energy upgrade to the PIP-II linac. The path includes the ability to instigate a new lower-energy, very high-power beam delivery system for experiments that can address much of the science program presented by the Booster Replacement Science Working Group. It also allows for the future possibility to go beyond 2.4 MW up to roughly 4 MW from the Main Injector.

43 PARTICLE ACCELERATORS↗

Scorpius Injector Solenoid Magnet Alignment and Characterization

The Scorpius project consists of a linear induction accelerator (LIA) that will be built in the coming few years at the Nevada test site. Large radius solenoid magnets are used to focus and transport the electron beam through the injector of Scorpius. Typically, solenoid magnets are wrapped directly around the beam pipe. They have an inner radius of about 8 cm and weight approximately 170 lbs. The magnets for the Scorpius injector are much larger, the largest of which extends to an inner radius of 45 cm with a weight of 709 lbs. In addition, the magnet measurement effort is typically time-consuming, taking several months to complete. With 5 different designs of 52 solenoid magnets to measure in under a year, this is a challenging undertaking that requires a systems solution. Due to the size, variability, and quantity of solenoid magnets, as well as the compressed schedule, the process of measuring the solenoid magnets before integrating them into their respective assemblies and shipping them to Nevada is a challenging problem. There is a need for the injector magnets to be characterized and aligned in preparation for the commissioning of Scorpius. This will ensure that the injector magnet installation meets requirements.

42 ENGINEERING↗

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field (Final Technical Report)

The project investigated 3-D linear instabilities of Harris current sheet under a realistic ion-to-electron mass ratio and a finite guide field. The excitation and 3-D properties of lower hybrid drift instability (LHDI), along with two new instabilities were demonstrated. The report presents the nonlinear properties of current sheet instabilities by employing Gyrokinetic electron and Fully kinetic ion (GeFi) particle simulation code. The GeFi simulation results are compared with fully kinetic PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ground-based microwave probing

Theoretical and experimental activities in remote sensing were conducted. The projects are: (1) the linear statistical inversion method, (2) vertical temperature profile determination from microwave emission measurements near 60 GHz, (3) radio path length correction determination from emission measurements at 20.6 GHz, and (4) liquid water content determination of thunderstorm cells by emission measurements at 10.7 GHz.

Westwater, E. R.↗

Radio structure of 3C147 determined by multi-element very long baseline interferometry

The radio structure of the quasar 3C147 has been determined from multi-baseline VLBI data at 609 MHz using both a conventional method and a technique which uses the 'closure' phase information to produce a good approximation to a synthesis map of the source. The structure is similar to the central part of M87, with a bright core and a linear 'jet' of a projected length of about 1.5 kpc which is concentrated in bright 'knots'.

Wilkinson, P. N.↗

Coded-aperture imaging in nuclear medicine

Coded-aperture imaging is a technique for imaging sources that emit high-energy radiation. This type of imaging involves shadow casting and not reflection or refraction. High-energy sources exist in x ray and gamma-ray astronomy, nuclear reactor fuel-rod imaging, and nuclear medicine. Of these three areas nuclear medicine is perhaps the most challenging because of the limited amount of radiation available and because a three-dimensional source distribution is to be determined. In nuclear medicine a radioactive pharmaceutical is administered to a patient. The pharmaceutical is designed to be taken up by a particular organ of interest, and its distribution provides clinical information about the function of the organ, or the presence of lesions within the organ. This distribution is determined from spatial measurements of the radiation emitted by the radiopharmaceutical. The principles of imaging radiopharmaceutical distributions with coded apertures are reviewed. Included is a discussion of linear shift-variant projection operators and the associated inverse problem. A system developed at the University of Arizona in Tucson consisting of small modular gamma-ray cameras fitted with coded apertures is described.

Smith, Warren E.↗

Project resource reallocation algorithm

A methodology for adjusting baseline cost estimates according to project schedule changes is described. An algorithm which performs a linear expansion or contraction of the baseline project resource distribution in proportion to the project schedule expansion or contraction is presented. Input to the algorithm consists of the deck of cards (PACE input data) prepared for the baseline project schedule as well as a specification of the nature of the baseline schedule change. Output of the algorithm is a new deck of cards with all work breakdown structure block and element of cost estimates redistributed for the new project schedule. This new deck can be processed through PACE to produce a detailed cost estimate for the new schedule.

Myers, J. E.↗

Spectroscopic quantification of projection effects in the SDSS redMaPPer galaxy cluster catalogue

ABSTRACT Projection effects, whereby galaxies along the line of sight to a galaxy cluster are mistakenly associated with the cluster halo, present a significant challenge for optical cluster cosmology. We use statistically representative spectral coverage of luminous galaxies to investigate how projection effects impact the low-redshift limit of the Sloan Digital Sky Survey (SDSS) redMaPPer galaxy cluster catalogue. Spectroscopic redshifts enable us to differentiate true cluster members from false positives and determine the fraction of candidate cluster members viewed in projection. Our main results can be summarized as follows: first, we show that a simple double-Gaussian model can be used to describe the distribution of line-of-sight velocities in the redMaPPer sample; secondly, the incidence of projection effects is substantial, accounting for ∼16 per cent of the weighted richness for the lowest richness objects; thirdly, projection effects are a strong function of richness, with the contribution in the highest richness bin being several times smaller than for low-richness objects; fourthly, our measurement has a similar amplitude to state-of-the-art models, but finds a steeper dependence of projection effects on richness than these models; and fifthly, the slope of the observed velocity dispersion–richness relation, corrected for projection effects, implies an approximately linear relationship between the true, three-dimensional halo mass and three-dimensional richness. Our results provide a robust, empirical description of the impact of projection effects on the SDSS redMaPPer cluster sample and exemplify the synergies between optical imaging and spectroscopic data for studies of galaxy cluster astrophysics and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Exponential integrators for non-linear diffusion

The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In this report, we present initial results of a proof of concept study that shows that exponential integrators can be an efficient alternative to traditional integration schemes.

42 ENGINEERING↗

Reduced basis approximations of parameterized dynamical partial differential equations via neural networks

Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. Here, we demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.

97 MATHEMATICS AND COMPUTING↗