Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗
DOE OSTI2025

Machine Learning based Correlation of the Mechanical Properties of Sub-sized and Standard-sized Specimens

Mechanical testing with sub-sized specimens is essential in the nuclear industry, offering the ability to conduct tests in confined spaces with lower irradiation and expediting material qualification. However, smaller specimens exhibit different material behavior across scales, a phenomenon known as the "specimen size effect". In this study, we compiled over 1,000 tensile testing records, covering 54 parameters such as material type, composition, manufacturing details, irradiation conditions, specimen dimensions, and tensile properties through a comprehensive literature review. We focus on correlating sub-sized and standard specimens’ tensile mechanical properties on SS316 alloy, which has the most extensive dataset available. We explore ML-based models and uncertainty quantification for tensile properties, analyze key factors influencing these properties, and compare the effectiveness of ML models with existing analytical methods in addressing the specimen size effect.

tensile properties↗
DOE OSTI2020

Evolving Trust in Robots: Specification Through Sequential and Comparative Meta-Analyses

Objective The objectives of this meta-analysis are to explore the presently available empirical findings on the antecedents of trust in robots and use this information to expand upon a previous meta-analytic review of the area. Background Human–robot interaction (HRI) represents an increasingly important dimension of our everyday existence. Currently, the most important element of these interactions is proposed to be whether the human trusts the robot or not. We have identified three overarching categories that exert effects on the expression of trust. These consist of factors associated with (a) the human, (b) the robot, and (c) the context in which any specific HRI event occurs. Method The current body of literature was examined and all qualifying articles pertaining to trust in robots were included in the meta-analysis. A previous meta-analysis on HRI trust was used as the basis for this extended, updated, and evolving analysis. Results Multiple additional factors, which have now been demonstrated to significantly influence trust, were identified. The present results, expressed as points of difference and points of commonality between the current and previous analyses, are identified, explained, and cast in the setting of the emerging wave of HRI. Conclusion The present meta-analysis expands upon previous work and validates the overarching categories of trust antecedent (human-related, robot-related, and contextual), as well as identifying the significant individual precursors to trust within each category. A new and updated model of these complex interactions is offered. Application The identified trust factors can be used in order to promote appropriate levels of trust in robots.

Behavioral Sciences↗
DOE OSTI2020

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimization of Lithium Content and Sintering Additives in Tape Cast Lithium Garnet Electrolyte Sheets

The lithium garnet (LLZO) remains one of the most promising electrolyte materials for producing solid state batteries. Tape casting is one method for producing the thin ceramic electrolyte layer required for a high-performance solid state battery. Processing of LLZO is complicated by several factors including lithium volatilization, abnormal grain growth and phase instability. Many of these issues are exacerbated in tape cast sheets by the sheets' higher surface area to volume ratio, when compared to pellet processing methods. Here, we report on an environmentally friendly aqueous tape casting process for LLZO using methylcellulose as a binder and compare final sheet properties to other, solvent-based, options. Solids loading, binder content, dispersant composition, and wetting agent type are varied to improve ease of casting and enhance properties of the as-cast tape. Lithium content in sintered LLZO sheets is controlled by the inclusion of excess Li 2 CO 3 as lithium source in the ceramic tape itself and by a novel method of producing a lithium-saturated furnace atmosphere. Sintering time and temperature are optimized, and rapid thermal processing is investigated as a method of reducing lithium loss. Al 2 O 3 and MgO are examined for use as sintering additives. Characterization of LLZO sheets is performed with XRD, EIS and SEM. Li 2 CO 3 and sintering additive levels are optimized for ionic conductivity and density. The optimization effort improved the conductivity and density of 100µm thick sintered sheets to greater than 3x10 -4 S/cm conductivity and 90% density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2023

Walking on the dark side: Anthropogenic factors limit suitable habitat for gray wolf ( Canis lupus ) in a large natural area covering Belarus and Ukraine

Due to successful conservation initiatives and legislations, the grey wolf (Canis lupus) is recolonising its historic range in Europe. However, wolves have never been extirpated across large areas in Eastern Europe but are often constrained to remote and inaccessible places due to centuries of persecution. This study aimed to identify the potentially suitable wolf habitats in Polesia, a massive cross-border lowland region extending over southern Belarus and northern Ukraine, which are often neglected in large carnivore studies at the continental scale. We hypothesized that anthropogenic rather than environmental factors govern wolf habitat suitability. We used a dataset of 4191 GPS locations obtained from radio-collared wolves (n = 26) and confirmed observations (n = 231) during 2014–2021 and applied maximum entropy method to estimate relative habitat suitability for wolves in Polesia. Artificial light at night (ALAN), proportion of cropland and tree cover were the most important factors affecting wolf habitat suitability. Road densities contributed poorly to predicting habitat suitability for wolves. Our models predicted a quarter of Polesia as suitable habitat and revealed priority areas connecting the important source populations in the Chornobyl Exclusion Zone in the east and the Bialowieza Forest in the west and thus essential for long-term wolf conservation. Our results provide the bases for effective, long-term wolf monitoring and management programs in both Belarus and Ukraine. However, national and transboundary wolf management in Polesia has been extremely challenging since 2022 due to the ongoing war and subsequent habitat degradation in this part of Europe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Uranium particle age dating, aggregation, and model age best estimators

We present important aspects of uranium particle age dating by Large-Geometry Secondary Ion Mass Spectrometry (LG-SIMS) that can introduce bias and increase model age uncertainties, especially for small, young, and/or low-enriched particles. This metrology is important for applications related to International Nuclear Safeguards. We explore influential factors related to model age estimation, including the effects of evolving surface chemistry on inter-element measurements of particles (e.g., Th and U), detector background, and aggregation methods using simulated and actual particle samples. We introduce a new model age estimator, called “mid68”, that supplements 95% confidence intervals, providing a “best estimate” and uncertainty about the most likely age. The mid68 estimator can be calculated using the Feldman and Cousins method or Bayesian methods and provides a value with a symmetric uncertainty that can be used for calculations and approximate aggregation of processed model age values when the raw data and correction factors are not available. For particles yielding low 230 Th counts amidst nonzero detector background, their underlying model age probability distributions are asymmetric, so the mid68 estimator provides additional robust information regarding the underlying model age likelihood. This study provides a comprehensive and timely examination of critical aspects of uranium particle age dating as more laboratories establish particle chronometry capabilities.

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

Predicting Thermoset Resin Cure Kinetics for Composite Processing

Predicting a material’s thermal properties is not only useful for designing parts but also for defining processing parameters. This is especially applicable to wet-filament winding operations and the production of composite parts using thermoset res ins. For these, the formation of cross-linking bonds is affected by time and temperature, but the degree of cross-linking can be modeled using phenomenological methods. In the current study, the Kissinger and Flynn/Wall/Ozawa methods were both used to calculate kinetic parameters and produced values that were in agreement with each other. The kinetic parameters activation energy ( E a ) and pre-exponential factor (𝐴) calculated from exothermic peak data increased between an epoxy and resole phenol resin system, with a cyanate ester having the highest. The successful application of both methods suggests applicability to future studies for developing cure recipes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Case Study of ISO/TS 6336-22 Micropitting Method: Preprint

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears — Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for a method to calculate the risk of micropitting in gear sets. It assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. A sliding parameter is used to adjust the film thickness in the contact zone to the region of highest sliding, which is where micropitting has been observed. The permissible specific film thickness is best determined through experience or testing, but there is an option to estimate it based on the lubricant’s failure load stage in FZG testing. This paper is a case study of how the method behaves when applied to several gear sets and applications that have experienced micropitting in the field.

gears↗
DOE OSTI2020

Case Study of ISO/TS 6336-22 Micropitting Calculation

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears - Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for a method to calculate the risk of micropitting in gear sets. It assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. A sliding parameter is used to adjust the film thickness in the contact zone to the region of highest sliding, which is where micropitting has been observed. The permissible specific film thickness is best determined through experience or testing, but there is an option to estimate it based on the lubricant’s failure load stage in FZG testing. This paper is a case study of how the method behaves when applied to several gear sets and applications that have experienced micropitting in the field.

ENGINEERING,WIND ENERGY↗
DOE OSTI2022

Multi-tool optimization for computer controlled optical surfacing

With the rapid development of precision technologies, the demand of high-precision optical surfaces has drastically increased. These optical surfaces are mainly fabricated with computer controlled optical surfacing (CCOS). In a CCOS process, a target surface removal profile is achieved by scheduling the dwell time for a set of machine tools. The optimized dwell time should be positive and smooth to ensure convergence to the target while considering CNC dynamics. The total run time of each machine tool is also expected to be balanced to improve the overall processing efficiency. In the past few decades, dwell time optimization for a single machine tool has been extensively developed. While the methods are applicable to multi-tool scenarios, they fail to consider the overall contributions of multiple tools simultaneously. In this paper, we conduct a systematic study on the strategies for multi-tool dwell time optimization and propose an innovative method for simultaneously scheduling dwell time for multiple tools for the first time. First, the influential factors to the positiveness and smoothness of dwell time solutions for a single machine tool are analyzed. The compensation strategies that minimize the residual while considering the CNC dynamics limit are then proposed. Afterwards, these strategies are extended to the proposed multi-tool optimization that further balances the run time of machine tools. Finally, the superiority of each strategy is carefully studied via simulation and experiment. The experiment is performed by bonnet polishing a 60 mm × 60 mm mirror with three tools of different diameters ( i.e., 12 mm, 8 mm, and 5 mm). The figure error of the mirror is reduced from 45.42 nm to 11.18 nm root mean square in 13.28 min. Moreover, the measured polishing result well coincides with the estimation, which proves the effectiveness of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparative Study of Wind Energy Potential Estimation Methods for Wind Sites in Togo and Benin (West Sub-Saharan Africa)

The characterization of wind speed distribution and the optimal assessment of wind energy potential are critical factors in selecting a suitable site for wind power plants (WPP). The Weibull distribution law has been used extensively to analyze the wind characteristics of candidate WPP sites, and to estimate the available and deliverable energy. This paper presents a comparative study of five wind energy resource assessment methods as they applied to the context of wind sites in West Sub-Saharan Africa. We investigated three numerical approaches, namely, the adaptive neuro-fuzzy inference system (ANFIS), the multilayer perceptron method (MLP), and support vector regression (SVR), to derive the distribution law of wind speeds and to optimally quantify the corresponding wind energy potential. Next, we compared these three approaches to two well-known Weibull distribution law-based methods: the empirical method of Justus (EMJ) and the maximum likelihood method (MLM). Case study results indicated that the neural network-based methods, ANFIS and MLP, yielded the most accurate distribution fits and wind energy potential estimates, and consequently, are the most recommended methods for the wind sites in Togo and Benin. The orders of magnitude of the root mean squared error (RMSE) in estimating the recoverable energy using ANFIS were, respectively, 10-4 and 10-5 for Lomé and Cotonou, while MLP achieved an RMSE order of magnitude of 10-3 for both sites.

17 WIND ENERGY↗
DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detectability of Varied Hybridization Scenarios Using Genome-Scale Hybrid Detection Methods

Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This phenomenon has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, as well as summary methods that predict individual hybridization events from a subset of taxa. However, a lack of empirical comparisons between methods – especially those pertaining to large networks with varied hybridization scenarios – hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson’s D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and HyDe, Patterson’s D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events between sets of three taxa. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow (γ). We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher relative false negative rates across all methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving extinct or unsampled ghost lineages, HyDe is the only method able to identify the direction of hybridization, distinguishing the source parental lineages from recipient hybrid lineages. Lastly, we test the methods on a dataset of ultraconserved elements from the bee subfamily Nomiinae, finding possible hybridization events between clades which correspond to regions of poor support in the species tree estimated in a previous study.

Bjorner, Marianne B.↗
DOE OSTI2022

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗