Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “micro scale models”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Saddle point approaches in strong field physics and generation of attosecond pulses

Attoscience is the emerging field that accesses the fastest electronic processes occurring at the atomic and molecular length scales with attosecond (1 as = 10 –18 s) time resolution having wide ranging physical, chemical, material science and biological applications. The quintessential and one of the most fundamental processes in this domain is the generation of phase locked XUV attosecond pulses. The theoretical approach to understand the process incorporates a fully quantum or semi classical or relativistic description of coherent charge dynamics in intense ultrashort electromagnetic fields driving a quantum system (an atom, a molecule, solid band gap materials or surface plasmas). Modelling of such physical and dynamical systems in science and also in many other branches often leads to equations represented in terms of complex multi-dimensional integrals. These integrals can often be solved using the stationary phase approximation, which leads to a series of equations identifying the points in the multi-dimensional space, having most significant contributions in their evaluation. These points are usually indicated as saddle points. The description of the dynamics of quantum mechanical or relativistic systems that results from such an approach enables near to classical physics intuitive perceptions of the processes under investigation. Thus, the saddle point methods are very powerful and valuable general theoretical tools to obtain asymptotic expressions of such solutions and help also to gain physical insights on the underlying phenomena. Such techniques developed in the past have been adapted to study the emission of as pulses by different physical systems and have been widely employed in calculating and estimating the response of matter to intense electromagnetic pulses on ultrafast time scales. Here we provide an extensive disposition of the saddle point approaches unifying their ubiquitous applications within the domain of attoscience valid for simple atomic to more complex condensed matter systems undergoing ultrafast dynamics and present current trends and advancements in the field. In this review we would delineate the methodology, present a synthesis of seminal works and describe the state of the art applications. Finally we also address ultrashort time dynamics of novel materials that have gained much attention recently, namely lower dimensional material systems and micro-plasma systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

HydraGNN v5.0

HydraGNN v5.0 expands the code base into a more portable, scalable, and flexible framework for scientific graph learning, with particular strength in atomistic machine-learning interatomic potentials and large-scale distributed training. The release adds Fully Sharded Data Parallel (FSDP) support alongside existing DDP and DeepSpeed paths, including FSDP-aware checkpointing and optimizer integration, and introduces a configurable multi-precision training workflow supporting FP32, BF16, and FP64 across GPUs and Intel XPUs. For atomistic modeling, HydraGNN v5.0 strengthens its MLIP capabilities through dynamic graph construction at every forward pass, energy-conserving force prediction via automatic differentiation, and per-atom energy loss formulations, while extending EGNN models to properly handle periodic boundary conditions. The release also broadens model expressiveness through graph-level attribute conditioning, adds new multi-task and model-parallel extensions such as MACE support and encoder/decoder branch optimization, and expands application coverage with integrated examples for datasets including OC25, Nabla2-DFT, QCML, Open Polymers 2026, and OPF. In parallel, HydraGNN v5.0 improves production readiness through performance optimizations for large-scale runs, stratified sampling and linear-regression preprocessing utilities, and tested installation scripts for DOE supercomputers including Frontier, Aurora, Perlmutter, and Andes. Overall, the release advances HydraGNN as a robust software platform for scalable graph neural networks across materials science, chemistry, and scientific machine learning workflows

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2024

The Effect of Interlayer Delay on the Heat Accumulation, Microstructures, and Properties in Laser Hot Wire Directed Energy Deposition of Ti-6Al-4V Single-Wall

Laser hot wire directed energy deposition (LHW-DED) is a layer-by-layer additive manufacturing technique that permits the fabrication of large-scale Ti-6Al-4V (Ti64) components with a high deposition rate and has gained traction in the aerospace sector in recent years. However, one of the major challenges in LHW-DED Ti64 is heat accumulation, which affects the part quality, microstructure, and properties of as-built specimens. These issues require a comprehensive understanding of the layerwise heat-accumulation-driven process–structure–property relationship in as-deposited samples. In this study, a systematic investigation was performed by fabricating three Ti-6Al-4V single-wall specimens with distinct interlayer delays, i.e., 0, 120, and 300 s. The real-time acquisition of high-fidelity thermal data and high-resolution melt pool images were utilized to demonstrate a direct correlation between layerwise heat accumulation and melt pool dimensions. The results revealed that the maximum heat buildup temperature of the topmost layer decreased from 660 °C to 263 °C with an increase to a 300 s interlayer delay, allowing for better control of the melt pool dimensions, which then resulted in improved part accuracy. Furthermore, the investigation of the location-specific composition, microstructure, and mechanical properties demonstrated that heat buildup resulted in the coarsening of microstructures and, consequently, the reduction of micro-hardness with increasing height. Extending the delay by 120 s resulted in a 5% improvement in the mechanical properties, including an increase in the yield strength from 817 MPa to 859 MPa and the ultimate tensile strength from 914 MPa to 959 MPa. Cooling rates estimated at 900 °C using a one-dimensional thermal model based on a numerical method allowed us to establish the process–structure–property relationship for the wall specimens. The study provides deeper insight into the effect of heat buildup in LHW-DED and serves as a guide for tailoring the properties of as-deposited specimens by regulating interlayer delay.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2025

Secondary electron emission for reticulated carbon foam surfaces using direct measurements and spectroscopic analysis

This study investigates secondary electron emission (SEE) characteristics of reticulated foams using direct measurements and analytical modeling. Total SEE was quantified, revealing suppression of up to 44% in carbon foam structures compared to planar graphite surfaces. An optimal geometric configuration was identified and supported by analytical models. SEE angular dependence experiments showed diverse behaviors: fiber-like behavior and directional dependence for pore and ligaments on the mm scale, with fuzz-like characteristics when the foam features are between 10–100 µm. Electron energy analyzer measurements showed that carbon foams preferentially suppress inelastic backscattered electrons (BSEs) more so than true secondary electrons (SEs). The analysis indicated a larger fraction of low-energy SE generation in foams compared to flat surfaces due to increased emission from curved fiber ligaments and tertiary SEs from high-energy BSEs. These findings have implications for design and optimization of materials with tailored electron emission properties for applications like plasma-facing components, spacecraft materials, and accelerator surfaces.

Auger↗
DOE OSTI2020

GIS Visualization of Transportation Energy Consumption

Transportation is witnessing unprecedented transformation, where emerging technologies are disrupting the way we travel. Be it sharing economy, or micro-mobility, the landscape of urban transportation is undergoing a much-needed paradigm shift. In order to capture and model these shifts, researchers need to be agile in their studies of the current and future transportation landscape. The need for agility in turn calls for sophisticated, intuitive, and reliable means to visualize and share transportation data. Travel is inherently spatial, so methods of visualizing transportation information should also be rooted in spatial analysis. Geographic Information Systems (GIS) are the premier technological framework to analyze and display spatial data. GIS tools have been used in the past to depict flow of vehicles, and display of network conditions (speed, congestion, etc.). However, in an ever-changing technological landscape, spanning advancements in vehicle as well as information systems, depicting vehicle flows falls short of providing a comprehensive picture of the impact these advancements have on travel related energy consumption. To address this issue, this research effort presents a web-based mapping application to visualize how travel related energy flows across a city. The application is being developed by National Renewable Energy Laboratory researchers to integrate various transportation energy consumption models within a GIS schema. This has the goal of enabling rapid analysis of energy impacts of the dramatically evolving transportation environment. The application is being developed using a Python framework for ArcGIS Online. Using road network data from the City of Columbus, Ohio and traffic data from CATT Laboratory's Regional Integrated Transportation Information System, transportation energy consumption will be modeled at a macro level. The application aims to provide highly accurate data while still maintaining the flexibility needed to adapt the model when new technologies arise. The methodology presented through this effort is expected to provide insights into how light-duty vehicles use energy on a large scale based on the road network they use. The application will also provide data exports to enable sharing and collaboration. " need for agility in turn calls for sophisticated, intuitive, and reliable means to visualize and share transportation data. Travel is inherently spatial, so methods of visualizing transportation information should also be rooted in spatial analysis. Geographic Information Systems (GIS) are the premier technological framework to analyze and display spatial data. GIS tools have been used in the past to depict flow of vehicles, and display of network conditions (speed, congestion, etc.). However, in an ever-changing technological landscape, spanning advancements in vehicle as well as information systems, depicting vehicle flows falls short of providing a comprehensive picture of the impact these advancements have on travel related energy consumption. To address this issue, this research effort presents a web-based mapping application to visualize how travel related energy flows across a city. The application is being developed by National Renewable Energy Laboratory researchers to integrate various transportation energy consumption models within a GIS schema. This has the goal of enabling rapid analysis of energy impacts of the dramatically evolving transportation environment. The application is being developed using a Python framework for ArcGIS Online. Using road network data from the City of Columbus, Ohio and traffic data from CATT Laboratory's Regional Integrated Transportation Information System, transportation energy consumption will be modeled at a macro level. The application aims to provide highly accurate data while still maintaining the flexibility needed to adapt the model when new technologies arise. The methodology presented through this effort is expected to provide insights into how light-duty vehicles use energy on a large scale based on the road network they use. The application will also provide data exports to enable sharing and collaboration.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Spatially resolved X-ray imaging and molecular characterisation of sulfur and iron in organic- and sulfur-rich hydrocarbon source mudstones

Sulfurisation of organic matter (OM) is a prominent preservation mechanism, however iron sulfide precipitation, particularly pyrite (FeS2), can counteract this mechanism. There is a dearth of high-resolution, spatially-resolved spectroscopic (redox) information on sulfur and iron inventories within organic-rich rocks that would improve our understanding of prevailing environmental conditions during deposition. Here, state-of-the-art synchrotron-based X-ray absorption and fluorescence analyses of key organic- and sulfur-rich mudstones demonstrate the potential of these techniques to non-destructively map and produce detailed spectroscopic information. Detailed high-resolution analyses (μm- to mm-scale) reveal the presence of widespread sulfurised OM in the Blackstone Band of the Kimmeridge Clay Formation, in line with a persistence of euxinia over a long temporal span and low reactive iron input, facilitating the preservation of OM through sulfurisation. In contrast, the presence of sulfurised OM was transitional in the Monterey Formation, consistent with fluctuating water column redox conditions, and is less significant in the Whitby Mudstone Formation, likely due to the high reactive iron concentrations outcompeting sulfurised OM formation. Analyses of sulfur species using model compounds further indicate that the Whitby Formation is strongly enriched in inorganic reduced sulfur minerals, while both the Kimmeridge Clay and Monterey Formations are dominated by organic sulfur species. These synchrotron-based observations improve our understanding of environmental conditions during the time of deposition of these mudstones and thus show great promise in the study of organic-rich sediments, especially in allowing their depositional settings to be more accurately reconstructed.

Kimmeridge Clay Formation↗
DOE OSTI2025

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗
DOE OSTI2021

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗
DOE OSTI2019

EGS Collab Earth Modeling: Integrated 3D Model of the Testbed

The EGS Collab project is conducting a series of stimulation and interwell flow tests in an intermediate scale (~10-20 m) testbed located on the 4850 level in the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the Poorman Formation from the western rib of the West Access Drift. Of the eight boreholes, one borehole is used as an injection/stimulation well, another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as monitoring wells. For 3D visualization of the testbed as well as to provide spatially accurate and consistent parameter data for geomechanical and geophysical process modeling tasks, we are using Leapfrog Software (Seequent Limited). For analysis of fractures and creation of conditioned stochastic discrete fracture network (CS-DFN) model of the testbed, we are using FracMan (Golder Associates Inc.). As of mid-May 2019, the EGS Collab Leapfrog database contains the general layout of the drifts/shafts at the 4850 and 4100 levels; generalized geologic framework model of a domain that spans both the 4850 and 4100 levels; the CS-DFN of the testbed; and layouts of the EGS Collab testbed boreholes, kISMET testbed boreholes, and other existing boreholes at the 4850 level. This information has been obtained from the SURF Vulcan database, detailed laser scan and point surveys of borehole collar locations, directional surveys of borehole trajectories, and mapping of features in the drifts, boreholes, and cores. Similarly, other types of data that have been incorporated into the Leapfrog database include static features of the testbed such as natural discontinuities (e.g., fractures, weeps), location of monitoring instruments, and stimulation notches. Baseline geophysical characterization results (e.g., seismic, ERT, temperature, etc.) are also imported into the Leapfrog database. The location and size of packers, various point-source signals detected during stimulation and flow tests (e.g., temperature anomalies, micro-earthquake events, etc.) are also compiled in the database. All feature locations are reported using the local Homestake Mine coordinate system. In this paper, we present some of the geostructural characteristics as well as the status of the testbed as represented by Leapfrog/FracMan database/visualizations.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

NGEE Arctic Soil Micro-warming Experiment Temperature Profiles, Council Road Mile Marker 71, Seward Peninsula, Alaska, 2017-2019

The reported soil temperature profile measurements (24 locations - 2 probes with 5 and 6 thermistors, respectively - Fig. 1 and Fig. 2) were initiated as part of a soil micro-warming experiment at the Council Road Mile Marker 71 Site (CN_MM71) in September 2017. From 2017 through August 2019, these were measurements of ambient pre-treatment plot conditions. This tussock tundra site (with underlying permafrost) experiences annual frost heaving that causes a vertical movement of the ground surface, hence causing some of the upper most temperature thermistors to be at surface level or above ground and recorded therefore air temperatures near the surface and not below ground soil temperature (see section on Quality Assurance). After August 2019, these measurements ended in preparation for transitioning to the experimental warming of individual plots. The reported temperature data are nominally 3-hour averages of measurements made at varying frequencies (3-hour maximum) with frequencies that depended upon power (solar) availability to operate the sensors (see section in documentation on Methods). The number of values and the standard deviation for each 3-hour average (where appropriate) are also provided. These measurements are located near the NGEE-Arctic CO2 and CH4 eddy covariance tower at the Council Road Site (US-NGC: NGEE Arctic Council, https://ameriflux.lbl.gov/sites/siteinfo/US-NGC#overview) and auxiliary data at https://doi.org/10.5440/1526749. This dataset contains four *.csv files and one *.pdf file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗
DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Association between soil organic carbon and calcium in acidic grassland soils from Point Reyes National Seashore, CA

Organo-mineral and organo-metal associations play an important role in the retention and accumulation of soil organic carbon (SOC). Recent studies have demonstrated a positive correlation between calcium (Ca) and SOC content in a range of soil types. However, most of these studies have focused on soils that contain calcium carbonate (pH > 6). To assess the importance of Ca-SOC associations in lower pH soils, we investigated their physical and chemical interaction in the grassland soils of Point Reyes National Seashore (CA, USA) at a range of spatial scales. Multivariate analyses of our bulk soil characterisation dataset showed a strong correlation between exchangeable Ca (Ca Exch ; 5–8.3 c.mol c kg –1 ) and SOC (0.6–4%) content. Additionally, linear combination fitting (LCF) of bulk Ca K-edge X-ray absorption near-edge structure (XANES) spectra revealed that Ca was predominantly associated with organic carbon across all samples. Scanning transmission X-ray microscopy near-edge X-ray absorption fine structure spectroscopy (STXM C/Ca NEXAFS) showed that Ca had a strong spatial correlation with C at the microscale. The STXM C NEXAFS K-edge spectra indicated that SOC had a higher abundance of aromatic/olefinic and phenolic C functional groups when associated with Ca, relative to C associated with Fe. In regions of high Ca-C association, the STXM C NEXAFS spectra were similar to the spectrum from lignin, with moderate changes in peak intensities and positions that are consistent with oxidative C transformation. Through this association, Ca thus seems to be preferentially associated with plant-like organic matter that has undergone some oxidative transformation, at depth in acidic grassland soils of California. Our study highlights the importance of Ca-SOC complexation in acidic grassland soils and provides a conceptual model of its contribution to SOC preservation, a research area that has previously been unexplored.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗