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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRS1974

Light beams and radiation patterns in random media

The purpose of this paper is to apply a numerical technique in solving problems involving light beams in a random medium. The technique starts by generating numerically sample media with prescribed statistical properties. Rays are then traced in these sample media, and the ray statistics are compiled. These statistics are utilized to find the mean square displacement and distribution of a beam. Problems on beam broadening and distortion of radiation patterns are considered. Fluctuations in wave amplitude and phase are also investigated. When possible, the numerical results are compared with the analytical results and experimental results.

Youakim, M. Y.↗
NASA NTRS1979

Endogenic craters on basaltic lava flows - Size frequency distributions

Circular crater forms, termed collapse depressions, which occur on many basalt flows on the earth have also been detected on the moon and Mars and possibly on Mercury and Io. The admixture of collapse craters with impact craters would affect age determinations of planetary surface units based on impact crater statistics by making them appear anomalously old. In the work described in the present paper, the techniques conventionally used in planetary crater counting were applied to the determination of the size range and size frequency distribution of collapse craters on lava flows in Idaho, California, and New Mexico. Collapse depressions range in size from 3 to 80 m in diameter; their cumulative size distributions are similar to those of small impact craters on the moon.

Greeley, R.↗
NASA NTRS1989

The area of isodensity contours in cosmological models and galaxy surveys

The contour crossing statistic, defined as the mean number of times per unit length that a straight line drawn through the field crosses a given contour, is applied to model density fields and to smoothed samples of galaxies. Models in which the matter is in a bubble structure, in a filamentary net, or in clusters can be distinguished from Gaussian density distributions. The shape of the contour crossing curve in the initially Gaussian fields considered remains Gaussian after gravitational evolution and biasing, as long as the smoothing length is longer than the mass correlation length. With a smoothing length of 5/h Mpc, models containing cosmic strings are indistinguishable from Gaussian distributions. Cosmic explosion models are significantly non-Gaussian, having a bubbly structure. Samples from the CfA survey and the Haynes and Giovanelli (1986) survey are more strongly non-Gaussian at a smoothing length of 6/h Mpc than any of the models examined. At a smoothing length of 12/h Mpc, the Haynes and Giovanelli sample appears Gaussian.

Ryden, Barbara S.↗
NASA NTRS1994

Statistical scaling relationships and size effects in the strength and creep rupture of fibrous composites

In this presentation we discuss a new theoretical model and supporting experimental results for the strength and lifetime in creep rupture of unidirectional, carbon fiber/epoxy matrix composites at ambient conditions. First we review the 'standard' Weibull/power-law methodology that has been standard practice. Then we discuss features of a recent model which build on the statistical aspects of fiber strength, micromechanical aspects of stress transfer around fiber breaks, and time-dependent creep of the matrix. The model is applied to 'microcomposites' consisting of seven fibers in a matrix for which strength and creep-rupture data are available. The model yields Weibull distributions in an envelope format for both strength and lifetime. The respective shape, scale and power-law parameters depend on such parameters as the Weibull shape parameter for fiber strength, the exponent for matrix creep, the effective load transfer length (which grows in time due to matrix creep) and the critical cluster size for failed fibers. The experimental results are consistent with the theory, though time-dependent debonding appears to be part of the failure process.

Phoenix, S. Leigh↗
NASA NTRS1994

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗
NASA NTRS1995

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.↗
NASA NTRS1993

Parameter estimation techniques based on optimizing goodness-of-fit statistics for structural reliability

New methods are presented that utilize the optimization of goodness-of-fit statistics in order to estimate Weibull parameters from failure data. It is assumed that the underlying population is characterized by a three-parameter Weibull distribution. Goodness-of-fit tests are based on the empirical distribution function (EDF). The EDF is a step function, calculated using failure data, and represents an approximation of the cumulative distribution function for the underlying population. Statistics (such as the Kolmogorov-Smirnov statistic and the Anderson-Darling statistic) measure the discrepancy between the EDF and the cumulative distribution function (CDF). These statistics are minimized with respect to the three Weibull parameters. Due to nonlinearities encountered in the minimization process, Powell's numerical optimization procedure is applied to obtain the optimum value of the EDF. Numerical examples show the applicability of these new estimation methods. The results are compared to the estimates obtained with Cooper's nonlinear regression algorithm.

Starlinger, Alois↗
NASA NTRS1995

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗
NASA NTRS1998

Design and Calibration of the X-33 Flush Airdata Sensing (FADS) System

This paper presents the design of the X-33 Flush Airdata Sensing (FADS) system. The X-33 FADS uses a matrix of pressure orifices on the vehicle nose to estimate airdata parameters. The system is designed with dual-redundant measurement hardware, which produces two independent measurement paths. Airdata parameters that correspond to the measurement path with the minimum fit error are selected as the output values. This method enables a single sensor failure to occur with minimal degrading of the system performance. The paper shows the X-33 FADS architecture, derives the estimating algorithms, and demonstrates a mathematical analysis of the FADS system stability. Preliminary aerodynamic calibrations are also presented here. The calibration parameters, the position error coefficient (epsilon), and flow correction terms for the angle of attack (delta alpha), and angle of sideslip (delta beta) are derived from wind tunnel data. Statistical accuracy of' the calibration is evaluated by comparing the wind tunnel reference conditions to the airdata parameters estimated. This comparison is accomplished by applying the calibrated FADS algorithm to the sensed wind tunnel pressures. When the resulting accuracy estimates are compared to accuracy requirements for the X-33 airdata, the FADS system meets these requirements.

Whitmore, Stephen A.↗
NASA NTRS2004

Physical Validation of TRMM TMI and PR Monthly Rain Products Over Oklahoma

The Tropical Rainfall Measuring Mission (TRMM) provides monthly rainfall estimates using data collected by the TRMM satellite. These estimates cover a substantial fraction of the earth's surface. The physical validation of TRMM estimates involves corroborating the accuracy of spaceborne estimates of areal rainfall by inferring errors and biases from ground-based rain estimates. The TRMM error budget consists of two major sources of error: retrieval and sampling. Sampling errors are intrinsic to the process of estimating monthly rainfall and occur because the satellite extrapolates monthly rainfall from a small subset of measurements collected only during satellite overpasses. Retrieval errors, on the other hand, are related to the process of collecting measurements while the satellite is overhead. One of the big challenges confronting the TRMM validation effort is how to best estimate these two main components of the TRMM error budget, which are not easily decoupled. This four-year study computed bulk sampling and retrieval errors for the TRMM microwave imager (TMI) and the precipitation radar (PR) by applying a technique that sub-samples gauge data at TRMM overpass times. Gridded monthly rain estimates are then computed from the monthly bulk statistics of the collected samples, providing a sensor-dependent gauge rain estimate that is assumed to include a TRMM equivalent sampling error. The sub-sampled gauge rain estimates are then used in conjunction with the monthly satellite and gauge (without sub- sampling) estimates to decouple retrieval and sampling errors. The computed mean sampling errors for the TMI and PR were 5.9% and 7.796, respectively, in good agreement with theoretical predictions. The PR year-to-year retrieval biases exceeded corresponding TMI biases, but it was found that these differences were partially due to negative TMI biases during cold months and positive TMI biases during warm months.

Fisher, Brad L.↗
NASA NTRS2005

Identification of a typical flight patterns

Method and system for analyzing aircraft data, including multiple selected flight parameters for a selected phase of a selected flight, and for determining when the selected phase of the selected flight is atypical, when compared with corresponding data for the same phase for other similar flights. A flight signature is computed using continuous- valued and discrete-valued flight parameters for the selected flight parameters and is optionally compared with a statistical distribution of other observed flight signatures, yielding a typicality scores for the same phase for other similar flights. A cluster analysis is optionally applied to the flight signatures to define an optimal collection of clusters. A level of atypicality for a selected flight is estimated, based upon an index associated with the cluster analysis.

Irving C Statler↗
NASA NTRS2016

Earth Observation System Flight Dynamics System Covariance Realism

This presentation applies a covariance realism technique to the National Aeronautics and Space Administration (NASA) Earth Observation System (EOS) Aqua and Aura spacecraft based on inferential statistics. The technique consists of three parts: collection calculation of definitive state estimates through orbit determination, calculation of covariance realism test statistics at each covariance propagation point, and proper assessment of those test statistics.

A-Train↗
NASA NTRS2017

Determination of the Limiting Magnitude

The limiting magnitude of an optical camera system is an important property to understand since it is used to find the completeness limit of observations. Limiting magnitude depends on the hardware and software of the system, current weather conditions, and the angular speed of the objects observed. If an object exhibits a substantial angular rate during the exposure, its light spreads out over more pixels than the stationary stars. This spreading causes the limiting magnitude to be brighter when compared to the stellar limiting magnitude. The effect, which begins to become important when the object moves a full width at half max during a single exposure or video frame. For targets with high angular speeds or camera systems with narrow field of view or long exposures, this correction can be significant, up to several magnitudes. The stars in an image are often used to measure the limiting magnitude since they are stationary, have known brightness, and are present in large numbers, making the determination of the limiting magnitude fairly simple. In order to transform stellar limiting magnitude to object limiting magnitude, a correction must be applied accounting for the angular velocity. This technique is adopted in meteor and other fast-moving object observations, as the lack of a statistically significant sample of targets makes it virtually impossible to determine the limiting magnitude before the weather conditions change. While the weather is the dominant factor in observing satellites, the limiting magnitude for meteors also changes throughout the night due to the motion of a meteor shower or sporadic source radiant across the sky. This paper presents methods for determining the limiting stellar magnitude and the conversion to the target limiting magnitude.

Kingery, Aaron↗
NASA NTRS2018

Follow-up Imaging of Disk Candidates from the Disk Detective Citizen Science Project: New Discoveries and False Positives in WISE Circumstellar Disk Surveys

The Disk Detective citizen science project aims to find new stars with excess 22 m emission from circumstellar dust in the All WISE data release from the Wide-field Infrared Survey Explorer. We evaluated 261 Disk Detective objects of interest with imaging with the Robo-AO adaptive optics instrument on the 1.5 m telescope at Palomar Observatory and with RetroCam on the 2.5 m du Pont Telescope at Las Campanas Observatory to search for background objects at 0 15–12 separations from each target. Our analysis of these data leads us to reject 7% of targets. Combining this result with statistics from our online image classification efforts implies that at most7.9%±0.2% of All WISE-selected infrared excesses are good disk candidates. Applying our false-positive rates to other surveys, we find that the infrared excess searches of McDonald et al. and Marton et al. all have false-positiverates >70%. Moreover, we find that all 13 disk candidates in Theissen & West with W4 signal-to-noise ratio >3are false positives. We present 244 disk candidates that have survived vetting by follow-up imaging. Of these,213 are newly identified disk systems. Twelve of these are candidate members of comoving pairs based on Gaia astrometry, supporting the hypothesis that warm dust is associated with binary systems. We also note the discovery of 22 m excess around two known members of the Scorpius–Centaurus association, and we identifyknown disk host WISEA J164540.79-310226.6 as a likely Sco-Cen member. Thirty of these disk candidates arecloser than 125 pc (including 26 debris disks), making them good targets for both direct-imaging exoplanetsearches.

Wide-field Infrared Survey Explorer↗
NASA NTRS2019

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

We present new Chandra constraints on the X-ray luminosity functions(XLFs)of X-ray binary(XRB)populations, as well as their scaling relations, for a sample of 38 nearby galaxies(D = 3.4–29 Mpc). Our galaxy sample is drawn primarily from the Spitzer Infrared Nearby Galaxies Survey(SINGS)and contains a wealth of Chandra(5.8 Ms total)and multiwavelength data, allowing for star formation rates(SFRs)and stellar masses(M⁎)to be measured on subgalactic scales. We divided the 2478 X-ray-detected sources into 21 subsamples in bins of specific SFR(sSFR ≡ SFR/M⁎)and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs(HMXBs), low-mass XRBs(LMXBs), and background sources from the cosmic X-ray background that respectively scale with SFR, M⁎, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above≈10^(38)erg/s, and elliptical galaxies that have relatively rich populations of globular clusters(GCs)show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies.

Bret D. Lehmer↗
NASA NTRS2021

Truncated ARQ Statistical Link Analysis for Dynamic Links

The future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non Gaussian and non linear effects such as atmospheric turbulence, scintillation, antenna mis-pointing, jitter, etc. These dynamic links thus will experience various degrees of fading loss, and some of these link disruptions cannot be effectively mitigated by forward error correction coding and/or interleaving. One effective way to ensure reliable communication is by using Automatic Repeat Request (ARQ) protocol, where the receiver acknowledges to the transmitter whether or not a data unit is successfully received. If a data unit is not successfully received (such as after a pre-set time-out), the transmitter would then re-transmit the lost data unit to the receiver. In a previous paper, we derived a statistical link analysis method of finding the optimal operating Signal-to-Noise Ratio (SNR) and estimating the latency of an ARQ scheme. In a more recent paper, we demonstrated the above method using the SNR distribution constructed from the Ka-band (32 GHz) flight data. To simplify the discussion, we considered the academic approach that the ARQ scheme allows for an infinite number of retransmissions. In this paper, we consider the more practical case of a truncated ARQ scheme, where there is a limit on the number of retransmissions. We derive the error probability, the optimal SNR setting, and the latency statistics of the correctly received frames of the truncated ARQ schemes. We first discuss the truncated ARQ link analysis principles using the Gaussian assumption for SNR distribution with a large variance. Next, we demonstrate the statistical truncated ARQ link analysis using the SNR distribution constructed from the Ka-band flight data. The results in this paper can be applied in the design of reliable communication systems such as the Consultative Committee for Space Data System (CCSDS) File Transfer Protocol (CFTP) and the Delay Tolerant Network (DTN).

Morabito, David↗
NASA NTRS2022

Quantifying urban, industrial, and background changes in NO2 during the COVID-19 lockdown period based on TROPOMI satellite observations

The COVID-19 lockdown had a large impact on anthropogenic emissions of air pollutants and particularly on nitrogen dioxide (NO 2 ). While the overall NO 2 decline over some large cities is well-established, understanding the details remains a challenge since multiple source categories contribute. In this study, a new method of isolation of three components (background NO 2 , NO 2 from urban sources, and NO 2 from industrial point sources) is applied to estimate the impact of the COVID-19 lockdown on each of them. The approach is based on fitting satellite data by a statistical model with empirical plume dispersion functions driven by a meteorological reanalysis. Population density and surface elevation data as well as coordinates of industrial sources were used in the analysis. The tropospheric NO 2 vertical column density (VCD) values measured by the Tropospheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor over 261 urban areas for the period from 16 March to 15 June 2020 were compared with the average VCD values for the same period in 2018 and 2019. While the background NO 2 component remained almost unchanged, the urban NO 2 component declined by −18 % to −28 % over most regions. India, South America, and a part of Europe (particularly, Italy, France, and Spain) demonstrated a −40 % to −50 % urban emission decline. In contrast, the decline over urban areas in China, where the lockdown was over during the analysed period, was, on average, only -4.4± 8 %. Emissions from large industrial sources in the analysed urban areas varied greatly from region to region from -4.8 ± 6 % for China to -40 ± 10 % for India. Estimated changes in urban emissions are correlated with changes in Google mobility data (the correlation coefficient is 0.62) confirming that changes in traffic were one of the key elements in the decline in urban NO2 emissions. No correlation was found between changes in background NO2 and Google mobility data. On the global scale, the background and urban components were remarkably stable in 2018, 2019, and 2021, with averages of all analysed areas all being within ±2.5 % and suggesting that there were no substantial drifts or shifts in TROPOMI data. The 2020 data are clearly an outlier: in 2020, the mean background component for all analysed areas (without China) was -6.0 % ± 1.2 % and the mean urban component was -26.7 ± 2.6 % or 20σ below the baseline level from the other years.

COVID-19↗