Engineering Papers⌕ Search

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.

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

At least 811 records · Page 45

Determination of Ionospheric Conductivities from UVI Intensity Ratios

Germany et al. examined the viability of using LBH emission ratios to infer auroral conductances, using modeled line brightnesses representative of the (then) planned UVI bandpasses. Here we extend that work by using actual UVI bandpass information instead of single modeled lines. Conductances are calculated using combined two stream and Field Line Interhemispheric Plasma (FLIP) model calculations to specify the ionosphere. The dependence of modeled conductances to UVI bandpass intensities is examined as a function of incident average energy, total energy flux, and changes in solar and magnetic activity levels. Potential parameterizations of conductance with UVI intensity ratios will be examined, as will error estimates of the conductance determinations.

Spann, James F., Jr.↗

Space Environments and Effects: Trapped Proton Model

An improved model of the Earth's trapped proton environment has been developed. This model, designated Trapped Proton Model version 1 (TPM-1), determines the omnidirectional flux of protons with energy between 1 and 100 MeV throughout near-Earth space. The model also incorporates a true solar cycle dependence. The model consists of several data files and computer software to read them. There are three versions of the mo'del: a FORTRAN-Callable library, a stand-alone model, and a Web-based model.

Huston, S. L.↗

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗

Determining the standoff distance of the bow shock: Mach number dependence and use of models

We explore the factors that determine the bow shock standoff distance. These factors include the parameters of the solar wind, as well as the size and shape of the obstacle. In this report we develop a semiempirical Mach number relation for the bow shock standoff distance in order to take into account the shock's behavior at low Mach numbers. This is done by determining which properties of the shock are most important in controlling the standoff distance and using this knowledge to modify the current Mach number relation. While the present relation has proven useful at higher Mach numbers, it has lacked effectiveness at the low Mach number limit. We also analyze the bow shock dependence upon the size and shape of the obstacle, noting that it is most appropriate to compare the standoff distance of the bow shock to the radius of curvature of the obstacle, as opposed to the distance from the focus of the object to the nose. Last, we focus our attention on the use of bow shock models in determining the standoff distance. We note that the physical behavior of the shock must correctly be taken into account, specifically the behavior as a function of solar wind dynamic pressure; otherwise, erroneous results can be obtained for the bow shock standoff distance.

Farris, M. H.↗

Investigation of Hypersonic Laminar Heating Augmentation in the Stagnation Region

Laminar stagnation region heating augmentation is investigated in the AEDC Tunnel 9 at Mach 10 by performing high frequency surface pressure and heat transfer measurements on the Orion CEV capsule at zero degree angle-of-attack for unit Reynolds numbers between 0.5 and 15 million per foot. Heating augmentation increases with Reynolds number, but is also model size dependent as it is absent on a 1.25-inch diameter model at Reynolds numbers where it reaches up to 15% on a 7-inch model. Heat transfer space-time correlations on the 7-inch model show that disturbances convect at the boundary layer edge velocity and that the streamwise integral scale increases with distance. Therefore, vorticity amplification due to stretching and piling-up in the stagnation region appears to be responsible for the stagnation point heating augmentation on the larger model. This assumption is reinforced by the f(exp -11/3) dependence of the surface pressure spectrum compared to the f(exp -1) dependence in the free stream. Vorticity amplification does not occur on the 1.25- inch model because the disturbances are too large. Improved free stream fluctuation measurements will be required to determine if significant vorticity is present upstream or mostly generated behind the bow shock.

Marineau, Eric C.↗

Linear and nonlinear responses to middle latitude surface temperature anomalies

Nonlinear responses to fixed and seasonally varying surface temperature anomalies in a two-level hemispheric time-dependent coupled atmosphere-surface mixed layer model are described. Linear stationary models that are equivalent to the nonlinear time-dependent model are used to analyze these responses. A model linearized around the climatological zonal state of the time dependent model and forced by anomalous surface temperatures does not provide a reasonable estimate for the anomalous reponses, which are considerably underestimated. Better responses are obtained when the anomalous stationary nonlinear eddy fluxes are included in a model linearized around the full climatology. However, this latter model is overly sensitive, and anomalous responses are a small residual balance to the forcing by the surface temperature anomalies and the anomalous transient eddy fluxes. To better understand these linear responses, an eigenanalysis of the climatological state is performed. Seasonal anomalies appear to be dominated by one characteristic pattern near resonance which can be associated with a slowly growing coupled atmosphere-ocean instability.

Roads, John O.↗

Eigenvalue assignment strategies in rotor systems

The work done to establish the control and direction of effective eigenvalue excursions of lightly damped, speed dependent rotor systems using passive control is discussed. Both second order and sixth order bi-axis, quasi-linear, speed dependent generic models were investigated. In every case a single, bi-directional control bearing was used in a passive feedback stabilization loop to resist modal destabilization above the rotor critical speed. Assuming incomplete state measurement, sub-optimal control strategies were used to define the preferred location of the control bearing, the most effective measurement locations, and the best set of control gains to extend the speed range of stable operation. Speed dependent control gains were found by Powell's method to maximize the minimum modal damping ratio for the speed dependent linear model. An increase of 300 percent in stable speed operation was obtained for the sixth order linear system using passive control. Simulations were run to examine the effectiveness of the linear control law on nonlinear rotor models with bearing deadband. The maximum level of control effort (force) required by the control bearing to stabilize the rotor at speeds above the critical was determined for the models with bearing deadband.

Youngblood, J. N.↗

A semi-empirical model for the complete orientation dependence of the growth rate for vapor phase epitaxy - Chloride VPE of GaAs

A semi-empirical model has been developed to determine the complete crystallographic orientation dependence of the growth rate for vapor phase epitaxy (VPE). Previous researchers have been able to determine this dependence for a limited range of orientations; however, our model yields relative growth rate information for any orientation. This model for diamond and zincblende structure materials is based on experimental growth rate data, gas phase diffusion, and surface reactions. Data for GaAs chloride VPE is used to illustrate the model. The resulting growth rate polar diagrams are used in conjunction with Wulff constructions to simulate epitaxial layer shapes as grown on patterned substrates. In general, this model can be applied to a variety of materials and vapor phase epitaxy systems.

Seidel-Salinas, L. K.↗

Preliminary Error Characterization and Parametric Error Model for the Automatic Dependent Surveillance - Contract Extended Projected Profile Message

A critical component of Trajectory-Based Operations (TBO) is the ability for a consistent and accurate 4-dimensional trajectory (4DT) to be shared and synchronized between airborne and ground systems as well as amongst various ground automation systems. The Aeronautical Telecommunication Network—Baseline 2 (ATN-B2) standard defines the Extended Projected Profile (EPP) trajectory that can be sent via Automatic Dependent Surveillance-Contract (ADS-C) from an aircraft to ground automation. The EPP trajectory message contains a representation of the reference trajectory from an aircraft’s Flight Management System (FMS). In this work, a set of scenarios were run in a high-fidelity aircraft and FMS simulation to perform an initial characterization of EPP trajectory errors under a given set of conditions. The parameters investigated were the route length, route type, wind magnitude error, wind direction error, and with and without a required time-of-arrival (RTA) constraint. In addition, linear regression was used to identify a set of EPP error models for the cross-track, vertical, and time errors.

Guerreiro, Nelson M.↗

Large-Eddy Simulations of Idealized Atmospheric Boundary Layers Using Nalu-Wind

Accurate prediction of wind-plant performance relies, in part, on properly characterizing the turbulent atmospheric boundary layer (ABL) flow in which wind turbines operate. Large-eddy simulation (LES) is a powerful tool for simulating ABLs because it resolves the largest, most energetic scales of three-dimensional turbulent motions. Yet LES predictions are well known to depend on modeling choices such as grid resolution, numerical discretization schemes, and closures for unresolved scales of turbulence. Here, we evaluate how these choices influence predictions of ABL winds using Nalu-Wind, a wind-specific fork of the open-source, generalized, unstructured, massively parallel flow solver NaluCFD/Nalu.

17 WIND ENERGY↗

Models for the Infrared Cavity of HH 46/47

We have modeled the limb-brightened cavity seen in the new Spitzer Space Telescope IR images of the southwest lobe of HH 46/47 as the bow shock driven by an outflow from a young, low-mass star. We present models in which the outflow is a perfectly collimated, straight jet, models in which the jet precesses, and finally a model in which the outflow takes the form of a latitude-dependent wind. We study cases in which the outflow moves into a constant-density cloud and into a stratified cloud. We find that the best agreement with the observed cavity is obtained for the precessing jet in a stratified cloud. However, the straight jet (traveling in a stratified cloud) also gives cavity shapes close to the observed one. The latitude-dependent wind model that we have computed gives cavity shapes that are substantially wider than the observed cavity. We therefore conclude that the cavity seen in the Spitzer observations of the southwest lobe of the HH 46/47 outflow do not seem to imply the presence of a latitude-dependent wind, as it can be modeled successfully with a perfectly collimated jet model.

Herbig Haro objects↗

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗

Using FU Orionis outbursts to constrain self-regulated protostellar disk models

One-dimensional, convective, vertical structure models and one dimensional time-dependent, radial diffusion models are combined to create a self-consistent picture in which FU Orionis outbursts occur in young stellar objects (YSOs) as the result of a large-scale, self-regulated, thermal ionization instability in the surrounding protostellar accretion disk. Although active accretion disks have long been postulated to be ubiqitous among low-mass YSOs, few constraints have until now been imposed on physical conditions in these disks. By fitting the results of time-dependent disk models to observed timescales of FU Orionis events, we estimate the magnitude of the effective viscous stress in the inner disk (r approximately less than 1 AU) to be, in accordance with an ad hoc 'alpha' prescription, the product of the local sound speed, pressure scale height, and an efficiency factor alpha of 10(exp -4) where hydrogen is neutral and 10(exp 3) where hydrogen is ionized. We hypothesize that all YSOs receive infall onto their outer disks which is steady (or slowly declining with time) and that FU Orionis outbursts are self-regulated, disk outbursts which occur only in systems which transport matter inward at a rate sufficiently high to cause hydrogen to be ionized in the inner disk. We estimate a critical mass flux of dm(sub crit)/dt = 5 x 10(exp 7) solar mass/yr independent of the magnitude of alpha for systems with one solar mass, three solar radius central objects. Infall accretion rates in the range of dm(sub in)/dt = 1-10) x 10(exp -6) solar mass/yr produce observed FU Orionis timescales consistent with estimates of spherical molecular cloud core collapse rates. Modeled ionization fronts are typically initiated near the inner edge of the disk and propogate out to a distance of several tens of stellar radii. Beyond this region, the disk transports mass steadily inward at the supplied constant infall rate. Mass flowing through the innermost disk annulus is equal to dm(sub in)/dt only in a time-averaged sense and is regulated by the ionization of hydrogen in the inner disk such that long intervals (approximately 1000 yr) of low-mass flux: (1-30) x 10(exp -8) solar mass/yr are punctuated by short intervals (approximately 100 yr) of high-mass flux: (1-30) x 10(exp -5) solar mass/yr. Timescales and mass fluxes derived for quiescent and outburst stages are consistent with estimates from observations of T Tauri and FU Orionis systems, respectively.

Bell, K. R.↗

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗