Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle filtering”

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 55 records · Page 3

Combining Remotely Sensed Evapotranspiration and an Agroecosystem Model to Estimate Center‐Pivot Irrigation Water Use at High Spatio‐Temporal Resolution

Abstract Estimating irrigation water use accurately is critical for sustainable irrigation and studying terrestrial water cycle in irrigated croplands. However, irrigation is not monitored in most places, and current estimations of irrigation water use has coarse spatial and/or temporal resolutions. This study aims to estimate irrigation water use at the daily and field scale through the proposed model‐data fusion framework, which is achieved by particle filtering with two configurations (concurrent, CON, and sequential, SEQ) by assimilating satellite‐based evapotranspiration (ET) observations into an advanced agroecosystem model, ecosys . Two types of experiments using synthetic and real ET observations were conducted to study the efficacy of the proposed framework for estimating irrigation water use at the irrigated fields in eastern and western Nebraska, United States. The experiments using synthetic ET observations indicated that, for two major sources of uncertainties of ET difference between observations and model simulations, which are bias and noise, noise had larger impacts on degrading the estimation performance of irrigation water use than bias. For the experiments using real ET observations, monthly and annual estimations of irrigation water use matched well with farmer irrigation records, with Pearson correlation coefficient ( r ) around 0.80 and 0.50, respectively. Although detecting daily irrigation records was very challenging, our method still gave a good performance with RMSE, BIAS, and r around 2.90, 0.03, and 0.4 mm/d, respectively. Our proposed model‐data fusion framework for estimating irrigation water use at high spatio‐temporal resolution could contribute to regional water management, sustainable irrigation, and better tracking terrestrial water cycle.

Environmental Sciences & Ecology↗

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 ↗

Knowing our Sun: data fusion for optimizing space weather forecasts

Outline: Space weather: model for predicting solar wind & polarity; Particle filtering: optimization with Monte Carlo; Simulation & observation: twin tests and real data; Back to Earth: independent work in inertial confinement fusion (ICF)

79 ASTRONOMY AND ASTROPHYSICS↗

Big Data For Operation and Maintenance Cost Reduction

The purpose of this research is to develop a first-of-a-kind framework for integrating Big Data capability into the daily activities of our current fleet of nuclear power plants. Big Data is traditionally defined as data sets with high volume, velocity, and heterogeneity, and the existing Big Data analytics capabilities are now widely popular in fields such as finance, weather, e-commerce, healthcare and sports. In the nuclear industry, while the volume and velocity of data may present computational challenges for existing analytics capabilities, data heterogeneity are seen to present the major challenge. This research project mainly focuses on incorporating the wide range of data heterogeneities in nuclear power plants into an integrated Big Data Analytics capability. The primary end-product of this project is a Big Data framework that is capable of dealing with the large volume and heterogeneity of the data found in nuclear power plants to extract timely and valuable information on equipment performance. The framework can generate system insights that are actionable relations between measurable impacts and the corresponding maintenance action plans and enable optimization of plant operation and maintenance based on the extracted information. The developed framework is capable of handling heterogeneous data including both image data and time-series sensor data. Specifically, this developed framework includes the following components. The first component is an overarching maintenance ontology which includes system insights required by maintenance optimization. The maintenance ontology interacts with other components in the developed framework. The second component handles Piping & Instrumentation Diagram (P&ID) data. It can be used to extract system components and their relations automatically from the P&IDs. This extracted information is stored in the first component, i.e., maintenance ontology, and is also used as input to the third component, i.e., a tool for generating the fault tree for the corresponding system. The generated fault tree in turn is stored in the ontology for assessing risk that is used as a criterion in maintenance policy optimization. The fourth component is a tool for inferring the parameters in the Markov degradation model for a nuclear system. It uses basic information from the ontology. The fifth component is a tool for assessing the degradation level using sensor measurement data, for example, pressure, flowrate. This tool can be used for determining corrective maintenance actions. The results obtained from components four and five are returned to the ontology. The sixth component of the framework is a tool for optimizing the maintenance policy for a nuclear system of interest. It takes certain basic information from the ontology, e.g., costs of maintenance actions and system failures, as input, and returns the optimal maintenance policy to the ontology. This tool can be used for determining predictive maintenance actions. A set of experiments have also been conducted to verify the algorithms developed in this project for nuclear system degradation monitoring. The experiments are based on four solenoid valves, similar to the ones used in nuclear power plants. The analyses based on the experimental data using two algorithms, i.e., the Randomized Window Decomposition (RWD) algorithm and the particle filtering algorithm, and the results are introduced in the report. The Big Data framework developed in this project can be used as a support tool in daily activities of plant operation and maintenance and will reduce current costs while maintaining or improving safety levels. Overall, the project will not only benefit existing reactors, however it will open new frontiers to realize the long overdue value of Big Data Analytics in the nuclear sphere.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗

Detecting Low-level Radiation Sources Using Border Monitoring Gamma Sensors

We consider a problem of detecting a low-level radiation source using a network of Gamma spectral sensors placed on the periphery of a monitored region. We propose a computationally light-weight, correlation-based method which is primarily intended for systems with limited computing capacity. Sensor measurements are combined at the fusion by first generating decisions at each time step and then taking their majority vote within a time widow. At each time step, decisions are generated using two strategies: (i) SUM method based on a threshold decision on a correlation statistic derived from measurements from all sensors, and (ii) OR method based on logical-OR of threshold decisions based on correlations statistics of individual sensor measurements. We derive analytical performance bounds for false alarm rates of SUM and OR methods, and show that their performance is enhanced by the temporal smoothing of majority vote within a time window. Using measurements from a test campaign, we generate a border monitoring scenario with twelve 2"x2" NaI Gamma sensors deployed on the periphery of 42m x 42m outdoor region. A Cs-137 source is moved in a straight-line across this region, starting several meters outside and finally moving away from it. We illustrate the performance of both correlation-based detection methods, and compare their performances with each other and with a particle filter method. Overall, under small false-alarm conditions, the OR fusion is found to produce better detection performance.

Sen, Satyabrata↗

Particle trapping, size-filtering, and focusing in the nonthermal plasma synthesis of sub-10 nanometer particles

Low-pressure nonthermal flowing plasmas are widely used for the gas-phase synthesis of nanoparticles and quantum dots of materials that are difficult or impractical to synthesize using other techniques. To date, the impact of temporary electrostatic particle trapping in these plasmas has not been recognized, a process that may be leveraged to control particle properties. Furthermore, we present experimental and computational evidence that, during their growth in the plasma, sub-10 nm silicon particles become temporarily confined in an electrostatic trap in radio-frequency excited plasmas until they grow to a size at which the increasing drag force imparted by the flowing gas entrains the particles, carrying them out of the trap. We demonstrate that this trapping enables the size filtering of the synthesized particles, leading to highly monodisperse particle sizes, as well as the electrostatic focusing of the particles onto the reactor centerline. Understanding of the mechanisms and utilization of such particle trapping will enable the design of plasma processes with improved size control and the ability to grow heterostructured nanoparticles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-Dimensional Analysis of Particle Distribution on Filter Layers inside N95 Respirators by Deep Learning

The global COVID-19 pandemic has changed many aspects of daily lives. Wearing personal protective equipment, especially respirators (face masks), has become common for both the public and medical professionals, proving to be effective in preventing spread of the virus. Nevertheless, a detailed understanding of respirator filtration-layer internal structures and their physical configurations is lacking. Here, we report three-dimensional (3D) internal analysis of N95 filtration layers via X-ray tomography. Using deep learning methods, we uncover how the distribution and diameters of fibers within these layers directly affect contaminant particle filtration. The average porosity of the filter layers is found to be 89.1%. Contaminants are more efficiently captured by denser fiber regions, with fibers <1.8 μm in diameter being particularly effective, presumably because of the stronger electric field gradient on smaller diameter fibers. This study provides critical information for further development of N95-type respirators that combine high efficiency with good breathability.

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↗

Characterization, Chemistry, and Particle Size Distribution of Fallout Particles Isolated from Filter Samples

Nuclear particulate fallout is the radioactive byproduct of a nuclear event formed by the mixture of proximate environmental materials with vaporized bomb debris. The debris can be transported into the atmosphere during cloud rise, raining out locally and globally constituting a radiation hazard. Questions remain on how entrainment of environmental material in the fireball affects fallout formation processes and radionuclide incorporation during cooling.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Convergence Analysis for an Online Data-Driven Feedback Control Algorithm

This paper presents convergence analysis of a novel data-driven feedback control algorithm designed for generating online controls based on partial noisy observational data. The algorithm comprises a particle filter-enabled state estimation component, estimating the controlled system’s state via indirect observations, alongside an efficient stochastic maximum principle-type optimal control solver. By integrating weak convergence techniques for the particle filter with convergence analysis for the stochastic maximum principle control solver, we derive a weak convergence result for the optimization procedure in search of optimal data-driven feedback control. Numerical experiments are performed to validate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Observations of supermicron-sized aerosols originating from biomass burning in southern Central Africa

During the 3 years of the ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) campaign, the NASA Orion P-3 was equipped with a 2D stereo (2D-S) probe that imaged particles with maximum dimension (D) ranging from 10 < D < 1280 µm. The 2D-S recorded supermicron-sized aerosol particles (SAPs) outside of clouds within biomass burning plumes during flights over the southeastern Atlantic off Africa's coast. Numerous SAPs with 10 < D < 1520 µm were observed in 2017 and 2018 at altitudes between 1230 and 4000 m, 1000 km from the coastline, mostly between 7–11° S. No SAPs were observed in 2016 as flights were conducted further south and further from the coastline. Number concentrations of refractory black carbon (rBC) measured by a single particle soot photometer ranged from 200 to 1200 cm –3 when SAPs were observed. Transmission electron microscopy images of submicron particulates, collected on Holey carbon grid filters, revealed particles with potassium salts, black carbon (BC), and organics. Energy-dispersive X-ray spectroscopy spectra also detected potassium, a tracer for biomass burning. These measurements provided evidence that the submicron particles originated from biomass burning. NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) 3 d back trajectories show a source in northern Angola for times when large SAPs were observed. Fire Information for Resource Management System (FIRMS) Moderate Resolution Imaging Spectroradiometer (MODIS) 6 active fire maps showed extensive biomass burning at these locations. Given the back trajectories, the high number concentrations of rBC, and the presence of elemental tracers indicative of biomass burning, it is hypothesized that the SAPs imaged by the 2D-S are examples of BC aerosol, ash, or unburned plant material.

54 ENVIRONMENTAL SCIENCES↗

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗

Laser-driven ionization mechanisms of aluminum for single particle aerosol mass spectrometry

Single particle aerosol mass spectrometry (SPAMS), an analytical technique for measuring the size and composition of individual micron-scale particles, is capable of analyzing atmospheric pollutants and bioaerosols much more efficiently and with more detail than conventional methods which require the collection of particles onto filters for analysis in the laboratory. Despite SPAMS’ demonstrated capabilities, the primary mechanisms of ionization are not fully understood, which creates challenges in optimizing and interpreting SPAMS signals. In this paper, we present a well-stirred reactor model for the reactions involved with the laser-induced vaporization and ionization of an individual particle. The SPAMS conditions modeled in this paper include a 248 nm laser which is pulsed for 8 ns to vaporize and ionize each particle in vacuum. The ionization of 1 μm, spherical Al particles was studied by approximating them with a 0-dimensional plasma chemistry model. The primary mechanism of absorption of the 248 nm photons was pressure-broadened direct photoexcitation to Al(y 2 D). Atoms in this highly excited state then undergo superelastic collisions with electrons, heating the electrons and populating the lower energy excited states. We found that the primary ionization mechanism is electron impact ionization of various excited state Al atoms, especially Al(y 2 D). Because the gas expands rapidly into vacuum, its temperature decreases rapidly. The rate of three-body recombination (e - + e - + Al + → Al + e - ) increases at low temperature, and most of the electrons and ions produced recombine within several μs of the laser pulse. The importance of the direct photoexcitation indicates that the relative peak heights of different elements in SPAMS mass spectra may be sensitive to the available photoexcitation transitions. We also discuss the effects of laser intensity, particle diameter, and expansion dynamics.

47 OTHER INSTRUMENTATION↗

A copula-based rank histogram ensemble filter

Serial ensemble filters implement triangular probability transport maps to reduce high-dimensional inference problems to sequences of state-by-state univariate inference problems. The univariate inference problems are solved by sampling posterior probability densities obtained by combining constructed prior densities with observational likelihoods according to Bayes' rule. Many serial filters in the literature focus on representing the marginal posterior densities of each state. However, rigorously capturing the conditional dependencies between the different univariate inferences is crucial to correctly sampling multidimensional posteriors. This work proposes a new serial ensemble filter, called the copula rank histogram filter (CoRHF), that seeks to capture the conditional dependency structure between variables via empirical copula estimates; these estimates are used to rigorously implement the triangular (state-by-state univariate) Bayesian inference. The success of the CoRHF is demonstrated on two-dimensional examples and the Lorenz'63 problem. A practical extension to the high-dimensional setting is developed by localizing the empirical copula estimation, and is demonstrated on the Lorenz'96 problem.

97 MATHEMATICS AND COMPUTING↗