SEARCH · Engineering Papers

Results for “dependability modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Convectively Generated Gravity Waves In The Tropical Stratosphere: Case Studies And Importance For The Circulation Of The Middle Atmosphere

The advent of high altitude aircraft measurements in the stratosphere over tropical convective systems has made it possible to observe the mesoscale disturbances in the temperature field that these systems excite. Such measurements show that these disturbances have horizontal scales comparable to those of the underlying anvils (about 50-100 km) with peak to peak theta surface variations of about 300-400 meters. Moreover, correlative wind measurements from the tropical phase of the Stratosphere-Troposphere Exchange Project (STEP) clearly show that these disturbances are gravity waves. We present two case studies of anvil-scale gravity waves over convective systems. Using steady and time-dependent linear models of gravity wave propagation in the stratosphere, we show: (1) that the underlying convective systems are indeed the source of the observed phenomena; and (2) that their generating mechanism can be crudely represented as flow over a time-dependent mountain. We will then discuss the effects gravity waves of the observed amplitudes have on the circulation of the middle atmosphere, particularly the quasi-biennial, and semiannual oscillations.

Pfister, Leonhard↗

NASA NTRS2023

Effects of Optimisation Parameters on Data-Driven Magnetofrictional Modelling of Active Regions

Context . The solar magnetic field plays an essential role in the formation, evolution, and dynamics of large-scale eruptive structures in the corona. The estimation of the coronal magnetic field, the ultimate driver of space weather, particularly in the ‘low’ and ‘middle’ corona, is presently limited due to practical difficulties. Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing the eruption dynamics. We present a detailed study of the effects of optimisation of the non-inductive electric field on the TMFM of AR 12473. Aims . We aim to study the effects of varying the non-inductive electric field on the data-driven coronal simulations, for two alternative parametrisations. By varying parameters controlling the strength of the non-inductive electric field, we wish to explore the changes in flux rope formation and their early evolution and other parameters, for instance, axial flux and magnetic field magnitude. Methods . We used the high temporal and spatial resolution cadence vector magnetograms from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). The non-inductive electric field component in the photosphere is critical for energising and introducing twist to the coronal magnetic field, thereby allowing unstable configurations to be formed. We estimated this component using an approach based on optimising the injection of magnetic energy. Results . Our data show that flux ropes are formed in all of the simulations except for those with the lower values of these optimised parameters. However, the flux rope formation, evolution and eruption time varies depending on the values of the optimisation parameters. The flux rope is formed and has overall similar evolution and properties with a large range of non-inductive electric fields needed to determine the non-inductive electric field component that is critical for energising and introducing twist to the coronal magnetic field. Conclusions . This study shows that irrespective of non-inductive electric field values, flux ropes are formed and erupted, which indicates that data-driven TMFM can be used to estimate flux rope properties early in their evolution without needing to employ a lengthy optimisation process.

A. Kumari↗

NASA NTRS1990

Auroral modeling of the 3371 A emission rate - Dependence on characteristic electron energy

An efficient two-stream auroral electron model is used to study the deposition of auroral energy and the dependence of auroral emission rates on characteristic energy. This model incorporates the concept of average energy loss to reduce the computation time. This simple two-stream model produces integrated emission rates that are in excellent agreement with the much more complex multistream model of Strickland et al. (1983) but disagrees with a recent study by Rees and Lummerzheim (1989) that indicates that the N2 second positive emission rate is a strongly decreasing function of the characteristic energy of the precipitating flux. These calculations reveal that a 10 keV electron will undergo approximately 160 ionizing collisions, with an average energy loss per collision of 62 eV before thermalizing. The secondary electrons are created with an average energy of 42 eV. When all processes including the backscattered escape fluxes are taken into account, the average energy loss per electron-ion pair is 35 eV in good agreement with laboratory results.

Richards, P. G.↗

DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗

DOE OSTI2023

Simulations of biomass compression-screw feeding using a compressible non-Newtonian constitutive model

There is global interest in the conversion of biomass into sustainable low-carbon-footprint fuels and chemicals as an alternative to non-renewable fossil feedstocks. Feeding biomass solids into pressurized reactors is one of the key steps in biomass conversion. Predicting mechanical failure and energy requirements for this step helps avoid upstream processing bottlenecks and enables efficient operation of a biorefinery. Here, in this work, we developed a predictive computational model for biomass screw feeders that capture the highly viscous, non-Newtonian and compressible behavior of biomass slurries. Biomass compressible behavior is formulated by an equation of state and the non-Newtonian rheology is represented by a density-dependent viscosity model. Experimental data from two compression screw-feeder systems are presented as a validation for our model. Our model successfully predicted the location of the compressed biomass “plug”, biomass flow rate, and the required torque at different operating conditions for the experimental conditions studied in this work.

09 BIOMASS FUELS↗

NASA NTRS1993

Verification of a two-dimensional infiltration model for the resin transfer molding process

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin was verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTI IM7/8HS preform and corn oil. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle.

Hammond, Vincent H.↗

DOE OSTI2022

An energy-dependent electro-thermal response model of CUORE cryogenic calorimeter

The Cryogenic Underground Observatory for Rare Events (CUORE) is the most sensitive experiment searching for neutrinoless double-beta decay (0νββ) in 130 Te. CUORE uses a cryogenic array of 988 TeO 2 calorimeters operated at ~10 mK with a total mass of 741 kg. To further increase the sensitivity, the detector response must be well understood. Here, we present a non-linear thermal model for the CUORE experiment on a detector-by-detector basis. We have examined both equilibrium and dynamic electro-thermal models of detectors by numerically fitting non-linear differential equations to the detector data of a subset of CUORE channels which are well characterized and representative of all channels. Further, we demonstrate that the hot-electron effect and electric-field dependence of resistance in NTD-Ge thermistors alone are inadequate to describe our detectors' energy-dependent pulse shapes. We introduce an empirical second-order correction factor in the exponential temperature dependence of the thermistor, which produces excellent agreement with energy-dependent pulse shape data up to 6 MeV. We also present a noise analysis using the fitted thermal parameters and show that the intrinsic thermal noise is negligible compared to the observed noise for our detectors.

47 OTHER INSTRUMENTATION↗

NASA NTRS1981

A multiloop generalization of the circle criterion for stability margin analysis

In order to provide a theoretical tool well suited for use in characterizing the stability margins (e.g., gain and phase margins) of multiloop feedback systems, multiloop input-output stability results generalizing the circle stability criterion are considered. Generalized conic sectors with 'centers' and 'radii' determined by linear dynamical operators are employed to enable an engineer to specify the stability margins which he desires as a frequency-dependent convex set of modeling errors (including nonlinearities, gain variations, and phase variations) which the system must be able to tolerate in each feedback loop without instability. The resulting stability criterion gives sufficient conditions for closed-loop stability in the presence of such frequency-dependent modeling errors, even when the modeling errors occur simultaneously in all loops.

Safonov, M. G.↗

NASA NTRS1987

A theoretical study of the lifetime and transport of large ionospheric density structures

A three-dimensional time-dependent ionospheric model was used to study the spatial and temporal evolution and transport of large-scale high-density ionospheric structures for a range of solar cycle, seasonal, and IMF conditions. Both density depletions and enhancements were considered. It was found that, depending on the IMF, horizontal plasma convection can cause an initial structure to break up into multiple structures of various sizes, to become stretched into elongated segments, or to remain as a single distorted structure. The lifetime of an F-region density structure depends on several factors, including its magnitude, the initial location where it was formed, the season, the solar cycle, and the convection pattern. For example, in summer, the effects of a large density structure can disappear in a few hours or last as long as nine hours, while in winter the effects can persist for 24 hours. The passage of perturbed plasma flux tubes through sunlit and auroral regions can significantly increase the lifetime of plasma enhancements.

Schunk, R. W.↗

NASA NTRS2009

Probabilistic Risk Assessment for Decision Making During Spacecraft Operations

Decisions made during the operational phase of a space mission often have significant and immediate consequences. Without the explicit consideration of the risks involved and their representation in a solid model, it is very likely that these risks are not considered systematically in trade studies. Wrong decisions during the operational phase of a space mission can lead to immediate system failure whereas correct decisions can help recover the system even from faulty conditions. A problem of special interest is the determination of the system fault protection strategies upon the occurrence of faults within the system. Decisions regarding the fault protection strategy also heavily rely on a correct understanding of the state of the system and an integrated risk model that represents the various possible scenarios and their respective likelihoods. Probabilistic Risk Assessment (PRA) modeling is applicable to the full lifecycle of a space mission project, from concept development to preliminary design, detailed design, development and operations. The benefits and utilities of the model, however, depend on the phase of the mission for which it is used. This is because of the difference in the key strategic decisions that support each mission phase. The focus of this paper is on describing the particular methods used for PRA modeling during the operational phase of a spacecraft by gleaning insight from recently conducted case studies on two operational Mars orbiters. During operations, the key decisions relate to the commands sent to the spacecraft for any kind of diagnostics, anomaly resolution, trajectory changes, or planning. Often, faults and failures occur in the parts of the spacecraft but are contained or mitigated before they can cause serious damage. The failure behavior of the system during operations provides valuable data for updating and adjusting the related PRA models that are built primarily based on historical failure data. The PRA models, in turn, provide insight into the effect of various faults or failures on the risk and failure drivers of the system and the likelihood of possible end case scenarios, thereby facilitating the decision making process during operations. This paper describes the process of adjusting PRA models based on observed spacecraft data, on one hand, and utilizing the models for insight into the future system behavior on the other hand. While PRA models are typically used as a decision aid during the design phase of a space mission, we advocate adjusting them based on the observed behavior of the spacecraft and utilizing them for decision support during the operations phase.

dynamic fault trees↗

NASA NTRS2008

Radiation Protection Effectiveness of Polymeric Based Shielding Materials at Low Earth Orbit

Correlations of limited ionizing radiation measurements onboard the Space Transportation System (STS; shuttle) and the International Space Station (ISS) with numerical simulations of charged particle transport through spacecraft structure have indicated that usage of hydrogen rich polymeric materials improves the radiation shielding performance of space structures as compared to the traditionally used aluminum alloys. We discuss herein the radiation shielding correlations between measurements on board STS-81 (Atlantis, 1997) using four polyethylene (PE) spheres of varying radii, and STS-89 (Endeavour, 1998) using aluminum alloy spheres; with numerical simulations of charged particle transport using the Langley Research Center (LaRC)-developed High charge (Z) and Energy TRaNsport (HZETRN) algorithm. In the simulations, the Galactic Cosmic Ray (GCR) component of the ionizing radiation environment at Low Earth Orbit (LEO) covering ions in the 1< or equals Z< or equals 28 range is represented by O'Neill's (2004) model. To compute the transmission coefficient for GCR ions at LEO, O'Neill's model is coupled with the angular dependent LaRC cutoff model. The trapped protons/electrons component of LEO environment is represented by a LaRC-developed time dependent procedure which couples the AP8min/AP8max, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment resulting from interaction of GCR ions with upper atmosphere is modeled through extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. With the validity of numerical simulations through correlation with PE and aluminum spheres measurements established, we further present results from the expansion of the simulations through the selection of high hydrogen content commercially available polymeric constituents such as PE foam core and Spectra fiber(Registered TradeMark) composite face sheet to assess their radiation shield properties as compared to generic PE.

Badavi, Francis F.↗

DOE OSTI2023

Insights into the Oxygen Evolution Reaction on Graphene-Based Single-Atom Catalysts from First-Principles-Informed Microkinetic Modeling

Single-atom transition metals embedded in nitrogen-doped graphene have emerged as promising electrocatalysts due to their high activity and low material cost. These materials have been shown to catalyze a variety of electrochemical reactions, but their active sites under reaction conditions remain poorly understood. Using first-principles density functional theory calculations, we develop a pH-dependent microkinetic model to evaluate the relative performance of transition metal catalysts embedded in fourfold N-substituted double carbon vacancies in graphene for the oxygen evolution reaction. We find that reaction pathways involving intermediates co-adsorbed on the metal site are preferred on all transition metals. These pathways lead to enhancements in catalytic activity and broaden the activity peak when compared with purely thermodynamics-based predictions. Furthermore, these findings demonstrate the importance of investigating reaction pathways on graphene-based catalysts and other two-dimensional (2D) materials that involve metal active centers decorated by spectator intermediate species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1986

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

NASA NTRS1990

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

NASA NTRS1998

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

NASA NTRS1989

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

NASA NTRSDate not supplied

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↗

NASA NTRS1980

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