Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock finding”

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.

442 records · Page 25

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

X-Ray Computed Tomography of Apollo Drive Tube 73001 as Part of the Apollo Next Generation Sample Analysis Program

Introduction: Starting in 2019, the Apollo Next Generation Sample Analysis (ANGSA) Program has enabled consortium studies of specially curated Apollo samples that were previously unstudied (or under studied). This program was intended as a participating scientist program for these Apollo samples, as well as a preparation for upcoming lunar sample return missions (e.g., Artemis). The work began with studies of unsealed core tube 73002 [1,2] that is the upper ~35 cm of an Apollo 17 double drive tube (~70 cm total depth) collected on the landslide deposit near Lara Crater at Station 3 in the Taurus Littrow Valley. More recently the program extended to the study of a variety of frozen Apollo 17 samples [3], as well as the gas extraction [4] and dissection [5] of 73001, the lower half of the station 3 double drive tube. Sample 73001 was sealed under vacuum on the Moon in a Core Sample Vacuum Container (CSVC). In this abstract we will examine the role of X-ray Computed Tomography (XCT) during the preliminary examination process for sealed core 73001, including: (1) engineering scans to aid in understanding the gas extraction process, (2) whole-core scanning prior to opening to inform extrusion and dissection work, and (3) individual particle scanning to characterize rock fragment lithologies for follow on studies. Methodology: Sample 73001 is a 33 cm long, 4 cm diameter regolith sample collected inside a lunar drive tube (~1 mm aluminum walls). That drive tube was sealed inside a 0.5 mm thick stainless steel (SS) CSVC. XCT scans for engineering purposes were done on the Nikon XTH 320 system at Johnson Space Center using the 225 kV multi-metal reflection source at 215 kV, 179 A, and a 38.49 m voxel size. Individual >4 mm particles separated from the core during processing (then triply sealed in Teflon bags) were also scanned at JSC using the 180 kV source at 90 kV, 33 A, and a 2.98 – 10.65 m voxel size. Whole-core scans were done at the University of Texas High-Resolution X-ray Computed Tomography Facility (UTCT) on the 225 kV reflection source on the North Star Imaging cabinet XCT system. These scans included: (1) a series of 9 overlapping super-resolution scans each covering a ~4 cm length of the tube at 190 kV, 180 µA, and a 12.9 µm voxel size and (2) a lower resolution continuous helical scan of the entire core at 190 kV, 180 µA, and a 51.8 µm voxel size. Progress and Results: Before piercing and extracting the gas from sample 73001, an XCT scan of the bottom portion of the CSVC was used to confirm the location of the Teflon cap on the inner drive tube and whether there was an “air” gap between the Teflon cap and the CSVC. This was to ensure the Teflon cap (which immobilized the soil inside the drive tube) was not accidentally pierced during gas extraction. Similarly, after piercing, the bottom and top portions of the CSVC were scanned in order to capture engineering knowledge about the results of the piercing process, as well as the metal knife edge vacuum seal (SS into an In-Ag alloy). Both scans will provide constraints on future work of this type, both for the design of tools used on the lunar surface during Artemis, as well as how we open future samples of this type. An unexpected finding from these engineering scans was that the device in the drive tube that immobilizes the regolith (the keeper) was not seated in the tube properly. This meant that (1) the drive tube could not be removed from the CSVC for the trip to UTCT, and (2) the procedure for opening and extruding the drive tube had to be modified. Had either of these things not been known prior to opening the CSVC, it could have led to an inability to XCT scan the whole core at high resolution and/or potential disruption of the core stratigraphy during extrusion. At UTCT, the entire length of the core was scanned at high resolution (12.9 microns per voxel). This scan serves multiple purposes: (1) A lower resolution (and uncorrected) version of these scans stitched together was used to help inform the processors of potential pitfalls during extrusion and dissection (voids, angular fragments at the edges of the sample, etc.); and (2) the full resolution corrected data will serve as the permanent in situ record of the stratigraphy of the sample and will enable future researchers to perform a variety of analyses. All 121 >4 mm particles separated during dissection pass 1 of sample 73001 have been individually scanned. These scans clearly show the lithology of each particle while keeping the particles in pristine condition (they are individually triply bagged in Teflon within our N2 purged processing cabinets). Because of the dust adhering to particle exteriors, it would otherwise be impossible to determine lithologies in a non-contaminating way. Thus far the types of lithologies seen in sample 73001 (e.g., regolith breccias, impact-melt breccias, agglutinates, and bas-alts) are similar to those previously identified in sample 73002 [2]. By the time of the meeting, all particles from all 3 dissection passes will have been scanned and statistics on the different lithologies in 73001 compiled. References: [1] Shearer et al. (2020) 51st LPSC, abstract 1181. [2] Zeigler et al. (2020) 51st LPSC, abstract 3023. [3] Kent et al. (2022) This Volume. [4] McDonald et al. (2022) European Lunar Symposium. [5] Gross et al (2022) This volume.

Moon↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no or some addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4 rocket engine fired on 1m2 area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon. The resulting computed values of plume impingement surface temperature, stagnation pressure, gas velocity, shear stress and heat flux were then matched as closely as possible in the Earth’s atmosphere in a sub-scale rocket engine GO2/CH4 test which was provided by Masten Space Systems in Mojave, California. The rocket engine was mounted on a test stand with vertical translation capabilities so that the landing operations of a lander could be simulated (Figure 1).The pad materials test coupons were placed at the surface of a large bin containing simulated lunar basalt regolith and subjected to a test firing as shown in Figure 2. The results of this testing will be presented with related findings, analysis and discussion.

Plume Surface Interaction↗

Identification And Survival Of Bacteriohopanepolyol In A Hot Spring Microbial Mat

The polar lipids of a hot spring microbial mat located in Yellowstone National Park were examined for the presence of bacteriohopanepolvols (BHP). BHP are a group of molecules consisting of a hopanoid (peotacyclic triterpene) linked via a n-alkyl polyhydroxylated chain to a variety of polar end groups. BHP have been isolated in varying amounts from phylogenetically diverse eubacterial groups including cyanobacteria, methanotrophs and the Rhodospirillaceae. The hopanoids are excellent biomarkers and have been detected in sedimentary rocks as old as 1.7 bya. In order to interpret the ancient organic record, it is important to understand the abundance, source and fate of such biomarker compounds in microbial mats. A 40 sq cm mat section was taken from a 52 to 55 C site in the effluent channel of Octopus Spring and was sampled vertically over approximately 16 mm. The first 5-6 mm was sectioned into a top green layer (310 mg dry weight) and several subjacent, deep orange layers (240 and 250 mg, respectively). The lower 10 mm of the mat was sectioned into two gelatinous orange layers containing a siliceous gritty material (260 and 440 mg) which increased with depth, and a bottom layer composed almost exclusively of siliceous sinter (4.1 g). The progressive decrease in total organic carbon from 45% in the top green layer to only 4% in the bottom layer reflects the observed increase in siliceous deposition. GC-MS analysis of the phospholipid and glycolipid fatty acids yielded predominantly saturated normal chain acids, n-15 to n-18, and iso-branched acids, i-15 to i-17. Small amounts of unsaturated fatty acids (16:1, two positional isomers of 18:1, and two cyclopropyl acids, C(sub 17) and C(sub 19)) were present mainly in the top layer. Esterified fatty acid which is a good index for intact cellular membrane, i.e. viable organisms, was highest in the top two layers (203 and 231 micro g/mg total lipid, respectively) and gradually decreased to 66 micro g/mg total lipid in the bottom layer. Small amounts of BHP were present in all six layers, however in this case, BHP was lowest in the top green and subjacent deep-orange layers (118 and 172 micro g/mg total lipid, respectively) and increased with depth reaching almost 400 micro g/mg in the bottom two layers. This data suggest that BHP are survivina the initial phase of mat degradation and may be preferentially enriched in any organic record of such thermal environments. The relatively low level of BHP in the top layer also suggests that cyanobacteria may not be the major source of BHP in this mat. Since Chloroflexus a major component of the deep-orange layer has been reported to lack BHP, this material may prove a valuable biomarker for some other mat inhabitant. Further isotopic characterization of this BHP should help resolve this finding.

Janke, Linda L.↗

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↗

Ar-39-Ar-40 Ages of Euerites and the Thermal History of Asteroid 4-Vesta

Eucrite meteorites are igneous rocks that derive from a large asteroid, probably 4 Vesta. Prior studies have shown that after eucrites formed, most were subsequently metamorphosed to temperatures up to equal to or greater than 800 C, and much later many were brecciated and heated by large impacts into the parent body surface. The uncommon basaltic, unbrecciated eucrites also formed near the surface but presumably escaped later brecciation, whereas the cumulate eucrites formed at depth where metamorphism may have persisted for a considerable period. To further understand the complex HED parent body thermal history, we determined new Ar-39-Ar-40 ages for nine eucrites classified as basaltic but unbrecciated, six eucrites classified as cumulate, and several basaltic-brecciated eucrites. Relatively precise Ar-Ar ages of two cumulate eucrites (Moama and EET87520) and four unbrecciated eucrites give a tight cluster at 4.48 +/1 0.01 Gyr. Ar-Ar ages of six additional unbrecciated eucrites are consistent with this age, within their larger age uncertainties. In contrast, available literature data on Pb-Pb isochron ages of four cumulate eucrites and one unbrecciated eucrite vary over 4.4-4.515 Gyr, and Sm-147 - Nd-143 isochron ages of four cumulate and three unbrecciated eucrites vary over 4.41-4.55 Gyr. Similar Ar-Ar ages for cumulate and unbrecciated eucrites imply that cumulate eucrites do not have a younger formation age than basaltic eucrites, as previously proposed. Rather, we suggest that these cumulate and unbrecciated eucrites resided at depth where parent body temperatures were sufficiently high to cause the K-Ar and some other chronometers to remain open diffusion systems. From the strong clustering of Ar-Ar ages at approximately 4.48 Gyr, we propose that these meteorites were excavated from depth in a single large impact event approximately 4.48 Gyr ago, which quickly cooled the samples and started the K-Ar chronometer. A large (approximately 460 km) crater postulated to exist on Vesta may be the source of these eucrites and of many smaller asteroids thought to be spectrally or physically associated with Vesta. Some Pb-Pb and Sm-Nd ages of cumulate and unbrecciated eucrites are consistent with the 4.48 Gyr Ar-Ar age, and the few older Pb-Pb and Sm-Nd ages may reflect isotopic closure prior to the large cratering event. One cumulate eucrite gives an Ar-Ar age of 4.25 Gyr; three additional cumulate eucrites give Ar-Ar ages of 3.4-3.7 Gyr; and two unbrecciated eucrites give Ar-Ar ages of approximately 3.55 Gyr. We attribute these younger ages to later impact heating. In addition, we find Ar-Ar impact-reset ages of several brecciated eucrites and eucritic clasts in howardites to fall in the range of 3.5-4.1 Gyr. Among these, Piplia Kalan, the first eucrite to show evidence for extinct 26 Al, was strongly impact heated approximately3.5 Gyr ago. When these data are combined with eucrite Ar-Ar ages in the literature, they confirm the previous suggestion that several large impact heating events occurred on Vesta over the time period approximately 4.1-3.4 Gyr ago. The onset of major impact heating may have occurred at similar times for both Vesta and the Moon, but impact heating appears to have persisted to a somewhat later time on Vesta compared to the Moon.

Bogard, Donald D.↗

Development and Field Testing Novel Natural Gas Surface Process Equipment for Replacement of Water as Primary Hydraulic Fracturing Fluid (Final Report)

Southwest Research Institute® (SwRI ® ), Schlumberger Technology Corporation (SLB), and Chevron Corporation ® (Chevron) developed a novel, optimized, and lightweight process for natural gas (NG) to replace water as a low-cost fracturing medium with a low environmental impact. Hydraulic fracturing is used to increase oil and NG production by injecting high-pressure fluid, primarily water, into a rock formation, which fractures the rock and releases trapped oil and NG. This method was developed to increase yield and make feasible production areas that would not otherwise be viable for large-scale oil and NG extraction using traditional drilling technologies. Since the fracturing fluid is composed of approximately 90% water, one of the principal drawbacks to hydraulic fracturing is its excessive water use and associated large environmental footprint. According to recent data collected at fracture sites within the United States, fracturing applications in North America can consume as much as 9 million gallons of water per well. During the fracturing process, some of the fracturing fluid is permanently lost and the portion that is recovered is contaminated by both fracturing chemicals and dissolved solids from the formation. The recovered water or flow-back represents a significant environmental challenge, as it must be treated before it can be reintroduced into the natural water system. Although there is some recycling for future fracturing, the majority of the flow-back water is hauled from the well site to a treatment facility or to an injection well for permanent underground disposal. To mitigate these issues, an optimized, lightweight and modular surface process using NG to replace a majority of the water was developed as a cost-effective and environmentally clean fracturing fluid. Using NG will result in significantly less consumption since the gas that is injected as a fracturing fluid will be mixed with the formation gas and extracted as if it were from the formation itself. This process will minimize the collection, waste, and treatment of large amounts of water and will reduce the environmental impact of transporting and storing the fracturing fluid. There are two major steps involved in utilizing NG as the primary fracturing medium: (1) increasing the supply pressure of NG to wellhead pressures suitable for fracturing and (2) mixing the required chemicals and proppant needed for the fracturing process at these elevated pressures. The second step (NG-proppant mixing at elevated pressures) still requires demonstration; but very similar processes have been demonstrated in the field with other gases such as nitrogen (N 2 ) and carbon dioxide (CO 2 ). However, the first step (a compact, on-site unit for generating high-pressure NG at costs feasible for fracturing) has not been developed and is currently not commercially available. Due to the inherent compressibility of NG, more energy is required to compress the gas than what is required to pump water (or other incompressible liquids) to the very high pressure required for downhole injection. This project aimed to develop a novel, hybrid method to overcome this challenge. The project accomplishments and findings are discussed in this report. Ultimately, the research and development efforts described herein demonstrate that fracturing with NG foam is a feasible alternative to using water.

02 PETROLEUM↗

Structural, Electrochemical, and (De)lithiation Mechanism Investigation of Cation-Disordered Rocksalt and Spinel Hybrid Nanomaterials in Lithium-Ion Batteries

Significant demand for lithium-ion batteries necessitates alternatives to Co- and Ni-based cathode materials. Cation-disordered materials using earth-abundant elements are being explored as promising candidates. Here, in this paper, we demonstrate a coprecipitation synthetic approach that allows direct preparation of disordered rocksalt Li 2.4 Fe 1.0 Ti 1.0 O 4.7 (r-LFTO·C) and spinel structured hybrid Li 0.5 Fe 1.0 Ti 0.9 O 3.2 ·C (s-LFTO·C) nanoparticles with a conformal conductive carbon coating. High-angle annular dark-field imaging coupled with electron energy loss spectroscopy mapping shows uniform Fe/Ti distribution with minor compositional variation among particles. Cation disorder was confirmed for both of the materials at an atomic level, with a short-range order more pronounced in r-LFTO·C. Operando X-ray absorption spectroscopy, ex situ hard X-ray photoelectron spectroscopy, ex situ soft X-ray absorption spectroscopy, and ex situ synchrotron X-ray diffraction were used to investigate (de)lithiation in the bulk and at the surface. Structurally, the r-LFTO·C demonstrated reversible partial Fe center migration between octahedral and tetrahedral sites during (de)lithiation. The r-LFTO·C evidenced that the redox of O was coincident with iron redox during initial electrochemical cycling, while iron redox dominated later cycling. In contrast, s-LFTO·C electrochemistry involved iron redox throughout the cycling process. The findings rationalize the differences in the electrochemistry where r-LFTO·C shows higher initial capacity yet poorer capacity retention over a voltage window where O redox can be accessed, while the s-LFTO·C shows lower initial capacity yet improved capacity retention.

25 ENERGY STORAGE↗

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗