Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical techniques”

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 703 records · Page 39

Shear dynamo, turbulence, and the magnetorotational instability

The formation, evolution, and detailed structure of accretion disks remain poorly understood, with wide implications across a variety of astrophysical disciplines. While the most pressing question – what causes the high angular momentum fluxes that are necessary to explain observations? – is nicely answered by the idea that the disk is turbulent, a more complete grasp of the fundamental processes is necessary to capture the wide variety of behaviors observed in the night sky. This thesis studies the turbulence in ionized accretion disks from a theoretical standpoint, in particular focusing on the generation of magnetic fields in these processes, known as dynamo. Such fields are expected to be enormously important, both by enabling the magnetorotational instability (which evolves into virulent turbulence), and through large-scale structure formation, which may transport angular momentum in different ways and be fundamental for the formation of jets. The central result of this thesis is the suggestion of a new large-scale dynamo mechanism in shear flows – the “magnetic shear-current effect” – which relies on a positive feedback from smallscale magnetic fields. As well as being a very promising candidate for driving field generation in the central regions of accretion disks, this effect is interesting because small-scale magnetic fields have historically been considered to have a negative effect on the large-scale dynamo, damping growth and leading to dire predictions for final saturation amplitudes. Given that small-scale fields are ubiquitous in plasma turbulence above moderate Reynolds numbers, the finding that they could instead have a positive effect in some situations is interesting from a theoretical and practical standpoint. The effect is studied using direct numerical simulation, analytic techniques, and novel statistical simulation methods. In addition to the dynamo, much attention is given to the linear physics of disks and its relevance to turbulence. This is studied using nonmodal stability theory, which both provides a highly intuitive connection between global domains and the commonly studied shearing box, and suggests that transient linear growth can often be more important than spectral instability. These realizations motivate the use of the quasi-linear models that are applied extensively throughout the turbulence and dynamo studies later in the thesis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupling Visualization, Simulation, and Deep Learning for Ensemble Steering of Complex Energy Models: Preprint

We have developed a framework for the exploration, design, and planning of energy systems that combines interactive visualization with machine-learning based approximations of simulations through a general purpose dataflow API. Our system provides a visual inter- face allowing users to explore an ensemble of energy simulations representing a subset of the complex input parameter space, and spawn new simulations to 'fill in' input regions corresponding to new enegery system scenarios. Unfortunately, many energy simula- tions are far too slow to provide interactive responses. To support interactive feedback, we are developing reduced-form models via machine learning techniques, which provide statistically sound esti- mates of the full simulations at a fraction of the computational cost and which are used as proxies for the full-form models. Fast com- putation and an agile dataflow enhance the engagement with energy simulations, and allow researchers to better allocate computational resources to capture informative relationships within the system and provide a low-cost method for validating and quality-checking large-scale modeling efforts.

97 MATHEMATICS AND COMPUTING↗

Modeling the joint statistics of satellite path XPD and attenuation

A new and simple technique for estimating from attenuation statistics alone the joint distribution of attenuation and cross-polarization impairments on a satellite-to-earth radio path is presented. It compares measurements with predictions for the CTS spacecraft 11.7 GHz right-hand circularly polarized downlink during the first six months of 1979.

Gaines, J. M.↗

Computational techniques for the analysis of small signals in high-statistics neutrino oscillation experiments

The current and upcoming generation of Very Large Volume Neutrino Telescopes – collecting unprecedented quantities of neutrino events – can be used to explore subtle effects in oscillation physics, such as (but not restricted to) the neutrino mass ordering. The sensitivity of an experiment to these effects can be estimated from Monte Carlo simulations. With the high number of events that will be collected, there is a trade-off between the computational expense of running such simulations and the inherent statistical uncertainty in the determined values. In such a scenario, it becomes impractical to produce and use adequately-sized sets of simulated events with traditional methods, such as Monte Carlo weighting. In this work we present a staged approach to the generation of expected distributions of observables in order to overcome these challenges. By combining multiple integration and smoothing techniques which address limited statistics from simulation it arrives at reliable analysis results using modest computational resources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Information and Statistics in Nuclear Experiment and Theory (ISNET)

As with all empirical sciences, nuclear physics operates in the virtuous cycle of the scientific method: observations inspire theoretical models; models lead to new predictions; predictions are tested in experiments; experiments lead to new observations; and so on. Evaluating what we are inferring, and how certain we are of it, is key to this process. These requirements, and a general interest in applying novel statistical, mathematical, and computational techniques, led to the formation of a dedicated research community entitled “Information and Statistics in Nuclear Experiment and Theory (ISNET)” (https://isnet-series.github.io/), which now includes more than 300 members. While the community’s interests lean toward nuclear theory, the unifying theme for this group is the inference of knowledge from data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cellular-based preemption system

A cellular-based preemption system that uses existing cellular infrastructure to transmit preemption related data to allow safe passage of emergency vehicles through one or more intersections. A cellular unit in an emergency vehicle is used to generate position reports that are transmitted to the one or more intersections during an emergency response. Based on this position data, the one or more intersections calculate an estimated time of arrival (ETA) of the emergency vehicle, and transmit preemption commands to traffic signals at the intersections based on the calculated ETA. Additional techniques may be used for refining the position reports, ETA calculations, and the like. Such techniques include, without limitation, statistical preemption, map-matching, dead-reckoning, augmented navigation, and/or preemption optimization techniques, all of which are described in further detail in the above-referenced patent applications.

Bachelder, Aaron D.↗

SU-E-T-497: Semi-Automated in Vivo Radiochromic Film Dosimetry Using a Novel Image Processing Algorithm

Purpose: To validate an automated image processing algorithm designed to detect the center of radiochromic film used for in vivo film dosimetry against the current gold standard of manual selection. Methods: An image processing algorithm was developed to automatically select the region of interest (ROI) in *.tiff images that contain multiple pieces of radiochromic film (0.5x1.3cm{sup 2}). After a user has linked a calibration file to the processing algorithm and selected a *.tiff file for processing, an ROI is automatically detected for all films by a combination of thresholding and erosion, which removes edges and any additional markings for orientation. Calibration is applied to the mean pixel values from the ROIs and a *.tiff image is output displaying the original image with an overlay of the ROIs and the measured doses. Validation of the algorithm was determined by comparing in vivo dose determined using the current gold standard (manually drawn ROIs) versus automated ROIs for n=420 scanned films. Bland-Altman analysis, paired t-test, and linear regression were performed to demonstrate agreement between the processes. Results: The measured doses ranged from 0.2-886.6cGy. Bland-Altman analysis of the two techniques (automatic minus manual) revealed a bias of -0.28cGy and a 95% confidence interval of (5.5cGy,-6.1cGy). These values demonstrate excellent agreement between the two techniques. Paired t-test results showed no statistical differences between the two techniques, p=0.98. Linear regression with a forced zero intercept demonstrated that Automatic=0.997*Manual, with a Pearson correlation coefficient of 0.999. The minimal differences between the two techniques may be explained by the fact that the hand drawn ROIs were not identical to the automatically selected ones. The average processing time was 6.7seconds in Matlab on an IntelCore2Duo processor. Conclusion: An automated image processing algorithm has been developed and validated, which will help minimize user interaction and processing time of radiochromic film used for in vivo dosimetry.

60 APPLIED LIFE SCIENCES↗

Calibration and Commissioning of the LSST

Understanding the nature of dark energy and dark matter remains one of the fundamental questions in physics today; impacting our understanding of particle physics, cosmology, and possibly theories of gravity. Given the scale and complexity of the next generation of cosmology experiments (e.g., the Rubin Observatory, the Euclid satellite mission, and the Roman space telescope) we are entering an era where statistical noise no longer determines the accuracy to which we can measure cosmological parameters. Our ability to control and correct for systematics will ultimately determine the scientific impact of these experiments. This award addressed the challenge of how we determine what limits the accuracy of our cosmological measures, what techniques are appropriate for measuring and calibrating the properties of galaxies to best constrain cosmological models, how to develop statistical techniques that are insensitive to systematic errors, and how to optimize survey strategies in order to minimize systematics while maximizing the speed at which an experiment can achieve its science objectives. In this final technical report for award DE-SC0011635 we describe a set of open-source frameworks that simulate the characteristics and properties of current and planned cosmology surveys and the application of these frameworks to the development of new methodologies for estimating the properties and distances to galaxies that are robust to noisy and incomplete data.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing the photon polarization in B → K*γ with conversion

We re-examine the possibility to measure the photon polarization in B → K*γ decays, via decays in which the photon subsequently undergoes nuclear conversion to a lepton pair. We obtain compact expressions for the full decay-plus-conversion amplitude. With these results we show that interference between the B → (K* → Kπ)γ decay and the γN → ℓ + ℓ – N conversion permits both the ratio and relative weak phase between the left- and right-handed photon amplitudes to be probed by an angular observable, constructed from the final state dilepton, kaon and pion kinematic configuration. Exploiting this technique will be experimentally challenging. However, we present special kinematic cuts that enhance the statistical power of this technique by an O(1) factor. Furthermore, we verify this effect and extract pertinent angular kinematic distributions with dedicated numerical simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adaptive statistical pattern classifiers for remotely sensed data

A technique for the adaptive estimation of nonstationary statistics necessary for Bayesian classification is developed. The basic approach to the adaptive estimation procedure consists of two steps: (1) an optimal stochastic approximation of the parameters of interest and (2) a projection of the parameters in time or position. A divergence criterion is developed to monitor algorithm performance. Comparative results of adaptive and nonadaptive classifier tests are presented for simulated four dimensional spectral scan data.

Gonzalez, R. C.↗

Inflation in random Gaussian landscapes

We develop analytic and numerical techniques for studying the statistics of slow-roll inflation in random Gaussian landscapes. As an illustration of these techniques, we analyze small-field inflation in a one-dimensional landscape. We calculate the probability distributions for the maximal number of e-folds and for the spectral index of density fluctuations n {sub s} and its running α {sub s} . These distributions have a universal form, insensitive to the correlation function of the Gaussian ensemble. We outline possible extensions of our methods to a large number of fields and to models of large-field inflation. These methods do not suffer from potential inconsistencies inherent in the Brownian motion technique, which has been used in most of the earlier treatments.

79 ASTRONOMY AND ASTROPHYSICS↗

Dosimetric difference amongst 3 techniques: TomoTherapy, sliding-window intensity-modulated radiotherapy (IMRT), and RapidArc radiotherapy in the treatment of late-stage nasopharyngeal carcinoma (NPC)

To investigate the dosimetric difference amongst TomoTherapy, sliding-window intensity-modulated radiotherapy (IMRT), and RapidArc radiotherapy in the treatment of late-stage nasopharyngeal carcinoma (NPC). Ten patients with late-stage (Stage III or IV) NPC treated with TomoTherapy or IMRT were selected for the study. Treatment plans with these 3 techniques were devised according to departmental protocol. Dosimetric parameters for organ at risk and treatment targets were compared between TomoTherapy and IMRT, TomoTherapy and RapidArc, and IMRT and RapidArc. Comparison amongst the techniques was done by statistical tests on the dosimetric parameters, total monitor unit (MU), and expected delivery time. All 3 techniques achieved similar target dose coverage. TomoTherapy achieved significantly lower doses in lens and mandible amongst the techniques. It also achieved significantly better dose conformity to the treatment targets. RapidArc achieved significantly lower dose to the eye and normal tissue, lower total MU, and less delivery time. The dosimetric advantages of the 3 techniques were identified in the treatment of late-stage NPC. This may serve as a guideline for selection of the proper technique for different clinical cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Statistical assessment of a unique time series analysis technique

Empirical tests of a particular least squares multiple prediction program were studied, as well as the relationship of this particular program to general least squares multiple prediction theory. Empirical evaluation and test of the program involved: (1) conversion of the program to run on an IBM 360/44 with replication of test results obtained on a Univac 1108; and (2) generation and analysis of Monte Carlo simulated data with the objective of comparison against results theoretically obtainable from spectrum analysis routines.

Benignus, V. A.↗

Analysis techniques for two-dimensional infrared data

In order to evaluate infrared detection and remote sensing systems, it is necessary to know the characteristics of the observational environment. For both scanning and staring sensors, the spatial characteristics of the background may be more of a limitation to the performance of a remote sensor than system noise. This limitation is the so-called spatial clutter limit and may be important for systems design of many earth application and surveillance sensors. The data used in this study is two dimensional radiometric data obtained as part of the continuing NASA remote sensing programs. Typical data sources are the Landsat multi-spectral scanner (1.1 micrometers), the airborne heat capacity mapping radiometer (10.5 - 12.5 micrometers) and various infrared data sets acquired by low altitude aircraft. Techniques used for the statistical analysis of one dimensional infrared data, such as power spectral density (PSD), exceedance statistics, etc. are investigated for two dimensional applicability. Also treated are two dimensional extensions of these techniques (2D PSD, etc.), and special techniques developed for the analysis of 2D data.

Winter, E. M.↗

Parameter estimation techniques based on optimizing goodness-of-fit statistics for structural reliability

New methods are presented that utilize the optimization of goodness-of-fit statistics in order to estimate Weibull parameters from failure data. It is assumed that the underlying population is characterized by a three-parameter Weibull distribution. Goodness-of-fit tests are based on the empirical distribution function (EDF). The EDF is a step function, calculated using failure data, and represents an approximation of the cumulative distribution function for the underlying population. Statistics (such as the Kolmogorov-Smirnov statistic and the Anderson-Darling statistic) measure the discrepancy between the EDF and the cumulative distribution function (CDF). These statistics are minimized with respect to the three Weibull parameters. Due to nonlinearities encountered in the minimization process, Powell's numerical optimization procedure is applied to obtain the optimum value of the EDF. Numerical examples show the applicability of these new estimation methods. The results are compared to the estimates obtained with Cooper's nonlinear regression algorithm.

Starlinger, Alois↗

A technique for generating correlated X-band weather degradation statistics

It is, known that the X-band noise temperature of DSN receivers can go up from 20 K to over 100 K or more, if the air is heavily laden with water vapor, although that is an uncommon occurrence. It is proposed that the DSN furnish flight projects relying on X-band degradation models, one for each DSN Complex. Such models would be in the form of a random process generator, say in an MBASIC program, that would permit the project to generate X-band degradation data with the right autocorrelations for periods of interest to the Projects. The autocorrelation modeling is especially important because bursts of degradation lasting several days can affect data storage and mission sequence design strategy. One approach which works if the degradation statistics obey a half-gaussian law. That is, the random variables are formed by taking the absolute values of another set of random variables, themselves having a Gaussian distribution. The technique then permits the half-gaussian random variables to have given one and two-step correlation coefficients.

Posner, E. C.↗

A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging

Monte Carlo simulations based on weighted simulation particles can solve a variety of population balance problems and allow thus to formulate a solution-framework for many chemical engineering processes. This study presents a novel concept for the calculation of coagulation rates of weighted Monte Carlo particles by introducing a family of transformations to non-weighted Monte Carlo particles. The tuning of the accuracy (named ‘stochastic resolution’ in this paper) of those transformations allows the construction of a constant-number coagulation scheme. Furthermore, a parallel algorithm for the inclusion of newly formed Monte Carlo particles due to nucleation is presented in the scope of a constant-number scheme: the low-weight merging. This technique is found to create significantly less statistical simulation noise than the conventional technique (named ‘random removal’ in this paper). Both concepts are combined into a single GPU-based simulation method which is validated by comparison with the discrete-sectional simulation technique. Two test models describing a constant-rate nucleation coupled to a simultaneous coagulation in 1) the free-molecular regime or 2) the continuum regime are simulated for this purpose.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying High Temperature Corrosion

Alloys designed for high temperature service can form a variety of surface oxides or scales based on the alloy composition, service temperature, gas environment, surface deposits, etc. The most severe environments, including molten salts and liquid metals can result in significant metal loss, void formation and/or pitting. Standard and evolving practices are discussed for measuring these types of degradationincluding measuring oxide scale thickness, internal oxidation and/or metal loss, correlating them to mass change and reporting the results. Manual and automated methodologies to measure oxide thickness and metal loss are reviewed and compared in order to provide a summary of techniques and identify the best practices for quantifying a variety of materials and damage mechanisms and producing statistically meaningful results. Such techniques create datasets useful for improved understanding of corrosion degradation mechanisms, potentially better predictive models and enhanced data analytics.

Su, Yi Feng↗