Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 217 records · Page 12

Estimation of high temperature low cycle fatigue on the basis of inelastic strain and strainrate

Fatigue life at elevated temperature can be predicted by introducing parametric values obtained from monotonic constitutive behavior into the Universal-Slopes Equation. For directionally solidified MAR-M200+HF at 975 C, these parameters are the maximum stress achievable under entirely plastic (time-independent) and purely creep (time-dependent) conditions and the corresponding inelastic strains, as well as the elastic modulus. For materials which exhibit plasticity/creep interaction, two more pairs of monotonic parameters must be evaluated for fatigue life prediction. This life-prediction method based on the Universal-Slopes Equation, resulted from a constitutive model characterizing monotonic and cyclic data as inelastic strainrate as a function of inelastic strain. Characterizing monotonic data is this way, permitted distinction between different material responses such as strain-hardening, strain-softening, and dynamic recovery effects. Understanding and defining the region of influence of each of these effects facilitated formulation of the constitutive model in relation to the mechanical and microstructural processes occurring in the material under cyclic loading.

Berkovits, A.↗

On reflection of Alfven waves in the solar wind

We have revisited the problem of propagation of toroidal and linear Alfven waves formulated by Heinemann and Olbert (1980) to compare WKB and non-WKB waves and their effects on the solar wind. They considered two solar wind models and showed that reflection is important for Alfven waves with periods of the order of one day and longer, and that non-WKB Alfven waves are no more effective in accelerating the solar wind than WKB waves. There are several recently published papers which seem to indicate that Alfven waves with periods of the order of several minutes should be treated as non-WKB waves and that these non-WKB waves exert a stronger acceleration force than WKB waves. The purpose of this paper is to study the origin of these discrepancies by performing parametric studies of the behavior of the waves under a variety of different conditions. In addition, we want to investigate two problems that have not been addressed by Heinemann and Olbert, namely, calculate the efficiency of Alfven wave reflection by using the reflection coefficient and identify the region of strongest wave reflection in different wind models. To achieve these goals, we investigated the influence of temperature, electron density distribution, wind velocity and magnetic field strength on the waves. The obtained results clearly demonstrate that Alfven wave reflection is strongly model dependent and that the strongest reflection can be expected in models with the base temperatures higher than 10(exp 6) K and with the base densities lower than 7 x 10(exp 7) cm(exp -3). In these models as well as in the models with lower temperatures and higher densities, Alfven waves with periods as short as several minutes have negligible reflection so that they can be treated as WKB waves; however, for Alfven waves with periods of the order of one hour or longer reflection is significant, requiring a non-WKB treatment. We also show that non-WKB, linear Alfven waves are always less effective in accelerating the plasma than WKB Alfven waves. Finally, it is evident from our results that the region of strongest wave reflection is usually located at the base of the models, and hence that interpretation of wave reflection based solely on the reflection coefficient can be misleading.

Krogulec, M.↗

Reflection of Alfven waves in the solar wind

We have revisited the problem of propagation of toroidal and linear Alfven waves formulated by Heinemann and Olbert (1980) to compare Wentzel-Kramers-Brillouin (WKB) and non-WKB waves and their effects on the solar wind. They considered two solar wind models and showed that reflection is important for Alfven waves with periods of the order of one day and longer and that non-WKB Alfven waves are no more effective in accelerating the solar wind than in WKB waves. There are several recently published papers that seem to indicate that Alfven waves with periods of the order of several minutes should be treated as non-WKB waves and that these non-WKB waves exert a stronger acceleration force than WKB waves. The purposse of this paper is to study the origin of these discrepancies by performing parametric studies of the behavior of the waves under a variety of different conditions. In addition, we want to investigate two problems that have not been addressed by Heinimann and Olbert, namely, calculate the efficieny of Alfven wave reflection by using the reflection coefficient and identfy the region of strongest wave reflection in different wind models. To achieve these goals, we investigate the influence of temperature, electron desity distribution, wind velocity, and magnetic field strength on te waves. The obtained results clearly demonstrate that Alfven wave reflection is strongly model dependent and that the strongest reflection can be expected in models with the base temperatures higher than 10(exp 6) K and with the base densities lower than 7 x 10(exp 7)/cu cm. In these models as well as in the models with lower temperatures and higher densities Alfven waves with periods as short as several minutes have negligible reflection so that they can be treated as WKB waves; however, for Alfven waves with periods of the order of one hour or longer reflection is significant, requiring a non-WKB treatment. We also show that non-WKB, linear Alfven waves are always less effective in accelerating the plasma than WKB Alfven waves. Finally, it is evident from our results that the region of strongest wave reflection is usually located at the base of the models and hence that interpretation of wave reflection based soley on the reflection coefficient can be misleading.

Krogulec, M.↗

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗

Parametric Study of Afterbody/nozzle Drag on Twin Two-dimensional Convergent-divergent Nozzles at Mach Numbers from 0.60 to 1.20

An investigation has been conducted in the Langley 16-Foot Transonic Tunnel to determine the effects of upper and lower external nozzle flap geometry on the external afterbody/nozzle drag of nonaxisymmetric two-dimensional convergent-divergent exhaust nozzles having parallel external sidewalls installed on a generic twin-engine, fighter-aircraft model. Tests were conducted over a Mach number range from 0.60 to 1.20 and over an angle-of-attack range from -5 to 9 deg. Nozzle pressure ratio was varied from jet off (1.0) to approximately 10.0, depending on Mach number.

Pendergraft, Odis C., Jr.↗

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↗

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗

Romie: A Domain-Independent Tool for Computer-Aided Robust Operations Management

Romie is a decision support tool based on AI's latest advances in the domain of robust scheduling. Unlike all its predecessors, the tool allows to (i) visually model the operational problem and context entirely (ii) optimize to find near-optimal schedules while taking uncertainty into account and deals with (iii) a combination of various {key performance indicators (KPIs). It comes with a web user interface. Part or all of the modelled activities may be associated to random variables describing their stochastic durations, in order to produce schedules that are robust w.r.t. temporal uncertainty. Hence, depending on the pursued KPIs, the schedules maximize a combination of the following terms: the probability of satisfying the problem constraints, the expected return/efficiency, the expected outcome quality, and even the operators' wellness by minimizing its expected extra-hours. Initially developed for spatial exploration and demonstration in the context of Mars analog missions, this versatile tool is here applied to operations management in both biotechnology manufacturing and robots parametrization in a cave exploration context.

Chien, Steve A.↗

The Impact of Inclination-dependent Attenuation on Ultraviolet Star Formation Rate Tracers

We examine and quantify how hybrid (e.g., UV+IR) star formation rate (SFR) estimators and the AFUV–β relation depend on inclination for disk-dominated galaxies using spectral energy distribution modeling that utilizes the inclination-dependent attenuation curves described in Doore et al. We perform this analysis on a sample of 133 disk-dominated galaxies from the CANDELS fields and 18 disk galaxies from the Spitzer Infrared Nearby Galaxies Survey and Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel samples. We find that both the hybrid SFR estimators and the AFUV–β relation present clear dependencies on inclination. To quantify this dependence in the hybrid SFR estimators, we derive an inclination and a far-UV–near-IR color-dependent parametric relation for converting observed UV and IR luminosities into SFRs. For the AFUV–β relation, we introduce an inclination-dependent component that accounts for the majority of the inclination dependence with the scatter of the relation increasing with inclination. We then compare both of these inclination-dependent relations to similar inclination-independent relations found in the literature. From this comparison, we find that the UV+IR correction factor and AFUV for our hybrid and AFUV–β relations, respectively, result in a reduction in the residual scatter of our sample by approximately a factor of 2. Therefore, we demonstrate that inclination must be considered in hybrid SFR estimators and the AFUV–β relation to produce more accurate SFR estimates in disk-dominated galaxies.

Keith Doore↗

Construction of Covariance Functions with Variable Length Fields

This article focuses on construction, directly in physical space, of three-dimensional covariance functions parametrized by a tunable length field, and on an application of this theory to reproduce the Quasi-Biennial Oscillation (QBO) in the Goddard Earth Observing System, Version 4 (GEOS-4) data assimilation system. These Covariance models are referred to as multi-level or nonseparable, to associate them with the application where a multi-level covariance with a large troposphere to stratosphere length field gradient is used to reproduce the QBO from sparse radiosonde observations in the tropical lower stratosphere. The multi-level covariance functions extend well-known single level covariance functions depending only on a length scale. Generalizations of the first- and third-order autoregressive covariances in three dimensions are given, providing multi-level covariances with zero and three derivatives at zero separation, respectively. Multi-level piecewise rational covariances with two continuous derivatives at zero separation are also provided. Multi-level powerlaw covariances are constructed with continuous derivatives of all orders. Additional multi-level covariance functions are constructed using the Schur product of single and multi-level covariance functions. A multi-level powerlaw covariance used to reproduce the QBO in GEOS-4 is described along with details of the assimilation experiments. The new covariance model is shown to represent the vertical wind shear associated with the QBO much more effectively than in the baseline GEOS-4 system.

Gaspari, Gregory↗

Measurements of DSD Second Moment Based on Laser Extinction

Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through rain (or dust in the rocket case) yields an estimate of the 2nd moment of the particle cloud, and rainfall drop size distribution (DSD) in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rainfall make in direct measurements of the DSD. Most common of these instruments are the rainfall rate gauge measuring the 1 1/3 th moment, (when using a D(exp 2/3) dependency on terminal velocity). Instruments that scatter microwaves off of hydrometeors, such as the WSR-880, vertical wind profilers, and microwave disdrometers, measure the 6th moment of the DSD. By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain, yield a measurement of the DSD 2nd moment, using the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required, depending on the intensity of the rainfall rate. For moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. A photo-detector could replace the camera, for automated installations. In order to spatially correlate the 2nd moment measurements to a collocated disdrometer or tipping bucket, the laser's beam path can be reflected multiple times using mirrors to restrict the spatial extent of the measurement. In cases where a disdrometer is not available, complete DSD estimates can be produced by parametric fitting of DSD model to the 2nd moment data in conjunction with tipping bucket data. In cases where a disdrometer is collocated, the laser extinction technique may yield a significant improvement to insitu disdrometer validation and calibration strategies

Lane, John E.↗

Contamination sensitivity of typical mirror coatings - A parametric study

Contamination budgeting for space optical systems basically starts at understanding the sensitivity of component performance, e.g. mirror reflectance, window transmittance, etc., to surface deposits. To evaluate contamination sensitivity for mirror coatings, eight types representative of those used in the vacuum ultraviolet, visible, and infrared were modeled assuming that the contaminant is uniformly deposited on the mirror surface. Parametric studies over a range of complex refractive indices combined with an examination of optical data available for several organic materials suggested division of the contaminant layer index into three categories, N = 1.5 + 0.1i, N = 1.5 + 0.5i, and N = 1.5 + 2.0i. Contaminant thickness sensitivity curves were then calculated for each of the selected mirror coatings. For comparative purposes, critical thicknesses for each type were extracted, assuming a reflectance loss of 10 percent was allowable. Critical thicknesses ranged from about 10 A to 1000 A depending on the specifics of coating design and spectral region.

Osantowski, J. F.↗

Effect of Scattering on the Heat Transfer Behavior of a Typical Semitransparent TBC Material on a Substrate

A parametric study was performed to examine the effects of isotropic scattering on the heat transfer in typical semitransparent thermal barrier coating on an opaque substrate. Some ceramic materials are semitransparent in the wavelength range where thermal radiation is important. Therefore, absorption, emission, and scattering of thermal radiation by the semitransparent layer and the emissivity of the substrate will affect the heat transfer and temperature in the layer. Scattering which depends on the structure of the semitransparent material can be used to increase the reflectivity of the layer and therefore decrease the heat transfer through the layer. A one dimensional model of a 1 mm thick semitransparent layer on a substrate was used in this study. The front of the semitransparent layer is heated by convection and radiation and the back of the substrate is cooled by radiation and convection. The coating is assumed to be semitransparent up to a cutoff wavelength and opaque at higher wavelengths. The absorption and scattering coefficient are constant over the wavelength range where the material is semitransparent. The absorption coefficient, scattering coefficients, width of the semitransparent band, and the bond coat emissivity are varied. Temperature profiles and heat flux through layers as a function of absorption and scattering coefficients are presented.

Spuckler, Charles M.↗

Autoignition of hydrogen injected transverse to a supersonic airstream

A parametric study of autoignition characteristics of transverse fuel injector patterns for a hydrogen-fueled, supersonic-combustion ramjet was conducted in direct-connect tests over a range of simulated hypersonic flight conditions. Autoignition limits, defined directly by visual observations and indirectly by measured wall pressure and temperature, are compared with a semi-empirical ignition limit model which accounts for first order effects on autoignition. It is shown that to initiate combustion in the immediate region of the fuel injector, ignition must occur in the separation region ahead of the fuel jet. Autoignition within this separation zone exhibits strong dependence on flow pressure and temperature, injector diameter, local boundary layer energy thickness, wall temperature, and on fluid dynamic interactions with adjacent downstream facing steps (such as inlet isolation steps). Secondary effects are also discussed. The empirical model is used to define limits of test temperature, pressure, and scale required for autoignition in a current subscale engine test program.

Mcclinton, C. R.↗

New Gravity Wave Treatments for GISS Climate Models

Previous versions of GISS climate models have either used formulations of Rayleigh drag to represent unresolved gravity wave interactions with the model resolved flow or have included a rather complicated treatment of unresolved gravity waves that, while being climate interactive, involved the specification of a relatively large number of parameters that were not well constrained by observations and also was computationally very expensive. Here, we introduce a relatively simple and computationally efficient specification of unresolved orographic and non-orographic gravity waves and their interaction with the resolved flow. We show comparisons of the GISS model winds and temperatures with no gravity wave parametrization; with only orographic gravity wave parameterization; and with both orographic and non-orographic gravity wave parameterizations to illustrate how the zonal mean winds and temperatures converge toward observations. We also show that the specifications of orographic and nonorographic gravity waves must be different in the Northern and Southern Hemispheres. We then show results where the non-orographic gravity wave sources are specified to represent sources from convection in the Intertropical Convergence Zone and spontaneous emission from jet imbalances. Finally, we suggest a strategy to include these effects in a climate dependent manner.

Geller, Marvin A.↗

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗

Experimental investigation of helicopter vibration reduction using rotor blade aeroelastic tailoring

A wind tunnel investigation has been conducted to parametrically investigate the effect of blade nonstructural mass on helicopter fixed- and rotating-system vibratory loads. The data were obtained using Mach- and aeroelastically-scaled model rotor blades which allowed for the addition of concentrated nonstructural masses at multiple points along the blade radius. Testing was conducted for advance ratios ranging from 0.10 to 0.35 for ten blade mass configurations. Three thrust levels were obtained at representative full-scale shaft angles for each configuration. Results indicate that proper placement of blade nonstructural mass can provide reductions in fixed-system vibratory loads, but that correct mass placement and the loads reduction realized are dependent upon flight condition. The data base obtained with this investigation provides a comprehensive set of fixed-system shears and moments, blade moment, and blade flap and lag response. The data set is well suited for use in the correlation and development of advanced rotorcraft analyses.

Wilbur, Matthew L.↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗