Microcomputed tomography (laboratory and synchrotron) of intact archeological human second metacarpa
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Abstract A substantial fraction of nitrogen (N) fertilizer applied in agricultural systems is not incorporated into crops and moves below the rooting zone as nitrate (NO 3 − ). Understanding mechanisms for soil N retention below the rooting zone and leaching to groundwater is essential for our ability to track the fate of added N. We used dual stable isotopes of nitrate (δ 15 N–NO 3 − and δ 18 O–NO 3 − ) and water (δ 18 O–H 2 O and δ 2 H–H 2 O) to understand the mechanisms driving nitrate leaching at three depths (0.8, 1.5, and 3.0 m) of an irrigated corn field sampled every 2 weeks from 2016 to 2020 in the southern Willamette Valley, Oregon, USA. Distinct periods of high nitrate concentrations with lower δ 15 N–NO 3 − values indicated that a portion of that nitrate was from recent fertilizer applications. We used a mixing model to quantify nitrate fluxes associated with recently added fertilizer N versus older, legacy soil N during these “fertilizer signal periods.” Nitrate leached below 3.0 m in these periods made up a larger proportion of the total N leached at that depth (∼52%) versus the two shallower depths (∼13%–16%), indicating preferential movement of recently applied fertilizer N through the deep soil into groundwater. Further, N associated with recent fertilizer additions leached more easily when compared to remobilized legacy N. A high volume of fall and winter precipitation may push residual fertilizer N to depth, potentially posing a larger threat to groundwater than legacy N. Optimizing fertilizer N additions could minimize fertilizer losses and reduce nitrate leaching to groundwater.
Water heating is the third largest electricity consumer in U.S. households, after space heating and cooling. Thus, water heaters represent a significant potential for reducing electricity consumption and associated CO2 emissions of residential buildings. To this end, this study proposes a model-free deep reinforcement learning (RL) approach that aims to minimize the electricity consumption and the CO2 emissions of a heat pump water heater without affecting user comfort. In this approach, a set of RL agents focusing on either electricity saving or emission reduction, with different look ahead periods, were trained using the deep Q-networks (DQN) algorithm and their performance was tested on different hot water usage and Marginal Operating Emissions Rate (MOER) profiles. The testing results showed that the RL agents that focus on electricity saving can save electricity in the range of 12–22% by operating the water heater with maximum heat pump efficiency and minimum electric element utilization. On the other hand, the RL agents that focus on emission reduction reduced emissions in the range of 18–37% by making use of the variable MOER values. These RL agents used the heat pump and/or an element when the MOER values are low due to the availability of renewable energy sources (e.g., solar and wind) and mostly avoided the periods of carbon-intensive periods. Overall, these results showed that the proposed RL approach can help minimize the electricity consumption and the CO2 emissions of a heat pump water heater without having any prior knowledge about the device.
Infrasound signals from open-vent volcanoes can contain information related to that volcano’s eruption dynamics and future activity. Studying a specific volcano’s acoustic signature during periods of relative quiescence is thus important for understanding potential warning signs of impending eruptions. Volcán Villarrica, located in southern Chile, has a long-lived active lava lake that produces continuous infrasound with spectral peaks centered around 1 Hz that vary by + / - ~ 0.2 Hz over day-to-week-long time scales. The infrasound frequency content has been shown to illuminate key volcanic properties such as eruption style and crater shape. Leading up to Villarrica’s most recent paroxysm in 2015, for instance, infrasound spectral changes coincided with a rise in Villarrica’s lava lake level. Quantifying and understanding the regular fluctuation in recorded infrasonic frequencies from Villarrica and other open-vent volcanoes is thus imperative during periods of relative calm. Our experiment entailed a week-long period of crater rim infrasound observations associated with stable, open-vent activity and revealed two independent source processes: spatter bursting events (relatively broadband infrasound) and lava lake-induced tremor (sustained signals peaked at 1 Hz). A comparison of these acoustic signals with results from 3-D finite-difference time-domain wave propagation models shows that while sound speed variability can influence Villarrica’s spectrum, it cannot explain the full extent of the observed frequency excursions. We conclude that source spectrum variability is primarily responsible for the frequency excursions. In conclusion, our work highlights the utility of data collected from a dense distribution of twenty infrasound sensors operating at the summit coupled with numerical modeling of sound radiation.
Abstract Evaluation of long-term temporal and spatial climatic change in mountainous regions is a critical challenge because of the interactive effects of multiple land and climatic factors and processes. Here we present the application of the statistical framework to the assessment of changes of climatic conditions, using data from 17 meteorological stations across the East River watershed near Crested Butte, Colorado, USA, and spanning the period from 1966 to 2021. The framework is developed based on (1) a time-series analysis of daily, monthly, and yearly averaged meteorological parameters (temperature, relative humidity, precipitation, wind speed, etc.), (2) evaluation and time series analysis of potential evapotranspiration (ET o ), actual evapotranspiration (ET), aridity index (AI), standard precipitation index (SPI) and standard precipitation-evapotranspiration index (SPEI), and (3) a temporal-spatial climatic zonation of the studied area based on the hierarchical clustering and PCA analysis of the SPEI, because the SPEI can be considered an integrative characteristic of the changes of climatic conditions. The Budyko model, with the application of the Penman–Monteith equation for the estimation of ET o , was used to determine the ET. The time series analysis of the AI is used to identify the periods with energy limited and water limited conditions. Hierarchical clustering of site locations for the three temporal segments of the SPEI showed a significant temporal-spatial shifts, indicating that dynamic climatic processes drive zonation patterns. Therefore, the watershed climatic zonation requires periodic re-evaluation based on the structural time series analysis of meteorological and water balance data.
Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).
Elucidating climatic impacts on stream nutrient export and stoichiometry will improve the understanding of forest carbon (C) storage in a warmer world. We constructed C, nitrogen (N), and phosphorus (P) budgets in four watersheds within a rain-snow transition site and another four within a higher-elevation, snow-dominated site, in California’s mixed-conifer zone. We used these two sites in a space-for-time substitution to assess the potential effects of warming on nutrient cycles in currently snow-dominated areas that will become more rain-dominated. During a non-drought period (water year (WY) 2004 – 2011), mean annual stream exports of C and N in particulate forms at the transition site were double that at the snow-dominated site, suggesting that sediment-associated nutrient losses may increase with warming. The transition site had 12% lower N but double the P content in mineral horizons, lower N:P mass ratios in organic horizons, and lower stream export of dissolved inorganic N than the snow-dominated site. Furthermore, these differences suggest that montane forests may receive lower supply of N relative to P with warming. In addition, given strong interests in forest thinning to increase drought resiliency, we examined changes in stream nutrient export after thinning and during a major drought period (WY 2013 – 2015). Stream exports of C, N, and P were similar between unthinned and thinned watersheds during drought, suggesting negligible thinning impacts on stream nutrient export during excessively dry periods. Taken together, our results suggest that as the climate warms, California’s montane forests may lose more nutrients through erosion and receive less N input than P.
Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.
Methane oxidation under periodic conditions and the oxygen storage capacity of a bilayer Pt/Pd/Al 2 O 3 over a Mn 0.5 Fe 2.5 O 4 spinel catalyst were studied before and after SO 2 exposure, and after simulated regeneration conditions. Prior to sulfur exposure, improvement in CH 4 oxidation conversion under periodic conditions compared to steady-state conditions was observed. After sulfur exposure at 100 °C, there was a loss in CH 4 oxidation performance and a loss of oxygen storage capacity of the spinel material. Here, the extent of regeneration from sulfur poisoning depends on the ability to induce the decomposition of sulfate species, and while all regeneration methods tested in this study did improve CH 4 conversion, regeneration methods under periodic conditions induced greater sulfur species desorption from the catalyst surface leading to improved CH 4 con version. Key regeneration parameters – temperature, feed composition, modulation amplitude and frequency – were optimized to induce S species decomposition and correlated to CH 4 oxidation activity recovery.
This study investigates the technical and economic feasibility of using high levels of solar energy penetration up to 400 MW into a smart grid system of 60,000 smart houses. A novel non-cooperative Stackelberg game is introduced that incorporates the profitability of the supply-side and helps in solving problems related to overgeneration and photovoltaic curtailment. The non-cooperative game is intended to find the optimal dynamic prices that would leverage distributed storage through the demand-side to stabilize the power grid operation. Ten cases are studied with five photovoltaic plant sizes and two battery designs. Here, a novel quantitative analysis of high levels of solar penetration as a percentage of the total electricity demand is introduced to evaluate the technical feasibility of the studied cases. To evaluate the economic viability of the proposed smart grid system, four metrics were used: the levelized cost of energy, the levelized cost of storage, the payback period, and the net present value. Two out of ten studied cases were concluded to be the most promising cases, one with a solar photovoltaic plant size of 200 MW and the other with 300 MW. The case with 300 MW solar plant is preferred as it paves the way for more solar energy deployment with a solar penetration percentage up to 67.78%. This case had a payback period of 10.72 years and a net present value of $\$51.44$ M for the solar plant and a payback period of 12.06 years and a net present value of $\$40.75$ M for the demand-side.
Joint and marginal distributions in the frequency, direction, and time domain are employed to demonstrate their value for wave energy resource characterization when full spectra are available. Insights gained through analysis of these distributions support wave energy converter concept design, operation and maintenance. Spatial trends in the wave energy resource and contributing wave energy systems along the continental shelf of the West Coast of the United States are investigated using the most recent two-dimensional wave spectra measurements at four buoys over an eleven year period (2008 to 2018). Resource hot spots and dominant resolved energy resource bands in the frequency-direction-time domain are delineated. Resource attributes, including frequency and directional spreading, and seasonal variability, are characterized using joint distributions and marginal distributions of wave power spectra. North Pacific westerly swells in the winter season, augmented by Aleutian low-pressure southwesterly swells, are the principal suppliers of the dominant resource and main drivers influencing resource attributes. In conclusion, the modification of these systems southward, especially the North Pacific westerly swells, explains the observed spatial resource trends. The dominant resource wave period shifts two seconds to higher wave periods, thirty degrees in the dominant direction band to a more northward orientation, and forward by one month.
Solar photovoltaic (PV) and wind energy systems are crucial for decarbonizing the electricity sector and achieving climate goals. However, these systems are weather-dependent, and ignoring the potential changes in their generation levels due to climate change could compromise achieving climate targets and meeting future electricity demand. This study evaluates the impact of climate change on the generation potential of wind and solar PV systems in India for three future periods, 2030 (2021-2040), 2050 (2041-2060), and 2070 (2061-2080) compared to the baseline year 2000 (1991-2010), under three emission scenarios: SSP245, SSP370, and SSP585. Solar PV generation levels consistently decline (up to 10 %) across all regions and scenarios. Wind energy shows more pronounced variability (-20 % to 30 %). The South and Southeastern regions of India show improvements in wind potential across all scenarios and time periods. This study also investigated the projected changes in the generation droughts of both energy systems. For solar PV, drought days increase across most regions (exceeding 500 days under SSP370 across the 20-year period). In contrast, wind energy sees a reduction in drought days, especially in parts of South and Southeast India (declines exceeding 50 days across different scenarios). For both energy systems, the patterns of generation drought and generation potential are similar, and indicate that Western and Northern India may be less favorable for the future expansion of solar PV and wind energy, respectively. These results highlight the need to account for the potential impacts in future capacity planning.
Actinium is an elusive element with untamed properties and represents a peculiar case in the periodic table, as its isotopes are all radioactive, the longest-lived one having only a 22-year half-life, and the availability of actinium isotopes remains very low (microgram level, at best), hindering research on its compounds. Despite being a natural element discovered more than 120 years ago, and despite an increasing interest in using one of its isotopes ( 225 Ac) for highly efficient cancer therapies, the chemistry of actinium is still largely unknown relative to other elements. Since Ac is the first element of the actinide series, it is accepted that its ion, Ac 3+ , is the most voluminous trivalent cation of the periodic table. However, the structural data available on Ac 3+ compounds are scarce and have mainly been collected in the 1940-1960's, when actinide chemistry was still in its infancy, and have not been put in perspective with the advances in the chemistry of other elements, making it difficult to accurately evaluate its actual size and coordination chemistry. Here, we review progress made on the chemistry of lanthanides and actinides and reevaluate the structural data published on Ac 3+ since the era of the Manhattan Project. The data are combined across different spectroscopic and characterization methods and presented in the context of periodic trends. When considering crystallographic data, solution chemistry results, and the nuclear properties of actinium isotopes, it appears that some structural parameters ascribed to the Ac 3+ ion may have been overestimated. This review can guide researchers interested in actinide sciences and those who are pursuing the development of actinium-based radiotherapies, from isotope production to clinical trials.
Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.
This paper provides an experimental and theoretical investigation of azimuthal instabilities in a lean, premixed combustor. It presents data from a high Reynolds number, industrial-scale facility with multiple fuel stages that enable variation of radial and azimuthal symmetries. The combustor exhibits a variety of behaviors – including standing, spinning, mixed waves, and intermittent behaviors, as well as periodic limit cycles and quasi-periodic oscillations – depending upon varying power and azimuthal symmetry. A second key objective of this paper is to present a weakly nonlinear framework that is used to understand and interpret these results, showing the influence of eigenmode degeneracy, background noise, and azimuthal non-uniformities in mean quantities such as temperature and flame-acoustic coupling. In particular, this analysis shows the importance of the linear instability growth/damping rate and azimuthal non-uniformities in controlling the relative prominence of standing, spinning, or mixed waves at the limit cycle, as well as the preference in the spinning direction. It also shows the conditions under which single frequency or quasi-periodic oscillations are dominant, and how these are more fundamentally influenced by azimuthal non-uniformities and background noise.
Optimal management and desalination of produced water is a major challenge for U.S. oil and gas development. Integrating rigorous desalination models into multi-period produced water network optimization problems presents several hurdles, which need to be tackled using advanced optimization strategies. Here, in this work, a novel multi-period produced water network formulation with separate solid and liquid flows is introduced to avoid singularities at zero flows. Rigorous steady state desalination models based on mechanical vapor recompression are embedded at the desalination sites in the network model. An integrated optimization formulation is developed to co-optimize the design of desalination units along with the operation of the network. Furthermore, a more robust approach based on the trust region filter method is developed to efficiently integrate complex desalination models into the multi-period planning problem. Both optimization approaches are demonstrated on a produced water network from the PARETO library (Drouven et al., 2022) using thermal desalination units. Our results show that while the TRF and integrated approaches have comparable solve times, the TRF approach has better performance reliability in terms of solver convergence. Furthermore, the optimal solution obtained by embedding rigorous models into the network is significantly different than when desalination costs are approximated using simple cost models, which motivates further research in this field.
Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.
The corrosion behavior of an Al 0.1 CoCrFeNi HEA in high temperature CO 2 at 850 °C for different exposure periods was investigated using different characterization methods. Further, weight gain increased with extended exposure periods, yet it remained lower than most conventional alloys. A double-layer oxide film consisting of Cr 2 O 3 and Al 2 O 3 formed after 100 and 500 h of exposure, while a single layer of Cr 2 O 3 developed after 1000 h of exposure. CoCrFeNi particles with an FCC structure were embedded in the oxide films irrespective of the exposure periods. Moreover, a transition zone containing numerous striped-like Al 2 O 3 oxides along the grain boundaries was developed underneath the oxide films.