Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle filter”

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 19 records

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗

Particle Filter Based Inference Testing

The primary intent of PAR-FIT (Particle Filter based Inference Testing) is to provide hard inductive evidence that a machine learning model is capable and proven for an individual test input. By examining training data used to form the underlying model functional correlation, an estimate of the reliability that a model will make the correct prediction can be made. The Sequential Probability Ratio Test is used to derive a qualitative evaluation for reliability based on hypothesis testing. The PAR-FIT framework achieves this by implementing a particle filter and the sequential probability ratio test algorithms on the machine learning model training data to determine relevancy of new individual test samples to the training dataset. The kernel function evaluates the local proximity and density of training data used to derive a prediction outcome. Particles are used to probabilistically determine which training data to evaluate for proximity. For test samples that are within a close proximity to and surrounded by multiple training data points, the evaluated reliability of the prediction is high. For test samples that are anomalies not represented by the training dataset, in low density data clusters, or are far from existing data points, the evaluated reliability is low as insufficient training evidence exists to suggest the model is capable of making the correct prediction. Sequential Probability Ratio Test is further used to determine when a hypothesis on whether a signal can be rejected or accepted for use. The ratio test collects sequence information from the particle filter to test whether the signal is anomalous or normal via hypothesis testing of the underlying distributions.

Chen, Edward [Idaho National Laboratory (INL), Ida↗

Applying the Kalman filter particle method to strange and open charm hadron reconstruction in the STAR experiment

Here, we apply KF Particle, a Kalman Filter package for secondary vertex finding and fitting, to strange and open charm hadron reconstruction in heavy-ion collisions in the STAR experiment. Compared to the conventional helix swimming method used in STAR, the KF Particle method improves the reconstructed Λ, Ω and $D$ 0 significance considerably. At the same time, we demonstrate that Monte Carlo simulation with the STAR detector responses can well reproduce the topological variable distributions reconstructed in real data using the KF Particle method, therefore retaining good control on the reconstruction efficiency uncertainties for strange and open charm hadrons measurements in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Ensemble Score Filter for Tracking High-Dimensional Nonlinear Dynamical Systems

We propose an ensemble score filter (EnSF) for solving high-dimensional nonlinear filtering problems with superior accuracy. A major drawback of existing filtering methods, e.g., particle filters or ensemble Kalman filters, is the low accuracy in handling high-dimensional and highly nonlinear problems. EnSF attacks this challenge by exploiting the score-based diffusion model, defined in a pseudo-temporal domain, to characterizing the evolution of the filtering density. EnSF stores the information of the recursively updated filtering density function in the score function, in stead of storing the information in a set of finite Monte Carlo samples (used in particle filters and ensemble Kalman filters). Unlike existing diffusion models that train neural networks to approximate the score function, we develop a training-free score estimation that uses mini-batch-based Monte Carlo estimator to directly approximate the score function at any pseudo-spatial-temporal location, which provides sufficient accuracy in solving high-dimensional nonlinear problems as well as saves tremendous amount of time spent on training neural networks. Another essential aspect of EnSF is its analytical update step, gradually incorporating data information into the score function, which is crucial in mitigating the degeneracy issue faced when dealing with very high-dimensional nonlinear filtering problems. High-dimensional Lorenz systems are used to demonstrate the performance of our method. EnSF provides surprisingly impressive performance in reliably tracking extremely high-dimensional Lorenz systems (up to 1,000,000 dimension) with highly nonlinear observation processes, which is a well-known challenging problem for existing filtering methods.

97 MATHEMATICS AND COMPUTING↗

Minimum entropy filtering for a single output non-Gaussian stochastic system using state transformation

This paper presents a novel filter design for the single-output stochastic non-linear systems subjected to non-Gaussian noises and the proposed assumptions. Based on a state transformation, the unmeasurable states of the systems can be estimated where non-linear terms in the systems have been eliminated. It has been shown that the estimation error is linearly dynamical regarding to the presented vector-valued filter gain which can be optimised by minimising the entropy-based performance criterion. In addition, the convergence of the presented algorithm is analysed in mean-square sense and a numerical example is given to verify the effectiveness of the presented filtering algorithm. Meanwhile, the extended Kalman filter, unscented particle filter and minimum entropy filter are given for the comparisons of the filtering performance. Following the presented framework, some extensions of the presented filtering algorithm are discussed to indicate the flexibility of the filter design. The contribution of this paper can be summarised as establishing a novel minimum entropy filtering framework which consists of model transformation, entropy optimisation and convergence analysis.

42 ENGINEERING↗

In-pixel integration of signal processing and AI/ML based data filtering for particle tracking detectors

We present the first physical realization of in-pixel signal processing with integrated AI-based data filtering for particle tracking detectors. Building on prior work that demonstrated a physics-motivated edge-AI algorithm suitable for ASIC implementation, this work marks a significant milestone toward intelligent silicon trackers. Our prototype readout chip performs real-time data reduction at the sensor level while meeting stringent requirements on power, area, and latency. The chip is taped-out in 28nm TSMC CMOS bulk process, which has been shown to have sufficient radiation hardness for particle experiments. This development represents a key step toward enabling fully on-detector edge AI, with broad implications for data throughput and discovery potential in high-rate, high-radiation environments such as the High-Luminosity LHC.

Parpillon, Benjamin [Fermilab; Illinois U., Chicag↗

Real-Time Radiological Source Term Estimation for Multiple Sources in Cluttered Environments

A particle filter algorithm is presented to estimate the position, strength, and cardinality of an unknown number of radioactive point sources in an obstacle-rich environment using count measurements. The algorithm addresses gaps in the prior literature by incorporating two novel elements. The first is a precomputation step in which local terrain and obstacle data is processed to compute attenuation kernels throughout the search area. This enables rapid estimation performance in obstacle-rich environments as measurements are gathered. The second novel feature is a dynamic particle allocation technique in which the number of particles is adjusted in real time to meet convergence goals. This feature allows the algorithm to scale more efficiently to scenarios with a larger number of sources. Furthermore, a series of computational experiments using simulated data demonstrates the algorithm’s performance in a cluttered environment with up to eight sources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Studies of optical and chemical properties of aged and fresh biomass burning absorbing aerosols for climate models

Particulate emissions from burning fuels from sub-Saharan Africa is currently poorly characterized. With increasing population and industrialization, emissions are expected to grow. Optical methods for measuring light absorption in the ultraviolet (UV) and visible (vis) regions of the spectrum are crucial for understanding the optical properties of BB aerosols (Bond and Bergstrom, 2006; Lack and Cappa, 2010; Moschos et al., 2021). One such instrument is the dual-spot AE33 model aethalometer is a common filter-based instrument that can be used to measure aerosol absorption properties (Hansen et al., 1984). This instrument has many figures of merit; however, the method's limitations include potential biases from scattering within the filter material, multiple scattering effects of the particles, filter loading (i.e. particle shadowing), and the need for wavelength-specific correction factors (Cλ), which depend on accurate intercomparison with other in situ techniques (Bond et al., 1999; Collaud Coen et al., 2010). This work aims to: 1. Optically characterize BB aerosol generated from African fuel 2. Correct the aethalometer at multiple wavelengths 3. Parametrize these correction factors. 4. Combusting the same fuels with two difference tube furnace systems assess the general applicability of this flexible means of producing BB aerosol. This was accomplished by conducting laboratory studies at Los Alamos National Laboratory (LANL) and North Carolina A&T State University (NCA&T), where aethalometer measurements were compared to in situ optical measurements. Mass cross sections of BB aerosol from these under-studied fuels will be determined along with their wavelength dependence.

09 BIOMASS FUELS↗

Research to Address Technical Barriers to Expanded Markets for Biodiesel and Biodiesel Blends (CRADA Final Report)

NREL and the Clean Fuels Alliance America (CFAA) will work cooperatively to assess the effects of biodiesel blends on the performance of modern diesel engines. This work will include research to understand the impact of biodiesel blends on the operation and durability of particle filters and NO x control sorbents/catalysts, to quantify the effect on emission control systems performance, and to understand effects on engine component durability. This research was performed at NREL Renewable Fuels and Lubricants (REFUEL) laboratory (an engine testing laboratory) as well as at third party labs paid directly by CFAA with NREL as part of the project management team. Also, research to develop appropriate ASTM standards for biodiesel quality and stability was conducted in NREL’s bench scale fuel chemistry laboratory. The cooperative project involved laboratory testing and research at NREL using biodiesel from a variety of sources and in collaboration with a broad range of other stakeholders. In addition, NREL will work with NBB to set up an Industrial Steering Committee to design the scope for the various tasks and to provide technical oversight to these projects. NREL and NBB will cooperatively communicate the study results to as broad an audience as possible. This research benefits the public by expanding markets for a domestically produced low-carbon intensity fuel for use in diesel engines.

09 BIOMASS FUELS↗

Orbitrap LC-MS Analysis of Nanoparticle Composition at the EPCAPE Mount Soledad site between 04 18 2023 and 06 14 2023

Weekly peak lists containing m/z, intensity, and assigned formula for filter samples, size selected for sub-100 nm particles. Filters were collected daily between 4/18/23 and 6/14/23, grouped based on calendar week for extraction, and analyzed via Thermo Scientific Q Exactive Plus Orbitrap LC-MS. Formulas were assigned to background-corrected peak lists and restricted to CHONS/CHONSNa atoms for the negative and positive modes respectively. Filters were grouped into calendar weeks 0-8 with dates provided in README text file.

54 ENVIRONMENTAL SCIENCES↗

Aerosol Microphysics and Chemical Measurements at Mt. Soledad and Scripps Pier during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) from February 2023 to February 2024 UCSD Library Collection

This dataset includes guest instrument measurements and other PI products for aerosol microphysics and chemical measurements collected at Mt. Soledad and Scripps Pier during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) from February 2023 to February 2024. The measurements include the following instruments at Mt. Soledad: High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS, Aerodyne), Scanning Electrical Mobility Spectrometer (SEMS, Brechtel Manufacturing Inc.), Aerodynamic Particle Sizer (APS, Droplet Measurements Technologies), Single Particle Soot Photometer (SP2, Drople Measurements Technologies), Meteorological Station (WXT520, Vaisala), Ozone (Teco), and trace gas proxies (Teledyne). In addition, the analyses of particle filters collected at Mt. Soledad for three dry-diameter size cuts (<1 micron, <0.5 micron, <0.18 micron) and at Scripps Pier for one dry-diameter size cut (<1 micron) by Fourier Transform Infrared (FTIR) and X-ray Fluorescence (XRF) are reported. A differential mobility analyzer operated as a scanning mobility particle sizer (SMPS, TSI Inc.), a printed particle optical spectrometer (POPS, Grimm), and a continuous flow diffusion cloud condensation nuclei (CCN, DMT) counter provide the mobility aerosol size distribution (30-360 nm), optical size distribution (150 - 6000 nm), size-resolved CCN distribution (30-360 nm) at 0.2, 0.4, 0.6, 0.8, and 1.0% supersaturation. Measurements are reported for both sampling from an isokinetic aerosol inlet and from a Counterflow Virtual Impactor (CVI, Brechtel Manufacturing Inc.). The data are available at the following link: https://library.ucsd.edu/dc/collection/bb0898306q

54 ENVIRONMENTAL SCIENCES↗

Proactive Frequency Stability Scheme via Bayesian Filters and Synchrophasors

Underfrequency (UF) load shedding schemes are traditionally implemented in two ways: One approach is based on manual load shedding, with system operators requesting loads to be shed ahead of anticipated stressful operating conditions. Manual load shedding is usually done through phone calls. The second method is automatic load shedding via underfrequency relays. Using static static settings, these schemes can be designed to operate in stages and drop previously identified loads. The main limitation of traditional load shedding schemes is that they are reactive and leave little room for optimized corrective actions. This work presents a proactive and automatic underfrequency load shedding solution for power systems. Measurements are captured via phasor measurement units (PMUs) at relatively low sampling rates of 30 Hz. These measurements are then processed by particle filters who predict the future state of the system's frequency. Based on these predictions excess load is determined and shed. Comparative case studies are performed in simulated environments. Easy-to-implement models, without hard-to-derive parameters, highlight potential aspects for real-life implementation.

Paramo, Gian↗

Microgrid Frequency Stability: A Proactive Scheme Based on Dynamic Predictions

The dynamic nature of microgrids introduces challenges in the context of frequency stability. This work presents a framework where the future state of microgrid frequency is predicted and corrective actions are optimized. Predictions are generated through Bayesian filters leveraging synchronized data acquired via PMUs. Taking a proactive approach makes it possible to optimize corrective actions considering dynamic system conditions. Testing is conducted via Matlab simulations. The performance of the solution presented in this work is compared to traditional load-shedding schemes and to predictive solutions found in literature. The results indicate that the proposed framework outperforms both. Some of the advantages of this framework include a reduction in amount of load dropped during compensation, the use of adaptive parameters which eliminates the need to simulate contingency conditions, and dynamic uncertainty quantification provided the particle filter.

Paramo, Gian↗

Flow Induced Vibration Studies in Pressurized Helium Gas Cooling Channels

Production of metastable Technetium-99 (Tc-99m), a radioactive tracer that emits gamma rays, is vital to the medical imaging community. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. The work presented in this report is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 through the irradiation of Mo-100 targets using an electron beam. The NMR target comprises a stack of approximately 60-70 Mo-100 disks with diameter 24 mm and thickness 0.74 mm held in stainless steel laminations, each separated using 0.25 mm thick stainless steel spacers. The symmetric target stack is housed in an Inconel vessel with two Inconel windows on either side. Two electron accelerators are used to produce 40 MeV, 3.16 µA electron beams each that penetrate the Inconel windows and irradiate the Mo-100 disks. Approximately 90% of the total 250 kW beam power is deposited in the NMR target during the irradiation process. During irradiation, pressurized helium gas flows through 0.25 mm thin gaps between the disks cooling the beam window, target disks, disk laminations and spacers. Both NMR and LANL have found during cold testing of the target system (no heat deposition) that the Mo-100 disks undergo significant mass loss and disk breakage due to vibrations induced by the flowing helium gas. The mass loss is not only undesirable due to monetary loss from reduced final quantities of Mo-99, but also due to the hazards associated with radioactive material trapped in the cooling lines and particle filters. The effect of flow rate and target geometry on the flow induced vibrations need to be quantified, and recommendations provided to minimize this mass loss. This work describes LANL’s experimental characterization of the flow induced vibrations and disk mass loss in a reduced scale set-up containing 10 Mo-100 disks. We use high speed imaging, displacement measurements and microphone measurements combined with signal processing to estimate the vibration frequency of each disk. The effect of disk thickness, target fit and duration of testing on the mass loss is described. We find that in the current configuration of NMR targets, the vibrations and mass loss on the first disk are minimized, while those in the adjacent disks are highest. The microphone and high-speed image data show that increased flow rates and increased duration of testing increases vibration frequency and mass loss. The mass loss is due to both disk rotation and back and forth motion. There are visible wear marks on the disks with the highest mass loss. We also note that the current NMR window gap reduces flow induced vibrations compared to the previous smaller gaps. Longer duration testing will provide more data and verification for the findings presented in this report. The work will be continued in FY 24.

43 PARTICLE ACCELERATORS↗

Multiple Changepoint Detection for Non‐Gaussian Time Series

ABSTRACT This article combines methods from existing techniques to identify multiple changepoints in non‐Gaussian autocorrelated time series. A transformation is used to convert a Gaussian series into a non‐Gaussian series, enabling penalized likelihood methods to handle non‐Gaussian scenarios. When the marginal distribution of the data is continuous, the methods essentially reduce to the change of variables formula for probability densities. When the marginal distribution is count‐oriented, Hermite expansions and particle filtering techniques are used to quantify the scenario. Simulations demonstrating the efficacy of the methods are given and two data sets are analyzed: 1) the proportion of home runs hit by Major League Baseball batters from 1920 to 2023 and 2) a six‐dimensional series of tropical cyclone counts from the Earth's basins of generation from 1980 to 2023. In the first series, beta marginal distributions are used to describe the proportions; in the second, Poisson marginal distributions seem appropriate.

Lund, Robert [Department of Statistics University ↗

Window Cooling Studies and Disk Vibration Testing on a Subset of Mo-100 Disks

Production of metastable Technetium-99 (Tc-99m), a radioactive tracer that emits gamma rays, is vital to the medical imaging community. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. The work presented in this report is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 through the irradiation of Mo-100 targets using an electron beam. The NMR target comprises a stack of approximately 60-70 Mo-100 disks with diameter 24 mm and thickness 0.74 mm held in stainless steel laminations, each separated using 0.25 mm thick stainless-steel spacers. The symmetric target stack is housed in an Inconel vessel with two Inconel windows on either side. Two electron accelerators are used to produce 40 MeV, 3.16 µA electron beams each that penetrate the Inconel windows and irradiate the Mo-100 disks. Approximately 90% of the total 250 kW beam power is deposited in the NMR target during the irradiation process, with a smaller percentage adding up to 2.2 kW of heat deposited on the Inconel window. During irradiation, pressurized helium gas flows through thin gaps between the disks cooling the beam window, target disks, disk laminations and spacers. Both NMR and LANL have found during cold testing of the target system (no heat deposition) that the Mo 100 disks undergo significant mass loss and disk breakage due to vibrations induced by the flowing helium gas. The mass loss is not only undesirable due to monetary loss from reduced final quantities of Mo-99, but also due to the hazards associated with radioactive material trapped in the cooling lines and particle filters. The effect of flow rate and target geometry on the flow induced vibrations need to be quantified, and recommendations provided to minimize this mass loss. LANL has previously also tested NMR’s Inconel beam window by heating the window, while flowing pressurized helium, using the average heat deposited on the window. However, the NMR beam is pulsed with a duty cycle of 12.5%, which introduces oscillation in temperature around the nominal 600 °C steady state value with each pulse. Available fatigue curves for Inconel are few, established for room temperature, and they are based on mechanical strain cycles not thermally induced strain as in the NMR target. The effect of pulsed beam heating on the Inconel window therefore needs to be quantified. This report details the experiments conducted to assess the factors that lead to mass loss in the NMR target disks as well as to understand the effect of a pulsed beam on NMR’s Inconel window. This work describes LANL’s experimental characterization of the flow induced vibrations and disk mass loss in a reduced scale set-up containing 5 to 10 Mo-100 disks. We use high speed imaging, displacement measurements and microphone measurements combined with signal processing to estimate the vibration frequency of each disk. The effect of disk thickness, target fit and duration of testing on the mass loss is described. We find that in the current configuration of NMR targets, the vibrations and mass loss on the first disk are minimized, while those in the adjacent disks are highest. The microphone and high-speed image data show that increased flow rates and increased duration of testing increases vibration frequency and mass loss. The mass loss is due to both disk rotation and back and forth motion. There are visible wear marks on the disks with the highest mass loss. We also note that the current NMR window gap reduces flow induced vibrations compared to the previous smaller gaps. Improved target holders significantly reduce disk mass loss to almost negligible quantities. This work finds that the larger window to first disk gap and improved target holder geometry should allow NMR to successfully conduct irradiations with minimal mass loss. The window tests were conducted to understand the effect of a pulsed beam on both the window longevity and to estimate the window temperature and displacement during pulsing. The experiments presented here were performed at significantly low power, due to the limitations of the induction heating system. The window temperature rose to approximately 73 °C with a significantly reduced power of 45 W without beam pulsing. With a 5 Hz pulse rate, 12.5% duty cycle, the window temperature remained constant at 26 °C. These experiments will be repeated with improved coil geometry and reported in upcoming journal papers.

42 ENGINEERING↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗