Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Statistical Aspects of North Atlantic Basin Tropical Cyclones During the Weather Satellite Era, 1960-2013

This Technical Publication (TP) is part 2 of a two-part study of the North Atlantic basin tropical cyclones that occurred during the weather satellite era, 1960-2013. In particular, this TP examines the inferred statistical relationships between 25 tropical cyclone parameters and 9 specific climate-related factors, including the (1) Oceanic Niño Index (ONI), (2) Southern Oscillation Index (SOI), (3) Atlantic Multidecadal Oscillation (AMO) index, (4) Quasi-Biennial Oscillation (QBO) index, (5) North Atlantic Oscillation (NAO) index of the Climate Prediction Center (CPC), (6) NAO index of the Climate Research Unit (CRU), (7) Armagh surface air temperature (ASAT), (8) Global Land-Ocean Temperature Index (GLOTI), and (9) Mauna Loa carbon dioxide (CO2) (MLCO2) index. Part 1 of this two-part study examined the statistical aspects of the 25 tropical cyclone parameters (e.g., frequencies, peak wind speed (PWS), accumulated cyclone energy (ACE), etc.) and provided the results of statistical testing (i.e., runs-testing, the t-statistic for independent samples, and Poisson distributions). Also, the study gave predictions for the frequencies of the number of tropical cyclones (NTC), number of hurricanes (NH), number of major hurricanes (NMH), and number of United States land-falling hurricanes (NUSLFH) expected for the 2014 season, based on the statistics of the overall interval 1960-2013, the subinterval 1995-2013, and whether the year 2014 would be either an El Niño year (ENY) or a non-El Niño year (NENY).

Wilson, Robert M.↗
NASA NTRS2016

On the Calculation of Uncertainty Statistics with Error Bounds for CFD Calculations Containing Random Parameters and Fields

This chapter discusses the ongoing development of combined uncertainty and error bound estimates for computational fluid dynamics (CFD) calculations subject to imposed random parameters and random fields. An objective of this work is the construction of computable error bound formulas for output uncertainty statistics that guide CFD practitioners in systematically determining how accurately CFD realizations should be approximated and how accurately uncertainty statistics should be approximated for output quantities of interest. Formal error bounds formulas for moment statistics that properly account for the presence of numerical errors in CFD calculations and numerical quadrature errors in the calculation of moment statistics have been previously presented in [8]. In this past work, hierarchical node-nested dense and sparse tensor product quadratures are used to calculate moment statistics integrals. In the present work, a framework has been developed that exploits the hierarchical structure of these quadratures in order to simplify the calculation of an estimate of the quadrature error needed in error bound formulas. When signed estimates of realization error are available, this signed error may also be used to estimate output quantity of interest probability densities as a means to assess the impact of realization error on these density estimates. Numerical results are presented for CFD problems with uncertainty to demonstrate the capabilities of this framework.

Uncertainty↗
NASA NTRS2016

Comparison of Atmospheric Delay Statistics from Deep Space Network Arrays and Nearby Test Interferometers

Several techniques have been explored and demonstrated that allow for greater data return on space-to-ground links. Among these techniques, arraying several smaller diameter dish antennas together is one method used in several arenas. These arrays can achieve larger effective area and gain than are available from a single larger antenna. This technique is routinely used by the NASA Deep Space Network (DSN) at 8.4 GHz where the incoming signals are much weaker than those experienced by the near-Earth satellite community. When considering arraying at much higher frequencies such as 32 GHz deep-space Ka-band, the phase alignment of the individual antenna signals is significantly disrupted by atmospheric turbulence. Since 2012, several downlink array demonstrations have been conducted using 32 GHz carrier signals emitted by the deep space probes Cassini and Kepler. Site test interferometers (STIs) that receive signals from geostationary satellites have been deployed at all three DSN tracking complexes for long-term monitoring of atmospheric delay fluctuations. In a previous DSN array demonstration study involving the Cassini spacecraft, it was shown that statistics of the adjusted STI phase fluctuations matched the statistics of concurrent array demonstration phase fluctuations. These adjustments accounted for differences in antenna separation, elevation angle and spacecraft frequencies. The STI antenna separations were about 200 m and the DSN antenna separations were about 300 m. These adjustments made use of the thick-layer turbulence model that was applicable to the Goldstone desert climate during the summer months for which the data were acquired. In this paper, we report on the results of additional array demonstrations involving the Kepler spacecraft and compare the adjusted STI phase fluctuations with those seen by a nearby two-element array of 34 m diameter antennas tracking Kepler’s 32 GHz signal at the Goldstone, California and Madrid, Spain DSN sites. We also discuss results from a demonstration using an array over a longer 12.5 km baseline. The Cassini and Kepler array demonstrations were found to validate the long term statistics acquired from several years of STI data as well as the models used to adjust the statistics for the conditions of an array. These statistics represent reliable estimates of the phase fluctuations that would be seen by an array tracking a deep space signal after applying appropriate adjustments for a given array configuration, elevation angle profile and observing frequency.

Morabito, David D.↗
NASA NTRS2023

Overview and Statistical Analysis of Boundary Layer Clouds and Precipitation Over the Western North Atlantic Ocean

Due to their fast evolution and large natural variability in macro- and microphysical properties, the accurate representation of boundary layer clouds in current climate models remains a challenge. One of the regions with large intermodel spread in the Coupled Model Intercomparison Project Phase 6 ensemble is the western North Atlantic Ocean. Here, statistically representative in situ measurements can help to develop and constrain the parameterization of clouds in global models. To this end, we performed comprehensive measurements of boundary layer clouds, aerosol, trace gases, and radiation in the western North Atlantic Ocean during the NASA Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) mission. In total, 174 research flights with 574 flight hours for cloud and precipitation measurements were performed with the HU-25 Falcon during three winter (February–March 2020, January–April 2021, and November 2021–March 2022) and three summer seasons (August–September 2020, May–June 2021, and May–June 2022). Here we present a statistical evaluation of 16 140 individual cloud events probed by the fast cloud droplet probe and the two-dimensional stereo cloud probe during 155 research flights in a representative and repetitive flight strategy allowing for robust statistical data analyses. We show that the vertical profiles of distributions of the liquid water content and the cloud droplet effective diameter (ED) increase with altitude in the marine boundary layer. Due to higher updraft speeds, higher cloud droplet number concentrations (N liquid ) were measured in winter compared to summer despite lower cloud condensation nucleus abundance. Flight cloud cover derived from statistical analysis of in situ data is reduced in summer and shows large variability. This seasonal contrast in cloud coverage is consistent with a dominance of a synoptic pattern in winter that favors conditions for the formation of stratiform clouds at the western edge of cyclones (post-cyclonic). In contrast, a dominant summer anticyclone is concomitant with the occurrence of shallow cumulus clouds and lower cloud coverage. The evaluation of boundary layer clouds and precipitation in the N liquid ED phase space sheds light on liquid, mixed-phase, and ice cloud properties and helps to categorize the cloud data. Ice and liquid precipitation, often masked in cloud statistics by a high abundance of liquid clouds, is often observed throughout the cloud. The ACTIVATE in situ cloud measurements provide a wealth of cloud information useful for assessing airborne and satellite remote-sensing products, for global climate and weather model evaluations, and for dedicated process studies that address precipitation and aerosol–cloud interactions.

Aerosol↗
DOE OSTI2013

Magnification Bias in Gravitational Arc Statistics

The statistics of gravitational arcs in galaxy clusters is a powerful probe of cluster structure and may provide complementary cosmological constraints. Despite recent progresses, discrepancies still remain among modelling and observations of arc abundance, specially regarding the redshift distribution of strong lensing clusters. Besides, fast "semi-analytic" methods still have to incorporate the success obtained with simulations. In this paper we discuss the contribution of the magnification in gravitational arc statistics. Although lensing conserves surface brightness, the magnification increases the signal-to-noise ratio of the arcs, enhancing their detectability. We present an approach to include this and other observational effects in semi-analytic calculations for arc statistics. The cross section for arc formation ({\sigma}) is computed through a semi-analytic method based on the ratio of the eigenvalues of the magnification tensor. Using this approach we obtained the scaling of {\sigma} with respect to the magnification, and other parameters, allowing for a fast computation of the cross section. We apply this method to evaluate the expected number of arcs per cluster using an elliptical Navarro--Frenk--White matter distribution. Our results show that the magnification has a strong effect on the arc abundance, enhancing the fraction of arcs, moving the peak of the arc fraction to higher redshifts, and softening its decrease at high redshifts. We argue that the effect of magnification should be included in arc statistics modelling and that it could help to reconcile arcs statistics predictions with the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Statistical Approaches for Pebble Bed Reactor Operations and Safeguards

The design of pebble bed reactors (PBRs) and their method of operation align more closely with statistical approaches used in manufacturing and process control than traditional safeguards statistical approaches. The reason is PBRs will employ a nondestructive assay (NDA) measurement (burnup measurement system [BUMS]) that is part of the fuel handling system supporting discharge decisions in addition to the reactor code which monitor performance. The integration of these two approaches provides the opportunity to monitor reactor performance statistically for both operations and safeguards in ways not achievable using other reactor designs. For light water reactors (LWRs), knowledge about reactor code performance in predicting irradiated special nuclear material (SNM) content historically was only achieved from special measurement campaigns or from fuel reprocessing. Conversely, through statistical comparison of the BUMS with the reactor code-predicted values, PBRs can achieve this in real time. The resulting SNM distribution is an indicator of reactor performance because factors such as transit time and path of the pebble fuel through the reactor determine the plutonium production and uranium depletion. By analyzing the predicted and measured values, opportunities exist to adjust operating parameters, fuel design, and other characteristics to optimize performance and fuel utilization. From a safeguards perspective, this approach also provides the information necessary to validate declared values and evaluate whether the reactor is being operated as expected. This paper outlines statistical approaches for PBRs that can be used to support both operations and safeguards.

Ball, Cory↗
DOE OSTI2016

Statistical and sensitivity analysis of failing rods in EPR LBLOCA

A statistical fuel failure analysis procedure has been developed, and is now applied in full scale to a large break loss-of-coolant accident (LBLOCA) scenario in EPR-type nuclear power plant. The goal is to statistically determine the number of failing fuel rods. It is also important to bring out the underlying causes for rod failures, and therefore a sensitivity analysis procedure for the preceding complex calculation chain is outlined and adopted. In order to produce a statistically reliable estimation on the number of failing rods, a large number of simulations are required. In the analysis, single rods are simulated with the coupled fuel performance - thermal hydraulics code FRAPTRAN-GENFLO. The statistically varied factors are divided into global and local by their range of influence. Transient thermal hydraulic and power history boundary conditions and model parameters are global and affect all the rods. Fuel rod manufacturing parameters and other rod position related factors are local and are related to a certain rod. The number of global scenarios is 59 as dictated by the Wilks' formula, and 1000 randomly sampled rods are simulated in each scenario. The thermal hydraulic and transient power boundary conditions for the fuel performance code are calculated with a statistical version of the system code APROS. The steady-state irradiation histories are simulated with FRAPCON, and the steady-state power histories used as boundary conditions for FRAPCON are obtained from SIMULATE calculations. Thus, a multistage calculation chain is required to consolidate the procedure. As an outcome, in the worst global scenario, 1.4% of the simulated rods failed. It can be concluded that the Finnish safety regulations, i.e. max. 10% of the rods allowed to fail, are met. In order to find out which input parameters have significance to the outcome of the analysis, a sensitivity analysis is done. At the moment, only the worst global scenario is considered but later on it is possible to include also the sensitivities of the global factors. Due to complexity of the existing data, first the relevant input parameters for the sensitivity analysis have to be specified. Data visualization with a cobweb graph is used for the screening. Then, selected sensitivity measures are calculated between the chosen input and output parameters. The sensitivity indices calculated are the Borgonovo's delta measure, the first order Sobol sensitivity index, and squared Pearson correlation coefficients. The first mentioned is a novelty in this context. As an outcome, the most relevant parameters with respect to the cladding integrity were determined to be the decay heat power during the transient, the thermal hydraulic conditions in the rod's location in the reactor, and the steady-state irradiation history of the rod as represented in this analysis by the rod burnup. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Memory access statistics monitoring

Systems, apparatuses, and methods related to memory access statistics monitoring are described. A host is configured to map pages of memory for applications to a number of memory devices coupled thereto. A first memory device comprises a monitoring component configured to monitor access statistics of pages of memory mapped to the first memory device. A second memory device does not include a monitoring component capable of monitoring access statistics of pages of memory mapped thereto. The host is configured to map a portion of pages of memory for an application to the first memory device in order to obtain access statistics corresponding to the portion of pages of memory upon execution of the application despite there being space available on the second memory device and adjust mappings of the pages of memory for the application based on the obtained access statistics corresponding to the portion of pages.

Roberts, David A.↗
DOE OSTI2023

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Statistic inversion of multi-zone transition probability models for aquifer characterization in alluvial fans

Understanding the heterogeneity arising from the complex architecture of sedimentary sequences in alluvial fans is challenging. This study develops a statistical inverse framework in a multi-zone transition probability approach for characterizing the heterogeneity in alluvial fans. An analytical solution of the transition probability matrix is used to define the statistical relationships among different hydrofacies and their mean lengths, integral scales, and volumetric proportions. A statistical inversion is conducted to identify the multi-zone transition probability models and estimate the optimal statistical parameters using the modified Gauss–Newton–Levenberg–Marquardt method. The Jacobian matrix is computed by the sensitivity equation method, which results in an accurate inverse solution with quantification of parameter uncertainty. We use the Chaobai River alluvial fan in the Beijing Plain, China, as an example for elucidating the methodology of alluvial fan characterization. The alluvial fan is divided into three sediment zones. In each zone, the explicit mathematical formulations of the transition probability models are constructed with optimized different integral scales and volumetric proportions. The hydrofacies distributions in the three zones are simulated sequentially by the multi-zone transition probability-based indicator simulations. Finally, the result of this study provides the heterogeneous structure of the alluvial fan for further study of flow and transport simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Multiparameter Methods are the New Field in Statistics

During the last decade statistical problems of a new kind have appeared more and more often, in which the dimensionality of observations is large, and the sample size is small or comparable with the dimensionality of observations. These problems may be characterized as multiparameter problems. The theory of multiparameter analysis, proposed by A.N.Kolmogorov and supported by the studies on the spectral theory of random matrices of V. A.Marchenko, L. A.Pastur, V. L.Girko, and others, revealed some specific phenomena appearing in statistics with a large number of weakly dependent variables. Particularly, there are stable relations between the principal parts of parameter set functions and the set of observed variables, which may be used to improve the statistical analysis. The use of these phenomena allows one to develop a thorough, always stable, and approximately optimal (irrespective of samples) versions of mostly used statistical procedures. The multiparameter versions of the discriminant analysis are of special interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Extension of the Concept of Invariance and Statistically Weakly Invariant Sets of Controllable Systems

We continue the study of statistically invariant and statistically weakly invariant sets with respect to controllable systems and differential inclusions launched by Prof. E. L. Tonkov. We examine properties of such statistical characteristics as the lower freq{sub *}(𝜑) and upper freq{sup *}(𝜑) relative frequencies of hitting a solution 𝜑(t) of a differential inclusion in a prescribed set. We obtain estimates and conditions of coincidence of these characteristics for functions whose difference tends to zero at infinity. We also present conditions of statistically weak invariance of a given set of a relatively controllable system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

An ensemble Kalman filter for statistical estimation of physics constrained nonlinear regression models

A central issue in contemporary science is the development of nonlinear data driven statistical–dynamical models for time series of noisy partial observations from nature or a complex model. It has been established recently that ad-hoc quadratic multi-level regression models can have finite-time blow-up of statistical solutions and/or pathological behavior of their invariant measure. Recently, a new class of physics constrained nonlinear regression models were developed to ameliorate this pathological behavior. Here a new finite ensemble Kalman filtering algorithm is developed for estimating the state, the linear and nonlinear model coefficients, the model and the observation noise covariances from available partial noisy observations of the state. Several stringent tests and applications of the method are developed here. In the most complex application, the perfect model has 57 degrees of freedom involving a zonal (east–west) jet, two topographic Rossby waves, and 54 nonlinearly interacting Rossby waves; the perfect model has significant non-Gaussian statistics in the zonal jet with blocked and unblocked regimes and a non-Gaussian skewed distribution due to interaction with the other 56 modes. We only observe the zonal jet contaminated by noise and apply the ensemble filter algorithm for estimation. Numerically, we find that a three dimensional nonlinear stochastic model with one level of memory mimics the statistical effect of the other 56 modes on the zonal jet in an accurate fashion, including the skew non-Gaussian distribution and autocorrelation decay. On the other hand, a similar stochastic model with zero memory levels fails to capture the crucial non-Gaussian behavior of the zonal jet from the perfect 57-mode model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A model and variance reduction method for computing statistical outputs of stochastic elliptic partial differential equations

We present a model and variance reduction method for the fast and reliable computation of statistical outputs of stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the hybridizable discontinuous Galerkin (HDG) discretization of elliptic partial differential equations (PDEs), which allows us to obtain high-order accurate solutions of the governing PDE; (2) the reduced basis method for a new HDG discretization of the underlying PDE to enable real-time solution of the parameterized PDE in the presence of stochastic parameters; and (3) a multilevel variance reduction method that exploits the statistical correlation among the different reduced basis approximations and the high-fidelity HDG discretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the reduced basis approximations. Furthermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the reduced basis approximations and the sizes of Monte Carlo samples to achieve a given error tolerance. We provide numerical examples to demonstrate the performance of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Statistically informed upscaling of damage evolution in brittle materials

The presence and growth of micro-cracks degrade the strength of brittle solids, greatly impacting the overall material response. Hence, the evolution of these micro-cracks must be accounted for in models describing the relationship between stress and strain so that accurate predictions of material failure can be made. The evolution of individual cracks and crack networks can be simulated with high-fidelity microscale models utilizing highly resolved meshes that can be computationally expensive to a point in which it limits the simulation scale. Hence, for many engineering applications that require simulations of large components, continuum-scale models, which cannot explicitly resolve individual cracks and thus lose important physical information, are required. In this work, we bridge these two scales by developing and implementing a continuum-scale effective moduli constitutive model that is informed by crack statistics generated from a high-fidelity model resolved using a finite-discrete element method (FDEM) implementation. Using statistical information describing the evolution of crack lengths and orientations, this model can capture the effects of brittle damage evolution without the need to resolve individual cracks. We have successfully captured the stress-strain behavior of the high-fidelity simulations using the statistics-based constitutive model shown through direct comparison of stress-strain curves. The curves match within error bars present in the strain-softening portions of the stress-strain curve of the high-fidelity results due to the statistical variation of the initial pre-existing crack network. In conclusion, the stress-strain curves are also compared to experimental results for similar loading conditions and show good qualitative agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Statistical Inference Over Persistent Homology Predicts Fluid Flow in Porous Media

Here, we statistically infer fluid flow and transport properties of porous materials based on their geometry and connectivity, without the need for detailed We summarize structure by persistent homology and then determines the similarity of structures using image analysis and statistics. Longer term, this may enable quick and automated categorization of rocks into known archetypes. We first compute persistent homology of binarized 3D images of material subvolume samples. The persistence parameter is the signed Euclidean distance from inferred material interfaces, which captures the distribution of sizes of pores and grains. Each persistence diagram is converted into an image vector. Here, we infer structural similarity by calculating image similarity. For each image vector, we compute principal components to extract features. We fit statistical models to features estimates material permeability, tortuosity, and anisotropy. We develop a Structural SIMilarity index to determine statistical representative elementary volumes.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Effects of forcing time scale on the simulated turbulent flows and turbulent collision statistics of inertial particles

In this paper, we study systematically the effects of forcing time scale in the large-scale stochastic forcing scheme of Eswaran and Pope [“An examination of forcing in direct numerical simulations of turbulence,” Comput. Fluids 16, 257 (1988)] on the simulated flow structures and statistics of forced turbulence. Using direct numerical simulations, we find that the forcing time scale affects the flow dissipation rate and flow Reynolds number. Other flow statistics can be predicted using the altered flow dissipation rate and flow Reynolds number, except when the forcing time scale is made unrealistically large to yield a Taylor microscale flow Reynolds number of 30 and less. We then study the effects of forcing time scale on the kinematic collision statistics of inertial particles. We show that the radial distribution function and the radial relative velocity may depend on the forcing time scale when it becomes comparable to the eddy turnover time. This dependence, however, can be largely explained in terms of altered flow Reynolds number and the changing range of flow length scales present in the turbulent flow. We argue that removing this dependence is important when studying the Reynolds number dependence of the turbulent collision statistics. The results are also compared to those based on a deterministic forcing scheme to better understand the role of large-scale forcing, relative to that of the small-scale turbulence, on turbulent collision of inertial particles. To further elucidate the correlation between the altered flow structures and dynamics of inertial particles, a conditional analysis has been performed, showing that the regions of higher collision rate of inertial particles are well correlated with the regions of lower vorticity. Regions of higher concentration of pairs at contact are found to be highly correlated with the region of high energy dissipation rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Monte Carlo simulations of time-of-flight PET with double-ended readout: calibration, coincidence resolving times and statistical lower bounds

Here, this paper demonstrates through Monte Carlo simulations that a practical positron emission tomograph with (1) deep scintillators for efficient detection, (2) double-ended readout for depth-of-interaction information, (3) fixed-level analog triggering, and (4) accurate calibration and timing data corrections can achieve a coincidence resolving time (CRT) that is not far above the statistical lower bound. One Monte Carlo algorithm simulates a calibration procedure that uses data from a positron point source. Annihilation events with an interaction near the entrance surface of one scintillator are selected, and data from the two photodetectors on the other scintillator provide depth-dependent timing corrections. Another Monte Carlo algorithm simulates normal operation using these corrections and determines the CRT. A third Monte Carlo algorithm determines the CRT statistical lower bound by generating a series of random interaction depths, and for each interaction a set of random photoelectron times for each of the two photodetectors. The most likely interaction times are determined by shifting the depth-dependent probability density function to maximize the joint likelihood for all the photoelectron times in each set. Example calculations are tabulated for different numbers of photoelectrons and photodetector time jitters for three 3 × 3 × 30 mm 3 scintillators: Lu 2 SiO 5 :Ce,Ca (LSO), LaBr 3 :Ce, and a hypothetical ultra-fast scintillator. To isolate the factors that depend on the scintillator length and the ability to estimate the DOI, CRT values are tabulated for perfect scintillator-photodetectors. For LSO with 4000 photoelectrons and single photoelectron time jitter of the photodetector J = 0.2 ns (FWHM), the CRT value using the statistically weighted average of corrected trigger times is 0.098 ns FWHM and the statistical lower bound is 0.091 ns FWHM. For LaBr 3 :Ce with 8000 photoelectrons and J = 0.2 ns FWHM, the CRT values are 0.070 and 0.063 ns FWHM, respectively. For the ultra-fast scintillator with 1 ns decay time, 4000 photoelectrons, and J = 0.2 ns FWHM, the CRT values are 0.021 and 0.017 ns FWHM, respectively. Lastly, the examples also show that calibration and correction for depth-dependent variations in pulse height and in annihilation and optical photon transit times are necessary to achieve these CRT values.

60 APPLIED LIFE SCIENCES↗