Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “angle control”

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 613 records · Page 34

Ring Pull Strain Analysis

Ring Pull Strain Analysis (RPSA) is a framework for analysis of Digital Image Correlation (DIC) data of gaugeless ring pull testing. It is meant for analysis of rings cut from tube cross sections and pulled in a displacement-controlled tensile test. This code takes load frame data, images take during testing, and csv files from the DIC analysis of those images and allows the user to analyze and plot data from the test. Some important features of the software include getting the strain values as a function of radial distance and angle, plotting of the strain around the ring, and stress-strain analysis of the load frame data.

Beck, Peter↗

Ring Pull Strain Analysis

Ring Pull Strain Analysis (RPSA) is a framework for analysis of Digital Image Correlation (DIC) data of gaugeless ring pull testing. It is meant for analysis of rings cut from tube cross sections and pulled in a displacement-controlled tensile test. This code takes load frame data, images take during testing, and csv files from the DIC analysis of those images and allows the user to analyze and plot data from the test. Some important features of the software include getting the strain values as a function of radial distance and angle, plotting of the strain around the ring, and stress-strain analysis of the load frame data. Additional capability is included with Ring Pull Coating Analysis (RPCA, also part of the software). Additional images are taken of samples with an outer coating and are used to correlate failure of the coating with strain values measured with DIC.

Beck, Peter↗

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Unusual Ferromagnetic Band Evolution and High Curie Temperature in Monolayer 1T‐CrTe2 on Bilayer Graphene

2D van der Waals ferromagnets hold immense promise for spintronic applications due to their controllability and versatility. Despite their significance, the realization and in-depth characterization of ferromagnetic materials in atomically thin single layers, close to the true 2D limit, has been scarce. Here, a successful synthesis of monolayer (ML) 1T-CrTe2 is reported on a bilayer graphene (BLG) substrate via molecular beam epitaxy. Using angle-resolved photoemission spectroscopy and magneto-optical Kerr effect measurements, that the ferromagnetic transition is observed at the Curie temperature (TC) of 150 K in ML 1T-CrTe2 on BLG, accompanied by unconventional temperature-dependent band evolutions. The spectroscopic analysis and first-principle calculations reveal that the ferromagnetism may arise from Goodenough-Kanamori super-exchange and double-exchange interactions, enhanced by the lattice distortion and the electron doping from the BLG substrate. These findings provide pivotal insight into the fundamental understanding of mechanisms governing 2D ferromagnetism and offer a pathway for engineering higher TC in 2D materials for future spintronic devices.

Park, Kyoungree↗

Data source authentication of synchrophasor measurement devices based on 1D-CNN and GRU

Synchrophasor measurement devices (SMDs) have been widely deployed to support real-time monitoring and control of power systems. In the meantime, data spoofing has emerged in recent years. Therefore, it is of great importance to study data authentication algorithms for detecting and defending the data spoofing effectively. Here, a one-dimensional convolutional neural network (1D-CNN) is utilized to extract temporal signatures hidden in frequency, voltage angle and amplitude data; then the gated recurrent unit (GRU) employs these temporal signatures for data source authentication. In case studies, the performances of different algorithms are tested in large-scale power systems with numerous SMDs for the first time, and comparisons among different algorithms show that the proposed algorithm can achieve a higher accuracy of data source authentication with a shorter time window.

47 OTHER INSTRUMENTATION↗

Polyphenols Weaken Pea Protein Gel by Formation of Large Aggregates with Diminished Noncovalent Interactions

Polyphenols are well-known native cross-linkers and gel strengthening agents for many animal proteins. However, their role in modifying plant protein gels remains unclear. In this study, multiple techniques were applied to unravel the influence of green tea polyphenols (GTP) on pea protein gels and the underlying mechanisms. We found that the elasticity and viscosity of pea protein gels decreased with increased GTP. The protein backbone became less rigid when GTP was present based on shortened T 1ρ H in relaxation solid-state NMR measurements. Electron microscopy and small-angle X-ray scattering showed that gels weakened by GTP possessed disrupted networks with the presence of large protein aggregates. Solvent extraction and molecular dynamic simulation revealed a reduction in hydrophobic interactions and hydrogen bonds among proteins in gels containing GTP. The current findings may be applicable to other plant proteins for greater control of gel structures in the presence of polyphenols, expanding their utilization in food and biomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗

Solid-State Sintering of Core-Shell Ceramic Powders Fabricated by Particle Atomic Layer Deposition

The properties of technical ceramics are highly dependent on their microstructure, which evolves during sintering. Sintering is the process by which ceramic parts are subjected to high temperatures to activate chemical diffusion and the consumption of porosity. During the initial stage of sintering, interparticle necks between neighboring particles form and subsequently increase in size, consuming porosity as the particle centers move closer together. To experimentally determine how this process depends on particle surface composition, particle atomic layer deposition (ALD) was used to deposit a thin film of amorphous aluminum oxide (Al2O3) onto yttria-stabilized tetragonal zirconia (3YSZ) particles, producing core-shell structured powders. The uniformity of the Al2O3 film was confirmed with transmission electron microscopy and energy dispersive spectroscopy. Scanning electron microscopy was used to observe microstructural evolution during sintering, and the dihedral angles of Al2O3 and 3YSZ grains were measured to determine the ratio of interfacial energies between the 3YSZ|3YSZ, 3YSZ|Al2O3, and Al2O3|Al2O3 interfaces. Analysis of the densification kinetics revealed that the initial stage of densification is dependent on the material at the surface of the particles (ie, the Al2O3 film) and is controlled by the diffusion of Al3+ cations through Al2O3. Once the Al2O3 film has coalesced, the sintering behavior is controlled by the densification of the core material (3YSZ). Thus, core-shell powders fabricated by particle ALD sinter by a two-step process where the kinetics are dependent on the material present at interparticle contacts.

alumina↗

Dataset_for_Molecular_Motion_Below_the_Glass_Transition_A_Solid-State_NMR_Study_of_Siloxane_Polymer_Dynamics Study

This dataset contains solid-state 1H and 13C NMR relaxometry data, differential scanning calorimetry (DSC) data, and size exclusion chromatography (SEC/GPC) data supporting the study of sub-glass-transition (sub-Tg) molecular dynamics in a composition- and sequence-controlled series of diphenyl-substituted polysiloxanes (PDMS, 14Ph, 33Ph, 50Ph, 67Ph, and 100Ph; 0–100% diphenylsiloxane content by mole).All solid-state NMR data were acquired on a 200 MHz Bruker Avance III HD spectrometer using a static 7 mm HX probe or a 4 mm HX probe under 4 kHz magic-angle spinning. Raw Bruker TopSpin experiment folders are included for: (1) variable-temperature 1H lineshape measurements used to determine linewidth (FWHM) as a function of temperature across the glass transition; (2) 1H T1 (saturation recovery with solid-echo detection), probing nanosecond-scale dynamics near the 1H Larmor frequency; (3) 1H T1rho (direct spin-lock, 62.5 kHz), probing microsecond-scale segmental dynamics; (4) 13C-detected Lee–Goldburg cross-polarization 1H T1rho (LGCPH T1rho) for 33Ph and 50Ph, resolving aromatic and aliphatic proton environments; and (5) 13C T1 relaxation for 33Ph and 50Ph. Differential scanning calorimetry data (TA Instruments DSC 25, −150 to +120 °C, up to +300 °C for 100Ph, 10 °C/min) are included for all six compositions and support the glass-transition temperatures in Table 1 and Figure 1. Size exclusion chromatography data (Agilent 1200 Series, PL-Gel 300 mixed-C column, THF mobile phase, polystyrene calibration standards) are included for the three synthesized copolymers (33Ph, 50Ph, 67Ph) and support the number-average molecular weights in Table 1. Processed data include per-composition relaxation-time summaries (Excel), curve-fitting and Bloembergen-Purcell-Pound (BPP) model analysis notebooks (Jupyter/Python), and Igor Pro (.pxp) master files used to generate the manuscript's figures.

Bloembergen-Purcell-Pound theory↗

Significance of the Direct Excitation Mechanism for High-Frequency Response of Premixed Flames to Flow Oscillations

Premixed flames are sensitive to flow disturbances, which can arise from acoustic or vortical fluctuations. For transverse instabilities, it is known that a dominant mechanism for flame response is “injector coupling,” whereby pressure oscillations associated with transverse waves excite axial flow disturbances. These axial flow disturbances then excite heat release oscillations. The objective of this paper is to consider another mechanism—the direct sensitivity of the unsteady heat release to transverse acoustic waves—and to compare its significance relative to the induced axial disturbances, in a linear framework. The rate at which the flame adds energy to the disturbance field is quantified using the Rayleigh criterion and evaluated over a range of control parameters, such as flame length and swirl number. The results show that radial modes induce heat release fluctuations that always add energy to the acoustic field, whereas heat release fluctuations induced by mixed radial-azimuthal modes can add or remove energy. These amplification rates are then compared to the flame response from induced axial fluctuations. For combustor-centered flames, these results show that the direct excitation mechanism has negligible amplification rates relative to the induced axial mechanism for radial modes. For transverse modes, the fact that the nozzle is located at a pressure node indicates that negligible induced axial velocity disturbances are excited; as such, the direct mechanism dominates. As a result, for flames that are not centered on pressure nodes, the direct mechanism for mixed modes dominates for certain nozzle locations and flame angles.

42 ENGINEERING↗

Understanding interfacial crystallization dynamics on carbon fiber reinforced polypropylene composite manufacturing

Reinforcing polymers with discontinuous fibers improves mechanical properties, such as strength and stiffness, and in some cases achieve isotropic properties, rendering them suitable for various engineering applications. Matrix materials are generally highly engineered thermosets (e.g. crosslinked epoxies), bonded to the fiber periphery by proprietary surface and sizing chemistries. Semicrystalline thermoplastic matrices are less utilized due to poor fiber-matrix bonding resulting in inefficient interfacial load-transfer in reinforced composites. However, flexibility with melt-processing or molding conditions can be leveraged to promote non-covalent interfacial bonding between matrix and fiber via crystallization of the matrix onto fiber surface. In the present study, we utilize a co-mingle chopped carbon and isotactic polypropylene fibers to form isotropic composites, tailoring interfacial immobilized matrix or interphase morphology to optimize performance through precise control of thermal processing/molding windows. Calorimetry and optical microscopy were employed to investigate the impact of carbon fiber at various volume fractions (10, 20, and 30 %) on isotactic polypropylene crystallization and mechanical performance. Variations in mechanical properties correspond to the structural evolution of the interfacial region and are correlated to underlying microstructural attributes using wide-angle X-ray scattering, thermal analysis, and low-field nuclear magnetic resonance spectroscopy. These results provide a practical framework for the manufacturing of thermoplastic matrix composites. In conclusion, the results presented provide a guide for the strategic optimization of interphase design, showcasing tailorable tensile strengths which outperform any isotactic polypropylene carbon fiber composites previously reported in literature.

36 MATERIALS SCIENCE↗

Snapshots of Ce 70 Toroid Assembly from Solids and Solution

Crystallization at the solid–liquid interface is difficult to spectroscopically observe and therefore challenging to understand and ultimately control at the molecular level. The Ce70-torroid formulated [Ce IV 70 (OH) 36 (O) 64 (SO 4 ) 60 (H 2 O) 10 ] 4– , part of a larger emerging family of M IV 70 -materials (M = Zr, U, Ce), presents such an opportunity. We elucidated assembly mechanisms by the X-ray scattering (small-angle scattering and total scattering) of solutions and solids as well as crystallizing and identifying fragments of Ce 70 by single-crystal X-ray diffraction. Fragments show evidence for templated growth (Ce 5 , [Ce 5 (O) 3 (SO 4 ) 12 ] 10- ) and modular assembly from hexamer (Ce 6 ) building units (Ce 13 , [Ce 13 (OH) 6 (O) 12 (SO 4 ) 14 (H 2 O) 14 ] 6– and Ce 62 , [Ce 62 (OH) 30 (O) 58 (SO 4 ) 58 ] 14- ). Ce 62 , an almost complete ring, precipitates instantaneously in the presence of ammonium cations as two torqued arcs that interlock by hydrogen boding through NH 4 + , a structural motif not observed before in inorganic systems. The room temperature rapid assemblies of both Ce 70 and Ce 62 , respectively, by the addition of Li + and NH 4 + , along with ion-exchange and redox behavior, invite exploitation of this emerging material family in environmental and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics of Neutrino Oscillations and Neutrino Cross Sections using the NOvA Experiment

NOvA is a long-baseline accelerator-based neutrino experiment based in the US. NOvA uses an intense neutrino beam produced at Fermilab’s accelerator complex to make physics measurements of neutrino oscillations, neutrino cross sections, and much more. For its physics goals, NOvA uses two functionally-identical detectors. The Near Detector (ND) is situated at Fermilab, 1 km from the neutrino target and the Far Detector (FD) is located at Ash River, MN, a distance of 810 km from the neutrino source. The ND receives a high statistics neutrino flux which gives a unique opportunity for high-precision neutrino cross-section measurements and is used as a control for the oscillation analyses. The FD is used to analyze the appearance and disappearance of the neutrinos arriving from the Fermilab. The purpose of the oscillation analysis is to understand the dominance of matter over antimatter in the universe, to resolve the ordering of neutrino masses, and to resolve the octant of the neutrino mixing angle theta23. In this talk, I will give an overview of the NOvA experiment. I will also talk about the status of the NOvA’s cross-section physics program and the latest results from the oscillation analyses. I will also talk about the future prospects of the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗

In-Plane Critical Magnetic Fields in Magic-Angle Twisted Trilayer Graphene

It has recently been shown that superconductivity in magic-angle twisted trilayer graphene survives to in-plane magnetic fields that are well in excess of the Pauli limit, and much stronger than the in-plane critical magnetic fields of magic-angle twisted bilayer graphene. The difference is surprising because twisted bilayers and trilayers both support the magic-angle flat bands thought to be the fountainhead of twisted graphene superconductivity. In this study, we show that the difference in critical magnetic fields can be traced to a C 2 M h symmetry in trilayers that survives in-plane magnetic fields, and also relative displacements between top and bottom layers that are not under experimental control at present. An gate electric field breaks the C 2 M h symmetry and therefore limits the in-plane critical magnetic field.

2-dimensional systems↗

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Hybrid-Learning Algorithm for Online Dynamic State Estimation in Multimachine Power Systems

With the increasing penetration of distributed generators in the smart grids, having knowledge of rapid real-time electromechanical dynamic states has become crucial to system stability control. Conventional Supervisory Control and Data Acquisition (SCADA)-based dynamic state estimation (DSE) techniques are limited by the slow sampling rates, while the emerging phasor measurement units (PMUs) technology enables rapid real-time measurements at network nodes. Using generator bus terminal voltages, we propose a hybrid-learning DSE (HL-DSE) algorithm to estimate the synchronous machine rotor angle and speed in real time. The HL-DSE takes the power system model into account and trains neuroestimators with real-time data in an online manner. Compared with traditional DSE methods, the HL-DSE overcomes limitations by using a data-driven approach in conjunction with the physical power system model. The time efficiency, accuracy, convergence, and robustness of the proposed algorithm are tested under noises and fault conditions in both small- and large-scale test systems. Simulation results show that the proposed HL-DSE is much more computationally efficient than widely used Kalman filter (KF)-based methods while maintaining comparable accuracy and robustness. In particular, HL-DSE is over 100 times faster than square-root unscented KF (SR-UKF) and 80 times faster than extended KF (EKF). The advantages and challenges of the HL-DSE are also discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Vapor-phase grafting of functional silanes on atomic layer deposited Al 2 O 3

Fundamental studies are needed to advance our understanding of selective adsorption in aqueous environments and develop more effective sorbents and filters for water treatment. Vapor-phase grafting of functional silanes is an effective method to prepare well-defined surfaces to study selective adsorption. In this investigation, we perform vapor phase grafting of five different silane compounds on aluminum oxide (Al2O3) surfaces prepared by atomic layer deposition. These silane compounds have the general formula L3Si–C3H6–X where the ligand, L, controls the reactivity with the hydroxylated Al2O3 surface and the functional moiety, X, dictates the surface properties of the grafted layer. We study the grafting process using in situ Fourier transform infrared spectroscopy and ex situ x-ray photoelectron spectroscopy measurements, and we characterize the surfaces using scanning electron microscopy, atomic force microscopy, and water contact angle measurements. Here, we found that the structure and density of grafted aminosilanes are influenced by their chemical reactivity and steric constraints around the silicon atom as well as by the nature of the anchoring functional groups. Methyl substituted aminosilanes yielded more hydrophobic surfaces with a higher surface density at higher grafting temperatures. Thiol and nitrile terminated silanes were also studied and compared to the aminosilane terminated surfaces. Uniform monolayer coatings were observed for ethoxy-based silanes, but chlorosilanes exhibited nonuniform coatings as verified by atomic force microscopy measurements.

36 MATERIALS SCIENCE↗

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↗