SEARCH · Engineering Papers

Results for “dependability modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Modeling Shock-Separated Flow with the Active Model Split

The predictive accuracy of typical hybrid RANS/LES models is limited by several shortcomings, such as modeled-stress depletion and a dependence on scalar grid measures. The "active model-split" (AMS) hybrid RANS/LES model was developed by Haering, Oliver, and Moser to address these shortcomings. In this model, the mean and fluctuating portions of the stress are modeled separately; this allows for a consistent treatment of the mean flow even in the presence of resolved fluctuations. Part of this hybrid model is an active forcing, which generates turbulent fluctuations in regions that support LES-type fluctuations. In this presentation, some recent results are shown using the active model-split for shock-separated flows. The effects of the active forcing are examined, including the effects on anisotropy and acoustics. Superior predictive accuracy over RANS is shown for some quantities of interest, such as the shock location.

RANS↗

NASA NTRS1984

Martian crater dark streak lengths - Explanation from wind tunnel experiments

The mechanism determining the length of dark erosional (type-Id) streaks near craters and hills on the surface of Mars is investigated using Mariner-9 images and the results of wind-tunnel studies (Greeley et al., 1974; Iversen et al., 1975 and 1976). An empirical model of time-dependent growth of crater-wake streaks in the downwind direction involving deflation by sand-grain saltation is developed and shown to give relatively good agreement with streak growth rates observed over a 38-day period of Mariner observations and with the streak-length/crater-diameter and streak-length/crater-height ratios analyzed statistically by Lee (1984). This time-dependent mechanism is thus considered a reasonable alternative (for dark streaks) to the blocking model proposed by Lee for both bright and dark streaks. Preliminary analysis of further wind-tunnel experiments (Iversen et al., 1982) supports the blocking model of bright-streak formation.

Iversen, J. D.↗

NASA NTRS1987

Development of control strategies for safe microburst penetration: A progress report

A single-engine, propeller-driven, general-aviation model was incorporated into the nonlinear simulation and into the linear analysis of root loci and frequency response. Full-scale wind tunnel data provided its aerodynamic model, and the thrust model included the airspeed dependent effects of power and propeller efficiency. Also, the parameters of the Jet Transport model were changed to correspond more closely to the Boeing 727. In order to study their effects on steady-state repsonse to vertical wind inputs, altitude and total specific energy (air-relative and inertial) feedback capabilities were added to the nonlinear and linear models. Multiloop system design goals were defined. Attempts were made to develop controllers which achieved these goals.

Psiaki, Mark L.↗

NASA NTRS2018

VLTI-GRAVITY Measurements of Cool Evolved Stars I. Variable Photosphere and Extended Atmosphere of the Mira Star R Peg

Context. Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Aim. Here we measure the variability of the oxygen-rich Mira star R Peg in near-continuum and molecular bands. Methods. We used near-infrared K-band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the Very Large Telescope Interferometer (VLTI). Results. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.24 m are steadily increasing with values of 8.7 0.1 mas, 9.4 0.1 mas, 9.8 0.1 mas, and 9.9 0.1 mas at visual phases of 0.15, 0.36, 0.45, 0.53, respectively. UD diameters at a bandpass around 2.05 m, dominated by water vapour, follow the near-continuum variability at larger UD diameters between 10.7 mas and 11.7 mas. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters, with values between 12.3 mas and 11.7 mas. Conclusions. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes in the extended atmosphere along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

photospheric radius↗

NASA NTRS2020

Development of a Non-Iterative Balance Load Prediction Algorithm for the NASA Ames Unitary Plan Wind Tunnel

A non-iterative load prediction algorithm for strain-gage balances was developed for the NASA Ames Unitary Plan Wind Tunnels that computes balance loads from the electrical outputs of the balance bridges and a set of state variables. A state variable could be, for example, a balance temperature difference or the bellows pressure of a flow-through balance. The algorithm directly uses regression models of the balance loads for the load prediction that were obtained by applying global regression analysis to balance calibration data. This choice greatly simplifies both implementation and use of the load prediction process for complex balance configurations as no load iteration needs to be performed. The regression model of a balance load is constructed by using terms from a total of nine term groups. Four term groups are derived from a Taylor Series expansion of the relationship between the load, gage outputs, and state variables. The remaining five term groups are defined by using absolute values of the gage outputs and state variables. Terms from these groups should only be included in the regression model if calibration data from a balance with known bi-directional outputs is analyzed. It is illustrated in detail how global regression analysis may be applied to obtain the coefficients of the chosen regression model of a load component assuming that no linear or massive near-linear dependencies between the regression model terms exist. Data from the machine calibration of a six-component force balance is used to illustrate both application and accuracy of the non-iterative load prediction process.

Ulbrich, Norbert M.↗

NASA NTRS2022

Assurance Equations: A Cost and Criticality Model for Optimizing Quality Assurance Surveillance

The cost of quality vs cost of failure correction has been a long-running topic of discussion within the Aerospace community. It leads directly to concepts of “risk tolerance”, and risk-based decision-making. It would be valuable if there was a way to compute the optimal investment in customer-executed quality assurance activities using defect significance with respect to performance objectives, the activities’ defect detection effectiveness, and the cost-penalty for late discovery of impactful defects. This optimization is particularly of interest to projects whose budget constraints significantly limit their risk management options.The cost to fix defects (i.e., failure correction) escalates as the project matures. There have been studies attempting to determine the relative cost of fixing defects discovered during various phases of a project life cycle with important implications, all of which suggest growth factors are large. The commonly referred to 1:10:100 rule represents a cost multiplier for repair/rework across the Design to Fab to Test hardware development phases. Cost premiums for QA activities also accumulate when they are treated as mandatory (due to schedule drag) or are performed later than their assigned phase.This paper describes the modeling of development phase -dependencies in the conduct of typical customer-executed quality assurance activities. Our initial modeling encompasses:• Distinct phases of the production lifecycle• Multiple kinds of Defects, each with some a-priori likelihood of being present• Each defect’s impact on performance Objectives for a type of hardware• The cost and efficacy of assurance techniques at detecting such Defects• The costs of fixing those Defects detected in a given phase of the production lifecycleThe model captures assurance activities’ abilities to Detect defects. Upon detection it is assumed that the Defect is immediately fixed. Defects that “escape” detection by some activity may thereafter be detected by a later activity, but by then the cost of fixing the Defect may have escalated. Defects are related to the performance Objectives they would detract from, were those Defects to remain present in the operating system.We have constructed and are exploring, a model that relates the importance of hardware system elements to mission objectives, the impact of types of Defects on those hardware types, the cost of customer-executed assurance activities (i.e., supplier controls) and their effectiveness towards reducing an impactful quality escape, and the cost of Defect correction across production phase. We describe the approach taken to select the key model aspects, why they are relevant to our NASA mission, and our efforts to populate it with relevant and contemporary data. We use a notional example to illustrate model design and function.

Plante, Jeannette↗

NASA NTRS1983

Application of an Uncoupled Elastic-plastic-creep Constitutive Model to Metals at High Temperature

A uniaxial, uncoupled constitutive model to predict the response of thermal and rate dependent elastic-plastic material behavior is presented. The model is based on an incremental classicial plasticity theory extended to account for thermal, creep, and transient temperature conditions. Revisions to he combined hardening rule of the theory allow for better representation of cyclic phenomenon including the high rate of strain hardening upon cyclic reyield and cyclic saturation. An alternative approach is taken to model the rate dependent inelastic deformation which utilizes hysteresis loops and stress relaxation test data at various temperatures. The model is evaluated and compared to experiments which involve various thermal and mechanical load histories on 5086 aluminum alloy, 304 stainless steel and Hastelloy-X.

Haisler, W. E.↗

NASA NTRS1997

Frozen-Flux Modelling for Epochs 1915 and 1980

The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.

OBrien, Michael S.↗

NASA NTRS1990

The wavelength dependence of broadband linear polarization in cool stars

A model is developed to investigate the effect of wavelength dependence of magnetically generated broadband linear polarization (BLP) in late-type stars as a function of line blanketing, observation filters, and other variables. Linear polarization is assumed to be a single 'average' line and the line opacity, magnetic area filling factor, and continuum intensity are incorporated to establish an equation for disk-generated BLP. The method is investigated by computing special models for the flare star BY Dra and the sun which are compared to previous calculations. The model is shown to agree with the calculations of BLP, and the magnetic BLP dependent on wavelengths varies significantly from the simple power law dependences on wavelength. The model sets lower limits to the area occupied by magnetic regions, and an upper limit to the expected BLP also results from the formulation.

Saar, Steven H.↗

NASA NTRS2004

Small-Scale Variability of Large Cloud Drops

Cloud droplet size distribution is one of the most fundamental subjects in cloud physics. Understanding of spatial distribution and small-scale fluctuations of cloud droplets is essential for both cloud physics and atmospheric radiation. For cloud physics, it relates to the coalescence growth of raindrops while for radiation, it has a strong impact on a cloud's radiative properties. Most of the existing cloud radiation and precipitation formation models assume that the mean number of drops with a given radius varies proportionally to volume. The analysis of microphysical data on liquid water drop sizes shows that, for sufficiently small volumes, the number is proportional to the drop size dependent power of the volume. For abundant small drops present, the exponent is 1 as assumed in the conventional approach. However, for rarer large drops, the exponents fall below unity. At small scales, therefore, the mean number of large drops decreases with volume at a slower rate than the conventional approach assumes, suggesting more large drops at these scales than conventional models account for; their impact is consequently underestimated. Size dependent models of spatial distribution of cloud drops that simulate the observed power laws show strong drop clustering, the more so the larger the drops. The degree of clustering is determined by the observed exponents. The strong clustering of large drops arises naturally from the observed power-law statistics. Current theories of photon-cloud interaction and warm rain formation will need radical revision in order to produce these statistics; their underlying equations are unable to yield the observed power law.

Marshak, Alexander↗

NASA NTRS1987

Atomic oxygen in the lower thermosphere

The 63-micron line due to thermospheric atomic oxygen O(P-3), using a far-infrared spectrometer on a balloon platform at 37 km altitude over Palestine, TX, on June 20, 1983. From measurements of the equivalent width of this line at two elevation angles, a weak angular dependence is found: the equivalent width increases by a factor of 1.5 + or - 0.3 as the angle decreases from +30 deg to +1 deg. Since the optical depth of the O(P-3) line is large, the measured line intensity cannot be directly converted to a column abundance. Instead, the measurements are interpreted in terms of radiative transfer through a 16-layer atmosphere extending to 200 km. A model atmosphere for summer at 30 deg N, with an exospheric temperature of 1300 K, including an assumed daytime atomic oxygen abundance profile constructed from recent chemical and dynamical models and a water vapor abundance profile constructed from recent experimental and model results is used. For this assumed O(P-3) vertical profile shape a multiplicative scaling factor of 0.8, with an altitude-dependent uncertainty is determined. In the best-determined layer the uncertainty in the multiplier is + or - 0.2 at 119 km. The model-dependent peak atomic oxygen density is 3.6 (+ or - 1.9) x 10 to the 11th/cu cm at an altitude of about 101 km.

Lin, Florence J.↗

NASA NTRS1986

Temperature dependent nonlinear metal matrix laminae behavior

An analytical method is described for computing the nonlinear thermal and mechanical response of laminated plates. The material model focuses upon the behavior of metal matrix materials by relating the nonlinear composite response to plasticity effects in the matrix. The foundation of the analysis is the unidirectional material model which is used to compute the instantaneous properties of the lamina based upon the properties of the fibers and matrix. The unidirectional model assumes that the fibers properties are constant with temperature and assumes that the matrix can be modelled as a temperature dependent, bilinear, kinematically hardening material. An incremental approach is used to compute average stresses in the fibers and matrix caused by arbitrary mechanical and thermal loads. The layer model is incorporated in an incremental laminated plate theory to compute the nonlinear response of laminated metal matrix composites of general orientation and stacking sequence. The report includes comparisons of the method with other analytical approaches and compares theoretical calculations with measured experimental material behavior. A section is included which describes the limitations of the material model.

Barrett, D. J.↗

NASA NTRSDate not supplied

Predicting Maximum Temperatures of a Li-ion Battery on a Simulated Flight Profile using a Model-based Prognostics

One of the challenges in using Li-ion packs in aeronautics is their safety, and thermal runaway (TR) is a major concern. The current engineering solutions to prevent a Li-ion pack from a catastrophic TR require additional mass and volume to isolate cells. The excess mass could be reduced by improving detection and, thus, preventing a TR event. One of the possible early warning indicators of a TR is crossing a threshold temperature. We have developed an approach, based on the Unscented Kalman Filter (UKF), to predict the likelihood of reaching the threshold temperature for simulated flight profiles. The current battery prognostics algorithms for aerospace predict state-of-charge (SOC) and end-of-life (EOL) [1]. We extended this two-level algorithm to predict the maximum temperature during discharge. The amount of heat generated in a cell depends on factors such as cell chemistry, cell packaging, total cycles, operating temperature, and abuse history [2]. Our semi-empirical thermal model depends on three phenomenological parameters which account for those factors. In addition, a two-parameter reduced-order model is developed to predict the temperature rise for short bursts of “random-walk” (RW) discharge current sequence, which simulates a flight's current-loading profile. The performance of these models on different datasets and types of current loading will be presented. To predict the maximal temperatures for future cycles we must estimate the evolution of thermal parameters as the batteries age. It is found that the parameters of the 3-parametric thermal model cannot be estimated only from the RW data. To address the issue, we will present two alternative approaches: i) expanding the datasets to include discharge profiles beyond RWs; ii) model reduction to a two-parametric model. The two approaches will be illustrated by an application to the cycling data from a commercial LG 18650 cell. References: 1. M. Daigle, C.S. Kulkarni, End-of-discharge and End-of-life Prediction in Lithium-ion Batteries with Electrochemistry-based Aging Models, in: AIAA Infotech @ Aerospace, American Institute of Aeronautics and Astronautics, San Diego, California, USA, 2016. 2. M. Börner, et. al, Correlation of aging and thermal stability of commercial 18650-type lithium ion batteries, Journal of Power Sources. 342 (2017) 382–392.

Thermal runaway↗

NASA NTRS1993

A global analysis of the ozone deficit in the upper stratosphere and lower mesosphere

The global measurements of temperature, ozone, water vapor, and nitrogen dioxide acquired by the Limb Infrared Monitor of the Stratosphere (LIMS), supplemented by a precomputed distribution of chlorine monoxide, are used to test the balance between odd oxygen production and loss in the upper stratosphere and lower mesosphere. An efficient photochemical equilibrium model, whose validity is ascertained by comparison with the results from a fully time-dependent one-dimensional model at selected latitudes, is used in the calculations. The computed ozone abundances are systematically lower than observations for May 1-7, 1979, which suggests, contrary to the conclusions of other recent studies, a real problem in model simulations of stratospheric ozone.

Eluszkiewicz, Janusz↗

NASA NTRS2018

3D Mapping of the Neutral X-Ray Absorption in the Local Interstellar Medium: The Gaia and XMM-Newton Synergy

We present a three-dimensional map of the hydrogen density distribution in the Galactic interstellar medium. The hydrogen-equivalent column densities were obtained from the Exploring the X-ray Transient and variable Sky project (EXTRAS) which provides equivalent NH values from X-ray spectral fits of observations within the XMM-Newton Data Release. EXTRAS include multiple fits for each source, allowing an accurate determination of the equivalent column densities, which depends on the continuum modelling of the spectra. A cross-correlation between the EXTRAS catalogue and the first Gaia Data Release was performed in order to obtain accurate parallax and distance measurements. We use a Bayesian method explained in Rezaei Kh. et al. (2017) in order to predict the most probable distribution of the density at any arbitrary point, even for lines of sight along which there are no initial observation. The resulting map shows small-scale density structures which could not have been modelled by using analytic density profiles. In this paper, we present a proof of concept of the kind of science possible with the synergy of these catalogues. However, given the systematic uncertainties connected to the source identification and to the dependence of NH on the spectral model, the present maps should be considered qualitatively at this point.

Gatuzz, Efrain↗

NASA NTRS1993

Heliocentric zoning of the asteroid belt by aluminum-26 heating

Variations in petrology among meteorites attest to a strong heating event early in solar system history, but the heat source has remained unresolved. Aluminum-26 has been considered the most likely high-energy, short-lived radionuclide (half-life 0.72 million years) since the discovery of its decay product - excess Mg-26 - in Allende CAI's. Furthermore, observation of relict Mg-26 in an achondritic clast and in feldspars within ordinary chondrites (3,4) provided strong evidence for live Al-26 in meteorite parent bodies and not just in refractory nebular condensates. The inferred amount of Al-26 is consistent with constraints on the thermal evolution of both ordinary and carbonaceous chondrite parent objects up to a few hundred kilometers in diameter. Meteorites can constrain the early thermal evolution of their parent body locations, provided that a link can be established between asteroid spectrophotometric signature and meteorite class. Asteroid compositions are heliocentrically distributed: objects thought to have experienced high metamorphic or even melting temperatures are located closer to the sun, whereas apparently unaltered or mildly heated asteroids are located farther away. Heliocentric zoning could be the result of Al-26 heating if the initial amount of the radionuclide incorporated into planetesimals was controlled by accretion time, which in turn varies with semimajor axis. Analytic expressions for planetary accretion may be integrated to given the time, tau, required for a planetesimal to grow to a specified radius: tau varies as a(sup n), where n = 1.5 to 3 depending on the assumptions about variations in the surface density of the planetesimal swarm. Numerical simulations of planetesimal accretion at fixed semimajor axis demonstrate that variations in accretion time among small planetesimals can be strongly nonlinear depending on the initial conditions and model assumptions. The general relationship with semimajor axis remains valid because it depends only on the initial orbit properties and distribution of the planesimal swarm. In order to demonstrate the basic dependence of thermal evolution on semimajor axis, we parameterized accretion time across the asteroid belt according to tau varies as a(sup n) and calculated the subsequent thermal history. Objects at a specified semimajor axis were assumed to have the same accretion time, regardless of size. We set the initial Al-26/Al-27 ratio = 6 x 10(exp -5) and treated n and tau(sub 0) at a(sub 0) = 3 AU as adjustable parameters. The thermal model included temperature-dependent properties of ice and rock (CM chondrite analog) and the thermodynamic effects of phase transitions.

Grimm, R. E.↗

NASA NTRS2002

Higher Order Modulation Intersymbol Interference Caused by Traveling-wave Tube Amplifiers

For the first time, a time-dependent, physics-based computational model has been used to provide a direct description of the effects of the traveling wave tube amplifier (TWTA) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves, Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry and operating characteristics of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept-amplitude and/or swept-frequency data. First, the TWT model using the three dimensional (3D) electromagnetic code MAFIA is presented. Then, this comprehensive model is used to investigate approximations made in conventional TWT black-box models used in communication system level simulations, To quantitatively demonstrate the effects these approximations have on digital signal performance predictions, including intersymbol interference (ISI), the MAFIA results are compared to the system level analysis tool, Signal Processing, Workstation (SPW), using high order modulation schemes including 16 and 64-QAM.

Kory, Carol L.↗

NASA NTRS2017

Effects of Adiabatic Heating on the High Strain Rate Deformation of Polymer Matrix Composites

Polymer matrix composites (PMCs) are increasingly being used in aerospace structures that are expected to experience complex dynamic loading conditions throughout their lifetime. As such, a detailed understanding of the high strain rate behavior of the constituents, particularly the strain rate, temperature, and pressure dependent polymer matrix, is paramount. In this paper, preliminary efforts in modeling experimentally observed temperature rises due to plastic deformation in PMCs subjected to dynamic loading are presented. To this end, an existing isothermal viscoplastic polymer constitutive formulation is extended to model adiabatic conditions by incorporating temperature dependent elastic properties and modifying the components of the inelastic strain rate tensor to explicitly depend on temperature. It is demonstrated that the modified polymer constitutive model is capable of capturing strain rate and temperature dependent yield as well as thermal softening associated with the conversion of plastic work to heat at high rates of strain. The modified constitutive model is then embedded within a strength of materials based micromechanics framework to investigate the manifestation of matrix thermal softening, due to the conversion of plastic work to heat, on the high strain rate response of a T700Epon 862 (T700E862) unidirectional composite. Adiabatic model predictions for high strain rate composite longitudinal tensile, transverse tensile, and in-plane shear loading are presented. Results show a substantial deviation from isothermal conditions; significant thermal softening is observed for matrix dominated deformation modes (transverse tension and in-plane shear), highlighting the importance of accounting for the conversion of plastic work to heat in the polymer matrix in the high strain rate analysis of PMC structures.

polymer matrix composites↗