Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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

Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns: Modeling Archive

This Modeling Archive is in support of a TES-SFA publication “Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns” (Craig et al., 2021). We ran and evaluated a multi-assumption soil organic carbon (SOC) model to investigate whether alternative assumptions regarding constraints on soil microbial biomass could lead to soil carbon saturation patterns. We developed this model in the Multi-Assumption Architecture and Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_Craig2021; Walker et al. 2018). Using MAAT, we embedded three alternative hypotheses in a microbially explicit three-pool SOC model: 1) the efficiency of mineral-associated SOC formation decreases as mineral-associated SOC approaches a maximum value (“Mineral saturation”), 2) the microbial biomass turnover rate increases with increasing microbial biomass (“Density-dependent turnover”), and 3) community carbon use efficiency decreases as microbial biomass increases toward an upper limit (“Density-dependent growth”). We ran a factorial combination of these hypotheses resulting in eight models for three different classes of model (linear decay, Michaelis-Menten decay, or reverse Michaelis-Menten decay), resulting in 24 models, 12 of which are presented or discussed in the related publication. Models were parameterized using values from previous studies with similar models (Wang et al. 2013, Wieder et al. 2014, Li et al. 2014, Georgiou et al. 2017, Hassink and Whitmore 1997) and ran to an approximate steady state (200 years) at six (6) different C input rates corresponding to 0.5, 1, 2, 4, 7, and 10 times the default input value. Further model details are available in the related publication. This archive contains output from three MAAT simulations, and scripts to run these simulations and process and plot the data. Simulations are labeled “lin”, “MM_highKm”, and “RMM_highKm” reflecting factorial runs for linear, Michealis-Menten, and reverse Michaelis-Menten models, respectively. This archive contains: • 3 R scripts prepended with “init_MAAT_” to initialize runs (1 for each simulation), • 1 csv file containing years over which to run simulations (“met_year.csv”), • 1 bash script (.bs) to run MAAT, • 6 XML files that are output from MAAT describing a run (2 for each simulation), • 3 model output csv files prepended by “out_” (1 for each simulation), and • 1 analysis R script for reproducing figures 3 and 4 in Craig et al. 2021. See included user guide (Craig_2021_modeling_archive_20210315.pdf) for file organization details.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

The Biotechnology Facility for International Space Station

The primary mission of the Cellular Biotechnology Program is to advance microgravity as a tool in basic and applied cell biology. The microgravity environment can be used to study fundamental principles of cell biology and to achieve specific applications such as tissue engineering. The Biotechnology Facility (BTF) will provide a state-of-the-art facility to perform cellular biotechnology research onboard the International Space Station (ISS). The BTF will support continuous operation, which will allow performance of long-duration experiments and will significantly increase the on-orbit science throughput.

NASA Center JSC↗
DOE OSTI2020

Synthesis of energetic material particles with controlled morphology

A surfactant-assisted self-assembly method can be used to crystallize energetic materials with controlled morphology. Microparticles of hexanitrohexaazaisowurtzitane (CL-20) formed by this method may have enhanced functional reproducibility due to their monodisperse nature, and decreased shock sensitivity due to their sub-2 μm particle size.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nanometer-Scale Correlations in Aqueous Salt Solutions

The intermediate-range (1–4 nm) correlation of cations, an- ions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, S(Q), a separate region at low Q (<1 °A-1) containing local diffraction maxima (pre-peaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the pre-peaks emerge from correlations arising from existence of small quasi close-packed lattice- like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the pre-peak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes. The authors were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated by Battelle for the U.S. Department of Energy under Contract No. DE-AC05-76RL01830. Computing resources were generously allocated by PNNL's Institutional Computing program. DOE Contract No. DEAC02- 06CH11357 supports operation of the Advanced Photon Source at Argonne National Laboratory.

Chemical physics, nucleation, Madelung potential↗
NASA NTRS1989

Droplet combustion experiment drop tower tests using models of the space flight apparatus

The Droplet Combustion Experiment (DCE) is an experiment that is being developed to ultimately operate in the shuttle environment (middeck or Spacelab). The current experiment implementation is for use in the 2.2 or 5 sec drop towers at NASA Lewis Research Center. Initial results were reported in the 1986 symposium of this meeting. Since then significant progress was made in drop tower instrumentation. The 2.2 sec drop tower apparatus, a conceptual level model, was improved to give more reproducible performance as well as operate over a wider range of test conditions. Some very low velocity deployments of ignited droplets were observed. An engineering model was built at TRW. This model will be used in the 5 sec drop tower operation to obtain science data. In addition, it was built using the flight design except for changes to accommodate the drop tower requirements. The mechanical and electrical assemblies have the same level of complexity as they will have in flight. The model was tested for functional operation and then delivered to NASA Lewis. The model was then integrated into the 5 sec drop tower. The model is currently undergoing initial operational tests prior to starting the science tests.

Haggard, J. B.↗
NASA NTRS2022

A Framework for Deep Learning Emulation of Numerical Models With a Case Study in Satellite Remote Sensing

Numerical models based on physics represent the state of the art in Earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning (DL), across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in Earth systems. A difficult test for DL-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A DL emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here, we examine, with a case study in satellite-based remote sensing, the hypothesis that DL approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the DL emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of DL and the growing desire for higher resolution simulations in the Earth sciences.

Bayesian Deep Learning↗
DOE OSTI2021

Community-Driven Methods for Open and Reproducible Software Tools for Analyzing Datasets from Atom Probe Microscopy

Atom probe tomography, and related methods, probe the three-dimensional architecture of a material. The software tools that microscopists use, and how these tools are connected into workflows, makes a substantial contribution to the accuracy and precision of such a material characterization experiment. Typically, we adapt methods from other communities like mathematics, data science, computational geometry, artificial intelligence, or scientific computing. We also realize that improving on research data management is a challenge when it comes to align with the FAIR data stewardship principles. Faced with this global challenge, we are convinced that collaborating is useful. Here, we report the results and challenges with an inter-laboratory call for developing test cases for several types of atom probe software tools. The results support why defining detailed recipes of software workflows and sharing these recipes is necessary and rewarding: Open source tools and (meta)data exchange can help to make our day-to-day data processing tasks become more efficient, the training of new users and knowledge transfer become easier, and assist us with automated quantification of uncertainties to gain access to substantiated results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vacancies in graphene: an application of adiabatic quantum optimization

Quantum annealers have grown in complexity to the point that quantum computations involving a few thousand qubits are now possible. In this paper, with the intentions to show the feasibility of quantum annealing to tackle problems of physical relevance, we used a simple model, compatible with the capability of current quantum annealers, to study the relative stability of graphene vacancy defects. By mapping the crucial interactions that dominate carbon-vacancy interchange onto a quadratic unconstrained binary optimization problem, our approach exploits the ground state as well as the excited states found by the quantum annealer to extract all the possible arrangements of multiple defects on the graphene sheet together with their relative formation energies. Furthermore, this approach reproduces known results and provides a stepping stone towards applications of quantum annealing to problems of physical–chemical interest.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Optimal adjustment sets for causal query estimation in partially observed biomolecular networks

Abstract Causal query estimation in biomolecular networks commonly selects a ‘valid adjustment set’, i.e. a subset of network variables that eliminates the bias of the estimator. A same query may have multiple valid adjustment sets, each with a different variance. When networks are partially observed, current methods use graph-based criteria to find an adjustment set that minimizes asymptotic variance. Unfortunately, many models that share the same graph topology, and therefore same functional dependencies, may differ in the processes that generate the observational data. In these cases, the topology-based criteria fail to distinguish the variances of the adjustment sets. This deficiency can lead to sub-optimal adjustment sets, and to miss-characterization of the effect of the intervention. We propose an approach for deriving ‘optimal adjustment sets’ that takes into account the nature of the data, bias and finite-sample variance of the estimator, and cost. It empirically learns the data generating processes from historical experimental data, and characterizes the properties of the estimators by simulation. We demonstrate the utility of the proposed approach in four biomolecular Case studies with different topologies and different data generation processes. The implementation and reproducible Case studies are at https://github.com/srtaheri/OptimalAdjustmentSet.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The Gut Microbiome and Alcoholic Liver Disease: Ethanol Consumption Drives Consistent and Reproducible Alteration in Gut Microbiota in Mice

Phenotypic health effects, both positive and negative, have been well studied in association with the consumption of alcohol in humans as well as several other mammals including mice. Many studies have also associated these same health effects and phenotypes to specific members of gut microbiome communities. Here we utilized a chronic plus binge ethanol feed model (Gao-binge model) to explore microbiome community changes across three independent experiments performed in mice. We found significant and reproducible differences in microbiome community assemblies between ethanol-treated mice and control mice on the same diet absent of ethanol. We also identified significant differences in gut microbiota occurring temporally with ethanol treatment. Peak shift in communities was observed 4 days after the start of daily alcohol consumption. We quantitatively identified many of the bacterial genera indicative of these ethanol-induced shifts including 20 significant genera when comparing ethanol treatments with controls and 14 significant genera based on temporal investigation. Including overlap of treatment with temporal shifts, we identified 25 specific genera of interest in ethanol treatment microbiome shifts. Shifts coincide with observed presentation of fatty deposits in the liver tissue, i.e., Alcoholic Liver Disease-associated phenotype. The evidence presented herein, derived from three independent experiments, points to the existence of a common, reproducible, and characterizable “mouse ethanol gut microbiome”.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2015

Characterizing Young Giant Planets with the Gemini Planet Imager: An Iterative Approach to Planet Characterization

After discovery, the first task of exoplanet science is characterization. However experience has shown that the limited spectral range and resolution of most directly imaged exoplanet data requires an iterative approach to spectral modeling. Simple, brown dwarf-like models, must first be tested to ascertain if they are both adequate to reproduce the available data and consistent with additional constraints, including the age of the system and available limits on the planet's mass and luminosity, if any. When agreement is lacking, progressively more complex solutions must be considered, including non-solar composition, partial cloudiness, and disequilibrium chemistry. Such additional complexity must be balanced against an understanding of the limitations of the atmospheric models themselves. For example while great strides have been made in improving the opacities of important molecules, particularly NH3 and CH4, at high temperatures, much more work is needed to understand the opacity of atomic Na and K. The highly pressure broadened fundamental band of Na and K in the optical stretches into the near-infrared, strongly influencing the spectral shape of Y and J spectral bands. Discerning gravity and atmospheric composition is difficult, if not impossible, without both good atomic opacities as well as an excellent understanding of the relevant atmospheric chemistry. I will present examples of the iterative process of directly imaged exoplanet characterization as applied to both known and potentially newly discovered exoplanets with a focus on constraints provided by GPI spectra. If a new GPI planet is lacking, as a case study I will discuss HR 8799 c and d will explain why some solutions, such as spatially inhomogeneous cloudiness, introduce their own additional layers of complexity. If spectra of new planets from GPI are available I will explain the modeling process in the context of understanding these new worlds.

exoplanets↗
NASA NTRS2019

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: An Application to XCO2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

Calle, Leonardo↗
NASA NTRS2019

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: an Application to Xco2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

segmentation algorithm↗
DOE OSTI2021

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparative analysis of plasmon modes in layered Lindhard metals and strange metals

The enigmatic strange metal remains one of the central unsolved problems of 21st century science. Understanding this phase of matter requires knowledge of the momentum- and energy-resolved dynamic charge susceptibility 𝜒⁡(𝑞,𝜔), especially at finite momentum. Inelastic electron scattering (EELS), performed in either transmission or reflection geometry, is a powerful probe of 𝜒⁡(𝑞,𝜔). For the prototypical strange metal Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝑥 , transmission- and reflection EELS, and infrared (IR) spectroscopy agree at 𝑞∼0, all revealing a highly damped plasmon near 1 eV. At larger 𝑞, however, EELS results show unresolved discrepancies. Since IR data are highly reproducible, it is advantageous to use IR data to calculate what the expected EELS response should be at modest 𝑞. Building on prior momentum-resolved reflection geometry M-EELS work [Chen et al., Phys. Rev. B 109, 045108 (2024)], we extend this approach to transmission EELS for finite stacks of metallic layers, comparing a “textbook” Lindhard metal to a strange metal. In the Lindhard case, the low-𝑞 response is dominated by long-lived, standing wave plasmon modes arising from interlayer Coulomb coupling, with in-plane dispersions that resemble the well-known Fetter modes of layered metals. This behavior depends only on the geometry and the long-range nature of the Coulomb interaction and is largely insensitive to layer details. At larger 𝑞, the response reflects the microscopic properties of individual layers. For the strange metal, calculations based on IR data predict a highly damped plasmon with weak dispersion and no distinct surface mode. While our results match IR and M-EELS at low 𝑞, they do not reproduce any published EELS spectra at large 𝑞, highlighting unresolved discrepancies that demand further experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗