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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The effect of metastable atomic levels in neutral tungsten for erosion and impurity transport applications

Long-lived metastable states in neutral tungsten are shown to potentially impact measurements of tungsten erosion from plasma facing components. Time-dependent collisional radiative modeling of neutral tungsten is used to analyze the role of these states in tungsten emission and ionization. The large number of non-quasistatic atomic states in neutral tungsten can take on the order of milliseconds to reach equilibrium, depending on plasma conditions, causing erosion measurements to be affected by the metastable populations. Previous measurements using the 400.88 nm tungsten emission line could be affected by these non-quasistatic metastable effects. Therefore, a scheme for measuring the relative metastable fractions is proposed through simultaneous observation of multiple ultraviolet spectral lines of neutral W. The accuracy of gross erosion measurements could potentially be increased by inclusion of these previously unconsidered metastable effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Explainable Machine Learning for Functional Data

Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).

Elastic Shape Analysis↗
DOE OSTI2023

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

SPT-3G D1: A Measurement of Secondary Cosmic Microwave Background Anisotropy Power

We report new measurements of millimeter-wave temperature power spectra in the angular multipole range $1700 \le \ell \le 11,000$ (wavelengths $13^\prime \gtrsim λ\gtrsim 2^\prime$). We use two years of data in three observing bands centered near 95, 150, and 220 GHz from the SPT-3G receiver on the South Pole Telescope that cover a 1646 deg$^2$ region of the Southern sky. Using the measured power spectra, we present constraints on the thermal and kinematic Sunyaev-Zel'dovich (SZ) effects, radio galaxies, and cosmic infrared background (CIB). We find that inferred SZ powers are dependent on the detailed modeling of the thermal SZ-CIB correlation, and to a lesser extent on the assumed angular dependence of the SZ spectra. We report constraints for simulation-based model templates as well as fits where the angular dependencies of the SZ and CIB power spectra are allowed to vary. In the latter case at $\ell=3000$, we find thermal SZ power at 143 GHz of $D_{3000}^{\rm tSZ} = 4.91\pm0.37\, μ{\rm K}^2$ and kinematic SZ power of $D_{3000}^{\rm kSZ} =1.75\pm0.86\, μ{\rm K}^2$. We use the measured kinematic SZ power to estimate the duration of reionization, noting that the reionization inferences are sensitive to the model choices and assumed level of homogeneous kinematic SZ power from the late-time universe. We find a 95% limit on the duration from an ionization fraction of 25% to 75% of $Δ^{50} z_{\rm re} <\,3.8$ based on a semi-analytic model, or a limit on the duration from an ionization fraction of 5% to 95% of $Δ^{90} z_{\rm re} <\,6.1$ based on the AMBER simulations.

Chaubal, P. [Melbourne U.]↗
DOE OSTI2020

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Analytical models of the strength and ductility of CNT reinforced metal matrix nano composites under elevated temperatures

Carbon nanotubes (CNTs) can greatly enhance the strength of metal matrix composites while resulting in ductility loss. This strength-ductility tradeoff dilemma always confines the development of material design and real-life applications. Here in this work, a temperature dependent strengthening analytical model is proposed for CNT reinforced metal matrix composites which considers three common strengthening mechanisms: Orowan looping effect, thermal expansion mismatch effect, and load bearing effect. The proposed model can predict composite material strength with different volume fractions of CNTs and under different temperatures. Combining the strengthening model with the stress based modified Mohr-Coulomb (sMMC) ductile fracture model, a ductility analytical model is then derived. This ductility analytical model includes the influences of temperature, multi-axial stress loading conditions, as well as the aforementioned three strengthening mechanisms. A good agreement has been achieved between literature published experimental data and analytical predictions for both composite material strength and ductility loss. The proposed two analytical models can provide a straightforward way to study the strength-ductility relationship for CNT reinforced metal matrix composites over a wide range of temperatures and different stress states, and then provide guidance on new material design and processing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗
DOE OSTI2022

Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution

The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The mechanism for acetate formation in electrochemical CO (2) reduction on Cu: selectivity with potential, pH, and nanostructuring

Nanostructured Cu catalysts have increased the selectivities and geometric activities for high value C–C coupled (C 2 ) products (ethylene, ethanol, and acetate) in the electrochemical CO (2) reduction reaction (CO (2) RR). The selectivity among the high-value C 2 products is also altered, where for instance the yield of acetate increases with alkalinity and is dependent on the catalyst morphology. The reaction mechanisms behind the selectivity towards acetate vs. other C 2 products remain controversial. In this work, we elucidate the reaction mechanism for acetate formation by using ab initio simulations, a coupled kinetic-transport model, and loading dependent experiments. We find that trends in acetate selectivity can be rationalized from variations in electrolyte pH and the local mass transport properties of the catalyst and not from changes in Cu's intrinsic activity. The selectivity mechanism originates from the transport of ketene, a stable (closed shell) intermediate, away from the catalyst surface into solution where it reacts to form acetate. While this type of mechanism has not yet been discussed in the CO (2) RR, variants of it may explain similar selectivity fluctuations observed for other stable intermediates like CO and acetaldehyde. Our proposed mechanism suggests that acetate selectivity increases with increasing pH, decreasing catalyst roughness and significantly varies with the applied potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)

The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Application of physiologically based pharmacokinetic modeling for sertraline dosing recommendations in pregnancy

Pregnancy is a period of significant change that impacts physiological and metabolic status leading to alterations in the disposition of drugs. Uncertainty in drug dosing in pregnancy can lead to suboptimal therapy, which can contribute to disease exacerbation. A few studies show there are increased dosing requirements for antidepressants in late pregnancy; however, the quantitative data to guide dose adjustments are sparse. We aimed to develop a physiologically based pharmacokinetic (PBPK) model that allows gestational-age dependent prediction of sertraline dosing in pregnancy. A minimal physiological model with defined gut, liver, plasma, and lumped placental-fetal compartments was constructed using the ordinary differential equation solver package, ‘mrgsolve’, in R. We extracted data from the literature to parameterize the model, including sertraline physicochemical properties, in vitro metabolism studies, disposition in nonpregnant women, and physiological changes during pregnancy. The model predicted the pharmacokinetic parameters from a clinical study with eight subjects for the second trimester and six subjects for the third trimester. Based on the model, gestational-dependent changes in physiology and metabolism account for increased clearance of sertraline (up to 143% at 40 weeks gestational age), potentially leading to under-dosing of pregnant women when nonpregnancy doses are used. The PBPK model was converted to a prototype web-based interactive dosing tool to demonstrate how the output of a PBPK model may translate into optimal sertraline dosing in pregnancy. Quantitative prediction of drug exposure using PBPK modeling in pregnancy will support clinically appropriate dosing and increase the therapeutic benefit for pregnant women.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Determination of energy-dependent neutron backgrounds using shadow bars

Understanding the neutron background is essential for determining the neutron yield from nuclear reactions. Here, the neutron backgrounds were determined for heavy-ion collision experiment using the shadow-bar method, where beams of 40,48 Ca at 56, 140 MeV/u impinged on targets of 58,64 Ni and 112,124 Sn. In the analysis presented here, brass shadow bars are placed in front of organic liquid scintillator neutron detectors to determine the energy-dependent neutron background fractions. The measurement of neutron spectra with and without shadow bars is important to determine the neutron background more accurately. The neutron background, along with its sources and systematic uncertainties, are explored with a focus on the impact of background models and their dependence on neutron energy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗