SEARCH · Search NASA

Results for “Reproducibility”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A crystal plasticity finite element model embedding strain-rate sensitivities inherent to deformation mechanisms: Application to alloy AZ31

We report the fundamental power-law relationship representing the flow rule in crystal visco-plasticity ensures uniqueness in the selection of slip systems accommodating imposed plastic strain-rates. The power-law relationship also introduces an artificially high strain-rate sensitivity in crystal plasticity simulations, unless a high value of the power-law exponent is used. However, the use of high values for the exponent is limited by numerical tractability. This paper presents a numerical method implemented in a crystal plasticity finite element (CPFE) model for embedding any value of the power-law exponent reflecting the true material strain-rate sensitivity. Importantly, the method does not increase computation time involved in the simulations. The enhanced CPFE model is used to interpret and predict a complex strain-rate sensitive response and microstructural evolution of AZ31 Mg alloy. Measured values of strain-rate sensitivity for slip and twinning modes are used in the simulations. Calculations show that the model successfully captures the phenomena pertaining to the effect of changing applied strain-rate on the mechanical response including flow stress and evolution of texture and twinning for a broad range of strain-rates ranging from 10 -3 s -1 to 10 3 s -1 and loading orientations in tension and compression. It is shown that such predictions are a consequence of not only relative amounts of slip and twinning activities driven by a set of accurately characterized hardening law parameters but also values of the strain-rate sensitivities inherent to individual deformation mechanisms. Besides, the model validates the measured strain-rate dependency of deformation mechanisms while accurately reproducing the mechanical data. Hence, the model can be used to verify and further refine or infer measured strain-rate sensitivity per deformation mechanism by reproducing experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

DOE OSTI2021

Prostate Bed Delineation Guidelines for Postoperative Radiation Therapy: On Behalf Of The Francophone Group of Urological Radiation Therapy

Prostate bed (PB) irradiation is considered the standard postoperative treatment after radical prostatectomy (RP) for tumors with high-risk features or persistent prostate-specific antigen, or for salvage treatment in case of biological relapse. Four consensus guidelines have been published to standardize practices and reduce the interobserver variability in PB delineation but with discordant recommendations. To improve the reproducibility in the PB delineation, the Francophone Group of Urological Radiotherapy (Groupe Francophone de Radiothérapie Urologique [GFRU]) worked to propose a new and more reproducible consensus guideline for PB clinical target volume (CTV) definition.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2022

A scalable and tunable thermoreversible polymer for 3D human pluripotent stem cell biomanufacturing

Human pluripotent stem cells (hPSCs) are an exciting and promising source to enable cell replacement therapies for a variety of unmet medical needs. Though hPSCs can be successfully derived into numerous physiologically relevant cell types, effective translation to the clinic is limited by challenges in scalable production of high-quality cells, cellular immaturity following the differentiation process, and the use of animal-derived components in culture. To address these limitations, we have developed a fully defined, reproducible, and tunable thermoreversible polymer for high-quality, scalable 3D cell production. Our reproducible synthesis method enables precise control of gelation temperature (24°C–32°C), hydrogel stiffness (100–4000 Pa), and the prevention of any unintended covalent crosslinking. After material optimization, we demonstrated hPSC expansion, pluripotency maintenance, and differentiation into numerous lineages within the hydrogel. Overall, this 3D thermoreversible hydrogel platform has broad applications in scalable, high-quality cell production to overcome the biomanufacturing burden of stem cell therapy.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Investigations of vertical wind variations at a mountain top in the Himalaya using Doppler Lidar observations and model simulations

Hight-resolved observations of vertical winds remain nearly non-existing over the Himalayas, despite of anticipated crucial role of vertical motions in transporting pollution across the Himalayan hms. The present study analyze the vertical wind observations from surface to 1 km above ground level over Manora Peak (29.4° N; 79.5° E; 1958 m ams!) in the Himalaya performed using a Doppler Lidar during the Ganges Valley Aerosol Experiment (GVAX). Vertical wind exhibited a pronounced diurnal variability at Manora Peak comprising of upward motions during the daytime (05-10 UT) and downward motions during nighttime typical of a mountainvalley system. Mean vertical wind speeds are observed to be varying from -0.8 to +0.8 ms- 1 during the study period with a variance of 0.1-1.5 m2s·2, which is attributed to the thermally driven turbulence. Mean vertical winds are observed to be stronger in the Doppler Lidar profiles above Manora Peak (-0_8 to 0.8 ms- 1) as compared to near surface measurements at this station using an ultrasonic anemometer ( -0.4 to 0.4 ms - 1 ), and low altitude stations in India. Daytime vertical wind speeds are observed to be higher during pre-monsoon (0.81 ms- 1), as compared to post-monsoon (0.24 ms- 1) and winter (0.33 ms- 1)_ Average Black Carbon (BC) concentrations are significantly higher during strong upward vertical winds, which indicates efficient transport of polluted air mass from low-altitude regions to the Himalaya. Weather Research and Forecasting (WRF) model reproduces the observed diurnal pattern in the vertical wind at the observation site however the model underestimates the variability.

Shukla, K K.↗

DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Magnetothermal properties of Ho 1-x Dy x Al 2 (x = 0, 0.05, 0.10, 0.15, 0.25 and 0.50) compounds

Magnetic and magnetocaloric properties of H o 1 - x D y x A l 2 compounds with x = 0 , 0.05 , 0.10 , 0.15 , 0.25 and 0.50 , modelled using a Hamiltonian that includes the exchange interactions between Ho-Dy, Ho-Ho and Dy-Dy ions in addition to the crystalline electric field and the Zeeman effects, have been compared with those determined experimentally. In order to reproduce experimentally observed global ferromagnetic ordering temperatures and spin reorientation transition temperatures as x Dy varies, the exchange interactions between Ho-Dy and Ho-Ho were set as free parameters and adjusted to match the experimental results. We demonstrate that heat capacity of polycrystalline materials in non-zero magnetic fields can be satisfactory reproduced by using the average of multiple magnetic field directions with respect to the crystallographic coordinate system, while reasonably good agreement between experimentally determined and theoretically predicted magnetocaloric effects can be achieved considering an average of only three field directions.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Effect of interatomic potential on the energetics of hydrogen and helium-vacancy complexes in bulk, or near surfaces of tungsten

Hydrogen (H) trapping by helium-vacancy (He$-$V) complexes in bulk and the near surface region of tungsten (W) have been investigated by molecular statics calculations that evaluate two different W$-$H interatomic potentials, which use the same W$-$He, He$-$He and He$-$H potentials. One of the W$-$H potentials is a bond-order potential (BOP) developed by Juslin et al., while the other is an embedding atom method (EAM) potential developed by Wang et al. Both potentials overestimate the H binding energies to He clusters in bulk W, as compared to DFT calculations, but properly predict the functional form of the H binding energies to He clusters with increasing number of He and H. The BOP simulations reveal that H binding energies to He x V complexes generally increase with increasing number of He. However, the EAM results indicate that the H binding energy as a function of number of He depends on the number of H, and the H binding energies change slightly at high He content. Compared with available DFT data, both BOP and EAM underestimate the H binding energies to He x V 2 H m complexes. Here, the BOP reproduces the He formation energy below a W surface, while the EAM potential better reproduces the H formation energy and the interactions between H and He$-$V complexes. Based on these comparisons, we determine that the EAM potential is more accurate than BOP for large-scale molecular dynamics simulations of W$-$He$-$H interactions. The EAM potential predicts that the difference in the average binding energies of H to stable He$-$V complexes near the W surface is less than 0.2 eV and the difference decreases with increasing He content. Thus, the EAM potential indicates that the effect of surfaces on H binding energies to large He$-$V complexes below the W surfaces can be ignored.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗

DOE OSTI2018

Synapse: Synthetic application profiler and emulator

Abstract Motivated by the need to emulate workload execution characteristics on high- performance and distributed heterogeneous resources, we introduce Synapse. Synapse is used as a proxy application (or "representative application") for real workloads, with the advantage that it can be tuned in different ways and dimensions, and also at levels of granularity that are not possible with real applications. Synapse has a platform-independent application profiler, and has the ability to emulate profiled workloads on a variety of resources. Experiments show that the automated profiling performed using Synapse captures an application's characteristics with high fidelity. The emulation of an application using Synapse can reproduce the application's execution behavior in the original run- time environment, and can also reproduce those behaviors on different run-time environments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Evaluation of electronic state of Cs-adsorbed clay minerals by NEXAFS analysis using DFT calculations

In this work, a combined analytical method of NEXAFS measurement and DFT-calculation was employed for the evaluation of Cs states in clay minerals. The measured Cs M 4,5 NEXAFS spectra of Cs halides were analyzed using the DFT-calculations in order to establish the analytical methods. The NEXAFS spectra of the Cs halides were well reproduced by incorporating the core-hole strength. The Cs M 4,5 NEXAFS spectrum of the clay minerals was well reproduced by the DFT-calculations including the major transitions and tail structures with the established method. Further evaluation of this spectrum by charge density analysis suggested that these major transitions and the tail structures likely reflect the bonding state and the local environment around the Cs atoms. Comparison of electronic states of Cs in the clay mineral with those in the Cs halides by DFT-calculations has shown that the interaction between Cs and the nearest-neighbor atom is largest in the clay mineral, which differs from those of Cs halides. This happened due to the energy level of Cs-5s and 5p, which was closer to that of O-2s and 2p than the s and p orbitals of other alkali metal and alkali earth metal elements.

36 MATERIALS SCIENCE↗

DOE OSTI2019

TENDL: Complete Nuclear Data Library for Innovative Nuclear Science and Technology

The TENDL library is now established as one of the major nuclear data libraries in the world, striving for completeness and quality of nuclear data files for all isotopes, evaluation methods, processing and applied performance. To reach this status, some basic principles have been applied which sets it apart from other libraries: reproducible dedicated evaluations when differential data are available, through determination of nuclear models implemented in TALYS and their parameters, completeness (with or without experimental data), format and processing standardization, automation of production and reproducibility. In this paper, we will outline how such an approach has become a reality, and recall some of the past successes since the first TENDL release in 2008. Next, we will demonstrate the performance of the latest TENDL releases for different application fields, as well as new approaches for uncertainty quantification based on Bayesian inference methods and possible differential and integral adjustments. Also, current limitations of the library performances due to modelling and needs for new and more precise experimental data will be outlined.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

An Automated Ultraclean Ion Exchange Separation Method for the Determinations of 232Th and 238U in Copper using Inductively Coupled Plasma Mass Spectrometry

This work presents a novel automated analytical method developed for high throughput ultrasensitive determinations of Th and U in copper using inductively coupled plasma mass spectrometry (ICP-MS). The method is based on the use of ultra clean sample preparation procedures including the use of a fully automated off-line chromatography system for the extraction and pre-concentration of the analytes prior to ICP-MS analysis. The separation system is equipped with a single reusable chromatographic column that, with the rigorously clean procedures developed herein, provided a low carryover, low background, and reproducible automated separation method. Isotope dilution methods were used for quantitation. Method Detection Limits (MDLs) of 3.7 and 9.4 fg·g-1 were obtained for the quantitation of 232Th and 238U in copper, respectively, corresponding to activities of 0.0148 and 0.116 microBq·kg-1 for 232Th and 238U, respectively. The analytical method provided high and reproducible tracer recoveries. The use of an automated system significantly reduced the extremely tedious active work time required for the chemist (by ca. 80%) relative to manually performing the separation. This method could be adapted to other critical sample matrices requiring the utmost in high throughput clean chemistry procedures for automated ultrasensitive analyses.

Arnquist, Isaac J.↗

DOE OSTI2021

UCGretina geant4 simulation of the GRETINA Gamma-Ray Energy Tracking Array

UCGretina, a geant4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton continuum. Both of these features are needed in order to accurately extract gamma-ray yields from spectra collected in in-beam gamma-ray spectroscopy measurements with beams traveling at v / c ≳0.3 at the National Superconducting Cyclotron Laboratory and the Facility for Rare Isotope Beams. In the process of developing the model, we determined that millimeter-scale layers of passive germanium surrounding the active volumes of the simulated crystals must be included in order to reproduce measured photopeak efficiencies. In conclusion, we adopted a simple model of effective passive layers and developed heuristic methods of determining passive-layer thicknesses by comparison of simulations and measurements for a single crystal and for the full array. Prospects for future development of the model are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

G4CMP: Condensed matter physics simulation using the G$\scriptsize{\mathrm{EANT}}$4 toolkit

G4CMP simulates phonon and charge transport in cryogenic semiconductor crystals using the Geant4 toolkit. The transport code is capable of simulating the propagation of acoustic phonons as well as electron and hole charge carriers. Herein processes for anisotropic phonon propagation, oblique charge-carrier propagation, and phonon emission by accelerated charge carriers are included. The simulation reproduces theoretical predictions and experimental observations such as phonon caustics, heat-pulse propagation times, and mean charge-carrier drift velocities. In addition to presenting the physics and features supported by G4CMP, this report outlines example applications from the dark matter and quantum information science communities. These communities are applying G4CMP to model and design devices for which the energy transported by phonons and charge carriers is germane to the performance of superconducting instruments and circuits placed on silicon and germanium substrates. The G4CMP package is available to download from GitHub: github.com/kelseymh/G4CMP.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Automated reactor physics analysis framework of High Flux Isotope Reactor low-enriched uranium silicide dispersion fuel designs

The High Flux Isotope Reactor (HFIR) is a versatile research reactor that provides one of the highest steady-state neutron fluxes of any reactor in the world. The HFIR reactor physics team investigated the conversion of the current 93 wt% highly enriched uranium U 3 O 8 -Al dispersion fuel to a 19.75% low-enriched uranium (LEU) U 3 Si 2 -Al dispersion fuel. The team continuously develops a Python module to streamline the analysis steps required for an LEU core design to ensure reproducible and agile design iteration. The Python module automates the data processing between analysis steps and automates the input perturbation for branch calculations and design changes. The automated framework has proven to significantly increase the efficiency and reproducibility of the reactor physics team to design High Flux Isotope Reactor (HFIR) LEU cores and thoroughly analyze performance metrics, safety metrics, and thermal safety margins. Consequently, the team can now respond rapidly to fuel fabrication engineer and thermal-hydraulic-structural analyst requests. Numerous combinations of LEU fuel designs are explored, of which two LEU fuel designs are presented here in this paper: a low density silicide design, and a high-density silicide design. Results show that both designs meet or exceed safety and performance metrics with exception for minor differences caused by the hardened spectrum from LEU.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗