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,261 records · Page 70

Assessment of Physics Models for Phase Transition Kinetics

The time dependence of phase diagrams and how to model rate dependent transitions remains one of the key unanswered questions in physics. When a material is loaded dynamically through equilibrium phase boundaries, it is the kinetics that determines the real time expression of a phase transition. Here we report the atomic and nanosecond-scale quantification of kinetics of shock-driven phase transition in multiple materials. We uniquely make use of a both a simple shock as well as shock-and-hold loading pathways compress different crystalline solids and induce structural phase transitions below melt. Coupling shock loading with time-resolved synchrotron x-ray diffraction (DXRD), we probe the structural transformations of these solids in the short-lived high pressure and temperature states generated. The novelty and power of using DXRD for the assessment of kinetics of phase transitions lies in the ability to discover and identify new phases and to examine kinetics without prior knowledge of a material's phase diagram. Our results provide a quantified expression and a physics model of kinetics of formation of high-pressure phases under shock loading: transition incubation time, evolution, completion time and crystallization rate.

36 MATERIALS SCIENCE↗

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analysis of MIR-18 results for physical and biological dosimetry: radiation shielding effectiveness in LEO

We compare models of radiation transport and biological response to physical and biological dosimetry results from astronauts on the Mir space station. Transport models are shown to be in good agreement with physical measurements and indicate that the ratio of equivalent dose from the Galactic Cosmic Rays (GCR) to protons is about 3/2:1 and that this ratio will increase for exposures to internal organs. Two biological response models are used to compare to the Mir biodosimetry for chromosome aberration in lymphocyte cells; a track-structure model and the linear-quadratic model with linear energy transfer (LET) dependent weighting coefficients. These models are fit to in vitro data for aberration formation in human lymphocytes by photons and charged particles. Both models are found to be in reasonable agreement with data for aberrations in lymphocytes of Mir crew members: however there are differences between the use of LET dependent weighting factors and track structure models for assigning radiation quality factors. The major difference in the models is the increased effectiveness predicted by the track model for low charge and energy ions with LET near 10 keV/micrometers. The results of our calculations indicate that aluminum shielding, although providing important mitigation of the effects of trapped radiation, provides no protective effect from the galactic cosmic rays (GCR) in low-earth orbit (LEO) using either equivalent dose or the number of chromosome aberrations as a measure until about 100 g/cm 2 of material is used.

Flight Experiment↗

Monte Carlo Simulation of THz Multipliers

Schottky Barrier diode frequency multipliers are critical components in submillimeter and Thz space based earth observation systems. As the operating frequency of these multipliers has increased, the agreement between design predictions and experimental results has become poorer. The multiplier design is usually based on a nonlinear model using a form of harmonic balance and a model for the Schottky barrier diode. Conventional voltage dependent lumped element models do a poor job of predicting THz frequency performance. This paper will describe a large signal Monte Carlo simulation of Schottky barrier multipliers. The simulation is a time dependent particle field Monte Carlo simulation with ohmic and Schottky barrier boundary conditions included that has been combined with a fixed point solution for the nonlinear circuit interaction. The results in the paper will point out some important time constants in varactor operation and will describe the effects of current saturation and nonlinear resistances on multiplier operation.

East, J.↗

Site‐selective fluorescence and spectroscopic properties of Yb‐doped lanthanum titanate glasses

Abstract Ytterbium‐doped lanthanum titanate glasses were prepared by levitation melting for the detailed characterization of the spectroscopic properties in the rare‐earth titanate glass host. Low‐temperature fluorescence spectroscopy reveals distinct site‐selectivity in both static and lifetime fluorescence measurements suggesting an absence of clustering as well as significant variation of local ytterbium environments. Typical site‐selectivity behavior of a shrinking Stark manifold with lower excitation energy is observed. At 77 K, both the mean emission frequency and the fluorescence lifetime initially increase as the excitation energy decreases from about 11100 to 10750 and then slightly decrease at lower excitation energy. Temperature‐dependent lifetime measurements between 77 and 420 K show a decreasing lifetime with increasing temperature and are well described by a two‐level thermal activation model. The temperature‐dependent fluorescence spectroscopy coupled with a room temperature white light absorption measurement allow the determination of the Stark energy levels of in lanthanum titanate glass as well as the calculation of the laser cross‐sections.

Materials Science↗

Incorporating concentration-dependent sediment microbial activity into methylmercury production kinetics modeling

We report in anoxic environments, anaerobic microorganisms carrying the hgcAB gene cluster can mediate the transformation of inorganic mercury (Hg(II)) to monomethylmercury (MMHg). The kinetics of Hg(II) transformation to MMHg in periphyton from East Fork Poplar Creek (EFPC) in Oak Ridge, TN have previously been modeled using a transient availability model (TAM). The TAM for Hg(II) methylation combines methylation/demethylation kinetics with kinetic expressions for processes that decrease Hg(II) and MMHg availability for methylation and demethylation (multisite sorption of Hg(II) and MMHg, Hg(II) reduction/Hg(0) oxidation). In this study, the TAM is used for the first time to describe MMHg production in sediment. We assessed MMHg production in sediment microcosms using two different sediment types from EFPC: a relatively anoxic, carbon-rich sediment with higher microbial activity (higher CO 2 production from sediment) and a relatively oxic, sandy, carbon-poor sediment with lower microbial activity (lower CO 2 production from sediment). Based on 16s rRNA sequencing, the overall microbial community structure in the two sediments was retained during the incubations. However, the hgcA containing methanogenic Euryarchaeota communities differed between sediment types and their growth followed different trajectories over the course of incubations, potentially contributing to the distinct patterns of MMHg production observed. The general TAM paradigm performed well in describing MMHg production in the sediments. However, the MMHg production and ancillary data suggested the need to revise the model structure to incorporate terms for concentration-dependent microbial activity over the course of the incubations. We modified the TAM to include Monod-type kinetics for methylation and demethylation and observed an improved fit for the carbon-rich, microbially active sediment. Overall our work shows that the TAM can be applied to describe Hg(II) methylation in sediments and that including expressions accounting for concentration-dependent microbial activity can improve the accuracy of the model description of the data in some cases.

54 ENVIRONMENTAL SCIENCES↗

Mean field model parameterization to recover coverage-dependent kinetics

Lateral interactions between adsorbates introduce coverage dependence into adsorption energies and activation barriers of surface reactions. Lattice-based kinetic Monte Carlo (kMC) simulations can capture these interactions quantitatively but are laborious to parameterize and solve. Mean field models are more tractable, but protocols to construct and parameterize them are unclear. Here we explore the ability of a coverage-aware mean-field model to map to a lattice-kMC model of a generic two-step reaction network, including quasi-equilibrated adsorption and rate-limiting dissociation steps. We derive expressions for mean-field and coverage-dependent adsorption energies and dissociation barriers and parameterize against lattice-kMC predictions. We show that the parameterized mean-field rates correlate with ground truth lattice-kMC results across a wide range of reaction conditions and identify regions where the mean field fails. The mean field model similarly captures kMC-derived rate-order, Arrhenius and Sabatier plots at a greatly reduced computational cost. Further, the results provide guidance for parameterizing mean-field models, benchmarked against explicit lattice-based approaches for incorporating the influence of coverage effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search for electroweak production of charginos and neutralinos in proton-proton collisions at $$ \sqrt{s} $$ = 13 TeV

A direct search for electroweak production of charginos and neutralinos is presented. Events with three or four leptons, with up to two hadronically decaying $\tau$ leptons, or two same-sign light leptons are analyzed. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions with a center of mass energy of 13 TeV, recorded with the CMS detector at the LHC. The results are interpreted in terms of several simplified models. These represent a broad range of production and decay scenarios for charginos and neutralinos. A parametric neural network is used to target several of the models with large backgrounds. In addition, results using orthogonal search regions are provided for all the models, simplifying alternative theoretical interpretations of the results. Depending on the model hypotheses, charginos and neutralinos with masses up to values between 300 and 1450 GeV are excluded at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Demonstrating the impact of ro-vibrationally excited H 2 on divertor detachment via population modelling

EDGE2D-EIRENE profiles characteristic for low recycling, high recycling and detached plasmas are postprocessed applying a vibrationally resolved collisional radiative (CR) model based on the Yacora solver to predict excited state densities and thus the molecular Lyman-, Werner- and Fulcher band emission intensities. These emission intensities are compared to predictions from the AMJUEL database standardly used in EIRENE. Lower Fulcher band and higher Lyman- and Werner band emission is predicted by the Yacora model in comparison to AMJUEL suggesting that inconsistent results from model based evaluation of measured spectra are expected depending on which model and molecular emission band is considered. The implications on detachment rate predictions are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Interstellar sodium lines and the two-component model of galactic H I regions.

The intensities and widths of the interstellar sodium lines are considered as possible discriminants between two different two-component models recently proposed for galactic H I regions. Depending on the model, the observed sodium lines are either formed entirely in the clouds or, alternatively, in both the clouds and the intercloud medium. Present observational evidence is insufficient to distinguish between these two possibilities.

Hobbs, L. M.↗

Cirrus clouds. I - A cirrus cloud model. II - Numerical experiments on the formation and maintenance of cirrus

A simplified cirrus cloud model is presented which may be used to investigate the role of various physical processes in the life cycle of a cirrus cloud. The model is a two-dimensional, time-dependent, Eulerian numerical model where the focus is on cloud-scale processes. Parametrizations are developed to account for phase changes of water, radiative processes, and the effects of microphysical structure on the vertical flux of ice water. The results of a simulation of a thin cirrostratus cloud are given. The results of numerical experiments performed with the model are described in order to demonstrate the important role of cloud-scale processes in determining the cloud properties maintained in response to larger scale forcing. The effects of microphysical composition and radiative processes are considered, as well as their interaction with thermodynamic and dynamic processes within the cloud. It is shown that cirrus clouds operate in an entirely different manner than liquid phase stratiform clouds.

Starr, D. OC.↗

Validation of clear-sky fluxes for tropical oceans from the Earth Radiation Budget Experiment

The existence and magnitude of a systematic bias in the clear-sky longwave fluxes from the Earth Radiation Budget Experiment (ERBE) is investigated. The bias is apparently introduced because the ERBE method for scene identification does not account for large zonal gradients in longwave absorption by water vapor. The ERBE fluxes are compared to fluxes calculated with a radiative transfer model from ship radiosonde measurements. The comparison is based upon an analysis of 5 yr of coincident satellite and radiosonde observations for equatorial ocean regions. The differences between the ERBE and model fluxes are examined as functions of sea surface temperature (SST) and relative humidity. The authors use height-mean relative humidity bar-RH as an index of atmospheric moisture. The average offset between and model ERBE fluxes ranges between +2 and +6 W/sq m for SSTs above 295 K, and the gradients with respect to SST are nearly identical. However, the difference between the model and ERBE depends significantly on the tropospheric relative humidity. ERBE fluxes exceed model fluxes for bar-RH above 70%, and the maximum offset of +9 to +12 W/sq m is consistent with previous estimates. There are also indications that the clear-sky fluxes for bar-RH below 25% may be underestimated by about 10-15 W/sq m. Since extreme values of height-mean humidity are relatively infrequent, the net bias introduced in the ERBE monthly mean clear-sky fluxes is generally less than the systematic error in estimates introduced in the ERBE monthly mean clear-sky fluxes is generally less than the systematic error in estimates of the instantaneous fluxes. These findings support earlier work on the coupling between SST and the atmospheric greenhouse effect, in particular the existence of a super greenhouse effect for oceans warmer than 300 K. Recent reports of much larger systematic differences are not supported by this analysis. The results indicate that comparison of general circulation model (GCM) and ERBE clear-sky longwave fluxes will depend explicitly on atmospheric humidity.

Collins, W. D.↗

CARES/Life Ceramics Durability Evaluation Software Enhanced for Cyclic Fatigue

The CARES/Life computer program predicts the probability of a monolithic ceramic component's failure as a function of time in service. The program has many features and options for materials evaluation and component design. It couples commercial finite element programs--which resolve a component's temperature and stress distribution--to reliability evaluation and fracture mechanics routines for modeling strength-limiting defects. The capability, flexibility, and uniqueness of CARES/Life have attracted many users representing a broad range of interests and has resulted in numerous awards for technological achievements and technology transfer. Recent work with CARES/Life was directed at enhancing the program s capabilities with regards to cyclic fatigue. Only in the last few years have ceramics been recognized to be susceptible to enhanced degradation from cyclic loading. To account for cyclic loads, researchers at the NASA Lewis Research Center developed a crack growth model that combines the Power Law (time-dependent) and the Walker Law (cycle-dependent) crack growth models. This combined model has the characteristics of Power Law behavior (decreased damage) at high R ratios (minimum load/maximum load) and of Walker law behavior (increased damage) at low R ratios. In addition, a parameter estimation methodology for constant-amplitude, steady-state cyclic fatigue experiments was developed using nonlinear least squares and a modified Levenberg-Marquardt algorithm. This methodology is used to give best estimates of parameter values from cyclic fatigue specimen rupture data (usually tensile or flexure bar specimens) for a relatively small number of specimens. Methodology to account for runout data (unfailed specimens over the duration of the experiment) was also included.

Nemeth, Noel N.↗

Application of a Temperature-Dependent Load Prediction Method to a RUAG Six-Component Block-Type Balance

Temperature-dependent data of a RUAG six-component block-type balance was analyzed to assess the accuracy of two load prediction methods for temperature-dependent balance data. The supplied data was prepared for the analysis by splitting it into calibration and check load data subsets. The first calibration data subset was obtained at a temperature of 294 Kelvin. The second calibration data subset was obtained at a temperature of 315 Kelvin. A subset of 38 points was extracted from the second data set and used as check loads so that the accuracy of the two load prediction methods could be tested. First, the Iterative Method in combination with an extended independent and dependent variable set was used for the balance load prediction. This approach fits electrical outputs as a function of loads and the temperature and, afterwards, constructs a load iteration scheme from the regression coefficients so that loads can be predicted from outputs and the temperature during a wind tunnel test. The Non-Iterative Method was also used for the load prediction. This alternate method can more easily be implemented in a data system as loads are directly fitted as a function of electrical outputs and the temperature. Analysis results for the axial force are only discussed in the paper as similar results were obtained for the other five load components. Results for both methods clearly show that the cross-product term constructed from either a primary gage load or a primary gage output and the temperature explains the majority of the temperature-dependent part of the predicted balance load. This term models the temperature dependent nature of the gage sensitivity. Therefore, it is recommended to apply primary gage loadings at different temperatures during a balance calibration whenever temperature effects need to be described. These loadings will contain information about the temperature-dependent nature of the gage sensitivities that can be quantified by related cross-product terms in regression models of the data.

strain-gage balance↗

Impulsive H-alpha diagnostics of electron-beam-heated solar flare model chromospheres

Time-dependent H-alpha profiles were computed for the dynamic model atmospheres of Fisher, Canfield, and McClymont (1985) simulating the effects of an intense impulsively initiated power-law beam of electrons incident on the chromosphere. The temporal response of H-alpha arises from three separate physical mechanisms, whose relative importance varies over the line profile. The fastest variations (typically less than 0.1 s for the explosive heating discussed here) arise from energy imbalance; these are apparent on chromospheric heating and cooling time scales and have their greatest amplitude at line center. Slower variations arise from ionization imbalance and are most apparent in the blue wing. The slowest variations arise from hydrodynamic effects and are related to the formation of a chromospheric condensation; these are most apparent in the red wing. These results provide a basis for the design and analysis of observations of H-alpha, in coordination with hard X-rays or microwaves, to test mechanisms of energy transport in solar flares.

Canfield, Richard C.↗

Position-dependent matter density probability distribution function

Here, we introduce the position-dependent probability distribution function (PDF) of the smoothed matter field as a cosmological observable. In comparison to the PDF itself, the spatial variation of the position-dependent PDF is simpler to model and has distinct dependence on cosmological parameters. We demonstrate that the position-dependent PDF is characterized by variations in the local mean density, and we compute the linear response of the PDF to the local density using separate universe N-body simulations. The linear response of the PDF to the local density field can be thought of as the linear bias of regions of the matter field selected based on density. We provide a model for the linear response, which accurately predicts our simulation measurements. We also validate our results and test the separate universe consistency relation for the local PDF using global universe simulations. We find excellent agreement between the two, and we demonstrate that the separate universe method gives a lower variance determination of the linear response.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗