Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Network”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Computational analysis of a pipe flow distributor for a thermocline based thermal energy storage system

In this paper, the overall efficiency of a concentrating solar power (CSP) plant depends on the effectiveness of thermal energy storage (TES) system. A single tank TES system consists of a thermocline region which produces the temperature gradient between hot and cold storage fluid by density difference. Preservation of this thermocline region in the tank during charging and discharging cycles depends on the uniformity of the velocity profile at any horizontal plane. Our objective is to maximize the uniformity of the velocity distribution using a pipe-network distributor by varying the number of holes, distance between the holes, position of the holes and number of distributor pipes. For simplicity, we consider that the diameter of the inlet, main pipe, the distributor pipes and the height and the width of the tank are constant. We use Hitec ® molten salt as the storage medium and the commercial software Gambit 2.4.6 and Fluent 6.3 for the computational analysis. We analyze the standard deviation in the velocity field and compare the deviations at different positions of the tank height for different configurations. Since the distance of the holes from the inlet and their respective arrangements affects the flow distribution throughout the tank; we investigate the impacts of rearranging the holes position on flow distribution. Impact of the number of holes and distributor pipes are also analyzed. We analyze our findings to determine a configuration for the best case scenario.

14 SOLAR ENERGY↗
DOE OSTI2014

Reverse Monte Carlo modeling in confined systems

An extension of the well established Reverse Monte Carlo (RMC) method for modeling systems under close confinement has been developed. The method overcomes limitations induced by close confinement in systems such as fluids adsorbed in microporous materials. As a test of the method, we investigate a model system of {sup 36}Ar adsorbed into two zeolites with significantly different pore sizes: Silicalite-I (a pure silica form of ZSM-5 zeolite, characterized by relatively narrow channels forming a 3D network) at partial and full loadings and siliceous Faujasite (which exhibits relatively wide channels and large cavities). The model systems are simulated using grand canonical Monte Carlo and, in each case, its structure factor is used as input for the proposed method, which shows a rapid convergence and yields an adsorbate microscopic structure in good agreement with that of the model system, even to the level of three body correlations, when these are induced by the confining media. The application to experimental systems is straightforward incorporating factors such as the experimental resolution and appropriate q-sampling, along the lines of previous experiences of RMC modeling of powder diffraction data including Bragg and diffuse scattering.

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

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
NASA NTRS2015

Overview of Global/Regional Models Used to Evaluate Tropospheric Ozone in North America

Ozone (O3) is an important greenhouse gas, toxic pollutant, and plays a major role in atmospheric chemistry. Tropospheric O3 which resides in the planetary boundary layer (PBL) is highly reactive and has a lifetime on the order of days, however, O3 in the free troposphere and stratosphere has a lifetime on the order of weeks or months. Modeling O3 mixing ratios at and above the surface is difficult due to the multiple formation/destruction processes and transport pathways that cause large spatio-temporal variability in O3 mixing ratios. This talk will summarize in detail the global/regional models that are commonly used to simulate/predict O3 mixing ratios in the United States. The major models which will be focused on are the: 1) Community Multi-scale Air Quality Model (CMAQ), 2) Comprehensive Air Quality Model with Extensions (CAMx), 3) Goddard Earth Observing System with Chemistry (GEOS-Chem), 4) Real Time Air Quality Modeling System (RAQMS), 5) Weather Research and Forecasting/Chemistry (WRF-Chem) model, National Center for Atmospheric Research (NCAR)'s Model for OZone And Related chemical Tracers (MOZART), and 7) Geophysical Fluid Dynamics Laboratory (GFDL) AM3 model. I will discuss the major modeling components which impact O3 mixing ratio calculations in each model and the similarities/differences between these models. This presentation is vital to the 2nd Annual Tropospheric Ozone Lidar Network (TOLNet) Conference as it will provide an overview of tools, which can be used in conjunction with TOLNet data, to evaluate the complex chemistry and transport pathways controlling tropospheric O3 mixing ratios.

Mixing Ratios↗
NASA NTRS2020

NeMO-Net – The Neural Multi-Modal Observation & Training Network for Global Coral Reef Assessment

We present NeMO-Net, the Srst open-source deep convolutional neural network (CNN) and interactive learning and training software aimed at assessing the present and past dynamics of coral reef ecosystems through habitat mapping into 10 biological and physical classes. Shallow marine systems, particularly coral reefs, are under significant pressures due to climate change, ocean acidification, and other anthropogenic pressures, leading to rapid, often devastating changes, in these fragile and diverse ecosystems. Historically, remote sensing of shallow marine habitats has been limited to meter-scale imagery due to the optical effects of ocean wave distortion, refraction, and optical attenuation. NeMO-Net combines 3D cm-scale distortion-free imagery captured using NASA FluidCam and Fluid lensing remote sensing technology with low resolution airborne and spaceborne datasets of varying spatial resolutions, spectral spaces, calibrations, and temporal cadence in a supercomputer-based machine learning framework. NeMO-Net augments and improves the benthic habitat classification accuracy of low-resolution datasets across large geographic ad temporal scales using high-resolution training data from FluidCam.NeMO-Net uses fully convolutional networks based upon ResNet and ReSneNet to perform semantic segmentation of remote sensing imagery of shallow marine systems captured by drones, aircraft, and satellites, including WorldView and Sentinel. Deep Laplacian Pyramid Super-Resolution Networks (LapSRN) alongside Domain Adversarial Neural Networks (DANNs) are used to reconstruct high resolution information from low resolution imagery, and to recognize domain-invariant features across datasets from multiple platforms to achieve high classification accuracies, overcoming inter-sensor spatial, spectral and temporal variations.Finally, we share our online active learning and citizen science platform, which allows users to provide interactive training data for NeMO-Net in 2D and 3D, integrated within a deep learning framework. We present results from the PaciSc Islands including Fiji, Guam and Peros Banhos 1 1 2 1 3 1 where 24-class classification accuracy exceeds 91%.

Chirayath, Ved↗
DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗
DOE OSTI2011

Alkali-free bioactive glasses for bone tissue engineering: A preliminary investigation

An alkali-free series of bioactive glasses has been designed and developed in the glass system CaO-MgO-SiO2-P2O5-CaF2 along diopside (CaMgSi2O6) – fluorapatite [Ca5(PO4)3F] – tricalcium phosphate (3CaO•P2O5) join. The silicate network in all the investigated glasses is predominantly coordinated in Q2 (Si) units while phosphorus tends to remain in orthophosphate (Q0) environment. The in vitro bioactivity analysis of glasses has been made by immersion of glass powders in simulated body fluid (SBF) while chemical degradation has been studied in Tris-HCl in accordance with ISO-10993-14. Some of the investigated glasses exhibit hydroxyapatite (HA) formation on their surface with in 1-12 h of their immersion in SBF solution. The sintering and crystallization kinetics of glasses has been investigated by differential thermal analysis (DTA) and hot-stage microscopy (HSM), respectively while the crystalline phase evolution in resultant glass-ceramics (GCs) has been studied in the temperature range of 800-900 oC using powder X-ray diffraction (XRD) and scanning electron microscope (SEM). The cell growth and osteogenic differentiation for glasses has been studied in vitro on sintered glass powder compacts using rat bone marrow mesenchymal stem cells. The as designed glasses are ideal candidates for their potential applications in bone tissue engineering in the form of bioactive glasses as well as glass/GC scaffolds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗
DOE OSTI2016

Tuning the tetrahedrality of the hydrogen-bonded network of water: Comparison of the effects of pressure and added salts

Experiments and simulations demonstrate some intriguing equivalences in the effect of pressure and electrolytes on the hydrogen-bonded network of water. Here, we examine the extent and nature of equivalence effects between pressure and salt concentration using relationships between structure, entropy, and transport properties based on two key ideas: first, the approximation of the excess entropy of the fluid by the contribution due to the atom-atom pair correlation functions and second, Rosenfeld-type excess entropy scaling relations for transport properties. We perform molecular dynamics simulations of LiCl–H{sub 2}O and bulk SPC/E water spanning the concentration range 0.025–0.300 molefraction of LiCl at 1 atm and pressure range from 0 to 7 GPa, respectively. The temperature range considered was from 225 to 350 K for both the systems. To establish that the time-temperature-transformation behaviour of electrolyte solutions and water is equivalent, we use the additional observation based on our simulations that the pair entropy behaves as a near-linear function of pressure in bulk water and of composition in LiCl–H{sub 2}O. This allows for the alignment of pair entropy isotherms and allows for a simple mapping of pressure onto composition. Rosenfeld-scaling implies that pair entropy is semiquantitatively related to the transport properties. At a given temperature, equivalent state points in bulk H{sub 2}O and LiCl–H{sub 2}O (at 1 atm) are defined as those for which the pair entropy, diffusivity, and viscosity are nearly identical. The microscopic basis for this equivalence lies in the ability of both pressure and ions to convert the liquid phase into a pair-dominated fluid, as demonstrated by the O–O–O angular distribution within the first coordination shell of a water molecule. There are, however, sharp differences in local order and mechanisms for the breakdown of tetrahedral order by pressure and electrolytes. Increasing pressure increases orientational disorder within the first neighbour shell while addition of ions shifts local orientational order from tetrahedral to close-packed as water molecules get incorporated in ionic hydration shells. The variations in local order within the first hydration shell may underlie ion-specific effects, such as the Hofmeister series.

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

Utah FORGE: Documentation on Discrete Fracture Network and Fracture Propagation Modelling

This dataset includes reports and a slide presentation on discrete fracture network (DFN) generation and hydraulic fracture modeling at the Utah FORGE site. It details the characterization of natural fractures using well log and core data, as well as stochastic modeling techniques. The reports describe simulations of hydraulic fracture propagation, fluid-mechanical interactions, and induced microseismicity. The dataset also includes history-matching of net pressure and analyses of fracture growth in naturally fractured geothermal reservoirs. The slides summarize key findings and future research directions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

FY20 Proxy App Suite Release: Report for ECP Proxy App Project Milestone ADCD-504-10

Version 4.0 of the ECP Proxy App Suite is practically unchanged from the previous release. The current set of proxies has proven useful for many aspects of benchmarking and co-design and we see little reason to alter the suite. Although there have been few changes to the ECP suite, the team has been hard at work in other areas. In the area of Machine Learning (ML) we have now created a separate proxy suite dedicated to this scientific applications of ML. The suite includes: miniGAN (Generative Adversarial Networks), miniRL (Reinforcement Learning), CRADL (inline inference), Cosmoflow-Benchmark (Convolutional Neural Network), and MLPerf-DeepCam (Climate Segmentation Benchmark). Section 3 contains more information about these proxies as well as the principles that are guiding the development of the suite. We have surveyed available proxies for several application domains including Computational Fluid Dynamics, Quantum Chemistry, Quantum Computing Simulation, Molecular Dynamics, Monte Carlo Transport, and Density Functional Theory to identify gaps in proxy coverage. Several new proxy apps are either already available or will be released soon to fill these gaps. Section 4 provides full details. Finally, section 5 reports on our continued collaboration with the ECP Continuous Integration (CI) effort to use proxy apps to help identify problems and roadblocks to cross-lab CI. We have also assisted the El Capitan Center of Excellence (COE) to stand up a CI system that will be used to test releases of HPE and AMD software stacks. We hope that the COE effort can serve as a model for ECP by showing how the Proxy App Team can work with AD and ST teams identify critical features, kernels, patterns, etc. and incorporate them into a CI system that will help ensure that Frontier and Aurora will provide those needed capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Deep Learning At Depth: Estimating subsurface parameters from geophysical monitoring data

Geophysical imaging techniques are a non-invasive way to image the subsurface and understand both subsurface solid (rock/soil) and fluid property distributions and their evolution in time. Inversions of the geophysical data, such as Electrical Resistance Tomography (ERT) data, are solved to estimate the subsurface property distributions, such as conductivity, and many inversion techniques smooth out sharp gradients in rock or fluid property distributions. Sharp gradients in subsurface properties tend to be present in situations with complex subsurface structures, which are common in many subsurface applications. We have successfully demonstrated that it is possible to inform, or constrain, inversions with neural networks trained on synthetic data with complex subsurface structures. Initial results suggest this process may be optimizable to yield property distributions that better represent the true property distributions than the same inversion process without the neural network constraint. Future work would optimize the neural network performance for this application and then apply the synthetic-data trained neural network to real data to understand the utility and performance of this technique for real data sets.

47 OTHER INSTRUMENTATION↗
NASA NTRS1976

A higher order panel method for general analysis and design applications in subsonic flow

A higher-order panel method is described for numerical solution of boundary-value problems relating to steady inviscid irrotational incompressible subsonic fluid flow in a domain. Both Neumann and Dirichlet boundary conditions are treated; two types of auxiliary conditions are used to remove the degrees of freedom that arise from specifying only the derivative of the perturbation velocity potential. Four general network types and two expansions of the induced potential kernel are employed in the numerical solution. Some results are presented which illustrate the modeling options and numerical characteristics of the method.

Johnson, F. T.↗
DOE OSTI2024

Changes in environmental and engineered conditions alter the plasma membrane lipidome of fractured shale bacteria

ABSTRACT Microorganisms that persist in fractured shale reservoirs cause several problems including secreting foul gases and forming biofilms. Current biocontrol measures often fail due to limited knowledge of their in situ activities. The plasma membrane protects the cell, mediates many of its critical functions, and responds to intracellular cues and ecological perturbations through physicochemical modifications. As such, it provides valuable insight into the physiological adaptation of microorganisms in disturbed environmental systems. Here, we (i) demonstrate how changes in salinity and hydraulic retention time (HRT) influence the plasma membrane intact polar lipid (IPL) chemistry of model bacterium, Halanaerobium congolense WG10, and mixed microbial consortia enriched from shale-produced fluids and (ii) elucidate adjustments in membrane IPL chemistry during biofilm growth relative to planktonic cells. We incubated H. congolense WG10 in chemostats under three salinities (7%, 13%, and 20% NaCl), operated under three HRTs (19.2, 24, and 48 h), and in drip flow biofilm reactors under the same salinity gradients. Also, mixed microbial consortia in produced fluids were enriched in triplicate chemostat vessels under three HRTs (19.2, 24, and 72 h) and biofilm reactors. Lipids were analyzed by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Our results show that phosphatidylglycerols, cardiolipins, and phosphatidylethanolamines were predominantly enriched in planktonic H. congolense WG10 cells grown at hypersalinity (20%) compared to optimum (13%). In addition, several zwitterionic phosphatidylcholines and phosphatidylethanolamines were higher in abundance during biofilm growth. These observations suggest that microbial adaptation and biofilm formation in fractured shale are enabled by strategic plasma membrane IPL chemistry adjustments. IMPORTANCE Microorganisms inadvertently introduced into the shale reservoir during fracturing face multiple stressors including brine-level salinities and starvation. However, some anaerobic halotolerant bacteria adapt and persist for long periods of time. They produce hydrogen sulfide, which sours the reservoir and corrodes engineering infrastructure. In addition, they form biofilms on rock matrices, which decrease shale permeability and clog fracture networks. These reduce well productivity and increase extraction costs. Under stress, microbes remodel their plasma membrane to optimize its roles in protection and mediating cellular processes such as signaling, transport, and energy metabolism. Hence, by observing changes in the membrane lipidome of model shale bacteria, Halanaerobium congolense WG10, and mixed consortia enriched from produced fluids under varying subsurface conditions and growth modes, we provide insight that advances our knowledge of the fractured shale biosystem. We also offer data-driven recommendations for improving biocontrol efficacy and the efficiency of energy recovery from unconventional formations.

03 NATURAL GAS↗
DOE OSTI2023

Convolutional Neural Networks for Problems in Transport Phenomena: A Theoretical Minimum

Convolutional neural network (CNN), a deep learning algorithm, has gained popularity in technological applications that rely on interpreting images (typically, an image is a 2D field of pixels). Transport phenomena is the science of studying different fields representing mass, momentum, or heat transfer. Some of the common fields are species concentration, fluid velocity, pressure, and temperature. Each of these fields can be expressed as an image(s). Consequently, CNNs can be leveraged to solve specific scientific problems in transport phenomena. Herein, we show that such problems can be grouped into three basic categories: (a) mapping a field to a descriptor (b) mapping a field to another field, and (c) mapping a descriptor to a field. After reviewing the representative transport phenomena literature for each of these categories, we illustrate the necessary steps for constructing appropriate CNN solutions using sessile liquid drops as an exemplar problem. If sufficient training data is available, CNNs can considerably speed up the solution of the corresponding problems. Finally, the present discussion is meant to be minimalistic such that readers can easily identify the transport phenomena problems where CNNs can be useful as well as construct and/or assess such solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

The bispectrum of cosmic string temperature fluctuations including recombination effects

We calculate the cosmic microwave background temperature bispectrum from cosmic strings, including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. We validate our model for the last scattering surface contribution by comparing the predicted temperature power spectrum with that obtained from a full Boltzmann code treatment applied to the Unconnected Segment Model of a string network. We find that including the last scattering contribution has only a small impact on the upper limit on the string tension resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely to be competitive with the power spectrum at any resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A GPU-accelerated package for simulation of flow in nanoporous source rocks with many-body dissipative particle dynamics

Mesoscopic simulations of hydrocarbon flow in source shales are concerning, in part due to the heterogeneous shale pores with sizes ranging from a few nanometers to a few micrometers. Additionally, the sub-continuum fluid–fluid and fluid–solid interactions in nano- to micro-scale shale pores, which are physically and chemically sophisticated, must be captured. To address those challenges, we present a GPU-accelerated package for simulation of flow in nano- to micro-pore networks with a many-body dissipative particle dynamics (mDPD) mesoscale model. Based on a fully distributed parallel paradigm, the code offloads all intensive workloads on GPUs. Other advancements, such as smart particle packing and no-slip boundary condition in complex pore geometries, are also implemented for the construction and the simulation of the realistic shale pores from 3D nanometer-resolution stack images. Our code is validated for accuracy and compared against the CPU counterpart for speedup. In our benchmark tests, the code delivers nearly perfect strong scaling and weak scaling (with up to 512 million particles) on up to 512 K20X GPUs on Oak Ridge National Laboratory’s (ORNL) Titan supercomputer. Moreover, a single-GPU benchmark on ORNL’s SummitDev and IBM’s AC922 suggests that the host-to-device NVLink can boost performance over PCIe by a remarkable 40%. Lastly, we demonstrate, through a flow simulation in realistic shale pores, that the CPU counterpart requires 840 Power9 cores to rival the performance delivered by our package with four V100 GPUs on ORNL’s Summit architecture. This simulation package enables quick-turnaround and high-throughput mesoscopic numerical simulations for exploring complex flow phenomena in nano- to micro-porous rocks with realistic pore geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗