Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time simulation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Applying Linear Optics from Closed Orbits Modulation for finding beam-based alignment of harmonic sextupoles

A fast and accurate beam-based alignment (BBA) method for harmonic sextupoles has been devel oped at NSLS-II using Linear Optics from Closed Orbit Modulation (LOCOM). The approach excites the beam with simultaneous sine-wave signals at two fast correctors, chosen with an appropriate phase advance to span the full betatron phase space. This strategy suppresses systematic errors from hysteresis, while additional errors from orbit drift and power-supply calibration are minimized by the short measurement time (a few minutes) and reliance solely on beam-based current-to-field conver sion of the sextupoles, with hysteresis explicitly included. Simulations indicate that Linear Optics from Closed Orbits (LOCO) combined with 0.5 mm local orbit bumps can resolve relative sextupole field offsets (∆k₂) with precision better than 10% of k₂, reflecting to beam-based alignment (BBA) accuracy finer than 50 µm. Moreover, employing machine-learning-optimized local orbit bumps en hances sextupole-induced quadrupole signals in a deterministic manner and maintains orbit stability under large sextupole strength variations (±40%). The proposed method is experimentally validated through proof-of-principle measurements at NSLS-II, demonstrating its potential as a fast, precise, and robust tool for harmonic sextupole alignment.

43 PARTICLE ACCELERATORS↗
NASA NTRS1972

Fast differential thermal analysis.

A study is conducted of the limits of time-dependent differential thermal analysis (DTA). Applications of DTA to the hysteresis phenomena of glasses are also reported. The computation of DTA sample temperatures is discussed, giving attention to feasible heating rates, transient effects to be measured, and the simulation of DTA responses to changes in thermal diffusivity similar to hysteresis phenomena in glasses. The reproducibility of temperature recording is checked with zone-refined organic chemicals.

Wunderlich, B.↗
DOE OSTI2025

Analysis and Validation of PMT s Waveforms in ICARUS LArTPC Using Monte Carlo Simulations

ICARUS (Imaging Cosmic and Rare Underground Signals) serves as the Far Detector in the Short Baseline Neutrino (SBN) program at Fermilab, playing a central role in investigating the potential existence of sterile neutrinos in the eV squared mass range. The detector consists of two large Liquid Argon Time Projection Chambers (LArTPCs) with a total capacity of 760 tons of liquid Argon. A key component of the system is its array of 360 Photo-Multiplier Tubes (PMTs), which detect the scintillation light produced by charged particles in liquid Argon; the fast scintillation signal enables accurate event timing, triggering, and reconstruction. Together with the TPC and CRT systems, the PMTs ensure precise interaction time measurements, which are crucial for distinguishing neutrino interactions from cosmic-ray backgrounds. ICARUS uses Hamamatsu R5912-MOD PMTs, optimized for cryogenic temperatures, with high quantum efficiency, excellent timing resolution, low dark current (around 10 nA at 1500 V), and broad spectral sensitivity (300–650 nm). These characteristics are crucial for the efficient detection of scintillation light. Analyzing the waveforms of PMT signals allows for a detailed comparison between experimental data and Monte Carlo simulations. This analysis is fundamental for improving the accuracy of neutrino event reconstruction, enhancing detector calibration, and optimizing the detector's performance for current and future operations.

Brio, V. [Catania U.] (ORCID:0009000088807391)↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Restricted Rotational Flexibility of the C5α-Methyl-Substituted Carbapenem NA-1-157 Leads to Potent Inhibition of the GES-5 Carbapenemase

Carbapenem antibiotics are used as a last-resort treatment for infections caused by multidrug-resistant bacteria. The wide spread of carbapenemases in Gram-negative bacteria has severely compromised the utility of these drugs and represents a serious public health threat. To combat carbapenemase-mediated resistance, new antimicrobials and inhibitors of these enzymes are urgently needed. Here, we describe the interaction of the atypically C5α-methyl-substituted carbapenem, NA-1-157, with the GES-5 carbapenemase. MICs of this compound against Escherichia coli, Klebsiella pneumoniae, and Acinetobacter baumannii producing the enzyme were reduced 4–16-fold when compared to MICs of the commercial carbapenems, reaching clinically sensitive breakpoints. When NA-1-157 was combined with meropenem, a strong synergistic effect was observed. Kinetic and ESI-LC/MS studies demonstrated that NA-1-157 is a potent inhibitor of GES-5, with a high inactivation efficiency of (2.9 ± 0.9) × 105 M –1 s –1 . Acylation of GES-5 by NA-1-157 was biphasic, with the fast phase completing within seconds, and the slow phase taking several hours and likely proceeding through a reversible tetrahedral intermediate. Deacylation was extremely slow (k 3 = (2.4 ± 0.3) × 10 –7 s –1 ), resulting in a residence time of 48 ± 6 days. MD simulation of the GES-5-meropenem and GES-5-NA-1-157 acyl-enzyme complexes revealed that the C5α-methyl group in NA-1- 157 sterically restricts rotation of the 6α-hydroxyethyl group preventing ingress of the deacylating water into the vicinity of the scissile bond of the acyl-enzyme intermediate. Furthermore, these data demonstrate that NA-1-157 is a potent irreversible inhibitor of the GES-5 carbapenemase.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

CaptureFOAM

CaptureFOAM provides utilities for fast and accurate direct air capture (DAC) contactor simulations, based on the OpenFOAM framework. At the core, CaptureFOAM is a suite of mass transfer and chemistry models built upon the OpenFOAM thin film model, which allows coarser mesh resolution and less constrictive time step requirements when compared to a traditional volume of fluid approach. This allows CaptureFOAM to conduct device-scale simulations for complex geometries including commercial packing materials and novel contactor shapes, matching experimental results in a wide variety of applications.

Kincaid, Kellis [Oak Ridge National Laboratory (OR↗
DOE OSTI2021

A Fast Algorithm for Scanning Transmission Electron Microscopy Imaging and 4D-STEM Diffraction Simulations

Scanning transmission electron microscopy (STEM) is an extremely versatile method for studying materials on the atomic scale. Many STEM experiments are supported or validated with electron scattering simulations. However, using the conventional multislice algorithm to perform these simulations can require extremely large calculation times, particularly for experiments with millions of probe positions as each probe position must be simulated independently. Recently, the plane-wave reciprocal-space interpolated scattering matrix (PRISM) algorithm was developed to reduce calculation times for large STEM simulations. Here, we introduce a new method for STEM simulation: partitioning of the STEM probe into “beamlets,” given by a natural neighbor interpolation of the parent beams. This idea is compatible with PRISM simulations and can lead to even larger improvements in simulation time, as well requiring significantly less computer random access memory (RAM). We have performed various simulations to demonstrate the advantages and disadvantages of partitioned PRISM STEM simulations. We find that this new algorithm is particularly useful for 4D-STEM simulations of large fields of view. We also provide a reference implementation of the multislice, PRISM, and partitioned PRISM algorithms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2021

NASA Engineering and Safety Center Technical Bulletin No. 21-02: Genesis Flight Mechanics Simulation

Genesis Flight Mechanics Simulation The NASA Engineering and Safety Center (NESC) consolidated and modernized a suite of legacy flight mechanics simulations, including the Flight Analysis and Simulation Tool (FAST), resulting in Genesis, a generic, multi-vehicle, variable-degree-of-freedom flight mechanics simulation for ascent, aerocapture, entry, descent, and landing (A2EDL) trajectory design. Genesis is more flexible, capable, and performant than FAST. It enables trajectory optimization and interactive trajectory generation. Its interoperability with Copernicus, an exo-atmospheric and interplanetary trajectory design tool, facilitates end-to-end trajectory optimization across all mission phases. Genesis is implemented in Julia, a new language for technical computing that combines the ease of use of scripting languages with the run-time performance of compiled languages.

EDL↗
DOE OSTI2024

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗
NASA NTRS1996

Construction of a Piezoresistive Neural Sensor Array

The construction of a piezoresistive - piezoelectric sensor (or actuator) array is proposed using 'neural' connectivity for signal recognition and possible actuation functions. A closer integration of the sensor and decision functions is necessary in order to achieve intrinsic identification within the sensor. A neural sensor is the next logical step in development of truly 'intelligent' arrays. This proposal will integrate 1-3 polymer piezoresistors and MLC electroceramic devices for applications involving acoustic identification. The 'intelligent' piezoresistor -piezoelectric system incorporates printed resistors, composite resistors, and a feedback for the resetting of resistances. A model of a design is proposed in order to simulate electromechanical resistor interactions. The goal of optimizing a sensor geometry for improving device reliability, training, & signal identification capabilities is the goal of this work. At present, studies predict performance of a 'smart' device with a significant control of 'effective' compliance over a narrow pressure range due to a piezoresistor percolation threshold. An interesting possibility may be to use an array of control elements to shift the threshold function in order to change the level of resistance in a neural sensor array for identification, or, actuation applications. The proposed design employs elements of: (1) conductor loaded polymers for a 'fast' RC time constant response; and (2) multilayer ceramics for actuation or sensing and shifting of resistance in the polymer. Other material possibilities also exist using magnetoresistive layered systems for shifting the resistance. It is proposed to use a neural net configuration to test and to help study the possible changes required in the materials design of these devices. Numerical design models utilize electromechanical elements, in conjunction with structural elements in order to simulate piezoresistively controlled actuators and changes in resistance of sensors. The construction of these devices may show significant improvement in ability to interrogate signals and in the control of effective compliance. This work focuses on the development a variety of series/parallel interconnected piezoresistive control elements for the neural sensing function.

Carlson, W. B.↗
DOE OSTI2020

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.

58 GEOSCIENCES↗
DOE OSTI2024

Overview of the KSTAR experiments toward fusion reactor

The Korean Superconducting Tokamak Advanced Research has been focused on exploring the key physics and engineering issues for future fusion reactors by demonstrating the long pulse operation of high beta steady-state discharge. Advanced scenarios are being developed with the goal for steady-state operation, and significant progress has been made in high ℓ i , hybrid and high beta scenarios with β N of 3. In the new operation scenario called fast ion regulated enhanced (FIRE), fast ions play an essential role in confinement enhancement. GK simulations show a significant reduction of the thermal energy flux when the thermal ion fraction decreases and the main ion density gradient is reversed by the fast ions in FIRE mode. Optimization of 3D magnetic field techniques, including adaptive control and real-time machine learning control algorithm, enabled long-pulse operation and high-performance ELM-suppressed discharge. Symmetric multiple shattered pellet injections (SPIs) and real-time disruption event characterization and forecasting are being performed to mitigate and avoid the disruptions associated with high-performance, long-pulse ITER-like scenarios. Finally, the near-term research plan will be addressed with the actively cooled tungsten divertor, a major upgrade of the NBI and helicon current drive heating, and transition to a full metallic wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mixed Ionic Electronic Conducting Quaternary Perovskites: Materials by Design for Solar Thermochemical Hydrogen

The innovative research conducted by Arizona State University and Princeton University in the project "Mixed Ionic-Electronic Conducting Quaternary Perovskites: Materials by Design for Solar Thermochemical Hydrogen" marks a significant stride forward in thermochemical water splitting. Through an intricate blend of computational design and experimental validation, the project delved into the promising potential of Mixed Ionic Electronic Conducting (MIEC) perovskites. These complex materials, characterized by their unique redox-active nature and adaptability in stoichiometry, present a promising frontier for efficient solar thermochemical hydrogen production. Firstly, the research enhanced the science by utilizing state-of-the-art computational methodologies to unravel the nuanced chemical potentials of MIEC perovskites. By simulating various off-stoichiometric scenarios and redox conditions, the team was able to predict material behaviors under diverse environmental conditions, a feat unachievable through conventional experimental methodologies alone. This approach not only fast-tracks the material screening process, significantly reducing the time from laboratory re-search to practical application, but also uncovers trends and correlations that are pivotal for future materials innovation. Regarding technical effectiveness, the project stands out in its economic feasibility. Traditional methods of materials discovery are often marred by high costs and extensive timeframes, owing to the iterative nature of experimental processes. However, by employing theoretical computations and validating these findings with targeted experiments, the project introduced a cost-effective paradigm for materials discovery and the first ever prediction, synthesis, and preliminary validation of a material solely from computational and theoretical considerations. This synergy between computation and experimentation expedites the discovery of optimal materials conducive to high-efficiency solar-to-hydrogen conversion processes. Furthermore, the public stands to benefit substantially from this research. The success of MIEC perovskites in solar thermochemical applications heralds a shift towards lower cost and lower electricity input for clean hydrogen production, hence potentially impacting climate and energy resilience. By improving the efficiency of solar-to-hydrogen conversions, the research paves the way for reduced dependency on fossil fuels, addressing the urgent global need for accessible and renewable energy sources. Moreover, the project's advancements contribute to scientific literacy in renewable energy technologies, empowering society through knowledge and spurring future innovations. In essence, this research project demonstrates significant progress in the realm of advanced water splitting through solar thermochemistry. Through its groundbreaking approaches in computational materials science and its implications for real-world applications, it holds the promise of a cleaner, more energy-resilient future.

08 HYDROGEN↗
NASA NTRS1975

Digital implementation of the TF30-P-3 turbofan engine control

The standard hydromechanical control modes for TF30-P-3 engine were implemented on a digital process control computer. Programming methods are described, and a method is presented to solve stability problems associated with fast response dynamic loops contained within the exhaust nozzle control. A modification of the exhaust nozzle control to provide for either velocity or position servoactuation systems is discussed. Transient response of the digital control was evaluated by tests on a real time hybrid simulation of the TF30-P-3 engine. It is shown that the deadtime produced by the calculation time delay between sampling and final output is more significant to transient response than the effects associated with sampling rate alone. For the main fuel control, extended update and calculation times resulted in a lengthened transient response to throttle bursts from idle to intermediate with an increase in high pressure compressor stall margin. Extremely long update intervals of 250 msec could be achieved without instability. Update extension for the exhaust nozzle control resulted in a delayed response of the afterburner light-off detector and exhaust nozzle overshoot with resulting fan oversuppression. Long update times of 150 msec caused failure of the control due to a false indication by the blowout detector.

Cwynar, D. S.↗
NASA NTRS1994

Optimizing tertiary storage organization and access for spatio-temporal datasets

We address in this paper data management techniques for efficiently retrieving requested subsets of large datasets stored on mass storage devices. This problem represents a major bottleneck that can negate the benefits of fast networks, because the time to access a subset from a large dataset stored on a mass storage system is much greater that the time to transmit that subset over a network. This paper focuses on very large spatial and temporal datasets generated by simulation programs in the area of climate modeling, but the techniques developed can be applied to other applications that deal with large multidimensional datasets. The main requirement we have addressed is the efficient access of subsets of information contained within much larger datasets, for the purpose of analysis and interactive visualization. We have developed data partitioning techniques that partition datasets into 'clusters' based on analysis of data access patterns and storage device characteristics. The goal is to minimize the number of clusters read from mass storage systems when subsets are requested. We emphasize in this paper proposed enhancements to current storage server protocols to permit control over physical placement of data on storage devices. We also discuss in some detail the aspects of the interface between the application programs and the mass storage system, as well as a workbench to help scientists to design the best reorganization of a dataset for anticipated access patterns.

Chen, Ling Tony↗
DOE OSTI2024

Thermal performance of phase change materials with anisotropic carbon fiber inserts

In thermal energy storage systems, phase change materials (PCMs) are widely used for thermal energy management. Most PCMs have low thermal conductivities, which limits the heat transfer rate within PCM and thus makes the phase-changing process very slow. However, thermal conductivities of PCMs can be altered by inserting targeted additives. We hypothesize that the size and shape of these additive inserts play a key role in thermal management efficiency. To this end, the impacts of carbon fiber (CF) inserts on the phase change behavior and consequent heat transfer efficiencies of inorganic and organic PCMs were investigated using experiments and simulations. Long, anisotropic CFs with high thermal conductivities formed continuous fast heat flux tunnels inside PCMs to enhance the heat transfer. Such CFs could extend the phase change fronts from the limited container-shaped interface to the larger surface of numerous CF inserts inside the PCM. These special CF inserts work with a new heat transfer mechanism, different from conventional small additives or long isotropic CF inserts. The thermal energy release rate increased by 2.5 times with 1 wt.% anisotropic CF inserts in inorganic PCM. However, CF inserts in liquid organic PCM hindered the natural convection and compromised the improved heat conduction. In conclusion, the lab-scale multiphysics simulations support these experimental observations and indicate that CF inserts have potential to enhance heat transfer in inorganic PCMs, but they are less effective in organic PCMs.

25 ENERGY STORAGE↗
DOE OSTI2021

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗