Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difference scaling”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Towards libraries of different ultrasmall amorphous silica cage topologies

Cage topology controlled at the nanometer length scale is expected to enable original functionalities, e.g., in catalysis and therapeutics. Cages are fundamental structural motifs of clathrates and their topological duals, Frank-Kasper phases, and are constituents of mesoporous silica frameworks. However, with one exception, they have not been synthesized as discrete particles. By varying reactant ratios of surfactant, oil, and silane, we report the discovery of 5(12)6(2) and 5(12)6(8) amorphous silica polyhedra, each coexisting with the previously identified 5(12) cage, using combined cryo-transmission electron microscopy (cryo-TEM) and single-particle reconstruction (SPR). Notably, the 5(12)6(8) polyhedron represents a topology not previously observed in mesoporous silica frameworks. Control over structural features is demonstrated, and insights into cage formation mechanisms are provided. Structural outcomes are summarized in a ternary morphology diagram alongside prior results, bridging serendipitous discovery and intentional design of silica cages displaying a level of control over silica polymerization rivaling nature.

Jang, Dong June [Department of Materials Science a↗
DOE OSTI2023

Validating Mixtures of 233 U, 235 U, and 239 Pu for the Sum-of-Fractions Method

The Sum-of-Fractions method is a technique used to assure that homogeneous mixtures of fissile and fissionable isotopes are below a minimum margin of k eff or reactivity. Current work by Pacific Northwest National Laboratory examines different mixtures of 233 U, 235 U, and 239 Pu to determine critical mass limits for mixtures of transuranic actinides lacking a validation basis. To provide a validation basis for these limits, the work presented here describes the results of a sensitivity and uncertainty analysis of various mixtures of these isotopes in various concentrations moderated and reflected by light water and polyethylene. The TSUNAMI-1D sequence in the SCALE code system was used to generate sensitivity coefficients for three different concentrations of mixtures of 233 U, 235 U, and 239 Pu. The TSUNAMI-IP sequence was then used for similarity assessment (c k ) with critical benchmark experiment sensitivity data files (SDFs) from the Oak Ridge National Laboratory Verified, Archived Library of Inputs and Data and the Nuclear Energy Agency SDF database. The VADER sequence in SCALE was used for statistical testing and to generate upper subcritical limits from the data to develop a basis for validating critical mass limits.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Indexed improvements for real-time trotter evolution of a (1 + 1) field theory using NISQ quantum computers

Today's quantum computers offer the possibility of performing real-time calculations for quantum field theory scattering processes motivated by high energy physics. In order to follow the successful roadmap which has been established for the calculation of static properties at Euclidean time, it is crucial to develop new algorithmic methods to deal with the limitations of current noisy intermediate-scale quantum (NISQ) devices and to establish quantitative measures of the progress made with different devices. In this paper, we report recent progress in these directions. We show that nonlinear aspects of the trotter errors allow us to take much larger step then suggested by low-order analysis. This is crucial to reach physically relevant time scales with today's NISQ technology. We propose to use an index averaging absolute values of the difference between the accurately calculated trotter evolution of site occupations and their actual measurements on NISQ machines (G index) as a measure to compare results that have been obtained from different hardware platforms. Using the transverse Ising model in one spatial dimension with four sites we apply this metric across several hardware platforms. We study the results including readout mitigation and Richardson extrapolations and show that the mitigated measurements are very effective based on the analysis of the trotter step size modifications. Furthermore, we discuss how this advance in the trotter step size procedures can improve quantum computing physics scattering results and how this technical advance can be applied to other machines and noise mitigation methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Scaling patch analysis of planar turbulent mixing layers

Proper scales for the mean flow and Reynolds shear stress in planar turbulent mixing layers are determined from a scaling patch analysis of the mean continuity and momentum equations. By seeking an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is determined as V ref =(dδ/dx)U ref , where dδ/dx is the growth rate of the mixing layer width and U ref =U h –U l is the difference between the velocity of the high speed stream U h and the velocity of the low speed stream U l . By seeking an admissible scaling for the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is determined as R uv,ref =U avg V ref =[$\frac{1}{2A_u}$$\frac{dδ}{dx}$]$U^{2}_{ref}$ where A u $_{=}^{def}$(U h –U l )/(U h +U l ) is the normalized velocity difference that emerges naturally in the admissible scaling of the mean momentum equation. Self-similar equations for the scaled mean transverse flow V* and Reynolds shear stress $R^{*}_{uv}$=R uv /R uv,ref are derived from the mean continuity and mean momentum equations. Finally, approximate equations for V* and $R^{*}_{uv}$ are developed and found to agree well with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

High throughput single cell multiomics platform [Abstract]

In this collaborative project, PNNL and Scienion will co-develop an integrated microfluidic technology to co-measure the transcriptome and proteome in single cells. The technology will enable us to efficiently separate proteins from mRNA transcripts between two microchips, barcode the molecules, and measure them separately with next-generation sequencing and mass spectrometry, respectively. Multicellular organisms contain diverse cell types and tremendous cell-to-cell heterogeneity that dictates a multitude of biological functions in both physiological and pathological environments. Even in the case of microbes, these genetically identical organisms can randomly differentiate into many subpopulations that assume different roles for the survival of the community. Bulk-scale measurements are insufficient to resolve such complexities. The development and applications of high throughput single-cell omics technologies have transformed our understanding of cellular heterogeneities and their differential responses to internal signaling events or external stimulations. Despite these advances, most current single-cell omic technologies provide information on only one type of biomolecule. Perse, such measurements provide incomplete information because the cell phenotype is determined by multiple layers of biomolecules and the interplay between genome, epigenome, transcriptome, and proteome. For example, mRNA abundance in one cell can not be precisely referred to the corresponding DNA and protein in other cells because of the potential subtle difference in genotype (e.g., somatic mutation or copy number variation) or phenotype (external microenvironment and cell-cell interactions). As such, parallel measurement of multiple biomolecules in the same single cells can offer unique advantages compared with measuring them separately in different single cells. Scienion is a world-leading biotech company focusing on precision liquid handling and its application in single-cell whole-genome sequencing and RNA sequencing. PNNL is the leading institution in ultrasensitive mass spectrometry, microfluidics, and untargeted single-cell proteomics (scProteomics). This collaboration will facilitate a unique fusion between scTranscriptomics capability at Scienion and scProteomics capabilities at PNNL to, for the first time, perform both untargeted transcriptomics and proteomics from the same single cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mapping foliar photosynthetic capacity in sub-tropical and tropical forests with UAS-based imaging spectroscopy: Scaling from leaf to canopy

Accurate understanding of the variability in foliar physiological traits across landscapes is critical to improve parameterization and evaluation of terrestrial biosphere models (TBMs) that seek to represent the response of terrestrial ecosystems to a changing climate. Numerous studies suggest imaging spectroscopy can characterize foliar biochemical and morphological traits at the canopy scale, but there is only limited evidence for retrieving canopy photosynthetic capacity (e.g., maximum carboxylation rate, V c,max and maximum electron transport rate, J max ). Moreover, the effect of canopy structure within forest communities on scaling up spectra-trait relationships from leaf to canopy level is not well known. To advance the spectra-trait approach and enable the estimation of key traits using remote sensing, we collected imaging spectroscopy data from an Unoccupied Aerial System (UAS) platform over two forest sites in China (a subtropical forest in Mt. Dinghu and a tropical rainforest in Xishuangbanna). At these sites, we also collected ground measurements of leaf spectra and traits, including biochemical (leaf nitrogen, phosphorus, chlorophyll, and water content), morphological (leaf mass per area, LMA) and physiological (V c,max25 and J max25 ) traits (n=135 tree-crowns from 42 species across two sites). Using a partial least-squares regression (PLSR) approach, we built and tested spectra-trait models with repeated cross-validation. The spectral models developed with leaf spectra were directly transferred to canopy spectra to evaluate the effect of canopy structure. Here we further applied canopy spectral models to map these traits at individual tree-crown scale. The results demonstrate that (1) UAS-based canopy spectra can be used to estimate V c,max (R 2 =0.55, nRMSE=11.79%), Jmax (R 2 =0.54, nRMSE=12.34%), and five additional foliar traits (R 2 =0.38-0.60, nRMSE=10.11-13.56%) at the tree-crown scale with demonstrated generalizability across two sites; (2) canopy structure strongly affects the spectratrait relationships from leaf to canopy level, but the effects vary considerably across foliar traits and cannot be well captured by the 4SAIL canopy radiative transfer model. UAS-based imaging spectroscopy maps large variability in all foliar traits (including physiological traits) with spatially explicit information, reproducing the field-observed inter- and intra-specific variations. These results demonstrate the capability of using UAS-based imaging spectroscopy for characterizing the variability of foliar physiological traits at individual tree-crown scale over forest landscapes and highlight the similar generalizability but different biophysical mechanisms underlying spectra-trait relationships at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Atomically Resolved Electrically Active Intragrain Interfaces in Perovskite Semiconductors

Deciphering the atomic and electronic structures of interfaces is key to developing state-of-the-art perovskite semiconductors. However, conventional characterization techniques have limited previous studies mainly to grain-boundary interfaces, whereas the intragrain-interface microstructures and their electronic properties have been much less revealed. Herein using scanning transmission electron microscopy, we resolved the atomic-scale structural information on three prototypical intragrain interfaces, unraveling intriguing features clearly different from those from previous observations based on standalone films or nanomaterial samples. These intragrain interfaces include composition boundaries formed by heterogeneous ion distribution, stacking faults resulted from wrongly stacked crystal planes, and symmetrical twinning boundaries. The atomic-scale imaging of these intragrain interfaces enables us to build unequivocal models for the ab initio calculation of electronic properties. Our results suggest that these structure interfaces are generally electronically benign, whereas their dynamic interaction with point defects can still evoke detrimental effects. This work paves the way toward a more complete fundamental understanding of the microscopic structure-property-performance relationship in metal halide perovskites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quality Guidelines for Energy System Studies: Capital Cost Scaling Methodology: Revision 4b Report

The National Energy Technology Laboratory (NETL) regularly updates legacy analysis with new studies and cases as the Department of Energy objectives change, technology performance improves, costs are reduced, regulations change, market drivers are established, fuel prices fluctuate, and any number of other relevant factors vary in the market. As legacy studies are updated by NETL, the underlying performance and cost of the cases presented changes, and as such, the methods for interpreting and scaling the cost estimates change. Therefore, it is important that NETL maintain public guidance documents associated with different sets of cost estimates that delineate how a specific set, based on report vintage and/or year published, should be scaled. This Quality Guidelines for Energy System Studies report, providing guidance on capital cost scaling, should generally be applied to NETL case costs included in the report “Cost and Performance Baseline for Fossil Energy Plants Volume 1: Bituminous Coal and Natural Gas to Electricity Revision 4b,” (NETL, 2025) or any cases derived from the cases presented in the referenced report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Examining Ohmic Losses in Fuel-Cell Catalyst Layers with Different Pt/C Ratios

A model of a fuel-cell cathode with distinct proton conductivities at the agglomerate and electrode scales was used to simulate polarization curves under conditions with facile oxygen transport. Here, the model was applied to electrodes containing agglomerates with different ratios of platinum to carbon in order to examine the hypothesis that intra-agglomerate proton conductivity is substantially lower than ionomer conductivity. The model predicts decreasing performance and increasing sensitivity to Pt/C ratio as intra-agglomerate conductivity declines relative to a constant through-electrode conductivity. Experiments on catalyst layers with Vulcan supports do not show the predicted trend indicating that mesoscale intra-agglomerate conductivity is not orders of magnitude lower than bulk ionomer conductivity.

08 HYDROGEN↗
DOE OSTI2020

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗
DOE OSTI2022

Bacterial community dynamics explain carbon mineralization and assimilation in soils of different land-use history

Soil dwelling microorganisms are key players in the terrestrial carbon cycle, driving both the degradation and stabilization of soil organic matter. Bacterial community structure and function vary with respect to land use; yet the ecological drivers of this variation remain poorly described and difficult to predict. For this work, we conducted a multi-substrate DNA-stable isotope probing experiment across cropland, old-field, and forest habitats to link carbon mineralization dynamics with the dynamics of bacterial growth and carbon assimilation. We tracked the movement of 13 C derived from five distinct carbon sources as it was assimilated into bacterial DNA over time. We show that carbon mineralization, community composition, and carbon assimilation dynamics all differed with respect to land use. We also show that microbial community dynamics affect carbon assimilation dynamics and are associated with soil DNA content. Soil DNA yield is easy to measure and may be useful in predicting microbial community dynamics linked to soil carbon cycling. Soil dwelling microorganisms are key players in the terrestrial carbon cycle, driving both the degradation and stabilization of soil organic matter. Microbial communities vary with respect to land use, but we still have an incomplete understanding of how variation in community structure links to variation in community function. DNA stable isotope probing (DNA-SIP) is a high-resolution method that can identify specific microbial taxa that assimilate carbon in situ. We conducted a large-scale multi-substrate DNA-SIP experiment to explore differences in bacterial activity across land-use regimes. We show that microbial community dynamics vary with land use, that these dynamics are linked to soil carbon cycling, and that they are associated with easily measured soil properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Subtask 2.6 – Optimization of Aerosol Mitigation Technology for Postcombustion CO 2 Capture

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion CO 2 capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery (EOR). The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the aerosols present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and aerosols provide nucleation sites where amine vapors can condense. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, aerosol components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to optimize aerosol mitigation technology for PCCC at a lignite-fired power plant. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on collection of aerosols at a low-rank coal-fired power station. Determine the impact of aerosols on the efficiency and degradation products of amine-based carbon capture systems fired with low-rank fuels. Work was conducted at Minnkota Power Cooperative’s (MPC’s) Milton R. Young (MRY) Station Unit 2 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. To gather initial data for sizing and specifying a WESP for this system, a temporary pilot-scale WESP was rented and installed on-site. Several different conditions were tested to examine the impact of flow rate, voltage, and current on WESP performance. The WESP was effective at removing large particulate (>200 nm) but caused an increase in fine particulate (<75 nm). Fine particulate material at the inlet and outlet of the WESP was collected, analyzed, and showed that crystalline sulfates carried over from the plant’s FGD unit were being converted to fine aerosols and SO 2 was being converted to SO 3 through the WESP. Additionally, the high moisture content of the flue gas stream at this sample location also contributed to an overall increase in aerosol mass under some of the test conditions. Using the results from the rented WESP, the project team installed a smaller-scale WESP upstream of the EERC’s slipstream CO 2 capture system. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then routed through the new WESP before passing to the CO 2 absorber columns. Fluor’s amine-based solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Aerosols and sulfur species were measured at multiple locations throughout the pilot-scale system. The inlet FGD and DCC removed much of the particulate matter and gaseous sulfur upstream of the WESP. With this configuration, the WESP achieved >95% particulate capture. The new WESP did not show any of the increases in SO 3 or other aerosol species that had been consistently observed with the larger-scale WESP installed immediately downstream of the plant’s full-scale FGD unit. Particulate samples captured and analyzed from the WESP inlet did not show any presence of crystalline sulfate materials, indicating that the pilot-scale FGD and DCC were efficient at reducing carryover from the plant’s full-scale FGD unit. A set of parametric tests were conducted on the new WESP to assess the impacts of flow rate, voltage, number of online WESP fields, and gas-phase sulfur content on aerosol and sulfur transformations. The results showed that the WESP performed similarly well at all sets of conditions. Sulfur and particulate matter exiting the WESP were further reduced through the absorber column as the amine-based solvent captured some of the residual contaminants. Solvent analysis showed that these species were slowly concentrating in the solvent over the duration of the test. When the WESP was taken offline and the sulfur slip through the FGD allowed to rise, the sulfate content in the solvent rose sharply, showing that the extra FGD and WESP were effective at reducing sulfate and cation uptake. This would be expected to extend amine-based solvent life by slowing the formation of heat-stable salts and other degradation products. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission and MPC.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

AI-Driven Frameworks for Characterizing Urban Energy Systems

We develop AI-driven frameworks to characterize urban energy systems with the goal of transforming planning by reducing the labor of model generation, scaling scenario exploration, and improving accuracy for localized analysis. The approach integrates top-down and bottom-up data to train different AI models that predict missing information and generate inputs and targeted scenarios for district-scale models. The result is a scalable framework that provides actionable insights for reliable and efficient planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

The role of electron-scale turbulence in the JET tokamak: experiments and modelling

Dedicated electron heat transport experiments have been carried out in L- and H-mode Deuterium plasmas of the JET-ILW tokamak to identify the amount of electron heat carried by electron-scale electron temperature gradient (ETG) modes. Ion cyclotron resonance heating at different positions has been used to probe the response of the electron temperature inverse gradient length R/L Te to changes in electron heat flux qe, while different amounts of neutral beam heating allowed to scan the ratio of ion to electron temperature T e /T i , which is a key parameter for the onset of ETGs. Results indicate a steepening of the normalized qe vs R/L Te curve above R/LT e ~ 8 for T e /T i ≤ 1, suggestive of the ETG onset. Ion-scale gyro-kinetic (GK) simulations match the ion heat flux and the low-R/L Te part of the qe curve, but do not reproduce such steepening at high R/L Te . Multi-scale GK simulations covering both ion and electron scales and including one impurity bundling light and heavy species indicate an ETG contribution only for R/L Te values larger than the experimental ones. Sensitivity studies of such result are difficult to achieve due to limitation in numerical resources. The quasi-linear TGLF model has been used for sensitivity studies. With the same bundled impurity as the GK multi-scale, TGLF shows the qe steepening at much larger R/L Te values than in experiment, but when using the real mix of light impurities neglecting the heavy impurities, TGLF gets closer to the experimental results. Profile simulations with TGLF including both light and heavy impurities show over-prediction of T e profiles and in some cases also of density, but good T i predictions, confirming issues with the model electron stiffness for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

SEAS Communication Engine: An Extensible, Flexible Wrapper for Co-Simulation Agents

When modeling and analyzing the power grid and other large scale systems, researchers often express scenarios as optimization problems and feed them into advanced software solvers. In order to allow multiple solvers to communicate with each other and share data from different domains, the National Renewable Energy Laboratory (NREL) and associated Department of Energy (DOE) labs have developed a software framework called the Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS). HELICS allows cosimulation via a collection of client libraries for different languages that can be called from the appropriate optimization software. However, these client libraries do not provide a higher level of abstraction beyond reading and writing data off of the shared HELICS bus. In this paper, we describe a new software library called the SEAS Communication Engine that exposes a higher-level API for running cosimulation problems. The SEAS Engine provides a class-based abstraction on top of the Python HELICS client, in order to allow users to implement their domain-specific cosimulations without needing to interact with core HELICS primitives. This will make adoption of HELICS and cosimulation in general easier, by exposing a simpler API. In the second part of the paper, we validate our library on a collection of different simulation examples, including the canonical IEEE 13 Bus Feeder. Lastly, we demonstrate using the SEAS Engine to directly call domain-specific code written in the Julia programming language. Our hope is that this will serve as a template for easily calling software in different programming languages via the SEAS Engine, thereby avoiding code duplication and complexity.

co-simulation↗
DOE OSTI2024

A network approach for multiscale catchment classification using traits

Abstract. The classification of river catchments into groups with similar biophysical characteristics is useful to understand and predict their hydrological behavior. The increasing availability of remote sensing and other large-scale geospatial datasets has enabled the use of advanced data-driven approaches to classify catchments using traits such as topography, geology, climate, land cover, land use, and human influence. Unsupervised clustering algorithms based on the Euclidean distance are commonly used for trait-based classification but are not suitable for highly dimensional data. In this study we present a new network-based method for multi-scale catchment classification, which can be applied to large datasets and used to determine the traits associated with different catchment groups. In this framework, two networks are analyzed in parallel: the first being where the nodes are traits and the second being where the nodes are catchments. In both cases, edges represent pairwise similarity, and a network cluster detection algorithm is used for the classification. The trait network is used to investigate redundancy in the trait data and to condense this information into a small number of interpretable categories. The catchments network is used to classify the catchments into clusters and to identify representative catchments for the different groups using the degree centrality metric. We apply this method to classify 9067 river catchments across the contiguous United States at both regional and continental scales using 274 non-categorical traits. At the continental scale, we identify 25 interpretable trait categories and 34 catchment clusters of sizes greater than 50. We find that catchments with similar trait categories are typically located in the same region, with different spatial patterns emerging among clusters dominated by natural and anthropogenic traits. We also find that the catchment clusters exhibit distinct hydrological behavior based on an analysis of streamflow indices. This network approach provides several advantages over traditional means of classification, including better separation of clusters, the use of alternate similarity metrics that are more suitable for highly dimensional data, and reducing redundancy in the trait information. The paired catchment–trait networks enable analysis of hydrological behavior using the dominant trait categories for each catchment cluster. The approach can be used at multiple spatial scales since the network topologies adjust automatically to reflect the trait patterns at the scale of investigation. Finally, the representative catchments identified as hub nodes in the network can be used to guide transferable observational and modeling strategies. The method is broadly applicable beyond hydrology for classification of other complex systems that utilize different types of trait datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations

ABSTRACT The spatial distribution of mono-abundance populations (MAPs, selected in [Fe/H] and [Mg/Fe]) reflect the chemical and structural evolution in a galaxy and impose strong constraints on galaxy formation models. In this paper, we use APOGEE data to derive the intrinsic density distribution of MAPs in the Milky Way, after carefully considering the survey selection function. We find that a single exponential profile is not a sufficient description of the Milky Way’s disc. Both the individual MAPs and the integrated disc exhibit a broken radial density distribution; densities are relatively constant with radius in the inner Galaxy and rapidly decrease beyond the break radius. We fit the intrinsic density distribution as a function of radius and vertical height with a 2D density model that considers both a broken radial profile and radial variation of scale height (i.e. flaring). There is a large variety of structural parameters between different MAPs, indicative of strong structure evolution of the Milky Way. One surprising result is that high-α MAPs show the strongest flaring. The young, solar-abundance MAPs present the shortest scale height and least flaring, suggesting recent and ongoing star formation confined to the disc plane. Finally we derive the intrinsic density distribution and corresponding structural parameters of the chemically defined thin and thick discs. The chemical thick and thin discs have local surface mass densities of 5.62 ± 0.08 and 15.69 ± 0.32 M⊙pc−2, respectively, suggesting a massive thick disc with a local surface mass density ratio between thick to thin disc of 36 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.

58 GEOSCIENCES↗