Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial heterogeneity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Continental United States may lose 1.8 petagrams of soil organic carbon under climate change by 2100

Abstract Aims High‐resolution information on soils’ vulnerability to climate‐induced soil organic carbon (SOC) loss can enable environmental scientists, land managers, and policy makers to develop targeted mitigation strategies. This study aims to estimate baseline and decadal changes in continental US surface SOC stocks under future emission scenarios. Location Continental United States. Time period 2014–2100. Methods We used recent SOC field observations ( n = 6,213 sites), environmental factors ( n = 32), and an ensemble machine learning (ML) approach to estimate baseline SOC stocks in surface soils across the continental United States at 100‐m spatial resolution, and decadal changes under the projected climate scenarios of Coupled Model Intercomparison Project Phase Six (CMIP6) earth system models (ESMs). Results Baseline SOC projections from ML approaches captured more than 50% of variability in SOC observations, whereas ESMs represented only 6–16% of observed SOC variability. ML estimates showed a mean total loss of 1.8 Pg C from US surface soils under the high‐emission scenario by 2100, whereas ESMs showed no significant change in SOC stocks with wide variation among ESMs. Both ML and ESM predictions agree on the direction of SOC change (net emissions or sequestration) across 46–51% of continental US land area. These differences are attributable to the high‐resolution site‐specific data used in the ML models compared to the relatively coarse grid represented in CMIP6 ESMs. Main conclusions Our high‐resolution estimates of baseline SOC stocks, identification of key environmental controllers, and projection of SOC changes from US land cover types under future climate scenarios suggest the need for high‐resolution simulations of SOC in ESMs to represent the heterogeneity of SOC. We found that the SOC change is sensitive to key soil related factors (e.g. soil drainage and soil order) that have not been historically considered as input parameters in ESMs, because currently more than 95% variability in the SOC of CMIP6 ESMs is controlled by net primary productivity, temperature, and precipitation. Using additional environmental factors to estimate the baseline SOC stocks and predict the future trajectory of SOC change can provide more accurate results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessing Multi-Dimensional Impacts of Achieving Sustainability Goals by Projecting the Sustainable Agriculture Matrix Into the Future

The concept of sustainability inherently spans multiple spatial scales, sectors, variables, and time horizons. This study links a recently developed method of assessing present-day agricultural sustainability across environmental, economic, and social dimensions with a process-based integrated assessment model, in order to allow forward-looking analysis of sustainability by region and scenario. The sustainable agriculture matrix estimates present-day agricultural sustainability at the national level using 18 indicator variables, of which this study estimates nine to the year 2100, using an enhanced version of the Global Change Analysis Model. Scenarios include a reference scenario, and scenarios that apply the following measures, both individually and in combination, that are thought to improve sustainability: yield intensification, transition toward more plant-based (“flexitarian”) diets, and economy-wide greenhouse gas emissions mitigation. The scenarios illustrate considerable complexity and tradeoffs inherent to efforts to improve agricultural sustainability in all regions globally. For example, yield intensification typically increases nitrogen pollution, flexitarian diets can reduce agricultural output, and greenhouse gas mitigation efforts may either increase deforestation or crowd out crop and livestock production due to consequent bioenergy demands. However, there is considerable inter-regional heterogeneity in the responses, and the importance of such secondary responses also differs by region. The analysis and post-processing methods developed in this study allow quantification and visualization of the absolute and relative magnitude of the tradeoffs between agricultural sustainability indicator variables across regions, time periods, and scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 2: Stratospheric and free-tropospheric response

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

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

Macroscale property assessment and indentation characteristics of thick section friction stir welded AA 5083

The spatial gradients in microstructure have been proven to be effective in optimizing the performance of components, and several techniques have been developed to achieve this microstructure. Nevertheless, the local microstructural modification methods inevitably result in the formation of gradient microstructure leading to a variation in the property. Herein this work systematically investigates the effect of gradient microstructure formed during thick section friction stir weld Al5083 alloy on its mechanical properties. A novel macro-indentation technique, profilometer-based indentation plastometry (PIP), is utilized to capture the variation in the mechanical properties in the feedstock (rolled plate) and weld region. The feedstock exhibited a gradient microstructure owing to different stress states during rolling. Additionally, the stirring of the tool resulted in a complex material flow pattern producing a microstructure gradient in the weld region. The aspect ratio of the grains was identified as the primary factor causing variation in indentation response in the feedstock. In contrast, variation in the grain size in the weld region resulted in property variation. Overall, the stir zone has a lower yield strength (158±4MPa) than the base material (173±3MPa). In addition, microstructure heterogeneity in the stir zone resulted in lower elongation than the base material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Patterns and controls on island‐wide aboveground biomass accumulation in second‐growth forests of Puerto Rico

Abstract Understanding the heterogeneity of biomass accumulation in second‐growth tropical forests following land use abandonment is important for informing ecosystem carbon models and forest restoration efforts. There is an urgent need for a broad sample of second‐growth forests to enhance our knowledge of carbon accumulation in human‐dominated landscapes, especially for older forests. Puerto Rico has predominantly second‐growth forests, ranging in age from approximately 25 to more than 80 years. We used an island‐wide sample of airborne lidar from the NASA Goddard Lidar, Hyperspectral, and Thermal (G‐LiHT) Airborne Imager collected on March 2017, forest inventory data, and data on forest age, precipitation, soils, and land use to estimate aboveground biomass stocks in moist and wet, second‐growth tropical forests. Biomass accumulation rates in Puerto Rico were lower, on average, than in other Neotropical forests. Median biomass across >16,700 ha of older second‐growth forests was 105 Mg ha −1 , and sampled biomass rarely surpassed 250 Mg ha −1 . Differences in biomass by age were large and persistent across different substrates and land uses, with a plateau in the pattern of island‐wide biomass accumulation after about 33 years. A spatial regression model showed that multiple factors were related to biomass accumulation, including time since abandonment, geologic substrate, past land use as coffee or pasture, precipitation, topographic wetness index, and slope. Our findings have important consequences for the total carbon storage and expected climate mitigation benefits of large‐scale reforestation efforts, and highlight the value of airborne lidar for quantifying biomass variability in complex tropical landscapes. Abstract in Spanish is available with online material.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Preliminary Study on TRISO Fuel Cross Section Generation

Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Molecular-Scale Characterization of Photoinduced Charge Separation in Mixed-Dimensional InSe–Organic van der Waals Heterostructures

Layered indium selenide (InSe) is an emerging two-dimensional semiconductor that has shown significant promise for high-performance transistors and photodetectors. The range of optoelectronic applications for InSe can potentially be broadened by forming mixed-dimensional van der Waals heterostructures with zero-dimensional molecular systems that are widely employed in organic electronics and photovoltaics. Here, we report the spatially resolved investigation of photoinduced charge separation between InSe and two molecules (C 70 and C 8 -BTBT) using scanning tunneling microscopy combined with laser illumination. We experimentally and computationally show that InSe forms type-II and type-I heterojunctions with C 70 and C 8 -BTBT, respectively, due to an interplay of charge transfer and dielectric screening at the interface. Laser-excited scanning tunneling spectroscopy reveals a ~0.25 eV decrease in the energy of the lowest unoccupied molecular orbital of C 70 with optical illumination. Furthermore, photoluminescence spectroscopy and Kelvin probe force microscopy indicate that electron transfer from InSe to C 70 in the type-II heterojunction induces a photovoltage that quantitatively matches the observed downshift in the tunneling spectra. In contrast, no significant changes are observed upon optical illumination in the type-I heterojunction formed between InSe and C 8 -BTBT. Density functional theory calculations further show that, despite the weak coupling between the molecular species and InSe, the band alignment of these mixed-dimensional heterostructures strongly differs from the one suggested by the ionization potential and electronic affinities of the isolated components. Self-energy-corrected density functional theory indicates that these effects are the result of the combination of charge redistribution at the interface and heterogeneous dielectric screening of the electron–electron interactions in the heterostructure. In addition to providing specific insight for mixed-dimensional InSe–organic van der Waals heterostructures, this work establishes a general experimental methodology for studying localized charge transfer at the molecular scale that is applicable to other photoactive nanoscale systems.

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

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2021

Reduction and Agglomeration of Supported Metal Clusters Induced by High-Flux X-ray Absorption Spectroscopy Measurements

Supported metal clusters are widely used in catalysis for many important reactions. To understand the catalytic properties, in situ/operando characterization techniques, such as X-ray absorption spectroscopy (XAS), provide essential details of the size, shape, and chemical composition of the cluster and the nature of the active sites. New-generation synchrotrons combined with focusing beamlines provide high-flux-density X-rays for improved detection sensitivity as well as higher time and spatial resolution. Understanding the effects of a high-flux-density X-ray beam on the catalyst during the actual measurement, whether XAS or another synchrotron-based technique, is crucial. This is especially important for in situ and operando studies where both the high flux density and reaction conditions can affect the catalyst structure. In this work, we investigated the effect of the flux density on rhodium clusters supported on Al 2 O 3 at two different beamlines: National Synchrotron Light Source II beamline 08-ID and Stanford Synchrotron Radiation Light Source (SSRL) beamline 4-1. We show that the higher flux density at beamline 08-ID causes the reduction of the highly dispersed RhO x /Al 2 O 3 catalyst, even at room temperature. Additionally, exposure to the higher flux density X-rays at beamline 08-ID during in situ reduction results in significant agglomeration of the Rh clusters. The final size of the Rh nanoparticles reduced at 310 °C is equivalent to that of particles formed after the reduction at 600–650 °C in the absence of the beam. Significant beam-induced reduction and agglomeration is also shown for Ni supported on beta zeolite during in situ reduction at an intermediate-flux-density beamline 9-3 at SSRL, indicating that beam-induced changes in heterogeneous catalysts could be common at intermediate- and high-flux-density beamlines. We provide precautions and recommendations for detecting and minimizing beam damage during in situ/operando XAS measurements.

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

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

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

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
DOE OSTI2021

Spatially-resolved lithiation dynamics from operando X-ray diffraction and electrochemical modeling of lithium-ion cells

In this study, energy dispersive X-ray diffraction is used to profile the time evolution of ordered Li x C 6 phases in solid electrodes of lithium-ion cells charged at rates between 0.2 and 4.7C (where 1C corresponds to full discharge in 1 h). The methods for quantifying lithium concentration in these phases from the acquired diffraction patterns are described. Compact expressions for time-dependent concentration gradients in the solid electrodes using orthogonal polynomial expansions are presented. Experimentally, these gradients persisted in lithiated graphite electrodes even after the cells rested at open-circuit for over 9 h. A multiphase electrochemical model of graphite intercalation captured many of the observed behaviors, including the progression of phase transitions and the persistent gradients at zero current. However, the magnitude of concentration gradients in both the oxide cathode and graphite anode is underestimated by the model, even at moderate currents.

25 ENERGY STORAGE↗
DOE OSTI2023

Highly Resolved Projections of Passenger Electric Vehicle Charging Loads for the Contiguous United States: Results From and Methods Behind Bottom-Up Simulations of County-Specific Household Electric Vehicle Charging Load (Hourly 8760) Profiles Projected Through 2050 for Differentiated Household and Vehicle Types

This report documents enhancements made to the TEMPO (Transportation Energy & Mobility Pathway Options TM ) model to project spatially, demographically, and temporally resolved national-scale EV charging load profiles and describes three scenarios and corresponding datasets created for the NREL demand-side grid (dsgrid) project in support of bulk power systems modeling. In brief, TEMPO was enhanced to disaggregate national and annual energy demand projections into household and county-level projections of passenger electric vehicle (EV) hourly charging load profiles (8760 profiles), accounting for consumer, travel, and temperature variations that impact EV energy demand. In alignment with NREL's forward-looking grid modeling, three scenarios for EV adoption covering 2020-2050 were created-- Annual Energy Outlook (AEO) Reference Case, Electrification Futures Study (EFS) High Electrification, and All EV Sales by 2035 --and associated datasets have been included in the dsgrid platform for public use.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

High methane concentrations in tidal salt marsh soils: Where does the methane go?

Abstract Tidal salt marshes produce and emit CH 4 . Therefore, it is critical to understand the biogeochemical controls that regulate CH 4 spatial and temporal dynamics in wetlands. The prevailing paradigm assumes that acetoclastic methanogenesis is the dominant pathway for CH 4 production, and higher salinity concentrations inhibit CH 4 production in salt marshes. Recent evidence shows that CH 4 is produced within salt marshes via methylotrophic methanogenesis, a process not inhibited by sulfate reduction. To further explore this conundrum, we performed measurements of soil–atmosphere CH 4 and CO 2 fluxes coupled with depth profiles of soil CH 4 and CO 2 pore water gas concentrations, stable and radioisotopes, pore water chemistry, and microbial community composition to assess CH 4 production and fate within a temperate tidal salt marsh. We found unexpectedly high CH 4 concentrations up to 145,000 μmol mol −1 positively correlated with S 2− (salinity range: 6.6–14.5 ppt). Despite large CH 4 production within the soil, soil–atmosphere CH 4 fluxes were low but with higher emissions and extreme variability during plant senescence (84.3 ± 684.4 nmol m −2 s −1 ). CH 4 and CO 2 within the soil pore water were produced from young carbon, with most Δ 14 C‐CH 4 and Δ 14 C‐CO 2 values at or above modern. We found evidence that CH 4 within soils was produced by methylotrophic and hydrogenotrophic methanogenesis. Several pathways exist after CH 4 is produced, including diffusion into the atmosphere, CH 4 oxidation, and lateral export to adjacent tidal creeks; the latter being the most likely dominant flux. Our findings demonstrate that CH 4 production and fluxes are biogeochemically heterogeneous, with multiple processes and pathways that can co‐occur and vary in importance over the year. This study highlights the potential for high CH 4 production, the need to understand the underlying biogeochemical controls, and the challenges of evaluating CH 4 budgets and blue carbon in salt marshes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spectral Induced Polarization-Biogeochemical Relationships for Remediation Amendment Monitoring

Geophysical tools such as electrical resistivity (ER) can indirectly monitor subsurface changes in response to remedial injections. These methods exhibit relatively low spatial resolution compared to sediment core characterization but are advantageous due to the ability to collect measurements non-intrusively over time across large volumes of the subsurface. Moreover, along with confirmatory groundwater or core sampling, geophysical methods can be used during active biogeochemical remedies to monitor short-term contaminant transformations and mobility, as well as part of an overall strategy for long-term monitoring of subsurface contaminated sites. The use of alternating current spectral induced polarization (SIP) provides significantly more information than conventional geophysical methods like direct current ER. SIP allows for monitoring changes in both solution and surface conductivity by separation of real and imaginary conductivity, respectively, as well as surface capacitance. In principle, SIP can measure indicators of remedy progression such as precipitation reactions that sequester contaminants, potentially providing a better indication of amendment delivery and reactivity as compared to conventional ER methods. However, multiple processes and material properties have overlapping (interacting) electrical responses across a range of frequencies. Hence, the purpose of this scoping study was to evaluate if SIP can be used to monitor (a) amendment delivery and (b) precipitation and reactivity of amendments under consideration at Hanford. SIP measurements were collected in fully saturated columns packed with sand and Hanford formation sediments containing (a) amendments that were highly conductive with significant capacitance (zero valent iron – ZVI, sulfur modified iron – SMI) and (b) amendments that exhibited low electrical conductivity with a small capacitance (calcite, apatite, bismuth). The sand was a quartz material with homogenous particle size that exhibited a relatively low surface conductivity. It was used as a control for comparison with the sediments from the Hanford Site, which have a greater surface conductivity because of their complex mineralogy and heterogeneous particle size distribution and may have complex interactions with amendments. The amendment mass fraction was varied to represent the different stages and subsurface locations associated with the delivery of a remedy. The primary objective was to identify the solution and solid surface changes associated with the delivery amendments and their secondary reactions within the subsurface. The figure below summarizes results for the amendments tested through this project. The SIP phase shift or imaginary conductivity change for the high conductivity amendments was more than 10 times that of the low conductivity amendments, highlighting the relative ease of detection of ZVI and SMI independent of the background signal from sand or Hanford sediments. The ZVI phase shift and imaginary conductivity changes occur primarily at high frequency (> 100 Hz) whereas SMI changes were at low frequency (0.01 to 10 Hz). SIP signals of SMI also increased over time and the maximum shifted to lower frequencies. The low conductivity amendments (calcite, abiotic and biotic apatite, bismuth subnitrate) exhibited relatively small phase and imaginary conductivity changes when added to sediments. The changes were above minimum detection limits (0.5 mrad for phase shift, 0.03 µS/cm for imaginary conductivity) for the highest concentration except for the commercial bismuth material. The lowest amendment concentration that can be detected is likely sediment specific, as minerals in sediments (clays, magnetite, Fe-oxides) have some capacitance and, therefore, exhibit a variable phase shift.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗