Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic parameter”

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

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

At least 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Probing Surface Plasmon Dynamics in Periodic Nanostructures through Ultrafast Electron Microscopy

Surface plasmon polaritons (SPPs) can be manipulated to localize and guide light in subwavelength distances, enabling them to find applications in a wide range of areas, from sensing to quantum computing. Among several methods of SPP excitation, periodic arrays of nano- and microstructures are of particular interest, as they enable engineering SPP properties through structural parameters. Here, in this study, using the photon-induced near-field electron microscopy (PINEM) technique, we investigated the mode formation, coupling, interference, and decay of SPPs in square and hexagonal arrays of circular nanoholes under both visible and near-infrared excitation. Polarization-resolved analysis revealed the key factors governing SPP localization and interference patterns, showing that the periodicity and symmetry of the array primarily determine the SPP interference patterns and their orientation, while pump polarization mainly modulates their intensity. Time-resolved PINEM measurements demonstrated the spatial dependence of the SPP temporal characteristics. In addition, cathodoluminescence (CL) spectroscopy was employed to examine the intrinsic plasmonic characteristics of the structure. Finite difference time domain (FDTD) simulations showed strong agreement with both PINEM and CL measurements on the spatial and spectral behavior of SPPs. Understanding the spatiotemporal dynamics of SPPs on nanostructures beyond the diffraction limit is crucial for optimizing plasmonic structures for advanced photonic and quantum technologies.

Plasmonics↗
DOE OSTI2020

Reconstruction of effective potential from statistical analysis of dynamic trajectories

The broad incorporation of microscopic methods is yielding a wealth of information on the atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here, we develop a method for the stochastic reconstruction of effective local potentials solely from observed structural data collected from molecular dynamics simulations (i.e., data analogous to those obtained via atomically resolved microscopies). Using the silicon vacancy defect in graphene as a model, we apply the statistical framework presented herein to reconstruct the free energy landscape from the calculated atomic displacements. Evidence of consistency between the reconstructed local potential and the trajectory data from which it was produced is presented, along with a quantitative assessment of the uncertainty in the inferred parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Multichannel dynamics in the OH+ n -butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Additionally, experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

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

Stone Content Influence on Land Surface Model Simulation of Soil Moisture and Evapotranspiration at Reynolds Creek Watershed

Considerable advancement in spatiotemporal resolution of remote sensing and ground-based measurements has enabled refinement of parameters used in land surface models for simulating surface water fluxes. However, land surface modeling capabilities are still inadequate for accurate representation of subsurface properties and processes, which continue to limit the accuracy of land surface model simulation. Our objective in this study was to examine the performance of the variously parameterized Noah-Multiphysics (Noah-MP) land surface model and to improve simulations of evapotranspiration (ET) and soil moisture dynamics in stony soils using verification from eddy covariance ET and in-situ soil moisture data. Meteorological and eddy covariance data from the Lower Sheep sub-catchment within the Reynolds Creek Experimental Watershed in Southwestern Idaho were employed to drive Noah-MP and to assess simulation results. We evaluated the performance of Noah-MP considering four different scenarios: (1) A one-layer soil profile with Noah-MP default soil hydraulic parameters, followed by five-layer soil profiles using -(2) Noah-MP default soil hydraulic parameters -(3) soil hydraulic parameters derived from a pedotransfer function using field observations, and -(4) hydraulic parameters from (3) which also accounted for stone content in each layer. Each modeling experiment was forced with the same set of initial conditions, atmospheric input and vegetation parameters. Our results indicated that enhanced representation of soil profile properties and stone content information noticeably improve the Noah-MP land surface model simulation of soil moisture content and evapotranspiration.

Land surface models, Soil hydrology, Noah-MP, Evap↗
DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
DOE OSTI2021

Seasonal changes in the drivers of water physico-chemistry variability of a small freshwater tidal river

Where rivers meet the sea, tides can exert a physical and chemical influence on the lower reaches of a river. How tidal dynamics in these tidal river reaches interact with upstream hydrological drivers such as storm rainfall, which ultimately determines the quantity and composition of material transferred from watersheds to estuaries, is currently unknown. We monitored a small freshwater tidal river in the Pacific Northwest, USA in high resolution over one year to evaluate the relative importance of tides versus upstream hydrological flows (i.e., base flow and precipitation events) on basic physico-chemical parameters (pH, dissolved oxygen, turbidity, specific conductivity, and temperature), and how these interactions relate to the downstream estuary. Tidal variability and diurnal cycles (i.e. solar radiation) dominated water physico-chemical variability in the summer, but the influence of these drivers was overshadowed by storm-driven sharp pulses in river physico-chemistry during the remainder of the year. Within such events, we found incidences of counterclockwise hysteresis of pH, counterclockwise hysteresis of dissolved oxygen, and clockwise hysteresis of turbidity, although systematic trends were not observed across events. The dominance of storm rainfall in the river’s physico-chemistry dynamics, and similar pulses of decreased pH observed in adjacent estuarine waters, suggest that the linkage between tidal streams and the broader system is variable throughout the year. High-frequency monitoring of tidal river biogeochemistry is therefore crucial to enable the assessment of how the relative strength of these drivers may change with future sea level rise and altered precipitation patterns to modulate biogeochemical dynamics across the land-ocean-atmosphere continuum.

aquatic↗
DOE OSTI2021

Scalable modular dynamic molten salt reactor system model with decay heat

Nodal dynamic models allow rapid reduced-order simulations of complex systems. This paper presents a publicly available system dynamic model of a molten salt reactor system with a topology of the Molten Salt Reactor Experiment (MSRE), in MATLAB/Simulink, building upon a model previously published and validated with MSRE data. It adds scaling of nominal power, dynamic representation of decay heat generation, a generic decay heat removal system, improved documentation, modularity, and test cases. One hour of system transient runs in less than 10 min on a laptop. This model is intended for first order engineering calculations, developing insights into system transient behavior, equipment sizing, parameter sensitivity studies. The model contains an easy to use toolkit which can be customized and expanded to describe any molten salt reactor system. Finally, scenarios involving off-normal transients are presented. Relation between shutdown reactivity insertions, decay heat removal, core recriticality, and safe shutdown conditions are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Synaptic Functionality and Neuromorphic Information Processing in Membrane Ion Channel Junctions

The human brain performs complex memory and computational tasks with high energy efficiency by regulating ion transport through membrane channels. These signaling mechanisms have been inspiring the development of nanofluidic memristors that emulate synaptic behavior. Here, in this study, we describe a membrane ion channel synapse (MICS), constructed from aqueous droplets linked by gramicidin A channels, that achieves neuromorphic functionality. MICS exhibits memristive ion transport with hysteretic current–voltage behavior arising from voltage-dependent channel formation and ion transport dynamics. MICS emulates a range of synaptic behaviors including associative learning. We further demonstrate its application in reservoir computing by performing handwritten digit classification and tic-tac-toe game and explore the system parameters that improve the computational performance. This droplet-based biomimetic synapse offers a potentially scalable and energy-efficient platform for next-generation neuromorphic computing systems.

Droplet interface bilayer↗
DOE OSTI2020

Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains

Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Automated model generation and parameter estimation of building energy models using an ontology-based framework

This study presents a methodology for automated model generation and parameter estimation of building energy models using semantic modeling and Bayesian estimation. Semantic modeling techniques are used to represent the system components and their interactions, facilitating the automatic generation of a simulation model from dynamic component models. The proposed approach is applied to a case study of a ventilation system where a simulation model is generated, calibrated, and assessed through different performance metrics. These metrics demonstrate the accuracy and reliability of both model point estimates and probabilistic prediction intervals across all model outputs. Overall, the proposed methodology offers a systematic and automated approach to model development and calibration in building energy systems, with potential applications in building performance analysis, monitoring, and optimization.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A Scale‐Adaptive Urban Hydrologic Framework: Incorporating Network‐Level Storm Drainage Pipes Representation

Abstract Below‐ground urban stormwater networks (BUSNs) significantly influence urban flood dynamics, yet their representation at the watershed or larger scales remains challenging. We introduce a scalable urban hydrologic framework that centers on a novel network‐level BUSN representation, balancing the needs for physical basis, parameter parsimony, and computational efficiency. Our framework conceptualizes an urban watershed into four interacting zones: hillslopes (natural), storm‐sewersheds (urban), a sub‐network channel (tributaries), and a main channel. We develop an innovative Graph Theory‐based algorithm to derive network‐level BUSN parameters from publicly available datasets, enabling efficient, scalable parameterization. We demonstrate this framework's applicability at nine representative watersheds in the Houston metropolitan region, USA, with urban imperviousness ranging from 0% to 64% and drainage areas ranging from 24 to 302 . Our model achieves satisfying computational efficiency, completing hourly time step simulations for 18 years in less than 5 sec per watershed on a standard PC. Validation against observed daily streamflow confirms that the model can capture small‐to‐large flood peaks and seasonal and annual water balance over these watersheds. Comparisons with the National Water Model show better performance in predicting flood peaks and overall water balance, underscoring the promises of our new framework for urban hydrologic modeling at large scales. Furthermore, analysis reveals nonlinear relationships between BUSNs' designed capacities and flood reduction effects. Our approach bridges the gap between detailed hydraulic and large‐scale hydrologic models, providing a valuable tool for urban flood prediction and management across broader spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Distributed Small-Signal Stability Conditions for Inverter-Based Unbalanced Microgrids

The proliferation of inverter-based generation and advanced sensing, controls, and communication infrastructure have facilitated accelerated deployment of microgrids. A coordinated network of microgrids can maintain reliable power delivery to critical facilities during extreme events. Low-inertia offered by the power-electronics–interfaced energy resources, however, can present significant challenges to ensuring stable operation of the microgrids. In this work, distributed small-signal stability conditions for inverter-based microgrids are developed that involve the droop-controller parameters and the network parameters (e.g. line impedances, loads). The distributed closed-form parametric stability conditions derived in this paper can be verified in a computationally efficient manner, facilitating reliable design and operations of networks of microgrids. Dynamic phasor models have been used to capture the effects of electromagnetic transients. Furthermore, numerical results are presented, along with PSCAD simulations, to validate the analytical stability conditions. Effects of design choices, such as the conductor types, and inverter sizes, on the small-signal stability of inverter-based microgrids are investigated to derive useful engineering insights.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Utah FORGE: Pump and Probe Test on a Mated Fracture Westerly Granite Sample

This dataset contains results from a pump and probe experiment conducted on a mated fracture Westerly Granite sample with a diameter of 1 inch and a height of 2 7/8 inches. The experiment was performed within an aluminum triaxial pressure vessel (TEMCO) to investigate the non-linear acoustic parameters of the sample under varying axial pressures. The confining pressure was maintained at 4 MPa, while the axial pressure was incrementally reduced from 17 MPa to 1 MPa in 1 MPa steps. 5 dynamic pressure oscillations of 0.5 MPa were applied over a 10 minute interval. The dataset includes pump controller data, linear variable differential transformer (LVDT) measurements, and acoustic data. The pump controller data tracks the confining and axial pressures, flow rates, and pump volumes. LVDT measurements provide detailed records of the axial displacement of the cell piston. Additionally, acoustic data captured by s-wave transducers is archived in a compressed file.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pilot Testing of a Modular Oxygen Production System Using Oxygen Binding Adsorbents (Final Report)

In this project, RTI and Air Liquide have focused on the development of innovative oxygen separation materials and technologies based on reversibly binding of oxygen to enable smaller and cheaper air separations. In alignment with that object the team has had multiple achievements to advance the technology including: Synthesized and characterized novel M-CoSalen (RTIO2Sorb-1) material which showed a dynamic oxygen adsorption capacity in excess of 1wt%. Developed synthesis routes for RTIO2Sorb-1 with commercially relevant techniques and produced batches of 0.25kg and over 4 kg of total synthesis. Synthesized and characterized extrudate forms of the RTIO2Sorb-1 with relevant crush strength and maintained dynamic oxygen adsorption capacity in near 1wt%. Developed techniques for forming structured beds in fiber shapes with conventional materials. Developed a model for O 2 sorption processes to analyze various bed configurations and cycle parameters. Studied classic O 2 solid sorbents reported in the literature to learn oxygen binding mechanisms. Design a 10 kg/day O 2 VPSA pilot system with 2 or 4-bed operation. Completed a techno-economic analysis based on the experimental results and modeling. The overall objective of the project was to design, fabricate, and test a modular O 2 production system and perform a techno-economic analysis (TEA) after testing to determine the cost-benefit of the advanced modular air separation system. The goals of this technology development project were to achieve a bed-size factor (BSF) of less than 600 lb-adsorbent/TPD O 2 (ton/day O 2 ) (as compared with the state-of-the-art BSF of 850), and O 2 -purity greater than 95% at a cost that is projected to be equivalent or lower than the current state of the art (SOTA), commercially available large-scale cryogenic air separation systems. To achieve these goals and objective, the project team executed on (1) oxygen binding adsorbent optimization and scale up, (2) adsorbent material formation process studies to form the adsorbent material into structured beds for rapid pressure swing adsorption (PSA) cycles with low pressure drop, fast mass transfer, and low attrition, (3) cycle development studies to optimize the PSA process, and (4) develop simulation tools for rapid cycle modeling and numerical evaluation/optimization.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

A New Vehicle-to-Vehicle Communication Technology Based on Binary Light Code

The proliferation of Connected and Autonomous Vehicles (CAVs) has necessitated the development of efficient, reliable communication methods between vehicles and their surroundings. This study presents a novel methodology, the VECTOR system, that addresses this need by converting dynamic vehicle data into a binary code and displaying it on an LED panel. The system involves a three-step process of data collection, polynomial fitting of velocity data, and encoding the polynomial parameters into binary form using a Cyclic Redundancy Check (CRC). The study also explores the practical application of this system by conducting experiments with modified Lincoln MKZ hybrids, demonstrating the feasibility and efficiency of this approach. The paper presents both the methodology and results of these experiments, further expanding upon the potential applications and implications for CAV technology. By comparing the final accuracy, we find this approach can achieve 80% of the original velocity data.

Ma, Ke↗
DOE OSTI2021

Spin Change of Asteroid 2012 TC4 Probably by Radiation Torques

Asteroid 2012 TC4 is a small (∼10 m) near-Earth object that was observed during its Earth close approaches in 2012 and 2017. Earlier analyses of light curves revealed its excited rotation state. We collected all available photometric data from the two apparitions to reconstruct its rotation state and convex shape model. We show that light curves from 2012 and 2017 cannot be fitted with a single set of model parameters; the rotation and precession periods are significantly different for these two data sets, and they must have changed between or during the two apparitions. Nevertheless, we could fit all light curves with a dynamically self-consistent model assuming that the spin states of 2012 TC4 in 2012 and 2017 were different. To interpret our results, we developed a numerical model of its spin evolution in which we included two potentially relevant perturbations: (i) gravitational torque due to the Sun and Earth and (ii) radiation torque, known as the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect. Despite our model simplicity, we found that the role of gravitational torques is negligible. Instead, we argue that the observed change of its spin state may be plausibly explained as a result of the YORP torque. To strengthen this interpretation, we verify that (i) the internal energy dissipation due to material inelasticity and (ii) an impact with a sufficiently large interplanetary particle are both highly unlikely causes of its observed spin state change. If true, this is the first case where the YORP effect has been detected for a tumbling body.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗