Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 685 records · Page 38

Research in solar physics: Some techniques for analyzing data from the ultraviolet spectrometer and polarimeter

Useful information for certain aspects of the analysis of data obtained by the Ultraviolet Spectrometer and Polarimeter (UVSP) on the Solar Maximum Mission (SMM) are contained. The meaning of the UVSPCO-ordinate system and the SMM roll, pitch, and yaw are described and the process for overlaying UVSP images is explained. The various computer programs that calculate the line of sight component of the SMM spacecraft velocity from the spacecraft to the Sun is described. The spacecraft velocity is used to correct or interpret the signal observed in USVP dopplergrams. A method of using the spacecraft velocity to calibrate UVSP dopplergrams and magnetograms, i.e., determine the width of the observed emission line is applied. The UVSP polarization analysis procedures are described and the expressions for the statistical uncertainties in various quantities obtained from UVSP measurements are given.

Henze, W., Jr.↗

Ring-diameter Ratios for Multi-ring Basins Average 2.0(0.5)D

The spacing of the concentric rings of planetary impact basins was studied. It is shown that a radial increment of x (sup 0.5) D, where x is about 2.0 and D = ring diameter, separates both (1) adjacent least-squares groups of rings and arcs of multi-ring basins on Mars, Mercury, and the Moon; and (2) adjacent rings of individual basins on the three planets. Statistics for ratios of ring diameters are presented, the first and most-applied parameter of ring spacing. It is found that ratios excluding rings flanking the main ring also have a mean spacing increment of about 2.0. Ratios including such rings, as for the least-squares groups, and (1) above, have a larger increment, averaging 2.1. The F-test indicates, that these spacings of basin ring locations, and mode of ring formation are controlled by the mechanics of the impact event itself, rather than by crustal properties.

Pike, R. J.↗

The statistical equilibrium of hydrogen and helium in a radiation field, with an application to interpreting Wolf-Rayet spectra

The interpretation of the intensities of the hydrogen and helium emission lines in O and Wolf-Rayet spectra in terms of the abundance of hydrogen relative to helium requires information regarding the distribution of hydrogen and helium atoms and ions over their several energy states. In addition, some estimate is needed regarding the transmission of the radiation through the stellar mantle. The present paper provides new information concerning the population of the energy levels of hydrogen and helium when statistical equilibrium occurs in the presence of a radiation field. The results are applied to an interpretation of the spectra of four Wolf-Rayet stars, taking into account the implications for interpreting the spectra of O stars, OB supergiants, and Be stars.

Bhatia, A. K.↗

Continuous fiber ceramic matrix composites for heat engine components

High strength at elevated temperatures, low density, resistance to wear, and abundance of nonstrategic raw materials make structural ceramics attractive for advanced heat engine applications. Unfortunately, ceramics have a low fracture toughness and fail catastrophically because of overload, impact, and contact stresses. Ceramic matrix composites provide the means to achieve improved fracture toughness while retaining desirable characteristics, such as high strength and low density. Materials scientists and engineers are trying to develop the ideal fibers and matrices to achieve the optimum ceramic matrix composite properties. A need exists for the development of failure models for the design of ceramic matrix composite heat engine components. Phenomenological failure models are currently the most frequently used in industry, but they are deterministic and do not adequately describe ceramic matrix composite behavior. Semi-empirical models were proposed, which relate the failure of notched composite laminates to the stress a characteristic distance away from the notch. Shear lag models describe composite failure modes at the micromechanics level. The enhanced matrix cracking stress occurs at the same applied stress level predicted by the two models of steady state cracking. Finally, statistical models take into consideration the distribution in composite failure strength. The intent is to develop these models into computer algorithms for the failure analysis of ceramic matrix composites under monotonically increasing loads. The algorithms will be included in a postprocessor to general purpose finite element programs.

Tripp, David E.↗

Determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution, appendix 4

The power spectrum for a stationary random process can be defined with the Wiener-Khintchine Theorem, which says that the power spectrum and the auto correlation function are a Fourier transform pair. To implement this theorem for signals that are discrete and of finite length we can use the Blackman-Tukey method. Blackman and Tukey (1958) show that a function w(tau), called a lag window, can be applied to the auto correlation estimates to obtain power spectrum estimates that are statistically stable. The Fourier transform of w(r) is called a spectral window. Typical choices for spectral windows show a distinct trade-off between the main lobe width and side lobe strength. A new idea for designing windows by taking linear combinations of the standard windows to produce hybrid windows was introduced by Smith (1985). We implement Smith's idea to obtain spectral windows with narrow main lobes and smaller (compared with typical windows) near side lobes. One of the main contributions of this thesis is that we show that Smith's problem is equivalent to a Quadratic Programming (QP) problem with linear equality and inequality constraints. A computer program was written to produce hybrid windows by setting up and solving the QP problem. We also developed and solved two variations of the original problem. The two variations involved changing the inequality constraints in both cases from non negativity on the combination coefficients to non negativity on the hybrid lag window itself. For the second variation, the window functions used to construct the hybrid window were changed to a frequency-variable set of truncated cosinusoids. A series of tests was run with the three computer programs to investigate the behavior of the hybrid spectral and lag windows. Emphasis was put on obtaining spectral windows with both relatively narrow main lobes and the lowest possible (for these algorithms) near side lobes. Some success was achieved for this goal. A 10 dB peak side lobe reduction over the rectangular spectral window without significant main lobe broadening was achieved. Also, average side lobe levels of -117 dB were reached at a cost of doubling the main lobe width (at the -3 dB point).

Ioup, George E.↗

Simulation of rotor blade element turbulence

A piloted, motion-based simulation of Sikorsky's Black Hawk helicopter was used as a platform for the investigation of rotorcraft responses to vertical turbulence. By using an innovative temporal and geometrical distribution algorithm that preserved the statistical characteristics of the turbulence over the rotor disc, stochastic velocity components were applied at each of twenty blade-element stations. This model was implemented on NASA Ames' Vertical Motion Simulator (VMS), and ten test pilots were used to establish that the model created realistic cues. The objectives of this research included the establishment of a simulation-technology basis for future investigation into real-time turbulence modeling. This goal was achieved; our extensive additions to the rotor model added less than a 10 percent computational overhead. Using a VAX 9000 computer the entire simulation required a cycle time of less than 12 msec. Pilot opinion during this simulation was generally quite favorable. For low speed flight the consensus was that SORBET (acronym for title) was better than the conventional body-fixed model, which was used for comparison purposes, and was determined to be too violent (like a washboard). For high speed flight the pilots could not identify differences between these models. These opinions were something of a surprise because only the vertical turbulence component on the rotor system was implemented in SORBET. Because of the finite-element distribution of the inputs, induced outputs were observed in all translational and rotational axes. Extensive post-simulation spectral analyses of the SORBET model suggest that proper rotorcraft turbulence modeling requires that vertical atmospheric disturbances not be superimposed at the vehicle center of gravity but, rather, be input into the rotor system, where the rotor-to-body transfer function severely attenuates high frequency rotorcraft responses.

Mcfarland, R. E.↗

Characterization of Sound Radiation by Unresolved Scales of Motion in Computational Aeroacoustics

Evaluation of the sound sources in a high Reynolds number turbulent flow requires time-accurate resolution of an extremely large number of scales of motion. Direct numerical simulations will therefore remain infeasible for the forseeable future: although current large eddy simulation methods can resolve the largest scales of motion accurately the, they must leave some scales of motion unresolved. A priori studies show that acoustic power can be underestimated significantly if the contribution of these unresolved scales is simply neglected. In this paper, the problem of evaluating the sound radiation properties of the unresolved, subgrid-scale motions is approached in the spirit of the simplest subgrid stress models: the unresolved velocity field is treated as isotropic turbulence with statistical descriptors, evaluated from the resolved field. The theory of isotropic turbulence is applied to derive formulas for the total power and the power spectral density of the sound radiated by a filtered velocity field. These quantities are compared with the corresponding quantities for the unfiltered field for a range of filter widths and Reynolds numbers.

Rubinstein, Robert↗

User's Guide for Monthly Vector Wind Profile Model

The background, theoretical concepts, and methodology for construction of vector wind profiles based on a statistical model are presented. The derived monthly vector wind profiles are to be applied by the launch vehicle design community for establishing realistic estimates of critical vehicle design parameter dispersions related to wind profile dispersions. During initial studies a number of months are used to establish the model profiles that produce the largest monthly dispersions of ascent vehicle aerodynamic load indicators. The largest monthly dispersions for wind, which occur during the winter high-wind months, are used for establishing the design reference dispersions for the aerodynamic load indicators. This document includes a description of the computational process for the vector wind model including specification of input data, parameter settings, and output data formats. Sample output data listings are provided to aid the user in the verification of test output.

Adelfang, S. I.↗

A TRMM-Calibrated Infrared Rainfall Algorithm Over Brazil

An improved version of the Convective/Stratiform Technique (CST) and its application over Brazil are presented. Keeping the major components of traditional CST, we have modified the scheme for classing convective and stratiform rainfall and tested various schemes to eliminate non-raining cirrus clouds. The parameters of the technique are calibrated using Tropical Rainfall Measuring Mission (TRMM) multi-sensor observations, including TMI, PR and Visible Infrared Scanner (VIRS). This technique takes advantage of the high temporal sampling rate of geosynchronous infrared channel and the better instantaneous rain observation of TRMM Microwave Imager (TMI) and PR. Moreover, sparse TMI rain estimates are used to dynamically calibrate the parameters of the technique to improve its performance. The technique is applied to make rainfall estimates on various time scales and to study rainfall statistics such as the distributions of rain intensity and storm duration over a four-month period beginning January 1999. The study period coincides with the TRMM/Large-Scale Biosphere-Atmosphere experiment in Amazonia (LBA) ground validation experiment, and observations from the LBA radar are used to validate the technique. Results show the mean diurnal cycle of precipitation over the LBA area, and are compared to radar data from the Tropical Oceans and Global Atmosphere (TOGA) radar. When the larger geographic region is considered, the analysis of the mean hourly estimates revealed the pronounced effects of rivers, topography and local circulations on the rainfall.

Negri, Andrew J.↗

Preliminary Evaluation of an Aviation Safety Thesaurus' Utility for Enhancing Automated Processing of Incident Reports

This document presents a preliminary evaluation the utility of the FAA Safety Analytics Thesaurus (SAT) utility in enhancing automated document processing applications under development at NASA Ames Research Center (ARC). Current development efforts at ARC are described, including overviews of the statistical machine learning techniques that have been investigated. An analysis of opportunities for applying thesaurus knowledge to improving algorithm performance is then presented.

Barrientos, Francesca↗

Vapor Hydrogen Peroxide as Alternative to Dry Heat Microbial Reduction

The Jet Propulsion Laboratory, in conjunction with the NASA Planetary Protection Officer, has selected vapor phase hydrogen peroxide (VHP) sterilization process for continued development as a NASA approved sterilization technique for spacecraft subsystems and systems. The goal is to include this technique, with appropriate specification, in NPG8020.12C as a low temperature complementary technique to the dry heat sterilization process. A series of experiments were conducted in vacuum to determine VHP process parameters that provided significant reductions in spore viability while allowing survival of sufficient spores for statistically significant enumeration. With this knowledge of D values, sensible margins can be applied in a planetary protection specification. The outcome of this study provided an optimization of test sterilizer process conditions: VHP concentration, process duration, a process temperature range for which the worst case D value may be imposed, a process humidity range for which the worst case D value may be imposed, and robustness to selected spacecraft material substrates.

vapor hydrogen peroxide↗

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, N-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-N-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations k 2 = 2k 1 = 0.2 h Mpc -1 by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the baryon acoustic oscillations (BAO) scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity (NG).

79 ASTRONOMY AND ASTROPHYSICS↗

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantitative analysis of Raman spectral parameters for carbon fibers: practical considerations and connection to mechanical properties

Although the literature on the Raman spectra of carbon fibers is vast, no consistent, robust predictive relationship between mechanical properties of carbon fibers and spectral parameters exists. This shortcoming is due to the use of numerous fitting functions to evaluate Raman spectra of carbon fibers and the inconsistencies in establishing the best fitting models in a statistically robust fashion. To address this gap, we present a comprehensive work on the Raman spectra of carbon fibers that combines a vast library of experimental data with a robust numerical analysis and a statistical evaluation of a wide range of suggested fitting models. This manuscript begins with a brief review of the commonly applied fitting models. Then, the Raman spectra of 32 commercially available polyacrylonitrile-based carbon fibers collected at excitation wavelengths 532, 633, and 785 nm are presented and the best fit for all fibers is evaluated based on several statistical criteria in conjunction with numerical calculations and physical arguments. The results suggest that high-performance fibers must be fit with at least five peaks, whereas high-tensile modulus fibers are best fit with at least six distinct peaks. Finally, we employ simultaneous fitting of the Raman spectra of specific fibers and wavelengths and demonstrate that strong correlations exist between mechanical properties and the D1 peak position and shape across the range of evaluated mechanical properties. Finally, we suggest straightforward improvements in fitting analysis procedures that can be implemented to increase coherency in the understanding of the underlying carbon fiber microstructure intuited from Raman spectroscopy.

36 MATERIALS SCIENCE↗

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

Emergence of robust anthropogenic increase of heat stress-related variables projected from CORDEX-CORE climate simulations

The information of when and where region-specific patterns in both mean and extreme temperatures leading to heat stress will emerge from the present-day climate variability is important to plan adaptation options, but to date studies on this issue still remain limited and fragmented. Here, we estimate the time of emergence (TOE) of temperature and wet-bulb temperature (Tw), a better indication of heat stress, using fine-scale, long-term regional climate model projections under the RCP2.6 and RCP8.5 scenarios across six different domains. Differently from previous studies, the TOE is determined using three methods applied on impact-relevant variables: two different signal-to-noise frameworks based on summer mean temperature and Tw and a statistical test to identify significant differences in daily extreme distributions. The TOE response to RCP2.6 and RCP8.5 with respect to the end of 20th century variability differs significantly regardless of which TOE metric is applied. For summer mean temperature, the land fraction reaching TOE is expected to exceed 90% by the 2050s under the RCP8.5, whereas the increase rate of land exposure to TOE tends to stagnate over time under the RCP2.6 so that more than 40% of land will not experience TOE by the end of the 21st century. Compared to temperature, the TOE of Tw is reached earlier in most of the wet tropics but is delayed in hot and dry regions because of the nonlinear response of Tw to humidity. For both temperature and Tw, the TOE appears earlier in regions with low baseline variability, such as in the tropics. Despite the uncertainties arising from the choice of TOE metrics, the vast majority of regions in Africa and southeast Asia experience TOE in the early 21st century under both the RCP2.6 and RCP8.5 scenarios, which stresses the urgent need for developing adequate adaptation strategies in these regions.

54 ENVIRONMENTAL SCIENCES↗