Search NASA⌕ Search

SEARCH · Search NASA

Results for “Normal Distribution”

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

A Statistical Simulation Approach to Safe Life Fatigue Analysis of Redundant Metallic Components

This paper introduces a dual active load path fail-safe fatigue design concept analyzed by Monte Carlo simulation. The concept utilizes the inherent fatigue life differences between selected pairs of components for an active dual path system, enhanced by a stress level bias in one component. The design is applied to a baseline design; a safe life fatigue problem studied in an American Helicopter Society (AHS) round robin. The dual active path design is compared with a two-element standby fail-safe system and the baseline design for life at specified reliability levels and weight. The sensitivity of life estimates for both the baseline and fail-safe designs was examined by considering normal and Weibull distribution laws and coefficient of variation levels. Results showed that the biased dual path system lifetimes, for both the first element failure and residual life, were much greater than for standby systems. The sensitivity of the residual life-weight relationship was not excessive at reliability levels up to R = 0.9999 and the weight penalty was small. The sensitivity of life estimates increases dramatically at higher reliability levels.

Matthews, William T.↗
NASA NTRS1960

Experimental Smoke and Electromagnetic Analog Study of Induced Flow Field About a Model Rotor in Steady Flight Within Ground Effect

Hovering and steady low-speed forward-flight tests were run on a 4-foot-diameter rotor at a ground height of 1 rotor radius. The two blades had a 2 to 1 taper ratio and were mounted in a see-saw hub. The solidity ratio was 0.05. Measurements were made of the rotor rpm, collective pitch, and forward-flight velocity. Smoke was introduced into the tip vortex and the resulting vortex pattern was photographed from two positions. Using the data obtained from these photographs, wire models of the tip vortex configurations were constructed and the distribution of the normal component of induced velocity at the blade feathering axis that is associated with these tip vortex configurations was experimentally determined at 450 increments in azimuth position from this electromagnetic analog. Three steady-state conditions were analyzed. The first was hovering flight; the second, a flight velocity just under the wake "tuck under" speed; and the third, a flight velocity just above this speed. These corresponded to advance ratios of 0, 0.022, and 0.030 (or ratios of forward velocity to calculated hovering induced velocity of approximately 0, 0.48, and 0.65), respectively, for the model test rotor. Cross sections of the wake at 450 intervals in azimuth angle as determined from the path of the tip vortex are presented graphically for all three cases. The nondimensional normal component of the induced velocity that is associated with the tip vortex as determined by an electromagnetic analog at 450 increments in azimuth position and at the blade feathering axis is presented graphically. It is shown that the mean value of this component of the induced velocity is appreciably less after tuck-under than before. It is concluded that this method yields results of engineering accuracy and is a very useful means of studying vortex fields.

Gray, Robin B.↗
NASA NTRS2000

Rolling Bearing Life Prediction-Past, Present, and Future

Comparisons were made between the life prediction formulas of Lundberg and Palmgren, Ioannides and Harris, and Zaretsky and full-scale ball and roller bearing life data. The effect of Weibull slope on bearing life prediction was determined. Life factors are proposed to adjust the respective life formulas to the normalized statistical life distribution of each bearing type. The Lundberg-Palmgren method resulted in the most conservative life predictions compared to Ioannides and Harris, and Zaretsky methods which produced statistically similar results. Roller profile can have significant effects on bearing life prediction results. Roller edge loading can reduce life by as much as 98 percent. The resultant predicted life not only depends on the life equation used but on the Weibull slope assumed, the least variation occurring with the Zaretsky equation. The load-life exponent p of 10/3 used in the American National Standards Institute (ANSI)/American Bearing Manufacturers Association (ABMA)/International Organization for Standardization (ISO) standards is inconsistent with the majority roller bearings designed and used today.

Zaretsky, E V↗
NASA NTRS1997

Effect of age and gender on sudomotor and cardiovagal function and blood pressure response to tilt in normal subjects

Normative data are limited on autonomic function tests, especially beyond age 60 years. We therefore evaluated these tests in a total of 557 normal subjects evenly distributed by age and gender from 10 to 83 years. Heart rate (HR) response to deep breathing fell with increasing age. Valsalva ratio varied with both age and gender. QSART (quantitative sudomotor axon-reflex test) volume was consistently greater in men (approximately double) and progressively declined with age for all three lower extremity sites but not the forearm site. Orthostatic blood pressure reduction was greater with increasing age. HR at rest was significantly higher in women, and the increment with head-up tilt fell with increasing age. For no tests did we find a regression to zero, and some tests seem to level off with increasing age, indicating that diagnosis of autonomic failure was possible to over 80 years of age.

Non-NASA Center↗
NASA NTRS2008

Tolerance and UQ4SIM: Nimble Uncertainty Documentation and Analysis Software

Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and variabilities is a necessary first step toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. The basic premise of uncertainty markup is to craft a tolerance and tagging mini-language that offers a natural, unobtrusive presentation and does not depend on parsing each type of input file format. Each file is marked up with tolerances and optionally, associated tags that serve to label the parameters and their uncertainties. The evolution of such a language, often called a Domain Specific Language or DSL, is given in [1], but in final form it parallels tolerances specified on an engineering drawing, e.g., 1 +/- 0.5, 5 +/- 10%, 2 +/- 10 where % signifies percent and o signifies order of magnitude. Tags, necessary for error propagation, can be added by placing a quotation-mark-delimited tag after the tolerance, e.g., 0.7 +/- 20% 'T_effective'. In addition, tolerances might have different underlying distributions, e.g., Uniform, Normal, or Triangular, or the tolerances may merely be intervals due to lack of knowledge (uncertainty). Finally, to address pragmatic considerations such as older models that require specific number-field formats, C-style format specifiers can be appended to the tolerance like so, 1.35 +/- 10U_3.2f. As an example of use, consider figure 1, where a chemical reaction input file is has been marked up to include tolerances and tags per table 1. Not only does the technique provide a natural method of specifying tolerances, but it also servers as in situ documentation of model uncertainties. This tolerance language comes with a utility to strip the tolerances (and tags), to provide a path to the nominal model parameter file. And, as shown in [1], having the ability to quickly mark and identify model parameter uncertainties facilitates error propagation, which in turn yield output uncertainties.

Kleb, Bil↗
NASA NTRS2011

Development of the ISS EMU SPEEDR

The Self Powered EVA EMU Data Recorder (SPEEDR) is an FPGA (Field-programmable gate array) based device designed to collect high-rate EMU (Extravehicular Mobility Unit) PLSS (Primary Life Support Subsystem) data for download at a later time. The existing EMU PLSS data down-link capability during EVA is one data packet every 2 minutes and is subject to bad packets or loss of signal. High-rate PLSS data is generated by the ECWS (Enhanced Caution and Warning System) but is not normally captured or distributed. Access to high-rate data will increase the capability of EMU anomaly resolution team to pinpoint issues remotely, saving crew time by reducing required call-down Q&A and on-orbit diagnostic activities. With no Shuttle flights post FY11, and potentially limited down-mass capability, the ISS crew and ground support personnel will have to be capable of on-orbit operations to maintain, diagnose, repair, and return to service EMU hardware, possibly through 2028. Collecting high-rate EMU PLSS data during both IVA (Intravehicular Activity) and EVA (Extravehicular Activity) operations will provide trending analysis for life extension and/or predictive performance. The SPEEDR concept has generated interest as a tool/technology that could be used for other ISS subsystems or future exploration-class space suits where hardware reliability/availability is critical and low/variable bandwidth may require "store then forward" methodology. Preliminary work in FY11 produced a functional prototype consisting of an FPGA evaluation board, custom memory/interface circuit board, and custom software. The SPEEDR concept includes a stand-alone battery that is recharged by a computer USB (Universal Serial Bus) port while data is being downloaded.

Bernard. Craig↗
NASA NTRS2012

Development of the Self-Powered Extravehicular Mobility Unit Extravehicular Activity Data Recorder

The Self-Powered Extravehicular Mobility Unit (EMU) Extravehicular Activity (EVA) Data Recorder (SPEEDR) is a field-programmable gate array (FPGA)-based device designed to collect high-rate EMU Primary Life Support Subsystem (PLSS) data for download at a later time. During EVA, the existing EMU PLSS data downlink capability is one data packet every 2 minutes and is subject to bad packets or loss of signal. Higher-rate PLSS data is generated by the Enhanced Caution and Warning System but is not normally captured or distributed. Access to higher-rate data will increase the capability of EMU anomaly resolution team to pinpoint issues remotely, saving crew time by reducing required call-down Q&A and on-orbit diagnostic activities. With no Space Shuttle flights post Fiscal Year 2011 (FY11), and potentially limited down-mass capability, the ISS crew and ground support personnel will have to be capable of on-orbit operations to maintain, diagnose, repair, and return to service EMU hardware, possibly through 2028. Collecting high-rate EMU PLSS data during both intravehicular activity (IVA) and EVA operations will provide trending analysis for life extension and/or predictive performance. The SPEEDR concept has generated interest as a tool/technology that could be used for other International Space Station subsystems or future exploration-class space suits where hardware reliability/availability is critical and low/variable bandwidth may require store then forward methodology. Preliminary work in FY11 produced a functional prototype consisting of an FPGA evaluation board, custom memory/interface circuit board, and custom software. The SPEEDR concept includes a stand-alone battery that is recharged by a computer Universal Serial Bus (USB) port while data are being downloaded.

Bernard, Craig↗
NASA NTRS2012

Coherent Detector for Near-Angle Scattering and Polarization Characterization of Telescope Mirror Coatings

A report discusses the difficulty of measuring scattering properties of coated mirrors extremely close to the specular reflection peak. A prototype Optical Hetero dyne Near-angle Scatterometer (OHNS) was developed. Light from a long-coherence-length (>150 m) 532-nm laser is split into two arms. Acousto-optic modulators frequency shift the sample and reference beams, establishing a fixed beat frequency between the beams. The sample beam is directed at very high f/# onto a mirror sample, and the point spread function (PSF) formed after the mirror sample is scanned with a pinhole. This light is recombined by a non-polarizing beam splitter and measured through heterodyne detection with a spectrum analyzer. Polarizers control the illuminated and analyzed polarization states, allowing the polarization dependent scatter to be measured. The bidirectional reflective or scattering distribution function is normally measured through use of a scattering goniometer instrument. The instrumental beam width (collection angle span) over which the scatterometer responds is typically many degrees. The OHNS enables measurement at angles as small as the first Airy disk diameter.

Macenka, Steven A.↗
NASA NTRS2012

Recent Upgrades to NASA SPoRT Initialization Datasets for the Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed several products for its National Weather Service (NWS) partners that can initialize specific fields for local model runs within the NOAA/NWS Science and Training Resource Center (STRC) Environmental Modeling System (EMS). In last year's NWA abstract on this topic, the suite of SPoRT products supported in the STRC EMS was presented, which includes a Sea Surface Temperature (SST) composite, a Great Lakes sea-ice extent, a Green Vegetation Fraction (GVF) composite, and NASA Land Information System (LIS) gridded output. This abstract and companion presentation describes recent upgrades made to the SST and GVF composites, as well as the real-time LIS runs. The Great Lakes sea-ice product is unchanged from 2011. The SPoRT SST composite product has been expanded geographically and as a result, the resolution has been coarsened from 1 km to 2 km to accommodate the larger domain. The expanded domain covers much of the northern hemisphere from eastern Asia to western Europe (0 N to 80 N latitude and 150 E to 10 E longitude). In addition, the NESDIS POES-GOES product was added to fill in gaps caused by the Moderate Resolution Imaging Spectroradiometer (MODIS) being unable to sense in cloudy regions, replacing the recently-lost Advanced Microwave Scanning Radiometer for EOS with negligible change to product fidelity. The SST product now runs twice per day for Terra and Aqua combined data collections from 0000 to 1200 UTC and from 1200 to 0000 UTC, with valid analysis times at 0600 and 1800 UTC. The twice-daily compositing technique reduces the overall latency of the previous version while still representing the diurnal cycle characteristics. The SST composites are available at approximately four hours after the end of each collection period (i.e. 1600 UTC for the nighttime analysis and 0400 UTC for the daytime analysis). The real-time MODIS GVF composite has only received minor updates in the past year. The domain was expanded slightly to extend further west, north, and east to improve coverage over parts of southern Canada. Minor adjustments were also made to the manner in which GVF is calculated from the distribution of maximum Normalized Difference Vegetation Index from MODIS. The presentation will highlight some examples of the substantial inter-annual change in GVF that occurred from 2010 to 2011 in the U.S. Southern Plains as a result of the summer 2011 drought, and the early vegetation green up across the eastern U.S. due to the very warm conditions in March 2012. Finally, the SPoRT LIS runs the operational Noah land surface model (LSM) in real time over much of the eastern half of the CONUS. The Noah LSM is continually cycled in real time, uncoupled to any model, and driven by operational atmospheric analyses over a long-term, multi-year integration. The LIS-Noah provides the STRC EMS with high-resolution (3 km) LSM initialization data that are in equilibrium with the operational analysis forcing. The Noah LSM within the SPoRT LIS has been upgraded from version 2.7.1 to version 3.2, which has improved look-up table attributes for several land surface quantities. The surface albedo field is now being adjusted based on the input real-time MODIS GVF, thereby improving the net radiation. Also, the LIS-Noah now uses the newer MODIS-based land use classification scheme (i.e. the International Biosphere-Geosphere Programme [IGBP]) that has a better depiction of urban corridors in areas where urban sprawl has occurred. STRC EMS users interested in initializing their LSM fields with high-resolution SPoRT LIS data should set up their model domain with the MODIS-IGBP 20-class land use database and select Noah as the LSM.

Case, Jonathan L.↗
NASA NTRS2017

An Application Programming Interface for Synthetic Snowflake Particle Structure and Scattering Data

The work by Kuo and colleagues on growing synthetic snowflakes and calculating their single-scattering properties has demonstrated great potential to improve the retrievals of snowfall. To grant colleagues flexible and targeted access to their large collection of sizes and shapes at fifteen (15) microwave frequencies, we have developed a web-based Application Programming Interface (API) integrated with NASA Goddard's Precipitation Processing System (PPS) Group. It is our hope that the API will enable convenient programmatic utilization of the database. To help users better understand the API's capabilities, we have developed an interactive web interface called the OpenSSP API Query Builder, which implements an intuitive system of mechanisms for selecting shapes, sizes, and frequencies to generate queries, with which the API can then extract and return data from the database. The Query Builder also allows for the specification of normalized particle size distributions by setting pertinent parameters, with which the API can also return mean geometric and scattering properties for each size bin. Additionally, the Query Builder interface enables downloading of raw scattering and particle structure data packages. This presentation will describe some of the challenges and successes associated with developing such an API. Examples of its usage will be shown both through downloading output and pulling it into a spreadsheet, as well as querying the API programmatically and working with the output in code.

reflectivity properties of snowflakes↗
NASA NTRSDate not supplied

Lower neck injury criteria for THOR and Hybrid III dummies in rear impact

The objectives of the study are to derive lower-neck-injury probability curves under rear impact loading from matched pair Hybrid III and THOR dummy tests. Methods: Twelve whole-body and 15 isolated head-neck rear-impact sled tests were conducted using the 2 dummies. They were positioned on a rigid seat that was attached to an acceleration sled. The dummies were positioned with the head parallel to the ground, torso against the seat back, and legs stretched such that there was no axial rotation. The acceleration pulse matched the previous in-house human cadaver tests. The 6-axis lower-neck load cell was used in both dummies. For the isolated rear-impact tests,the head-neck was excised from the dummies, and the lower-neck load cell was mounted to the top of the sled with the head parallel to the ground and such that the acceleration vector was in the rear impact mode. Lower-neck loads and lower-neck-injury criteria (LNij) were obtained using the load cell data and survival analysis to develop injury-assessment-risk values and curves for both dummies. The LNij criteria were determined for both dummies using the intercept value corresponding to the 90% probability level for the forces and moments. Results: The log normal and Weibull distributions were the optimal distributions for the Hybrid III and THOR devices. At the 50% risk level, the mean LNij of 1.1 and NCIS was 0.77 for the Hybrid III, and 1.5 NCIS was 0.30 for the THOR device. The quality indices were in the fair range and good range for the 2 dummies, respectively, at this risk level. Conclusions: The lower-neck-based IARCs and IARVs serve as the first dataset for injury assessments, and the THOR may be a better test device for assessing injures in rear impact environments.

Narayan Yoganandan↗
NASA NTRS1994

Stagnation Region Heat Transfer: The Influence of Turbulence Parameters, Reynolds Number and Body Shape

The effect of velocity gradient on stagnation region heat transfer augmentation by free stream turbulence was investigated. Heat transfer was measured in the stagnation region of four models with elliptical leading edges with ratios of major to minor axes of 1:1, 1.5:1, 2.25:1, and 3:1. Four geometrically similar, square bar, square mesh, biplane grids were used to generate free stream turbulence with different intensities and length. Heat transfer measurements were made for the following ranges of parameters: Reynolds number, based on leading edge diameter, 37,000 to 228,000; dimensionless leading edge velocity gradient, 1.20 to 1.80; turbulence intensity, 1.1 to 15.9%; and length scale to leading edge diameter ratio, 0.05 to 0.30. Stagnation point heat transfer augmentation by free stream turbulence can be predicted using a modified version of a previously developed correlation for a circular leading edge. Heat transfer augmentation was independent of body shape at the stagnation point. The heat transfer distribution down-stream from the stagnation point can be predicted using the normalized laminar heat transfer distribution.

Vanfossen, G. James↗
NASA NTRS1976

Radiance, polarization, and ellipticity of the radiation in the earth's atmosphere

The complete radiation field including polarization is calculated for a model of the real atmosphere by the matrix operator method. The radiance, direction and amount of polarization, and ellipticity are obtained at the top and bottom of the atmosphere for three values of the surface albedo (0; 0.15 0.90) and five solar zenith angles. Scattering and absorption by molecules (including ozone) and by aerosols are taken into account together with the variation of the number density of these substances with height. All results are calculated for both a normal aerosol number and a distribution which is one-third of the normal amount at all heights. The calculated values show general qualitative agreement with the available experimental measurements. The position of the neutral points of the polarization in the principal plane is a sensitive indicator of the characteristics of the aerosol particles in the atmosphere, since it depends on the sign and value of the single scattered polarization for scattering angles around 20 deg and 160 deg for transmitted and reflected photons respectively.

Hitzfelder, S. J.↗
NASA NTRS1978

Coronal gas in the Galaxy. II - A statistical analysis of O VI absorptions

This paper deals with general inferences about the low-density phase of interstellar gas having temperatures well above 200,000 K (the 'coronal gas') which can be drawn from O VI absorption data for 72 stars. Attention is given to the behavior of radial velocities, possible evidence for circumstellar O VI, the space distribution of the O VI gas, crowding of normal interstellar gas, and temperature distributions for the coronal gas. A model is adopted in which the coronal gas is contained within randomly distributed and nonoverlapping parcels, each with a size, pressure, and internal temperature distribution that do not vary markedly from one unit to the next. It is shown that the one-dimensional velocity dispersion for O VI regions (26 km/s) is substantially higher than the value for ordinary interstellar clouds (6.4 km/s).

Jenkins, E. B.↗
NASA NTRS1992

A statistical model for evaluating GOPEX uplink performance

This article describes a statistical model to analyze the signal intensity received at the solid-state imaging (SSI) camera of the Galileo optical communications system from an Earth-based transmitter (GOPEX) demonstration. The analytical model assumes that the optical beam possesses a Gaussian profile and the communication channel has a log-normal scattering characteristic. The atmospheric-induced jitter is modelled as two independent zero mean Gaussian random variables. By modelling the system parameters as a set of independent and identically distributed (iid) random variables, the combined impact of uncertainties due to system parameters and the turbulent atmosphere can be approximated by a log-normal distributed signal intensity at the spacecraft. A Monte-Carlo software simulation package was also developed to compute the confidence interval probabilities for general optical beam profiles. Numerical results show that the approximation is valid for a wide range of operation scenarios.

Kiesaleh, K.↗
NASA NTRS2007

Self-consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves. 2. Wave Induced Ring Current Precipitation and Thermal Electron Heating

This paper continues presentation and discussion of the results from our new global self-consistent theoretical model of interacting ring current ions and propagating electromagnetic ion cyclotron waves [Khazanov et al., 2006]. To study the effects of electromagnetic ion cyclotron wave propagation and refraction on the wave induced ring current precipitation and heating of the thermal plasmaspheric electrons, we simulate the May 1998 storm. The main findings after a simulation can be summarized as follows. Firstly, the wave induced ring current precipitation exhibits quite a lot of fine structure, and is highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 x 10(exp 6) (cm(raised dot) s(raised dot) sr(raised dot) (sup -1)) are observed during the maill and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not connected to the most intense waves in simple manner. The characteristics of the wave power spectral density distribution over the wave normal angle are extremely crucial for the effectiveness of the ring current ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from RAM [Kozyra et al., 1997a] reveals that although we observe a qualitative agreement between the localizations of the wave induced precipitations in the models, there is no quantitative agreement between the magnitudes of the fluxes. The quantitative differences are mainly due to a qualitative difference between the characteristics of the wave power spectral density distributions over the wave normal angle in RAM and in our model. Thirdly, the heat fluxes to plasmaspheric electrons caused by Landau resonate energy absorption from electromagnetic ion cyclotron waves are observed in the postnoon-premidnight MLT sector, and can reach the magnitude of 10(exp 11) eV/(cm(sup 2)(raised dot)s). The Coulomb energy degradation of the RC H(+) and O(+) ions maximizes at about 10(exp 11) (eV/(cm(sup 2) (raised dot) s), and typically leads to electron energy deposition rates of about 2(raised dot) 10(exp 10) (eV/(cm(sup 2)(raised dot)s) which are observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the ring current with the plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, ring current, and topside ionosphere [Gurgiolo et al., 2005]. Finally, the wave induced intense electron heating has a structure of the spot-like patches along the most enhanced density gradients in the plasmasphere boundary layer and can be a possible driver to the observed but still not explained small-scale structures of enhanced emissions in the stable auroral red arcs.

Khazanov, G. V.↗
NASA NTRS2022

Structural and chemical complexity of minerals: An update

The complexities of chemical composition and crystal structure are fundamental characteristics of minerals that have high relevance to the understanding of their stability, occurrence and evolution. This review summarises recent developments in the field of mineral complexity and outlines possible directions for its future elaboration. The database of structural and chemical complexity parameters of minerals is updated by H-correction of structures with unknown H positions and the inclusion of new data. The revised average complexity values (arithmetic means) for all minerals are 3.54(2) bits/atom and 345(10) bits/cell (based upon 4443 structure reports). The distributions of atomic information amounts, chemIG and strIG, versus the number of mineral species fit the normal modes, whereas the distributions of total complexities, chemIG,total and strIG,total, along with numbers of atoms per formula and per unit cell are log normal. The three most complex mineral species known today are ewingite, morrisonite and ilmajokite, all either discovered or structurally characterised within the last five years. The most important complexity-generating mechanisms in minerals are: (1) the presence of isolated large clusters; (2) the presence of large clusters linked together to form three-dimensional frameworks; (3) formation of complex three-dimensional modular frameworks; (4) formation of complex modular layers; (5) high hydration state in salts with complex heteropolyhedral units; and (6) formation of ordered superstructures of relatively simple structure types. The relations between symmetry and complexity are considered. The analysis of temporal dynamics of mineralogical discoveries since 1875 with the step of 25 years show the increasing chemical and structural complexities of human knowledge of the mineral kingdom in the history of mineralogy. In the Earth’s history, both diversity and complexity of minerals experience dramatic increases associated with the formation of Earth’s continental crust, initiation of plate tectonics and the Great Oxidation event.

complexity↗
DOE OSTI2026

Machine Learning-Based Anomaly Detection for PMT Data Quality Monitoring in the SBN and DUNE

Maintaining high-quality detector data is essential for achieving the scientific objectives of the Short-Baseline Neutrino (SBN) Program at Fermilab. Current data quality monitoring (DQM) procedures rely primarily on threshold-based metrics and manual inspection of detector monitoring plots, making the detection of subtle or gradually developing anomalies both time-consuming and dependent on expert interpretation. This project developed and evaluated a machine-learning workflow for automatically identifying anomalous photomultiplier tube (PMT) channels in the Short-Baseline Near Detector (SBND) using optical-hit amplitude data. A Python-based analysis program was developed to process ROOT files, extract statistical features describing individual PMT amplitude distributions, and generate feature vectors for anomaly detection. These features were used to train an Isolation Forest model using data representing normal detector operation. The trained model was subsequently applied to independent detector runs to identify channels exhibiting statistically unusual behavior relative to the learned reference response. To support expert interpretation, the workflow generated complementary diagnostic products, including anomaly score distributions, normalized amplitude comparisons, decision-tree visualizations, and principal component analysis (PCA) projections. This project demonstrated the feasibility of integrating unsupervised machine learning into detector data-quality monitoring and developed a complete workflow for automated PMT performance assessment to aid expert-driven review. Beyond its technical contributions, the VFP appointment fostered a research collaboration between Aurora University and Fermilab and provided direct workforce development benefits by training the visiting faculty member in detector-scale machine-learning methods that are now being incorporated into undergraduate coursework and research. The methodology developed here provides a foundation for future applications to ProtoDUNE and other liquid argon time projection chamber (LArTPC) detectors, contributing to ongoing efforts to improve detector reliability, reduce manual monitoring requirements, and enable scalable data quality monitoring for future large-scale neutrino experiments, including the Deep Underground Neutrino Experiment (DUNE).

Colón Santana, Juan A. [Unlisted, US, IL]↗