Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistically designed experiments”

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 793 records · Page 44

Dynamic performance of an aero-assist spacecraft - AFE

Dynamic performance of the Aero-assist Flight Experiment (AFE) spacecraft was investigated using a high-fidelity 6-DOF simulation model. Baseline guidance logic, control logic, and a strapdown navigation system to be used on the AFE spacecraft are also modeled in the 6-DOF simulation. During the AFE mission, uncertainties in the environment and the spacecraft are described by an error space which includes both correlated and uncorrelated error sources. The principal error sources modeled in this study include navigation errors, initial state vector errors, atmospheric variations, aerodynamic uncertainties, center-of-gravity off-sets, and weight uncertainties. The impact of the perturbations on the spacecraft performance is investigated using Monte Carlo repetitive statistical techniques. During the Solid Rocket Motor (SRM) deorbit phase, a target flight path angle of -4.76 deg at entry interface (EI) offers very high probability of avoiding SRM casing skip-out from the atmosphere. Generally speaking, the baseline designs of the guidance, navigation, and control systems satisfy most of the science and mission requirements.

Chang, Ho-Pen↗

Inland and Near Shore Water Profiles Derived from the High Altitude Multiple Altimeter Beam Experimental Lidar (MABEL)

The Advanced Topographic Laser Altimeter System (ATLAS) on the Ice, Cloud, and Land Elevation Satellite (ICESat-2) mission is a six beam, low energy, high repetition rate, 532 nm laser transmitter with photon counting detectors. Although designed primarily for detecting height changes in icecaps, sea ice and vegetation, the polar-orbital satellite will observe global surface water during its designed three year life span, including inland water bodies, coasts, and open oceans. In preparation for the mission, an ICESat-2 prototype or the Multiple Altimeter Beam Experimental Lidar (MABEL), was built and flown on high altitude aircraft experiments over a range of inland and near-shore targets. The purpose was to test the ATLAS concept and to provide a database for developing an algorithm that detects along track surface water height and light penetration under a range of atmospheric and water conditions. The current analysis examines the datasets of three MABEL transects observed from 20 km above ground of coastal and inland waters conducted in 2012 and 2013. Transects ranged from about 2 to 12 km in length and included the middle Chesapeake Bay, the near shore Atlantic coast at Virginia Beach, and Lake Mead. Results indicate MABEL's high capability for retrieving surface water height statistics with a mean height precision of approximately 5-7 cm per 100m segment length. Profiles of attenuated subsurface backscatter, characterized using a Signal to Background Ratio written in Log10 base, or LSBR0, were observed over a range of 1.3 to 9.3 meters depending on water clarity and atmospheric background. Results indicate that observable penetration depth, although primarily dependent on water properties, was greatest when solar background rate was low. Near shore bottom reflectance was detected only at the Lake Mead site down to maximum of 10 m under a clear night sky and low turbidity of approximately 1.6 Nephelometric Turbidity Units (NTU). The overall results suggest that the feasibility of retrieving operational surface water height statistics from space-based photon counting systems such as ATLAS is very high for resolutions down to about 100m, even in partly cloudy conditions. The capability to observe subsurface backscatter profiles is achievable but requires much longer transects of several hundreds of meters.

Topographic↗

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Radio occultation studies of the Venus atmosphere with the Magellan spacecraft. 1: Experimental description and performance

While primarily designed for radar studies of the Venus surface, the high effective isotropic radiated power (EIRP) from the Magellan spacecraft makes it an ideal transmitter for use in radio occultation measurements of the refractivity and absorptivity of the Venus atmosphere. Such experiments have been conducted involving transmissions at 2.3 GHz and 8.4 GHz (13 cm and 3.6 cm, respectively), during spacecraft ingress. Since the stability of the spacecraft transmitter is critical for accurately determining the Doppler shift and amplitude attenuation created as the ray penetrates the atmosphere, the spacecraft transmitter was locked to a 2.1 GHz uplink from a 70-meter DSN station which also received the signals. Because of the high directivity of the spacecraft antenna, and the significant ray bending in the deep Venus atmosphere, a spacecraft tracking maneuver was designed to keep the spacecraft antenna pointed in the direction of the refracted ray path back to Earth. This tracking maneuver, plus the high EIRP of the Magellan transmitter has yielded 3.6 cm refractivity and absorptivity profiles down to the 35 km altitude and 13 cm profiles down to the altitude of critical refraction (approximately 33 km). The statistical uncertainties in the derived profiles are significantly lower than those previously obtained, resulting in extremely accurate profiles of H2SO4 (g) abundance as discussed in an accompanying paper.

Steffes, Paul G.↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessing Interpersonal and Communication Skills in Radiation Oncology Residents: A Pilot Standardized Patient Program

Purpose: There is a lack of data for the structured development and evaluation of communication skills in radiation oncology residency training programs. Effective communication skills are increasingly emphasized by the Accreditation Council for Graduate Medical Education and are critical for a successful clinical practice. We present the design of a novel, pilot standardized patient (SP) program and the evaluation of communication skills among radiation oncology residents. Methods and Materials: Two case scenarios were developed to challenge residents in the delivery of “bad news” to patients: one scenario regarding treatment failure and the other regarding change in treatment plan. Eleven radiation oncology residents paired with 6 faculty participated in this pilot program. Each encounter was scored by the SPs, observing faculty, and residents themselves based on the Kalamazoo guidelines. Results: Overall resident performance ratings were “good” to “excellent,” with faculty assigning statistically significant higher scores and residents assigning lower scores. We found inconsistent inter rater agreement among faculty, residents, and SPs. SP feedback was also valuable in identifying areas of improvement, including more collaborative decision making and less use of medical jargon. Conclusions: The program was well received by residents and faculty and regarded as a valuable educational experience that could be used as an annual feedback tool. Poor inter rater agreement suggests a need for residents and faculty physicians to better calibrate their evaluations to true patient perceptions. High scores from faculty members substantiate the concern that resident evaluations are generally positive and nondiscriminating. Faculty should be encouraged to provide honest and critical feedback to hone residents' interpersonal skills.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Intrinsic Kinetics of Polyethylene Terephthalate Pyrolysis via Micropyrolysis and Multivariate Chromatographic Analysis

This study provides an in-depth investigation of the primary decomposition of polyethylene terephthalate (PET) via pyrolysis, employing an experimental-analytic workflow that integrates design of experiments (DoE), micropyrolysis coupled with comprehensive two-dimensional gas chromatography (GC×GC), and multivariate data analysis to verify intrinsic kinetic conditions and elucidate evolving product distributions for mapping key reaction pathways. Peaks that could not be identified using commercial spectral libraries were assigned using Mass Frontier simulations, enabling the identification of divinyl terephthalate, ethyl vinyl terephthalate, and 2-(benzoyloxy)ethyl vinyl terephthalate. A polar×polar (non-orthogonal) column set tailored for the detection of carboxylic acids enhanced the quantification of benzoic acid, 4-vinylbenzoic acid, 4-ethylbenzoic acid, and methylbenzoic acid by up to 6-fold relative to an orthogonal column combination (non-polar×mid-polar). Moreover, pyrolysis variables were systematically evaluated using a Box- Behnken design (BBD), encompassing pyrolysis temperature (500−600 °C), sample weight (50−150 μg), and carrier gas flow rate (100−300 mL min −1 ). Among these, pyrolysis temperature was the only statistically significant factor influencing product yields, ranging from 58.78 to 84.26 wt %. In contrast, neither the sample weight nor the carrier gas flow rate had a significant effect on product yields within the evaluated experimental space. At 600 °C, the major pyrolysis products were benzoic acid (up to 20.20 ± 1.46 wt %) and CO 2 (up to 21.28 ± 1.46 wt %), which can be produced through decarboxylation reactions. These findings underscore the critical importance of selecting appropriate analytical columns for the accurate quantification of heteroatomcontaining products such as carboxylic acids, which may otherwise be underestimated or undetected due to their reactivity with the stationary phase of non-polar and mid-polar columns, as well as other GC components. They also highlight the importance of selecting pyrolysis conditions for investigating the primary decomposition of PET under an isothermal kinetically limited regime.

aromatic compounds↗

Monitoring of International Space Station Telemetry Using Shewhart Control Charts

Shewhart control charts have been established as an expedient method for analyzing dynamic, trending data in order to identify anomalous subsystem performance as soon as such performance would exceed a statistically established baseline. Additionally, this leading indicator tool integrates a selection methodology that reduces false positive indications, optimizes true leading indicator events, minimizes computer processor unit duty cycles, and addresses human factor concerns (i.e., the potential for flight-controller data overload). This innovation leverages statistical process control, and provides a relatively simple way to allow flight controllers to focus their attention on subtle system changes that could lead to dramatic off-nominal system performance. Finally, this capability improves response time to potential hardware damage and/or crew injury, thereby improving space flight safety. Shewhart control charts require normalized data. However, the telemetry from the ISS Early External Thermal Control System (EETCS) was not normally distributed. A method for normalizing the data was implemented, as was a means of selecting data windows, the number of standard deviations (Sigma Level), the number of consecutive points out of limits (Sequence), and direction (increasing or decreasing trend data). By varying these options, and treating them like dial settings, the number of nuisance alerts and leading indicators were optimized. The goal was to capture all leading indicators while minimizing the number of nuisances. Lean Six Sigma (L6S) design of experiment methodologies were employed. To optimize the results, Perl programming language was used to automate the massive amounts of telemetry data, control chart plots, and the data analysis.

Fitch, Jeffery T.↗

Inland and Near-Shore Water Profiles Derived from the High-Altitude Multiple Altimeter Beam Experimental Lidar (MABEL)

The Advanced Topographic Laser Altimeter System (ATLAS) on the Ice, Cloud, and Land Elevation Satellite (ICESat-2) mission is a six beam, low energy, high repetition rate, 532-nanometer laser transmitter with photon counting detectors. Although designed primarily for detecting height changes in ice caps, sea ice, and vegetation, the polar-orbiting satellite will observe global surface water during its designed three-year life span, including inland waterbodies, coasts, and open oceans. In preparation for the mission, an ICESat-2 prototype, the Multiple Altimeter Beam Experimental Lidar (MABEL), was built and flown on high-altitude aircraft experiments over a range of inland and near-shore targets. The purpose was to test the ATLAS concept and to provide a database for developing an algorithm that detects along track surface water height and light penetration under a range of atmospheric and water conditions. The current analysis examines the data sets of three MABEL transects observed from 20 kilometers above ground of coastal and inland waters conducted in 2012 and 2013. Transects ranged from about 2 to 12 kilometers in length and included the middle Chesapeake Bay, the near-shore Atlantic coast at Virginia Beach, and Lake Mead. Results indicate MABEL's high capability for retrieving surface water height statistics with a mean height precision ofapproximately 5-7 centimeters per 100-meter segment length. Profiles of attenuated subsurface backscatter, characterized using a Signal to Background Ratio written in Log10 base, or LSBR (sub 0), were observed over a range of 1.3 to 9.3 meters, depending on water clarity and atmospheric background. Results indicate that observable penetration depth, although primarily dependent on water properties, was greatest when the solar background rate was low. Near-shore bottom reflectance was detected only at the Lake Mead site down to a maximum of 10 meters under a clear night sky and low turbidity of approximately 1.6 Nephelometric Turbidity Units (NTU). The overall results suggest that the feasibility of retrieving operational surface water height statistics from space-based photon counting systems such as ATLAS is very high for resolutions down to about 100 meters, even in partly cloudy conditions. The capability to observe subsurface backscatterprofiles is achievable but requires much longer transects of several hundreds of meters.

inland↗

The DESI-Lensing Mock Challenge: large-scale cosmological analysis of 3x2-pt statistics

The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of algorithms for understanding the temporal and spatial variability of the earth's radiation balance

A brief description is given of how temporal and spatial variability in the earth's radiative behavior influences the goals of satellite radiation monitoring systems and how some previous systems have addressed the existing problems. Then, results of some simulations of radiation budget monitoring missions are presented. These studies led to the design of the Earth Radiation Budget Experiment (ERBE). A description is given of the temporal and spatial averaging algorithms developed for the ERBE data analysis. These algorithms are intended primarily to produce monthly averages of the net radiant exitance on regional, zonal, and global scales and to provide insight into the regional diurnal variability of radiative parameters such as albedo and long-wave radiant exitance. The algorithms are applied to scanner and nonscanner data for up to three satellites. Modeling of dialy shortwave albedo and radiant exitance with satellite samling that is insufficient to fully account for changing meteorology is discussed in detail. Studies performed during the ERBE mission and software design are reviewed. These studies provide quantitative estimates of the effects of temporally sparse and biased sampling on inferred diurnal and regional radiative parameters. Other topics covered include long-wave diurnal modeling, extraction of a regional monthly net clear-sky radiation budget, the statistical significance of observed diurnal variability, quality control of the analysis, and proposals for validating the results of ERBE time and space averaging.

Brooks, D. R.↗

Developing a Model-based Capability to Analyze Requirements for the Climate Observing System

Models are foundational for estimating states of the earth's climate system, both as tools to extrapolate information in time and space, and as observation 'operators' used to relate what is analyzed and predicted to what is observed. Expanding the simulation approach further, observing system simulation experiments (OSSEs) are designed to mimic the complete process of analyzing the climate state by replacing real observations with entirely simulated ones determined from a model-based depiction of nature. OSSEs provide a framework to 'fly' simulated satellite instruments through a synthetic atmosphere and investigate the trade-spaces of measurements for various satellite configurations and sampling strategies, and assess their measurement impact on modeling and forecasting capabilities. Such a tool is a crucial but as yet unfulfilled need for future mission selection and design. The components of a state-of-the-art OSSE system are being assembled at the Global Modeling and Assimilation Office (GMAO, Code 610.1) at NASA/GSFC, leveraging on the GMAO's existing modeling and data assimilation infrastructure for numerical weather prediction (NWP). The OSSE framework is based on the GMAO's Goddard Earth Observing System atmospheric general circulation model, version 5 (GEOS-5) and the Gridpoint Statistical Interpolation (GSI) observational analysis scheme, combined with the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model developed by the Atmospheric Chemistry and Dynamics Branch (Code 613.3). This system is an evolving, key component of Goddard's planned development of an Integrated Earth System Analysis (IESA) capability, which will bring together into a single, fully interactive system Goddard's modeling and assimilation efforts in atmosphere, ocean and chemistry and aerosols to provide a comprehensive analysis and prediction system for weather and climate In addition to providing a state-of-the-art capability for assimilating current observation types, GEOS-5, and the future IESA, provide the capability to identify the need for, and assess the potential impact of, future observing systems under consideration for improving weather and climate prediction.

Gelaro, Ronald↗

Aerobic Capacity Following Long Duration International Spaces Station (ISS) Missions: Preliminary Results

Maximum oxygen uptake (VO2max) is reduced immediately following space flights lasting <15 d, but has not been measured following long-duration missions. The purpose of this study is to measure VO2max and maximum work rate (WRmax) data from astronauts following ISS flights (91 to 188 d). Methods: Five astronauts [3 M, 2 F: 47+/-6 yr, 174+/-6 cm, 71.9+/-10.9 kg (mean +/- SD)] have participated in the study. Subjects performed upright cycle exercise tests to symptom-limited maximum. An initial test was done approx.270 d before flight to establish work rates for subsequent tests. Subsequent tests, conducted approx.45 d before flight and repeated on the first or second day (R+1/2) and at approx.10 d (R+10) following landing, consisted of 3 5 min stages designed to elicit 25%, 50%, and 75% of preflight VO2max, followed by 25 W(dot)/min increases. VO2, WR, and heart rate (HR) were measured using the ISS Portable Pulmonary Function System [Damec, Odense, DK]. Descriptive statistics are reported. Results: On R+1/2 mean VO2max decreased compared to preflight (Pre: 2.98+/-0.99, R+1/2: 2.63+/-0.56 L(dot)/min); 4 of 5 subjects demonstrated a loss of > 6%. WRmax also decreased on R+1/2 compared to preflight (Pre: 245+/-69, R+1/2: 210+/-45 W). On R+10, VO2max was 2.86+/-0.62 L(dot)/min, with 2 subjects still demonstrating a loss of > 6% from preflight. WRmax on R+10 was 240+/-49 W. HRmax did not change from pre to post-flight. Conclusions: These preliminary results, from the first 5 of 12 planned subjects of an ongoing ISS study, suggest that the majority of astronauts will experience a decrease in VO2max after long-duration space-flight. Interestingly, the two astronauts with the highest preflight VO2max had the greatest loss on R+1/2, and the astronaut with the lowest preflight VO2max increased by 13%. Thus, maintenance of VO2max may be more difficult in astronauts who have a high aerobic capacity, perhaps requiring more intense in-flight exercise countermeasure prescriptions.

Moore, Alan D.↗

Mitigating the impact of the DESI fiber assignment on galaxy clustering

In this paper, we present a simple strategy to mitigate the impact of an incomplete spectroscopic redshift galaxy sample as a result of fiber assignment and survey tiling. The method has been designed for the Dark Energy Spectroscopic Instrument (DESI) galaxy survey but may have applications beyond this. We propose a modification to the usual correlation function that nulls the almost purely angular modes affected by survey incompleteness due to fiber assignment. Predictions of this modified statistic can be calculated given a model of the two point correlation function. The new statistic can be computed with a slight modification to the data catalogues input to the standard correlation function code and does not incur any additional computational time. Lastly, we show that the spherically averaged baryon acoustic oscillation signal is not biased by the new statistic.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of Long-Term Quality Control Data for a 137 Cs Dosimetry Calibration Source

Strict quality assurance programs are required for many radiological applications, but these seldom exist for verifying dosimetry calibration sources. After initial characterization of a dosimetry calibration facility, quality control procedures are recommended to ensure the early detection of any changes or malfunctions. These also result in refined knowledge about average dose rate and experimental variations in dose delivery. This paper describes the implementation of a phase I quality control protocol for a 137 Cs dosimetry calibration source and includes an analysis of the resulting data collected over a 24-mo period. During this time, substantial data was collected to establish trial control limits. Air kerma rate measurements were obtained using an ion chamber and were adjusted for decay, corrected for ambient temperature, pressure and humidity, and then analyzed using quality control charts. Three variations of rational subgrouping methods were used in order to find assignable causes of error, and Nelson's Rules were followed to detect any non-random statistical variations. Measurements were subgrouped according to same-day measurements in order to detect positional errors as well as atmospheric correction errors. Additionally, measurements were subgrouped according to analogous experimental setups in order to detect failure in equipment or incorrect settings. Both were analyzed using the X-bar and R chart method. Similarly, individuals and moving ranges charts were used to carefully examine each position in order to observe any situational errors that may occur which include timing, positional, or interference errors. Each method was successful in identifying unique out-of-control data points that occurred during the phase I application of forming control limits. Furthermore, over the 24-mo period, enough data points were deemed in-control to establish reliable trial limits. Future experiments will include the phase II application of gaining more reliable measurements in order to fine-tune the limits, as well as performing a designed experiment, where variables are purposefully changed in order to test the variation of the data.

61 RADIATION PROTECTION AND DOSIMETRY↗

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metalearners for estimating heterogeneous treatment effects using machine learning

There is growing interest in estimating and analyzing heterogeneous treatment effects in experimental and observational studies. We describe a number of metaalgorithms that can take advantage of any supervised learning or regression method in machine learning and statistics to estimate the conditional average treatment effect (CATE) function. Metaalgorithms build on base algorithms—such as random forests (RFs), Bayesian additive regression trees (BARTs), or neural networks—to estimate the CATE, a function that the base algorithms are not designed to estimate directly. We introduce a metaalgorithm, the X-learner, that is provably efficient when the number of units in one treatment group is much larger than in the other and can exploit structural properties of the CATE function. For example, if the CATE function is linear and the response functions in treatment and control are Lipschitz-continuous, the X-learner can still achieve the parametric rate under regularity conditions. We then introduce versions of the X-learner that use RF and BART as base learners. In extensive simulation studies, the X-learner performs favorably, although none of the metalearners is uniformly the best. In two persuasion field experiments from political science, we demonstrate how our X-learner can be used to target treatment regimes and to shed light on underlying mechanisms. A software package is provided that implements our methods.

42 ENGINEERING↗