Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex mixtures”

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.

750 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A structural underpinning of the lower critical solution temperature (LCST) behavior behind temperature-switchable liquids

In this work, we use state-of-the-art X-ray scattering and molecular dynamics to analyze amine-water mixtures that show the unusual lower critical solution temperature (LCST) behavior. The goal is to provide direct experimental evidence for the entropy-lowering molecular cluster formation hypothesized as necessary for LCST behavior. Differential wide-angle X-ray scattering and pair distribution analysis and small-angle X-ray scattering measurements were combined with molecular modeling and liquid-liquid equilibrium measurements, revealing direct experimental evidence for the hypothesis. However, the response of the amine phase to accommodating water is even more subtle than the simple hypothesis suggests, with the formation of robust nanoscale reverse micelles. The techniques developed in this paper can be expected to yield insights in the use of temperature-switchable liquids in solvent extraction and other separations, and the stabilization of organelles in living cells that do not have physical membranes but do require compositional gradients to operate.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the role of earth alkaline doping in perovskite-based methane dry reforming catalysts

To elucidate the role of earth alkaline doping in perovskite-based dry reforming of methane (DRM) catalysts, we embarked on a comparative and exemplary study of a Ni-based Sm perovskite with and without Sr doping. While the Sr-doped material appears as a structure-pure Sm 1.5 Sr 0.5 NiO 4 Ruddlesden Popper structure, the undoped material is a NiO/monoclinic Sm 2 O 3 composite. Hydrogen pre-reduction or direct activation in the DRM mixture in all cases yields either active Ni/Sm 2 O 3 or Ni/Sm 2 O 3 /SrCO 3 materials, with albeit different short-term stability and deactivation behavior. The much smaller Ni particle size after hydrogen reduction of Sm 1.5 Sr 0.5 NiO 4 , and of generally all undoped materials stabilizes the short and long-term DRM activity. Carbon dioxide reactivity manifests itself in the direct formation of SrCO 3 in the case of Sm 1.5 Sr 0.5 NiO 4 , which is dominant at high temperatures. For Sm 1.5 Sr 0.5 NiO 4 , the CO : H 2 ratio exceeds 1 at these temperatures, which is attributed to faster direct carbon dioxide conversion to SrCO 3 without catalytic DRM reactivity. As no Sm 2 O 2 CO 3 surface or bulk phase as a result of carbon dioxide activation was observed for any material – in contrast to La 2 O 2 CO 3 – we suggest that oxy-carbonate formation plays only a minor role for DRM reactivity. Rather, we identify surface graphitic carbon as the potentially reactive intermediate. Graphitic carbon has already been shown as a crucial reaction intermediate in metal-oxide DRM catalysts and appears both for Sm 1.5 Sr 0.5 NiO 4 and NiO/monoclinic Sm 2 O 3 after reaction as crystalline structure. It is significantly more pronounced for the latter due to the higher amount of oxygen-deficient monoclinic Sm 2 O 3 facilitating carbon dioxide activation. Despite the often reported beneficial role of earth alkaline dopants in DRM catalysis, we show that the situation is more complex. In our studies, the detrimental role of earth alkaline doping manifests itself in the exclusive formation of the sole stable carbonated species and a general destabilization of the Ni/monoclinic Sm 2 O 3 interface by favoring Ni particle sintering.

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

IDAES-PSE 2.4.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications.. Deprecations • Convergence Analysis tool (idaes/core/util/convergence): deprecated in favor of new Parameter Sweep tools. To be removed in v3.0.0. New Beta Capabilities • Parameter Sweep Tool (idaes.core.util.parameter_sweep) o A new API for defining and performing parameter sweep studies on IDAES models has been developed • Diagnostics Tools (idaes.core.util.model_diagnostics) o New methods for identifying duplicate variables and constraints have been added to the diagnostics toolbox o New tools for detecting ill conditioning in Jacobians have been developed and are available in the model_diagnostics module. These provide alternatives to the existing DegeneracyHunter toolbox, and will eventually be merged with this capability, but initial working versions have been provided as beta capabilities for interested users o IpoptConvergenceAnalysis (replaces deprecated Convergence Analysis tool):  A new tool for performing convergence analysis studies that leverages the new Parameter Sweep tools has been developed. This tool allows users to define the input parameters to their model and sampling methods for these (leveraging Pysmo's sampling tools) and to then solve their model across the sampled domains and return a summary of the solver performance (IPOPT only) Improved Models • Thickener model (idaes.models.unit_models.solid_liquid.thickener) o Improved model to include predictive correlations for unit sizing based on settling velocity measurements (steady-state only) • Modular Property Packages o Added general support for calculating critical properties of mixtures using defined Equation of State modules. New API defined for Equation of State modules in order to define the necessary constraints for calculating critical properties (most EoS modules DO NOT support calculation of critical properties (yet)) o Added new methods to Cubic Equation of State module to support calculation of critical properties

DiagnosticsToolbox↗
DOE OSTI2020

Development of a Testbed for Pipeline Flushing - 20425

Pipelines which carry high-level waste within the DOE complex should be properly flushed in the cases where stationary or moving beds of solid sediment occur, or lines are prone to hydrogen gas buildup. Present guidelines establish a minimum for the amount of water (flush volume) and flush velocity values used for post-transfer flushing operations to achieve satisfactory cleaning of pipelines. However, further studies are needed to find optimal operating velocity modes/values in order to minimizing the flush volume and consequent downstream waste. These efforts are significantly helpful to DOE waste remediation sites by preserving tank storage, preventing additional waste processing, and minimizing dilution and changes in waste chemistry. An experimental test loop was recently developed for study of non-Newtonian slurry flushing at the Florida International University (FIU). This loop was designed to create sediment beds of various materials and bed heights and investigate parameters that effect the efficiency of flushing operations. The objective is to find flush velocity values/modes which lead to satisfactory cleaning of transport lines with a minimum amount of water usage. A 165 ft experimental test loop made of 3-inch carbon steel pipes was constructed and used in initial testing campaigns. This pipe loop was equipped with real-time monitoring instruments to record variation of pressure, density, mass flow rate, and temperature during testing. Post-flush in-situ no-flow evaluations were performed using ultrasonic, endoscopy, and visual inspection. Additional post-flush evaluations were conducted using filtration and density monitoring of residuals within a designated circulations loop. These provisions enable comparison of pipeline cleanness before the start and after termination of flushing operations for efficiency evaluation. This paper presents efforts associated with flushing of kaolin-water mixtures at different concentrations and with various flush volume and flush velocity values/modes. In initial testing, efforts focused on creating repeatable sediment beds inside the pipeline in both fully-flooded (no pre-flush drainage) and gravity-drained (pre-flush drained) systems. Bed characterization (i.e., measurement of sediment height and solids concentration) was conducted to ensure about consistency of initial conditions between tests. Preliminary results for flushing study of 10% vol. kaolin-water mixtures starting from a fully-flooded initial condition showed that flushing with 1-line volume (65 gallons) could remove up to 72 percent of solids from the pipeline. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Rhenium-Sulfido and -Dithiolato Corroles: Reflections on Chalcophilicity

The high-temperature (~180 °C) reaction between free-base meso-triarylcorroles and Re 2 (CO) 10 , followed by exposure to PCl 3 and thiols (or elemental sulfur), affords rhenium-sulfido (ReS) corroles in 67-76% yields. The use of shorter reaction times, lower temperatures (~130 °C), and a dithiol (e.g., ethane-1,2-dithiol) also allows the isolation of rhenium-dithiolato corroles, presumptive intermediates on the path to ReS corroles. The ReS corroles exhibit high thermal stability and two reversible oxidations and reductions in their cyclic voltammograms, with redox potentials nearly identical to those observed for analogous ReO corroles. The electrochemical HOMO-LUMO gaps of the complexes, at 2.2 eV, are consistent with ligand-centered oxidation and reduction. The UV-vis spectra of ReS corroles, on the other hand, differ significantly from those of their ReO counterparts. Scalar-relativistic DFT calculations suggest that this difference reflects low-energy LUMO+2 and LUMO+3 levels, consisting of Re-S π-antibonding interactions; the ReO corroles, in contrast, exhibit a larger LUMO+1/LUMO+2 gap, as expected for a relatively classical Gouterman-type metalloporphyrin analogue. The high stability of ReS corroles is consistent with geochemists' view of rhenium as a moderately chalcophilic element (i.e., one that partitions into sulfide melts) as well as with a recent quantitative analysis of thiophilicity, which indicates that rhenium's oxophilicity and thiophilicity are essentially evenly balanced.

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

Quantification of Residual Water in Spent Fuel Dry Storage Canisters Using Mass Spectrometry

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal, such as fuel degradation; cladding corrosion, embrittlement, or breaching; and the creation of a flammable environment via radiolysis of water. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents details on the quantification of residual water in the Advanced Drying Cycle Simulator (ADCS), an apparatus built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. The ADCS is outfitted with thermocouples to measure the thermal response of the ADCS to simulated decay heats and internal helium fill pressures relevant to commercial drying procedures. The ADCS is also instrumented with pressure transducers to measure the pressures and vacuum levels observed during simulated commercial drying. The most unique instrumentation used for quantifying residual water in the ADCS is a Hiden Analytical HPR-30 mass spectrometer (MS), which measures gas compositions of the ADCS internal free volume, based on partial pressures calculated from relative proportions of gas molecules detected by the MS. This report details the methodology used to implement MS measurements in quantifying residual water in the ADCS. This methodology includes the calibration of the HPR-30 MS to a Buck Research Instruments CR-4 chilled mirror hygrometer, which itself is calibrated to a NIST-traceable standard. Data collected by both the MS and the chilled mirror hygrometer from water/helium mixtures ranging from 150 to 500,000 ppmv water in helium were used to generate calibration curves, establishing a source of verification of MS measured water contents. Details regarding water content measurement uncertainties are included in this report, defining the accuracy and verifiability of the HPR-30 MS in measuring residual water content in simulated dry storage canister environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Standoff Detection of Chemical Plumes Using Swept-ECQCL Remote Detection Technology: Results from EMRTC Complex Terrain Dual Tracer Experiment

As part of the Complex Terrain Dual Tracer Experiments performed in support of the LYNM program, a swept-wavelength external cavity quantum cascade laser (swept-ECQCL) remote detection system was deployed to measure chemical plume propagation with high temporal resolution. The swept-ECQCL remote detection technology was initially developed at PNNL, and is now under further development at Opticslah, LLC under SBIR funding from DOE (DOE SBIR Phase I Project DE-SC0019855, Remote Detection Technologies, Program Manager Chris Ramos). The swept-ECQCL technology provides a single-frequency infrared laser source which is rapidly scanned over a large wavelength range to measure the infrared absorption of gases, liquids, or solids. Analysis of the measured infrared absorption spectrum is used to detect, identify, and quantify the chemical species which are present. While most previous laser-based techniques used small wavelength tuning ranges to measure only a few isolated spectral lines of a few species, the swept-ECQCL provides the large tuning range needed to measure multiple species simultaneously, and to detect broad absorption features from large molecules or condensed phase materials. In addition, the swept-ECQCL systems have high enough spectral resolution and sensitivity to measure narrow lines of small molecules and determine isotope ratios. The remote detection system for gas plumes is based on the swept-ECQCL technology, and consists of a swept-ECQCL source, infrared detector, control electronics, and beam transmitting and alignment optics. The remote detection system directs the infrared output of the swept-ECQCL to a target retro reflector typically located at ~10-10,000 m distance. The swept-ECQCL wavelength is varied continuously over a large wavelength range (>1000 nm) at rates up to 1000 Hz. The swept ECQCL light reflected from the remote target is collected and focused onto an infrared photodetector. By measuring the detector signal continuously as the swept-ECQCL wavelength is varied, a time-series of absorption spectra are measured for the gases along the beam line of-sight. Based on analysis and fitting of the measured absorption spectra, the time-dependent chemical concentrations of species are determined. Multiple species and mixtures can be detected simultaneously, for all species with absorption features above the sensor noise floor. The sensor output is thus a continuous, high-speed record of identified chemicals along with their concentrations and isotope ratios of interest. The high-speed operation enables low-noise operation in conditions of atmospheric turbulence and allows tracking of fluctuating concentrations in transient chemical plumes as needed for reliable remote detection applications. All components operate with minimal cooling requirements and no cryogens are needed. In addition, the output infrared light is invisible with an intensity below the maximum permissible exposure (MPE) threshold of 100 mW/cm2 for these infrared wavelengths and thus eye safe. The measurement concept for plume detection using the swept-ECQCL remote detection system is as follows. The ECQCL beam is directed toward a remote retro-reflector, which reflects the beam back to the ECQCL system where it is focused onto an infrared photodetector. The line-of-sight from the ECQCL system to the retro-reflector defines the measurement path. The detected infrared signal is recorded continuously and detects the absorption of molecular species along the measurement path. Before a plume release, the system measures absorption from atmospheric constituents, which in this LWIR band consist primarily of lines from H 2 O and CO 2 , with additional weaker absorption from N 2 O, NH 3 , and O 3 possibly contributing. The system then measures any changes in absorption due to the tracer compounds when they enter the measurement path. The shape of absorption features is used to distinguish and identify the different tracer compounds, and the strength of the absorption is used to determine the path-integrated concentration – also called the column density.

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

High-resolution rovibrational spectroscopy of trans -formic acid in the v 1 OH stretching fundamental: Dark state coupling and evidence for charge delocalization dynamics

High-resolution infrared (IR) reduced-Doppler absorption spectra of jet-cooled gas phase trans-formic acid in the v 1 OH stretching fundamental region are reported for the first time, obtained by supersonically expanding trans- formic acid/Ar mixtures through a slit jet nozzle source and rotationally cooling to T rot ≈ 10.9(5) K, with ab- sorption signals recorded by high-resolution difference-frequency IR absorption spectroscopy. Two a/b-type rovibrational bands of comparable intensity, one ~10-fold weaker b-type band, and one ~6-fold weaker a-type band are observed, with vibrational band origins at 3570.493(5), 3566.793(5), 3560.032(9), and 3534.6869(2) cm –1 , respectively. Based on previous Raman jet spectroscopic work by Nejad and Sibert [A. Nejad, E.L. Sibert III, The Raman jet spectrum of trans-formic acid and its deuterated isotopologs: Combining theory and experi- ment to extend the vibrational database, J. Chem. Phys. 154(6) (2021) 064301.], these four rovibrational bands have been assigned to v 1 , (v 2 + v 7 ), (v 6 + 2v 7 + 2v 9 ), and 2v 3 , respectively. Specifically, two of the three upper dark states (2 1 7 1 (a') and 6 1 7 2 9 2 (a')) are close enough to the “bright” 1 1 (a') state to facilitate strong anharmonic resonance interactions, which results in intensity mixing into the two zero-order bands that would otherwise be “dark”. Furthermore, our high-resolution spectral analysis reveals that there are local rotational crossings be- tween these zero-order 1 1 and 2 1 7 1 states resulting in extra lines (i.e., some upper levels in the nominally v 1 band have majority zero-order 2 1 7 1 state character). This motivates development of a 3 coupled state (1 1 , 2 1 7 1 , and 6 1 7 2 9 2 ) picture to aid in the spectral analysis, which is able to match all 3 observed band origins and relative band intensities, as well as indicate the necessity of multistate (> 2) coupling. Though limited by range of J and Ka levels (J’ ≤ 9 and K a ’ ≤ 3) populated at supersonic jet temperatures, this work offers first precision spec- troscopic analysis of trans-formic acid in the v 1 OH stretch region, which should aid in assignment of the more complete yet highly congested room temperature FTIR spectra. Lastly, and in sharp contrast to the spectral complexity in the three predominantly b-type bands, the lone a-type 2v 3 rovibrational band at 3534.6869(2) cm –1 is well described by a simple, rigid asymmetric top Hamiltonian.

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

Synthesis of Hf 6 Ta 2 O 17 superstructure via spark plasma sintering for improved oxidation resistance of multi-component ultra-high temperature ceramics

Ultra-high temperature ceramics (UHTCs) have shown aspiration to overcome challenges in the thermal protection system (TPS) by designing new materials referred to as multi-component UHTCs (MC-UHTCs) in the compositional space. MC-UHTCs have shown remarkable improvement in oxidation resistance due to the formation of the Hf6Ta2O17 superstructure during plasma exposure. Herein, the Hf 6 Ta 2 O 17 superstructure is synthesized via a solid-state reaction between HfO 2 and Ta 2 O 5 powder mixtures during spark plasma sintering (SPS). The compositions chosen are 50 vol% of HfO 2 -50 vol% of Ta 2 O 5 (50HO-50TO) and 70 vol% of HfO 2 -30 vol% of Ta 2 O 5 (70HO-30TO). The phase quantification via Rietveld analysis showed Hf 6 Ta 2 O 17 as a principal phase with some residual Ta 2 O 5 phase in both the samples. The high-temperature thermal stability of the samples was evaluated using high-velocity plasma jet exposure for up to 3 min. 50HO-50TO was able to withstand the intense plasma condition, which is attributed to the higher content of the Hf 6 Ta 2 O 17 phase (~84%) and lower strain in the Ta 2 O 5 phase. The augmentation in the Hf 6 Ta 2 O 17 phase to 94.7% (in 50HO-50TO) post plasma exposure has been attributed to the invariant transformation from a liquid state to Hf 6 Ta 2 O 17 at temperatures >2500 °C during testing. The mechanical integrity is elucidated from the insignificant change in the hardness ~13.3 GPa before and 11.2 GPa after plasma exposure of the 50HO-50TO sample. As a result, the Hf 6 Ta 2 O 17 superstructure's thermo-mechanical stability suggests developing novel oxidation-resistant MC-UHTCs in compositional space for reusable space vehicle applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Curious Case of [AnH(NR 2 ) 3 ] (An = Th, U; R = SiMe 3 ): Two Monomeric Actinide Hydrides Revisited

The reaction of AnCl 4 (DME) x (An = Th, x = 2; An = U, x = 0) with 4 equiv of NaNR 2 (R = SiMe 3 ) in THF at 65 °C results in the formation of [An{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (An = U, 1; An = Th, 2), and not the reported monomeric actinide hydrides, [AnH(NR 2 ) 3 ], as expected. Furthermore, both complexes 1 and 2 were characterized by X-ray crystallography. Surprisingly, their unit cell parameters are remarkably close to those reported for [AnH(NR 2 ) 3 ], suggesting that the original crystals of [AnH(NR 2 ) 3 ] were, in fact, [An{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ], but were misidentified. Reduction of 1 with 1.1 equiv of KC 8 in THF, in the presence of 1 equiv of 2.2.2-cryptand, results in the formation of [K(2.2.2-cryptand)][U{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (3) in good yield. Likewise, the reaction of 1 with 1 equiv of bis(diisopropylamino)cyclopropenylidene (BAC) results in the formation of the BAC adduct, [(BAC)U{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (4), in moderate yield. Finally, the addition of H 2 (10 bar) to 2 in C 6 D 6 at room temperature results in the formation of the targeted monomeric hydride, [ThH(NR 2 ) 3 ], in 32% yield, according to integrations against an internal standard. However, removal of the H 2 atmosphere results in rapid reformation of 2. In contrast, the addition of H 2 (10 bar) to 1 in C 6 D 6 at room temperature results in no apparent reaction.

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

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
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↗