Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution parameters”

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 595 records · Page 33

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: A New Type of Mass-Flux Model - Part 1

Higher-order closure (HOC) models have been proposed for parameterization of the turbulent planetary boundary layer (PBL). HOC models must include closures for higher-order moments (e.g., fourth moments in third-order closure models), for pressure terms, and for dissipation terms. Mass-flux closure (MFC) models have been proposed for parameterization of cumulus convection and, more recently, the convective PBL. MFC models include closures for lateral mass exchanges and for pressure terms (which are usually ignored). The authors developed a new kind of model that combines HOC and MFC, which they hope will be useful for the parameterization of both the PBL and cumulus convection, in a unified framework. Such a model is particularly well suited to regimes in which the PBL turbulence and the cumulus convection are not well separated, for example, the broken stratocumulus and shallow cumulus regimes. The model makes use of an assumed joint probability distribution for the variables of interest, and the equations typically used in HOC models can be derived by integrating over the distribution. Accordingly, the model is called Assumed-Distribution Higher-Order Closure (ADHOC). The prognostic variables of ADHOC are the mean state, the second and third moments of the vertical velocity, and the vertical fluxes of other quantities of interest. All of the parameters of the distribution can be determined from the predicted moments, thereafter the joint distribution is effectively known, and so any and all moments can be constructed as needed. In this way, the usual closure problem of "higher moments" is avoided. The pressure-term parameterizations previously developed for HOC models are used to predict the convective fluxes and the moments of the vertical velocity. In companion papers, parameterizations of lateral mass exchanges and subplume-scale fluxes are presented, and then ADHOC is applied to several observationally based tropical, subtropical, and dry convective boundary layers.

Lappen, Cara-Lyn↗

Constraining the Neutron Star Mass–Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes

We describe the X-ray pulse profile models we use and how we use them to analyze Neutron Star Interior Composition Explorer(NICER)observations of rotation-powered millisecond pulsars to obtain information about the mass–radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. and Riley et al. and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. and Riley et al. by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass–radius plane that are statistically consistent with one another, and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar’s spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in the estimates made using this technique, some of which may warrant further investigation.

Slavko Bogdanov↗

Lifetime Reliability Evaluation of Structural Ceramic Parts with the CARES/LIFE Computer Program

The computer program CARES/LIFE calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. This program is an extension of the CARES (Ceramics Analysis and Reliability Evaluation of Structures) computer program. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker equation. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled using either the principle of independent action (PIA), Weibull's normal stress averaging method (NSA), or Batdorf's theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. Two example problems demonstrating cyclic fatigue parameter estimation and component reliability analysis with proof testing are included.

Nemeth, Noel N.↗

Time-dependent reliability analysis of ceramic engine components

The computer program CARES/LIFE calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. This program is an extension of the CARES (Ceramics Analysis and Reliability Evaluation of Structures) computer program. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing either the power or Paris law relations. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. Two example problems demonstrating proof testing and fatigue parameter estimation are given.

Nemeth, Noel N.↗

Durability evaluation of ceramic components using CARES/LIFE

The computer program CARES/LIFE calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. This program is an extension of the CARES (Ceramics Analysis and Reliability Evaluation of Structures) computer program. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker equation. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. Application of this design methodology is demonstrated using experimental data from alumina bar and disk flexure specimens which exhibit SCG when exposed to water.

Nemeth, Noel N.↗

Optimum design of structures subject to general periodic loads

A simplified version of Icerman's problem regarding the design of structures subject to a single harmonic load is discussed. The nature of the restrictive conditions that must be placed on the design space in order to ensure an analytic optimum are discussed in detail. Icerman's problem is then extended to include multiple forcing functions with different driving frequencies. And the conditions that now must be placed upon the design space to ensure an analytic optimum are again discussed. An important finding is that all solutions to the optimality condition (analytic stationary design) are local optima, but the global optimum may well be non-analytic. The more general problem of distributing the fixed mass of a linear elastic structure subject to general periodic loads in order to minimize some measure of the steady state deflection is also considered. This response is explicitly expressed in terms of Green's functional and the abstract operators defining the structure. The optimality criterion is derived by differentiating the response with respect to the design parameters. The theory is applicable to finite element as well as distributed parameter models.

Reiss, Robert↗

Effects of anisotropic electron-ion interactions in atomic photoelectron angular distributions

The photoelectron asymmetry parameter beta in LS-coupling is obtained as an expansion into contributions from alternative angular momentum transfers j sub t. The physical significance of this expansion of beta is shown to be that: (1) the electric dipole interaction transfers to the atom a charcteristic single angular momentum j sub t = sub o, where sub o is the photoelectron's initial orbital momentum; and (2) angular momentum transfers indicate the presence of anisotropic interaction of the outgoing photoelectron with the residual ion. For open shell atoms the photoelectron-ion interaction is generally anisotropic; photoelectron phase shifts and electric dipole matrix elements depend on both the multiplet term of the residual ion and the total orbital momentum of the ion-photoelectron final state channel. Consequently beta depends on the term levels of the residual ion and contains contributions from all allowed values of j sub t. Numerical calculations of the asymmetry parameters and partial cross sections for photoionization of atomic sulfur are presented.

Dill, D.↗

Parameter Estimation in Astronomy with Poisson-Distributed Data: The (CHI)2(gamma) Statistic - 1

Applying the standard weighted mean formula, [Sigma (sub i)n(sub i)ssigma(sub i, sup -2)], to determine the weighted mean of data, n(sub i), drawn from a Poisson distribution, will, on average, underestimate the true mean by approx. 1 for all true mean values larger than approx.3 when the common assumption is made that the error of the i th observation is sigma(sub i) = max square root of n(sub i), 1).This small, but statistically significant offset, explains the long-known observation that chi-square minimization techniques which use the modified Neyman'chi(sub 2) statistic, chi(sup 2, sub N) equivalent Sigma(sub i)((n(sub i) - y(sub i)(exp 2)) / max(n(sub i), 1), to compare Poisson - distributed data with model values, y(sub i), will typically predict a total number of counts that underestimates the true total by about 1 count per bin. Based on my finding that weighted mean of data drawn from a Poisson distribution can be determined using the formula [Sigma(sub i)[n(sub i) + min(n(sub i), 1)](n(sub i) + 1)(exp -1)] / [Sigma(sub i)(n(sub i) + 1)(exp -1))], I propose that a new chi(sub 2) statistic, chi(sup 2, sub gamma) equivalent, should always be used to analyze Poisson- distributed data in preference to the modified Neyman's chi(exp 2) statistic. I demonstrated the power and usefulness of,chi(sub gamma, sup 2) minimization by using two statistical fitting techniques and five chi(exp 2) statistics to analyze simulated X-ray power - low 15 - channel spectra with large and small counts per bin. I show that chi(sub gamma, sup 2) minimization with the Levenberg - Marquardt or Powell's method can produce excellent results (mean slope errors approx. less than 3%) with spectra having as few as 25 total counts.

Mighell, Kenneth J.↗

CARES/LIFE Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program

This manual describes the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction (CARES/LIFE) computer program. The program calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. CARES/LIFE is an extension of the CARES (Ceramic Analysis and Reliability Evaluation of Structures) computer program. The program uses results from MSC/NASTRAN, ABAQUS, and ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker law. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled by using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. The probabilistic time-dependent theories used in CARES/LIFE, along with the input and output for CARES/LIFE, are described. Example problems to demonstrate various features of the program are also included.

Nemeth, Noel N.↗

The CARMENES Search for Exoplanets Around M Dwarfs: Revisiting the GJ 581 Multi-Planetary System With New Doppler Measurements From CARMENES, HARPS, and HIRES

Context.GJ 581 is a nearby M dwarf known to host a packed multiple planet system composed of two super-Earths and a Neptune-mass planet. We present new orbital analyses of the GJ 581 system, utilizing recent radial velocity (RV) data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Aims.Our aim was to analyze the post-discovery spectroscopic data of GJ 581, which were obtained with CARMENES. In addition, we used publicly available HIRES and HARPS spectroscopic data to seek evidence of the known and disputed exoplanets in this system. We aimed to investigate the stellar activity of GJ 581 and update the planetary system’s orbital parameters using state-of-the-art numerical models and techniques. Methods.We performed a periodogram analysis of the available precise CARMENES, HIRES, and HARPS RVs and of stellar activity indicators. We conducted detailed orbital analyses by testing various orbital configurations consistent with the RV data. We studied the posterior probability distribution of the parameters fit to the data and we explored the long-term stability and overall orbital dynamics of the GJ 581 system. Results. We refined the orbital parameters of the GJ 581 system using the most precise and complete set of Doppler data available. Consistent with the existing literature, our analysis confirms that the system is unequivocally composed of only three planets detectable in the present data, dismissing the putative planet GJ 581 d as an artifact of stellar activity. Our N-body fit reveals that the system’s inclination is i=47.0+14.6−13.0deg, which implies that the planets could be up to 30% more massive than their previously reported minimum masses. Furthermore, we report that the GJ 581 system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances.

Satellites↗

IMF orientation, solar wind velocity, and Pc 3-4 signals - A joint distribution

Separate studies using the same micropulsation data base in the period range 10-150 s have shown earlier that signal levels recorded during September, October, and November 1969 at Calgary correlated positively with both solar-wind alignment of the IMF and solar-wind speed, but each correlation contained enough scatter to allow for the influence of the other factor. In this report, joint correlations of velocity and field direction with parameters representing hourly distributions rather than minima of IMF orientation angle display the relative effect of the two agents on magnetic pulsation signal levels. The joint correlations reduce the overall scatter and show that solar-wind speeds above 200-300 km/s and angles between the IMF and the sun-earth line of less than 50-60 deg are associated with enlarged magnetic pulsation amplitudes. These threshold effects tend to support both the bow-shock origin and the Kelvin-Helmholtz amplification of daytime signal transients in the Pc 3, 4 period ranges.

Greenstadt, E. W.↗

Bulk properties and velocity distributions of water group ions at Comet Halley - Giotto measurements

In the region upstream of Comet Halley, pickup heavy ions of cometary origin were directly observed by the implanted ion spectrometer on Giotto. Diffusion of this population in pitch angle and in energy, during the approach to the comet and on the outbound leg is discussed. The two data sets are compared and qualitative ideas on scattering timescales are inferred. In addition the bulk parameters of these distributions have been computed and a comparison of the observed speed in the solar wind frame and the observed density with expectations is presented. Pitch angle scattering occurs more slowly than expected with filled shells appearing at 2,500,000 km, and significant energy diffusion does not occur until the bow shock region. Also the shell distributions downstream of the shock flow at the bispherical bulk speed (related to the Alfven speed) along the magnetic field with respect to the solar wind in accordance with conservation of energy between the pickup ions and the wave turbulence.

Coates, A. J.↗

Linear Estimation of Particle Bulk Parameters from Multi-Wavelength Lidar Measurements

An algorithm for linear estimation of aerosol bulk properties such as particle volume, effective radius and complex refractive index from multiwavelength lidar measurements is presented. The approach uses the fact that the total aerosol concentration can well be approximated as a linear combination of aerosol characteristics measured by multiwavelength lidar. Therefore, the aerosol concentration can be estimated from lidar measurements without the need to derive the size distribution, which entails more sophisticated procedures. The definition of the coefficients required for the linear estimates is based on an expansion of the particle size distribution in terms of the measurement kernels. Once the coefficients are established, the approach permits fast retrieval of aerosol bulk properties when compared with the full regularization technique. In addition, the straightforward estimation of bulk properties stabilizes the inversion making it more resistant to noise in the optical data. Numerical tests demonstrate that for data sets containing three aerosol backscattering and two extinction coefficients (so called 3 + 2 ) the uncertainties in the retrieval of particle volume and surface area are below 45% when input data random uncertainties are below 20 %. Moreover, using linear estimates allows reliable retrievals even when the number of input data is reduced. To evaluate the approach, the results obtained using this technique are compared with those based on the previously developed full inversion scheme that relies on the regularization procedure. Both techniques were applied to the data measured by multiwavelength lidar at NASA/GSFC. The results obtained with both methods using the same observations are in good agreement. At the same time, the high speed of the retrieval using linear estimates makes the method preferable for generating aerosol information from extended lidar observations. To demonstrate the efficiency of the method, an extended time series of observations acquired in Turkey in May 2010 was processed using the linear estimates technique permitting, for what we believe to be the first time, temporal-height distributions of particle parameters.

linear estimation↗

Seismic properties of the crust and uppermost mantle of North America

Seismic refraction profiles for the North American continent were compiled. The crustal models compiled data on the upper mantle seismic velocity (P sub n), the crustal thickness (H sub c) and the average seismic velocity of the crystalline crust (V sub p). Compressional wave parameters were compared with shear wave data derived from surface wave dispersion models and indicate an average value for Poisson's ratio of 0.252 for the crust and of 0.273 for the uppermost mantle. Contour maps illustrate lateral variations in crustal thickness, upper mantle velocity and average seismic velocity of the crystalline crust. The distribution of seismic parameters are compared with a smoothed free air anomaly map of North America and indicate that a complidated mechanism of isostatic compensation exists for the North American continent. Several features on the seismic contour maps also correlate with regional magnetic anomalies.

Braile, L. W.↗

Improvement of Stratospheric Aerosol Extinction Retrieval from OMPS/LP Using a New Aerosol Model

The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) has been flying on the Suomi National Polar-orbiting Partnership (S-NPP) satellite since October 2011. It is designed to produce ozone and aerosol vertical profiles at 2-kilometer vertical resolution over the entire sunlit globe. Aerosol extinction profiles are computed with Mie theory using radiances measured at 675 nanometers. The operational Version 1.0 (V1.0) aerosol extinction retrieval algorithm assumes a bimodal lognormal aerosol size distribution (ASD) whose parameters were derived by combining an in situ measurement of aerosol microphysics with the Stratospheric Aerosol and Gas Experiment (SAGE II) aerosol extinction climatology. Internal analysis indicates that this bimodal lognormal ASD does not sufficiently explain the spectral dependence of LP-measured radiances. In this paper we describe the derivation of an improved aerosol size distribution, designated Version 1.5 (V1.5), for the LP retrieval algorithm. The new ASD uses a gamma function distribution that is derived from Community Aerosol and Radiation Model for Atmospheres (CARMA)-calculated results. A cumulative distribution fit derived from the gamma function ASD gives better agreement with CARMA results at small particle radii than bimodal or unimodal functions. The new ASD also explains the spectral dependence of LP-measured radiances better than the V1.0 ASD. We find that the impact of our choice of ASD on the retrieved extinctions varies strongly with the underlying reflectivity of the scene. Initial comparisons with collocated extinction profiles retrieved at 676 nanometers from the SAGE III instrument on the International Space Station (ISS) show a significant improvement in agreement for the LP V1.5 retrievals. Zonal mean extinction profiles agree to within 10 percent between 19 and 29 kilometers, and regression fits of collocated samples show improved correlation and reduced scatter compared to the V1.0 product. This improved agreement will motivate development of more sophisticated ASDs from CARMA results that incorporate latitude, altitude and seasonal variations in aerosol properties.

Ozone Mapping Profiler Suite Limb Profiler (OMPS/L↗

Technique for Solving Electrically Small to Large Structures for Broadband Applications

Fast iterative algorithms are often used for solving Method of Moments (MoM) systems, having a large number of unknowns, to determine current distribution and other parameters. The most commonly used fast methods include the fast multipole method (FMM), the precorrected fast Fourier transform (PFFT), and low-rank QR compression methods. These methods reduce the O(N) memory and time requirements to O(N log N) by compressing the dense MoM system so as to exploit the physics of Green s Function interactions. FFT-based techniques for solving such problems are efficient for spacefilling and uniform structures, but their performance substantially degrades for non-uniformly distributed structures due to the inherent need to employ a uniform global grid. FMM or QR techniques are better suited than FFT techniques; however, neither the FMM nor the QR technique can be used at all frequencies. This method has been developed to efficiently solve for a desired parameter of a system or device that can include both electrically large FMM elements, and electrically small QR elements. The system or device is set up as an oct-tree structure that can include regions of both the FMM type and the QR type. The system is enclosed with a cube at a 0- th level, splitting the cube at the 0-th level into eight child cubes. This forms cubes at a 1st level, recursively repeating the splitting process for cubes at successive levels until a desired number of levels is created. For each cube that is thus formed, neighbor lists and interaction lists are maintained. An iterative solver is then used to determine a first matrix vector product for any electrically large elements as well as a second matrix vector product for any electrically small elements that are included in the structure. These matrix vector products for the electrically large and small elements are combined, and a net delta for a combination of the matrix vector products is determined. The iteration continues until a net delta is obtained that is within the predefined limits. The matrix vector products that were last obtained are used to solve for the desired parameter. The solution for the desired parameter is then presented to a user in a tangible form; for example, on a display.

Jandhyala, Vikram↗

Gravitational lens models of arcs in clusters

It is now well established that the luminous arcs discovered in clusters of galaxies, in particular those in Abell 370 and Cluster 2244-02, are produced by gravitational lensing of background sources. The arcs are modeled and constraints are placed on the distribution of the mass in the clusters and the shape and size of the sources. The models require, as expected, a large amount of dark matter in the clusters and a mass-to blue-light ratio for the cluster which exceeds 100 solar mass/solar luminosity and could be as high as 1000 solar mass/solar luminosity depending on cosmological parameters and the distribution of the dark matter. Furthermore, it is found that in the case of the arc in A370 the dark matter must have a different distribution than the luminous galaxies, while for the arc in Cl 2244 the dark matter can have a distribution similar to that of the light matter (galaxies) or a separate distribution.

Bergmann, Anton G.↗

Electrostatic waves due to field-aligned electron beams in the low-latitude boundary layer

Mass-resolved ion, electron, and plasma wave data obtained from several low-latitude boundary layer (LLBL) crossings by the AMPTE CCE satellite are analyzed. The data clearly separate the LLBL from the adjacent magnetosheath and magnetosphere. Attention was focused on wave-particle interactions involving electrons. Electron beams were found to be present in the LLBL during the southward interplanetary magnetic field, along with a simultaneous enhancement of electrostatic waves with parallel polarization. Linear theory analysis shows that for plasma conditions in the LLBL, electron beams are unstable to electrostatic waves that propagate parallel to the local magnetic field, in agreement with observations. A numerical simulation study of the beam-plasma interaction in the LLBL shows that the instability saturates by thermalization of the beam but that a beamlike structure can still remain in the electron distribution for certain initial parameters. It is suggested that peaks in the electron velocity distribution function may be found in the LLBL away from the beam source region.

Peroomian, V.↗