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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

F region above Kauai - Measurement, model, modification

Results of the Lagopedo II experiment conducted from Kauai, Hawaii to investigate the ionospheric modification that occurs when rocket combustion products are introduced into the O(+)-rich F region are presented. The experiment involved the detonation of a chemical explosion in the F2 peak accompanied by rocket-borne measurements of ion composition and electron content in the vicinity of the explosion. The experimental data is found to be in good agreement with the predictions of a model of the nighttime ion densities in the midlatitude laminar ionosphere, with the exception of N2(+) densities before the explosion. H2O(+) and H3O(+) currents produced by considerable H2O outgassing from the rocket are used to determine a H3O(+)/H2O(+) dissociative recombination rate averaging 1.6 to 1.08, depending on model assumptions. At the time of the explosion, an ionic void 1 km in radius is observed, the boundary of which is characterized by a steep gradient in ionic densities. Evidence of variations in the concentrations of ambient ion species, new reactant species and ionic depletion by sweeping is also obtained.

Johnson, C. Y.↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
DOE OSTI2021

Direct driving of electronic and phononic degrees of freedom in a honeycomb bilayer with infrared light

Here, we study theoretically AB-stacked honeycomb bilayers driven by light in resonance with an infrared phonon within a tight-binding description. We characterize the phonon properties of honeycomb bilayers with group theory and construct an electronic time-dependent tight-binding model for the system following photoexcitation in resonance with an infrared phonon. We adopt an atomically adiabatic approximation [V. Mohanty and E. J. Heller, Proc. Natl. Acad. Sci. 116, 18316 (2019)] to describe classically vibrating nuclei, but obtain the Floquet quasienergy spectrum associated with the time-dependent model exactly. We introduce a general scheme to disentangle the complex low-frequency Floquet spectrum to elucidate the relevant Floquet bands. As a prototypical example, we consider bilayer graphene. We find that light in the low-frequency regime can induce a band gap in the quasienergy spectrum in the vicinity of the $K$ points even if it is linearly polarized, in contrast with the expectations within the Born-Oppenheimer approximation and the high-frequency regime. Finally, we analyze the diabaticity of the driven electron and driven phonon processes and find contrasting effects on the autocorrelation functions at the same driving frequency: Driven phonons preserve the character of the initial state while driven electrons exhibit strong deviations within a few drive cycles. The procedure outlined here can be applied to other materials to describe the combined effects of low-frequency light on phonons and electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Modeling of high-latitude currents in a substorm

The currents and electric fields in the area covered by the poleward set of field-aligned currents are computed by means of a quantitative model, in which the high-latitude band considered carries most of the poleward electrojet, and lies poleward of the region covering the inner magnetosphere and corresponding ionosphere. The Birkeland current and aurorally enhanced conductivity are assumed to be uniformly distributed across the band. The time-dependent conductivity model used is based on electron fluxes and mean energies measured from the S3-2 satellite. Joule heating of the upper atmosphere was found to be about 2 x 10 to the 11th W during the substorm period. Model values for the strength of the electrojet and the amount of Joule heating agree to within about 20% with results based on a simple Cowling conductivity band.

Karty, J. L.↗
DOE OSTI2022

Uncertainty reduction in residual stress measurements by an optimised inverse solution using nonconsecutive polynomials

Many destructive methods for measuring residual stresses such as the slitting method require an inverse analysis to solve the problem. The accuracy of the result as well as an uncertainty component (the model uncertainty) depends on the basis functions used in the inverse solution. The use of a series expansion as the basis functions for the inverse solution was analysed in a previous work for the particular case where functions orders grew consecutively. The present work presents a new estimation of the model uncertainty and a new improved methodology to select the final basis functions for the case where the basis is composed of polynomials. Including nonconsecutive polynomial orders in the basis generates a larger space of possible solutions to be evaluated and allows the possibility to include higher-order polynomials. The paper includes a comparison with two other inverse analyses methodologies applied to synthetically generated data. With the new methodology, the final error is reduced and the uncertainty estimation improved.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hybrid model to evaluate the frequency-dependent leakage inductance of partially-filled transformers

The leakage inductance of a transformer designed for a power electronic converter can drop significantly as the switching frequency is increased due to skin and proximity effects. Although the magnetic image method-based double-2-D model can predict the low-frequency leakage inductance of a partially-filled transformer with sufficient accuracy, it is inherently a frequency-independent model. While Dowell’s 1-D model uses frequency-dependent relations to account for both skin and proximity effects, its accuracy is severely affected by the assumed winding geometry. In this paper, a semianalytical hybrid model is proposed that uses superposition to combine a modified Dowell’s model with the double-2-D model to predict the true leakage inductance of partially-filled transformers at any given frequency. All three conductor types—round, foil, and litz wire—are modeled and analyzed. The quasi-2-D model is further investigated on a variable inductance transformer (VIT) whose winding geometry can be modified mechanically to vary its leakage inductance. With less than 5 % error throughout, the semi-analytically evaluated leakage inductances are in excellent agreement with the finite element method (FEM) simulated and experimentally measured leakage inductances.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
NASA NTRS2000

Evaluation of the Regional Atmospheric Modeling System in the Eastern Range Dispersion Assessment System

The Applied Meteorology Unit is conducting an evaluation of the Regional Atmospheric Modeling System (RAMS) contained within the Eastern Range Dispersion Assessment System (ERDAS). ERDAS provides emergency response guidance for operations at the Cape Canaveral Air Force Station and the Kennedy Space Center in the event of an accidental hazardous material release or aborted vehicle launch. The prognostic data from RAMS is available to ERDAS for display and is used to initialize the 45th Range Safety (45 SW/SE) dispersion model. Thus, the accuracy of the 45 SW/SE dispersion model is dependent upon the accuracy of RAMS forecasts. The RAMS evaluation task consists of an objective and subjective component for the Florida warm and cool seasons of 1999-2000. The objective evaluation includes gridded and point error statistics at surface and upper-level observational sites, a comparison of the model errors to a coarser grid configuration of RAMS, and a benchmark of RAMS against the widely accepted Eta model. The warm-season subjective evaluation involves a verification of the onset and movement of the Florida east coast sea breeze and RAMS forecast precipitation. This interim report provides a summary of the RAMS objective and subjective evaluation for the 1999 Florida warm season only.

Case, Jonathan↗
NASA NTRS2005

Small-Scale Drop-Size Variability: Empirical Models for Drop-Size-Dependent Clustering in Clouds

By analyzing aircraft measurements of individual drop sizes in clouds, it has been shown in a companion paper that the probability of finding a drop of radius r at a linear scale l decreases as l(sup D(r)), where 0 less than or equals D(r) less than or equals 1. This paper shows striking examples of the spatial distribution of large cloud drops using models that simulate the observed power laws. In contrast to currently used models that assume homogeneity and a Poisson distribution of cloud drops, these models illustrate strong drop clustering, especially with larger drops. The degree of clustering is determined by the observed exponents D(r). The strong clustering of large drops arises naturally from the observed power-law statistics. This clustering has vital consequences for rain physics, including how fast rain can form. For radiative transfer theory, clustering of large drops enhances their impact on the cloud optical path. The clustering phenomenon also helps explain why remotely sensed cloud drop size is generally larger than that measured in situ.

Marshak, Alexander↗
NASA NTRS1989

A magnetohydrodynamical model of the inner coma of Comet Halley

The magnetometer onboard the Giotto spacecraft observed a diamagnetic cavity surrounding the nucleus of Comet Halley. The location of the boundary of this diamagnetic cavity is determined by a balance between an inward magnetic pressure gradient force and an outward ion-neutral drag force, associated with collisions between the outwardly flowing neutrals and the stagnated ions. A one-dimensional time-dependent MHD model has been developed for the inner coma of Comet Halley, and includes ion-neutral collisions, photochemical production and loss of plasma, and finite conductivity. This model is used to investigate the plasma dynamics in the vicinity of the diamagnetic cavity boundary surface. A narrow transition layer with enhanced plasma density is shown to exist just inside the boundary, although a full understanding of this layer will require a two- or three-dimensional MHD model. The flux of cometary ions into this shocklike layer is removed by electron-ion recombination. The thickness of this layer is determined by the Mach number of the incident flow.

Cravens, T. E.↗
NASA NTRS1998

The Impact of FU Orionis Outbursts and the Solar Nebula

Protostellar systems are variable on many timescales. One of the most dramatic forms of variability known to occur in low mass stellar systems is the FU Orionis outburst (Herbig 1977). Throughout a typical outburst lasting several decades, system luminosities may be a hundred times what is typical of the quiesent state. FU Orionis outburst events are thought to have significant impact on the thermal structure of the protosolar nebula; their existence has been used to explain features in the meteoritic record from thermally induced homogenization to the formation of chondrules. Until recently, the magnitude of the likely effect from such outbursts has been largely speculative due to the lack of a detailed understanding of the outburst mechanism. Recent numerical models (Bell\& Lin 1994) have demonstrated the viability of the observational hypothesis (Hartmann\& Kenyon 1985) that the radiation observed during outburst is emitted by a luminous circumstellar disk transporting mass at a thousand times the quiesent rate. Light curves and color and line width evolution observed in FU Orionis systems are naturally explained by time dependent outbursting model disks (Bell et al. 1995). The radial temperature structure and shape of the disk during outburst derived from these models may be used to calculate the outburst's expected impact on primitive material at various radii throughout the disk. In this review, we will begin by discussing what is known about the FU Orionis outburst phenomenon from recent observations and theory including statistically deduced outburst timescales and observed peak temperatures. Unless covered by another author, we will discuss the evidence which suggests that outburst radiation is emitted by a circumstellar disk rather than by the star and will briefly review the thermal instability as a mechanism for outburst. We will then report on recent work which investigates the likely heating of solar nebula material due to FU Orionis outbursts including the following effects: (1) heating of the planet forming region by direct radiation from the hot inner nebula; (2) heating by the diffuse radiation field of a coccooning envelope; and (3) time-dependent penetration of the increased luminosity from the above sources into the optically thick nebula. Some of this work is currently in progress. The potential effects on condensation and migration in the nebula and the thermal processing of solids will be evaluated.

Bell, Robbins↗
DOE OSTI2021

The Significance of Aggregation Methods in Functional Group Modeling

The growth of forests and the feedbacks between forests and environmental changes are central issues in the planetary carbon cycle, global climate change, and basic plant ecology. A challenge to understanding both growth and feedbacks from local to global scales is that many critical metabolic processes vary among species. An innovation in solving this challenge is the recognition that species can be lumped into “functional groups” based on metabolic similarity, and these functional groups can then be studied in computational models that simulate ecosystem function. Despite the vast resources devoted to functional group studies and the progress made by them, an important logical and biological question has not been formally addressed, “How do the groupings alter the results of modeling studies?” To what extent do modeling results depend on the choices made in aggregating taxa into functional groups. Here, we consider the effects of using different aggregation strategies in simulating the carbon dynamics of a deciduous forest. Understanding the impacts that aggregation strategy has on efforts to simulate regional-to-global-scale forest dynamics offers insights into both ecosystem regulation and model function and addresses this central problem in the study of carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Ionospheric response to the sustained high geomagnetic activity during the March '89 great storm

A simulation was conducted to model the high-latitude ionospheric to the sustaied level of high geomagnetic activity for the great magnetic storm period of March 13-14, 1989. The geomagnetic and solar activity indices and the Defense Meterological Satellite Program (DMSP) F8 and F9 satellite data for particle precipitation and high-latitude convection were used as inputs to a time-dependent ionospheric model (TDIM). The results of the TDIM were compared to both DMSP plasma density data and ground-based total electron content (TEC) measurements for the great storm period as well as with earlier storm observations. The comparisons show that the overall structure of the high-latitude ionosphere was dominated by an increased convection speed within the polar cap that led to increased ion temperatures. In turn, this enhanced the NO(+) density, raised the atomic-to-molecular ion transition height to over 300 km, decreased N(sub m)F(sub 2), increased h(sub m)F(sub 2), and in places either increased n(sub e) at 800 km or slightly decreased it. The morphology of the ionosphere under these extreme conditions was considerably different than that modeled for less distributed intervals. These differences included the character of the dayside tongue of ionization that no longer extended deep into the polar cap. Instead, as a result of the ion heating and consequent reduction in N(sub m)F(sub 2), a large polar hole occupied much of the polar region. This polar hole extended beyond the auroral oval and merged with the night sector midatitude trough. The limitaions associated with the applicability of the TDIM to the geomagnetic conditions present on March 13 and 14 are discussed. The primary limitations of the TDIM derive from the limited temporal resolution of the model input parameters and the lack of suitably dynamic thermospheric specification for the great storm conditions. These limitations leads to midlatitude ionospheric storm phases that do no follow those observed.

Sojka, J. J.↗
NASA NTRS1994

Patches in the polar ionosphere: UT and seasonal dependence

The seasonal and Universal Time (UT) dependencies of patches in the polar ionosphere are simulated using the Utah State University time dependent ionospheric model (TDIM). Patch formation is achieved by changing the plasma convection pattern in response to temporal changes in the interplanetary magnetic field (IMF) B(sub y) component during periods of southwrd IMF. This mechanism redirects the plasma flow from the dayside high-density region, which is the source of the tongue of ionization (TOI) density feature, through the throat and leads to patches, rather than a continuous TOI. The model predicts that the patches are absent at winter solstice (northern hemisphere) between 0800 and 1200 UT and that they have their largest seasonal intensity at winter solstice between 2000 and 2400 UT. Between winter solstice and equinox, patches are strong and present all day. Patches are present in summer as well, although their intensity is only tens of percent above the background density. These winter-to-equinox findings are also shown to be consistent with observations. The model was also used to predict times at which patch observations could be performed to determine the contributions from other patch mechanisms. This observational window is +/- 20 days about winter solstice between 0800 and 1200 UT in the northern hemisphere. In this observational window the TOI is either absent or reduced to a very low density. Hence the time dependent electric field mechanism considered in this study does not produce patches, and if they are observed, then they must be due to some other mechanism.

Sojka, J. J.↗
DOE OSTI2025

Systematic uncertainties in the measurement of the neutron lifetime using the Lunar Prospector neutron spectrometer

The lifetime of free neutrons measured in the laboratory has a longstanding disparity of ≈9 s. A space-based technique has recently been proposed to independently measure the neutron lifetime using interactions between the galactic cosmic rays and a low atmosphere planetary body. This technique has not produced competitive results yet due to constraints of nonoptimized data that contain large systematic errors. We use data from the neutron spectrometer on-board NASA's Lunar Prospector, and study two large systematics in the measurement of neutron lifetime: the lunar subsurface temperature and the lunar surface composition. We use the HeCd and HeSn neutron spectrometer data when the spacecraft was in a highly elliptical orbit during the orbit insertion period. We report the neutron lifetime using four different models that each have different choices of surface temperature and composition. The 5° [Prettyman et al ., J. Geophys. Res.: Planets 111, 2005JE002656 (2006)] and 2° [Wilson et al., Phys. Rev. C 104, 045501 (2021)] rebinned maps result in 777.6±11.7 s and 739.6±10.8 s, respectively. For the 20° map (Prettyman et al., 2006), constant equatorial and a latitude-dependent temperature model result in 738.6±10.8 s and 767.3±11.2 s, respectively. Increasing the complexities of the models accounting for the systematic effects increase the measured lifetime. However, the reported measurements are not competitive with the laboratory results due to large unaccounted systematics resulting from nonoptimized measurements and modeling assumptions. This work serves as a study of systematic uncertainties for future neutron lifetime measurements using the space-based technique. We estimate the effect on the lifetime from the choice of temperature model to be to be 28.7 ±15.5 s, and choice of compositional map (for 20° and 5° maps) to be 10.3 ±12.2 s.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Breakdown of vertically propagating two-dimensional gravity waves forced by orography

The propagation of orographic gravity waves into an atmosphere with exponentially decreasing density is simulated with a two-dimensional, nonlinear, time-dependent numerical model. After the stationary wave is established over the mountain, the model predicts that wave breaking causes a large reduction of the vertical momentum flux in the flow, not only at levels where wave breaking is present, but also far below the lowest occurrence of overturning. More than half of the decrease in momentum flux is explained by the presence of large amplitude, downward propagating waves, which are generated in regions of wave breaking. The downward propagating waves appear almost simultaneously with overturning, and have nonzero phase speeds, suggesting a strongly nonlinear generation mechanism that depends on local wave properties. The generation of these downward propagating waves is a robust process, insensitive to mountain height, mountain width, or density scale height. These results have important implications for observational studies of orographically generated waves as well as for schemes that seek to parameterize the effects of orogoraphy in large-scale models.

Bacmeister, Julio T.↗
DOE OSTI2020

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1983

Ship waves and lee waves

Dynamics analogous to those of surface ship waves on water of finite depth are noted for the three-dimensional trapped lee wave modes produced by an isolated obstacle in a stratified fluid. This vertical trapping of wave energy is modeled by uniform upstream flow and stratification, bounded above by a rigid lid, and by a semiinfinite fluid of uniform stability whose wind velocity increases exponentially with height, representing the atmosphere. While formal asymptotic solutions are produced, limited quantitative usefulness is obtained through them because of the limitations of the approximations and the infinity of modes in the solution. Time-dependent numerical models are accordingly developed for both surface ship waves and internal and atmospheric ship waves, yielding a variety of results.

Sharman, R. D.↗