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 685 records · Page 38

Self-consistent modelling of the polar thermosphere and ionosphere to magnetospheric convection and precipitation (invited review)

It has recently been demonstrated that the dramatic effects of plasma precipitation and convection on the composition and dynamics of the polar thermosphere and ionosphere include a number of strong interactive, or feedback, processes. To aid the evaluation of these feedback processes, a joint three dimensional time dependent global model of the Earth's thermosphere and ionosphere was developed in a collaboration between University College London and Sheffield University. This model includes self consistent coupling between the thermosphere and the ionosphere in the polar regions. Some of the major features in the polar ionosphere, which the initial simulations indicate are due to the strong coupling of ions and neutrals in the presence of strong electric fields and energetic electron precipitation are reviewed. The model is also able to simulate seasonal and Universal time variations in the polar thermosphere and ionospheric regions which are due to the variations of solar photoionization in specific geomagnetic regions such as the cusp and polar cap.

Rees, D.↗

Discrete event simulation tool for analysis of qualitative models of continuous processing systems

An artificial intelligence design and qualitative modeling tool is disclosed for creating computer models and simulating continuous activities, functions, and/or behavior using developed discrete event techniques. Conveniently, the tool is organized in four modules: library design module, model construction module, simulation module, and experimentation and analysis. The library design module supports the building of library knowledge including component classes and elements pertinent to a particular domain of continuous activities, functions, and behavior being modeled. The continuous behavior is defined discretely with respect to invocation statements, effect statements, and time delays. The functionality of the components is defined in terms of variable cluster instances, independent processes, and modes, further defined in terms of mode transition processes and mode dependent processes. Model construction utilizes the hierarchy of libraries and connects them with appropriate relations. The simulation executes a specialized initialization routine and executes events in a manner that includes selective inherency of characteristics through a time and event schema until the event queue in the simulator is emptied. The experimentation and analysis module supports analysis through the generation of appropriate log files and graphics developments and includes the ability of log file comparisons.

Malin, Jane T.↗

Accelerated Strength Testing of Thermoplastic Composites

Constant ramp strength tests on unidirectional thermoplastic composite specimens oriented in the 90 deg. direction were conducted at constant temperatures ranging from 149 C to 232 C. Ramp rates spanning 5 orders of magnitude were tested so that failures occurred in the range from 0.5 sec. to 24 hrs. (0.5 to 100,000 MPa/sec). Below 204 C, time-temperature superposition held allowing strength at longer times to be estimated from strength tests at shorter times but higher temperatures. The data indicated that a 50% drop in strength might be expected for this material when the test time is increased by 9 orders of magnitude. The shift factors derived from compliance data applied well to the strength results. To explain the link between compliance and strength, a viscoelastic fracture model was investigated. The model, which used compliance as input, was found to fit the strength data only if the critical fracture energy was allowed to vary with temperature reduced stress rate. This variation in the critical parameter severely limits its use in developing a robust time-dependent strength model. The significance of this research is therefore seen as providing both the indication that a more versatile acceleration method for strength can be developed and the evidence that such a method is needed.

Reeder, J. R.↗

Modelling ionospheric density structures

Large-scale density structures are a common feature in the high-latitude ionsphere. The structures were observed in the dayside cusp, polar cap, and nocturnal auroral region over a range of altitudes, including the E-region, F-region and topside ionosphere. The origins, lifetimes and transport characteristics of large-scale density structures were studied with the aid of a three-dimensional, time-dependent ionospheric model. Blob creation due to particle precipitation, the effect that structured electric fields have on the ionosphere, and the lifetimes and transport characteristics of density structures for different seasonal, solar cycle, and interplanetary magnetic field (IMF) conditions were studied. The main conclusions drawn are: (1) the observed precipitation energy fluxes are sufficient for blob creation if the plasma is exposed to the precipitation for 5 to 10 minutes; (2) structured electric fields produce structured electron densities, ion temperatures, and ion composition; (3) the lifetime of an F-region density structure depends on several factors, including the initial location where it was formed, the magnitude of the perturbation, season, solar cycle and IMF; and (4) depending on the IMF, horizontal plasma convection can cause an initial structure to break up into multiple structures of various sizes, remain as a single distorted structure, or become stretched into elongated segments.

Schunk, R. W.↗

The Land Surface in Current and Planned MERRA Reanalysis Products

Current global atmospheric reanalysis products such as the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5), the NASA Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2), and the Japanese Reanalysis for Three Quarters of a Century (JRA-3Q) provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff, that are widely used in research and applications. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture, soil temperature, snow conditions, and screen-level air temperature and humidity from satellite observations and in situ measurements. In this presentation, we review the land surface modeling and data assimilation components of the suite of current and planned MERRA reanalysis products. In addition to MERRA-2, we will discuss the latest NASA reanalysis, MERRA for the 21st century (M21C), which is currently under production, as well as the development and planning of the next version of the MERRA reanalysis, tentatively labeled MERRA-3. In MERRA-2, observations-based precipitation data products are used to correct the precipitation falling on the land surface. Outside of the high-latitudes and Africa, the daily, 0.5-degree, gauge-based Climate Prediction Center (CPC) Unified (CPCU) product is used. In Africa, the pentad, 2.5-degree, satellite- and gauge-based CPC Merged Analysis of Precipitation (CMAP) product is used. Poleward of 62.5 degrees latitude, the land surface sees the precipitation generated by the atmospheric model in the cycling data assimilation system. This configuration provides improved soil moisture estimates compared to those of the original (version 1) MERRA estimates, which did not benefit from the use of precipitation observations. Moreover, the use of precipitation observations facilitates a seamless spin-up of the land surface initial conditions across the MERRA-2 production streams. The use of a gauge-only precipitation product in MERRA-2 across much of the globe, however, adversely impacts the quality of the MERRA-2 land surface estimates in regions with poor gauge coverage, including most of South America and Australia. Therefore, the forthcoming M21C reanalysis uses satellite- and gauge-based precipitation from the Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement Mission (IMERG). This change results in significant improvements in the quality of the M21C soil moisture estimates in the Southern Hemisphere compared to those from MERRA-2. Planning for MERRA-3 focuses on the assimilation of soil moisture observations from the Soil Moisture Active Passive (SMAP) mission and the Advanced Scatterometer (ASCAT), along with snow cover area fraction observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to further improve the quality of the land surface estimates from the reanalysis. As a first step towards the assimilation of land surface observations in MERRA-3, the offline (land-only) M21C-Land reanalysis is currently under development as a supplemental M21C product that includes the assimilation of SMAP, ASCAT, and MODIS observations. Preliminary results from M21C and M21C-Land will be discussed in the context of MERRA-2 and plans for MERRA-3.

Rolf Reichle↗

Comparison of the Nimbus-4 BUV ozone data with the Ames two-dimensional model

A comparison is made of the first two years of Nimbus 4 backscattered ultraviolet (BUV) ozone measurements with the predictions of the Ames two dimensional model. The ozone observations used consist of the mixing ratio on the 1, 2, 5, and 10 mb pressure surfaces. The data are zone and time averaged to obtain seasonal means for 1970 and 1971 and are found to show strong and repeatable meridional and seasonal dependencies. The model used for comparison with the observations extends from 80 N to 80 S latitude and from altitudes of 0 to 60 km with 5 deg horizontal grid spacing and 2.5 km vertical grid spacing. Chemical reaction and photolysis rates are diurnally averaged and the photodissociation rates are corrected for the effects of scattering. The large altitude, latitude, and seasonal changes in the ozone data agree with the model predictions. Model predictions of the sensitivity of the comparisons to changes in the assumed mixing ratios of water vapor, odd nitrogen, and odd chlorine, as well as to changes in the ambient temperature and transport parameters are also shown.

Borucki, W. J.↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗