Search NASA⌕ Search

SEARCH · Search NASA

Results for “Natural Fractures”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling Validation of Mechanical and Thermal Borehole Breakouts From Polyaxial Laboratory Tests to Determine Maximum Horizontal in Situ Stress

To accurately assess rock behavior in the deep subsurface, it is necessary to measure the in situ stress directions and magnitudes. The current methods for measuring the in situ stress state in the deep subsurface primarily include hydraulic fracturing tests (i.e., the minimum horizontal stress) or occasionally observing existing compressive borehole breakouts (i.e., the maximum horizontal stress) that have occurred naturally from drilling. If there are no existing compressive breakouts in a borehole, the maximum horizontal in situ stress cannot be estimated with much confidence. In response to this data gap, a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. In support of developing the thermal breakout technology, polyaxial laboratory tests have been performed on small-scale boreholes within rock blocks where mechanically- and thermally-induced borehole breakouts have been created. Numerical models along with the principal of superposition were created and used to analyze the polyaxial laboratory tests to predict the maximum horizontal stress, given the same data that would be obtained in an actual subsurface borehole thermal breakout test. Multiple failure criteria were used to evaluate the best prediction of the breakout onset and maximum horizontal stress. The maximum horizontal stress predictions were compared to the actual maximum horizontal stress applied at breakout using acoustic events recorded from emission sensors in the polyaxial tests. The results showed consistent results that could be used to refine the modeling approach and failure criterion that are used to make the maximum horizontal stress predictions. This study provided insight and validation for the thermal breakout stress measurement concept.

Jones, Matt↗
DOE OSTI2023

GeoThermalCloud for EGS – An Open-source, User-friendly, Scalable AI Workflow for Modeling Enhanced Geothermal Systems

Enhanced Geothermal Systems (EGS) offer a vast potential to expand the use of geothermal energy. Heat is extracted from this engineered system by injecting relatively cold water into subsurface fractures, which are in contact with hot dry rock, and brought back to surface through production wells. Creating EGS requires improving the natural permeability of hot crystalline rocks. In this short conference paper, we present a reproducible workflow for modeling EGS. Our workflow called the GeoThermalCloud (GTC) for EGS, leverages recent advances in machine learning, deep learning, and high-performance computing. This GTC framework is currently being made open-source, user-friendly, and reproducible through python scripts as well as Google Colab/Jupyter Notebooks. This GTC for EGS modeling scripts are made available at https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS and will constantly be updated to cater for geothermal community. Current GTC framework provides scripts to train deep learning (DL) models for techno-economics and data worth analysis. The Geothermal Design Tool (https://github.com/GeoDesignTool/GeoDT.git), a fast and simplified multi-physics solver, is used to develop a database for training DL models. This short paper provides details on the scripts to curate, process, and train DL models. The scripts can easily be modified to train on databases generated by other popular open-source simulators such as PFLOTRAN, STOMP, TOUGH, and GEOSX or commercial software such as ResFrac and COMSOL.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Effects of hydrogen partial pressure on crack initiation and growth rate in vintage X52 steel

There is an increasing interest in ensuring compatibility of existing natural gas infrastructure with conveyance of hydrogen, especially for vintage (pre-1970s) pipeline steels. While hydrogen is known to affect the fatigue and fracture resistance of these materials, there is a lack of data on its effects on crack initiation at low hydrogen partial pressure. To address this gap, the current work presents a series of fatigue tests on vintage API Grade X52 steel in pure nitrogen, air and gaseous hydrogen at partial pressures of 1 and 207 bar. Circumferential notch tension (CNT) specimens with direct current potential difference (DCPD) measurements are used to get S–N curves for crack initiation and failure, and compact tension (CT) tests to assess crack growth rates. Results indicate that air accelerates crack initiation, low hydrogen partial pressure affects the crack growth rate, and higher pressures have the largest impact on both. Surprisingly, crack initiation is faster in air than at 1 bar hydrogen partial pressure. The CNT tests also reveal that the orientation of the banded microstructure has no effect on the crack initiation site, but it influences its growth behavior. Moreover, fractography shows that hydrogen also promotes instances of quasi-cleavage and intergranular fracture, linked to embrittlement of the material. Here, these results demonstrate that, while hydrogen may impact the entirety of the fatigue process, its effects on the early stages of damage accumulation become relevant only at larger pressures. More comprehensive studies of crack initiation in hydrogen atmospheres will be critical to future codes and standards for hydrogen infrastructure.

Circumferential notch tension↗
DOE OSTI2020

The emergence of small-scale self-affine surface roughness from deformation

Most natural and man-made surfaces appear to be rough on many length scales. There is presently no unifying theory of the origin of roughness or the self-affine nature of surface topography. One likely contributor to the formation of roughness is deformation, which underlies many processes that shape surfaces such as machining, fracture, and wear. Using molecular dynamics, we simulate the biaxial compression of single-crystal Au, the high-entropy alloy Ni 36.67 Co 30 Fe 16.67 Ti 16.67 , and amorphous Cu 50 Zr 50 and show that even surfaces of homogeneous materials develop a self-affine structure. By characterizing subsurface deformation, we connect the self-affinity of the surface to the spatial correlation of deformation events occurring within the bulk and present scaling relations for the evolution of roughness with strain. These results open routes toward interpreting and engineering roughness profiles.

42 ENGINEERING↗
DOE OSTI2022

Salinity and hydraulic retention time induce membrane phospholipid acyl chain remodeling in Halanaerobium congolense WG10 and mixed cultures from hydraulically fractured shale wells

Bacteria remodel their plasma membrane lipidome to maintain key biophysical attributes in response to ecological disturbances. For Halanaerobium and other anaerobic halotolerant taxa that persist in hydraulically fractured deep subsurface shale reservoirs, salinity, and hydraulic retention time (HRT) are important perturbants of cell membrane structure, yet their effects remain poorly understood. Membrane-linked activities underlie in situ microbial growth kinetics and physiologies which drive biogeochemical reactions in engineered subsurface systems. Hence, we used gas chromatography–mass spectrometry (GC–MS) to investigate the effects of salinity and HRT on the phospholipid fatty acid composition of H. congolense WG10 and mixed enrichment cultures from hydraulically fractured shale wells. We also coupled acyl chain remodeling to membrane mechanics by measuring bilayer elasticity using atomic force microscopy (AFM). For these experiments, cultures were grown in a chemostat vessel operated in continuous flow mode under strict anoxia and constant stirring. Our findings show that salinity and HRT induce significant changes in membrane fatty acid chemistry of H. congolense WG10 in distinct and complementary ways. Notably, under nonoptimal salt concentrations (7% and 20% NaCl), H. congolense WG10 elevates the portion of polyunsaturated fatty acids (PUFAs) in its membrane, and this results in an apparent increase in fluidity (homeoviscous adaptation principle) and thickness. Double bond index (DBI) and mean chain length (MCL) were used as proxies for membrane fluidity and thickness, respectively. These results provide new insight into our understanding of how environmental and engineered factors might disrupt the physical and biogeochemical equilibria of fractured shale by inducing physiologically relevant changes in the membrane fatty acid chemistry of persistent microbial taxa. GRAPHICAL ABSTRACT Salinity significantly alters membrane bilayer fluidity and thickness in Halanaerobium congolense WG10.

03 NATURAL GAS↗
DOE OSTI2020

Sapphire advanced mitigation process: wet etch to expose sub-surface damage and increase laser damage resistance and mechanical strength

In this study, a novel, to the best of our knowledge, method of wet chemical etching of sapphire workpieces (such as optics, wafers, windows, and cones), called the sapphire advanced mitigation process (or sapphire AMP), has been developed that exposes sub-surface mechanical damage created during the optical fabrication process and significantly enhances the surface laser damage resistance (> 2×) and mechanical strength (up to ~ 2.6×). Sapphire AMP involves first treating the workpiece with a mixture of sulfuric and phosphoric acid ([H 2 SO 4 ]:[H 3 PO 4 ]= 1 : 3) at 220°C, followed with phosphoric acid at 160°C, then with sodium hydroxide base (NaOH) and surfactant at 40°C, and finally with a high-pressure deionized water spray rinse. Sapphire AMP has been demonstrated on both A- and C-plane sapphire workpieces. The mechanism of this etch process involves the reaction of the sapphire (Al 2 O 3 ) surface with sulfuric acid (H 2 SO 4 ) forming aluminum sulfate [Al2(SO4)3], which has low solubility. The high phosphoric acid content in the first and second steps of sapphire AMP results in the efficient conversion of Al 2 (SO 4 ) 3 to aluminum phosphate (AlPO 4 ), which is very soluble, greatly reducing reaction product redeposition on the workpiece surface. Sapphire AMP is shown to expose sub-surface mechanical damage on the sapphire surface created during the grinding and polishing processes, whose etched morphology has either isotropic or anisotropic evolution depending on the nature of the initial surface damage. Sapphire AMP was also designed to remove the key known surface, laser absorbing precursors (namely, foreign chemical impurities, the fracture surface layer of preexisting sub-surface damage, and reaction product or foreign species redeposition or precipitation). Static and sliding indention induced surface microfractures on sapphire are shown after sapphire AMP to have a significant decrease in the fast photoluminescence intensity (a known metric for measuring the degree of laser damaging absorbing precursors). In addition, the onset of laser damage (at 351 nm 3 ns) on sapphire AMP treated workpieces was shown to increase in fluence from ~ 4 to >9.5J/cm 2 . Finally, biaxial ball-on-ring mechanical tests on sapphire disks showed an increase in the failure stress from 340 MPa (with pre-existing 28 µm flaws) to ~ 900MPa after sapphire AMP, which is attributed to the blunting of the surface microfractures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The mechanics of brittle granular materials with coevolving grain size and shape

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Toughness of Interlocking Metasurfaces

Interlocking metasurfaces (ILMs) are arrays of autogenous latching unit cells patterned across a surface. These create structural joints similar to bioinspired suture joints but patterned over a 2D surface rather than a 1D seam. This enables ILMs to be an alternative to conventional joining technologies such as bolts, welds, and adhesives. However, compared to conventional joining methods, relatively little is known of the engineering considerations for designing structural ILMs. Herein, the interfacial toughness of an archetypal ILM is examined for the first time. Under the conditions studied here, the ILM is substantially tougher than the material from which it is made, in this case, exhibiting up to a 50% increase in interfacial crack initiation energy over the solid base material, a photocured 3D‐printed polymer. Through experimental tests using in‐situ digital image correlation along with complementary computational analyses, the mechanism of toughening in the ILM structure and the origins of toughness anisotropy are revealed. The increase in toughness is associated with cross‐cell interactions, that is, load‐sharing across unit cells, which give rise to a finite process zone length with different effective material properties. In this way, ILM toughening is analogous to crack blunting in ductile materials or fiber bridging in composites; yet here, the ILM is composed of a single‐phase base material and so the architected toughening is geometric in nature and hence amenable to future topological optimization.

fracture↗
DOE OSTI2024

Material selection and manufacturing for high‐temperature heat exchangers: Review of state‐of‐the‐art development, opportunities, and challenges

Abstract Many energy systems demand heat transfer at high temperatures to keep up with high demand for power, so high‐temperature material that can perform and last under these harsh conditions is needed for heat exchangers. The engineering requirements for these high‐temperature heat exchanger material call for high thermal conductivity, high resistance to fracture, high resistance to creep deformation, environmental stability in environments associated with the application, and high modulus of elasticity while maintaining low cost to make and maintain. Naturally, ceramics are a good solution for this endeavor. In the past, high‐temperature heat exchangers made from ceramics have been used. We provide examples of ceramics in relevant heat exchange applications and provide motivation where additive manufacturing (AM) can improve efficiency. AM for the relevant material is under development, and we provide insight on the AM of ceramic materials and examples of AM heat exchangers keeping cost in mind. The motivation of the review paper is to provide a framework for material and manufacturing selection for high‐temperature heat exchangers for AM to keep up with the demand for better efficiency, better material, better manufacturing, and cost moving forward with AM technology in high‐temperature ceramic heat exchangers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stress Corrosion Cracking of NaSICON Membranes in Aqueous Electrolytes for Redox-Flow Batteries

The sodium super-ionic conductor (NaSICON) has versatile applications as a ceramic electrolyte for energy storage, where it can serve as an impermeable separator in solid-state batteries and redox flow systems. In particular, NaSICON systems have been proposed to be relatively stable in contact with water, making them compatible with aqueous battery chemistries. However, owing to their brittle nature and metal oxide constituents, stress-corrosion cracking (SCC) is an important failure mechanism that has not been previously explored. In this study, we assess the fracture toughness of NaSICON membranes in contact with aqueous solutions that are relevant to redox flow systems. Microindentation was performed to generate visible surface cracks and residual stress, which were observed to grow in length after exposure to aqueous solutions. This allows for a quantitative measurement of fracture toughness, which decreases after exposure to water. To contextualize these results, we develop a simplified model of the fracture behavior in aqueous redox-flow batteries that incorporate NaSICON membranes, illustrating the importance of SCC in cell design. This work provides quantitative insights into SCC as a failure mode in NaSICON, enhancing our understanding of the chemo-mechanical behavior of ceramic electrolytes in contact with aqueous solutions.

Zhang, Mengyao↗
DOE OSTI2020

Effects of HFIR neutron irradiation on fracture toughness properties of standard and Ni-doped F82H

F82H is the Japanese reference reduced-activation ferritic-martensitic (RAFM) steel for fusion blanket applications. The harsh environment of a fusion reactor, such as neutron irradiation and He/H damage, can result in significant degradation of F82H fracture toughness. Therefore, understanding the fracture toughness behavior of F82H in the fusion environment is critical to ensure the long-term safe operation of the fusion reactor. Here, we summarize seven irradiation campaigns of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL) covering five variants of F82H steels, including F82H IEA, F82H Mod3, F82H doped with 1.4% natural Ni, F82H doped with 1.4% 58 Ni, and F82H doped with 1.4% 60 Ni. The irradiation temperatures covered the range from 220 °C to 530 °C and the neutron irradiation dose spanned 4 dpa to 70 dpa. The effects of neutron irradiation temperature, dose, materials composition, Ni doping, and He production on F82H fracture toughness are discussed. Our results showed that irradiation embrittlement monotonically decreased with increasing irradiation temperature until 400 °C for F82H IEA and F82H Mod3. F82H Mod3 showed better fracture toughness than F82H IEA both before and after neutron irradiation. We determined that 1.4% Ni alloying can be applied to F82H for simulating He effect in a fission reactor without jeopardizing the fracture toughness of the material. However, more studies are needed to understand the effect of high dose (>20 dpa) and He production on F82H fracture toughness.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insights into the soft brittle-to-ductile transition from discrete dislocation dynamics

The Brittle-to-ductile transition (BDT) in body centered cubic metals exhibits a soft transition wherein the fracture toughness gradually rises to before the onset of ductility. The resultant brittle-to-ductile transition temperature can be described with an Arrhenius relationship whose activation energy is related to plasticity in the material. To provide further insight into the nature of the BDTT, in this work we utilized a discrete dislo- cation dynamics model with a crack to simulate the BDT and how it depends on the thermally activated nature of plasticity. The interrelationship between the BDT activation energy and the dislocation mobility parameters were determined via the calculation of first order sensitivity coefficients. This analysis allows us to demonstrate that the activation energy for the BDT is directly related to the activation energy for plasticity through an effective stress that defines this relationship. This effective stress physically is the average stress on the dislocations that move out of the crack. Lastly, we are able to show that this effective stress is dictated by the low temperature fracture toughness or cleave energy of the material and the source position, the latter of which can be affected by processing. Collectively, these results provide new insight into what controls the thermal activation of the BDT and what are the important parameters to control it.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗