Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Image Characterization of Reactions Generated by an Aeroengine Micromixing Injector for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blend

An aeroengine micromixing injector, originally designed for lean direct injection of jet fuel, was adapted to work with hydrogen, natural gas, and any blend in between. The ultimate goal of the design was to achieve low NOx emissions when operated on pure hydrogen; to better characterize its performance, flame imaging diagnostics was implemented. Three types of cameras mounted parallel to the injector base and pointing towards the quartz combustor cylinder with the same angle were operated simultaneously: Nikon D90, Dynacolor FB-N9-U, and Phantom v7.1, to capture the visible spectrum, OH* chemiluminescence, and dynamic behavior of the flames, respectively. The current work also includes the overlapping of visible spectrum and OH* chemiluminescence, as an effort to qualitatively define the heat release over the flame area, i.e., UV over visible domain. Additionally, a z-type schlieren configuration was employed to reinforce the presence of some small-scale details occurring near the injector outlet ports. Using the airbox pressure drop, preheat temperature, fuel composition and flame temperature as the factors of study, a Box Behnken model was designed yielding into a 16-points matrix for the operability region, at atmospheric pressure. The main responses were extracted from the OH* images, because of its correlation with the flame heat release, these being the flame area and its average brightness, the heat release area and its center of gravity and leading edge. To assess the importance of the main factors and their interactions, an analysis of variance (ANOVA) was performed for each response, considering as significant each contributor with a p-value below 0.05. For these five responses it was found that all four main factors need to be included into the model (most of them because its own significance, others to ensure the model hierarchy), as well as some common interactions as the preheat and flame temperatures or pressure drop and fuel composition, and quadratic terms from fuel composition or flame temperature. Three of the five models presented R2 and coefficient of variance (C.V.) indicators around 0.97 and 5%, respectively; the leading edge ANOVA yielded into R2=0.78 and C.V.=18.88%; and the heat release area showed R2=0.91 and C.V.=24.13%. Lastly, the confinement ratio effect of the combustor over the flame structure was quantified, for three different conditions: the original 80 mm inner diameter and 200 mm length cylinder was compared against a shorter 150 mm tube (same inner diameter) and a narrower 47 mm inner diameter one (same 200 mm length). The length had a smaller effect than the inner diameter, even though both factors were not strong enough to make a significant deviation on most of the flame parameters: the difference was within the repeatability margin of error for the brightness, flame area, and center of gravity and leading edge of the heat release area.

Imaging diagnostics, flame characterization, OH* c↗
DOE OSTI2025

Nonlocal nucleon matrix elements in the rest frame

Extracting parton structure from lattice quantum chromodynamics (QCD) calculations requires studying the coordinate scale 𝑧 3 dependence of the matrix elements of bilocal operators. The most significant contribution comes from the 𝑧 3 dependence induced by ultraviolet (UV) renormalization of the Wilson line. We demonstrate that the next-to-leading order perturbative calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice QCD rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. The residual discrepancies can be modeled by a leading effect from the structure of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Direct quarkonium production in DIS from a joint CGC and NRQCD framework

We compute the differential cross section for direct quarkonium production in high-energy electron-nucleus collisions at small 𝑥. Our computation is performed within the nonrelativistic QCD factorization formalism that separates the calculation into short distance coefficients and long distance matrix elements that depend on the color and spin of the state. We obtain the short distance coefficients of the production of the heavy quark pair within the framework of the color glass condensate effective field theory, which resums coherent multiple interactions of the heavy quark pair with the nucleus to all orders. Our results are expressed as the convolution of perturbatively calculable functions with multipoint lightlike Wilson line correlators. In the correlation limit, we establish the correspondence between our color glass condensate formulation with calculations employing the transverse momentum dependent (TMD) framework. We extend this correspondence by resumming kinematic power corrections within the improved TMD framework, which interpolates between the TMD formalism and 𝑘 ⊥ -factorization formalism. We present a detailed numerical analysis, focusing on 𝐽/𝜓 production in the kinematics accessible at the future Electron-Ion Collider, highlighting the importance of genuine higher-order saturation contributions when the electron collides with a large nucleus. Our results are also valid in the photoproduction limit where we expect the largest contribution from genuine higher-order saturation contributions which could be accessed in ultraperipheral collisions of relativistic heavy ions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2021

Comparison on the properties of ITZs in fly ash-based geopolymer and Portland cement concretes with equivalent flowability

Highlights: • Modelled ITZs can effectively avoid the effects caused by different factors and test locations of ITZ. • ITZ properties of PC and geopolymer concretes are compared based on equivalent flowability of binders. • Geopolymer ITZ exhibits stronger bonding than that of PC, and higher mechanical properties than geopolymer paste. • Strong ITZs in geopolymer concrete are attributed to the layer of gel-rich paste with dense microstructures. • Elastic modulus and hardness properties of ITZs in PC concrete are increased due to the CH crystals. This paper aims to compare the properties of interfacial transition zones (ITZs) in Portland cement (PC) concrete and geopolymer concrete. Portland cement and geopolymer pastes were designed with the equivalent flowability to provide similar mix and casting condition of ITZs. Two types of modelled ITZs were prepared to facilitate the nanoindentations across ITZs, microstructural characterization, and comparison on the properties of ITZs with less influential factors. The results showed that the interfacial bonding of ITZs between geopolymer matrix and aggregate is relatively stronger than the counterpart in the PC concrete. There is a high amount of crystalline hydration products in the ITZs of PC concrete, but a layer gel-rich paste with denser microstructures in the ITZs of geopolymer concrete. Additionally, the interface morphology and nanoindentation analysis indicate that the property of ITZs in the modelled geopolymer concrete is not poorer than that of the corresponding geopolymer paste.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Challenges and physiological implications of wood feeding in termites

Termites are fascinating insects for a number of reasons, one of which being their specialization on diets of wood lignocellulose. The goal of this review is to consider stress-inducing characteristics of wood and apparent molecular-physiological adaptations in termite guts to overcome these stressors. Defensive factors present in wood include extractive secondary plant metabolites, lignin and related phenolics, crystalline cellulose, and low nitrogen content. Molecular-physiological adaptations of the termite gut to deal with these factors include robust detoxification and antioxidant machinery, the production of a peritrophic matrix and a wide range of cellulases from host and symbiotic sources, and creation of niches available to nitrogen-fixing bacterial symbionts. Considering termite gut physiology and symbioses in the context of stress-response has applied implications. Overall, these outcomes can include development of efficient biomass breakdown strategies, protection of microbes during industrial processing applications, and safeguarding wooden structures from termite damage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Probing quantum entanglement with generalized parton distributions at the Electron-Ion Collider

Within the collinear factorization framework based on generalized parton distributions, we calculate the spin density matrix of exclusively produced quark and antiquark pairs $𝑢\bar{𝑢}$, $𝑑⁢\bar{d}$, $𝑠⁢\bar{𝑠}$, $𝑐\bar{⁢𝑐}$, $𝑏\bar{𝑏}$ in electron-proton scattering. The presence of both real and imaginary parts in the scattering amplitudes leads to a rich pattern of entanglement between the quark and the antiquark. We map out kinematical regions where the pairs exhibit entanglement, Bell nonlocality, and nonstabilizerness. We also predict that massive quarks and antiquarks are transversely polarized, similar to the well-known transverse hyperon polarization in unpolarized collisions. In strangeness, charm, and bottom productions, the polarization can reach 50%–80% in certain kinematic regions in the low-energy runs of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Carbon Storage Technical Viability Approach (CS TVA): An Integrated Approach for Feasibility and Data Resource Assessment

There is currently a poor understanding and lack of workflow to understand the technical viability of carbon storage spatially. To address this gap, the multi-faceted Carbon Storage Technical Viability Approach (CS TVA) is being developed to incorporate CO2 storage resources, environmental and socio-economic justice (EJ/SJ) factors to enable more comprehensive assessments. The CS TVA includes a (1) matrix framework, (2) an integrated and labeled database, (3) a data availability assessment workflow, and (4) spatial data availability assessment results. This approach leverages spatial and data science analytics to communicate data density, uncertainty, and gaps. The workflow can be applied in whole or in part, based on user needs.

Rodriguez, Neyda Cordero↗
DOE OSTI2020

Ultra-wideband mutual coupling compensation of active electronically scanned arrays in multi-channel radar systems

Technologies pertaining to the calibration of ultra-wide-band active-electronically-scanned arrays to compensate for transmit and receive mutual coupling between array elements and/or non-ideal isolated channels are described herein. A plurality of near-field measurements are taken by way of probe or probes that may be moved among a plurality of positions aligned with respective elements in the array. For each position of the probe, each of the elements of the array is stimulated to transmit or receive a calibration signal to or from the probe, respectively. Frequency-domain transfer functions are computed from the received signals by the probe or the element for each of the array elements in each of the positions of the probe. The inverse of the matrix of transfer functions comprise frequency-domain transmit and receive correction factors that are used to modify desired array inputs/outputs such that the modified signals correct for mutual coupling between elements in the array.

42 ENGINEERING↗
DOE OSTI2025

Environmental matrix and moisture influence soil microbial phenotypes in a simplified porous media incubation

Soil moisture and porosity regulate microbial metabolism by influencing factors, such as system chemistry, substrate availability, and soil connectivity. However, accurately representing the soil environment and establishing a tractable microbial community that limits confounding variables is difficult. Here, we use a reduced-complexity microbial consortium grown in a glass bead porous media amended with chitin to test the effects of moisture and a structural matrix on microbial phenotypes. Leveraging metagenomes, metatranscriptomes, metaproteomes, and metabolomes, we saw that our porous media system significantly altered microbial phenotypes compared with the liquid incubations, denoting the importance of incorporating pores and surfaces for understanding microbial phenotypes in soils. These phenotypic shifts were mainly driven by differences in expression of Streptomyces and Ensifer, which included a significant decrease in overall chitin degradation between porous media and liquid. Our findings suggest that the success of Ensifer in porous media is likely related to its ability to repurpose carbon via the glyoxylate shunt amidst a lack of chitin degradation byproducts while potentially using polyhydroxyalkanoate granules as a C source. We also identified traits expressed by Ensifer and others, including motility, stress resistance, and carbon conservation, that likely influence the metabolic profiles observed across treatments. Together, these results demonstrate that porous media incubations promote structure-induced microbial phenotypes and are likely a better proxy for soil conditions than liquid culture systems. Furthermore, they emphasize that microbial phenotypes encompass not only the multi-enzyme pathways involved in metabolism but also include the complex interactions with the environment and other community members.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A DE1 BINDING FACTOR 1–GLABRA2 module regulates rhamnogalacturonan I biosynthesis in Arabidopsis seed coat mucilage

Abstract The mucilage surrounding hydrated Arabidopsis thaliana seeds is a specialized extracellular matrix composed mainly of the pectic polysaccharide rhamnogalacturonan I (RG-I). Although, several genes responsible for RG-I biosynthesis have been identified, the transcriptional regulatory mechanisms controlling RG-I production remain largely unknown. Here we report that the trihelix transcription factor DE1 BINDING FACTOR 1 (DF1) is a key regulator of mucilage RG-I biosynthesis. RG-I biosynthesis is significantly reduced in loss-of-function mutants of DF1. DF1 physically interacts with GLABRA2 (GL2) and both proteins transcriptionally regulate the expression of the RG-I biosynthesis genes MUCILAGE MODIFIED 4 (MUM4) and GALACTURONOSYLTRANSFERASE-LIKE5 (GATL5). Through chromatin immunoprecipitation-quantitative PCR and transcriptional activation assays, we uncover a cooperative mechanism of the DF1–GL2 module in activating MUM4 and GATL5 expression, in which DF1 binds to the promoters of MUM4 and GATL5 through interacting with GL2 and facilitates the transcriptional activity of GL2. The expression of DF1 and GL2 is directly regulated by TRANSPARENT TESTA GLABRA2 (TTG2) and, in turn, DF1 directly represses the expression of TTG2. Taken together, our data reveal that the transcriptional regulation of mucilage RG-I biosynthesis involves a regulatory module, comprising DF1, GL2, and TTG2.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Nucleon axial and pseudoscalar form factors using twisted-mass fermion ensembles at the physical point

We compute the nucleon axial and pseudoscalar form factors using three N f = 2 + 1 + 1 twisted-mass fermion ensembles with all quark masses tuned to approximately their physical values. The values of the lattice spacings of these three physical point ensembles are 0.080, 0.068, and 0.057 fm and spatial sizes 5.1, 5.44, and 5.47 fm, respectively, yielding m π L > 3.6 . Convergence to the ground-state matrix elements is assessed using multistate fits. We study the momentum dependence of the three form factors and check the partially conserved axial-vector current (PCAC) hypothesis and the pion pole dominance (PPD). We show that in the continuum limit, the PCAC and PPD relations are satisfied. We also show that the Goldberger-Treimann relation is approximately fulfilled and determine the Goldberger-Treiman discrepancy. Our final results are g A = 1.245 ( 28 ) ( 14 ) for the nucleon axial charge, ⟨ r A 2 ⟩ = 0.339 ( 48 ) ( 06 ) fm 2 for the axial radius, g π N N ≡ lim Q 2 → − m π 2 G π N N ( Q 2 ) = 13.25 ( 67 ) ( 69 ) for the pion-nucleon coupling constant, and G P ( 0.88 m μ 2 ) ≡ g P * = 8.99 ( 39 ) ( 49 ) for the induced pseudoscalar form factor at the muon capture point. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2022

Decoherence factor as a convolution: an interplay between a Gaussian and an exponential coherence loss

Abstract We identify and investigate the origin and nature of the transition between Gaussian and exponential forms of decoherence: the decoherence factor (that controls the time dependence of the off-diagonal terms of the density matrix expressed in the pointer basis representation) is the convolution of the Fourier transforms of the spectral density and of the overlap (between the eigenstates the environment with and without couplings to the system). Spectral density alone tends to lead to the (approximately) Gaussian decay of coherence while the overlap alone results in a (largely) exponential decay. We show that these two contributions combine as a convolution, their relative importance controlled by the strength of the system-environment coupling. The resulting decoherence factor in the strong and weak coupling limits leads to predominantly Gaussian or exponential decay, respectively, as is demonstrated with two paradigmatic examples of decoherence—a spin-bath model and the quantum Brownian motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Renormalization and scale evolution of the soft-quark soft function

Soft functions defined in terms of matrix elements of soft fields dressed by Wilson lines are central components of factorization theorems for cross sections and decay rates in collider and heavy-quark physics. While in many cases the relevant soft functions are defined in terms of gluon operators, at subleading order in power counting soft functions containing quark fields appear. We present a detailed discussion of the properties of the soft-quark soft function consisting of a quark propagator dressed by two finite-length Wilson lines connecting at one point. This function enters in the factorization theorem for the Higgs-boson decay amplitude of the h → γγ process mediated by light-quark loops. We perform the renormalization of this soft function at one-loop order, present a conjecture for its two-loop anomalous dimension and discuss solutions to its renormalization-group evolution equation in momentum space, in Laplace space and in the “diagonal space”, where the evolution is strictly local in the momentum variable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Chlorophylls as primary electron acceptors in reaction centers

Photosynthetic Reaction Centers (RCs) can be considered blueprints for highly efficient energy transfer. Embedded with an array of cofactors, including (bacterio)chlorophyll ((B)Chl) and (B)pheophytin ((B)Pheo) molecules, RCs function with a high quantum yield that spans a wide spectral range. Understanding the principles that underlie their function can influence the design of the next generation of artificial photosynthetic devices. We are particularly interested in the factors that influence the early stages of light-driven charge separation in RCs. With the recent publication of several highly anticipated RC structures and advanced computational methods available, it is possible to probe both the geometric and electronic structures of an array of RCs. In this chapter, we review the electronic and geometric structures of the (B)Chl and (B)Pheo primary electron acceptors from fiveRCs, comprising both Type I and Type II RCs and representing both heterodimeric and homodimeric systems. We showcase the dimeric A0•- state of Type I RCs, whereby the unpaired electron is delocalized, to various extents, over two (B)Chl molecules (B)Chl2 and (B)Chl3. This delocalization is controlled by several factors, including the structure of the (B)Chls, interactions with the surrounding protein matrix, and the orientation and distances of the cofactors themselves. In contrast, the primary acceptors of Type II RCs are entirely monomeric, with electron density residing solely on the (B)Pheo. We compare the natural design of the primary acceptors of the Type I and Type II RCs from both an evolutionary and application based perspective.

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

Statistical Multiobjective Optimization of Thiospinel CoNi 2 S 4 Nanocrystal Synthesis via Design of Experiments

Thiospinels, such as CoNi 2 S 4 , are showing promise for numerous applications, including as catalysts for the hydrogen evolution reaction, hydrodesulfurization, and oxygen evolution and reduction reactions; however, CoNi 2 S 4 has not been synthesized as small, colloidal nanocrystals with high surface-area-to-volume ratios. Traditional optimization methods to control nanocrystal attributes such as size typically rely upon one variable at a time (OVAT) methods that are not only time and labor intensive but also lack the ability to identify higher-order interactions between experimental variables that affect target outcomes. Herein, we demonstrate that a statistical design of experiments (DoE) approach can optimize the synthesis of CoNi 2 S 4 nanocrystals, allowing for control over the responses of nanocrystal size, size distribution, and isolated yield. After implementing a 2 5–2 fractional factorial design, the statistical screening of five different experimental variables identified temperature, Co:Ni precursor ratio, Co:thiol ratio, and their higher-order interactions as the most critical factors in influencing the aforementioned responses. Second-order design with a Doehlert matrix yielded polynomial functions used to predict the reaction parameters needed to individually optimize all three responses. A multiobjective optimization, allowing for the simultaneous optimization of size, size distribution, and isolated yield, predicted the synthetic conditions needed to achieve a minimum nanocrystal size of 6.1 nm, a minimum polydispersity (σ/$\bar{d}$) of 10%, and a maximum isolated yield of 99%, with a desirability of 96%. The resulting model was experimentally verified by performing reactions under the specified conditions. Furthermore, our work illustrates the advantage of multivariate experimental design as a powerful tool for accelerating control and optimization in nanocrystal syntheses.

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

Fracture-matrix fluid exchange in oil-bearing unconventional mudstones

The poromechanical properties of unconventional reservoir materials are in large part dictated by their mineralogy. Since these properties govern the response to stress experienced during hydraulic fracturing, fluid production, and fluid injection, they play a central role in the formation of microcracks or bedding delaminations which ultimately dominate mass transport. In this work we study access to the porosity of end member unconventional reservoir materials, where the end members are predominantly dictated by carbonate content. Access to the porosity is quantified using state of the art 3D x-ray computed tomography coupled with physics informed data analytics. Xenon gas, which attenuates x-rays, provides a spatiotemporal map of access to the porosity. The accessible porosity is quantified over a range of net confining stress relevant to the manmade disturbances listed above. These experiments demonstrate that heavily carbonated mudstones are nearly impermeable at the core (~ cm) scale, while carbonate free analogues afford better access to the microstructure. Consistent with previous qualitative 2D radiographs, access to the interior of the clastic mudstones is first observed along planar microcracks, followed by slow penetration into the surrounding matrix. Physics informed data analytics of the 3D tomography measurements presented here show that these microcracks do not permit uniform access to the adjacent rock matrix. In addition, variation of the effective pressure elucidates the mechanisms that govern fracture/matrix fluid exchange. Under conditions consistent with hydrocarbon production fluid accumulates in the immediate vicinity of the nearest microcrack. While there is clear evidence that, as intended, part of this accumulation is from the more distant matrix, fluid is also squeezed out of the microcrack. The fluid build-up at the microcrack indicates that migration out of the rock is hindered by the coupled poroelastic response of the microcrack and adjacent rock matrix. We show that these mechanisms ultimately account for the meager oil recovery factors realized in practice. These insights have implications for making reservoir scale predictions based on core scale observations, and provide a basis for devising new asset development techniques to access more porosity, and enhance fluid extraction. Finally, these findings shed light on key features and mechanisms that govern shale storage capacity, with relevance to other important industrial processes, such as geologic CO 2 storage.

58 GEOSCIENCES↗