Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling 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 631 records · Page 35

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Correcting correlation functions for redshift-dependent interloper contamination

ABSTRACT The construction of catalogues of a particular type of galaxy can be complicated by interlopers contaminating the sample. In spectroscopic galaxy surveys this can be due to the misclassification of an emission line; for example in the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift [O ii] emitters may make up a few per cent of the observed Ly α emitter (LAE) sample. The presence of contaminants affects the measured correlation functions and power spectra. Previous attempts to deal with this using the cross-correlation function have assumed sources at a fixed redshift, or not modelled evolution within the adopted redshift bins. However, in spectroscopic surveys like HETDEX, where the contamination fraction is likely to be redshift dependent, the observed clustering of misclassified sources will appear to evolve strongly due to projection effects, even if their true clustering does not. We present a practical method for accounting for the presence of contaminants with redshift-dependent contamination fractions and projected clustering. We show using mock catalogues that our method, unlike existing approaches, yields unbiased clustering measurements from the upcoming HETDEX survey in scenarios with redshift-dependent contamination fractions within the redshift bins used. We show our method returns autocorrelation functions with systematic biases much smaller than the statistical noise for samples with at least as high as 7 per cent contamination. We also present and test a method for fitting for the redshift-dependent interloper fraction using the LAE–[O ii] galaxy cross-correlation function, which gives less biased results than assuming a single interloper fraction for the whole sample.

79 ASTRONOMY AND ASTROPHYSICS↗

THE 60 MONTH ALL-SKY BURST ALERT TELESCOPE SURVEY OF ACTIVE GALACTIC NUCLEUS AND THE ANISOTROPY OF NEARBY AGNs

Surveys above 10 keV represent one of the best resources to provide an unbiased census of the population of active galactic nuclei (AGNs). We present the results of 60 months of observation of the hard X-ray sky with Swift/Burst Alert Telescope (BAT). In this time frame, BAT-detected (in the 15-55 keV band) 720 sources in an all-sky survey of which 428 are associated with AGNs, most of which are nearby. Our sample has negligible incompleteness and statistics a factor of {approx}2 larger over similarly complete sets of AGNs. Our sample contains (at least) 15 bona fide Compton-thick AGNs and 3 likely candidates. Compton-thick AGNs represent {approx}5% of AGN samples detected above 15 keV. We use the BAT data set to refine the determination of the log N-log S of AGNs which is extremely important, now that NuSTAR prepares for launch, toward assessing the AGN contribution to the cosmic X-ray background. We show that the log N-log S of AGNs selected above 10 keV is now established to {approx}10% precision. We derive the luminosity function of Compton-thick AGNs and measure a space density of 7.9{sup +4.1}{sub -2.9} Multiplication-Sign 10{sup -5} Mpc{sup -3} for objects with a de-absorbed luminosity larger than 2 Multiplication-Sign 10{sup 42} erg s{sup -1}. As the BAT AGNs are all mostly local, they allow us to investigate the spatial distribution of AGNs in the nearby universe regardless of absorption. We find concentrations of AGNs that coincide spatially with the largest congregations of matter in the local ({<=}85 Mpc) universe. There is some evidence that the fraction of Seyfert 2 objects is larger than average in the direction of these dense regions.

79 ASTRONOMY AND ASTROPHYSICS↗

The correlation between metallicity and debris disk mass

We find that the initial dust masses in planetary debris disks are correlated with the metallicities of their central stars. We compiled a large sample of systems, including Spitzer, the Herschel DUNES and DEBRIS surveys, and WISE debris disk candidates. We also merged 33 metallicity catalogs to provide homogeneous [Fe/H] and σ{sub [Fe/H]} values. We analyzed this merged sample, including 222 detected disks (74 warm and 148 cold) around a total of 187 systems (some with multiple components) and 440 disks with only upper limits (125 warm and 315 cold) around a total of 360 systems. The disk dust masses at a common early evolutionary point in time were determined using our numerical disk evolutionary code, evolving a unique model for each of the 662 disks backward to an age of 1 Myr. We find that disk-bearing stars seldom have metallicities less than [Fe/H]=−0.2 and that the distribution of warm component masses lacks examples with large mass around stars of low metallicity ([Fe/H]<−0.085). Previous efforts to find a correlation have been largely unsuccessful; the primary improvements supporting our result are (1) basing the study on dust masses, not just infrared excess detections; (2) including upper limits on dust mass in a quantitative way; (3) accounting for the evolution of debris disk excesses as systems age; (4) accounting fully for the range of uncertainties in metallicity measurements; and (5) having a statistically large enough sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of sampling plan options for tank 16H from the perspective of statistical uncertainty

This report develops a concentration variability model for Tank 16H in order to compare candidate sampling plans for assessing the concentrations of analytes in the residual material in the annulus and on the floor of the primary vessel. A concentration variability model is used to compare candidate sampling plans based on the expected upper 95% confidence limit (UCL95) for the mean. The result is expressed as a rank order of candidate sampling plans from lowest to highest expected UCL95, with the lowest being the most desirable from an uncertainty perspective.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Analysis of Biomolecular Force Fields for Simulations of Polyglutamine in Solution

Polyglutamine (polyQ) peptides are a useful model system for biophysical studies of protein folding and aggregation, both for their intriguing aggregation properties and their own relevance to human disease. The genetic expansion of a polyQ tract triggers the formation of amyloid aggregates associated with nine neurodegenerative diseases. Several clearly identifiable and separable factors, notably the length of the polyQ tract, influence the mechanism of aggregation, its associated kinetics, and the ensemble of structures formed. Atomistic simulations are well positioned to answer open questions regarding the thermodynamics and kinetics of polyQ folding and aggregation. The additional, explicit representation of water permits deeper investigation of the role of solvent dynamics, and it permits a direct comparison of simulation results with infrared spectroscopy experiments. The generation of meaningful simulation results hinges on satisfying two essential criteria: achieving sufficient conformational sampling to draw statistically valid conclusions, and accurately reproducing the intermolecular forces that govern system structure and dynamics. In this work, we examine the ability of 12 biomolecular force fields to reproduce the properties of a simple, 30-residue polyQ peptide (Q30) in explicit water. In addition to secondary and tertiary structure, we consider generic structural properties of polymers that provide additional dimensions for analysis of the highly degenerate disordered states of the molecule. We find that the 12 force fields produce a wide range of predictions. We identify AMBER ff99SB, AMBER ff99SB*, and OPLS-AA/L to be most suitable for studies of polyQ folding and aggregation.

97 MATHEMATICS AND COMPUTING↗

Modeling Freedom From Progression for Standard-Risk Medulloblastoma: A Mathematical Tumor Control Model With Multiple Modes of Failure

Purpose: As pediatric medulloblastoma (MB) is a relatively rare disease, it is important to extract the maximum information from trials and cohort studies. Here, a framework was developed for modeling tumor control with multiple modes of failure and time-to-progression for standard-risk MB, using published pattern of failure data. Methods and Materials: Outcome data for standard-risk MB published after 1990 with pattern of relapse information were used to fit a tumor control dose-response model addressing failures in both the high-dose boost volume and the elective craniospinal volume. Estimates of 5-year event-free survival from 2 large randomized MB trials were used to model the time-to-progression distribution. Uncertainty in freedom from progression (FFP) was estimated by Monte Carlo sampling over the statistical uncertainty in input data. Results: The estimated 5-year FFP (95% confidence intervals [CI]) for craniospinal doses of 15, 18, 24, and 36 Gy while maintaining 54 Gy to the posterior fossa was 77% (95% CI, 70%-81%), 78% (95% CI, 73%-81%), 79% (95% CI, 76%-82%), and 80% (95% CI, 77%-84%) respectively. The uncertainty in FFP was considerably larger for craniospinal doses below 18 Gy, reflecting the lack of data in the lower dose range. Conclusions: Estimates of tumor control and time-to-progression for standard-risk MB provides a data-driven setting for hypothesis generation or power calculations for prospective trials, taking the uncertainties into account. The presented methods can also be applied to incorporate further risk-stratification for example based on molecular biomarkers, when the necessary data become available.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Serum Proteome Signature of Radiation Response: Upregulation of Inflammation-Related Factors and Downregulation of Apolipoproteins and Coagulation Factors in Cancer Patients Treated With Radiation Therapy—A Pilot Study

Purpose: Ionizing radiation affects the proteome of irradiated cells and tissue, yet data concerning changes induced during radiation therapy (RT) in human blood are fragmentary and inconclusive. We aimed to identify features of serum proteome and associated processes involved in response to partial body irradiation during cancer treatment. Methods and Materials: Twenty patients with head and neck squamous cell cancer (HNSCC) and 20 patients with prostate cancer received definitive intensity modulated RT. Blood samples were collected before RT, just after RT, and 1 month after the end of RT. Complete serum proteome was analyzed in individual samples, using a shotgun liquid chromatography-tandem mass spectrometry approach which allowed identification of approximately 450 proteins. Approximately 100 unique proteins were quantified in all samples after exclusion of immunoglobulins, and statistical significance of differences among consecutive samples was assessed. Processes associated with quantified proteins and their functional interactions were predicted using gene ontology tools. Results: RT-induced changes were marked in the HNSCC patient group: 22 upregulated and 33 downregulated proteins were detected in post-RT sera. Most of the changes reversed during follow-up, yet levels of some proteins remained affected 1 month after the end of RT. RT-upregulated proteins were associated with acute phase, inflammatory response, and complement activation. RT-downregulated proteins were associated with transport and metabolism of lipids (plasma apolipoproteins) and blood coagulation. RT-induced changes were much weaker in prostate cancer patients, which corresponded to differences in acute radiation toxicity observed in both groups. Nevertheless, general patterns of RT-induced sera proteome changes were similar in both of the groups of cancer patients. Conclusions: In this pilot study, we proposed to identify a molecular signature of radiation response, based on specific features of serum proteome. The signature included upregulation of factors involved in acute or inflammatory response but also downregulation of plasma apolipoproteins and factors involved in blood coagulation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Compressive sampling of polynomial chaos expansions: Convergence analysis and sampling strategies

Sampling orthogonal polynomial bases via Monte Carlo is of interest for uncertainty quantification of models with random inputs, using Polynomial Chaos (PC) expansions. It is known that bounding a probabilistic parameter, referred to as coherence, yields a bound on the number of samples necessary to identify coefficients in a sparse PC expansion via solution to an ℓ{sub 1}-minimization problem. Utilizing results for orthogonal polynomials, we bound the coherence parameter for polynomials of Hermite and Legendre type under their respective natural sampling distribution. In both polynomial bases we identify an importance sampling distribution which yields a bound with weaker dependence on the order of the approximation. For more general orthonormal bases, we propose the coherence-optimal sampling: a Markov Chain Monte Carlo sampling, which directly uses the basis functions under consideration to achieve a statistical optimality among all sampling schemes with identical support. We demonstrate these different sampling strategies numerically in both high-order and high-dimensional, manufactured PC expansions. In addition, the quality of each sampling method is compared in the identification of solutions to two differential equations, one with a high-dimensional random input and the other with a high-order PC expansion. In both cases, the coherence-optimal sampling scheme leads to similar or considerably improved accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dosimetric comparison of rectal-sparing capabilities of rectal balloon vs injectable spacer gel in stereotactic body radiation therapy for prostate cancer: lessons learned from prospective trials

This study aimed to compare the rectal-sparing capabilities of rectal balloons vs absorbable injectable spacer gel in stereotactic body radiation therapy (SBRT) for prostate cancer. Patient samples included in this analysis were obtained from 2 multi-institutional prospective trials of SBRT for prostate cancer using a rectal balloon (n = 36 patients) and injectable spacer gel (n = 36). Treatment prescription dose was 45 Gy in 5 fractions in 42 patients; for equal comparison, the remaining 30 patients were rescaled to 45 Gy from 47.5 Gy prescription (n = 6) and 50 Gy prescription (n = 24). The median prostate volumes and body mass index in the 2 patient samples were not statistically significantly different (p = 0.67 and 0.45, respectively), supporting anatomic similarity between cohorts. The injectable spacer gel achieved dosimetric superiority over the rectal balloon with respect to the maximum dose to the rectum (42.3 vs 46.2 Gy, p vs 35.1 Gy, p Gy ), V40{sub Gy}, and V30{sub Gy} (0.3 vs 1.7 cc, 1 vs 5.4 cc, and 4.1 vs 9.6 cc, respectively; p Gy of the rectum or the dose to 50% of the rectal circumference (p = 0.29 and 0.06, respectively). The V18.3{sub Gy} of the bladder was significantly larger with the rectal balloon (19.9 vs 14.5 cc, p = 0.003). In this analysis of patients enrolled on 2 consecutive multi-institutional prospective trials of SBRT for prostate cancer, the injectable spacer gel outperformed the rectal balloon in the majority of the examined and relevant dosimetric rectal-sparing parameters. The rectal balloon did not outperform the injectable spacer gel in any measured rectal dose parameter.

61 RADIATION PROTECTION AND DOSIMETRY↗

A novel in situ trigger combination method

Searches for rare physics processes using particle detectors in high-luminosity colliding hadronic beam environments require the use of multi-level trigger systems to reject colossal background rates in real time. In analyses like the search for the Higgs boson, there is a need to maximize the signal acceptance by combining multiple different trigger chains when forming the offline data sample. In such statistically limited searches, datasets are often amassed over periods of several years, during which the trigger characteristics evolve and system performance can vary significantly. Reliable production cross-section measurements and upper limits must take into account a detailed understanding of the effective trigger inefficiency for every selected candidate event. We present as an example the complex situation of three trigger chains, based on missing energy and jet energy, that were combined in the context of the search for the Higgs (H) boson produced in association with a $W$ boson at the Collider Detector at Fermilab (CDF). We briefly review the existing techniques for combining triggers, namely the inclusion, division, and exclusion methods. We introduce and describe a novel fourth in situ method whereby, for each candidate event, only the trigger chain with the highest a priori probability of selecting the event is considered. We compare the inclusion and novel in situ methods for signal event yields in the CDF $WH$ search. This new combination method, by virtue of its scalability to large numbers of differing trigger chains and insensitivity to correlations between triggers, will benefit future long-running collider experiments, including those currently operating on the Large Hadron Collider.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

ParNCL and ParGAL: Data-parallel Tools for Postprocessing of Large-scale Earth Science Data

Earth science high-performance applications often require extensive analysis of their output in order to complete the scien- tific goals or produce a visual image or animation. Often this analysis cannot be done in situ because it requires calculating time-series statistics from state sampled over the entire length of the run or analyzing the relationship between similar time series from previous simulations or observations. Many of the tools used for this postprocessing are not themselves high- performance applications, but the new Parallel Gridded Analysis Library (ParGAL) provides high-performance data-parallel versions of several common analysis algorithms for data from a structured or unstructured grid simulation. The library builds on several scalable systems, including the Mesh Oriented DataBase (MOAB), a library for representing mesh data that sup- ports structured, unstructured finite element, and polyhedral grids; the Parallel-NetCDF (PNetCDF) library; and Intrepid, an extensible library for computing operators (such as gradient, curl, and divergence) acting on discretized fields. We have used ParGAL to implement a parallel version of the NCAR Command Language (NCL) a scripting language widely used in the climate community for analysis and visualization. The data-parallel algorithms in ParGAL/ParNCL are both higher performing and more flexible than their serial counterparts.

97 MATHEMATICS AND COMPUTING↗

Evidence of Skewness and Sub-Gaussian Character in Temperature-Dependent Distributions of One Million Electronic Excitation Energies in PbS Quantum Dots

Obtaining statistical distributions by sampling a large number of conformations is vital for an accurate description of temperature-dependent properties of chemical systems. However, constructing distributions with 105??106 samples is computationally challenging because of the prohibitively high computational cost of performing rst-principles quantum-mechanical calculations. In this work, we present a new technique called the Eective Stochastic Potential Conguration Interaction Singles (ESP-CIS) method to obtain excitation energies. The ESP-CIS method uses random matrix theory for the construction of an eective stochastic representation of the Fock operator and combines it with linear-response theory. Here, we perform excitation calculations on PbS quantum dots (0.75-1.75 nm) at temperatures of 200-400 K. Results from a total of 15 million excitation energy calculations (1 million samples for each of the three dots at ve dierent temperatures) revealed the distributions to be sub-Gaussian in nature with negative skewness. The distributions progressively became red-shifted with increasing temperature. The study shows that random distortions in structure from the minimum-energy geometry of these quantum dots are more likely to lower the excitation energies.

Scher, Jeremy A.↗

Dependence of the critical temperature in overdoped copper oxides on superfluid density

The physics of underdoped copper-oxide superconductors, including the pseudogap, spin and charge ordering, and their relation to superconductivity 1-3 , is intensely debated. The overdoped side is perceived as simpler, with strongly-correlated fermion physics evolving smoothly into the conventional Bardeen-Cooper-Schrieffer (BCS) behavior. Pioneering studies on a few overdoped samples 4-11 indicated that the superfluid density was much smaller than expected, but this was attributed to pair-breaking, disorder, and phase separation. Here, we test this conjecture by studying how the magnetic penetration depth λ and the phase stiffness ρs depend on temperature and doping, scanning densely the entire overdoped side of the La 2-x Sr x CuO 4 (LSCO) phase diagram. We have measured the absolute values of λ and ρs to the accuracy of ±1% in thousands of cuprate samples; the large statistics reveals clear trends and intrinsic properties. The films are quite homogeneous; variations in the critical temperature (T c ) within a film are very small (< 1 K). At every doping, ρs(T) decreases linearly with temperature. The T c (ρ s0) dependence is linear but with an offset, (T c - T 0 ) ∝ ρs0 where T0 ≈ 7 K, except very close to the origin where Tc ∝ √ρ s0. This scaling law defies the standard BCS description, posing a challenge to theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Delayed Slater determinant update algorithms for high efficiency quantum Monte Carlo

Within ab initio Quantum Monte Carlo simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunction. Each Monte Carlo step requires finding the determinant of a dense matrix. This is most commonly iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the inverse. The overall computational cost is therefore formally cubic in the number of electrons or matrix size. To improve the numerical efficiency of this procedure, we propose a novel multiple rank delayed update scheme. This strategy enables probability evaluation with application of accepted moves to the matrices delayed until after a predetermined number of moves, K . The accepted events are then applied to the matrices en bloc with enhanced arithmetic intensity and computational efficiency via matrix-matrix operations instead of matrix-vector operations. Here this procedure does not change the underlying Monte Carlo sampling or its statistical efficiency. For calculations on large systems and algorithms such as diffusion Monte Carlo where the acceptance ratio is high, order of magnitude improvements in the update time can be obtained on both multi- core CPUs and GPUs.

97 MATHEMATICS AND COMPUTING↗

CHARACTERIZING THE FORMATION HISTORY OF MILKY WAY LIKE STELLAR HALOS WITH MODEL EMULATORS

We use the semi-analytic model ChemTreeN, coupled to cosmological N-body simulations, to explore how different galaxy formation histories can affect observational properties of Milky Way like galaxies' stellar halos and their satellite populations. Gaussian processes are used to generate model emulators that allow one to statistically estimate a desired set of model outputs at any location of a p-dimensional input parameter space. This enables one to explore the full input parameter space orders of magnitude faster than could be done otherwise. Using mock observational data sets generated by ChemTreeN itself, we show that it is possible to successfully recover the input parameter vectors used to generate the mock observables if the merger history of the host halo is known. However, our results indicate that for a given observational data set, the determination of 'best-fit' parameters is highly susceptible to the particular merger history of the host. Very different halo merger histories can reproduce the same observational data set, if the 'best-fit' parameters are allowed to vary from history to history. Thus, attempts to characterize the formation history of the Milky Way using these kind of techniques must be performed statistically, analyzing large samples of high-resolution N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

WEAK HARD X-RAY EMISSION FROM TWO BROAD ABSORPTION LINE QUASARS OBSERVED WITH NuSTAR: COMPTON-THICK ABSORPTION OR INTRINSIC X-RAY WEAKNESS?

We present Nuclear Spectroscopic Telescope Array (NuSTAR) hard X-ray observations of two X-ray weak broad absorption line (BAL) quasars, PG 1004+130 (radio loud) and PG 1700+518 (radio quiet). Many BAL quasars appear X-ray weak, probably due to absorption by the shielding gas between the nucleus and the accretion-disk wind. The two targets are among the optically brightest BAL quasars, yet they are known to be significantly X-ray weak at rest-frame 2-10 keV (16-120 times fainter than typical quasars). We would expect to obtain Almost-Equal-To 400-600 hard X-ray ({approx}> 10 keV) photons with NuSTAR, provided that these photons are not significantly absorbed (N{sub H} {approx}< 10{sup 24} cm{sup -2}). However, both BAL quasars are only detected in the softer NuSTAR bands (e.g., 4-20 keV) but not in its harder bands (e.g., 20-30 keV), suggesting that either the shielding gas is highly Compton-thick or the two targets are intrinsically X-ray weak. We constrain the column densities for both to be N{sub H} Almost-Equal-To 7 Multiplication-Sign 10{sup 24} cm{sup -2} if the weak hard X-ray emission is caused by obscuration from the shielding gas. We discuss a few possibilities for how PG 1004+130 could have Compton-thick shielding gas without strong Fe K{alpha} line emission; dilution from jet-linked X-ray emission is one likely explanation. We also discuss the intrinsic X-ray weakness scenario based on a coronal-quenching model relevant to the shielding gas and disk wind of BAL quasars. Motivated by our NuSTAR results, we perform a Chandra stacking analysis with the Large Bright Quasar Survey BAL quasar sample and place statistical constraints upon the fraction of intrinsically X-ray weak BAL quasars; this fraction is likely 17%-40%.

79 ASTRONOMY AND ASTROPHYSICS↗