Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear filters”

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.

194 records · Page 11

A Multiwavelength Search for Intrinsic Linear Polarization in Wolf–Rayet Winds

Wolf–Rayet stars have strong, hot winds, with mass-loss rates at least a factor of 10 greater than their O-star progenitors, although their terminal wind speeds are similar. In this paper, we use the technique of multiband linear polarimetry to extract information on the global asymmetry of the wind in a sample of 47 bright Galactic WR stars. Our observations also include time-dependent observations of 17 stars in the sample. The path to our goal includes removing the dominating component of wavelength-dependent interstellar polarization (ISP), which normally follows the well-known Serkowski law. We include a wavelength-dependent ISP position angle parameter in our ISP law and find that 15 stars show significant results for this parameter. We detect a significant component of wavelength-independent polarization due to electron scattering in the wind for 10 cases, with most WR stars showing none at the ∼0.05% level precision of our data. The intrinsically polarized stars can be explained with binary interaction, large-scale wind structure, and clumping. We also found that 5 stars out of 19 observed with the Strömgren b filter (probing the complex λ4600–4700 emission-line region) have significant residuals from the ISP law and propose that this is due to wind asymmetries. We provide a useful catalog of ISP for 47 bright Galactic WR stars and upper limits on the possible level of intrinsic polarization.

79 ASTRONOMY AND ASTROPHYSICS↗

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗

Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment

We present a detection of 21 cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment. Radio observations acquired over 102 nights are used to construct maps that are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood ratio test, yields a detection significance of 7.1σ (LRG), 5.7σ (ELG), and 11.1σ (QSO). These are the first 21 cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (H I ), defined as ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}\equiv {10}^{3}\,{{\rm{\Omega }}}_{{\rm{H}}\,{\rm\small{I}}}\left({b}_{{\rm{H}}\,{\rm\small{I}}}+\langle \,f{\mu }^{2}\rangle \right)$, where Ω H I is the cosmic abundance of H I , b H I is the linear bias of H I , and $\langle$fμ 2 $\rangle$ = 0.552 encodes the effect of redshift-space distortions at linear order. We find ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.51}_{-0.97}^{+3.60}$ for LRGs (z = 0.84), ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={6.76}_{-3.79}^{+9.04}$ for ELGs (z = 0.96), and ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.68}_{-0.67}^{+1.10}$ for QSOs (z = 1.20), with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and we find a nonzero bias Δ v = - 66 ± 20 km s -1 for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at z = 1.30 producing the highest-redshift 21 cm intensity mapping measurement thus far.

79 ASTRONOMY AND ASTROPHYSICS↗

LENS: Learning Enabled Network Synthesis

RTRC and UMD have developed novel machine learning based methods under the ARPA-E DIFFERENTIATE program for rapid acceleration of hypothesis generation in complex architecture design spaces involving both discrete choices of component inclusion and interconnection and continuous parametric decisions. The project named Learning Enabled Network Synthesis (LENS) further demonstrated the developed methods on challenging electrical power converter design problems by identifying the most suitable circuit topologies and simultaneously selecting the most appropriate components to achieve optimized design of power converter with improved performances. We demonstrated that LENS could enable exploration of very large design space of circuit topologies and components by addressing the limitations of conventional design process in non-linear, high switching speed, multi-dimensional power converter design and optimization. The key innovation developed in LENS is the seamless integration of statistical learning and logical reasoning techniques and building on the individual strengths of these techniques for rapid hypothesis discovery. The main component of LENS comprises of: 1) Graph Reasoning Engine (GRE) to enforce composition rules that rapidly reject all discrete architectures that are composed incorrectly and generates an adaptive database of feasible designs which can be used by ML modules, 2) Graph Generative Learning module which is a deep neural network based generative model for graph architectures which can enable design space exploration beyond the dataset generated by the GRE, 3) Graph Reduced Order Model (ROM) for graph domains for accelerating computation of output metrics, and 4) Active learning and Rule Discovery module for sample efficient learning and extracting logical rules from the learned ML models which will be integrated in the GRE to enhance the filtering effectiveness. LENS approach can be applied to any design domains where designs can be represented as multi-attribute graphs. The LENS team integrated the various technical innovations listed above into an optimization pipeline and exercised the optimization pipeline on the converter design problem. The LENS project demonstrated that the developed AI/ML technologies can be used to generate novel converter circuits >45x faster than experts on chosen use-cases. This can enable faster design space exploration and identification of new designs which are not considered by experts due to the increasing design space complexity. This has significant potential impact on the public and energy needs of the country. It is currently estimated that 30% of all electrical powers generated passes through power converters. The future estimate is that 80% of all power generated would be passing through converters. LENS fills a critical gap in this space since by accelerating the design process the designers would be able to generate more efficient converters which can lead to significant energy savings for the country.

42 ENGINEERING↗

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗

Machine-Protection System for the PIP-II Linear Accelerator: Architecture and Data Processing Schemes

PIP-II at Fermilab features a brand-new, 800 MeV leading-edge Linear Accelerator (Linac) that will enable the Fermilab complex to deliver more than a megawatt of beam power to the Long Base-line Neutrino Facility (LBNF), requiring a robust Machine Protection System (MPS) to prevent beam-induced damage to delicate cryomodules and vacuum components. Signals in various accelerator components, in both analog and digital format, are collected and processed to assest the operating conditions of the accelerator. In the situation a potential damage is foreseen, the beam permit signal will be drop. A critical element of this system is the fast Analog Machine Protection System, a high-bandwidth platform designed for real-time beam loss monitoring. The system utilizes a modular, FPGA-based architecture to digitize signals from beam-sensing devices, such as AC Current Transformers (ACCTs), non-invasive Ring Pickups (RPUs) and beam scrapers. To eliminate both high frequency as well as 60 Hz noise that might be picked up in the industrial environment, data processing schemes beyond regular digital filtering are planned and implemented.

Wu, Jinyuan [Fermilab] (ORCID:0000000344329521)↗

Identification of engine oil-derived ash nanoparticles and ash formation process for a gasoline direct-injection engine

Engine oil-derived ash particles emitted from internal combustion (IC) engines are unwanted by-products, after oil is involved in in-cylinder combustion process. Since they typically come out together with particulate emissions, no detail has been reported about their early-stage particles other than agglomerated particles loaded on aftertreatment catalysts and filters. To better understand ash formation process during the combustion process, here differently formulated engine oils were dosed into a fuel system of a gasoline direct injection (GDI) engine that produces low soot mass emissions at normal operating conditions to increase the chances to find stand-alone ash particles separated from soot aggregates in the sub-20-nm size range. In addition to them, ash/soot aggregates in the larger size range were examined using scanning transmission electron microscopy (STEM)-X-ray electron dispersive spectroscopy (XEDS) to present elemental information at different sizes of particles from various oil formulations. The STEM-XEDS results showed that regardless of formulated oil type and particle size, Ca, P and C were always contained, while Zn was occasionally found on relatively large particles, suggesting that these elements get together from an early stage of particle formation. The S, Ca and P K-edge X-ray absorption near edge structure (XANES) analyses were performed for bulk soot containing raw ash. The linear combination approach & cross-checking among XANES results proposed that Ca 5 (OH)(PO 4 ) 2 , Ca 3 (PO 4 ) 2 and Zn 3 (PO 4 ) 2 are potentially major chemical compounds in raw ash particles, when combined with the STEM-XEDS results. Despite many reports that CaSO 4 is a major ash chemical when ash found in DPF/GFP systems was examined, it was observed to be rarely present in raw ashes using the S K-edge XANES analysis, suggesting ash transformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trace Chemical Detection Using Intercalated MXenes as a Signal Enhancing Substrate in Optical Probes

MXenes are 2D materials composed of layered transition metal nitrides or carbides. These materials are synthesized by HF exfoliation from MAX phases (Ti{sub 3}AlC{sub 2}). The 2D nanomaterial was synthesized by the removal of the 'A' element, resulting in a Mxene product (Ti{sub 3}C{sub 2}). MXenes have the general formula M{sub n+1}X{sub n}T{sub x}, where M is an early transition metal, X is Carbon and/or Nitrogen, and T accounts for surface terminated functional groups such as Fluoride, hydroxyl, and oxygen. These materials have very unique properties, similar to graphene, that allows them to be applied in a variety of trace detection techniques including surface-enhanced Raman spectroscopy (SERS). MXenes have also been demonstrated to selectively uptake uranyl ion, UO{sub 2}{sup 2+}. If this property can be combined with SERS or fluorescence detection, it may be possible to use MXenes as the basis for an alternative method to kinetic phosphorescence analysis (KPA) for trace uranyl measurements. Objectives: To confirm that MXene Nano materials are suitable substrates for SERS and sensor development by enhancing Raman signaling. To determine if certain MXene preparation methods yield materials that are more suitable for trace sensing methods. To determine uranyl uptake properties of these MXene materials and test them for analytical signals. Sample Preparation: Preparation of Ti{sub 3}C{sub 2}MXene (at FSU). MXenes were prepared by etching Al from Ti{sub 3}AlC{sub 2} (MAX phase)material. Two etching techniques yield different MXene products: LiF/HCl: Milder reaction, larger MXene flakes. HF: Harsher reaction, smaller flakes, larger layer separation. Products washed to remove etchant, vacuum filtered, and dried. Dried MXene flakes are air-stable. Film preparation for sensor testing (at SRNL): Suspend powder in diH{sub 2}O, purge with Ar, sonicate for 30 min. Centrifuge and collect supernate with suspended particles. Observed LiF-etched Mxene yielded a higher density of particles and darker collected solution. Drop-cast (4 ml) supernate onto slides and dried with Ar. For Rhodamine B (RhB) testing, drop-cast 4 ml drops onto Mxene spots and dried with Ar. Scanning Electron Microscopy conditions: 10 kV Beam energy, high vacuum; Working distance of 8 mm; beam penetration depth appx. 4 microns, beam spot size appx. 2 nanometers. Results: Detection of aluminum correlates with bright spots on image. Presence of aluminum shows that LiF/HCl etching was less thorough than HF etching. Trace Cl detection in LiF images suggests incomplete rinsing. HF has smaller feature size, more layer structure, and increased homogeneity, consistent with expectations. Macroscopic Raman spectroscopy measurements: 532 nm excitation, ∼50 mW with a ∼100 micron spot size (InPhotonics RPB probe). Kaiser Optical Holospec f/1.8 spectrometer with cooled (-60 deg.C) Andor iDus OE420 CCD. LiF 1x supernate showed good signal for trace measurements of Rhodamine B. HF and 1/4x LiF supernates showed little Mxene or Rhodamine B signal. Low deposition densities led to excess background signal from glass slides. For LiF film, response is linear with Rhodamine B concentration over range tested. Will retest with Raman microscope (∼1 micron spot size) to characterize SERS of more dilute LiF and HF etched Mxenes. Conclusions: The LiF etched material was more suitable for macroscopic SERS measurements because it was more concentrated, resulting in a thicker film than the HF etched Mxene and diluted LiF sample. However, the other materials may give greater SERS enhancements, which we hope to determine from measurements with the Raman microscope. From characterization with SEM we concluded that the HF etched Mxene is more uniform/homogenous and has smaller particle size than the LiF etched Mxene. There is still aluminum present in both samples indicating that etching wasn't complete, but the removal of the aluminum was more efficient in the HF method. Path Forward: Observe SERS with Raman microscopy, to obtain better signals for the more diluted samples and be able to compare enhancement effects for the different MXenes. Characterize uranyl sorption into MXene films and test Raman and fluorescence signals. Revisit the etching conditions to improve removal of aluminum. FSU and SRNL will continue to collaborate to create and characterize different Mxene materials and test their usefulness for sensor applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗

Experimental Safe Extremum Seeking for Accelerators

Here, we demonstrate the recent designs of safe extremum seeking (Safe ES) on the 1-km-long charged particle accelerator at the Los Alamos Neutron Science Center (LANSCE). Safe ES is a modification of extremum seeking (ES) which, in addition to minimizing an analytically unknown cost, also employs a safety filter based on an analytically unknown control barrier function (CBF) safety metric. Tuning is necessitated by accelerators being large complex systems, with many drifting parameters due to thermal effects and degradation. At the same time, safe operation (the maintenance of state constraints) is crucial, as damage brings astronomical costs, both financially and in operation downtime. Our measured (but analytically unknown) safety metric is the beam current. We perform multivariable Safe ES on three accelerator applications, in which we adapt 4, 6, and 3 magnet strength parameters, respectively. Two of the three applications are for validated simulation models of beamlines at LANSCE: the first for the proton radiography (pRad) beamline of 800-MeV protons for spot size tuning; the second on a high-performance code, HPSim, for tuning the low-energy beam transport (LEBT) region that contains a beam of 750-keV protons. The third is an experimental tuning of the steering magnets in the LEBT at LANSCE.

43 PARTICLE ACCELERATORS↗

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

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↗

Using the Gerchberg-Saxton algorithm to reconstruct nonmodulated pyramid wavefront sensor measurements

Adaptive optics (AO) is a technique to improve the resolution of ground-based telescopes by correcting, in real-time, optical aberrations due to atmospheric turbulence and the telescope itself. With the rise of Giant Segmented Mirror Telescopes (GSMT), AO is needed more than ever to reach the full potential of these future observatories. One of the main performance drivers of an AO system is the wavefront sensing operation, consisting of measuring the shape of the above mentioned optical aberrations. Aims. The non-modulated pyramid wavefront sensor (nPWFS) is a wavefront sensor with high sensitivity, allowing the limits of AO systems to be pushed. The high sensitivity comes at the expense of its dynamic range, which makes it a highly non-linear sensor. We propose here a novel way to invert nPWFS signals by using the principle of reciprocity of light propagation and the Gerchberg-Saxton (GS) algorithm. We test the performance of this reconstructor in two steps: the technique is first implemented in simulations, where some of its basic properties are studied. Then, the GS reconstructor is tested on the Santa Cruz Extreme Adaptive optics Laboratory (SEAL) testbed located at the University of California Santa Cruz. This new way to invert the nPWFS measurements allows us to drastically increase the dynamic range of the reconstruction for the nPWFS, pushing the dynamics close to a modulated PWFS. The reconstructor is an iterative algorithm requiring heavy computational burden, which could be an issue for real-time purposes in its current implementation. However, this new reconstructor could still be helpful in the case of many wavefront control operations. In conclusion, this reconstruction technique has also been successfully tested on the Santa Cruz Extreme AO Laboratory (SEAL) bench where it is now used as the standard way to invert nPWFS signal.

47 OTHER INSTRUMENTATION↗