Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

How Many Elements Matter?

Some studies of stars' multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and $\mathrm{log}\,g$ minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [X/H] ≲ 0.02 dex, consistent with APOGEE's reported statistical uncertainties of ~0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

79 ASTRONOMY AND ASTROPHYSICS↗

Single molecule characterization of anomalous transport in a thin, anisotropic film

Here, the diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayers (PEMs) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film’s matrix.

36 MATERIALS SCIENCE↗

A numerical method for relating two- and three-dimensional pressure distributions on transonic wings

This paper presents a preliminary design method for determining a wing's design pressure distribution and geometry based on airfoil normal Mach numbers and airfoil loading. In this method, the perturbation velocities in supercritical regions are computed from airfoil transonic normal Mach numbers and include the influence of local sweep, taper, and three-dimensional induced velocities, so that the appearance and strength of shocks can be expected to resemble those of the airfoil. The velocities in subcritical wing regions are scaled first with simple sweep theory, and then to achieve the desired load distribution. The method was applied to the design of an oblique flying wing, using a linear potential method. The required wing area could be reduced by 14 percent using this method rather than simple sweep theory.

Kroo, Ilan↗

Microstructural Effects on SCC Initiation PWR Primary Water Cold-Worked Alloy 600

SCC initiation behavior of one mill annealed alloy 600 plate heat was investigated in simulated PWR primary water under constant load at yield stress with in-situ direct current potential drop (DCPD) monitoring for crack initiation. Twelve specimens were tested at similar cold work levels among which three showed much shorter SCC initiation times (<400 hrs) than the others (>1200 hrs). Post-test examinations revealed that these three specimens all feature an inhomogeneous microstructure where the primary crack always nucleated along the boundary of large elongated grains protruding normally into the gauge. In contrast, such microstructure was either not observed or did not extend deep enough into the gauge in the other specimens exhibiting ~3-6X longer initiation times. In order to better understand the role of this microstructural inhomogeneity in SCC initiation, high-resolution microscopy was performed to compare carbide morphology and strain distribution between the long grains and normal grains, and their potential effects on SCC initiation are discussed in this paper.

Alloy 600↗

Use of collateral information to improve LANDSAT classification accuracies

Methods to improve LANDSAT classification accuracies were investigated including: (1) the use of prior probabilities in maximum likelihood classification as a methodology to integrate discrete collateral data with continuously measured image density variables; (2) the use of the logit classifier as an alternative to multivariate normal classification that permits mixing both continuous and categorical variables in a single model and fits empirical distributions of observations more closely than the multivariate normal density function; and (3) the use of collateral data in a geographic information system as exercised to model a desired output information layer as a function of input layers of raster format collateral and image data base layers.

Strahler, A. H.↗

Comprehensive study of the X-ray flares from gamma-ray bursts observed by Swift

X-ray flares are generally supposed to be produced by later activities of the central engine, and may share a similar physical origin with the prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by Swift from 2005 April to 2015 March. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior L{sub p}∝T{sub peak,z}{sup −1.27}. The flare duration increases with peak time. The 0.3–10 keV isotropic energy of the distribution of X-ray flares is a log-normal peaked at 10{sup 51.2} erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times, and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models based on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be dominated by Poynting flux.

79 ASTRONOMY AND ASTROPHYSICS↗

An analytical elastic plastic contact model with strain hardening and frictional effects for normal and oblique impacts

Impact between metallic surfaces is a phenomenon that is ubiquitous in the design and analysis of mechanical systems. We found that to model this phenomenon, a new formulation for frictional elastic–plastic contact between two surfaces is developed. The formulation is developed to consider both frictional, oblique contact (of which normal, frictionless contact is a limiting case) and strain hardening effects. The constitutive model for normal contact is developed as two contiguous loading domains: the elastic regime and a transitionary region in which the plastic response of the materials develops and the elastic response abates. For unloading, the constitutive model is based on an elastic process. Moreover, the normal contact model is assumed to only couple one-way with the frictional/tangential contact model, which results in the normal contact model being independent of the frictional effects. Frictional, tangential contact is modeled using a microslip model that is developed to consider the pressure distribution that develops from the elastic–plastic normal contact. This model is validated through comparisons with experimental results reported in the literature, and is demonstrated to be significantly more accurate than 10 other normal contact models and three other tangential contact models found in the literature.

36 MATERIALS SCIENCE↗

Monte Carlo lattice models for adsorbed polymer conformation

The adhesion between a polymer film and a metal surface is of great technological interest. However, the prediction of adhesion and wear properties of polymer coated metals is quite difficult because a fundamental understanding of the polymer surface interaction does not yet exist. A computer model for the conformation of a polymer molecule adsorbed on a surface is discussed. The chain conformation is assumed to be described by a partially directed random walk on a three dimensional simple cubic lattice. An attractive surface potential is incorporated into the model through the use of a random walk step probability distribution that is anisotropic in the direction normal to the attractive surface. The effects of variations in potential characteristics are qualitatively included by varying both the degree of anisotropy of the step distribution and the range of the anisotropy. Polymer conformation is characterized by the average end to end distance, average radius of gyration, and average number of chain segments adsorbed on the surface.

Good, B. S.↗

Rates and properties of Type Ia supernovae in galaxy clusters within the dark energy survey

We identify 66 photometrically classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) that have occurred within red-sequence selected galaxy clusters. We compare light-curve and host galaxy properties of the cluster SNe to 1024 DES SNe Ia located in field galaxies, the largest comparison of two such samples at high redshift (z > 0.1). We find that cluster SN light curves decline faster than those in the field (97.7 per cent confidence). However, when limiting these samples to host galaxies of similar colour and mass, there is no significant difference in the SN light-curve properties. Motivated by previous detections of a higher-normalized SN Ia delay-time distribution in galaxy clusters, we measure the intrinsic rate of SNe Ia in cluster and field environments. We find the average ratio of the SN Ia rate per galaxy between high-mass (⁠|$10\le \log \mathrm{(\mathit{ M}_{*}/{\rm \mathit{ M}}_{\odot })} \le 11.25$|⁠) cluster and field galaxies to be 0.594 ± 0.068. This difference is mass-dependent, with the ratio declining with increasing mass, which suggests that the stellar populations in cluster hosts are older than those in field hosts. We show that the mass-normalized rate (or SNe per unit mass) in massive–passive galaxies is consistent between cluster and field environments. Additionally, both of these rates are consistent with rates previously measured in clusters at similar redshifts. We conclude that in massive–passive galaxies, which are the dominant hosts of cluster SNe, the cluster delay-time distribution is comparable to the field.

79 ASTRONOMY AND ASTROPHYSICS↗

Further Investigation of Employing Binomial Small-Sample Estimator in M-LROM

Reduced Order Modeling (ROM) is essential to enable the identification of the so-called active parameter subspace for high-fidelity models that are employed in computationally intensive analyses such as uncertainty characterization and data assimilation. Given the complexity of models for real-world application, brute force methods to identify the active subspace may not be computationally efficient to realize the value of reduction, since in realistic problems, the full scale model must be executed few hundred times, which may be computationally taxing. Earlier work has introduced a new approach, which represents a further development of our previous work on ROM, referred to as multi-level ROM (M-LROM). In M-LROM, approximate models are employed to identify the active subspace instead of the full scale model. In this summary, we present our recent results on M-LROM, specifically in regard to the construction of reliable error metrics that can bound the errors resulting from the reduction. The idea is to employ a lower fidelity model, i.e., an approximation of the full scale or high fidelity model, to build the active subspace, and via the proposed error bound, the quality of the active subspace can be determined by calculating an upper bound on the errors resulting from the reduction. For this approach to be practical, the error bound must be realistic, implying that it should be close to the actual errors expected, and not much higher. Therefore, an important goal of this summary is to compare different error bounds in terms of the relation between their estimated bounds to the actual errors. In particular, we compare the proposed error bound which numerically construct an estimator using binomial distribution to the widely used small-sample norm estimators. Unlike most probabilistic error estimation methods, the proposed bound employs a binomial distribution to construct the estimator whereas in other conventional methods, the distribution employed is typically the standard Gaussian, normalized over a unit hypersphere as initially proposed by the seminal paper of J. Dixon on random matrix theory. Previous publications provided evidence that the binomial distribution calculates more realistic bounds as compared to the Gaussian distribution, which tends to be conservative, implying that the bound is much bigger than the actual errors, which diminishes the value of the reduction, as now must increase the size of the active subspace to reach the desired upper bound on the reduction errors. This work has provided a rigorous approach by which error bounds on the reduction errors can be established without having to execute the high fidelity model during the construction of the active subspace. The high fidelity model is only need to test the error bound using small number of executions. This methodology paves the way for a more efficient application of ROM to complex engineering models, where standard ROM techniques may be too taxing to optimally realize their value in practical situations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparative Prospective Study of Load Distribution Projection Among Patients with Vertebral Fractures Treated with Percutaneous Vertebroplasty and a Control Group of Healthy Volunteers

Purpose: Through a prospective comparison of patients with vertebral fractures and normal population, we illustrate effect of percutaneous vertebroplasty (PV) upon projection of load distribution changes. Methods: Vertebroplasty group (36 symptomatic patients with osteoporotic vertebral fractures) was evaluated on an electronic baropodometer registering projection of weight bearing areas on feet. Load distribution between right and left foot (including rear-front of the same foot) during standing and walking was recorded and compared before (group V1) and the day after (group V2) PV. Control group (30 healthy asymptomatic volunteers-no surgery record) were evaluated on the same baropodometer. Results: Mean value of load distribution difference between rear-front of the same foot was 9.45 ± 6.79 % (54.72–45.28 %) upon standing and 14.76 ± 7.09 % (57.38–42.62 %) upon walking in the control group. Respective load distribution values before PV were 16.52 ± 11.23 and 30.91 ± 19.26 % and after PV were 10.08 ± 6.26 and 14.25 ± 7.68 % upon standing and walking respectively. Mean value of load distribution variation between the two feet was 6.36 and 14.6 % before and 4.62 and 10.4 % after PV upon standing and walking respectively. Comparison of load distribution variation (group V1–V2, group V1-control group) is statistically significant. Comparison of load distribution variation (group V2-control group) is not statistically significant. Comparison of load distribution variation among the two feet is statistically significant during walking but not statistically significant during standing. Conclusions: There is a statistically significant difference when comparing load distribution variation prior vertebroplasty and that of normal population. After vertebroplasty, this difference normalizes in a statistically significant way. PV is efficient on equilibrium-load distribution improvement as well.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Effects of the guide field on electron distribution functions in the diffusion region of asymmetric reconnection

Crescent electron distribution functions in the electron diffusion region in asymmetric reconnection with a guide field (strength less than 50% of the upstream magnetic field) are investigated by theory and fully kinetic simulations. Electron motion in the electron diffusion region is a combination of meandering across the current sheet and gyration around the guide field. Combining the two canonical momentum conservation laws in a one-dimensional model with the energy conservation law, we derive the shape of crescent electron distribution functions in both velocity planes: perpendicular to the magnetic field and coplanner with the magnetic field. There are two major effects of the guide field on crescent distribution functions: (1) widening the opening angle of the crescent in a reduced distribution function and (2) the linear cut-off of a distribution function nonparallel to the magnetic field. We also discuss a translational mapping (perpendicular to the linear cut-off) of distribution functions between the X-line and another position. In addition, considering the gyration due to the normal magnetic field, we predict the shape of distribution functions in general positions in the electron diffusion region. Finally, the predictions will be useful to identify features of magnetic reconnection in a guide field by satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cumulative slant path rain attenuation associated with COMSTAR beacon at 28.56 GHz for Wallops Island, Virginia

Yearly, monthly, and time of day fade statistics are presented and characterized. A 19.04 GHz yearly fade distribution, corresponding to a second COMSTAR beacon frequency, is predicted using the concept of effective path length, disdrometer, and rain rate results. The yearly attenuation and rain rate distributions follow with good approximation log normal variations for most fade and rain rate levels. Attenuations were exceeded for the longest and shortest periods of times for all fades in August and February, respectively. The eight hour time period showing the maximum and minimum number of minutes over the year for which fades exceeded 12 db were approximately between 1600 to 2400, and 0400 to 1200 hours, respectively. In employing the predictive method for obtaining the 19.04 GHz fade distribution, it is demonstrated theoretically that the ratio of attenuations at two frequencies is minimally dependent of raindrop size distribution providing these frequencies are not widely separated.

Goldhirsh, J.↗

SU-E-T-227: Least Restrictive Assignment of Dose-Distance Differences (LRAD) in IMRT On a Conventional Linac

Purpose: To develop a method to compare dose distributions based on the least restrictive of dose and distance differences, and to evaluate the method for IMRT QA on a conventional linac. Methods: Lrad is defined as the least restrictive of the dose difference at each point in the distributions being compared (%D), with dose normalized by the same factor, and the distance to agreement (DTA in millimeters), limited to directions parallel and perpendicular to MLC motion. If %D<%Do, Lrad=%D. Otherwise, if L=%D/DTA>Lo, Lrad=DTA, else Lrad=%D. Lrad fails at any point when Lrad>DTAo or Lrad>%Do. We compared LRAD (Lo=1%/mm, DTAo=2mm, %Do=2%) with Gamma index tool (2%/2mm pass/fail criteria) for IMRT QA with 2D detector arrays. Dose distribution measured for 6 MV and 18 MV x-rays for various field sizes and shapes, including 6 IMRT plans, were compared to measurements made with output, beam flatness, beam symmetry, MLC leaf positions, and detector array translation and rotation adjusted at up to twice published acceptable tolerance. Results: LRAD was more effective than Gamma in detecting a 2% drift in output (dose per monitor unit), 1–2 mm drifts in MLC leaf position for 1–8 leaves in the field, and beam flatness and symmetry drifts that resulted in a 3% change in the off-axis ratio. LRAD and Gamma were equally effective for larger output changes and for detector translations and rotations. For all ∼400 measurements performed, LRAD detected 78% of the error and Gamma detected 52% of the errors. By using a color scale to display dose differences and grey scale for spatial errors, the reason for failure was self evident, whether it be beam or mechanical changes. Conclusion: LRAD proved to be advantageous in its ability to detect and distinguish between realistic changes in beam and mechanical settings in comparison to the Gamma Index.

60 APPLIED LIFE SCIENCES↗

Modeling energy balance while correcting for measurement error via free knot splines

Measurements of energy balance components (energy intake, energy expenditure, changes in energy stores) are often plagued with measurement error. Doubly-labeled water can measure energy intake (EI) with negligible error, but is expensive and cumbersome. An alternative approach that is gaining popularity is to use the energy balance principle, by measuring energy expenditure (EE) and change in energy stores (ES) and then back-calculate EI. Gold standard methods for EE and ES exist and are known to give accurate measurements, albeit at a high cost. We propose a joint statistical model to assess the measurement error in cheaper, non-intrusive measures of EE and ES. We let the unknown true EE and ES for individuals be latent variables, and model them using a bivariate distribution. We try both a bivariate Normal as well as a Dirichlet Process Mixture Model, and compare the results via simulation. Our approach, is the first to account for the dependencies that exist in individuals’ daily EE and ES. We employ semiparametric regression with free knot splines for measurements with error, and linear components for error free covariates. We adopt a Bayesian approach to estimation and inference and use Reversible Jump Markov Chain Monte Carlo to generate draws from the posterior distribution. Based on the semiparameteric regression, we develop a calibration equation that adjusts a cheaper, less reliable estimate, closer to the true value. Along with this calibrated value, our method also gives credible intervals to assess uncertainty. A simulation study shows our calibration helps produce a more accurate estimate. Furthermore, our approach compares favorably in terms of prediction to other commonly used models.

60 APPLIED LIFE SCIENCES↗

Communication: Ultraviolet photodissociation dynamics of the simplest Criegee intermediate CH{sub 2}OO

The velocity and angular distributions of O {sup 1}D photofragments arising from UV excitation of the CH{sub 2}OO intermediate on the B {sup 1}A′←X {sup 1}A′ transition are characterized using velocity map ion imaging. The anisotropic angular distribution yields the orientation of the transition dipole moment, which reflects the π*←π character of the electronic transition associated with the COO group. The total kinetic energy release distributions obtained at several photolysis wavelengths provide detail on the internal energy distribution of the formaldehyde cofragments and the dissociation energy of CH{sub 2}OO X {sup 1}A′ to O {sup 1}D + H{sub 2}CO X {sup 1}A{sub 1}. A common termination of the total kinetic energy distributions, after accounting for the different excitation energies, gives an upper limit for the CH{sub 2}OO X {sup 1}A′ dissociation energy of D{sub 0}≤ 54 kcal mol{sup −1}, which is compared with theoretical predictions including high level multi-reference ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak (R = 0.9 m, B t = 5 T, and I p = 2 MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric fields and secondary emission near a dielectric-metal interface

Dielectric surface-charge distributions near a metal-dielectric interface in vacuum depend upon secondary emission processes in the presence of normal and tangential components of electric field. From measured charge distributions created by exposing a specimen of fluorinated ethylene-propylene to monoenergetic electron fluxes, it has been possible to calculate potentials and fields on and near the dielectric surface. The effect of the normal electric field upon secondary emission is measured directly, and the effect of the tangential field is inferred from the charge-distribution data. The critical point (unity crossover) for secondary emission is shifted by the application of fields, so that it occurs at much higher primary energies than normally. Primary beams having energies up to 20 keV are used, and surface fields are as high as 20 kV/mm.

Robinson, J. W.↗