Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient sampling”

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 OSTI2023

Salinity exposure affects lower-canopy specific leaf area of upland trees in a coastal deciduous forest

Sea level rise and increasing storm surges are likely to affect the canopy physiology, ecology, and structure of coastal forests, even well in advance of tree mortality. Laboratory and greenhouse studies have documented that saltwater exposure can trigger changes in leaf-level physiology and morphology, but few in situ studies have examined how tree-specific leaf area (SLA), the ratio of leaf area to mass and a crucial trait and model parameter, is affected by saline soils. We conducted an observational study of SLA in a mid-Atlantic (USA) coastal deciduous forest, taking advantage of a natural gradient in salinity along a tidal creek. Measured SLA of the 239 trees and seven species sampled ranged from Carya glabra (N = 6 trees, mean SLA = 277.9±36.3cm 2 /g) to Fagus grandifolia (N=60, 321.9±62.9cm 2 g); as expected, trees species and canopy position (sun versus shade) significantly affected SLA. For trees (N=100) directly exposed to the tidal creek, salinity was highly significant after accounting for species (P<0.001), with trees in the lower reaches of the creek having lower SLA. Leaf area index (LAI), computed from SLA and litter traps, ranged from 4.8 to 15.8 and was inversely related to salinity exposure; the spatial variability in leaf litter production contributed much more to LAI uncertainty than did SLA variability. These in situ results are correlative but consistent with the hypothesis, based on previous greenhouse studies, that the stress of chronic salinity exposure changes species’ leaf morphology. Our findings are useful for understanding the growing effects of saltwater intrusion into upland forests, as well as parameterizing and testing ecosystem-scale models simulating forest stressors and disturbances at the terrestrial-aquatic interface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impacts of Degradation on Water, Energy, and Carbon Cycling of the Amazon Tropical Forests

Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high-resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED-2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED-2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data used in: "Salinity exposure affects lower-canopy specific leaf area of upland trees in a coastal deciduous forest"

We conducted an observational study of SLA in a mid-Atlantic (USA) coastal deciduous forest, taking advantage of a natural gradient in salinity along a tidal creek. Measured SLA of the 239 trees and seven species sampled ranged from Carya glabra (N = 6 trees, mean SLA = 277.9 ± 36.3 cm2 g-1) to Fagus grandifolia (N = 60, 321.9 ± 62.9 cm2 g-1). This work is in press at Forest Ecology and Management.This dataset has two files, "sla_data.csv" (the data) and "sla_metadata.txt" (the metadata). Both are plain-text files and can be read by most software such as text editors, R, Python, Excel, etc. The data file is comma-separated values with 11 columns and 239 rows of data. The metadata file has 11 rows, one for each data column, and giving a plain-text description for each column and any associated units.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An Online Dynamic Amplitude-Correcting Gradient Estimation Technique to Align X-ray Focusing Optics

High-brightness X-ray pulses, as generated at synchrotrons and X-ray free electron lasers (XFELs), are used in a variety of scientific experiments. At these facilities, measurements often require optical equipment, e.g Compound Refractive Lenses (CRLs) to be precisely aligned and focused. The lateral alignment of CRLs to a beamline requires precise positioning along four axes: two translational, and the two rotational. At a synchrotron, alignment is often accomplished manually. However, XFEL beamlines present a beam brightness that fluctuates stochastically, making manual alignment a time-consuming endeavor. Automation using simplex or classic stochastic descent often fails, given the errant gradient estimates. Herein we present a dynamic-amplitude correction to the usual gradient based on the combination of a generalized finite difference stencil and a time-dependent sampling pattern. Intensity is recorded periodically, then used to normalize numerical derivatives against fluctuations. Error expectation is analyzed, and efficacy is demonstrated on classic benchmarks. We provide a proof of concept by laterally aligning optics on a simulated XFEL beamline using data recorded at both synchrotron and XFEL facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

BEYONDPLANCK II. CMB mapmaking through Gibbs sampling

We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Temperature distribution in a laser-heated diamond anvil cell as described by finite element analysis

Finite element analysis (FEA) is a powerful tool for numerically solving partial differential equations over complex geometries and is thus useful for analyzing heat transport in laser-heated diamond anvil cell (LHDAC) experiments. Our models expand on previously published simulations by calculating the volume-averaged temperatures of both the sample and insulation/pressure media under steady-state heating to determine the thermal pressure of the hot sample. Our goal is to produce an accurate relationship between the measured surface temperature of the absorbing sample and the temperature of the transparent insulating media, which is used to determine thermal pressure but susceptible to steep temperature gradients. We find that in doing so, our FEA models of temperature within the pressure/insulation media can differ from simplified estimates of temperature gradients by more than a factor of 2. We also explore temperature-dependent and temperature-independent thermal conductivity models and find that the volume-averaged temperatures differ by up to a factor of 1.3, forcing the predicted thermal pressures determined to also differ by up to a factor of 1.5 at a temperature of 2000 K at 50 GPa for neon. Higher temperatures exacerbate this difference. We also find that unintentional asymmetric sample insertion and sample heating, which are common in LHDAC experiments, do not have a first-order effect on volume-averaged temperatures. The FEA models, available in both Python and FlexPDE, are versatile across different sample geometries, materials, and heat source laser shapes.

Farah, Frederick↗
DOE OSTI2025

Multi-fidelity is the new annealing: Gradient-free learning of posterior densities via transport maps

To tackle concentrated, multi-modal Bayesian inference problems, we propose using an annealed importance sampling procedure. To do this, we form a sequence of annealed distributions and employ transport maps to act as a surrogate of each distribution. This process is demonstrated on a few examples, favorably showing its potential for efficiently parallelizing the process of PDE evaluations and allowing for surrogates that we can sample from exactly.

van Bloemen Waanders, Bart G [Sandia National Labo↗
DOE OSTI2023

Efficient Subset Simulation using Hamiltonian Neural Network enhanced Markov Chain Monte Carlo Methods

The Monte Carlo method delivers an unbiased estimate of the probability of failure. However, the variance of the estimate depends on the number of evaluated samples. This number must be very large for estimations of a low probability of failure. If the evaluation of each sample is computationally expensive, the crude Monte Carlo simulation strategy is impracticable. Therefore, subset simulations are used to reduce the required number of evaluations. Subset simulations require a Markov Chain Monte Carlo sampler, such as the random walk Metropolis-Hastings algorithm. The algorithm, however, struggles with sampling in low-probability regions, especially if they are narrow. As a consequence, advanced Markov Chain Monte Carlo simulations have been developed. In particular, the Hamiltonian Monte Carlo method explores the target distribution rapidly. Driven by the idea of Hamiltonian dynamics, this sampler provides a non-random walk through the target distribution. The incorporation of subset simulation and Hamiltonian Monte Carlo methods has shown promising results for reliability analysis. One downside of the Hamiltonian Monte Carlo method is that gradient evaluations are computationally expensive, especially when dealing with high-dimensional problems and evaluating long trajectories. We show that integrating Hamiltonian neural networks in Hamiltonian Monte Carlo simulations significantly speeds up the sampling task. Furthermore, the enhancement of adaptive trajectory length within the Hamiltonian Monte Carlo results in the efficient proposal of the following states. Based on this recent enhancement, we provide a fast sampling strategy for subset simulations using Hamiltonian neural networks to replace the evaluation of the gradient and significantly speed up the Hamiltonian Monte Carlo simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First Measurement of Z Opacity Sample Evolution near Solar Interior Conditions Using Time-Resolved Spectroscopy

Opacity model differences with Fe opacity measurements at high temperature (T>1⁢8⁢0 eV ) and high electron density (𝑛 𝑒 >3×10 22 cm −3 ) at nearly solar interior conditions have remained unresolved [Bailey et al., Nature 517, 56 (2015) and Nagayama et al., Phys. Rev. Lett. 122, 235001 (2019)]. Systematic errors from temporal gradients are one hypothesis for the discrepancy. Past data recorded on x-ray film provided spectral measurements over a time determined by the 2.8-ns backlighter duration. Here, we present the first measurements of opacity sample temporal evolution using novel hCMOS ultrafast x-ray camera technology. The measured conditions, measured backlighter time history, and modeled opacities are used to show that temporal gradients do not resolve the model-data discrepancy. The methods demonstrated provide potential advantages, including opacities at more extreme conditions, spectral line shift measurements, and improved measurements at other facilities.

atomic spectra↗
DOE OSTI2022

Label-Free Profiling of up to 200 Single-Cell Proteomes per Day Using a Dual-Column Nanoflow Liquid Chromatography Platform

Single-cell proteomics (SCP) has great potential to advance biomedical research and personalized medicine. The sensitivity of such measurements increases with low-flow separations (<100 nL/min) due to improved ionization efficiency, but the time required for sample loading, column washing, and regeneration in these systems can lead to low measurement throughput and inefficient utilization of the mass spectrometer. Herein, we developed a two-column liquid chromatography (LC) system that dramatically increases the throughput of label-free SCP using two parallel subsystems to multiplex sample loading, online desalting, analysis, and column regeneration. The integration of MS1-based feature matching increased proteome coverage when short LC gradients were used. The high-throughput LC system was reproducible between the columns, with a 4% difference in median peptide abundance and a median CV of 18% across 100 replicate analyses of a single-cell-sized peptide standard. An average of 621, 774, 952, and 1622 protein groups were identified with total analysis times of 7, 10, 15, and 30 min, corresponding to a measurement throughput of 206, 144, 96, and 48 samples per day, respectively. When applied to single HeLa cells, we identified nearly 1000 protein groups per cell using 30 min cycles and 660 protein groups per cell for 15 min cycles. Finally, we explored the possibility of measuring cancer therapeutic targets with a pilot study comparing the K562 and Jurkat leukemia cell lines. This work demonstrates the feasibility of high-throughput label-free single-cell proteomics.

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

Understanding the role of segmentation on process-structure–property predictions made via machine learning

Here, the present study investigated the effect of porosity surface determination methods on performance of machine learning models used to predict the tensile properties of AlSi10Mg processed by laser powder bed fusion from micro-computed tomography data. Machine learning models applied in this work include support vector machines, neural networks, decision trees, and Bayesian classifiers. The effects of isosurface thresholding and local gradient approaches for porosity segmentation, as well as image filtering schemes, on model precision were evaluated for samples produced under differing levels of global energy density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Soil Particulate Organic Matter is Related to Ericoid Mycorrhizal Shrubs, not Ectomycorrhizal Trees, in a Temperate Forest

Mycorrhizal associations are key drivers of soil biogeochemistry, but previous studies have focused almost exclusively on ectomycorrhizal (EcM) and arbuscular mycorrhizal (AM) associations. Ericoid mycorrhizal (ErM) shrubs frequently occur in forest understories and are expanding in response to disturbance, but are rarely considered in biogeochemical frameworks. We investigated the relationships of understory ErM shrubs and overstory trees on carbon (C) and nitrogen (N) in soil organic matter fractions in a southern Appalachian temperate forest. We sampled the 0–10 cm mineral soil layer from 43 plots at the Coweeta Hydrologic Laboratory, across gradients in overstory EcM dominance and understory ErM shrub biomass. Soil C:N ratios increased with both increasing EcM dominance and increasing ErM shrub biomass. However, total particulate organic matter (POM) C, and the proportion of C and N held in POM increased with increasing ErM shrub biomass, but not with increasing EcM dominance. In contrast, mineral-associated organic matter (MAOM) C and N were negatively associated with EcM dominance, but were not related to ErM shrubs. Our findings suggest that ErM shrubs facilitate POM formation while AM trees promote MAOM formation. Because ErM shrub biomass represents a small fraction of total forest biomass, our work provides evidence that ErM shrubs have an outsized effect on soil organic matter, which advocates for their inclusion in mechanistic studies and biogeochemical frameworks.

Bonilla, Kayla A. [University of Georgia, Athens, ↗
DOE OSTI2022

Strain-rate dependent deformation mechanisms in single-layered Cu, Mo, and multilayer Cu/Mo thin films

Here, strain-rate sensitivity and rate-dependent hardness, over a range of 10 -2 to 10 2 s -1 , of sputter-deposited single-layered Cu, Mo, and 5 nm Cu/ 5 nm Mo, and 100 nm Cu/ 100 nm Mo multilayer films with a total film thickness of 5 μm were measured using nanoindentation. The plastic zone underneath the nanoindents was characterized via cross-sectional transmission electron microscopy (XTEM). The multilayer films exhibited enhanced hardness but slightly reduced strain-rate sensitivity with decreasing layer thickness from 100 nm to 5 nm. Only the 5 nm Cu/ 5 nm Mo multilayer film exhibited shear bands underneath the nanoindents, and the size of the shear bands increased with increasing strain rate. In contrast, the 100 nm Cu/ 100 nm Mo multilayer film exhibited material pile-up around the indents and significant nanotwinning within Cu grains. The effect of strain rate and layer thickness on the hardness and strain rate sensitivity of the multilayer thin films is interpreted using a modified confined layer slip (CLS) model. The reduced rate sensitivity at 5 nm as compared to 100 nm correlates with abundant growth nanotwins in the Cu grains in 100 nm and formation of shear bands in 5 nm multilayers. In single layer films, a substructure with a high density of dislocations was observed consistent with the plastic strain gradient in the indent plastic zone. No evidence of deformation twins was noted in any of the samples.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Real-time breath analysis towards a healthy human breath profile

Abstract The direct analysis of molecules contained within human breath has had significant implications for clinical and diagnostic applications in recent decades. However, attempts to compare one study to another or to reproduce previous work are hampered by: variability between sampling methodologies, human phenotypic variability, complex interactions between compounds within breath, and confounding signals from comorbidities. Towards this end, we have endeavored to create an averaged healthy human ‘profile’ against which follow-on studies might be compared. Through the use of direct secondary electrospray ionization combined with a high-resolution mass spectrometry and in-house bioinformatics pipeline, we seek to curate an average healthy human profile for breath and use this model to distinguish differences inter- and intra-day for human volunteers. Breath samples were significantly different in PERMANOVA analysis and ANOSIM analysis based on Time of Day, Participant ID, Date of Sample, Sex of Participant, and Age of Participant ( p < 0.001). Optimal binning analysis identify strong associations between specific features and variables. These include 227 breath features identified as unique identifiers for 28 of the 31 participants. Four signals were identified to be strongly associated with female participants and one with male participants. A total of 37 signals were identified to be strongly associated with the time-of-day samples were taken. Threshold indicator taxa analysis indicated a shift in significant breath features across the age gradient of participants with peak disruption of breath metabolites occurring at around age 32. Forty-eight features were identified after filtering from which a healthy human breath profile for all participants was created.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Training Restricted Boltzmann Machines With a D-Wave Quantum Annealer

Restricted Boltzmann Machine (RBM) is an energy-based, undirected graphical model. It is commonly used for unsupervised and supervised machine learning. Typically, RBM is trained using contrastive divergence (CD). However, training with CD is slow and does not estimate the exact gradient of the log-likelihood cost function. In this work, the model expectation of gradient learning for RBM has been calculated using a quantum annealer (D-Wave 2000Q), where obtaining samples is faster than Markov chain Monte Carlo (MCMC) used in CD. Training and classification results of RBM trained using quantum annealing are compared with the CD-based method. The performance of the two approaches is compared with respect to the classification accuracies, image reconstruction, and log-likelihood results. The classification accuracy results indicate comparable performances of the two methods. Image reconstruction and log-likelihood results show improved performance of the CD-based method. It is shown that the samples obtained from quantum annealer can be used to train an RBM on a 64-bit “bars and stripes” dataset with classification performance similar to an RBM trained with CD. Though training based on CD showed improved learning performance, training using a quantum annealer could be useful as it eliminates computationally expensive MCMC steps of CD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spatial variation in soil microbial processes as a result of woody encroachment depends on shrub size in tallgrass prairie

Aims: As woody plants encroach into grassland ecosystems, we expect that altered plant-soil interactions will lead to changes in the microbial processes that affect carbon storage and nutrient cycling. Specifically, this research aimed to address how (1) soil chemistry, (2) microbial nutrient demand, and (3) the rate and source of potential soil C mineralization vary spatially under individual woody shrubs of varying size within a mesic grassland. Methods: Here, we collected soil samples from the center, the midpoint between the center and edge, the edge, and the shrub-grass ecotone of multiple Cornus drummondii shrubs across a shrub-size gradient in infrequently burned tallgrass prairie. Results: We found total soil carbon and total soil nitrogen increased with shrub size in every location but the edge. Microbial demand for nitrogen also increased as shrubs increased in size. Across all shrub sizes and sampling locations, potential soil carbon mineralization rates were higher when microbes broke down proportionally more shrub-derived (C 3 ) organic matter than grass-derived (C 4 ) organic matter. Conclusions: Our results suggest that the spatio-temporal context of woody encroachment is critical for understanding its impact on belowground microbial processes. In this ecosystem, a longer period of occupancy by woody plants increases potentially mineralizable carbon.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Densification, microstructure, and mechanical properties of Mo–30W alloys fabricated from conditioned powders

Refractory alloys, such as molybdenum-based systems, are attracting growing interest for applications in extreme environments, such as in the nuclear and aerospace industries. Recent advances in sintering technologies, coupled with mechanical alloying, have enabled the tailored design of these alloys by leveraging powder characteristics to control final microstructures and mechanical properties. In this study, Mo-30W alloys were fabricated using electric field-assisted sintering (EFAS) from ball-milled powders with and without hydrogen treatment to investigate the influence of surface oxides on material properties and sintering behavior. The results revealed that samples processed from as-ball-milled powder contained a high density of oxides within the microstructure, whereas oxide presence was significantly reduced in samples fabricated from hydrogen-treated powders. Interestingly, the two powder types led to opposite trends in grain size distribution: samples from untreated powders exhibited grain refinement from sample periphery to the center, while samples from hydrogen-treated powders showed grain coarsening toward the center. This behavior is attributed to temperature gradients present during sintering due to electrical percolation pathway differences during Joule heating. The powder surface oxides may have influenced the temperature distribution and grain evolution. Microhardness profiles measured along both axial and thickness directions were consistent with the grain size distribution. Furthermore, oxide films on powder surfaces have delayed densification by hindering particle necking and atomic diffusion during sintering.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Self-Sealing Mafic Sills for Carbon and Hydrogen Storage

Tabular igneous intrusions (sills) are common features in sedimentary basins and have the potential to be useful seals for fluids (e.g., CO 2 or H 2 ) in geological storage scenarios, and may be important as the need for geologic carbon sequestration and alternative fuel storage increases. This is advantageous in regions without ready access to large-volume reservoirs in depleted hydrocarbon plays and saline aquifers, like the northeastern United States. Moreover, geological H 2 storage requires more demanding conditions than typically associated with oil and gas. An enhanced seal will have the properties of a traditional seal - competent, low permeability, and laterally extensive - with the addition of being able to self-seal pre-existing and induced fractures. Self-sealing will occur along fluid pathways like fractures where CO 2 , water, and minerals like plagioclase, olivine, and pyroxene react together. Dolerite sills from the Gettysburg Basin, Pennsylvania, have remarkably low permeability and homogeneous compositions that include minerals that will readily react with CO 2 dissolved in water. Here, we characterize the physical properties and chemical gradients within several mafic sills cored in five boreholes. In addition, preliminary CO 2 -reaction experiments on dolerite samples demonstrated rapid carbonate mineralization.

58 GEOSCIENCES↗