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,117 records · Page 62

A new portable random number generator wrapper library

Random number generator is an important component of many scientific projects. Many projects are written using programming models (like OpenMP and SYCL) to target different architectures. However, some programming models do not provide a random number generator. In this work, we introduce our random number generator wrapper. It is a header-only library that supports three distributions of random numbers: uniform, normal, and poisson. On the GPU backend, it wraps the cuRAND and rocRAND library, and supports various random number engines. It also wraps random123, a counterbased random number generator, on both CPU and GPU. With this library, we can generate random numbers with a few lines of code and target both GPU and multi-thread CPU with the same code. We also investigate the performance and scalability of this wrapper on different architectures with different engines and the number of cores.

97 MATHEMATICS AND COMPUTING↗

Non-local thermodynamic equilibrium 1.5D modeling of red giant stars

Spectra for two-dimensional (2D) stars in the 1.5D approximation are created from synthetic spectra of one-dimensional (1D) non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star while varying the temperature difference between the two 1D component models (ΔT {sub 1.5D}) and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred T {sub eff} values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the T {sub eff} of the 1.5D stars: UBVRI photometric colors, absolute surface flux spectral energy distributions (SEDs), relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred T {sub eff} value increases with increasing ΔT {sub 1.5D}. In all cases, the inferred T {sub eff} value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of ΔT {sub 1.5D} within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between T {sub eff} values derived from NLTE and LTE modeling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled and largely depends on the observable quantity being fit.

79 ASTRONOMY AND ASTROPHYSICS↗

Light-front spin-dependent spectral function and nucleon momentum distributions for a three-body system

Poincare covariant definitions for the spin-dependent spectral function and for the momentum distributions within the light-front Hamiltonian dynamics are proposed for a three-fermion bound system, starting from the light-front wave function of the system. The adopted approach is based on the Bakamjian–Thomas construction of the Poincaré generators, which allows one to easily import the familiar and wide knowledge on the nuclear interaction into a light-front framework. The proposed formalism can find useful applications in refined nuclear calculations, such as those needed for evaluating the European Muon Collaboration effect or the semi-inclusive deep inelastic cross sections with polarized nuclear targets, since remarkably the light-front unpolarized momentum distribution by definition fulfills both normalization and momentum sum rules. As a result, also shown is a straightforward generalization of the definition of the light-front spectral function to an A -nucleon system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of $\pi^-$ production in $^7$Be+$^9$Be collisions at beam momenta from 19$A$ to 150$A$GeV/$c$ in the NA61/SHINE experiment at the CERN SPS

The NA61/SHINE collaboration studies at the CERN Super Proton Synchrotron (SPS) the onset of deconfinement in hadronic matter by the measurement of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra and mean multiplicities of $\pi^{-}$ mesons produced in the 5% most central $^7$Be + $^9$Be collisions at beam momenta of 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$ obtained by the so-called $h^-$ method which does not require any particle identification. The shape of the transverse mass spectra differs from the shapes measured in central Pb + Pb collisions and inelastic p + p interactions. The normalized width of the rapidity distribution decreases with increasing collision energy and is in between the results for inelastic nucleon-nucleon and central Pb + Pb collisions. The mean multiplicity of pions per wounded nucleon in central $^7$Be + $^9$Be collisions is close to that in central Pb + Pb collisions up to 75$A$GeV/$c$. However, at the top SPS energy the result is close to the one for inelastic nucleon-nucleon interactions. The results are discussed in the context of predictions for the onset of deconfinement at the CERN SPS collision energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multilevel Sequential Monte Carlo Samplers for Normalizing Constants

This article considers the sequential Monte Carlo (SMC) approximation of ratios of normalizing constants associated to posterior distributions which in principle rely on continuum models. Therefore, the Monte Carlo estimation error and the discrete approximation error must be balanced. A multilevel strategy is utilized to substantially reduce the cost to obtain a given error level in the approximation as compared to standard estimators. Two estimators are considered and relative variance bounds are given. The theoretical results are numerically illustrated for two Bayesian inverse problems arising from elliptic partial differential equations (PDEs). The examples involve the inversion of observations of the solution of (i) a 1-dimensional Poisson equation to infer the diffusion coefficient, and (ii) a 2-dimensional Poisson equation to infer the external forcing.

97 MATHEMATICS AND COMPUTING↗

QuantumNAT: quantum noise-aware training with noise injection, quantization and normalization

Parameterized Quantum Circuits (PQC) are promising towards quantum advantage on near-term quantum hardware. However, due to the large quantum noises (errors), the performance of PQC models has a severe degradation on real quantum devices. Take Quantum Neural Network (QNN) as an example, the accuracy gap between noise-free simulation and noisy results on IBMQ-Yorktown for MNIST-4 classification is over 60%. Existing noise mitigation methods are general ones without leveraging unique characteristics of PQC; on the other hand, existing PQC work does not consider noise effect. To this end, we present QuantumNAT, a PQC-specific framework to perform noise-aware optimizations in both training and inference stages to improve robustness. We experimentally observe that the effect of quantum noise to PQC measurement outcome is a linear map from noise-free outcome with a scaling and a shift factor. Motivated by that, we propose post-measurement normalization to mitigate the feature distribution differences between noise-free and noisy scenarios. Furthermore, to improve the robustness against noise, we propose noise injection to the training process by inserting quantum error gates to PQC according to realistic noise models of quantum hardware. Lastly, post-measurement quantization is introduced to quantize the measurement outcomes to discrete values, achieving the denoising effect. Extensive experiments on 8 classification tasks using 6 quantum devices demonstrate that QuantumNAT improves accuracy by up to 43%, and achieves over 94% 2-class, 80% 4-class, and 34% 10-class classification accuracy measured on real quantum computers. The code for construction and noise-aware training of PQC is available in the TorchQuantum library.

Wang, Hanrui↗

Correction of confidence intervals in excess relative risk models using Monte Carlo dosimetry systems with shared errors

Introduction: Measurement errors are ubiquitous in epidemiological studies, especially for doses arising from environmental and occupational exposures, which are difficult to assess. Correction methods for independent measurement errors, have been comprehensively described elsewhere, including regression calibration, simulation extrapolation (SIMEX), etc. However, when exposure estimates are constructed using complicated physical and biological models, uncertainties in the dosimetry system can become very complex. In such systems, uncertainties in some parameters may affect a large group of study participants simultaneously, and hence are “shared” as opposed to independent dosimetric uncertainties. Representing the structure of the uncertainty in the dosimetry by providing multiple realizations of possible dose has been described for several studies. The goal of this paper is to implement, and explore the properties of, a generalized estimating equation-based approach to shared uncertainties based on principles described in an earlier paper. Methods: The effects of shared and unshared dosimetry error on parameter estimates and confidence intervals in a generalized excess relative risk model are examined. A post-hoc correction method is implemented that constructs corrected confidence intervals based on an underlying variance structure described in an earlier paper. In this correction method it is assumed that the estimated slope parameter in a linear excess relative risk model is distributed as a mixture of normal and lognormal components. A simulation study of the behavior of this correction was developed based on multiple realizations provided by the Monte Carlo dosimetry system used for the plutonium exposures in the Mayak Worker Cohort and a risk model with close similarity to that used in analyses of lung cancer data in that cohort. Results: In the simulation study we found that shared dose errors severely degraded the coverage probabilities of confidence limits for the slope parameter in the linear excess relative risk model relating uncertain exposure to lung cancer risk. Using our correction method to adjust the confidence intervals provided much improved coverage probabilities for this parameter. A smaller amount of degradation in confidence interval coverage for other parameters (e.g. those in the baseline and dose modification terms) was noted, and the proposed correction method also improved coverage for these parameters as well. Discussion: Compared to other recently described methods our proposed correction method relies on relatively familiar software to fit dose response models. The basic modeling philosophy described here first uses existing Poisson regression methods to estimate parameters in the excess relative risk model; in this first stage mean dose (computed from the multiple realizations) is used as the independent variable in the regression. Once the model is fit we then correct (widen) the standard errors for all parameters using the proposed correction method. The simulation studies indicate that the methods described here are an important contribution to the available tools for dealing with shared measurement errors in complex risk analysis.

61 RADIATION PROTECTION AND DOSIMETRY↗

Relativistic Alfvén Turbulence at Kinetic Scales

Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close to k −3 or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

Description and calibration of the Langley 6- by 19-inch transonic tunnel

A description and calibration is presented of the Langley 6- by 19-inch transonic tunnel which is a two-dimensional facility with top and bottom slotted walls used for testing two-dimensional airfoil sections. Basic tunnel-empty Mach number distributions and schlieren flow photographs as well as integrated normal-force coefficients, pitching-moment coefficients, surface-pressure distributions, and schlieren flow photographs of an NACA 0012 airfoil calibration model are presented. The Mach number capability of the facility is from 0.5 to about 1.1 with a corresponding Reynolds number range of 1.5 million to 3 million based on a 4.0-in. model chord. Comparisons of experimental results from the tests with previous data are also presented.

Ladson, C. L.↗

The optical oblateness of Mars

Two yellow and two blue photometric plates of Mars obtained on the near-opposition date of November 21/22, 1958 were studied with a microdensitometer. The plates were provided with a photometric tube scale. The measurements were converted into absolute normalized intensities with the aid of integrated magnitudes of the planet adopted from other sources. Exposures, coordinates and isophotic points are discussed together with the normalization of isophotic radii, standard distributions of brightness, reduced radii, atmospheric levels, and apparent oblateness.

Opik, E. J.↗

A high order theory for uniform and laminated plates

A theory of plate deformation is derived which accounts for the effects of transverse shear deformation, transverse normal strain, and a nonlinear distribution of the in-plane displacements with respect to the thickness coordinate. The theory is compared with lower order plate theories through application to a particular problem involving a plate acted upon by a sinusoidal surface pressure. Comparison is also made with exact elasticity solution of this problem. It is found that when the ratio of the characteristic length of the load pattern to the plate thickness is of the order of unity, lower order theories are inadequate and the present high order theory is required to give meaningful results. Results are given for the bending of symmetric cross-ply and angle-ply laminates. Comparison with exact elasticity solutions indicates that the present plate theory is sufficiently accurate for predicting the behavior of thick laminates.

Lo, K. H.↗

Analysis of delamination in unidirectional and crossplied fiber composites containing surface cracks

A two-dimensional hybrid stress finite element analysis is described which was used to study the local stress field around delamination cracks in composite materials. The analysis employs a crack tip singularity element which is embedded in a matrix interlayer between plies of the laminate. Results are given for a unidirectional graphite/epoxy laminate containing a delamination emanating from a surface crack through the outside ply. The results illustrate several aspects of delamination cracks: (1) the localization of the singular stress domain within the interlayer; (2) the local concentration of stress in the ply adjacent to the crack; (3) the nature of the transverse normal and interlaminar shear stress distributions; and (4) the relative magnitudes of K sub 1 and K sub 2 associated with the delamination. A simple example of the use of the analysis in predicting delamination crack growth is demonstrated for a glass/epoxy laminate. The comparisons with experimental data show good agreement.

Wang, S. S.↗

Model aerodynamic test results for two variable cycle engine coannular exhaust systems at simulated takeoff and cruise conditions. Comprehensive data report. Volume 3: Graphical data book 1

A graphical presentation of the aerodynamic data acquired during coannular nozzle performance wind tunnel tests is given. The graphical data consist of plots of nozzle gross thrust coefficient, fan nozzle discharge coefficient, and primary nozzle discharge coefficient. Normalized model component static pressure distributions are presented as a function of primary total pressure, fan total pressure, and ambient static pressure for selected operating conditions. In addition, the supersonic cruise configuration data include plots of nozzle efficiency and secondary-to-fan total pressure pumping characteristics. Supersonic and subsonic cruise data are given.

Nelson, D. P.↗

Vacuum UV and soft X-ray optical emissions from electron impact on metals

Spectrograms from 2200 to 20 A of optical emissions produced by a simple lamp using fast (3-10-keV) electron impact on tungsten and tantalum targets are presented. The lamp dissipated up to 5 kW/sq cm on the target surface, simultaneously cleaning the surface and generating bremsstrahlung in the soft X-ray region with 3-keV electron impact and longer wavelength transition radiation at the higher energy. In both the bremsstrahlung and transition radiation modes, the lamp emissions were of sufficient strength and reproducibility to qualify the two mechanisms as potential processes for use in a secondary standard light source. Calculated spectral distributions of transition radiation for normally incident electrons on tungsten, tantalum, and aluminum targets, with a viewing direction 45 deg from the normal, are also shown.

Hughes, R. H.↗

Three-dimensional analysis of 0/90s and 90/0s laminates with a central circular hole

Stress distributions were calculated near a circular hole in laminates, using a three dimensional finite element analysis. These stress distributions were presented three ways: through the thickness at the hole boundary, along radial lines at the 0/90 and 90/0 interfaces, and around the hole at these interfaces. The interlaminar normal stress, and the shear stress, distributions had very steep gradients near the hole boundary, suggesting interlaminar stress singularities. The largest compressive stress occurred at about 60 deg from the load axis. A simple procedure was introduced to calculate interlaminar stresses near the hole boundary. It used stresses calculated by an exact two dimensional analysis of a laminate with a hole as input to a quasi three dimensional model. It produced stresses that agreed closely with those from the three dimensional finite element model.

Raju, I. S.↗

Surface topography-connections between lubrication and failure initiation

The characteristics of the initial surface topography is intimately connected to the machining process by which it is produced. Both processes create a near surface region of residual stresses, microstructure and hardness that is different from the bulk. The material properties in this rather undefined region can significantly influence the mode of failure, such as wear, scuffing or fatigue, as well as the degree of failure resistance for a given material. Under full film elastohydrodynamic (EHD) conditions, where shear is accomodated within a relatively thick lubricant film, the normal and shear stresses are distributed uniformly over the near surface region, and the surface topography has little influence on the lubrication or failure process. Under more typical conditions where surface roughness and lubricant film thickness are of the same order of magnitude, the surface topography not only emerges as an important parameter in failure initiation, but it also becomes intimately involved in the lubrication process itself.

Wedeven, L. D.↗

Wave propagation into the middle atmosphere

Recent observations of various types of waves propagating into the middle atmosphere are reviewed. Emphasis is made on the excitation processes in the lower atmosphere and their vertical propagation through the background flow as a function of the latitude, height and season. The following subjects are discussed: (1) Vertical propagation of quasi-stationary forced Rossby waves into the winter stratosphere in connection with the sudden warming; (2) Spectral distribution and seasonal characteristics of normal mode (free) Rossby waves and the asymmetry of the Northern and Southern Hemispheres; and (3) Seasonal variation of internal gravity waves in the middle atmosphere. Further discussions are presented for future studies based on accumulated observational data during the MAP period.

Hirota, I.↗

Return stroke velocities and currents using a solid state silicon detector system

A small, portable device has been developed to measure return stroke velocities. With the device, velocities from 135 strokes that consist of 92 natural return strokes and 43 triggered return strokes have been analyzed. The average return stroke velocity for longer channels, greater than 500 meters, is 1.2 + or - 0.3 x 10 to the 8th m/s for both natural and triggered return strokes. For shorter channel lengths, less than 500 m, natural lightning has a statistically higher average return stroke velocity of 1.9 + or - 0.7 x 10 to the 8th m/s than triggered lightning with an average return stroke velocity of 1.4 + or - 0.4 x 10 to the 8th m/s. Using the transmission line model of the return stroke, natural lightning has a peak current distribution that is log-normal with a median value of 19 kA. Return stroke velocities and currents were determined for two distant single stroke natural positive cloud-to-ground flashes. The velocities were 1.0 and 1.7 x 10 to the 8th ms/s while the estimated peak current for each positive flash was over 125 kA.

Mach, Douglas M.↗