Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Commissioning of a modulated pulse-power magnetron sputtering system for depositing niobium thin films

A modulated pulse-power magnetron sputtering system (MPPMS) was commissioned and used to deposit niobium on copper. MPPMS utilizes a modulated pulse train with high peak-to-average energy ratios, thereby increasing the energy density leading to the production of the energetic ions. The MPPMS modulator was designed and constructed based on large voltage/current IGBT modules, which can be externally modulated, produce bi-polar waveforms, and are readily available. A high voltage power supply provides moderate current to the IGBT module to create the precise pulse structure. The MPPMS system was used to deposit ~1 mm thick, niobium film on a copper substrate held at ~77 K to minimize Nb atom surface diffusion and maximize intrinsic film energy. Finally, the Nb film had an average grain size of ~8 nm.

42 ENGINEERING↗
DOE OSTI2021

Environmental life cycle impacts of small wastewater treatment plants: Design recommendations for impact mitigation

The objective of this study was to quantify potential mitigation of environmental impacts from the operation and construction of wastewater treatment plants (WWTP) from implementing specific design recommendations. The study investigated small WWTPs, many of which are serving slow growing or declining populations. Life Cycle Assessment methodology was used to evaluate and compare the inventory and environmental impacts of nine small WWTP case studies. Detailed inventory data was collected from the facilities’ engineering design plans and utility bills. One recommended practice was to avoid significant overdesign by planning for no lower than a 75% capacity utilization by the facilities’ end-of-life. A theoretical correction to a 75% capacity utilization was estimated to mitigate 0.4% of lifetime electricity usage and 1% of secondary process concrete for every 1% reduction in design average flow rate. Relatedly, a 0.4% mitigation in the Carcinogenic and Global Warming impacts could be achieved for every 1% reduction in design average flow toward a 75% capacity utilization. Other suggested practices were focused on conveyance, namely, to minimize non-process facility area and to use polyvinyl chloride pipe instead of ductile iron pipe where possible. Furthermore, the latter practice was estimated to mitigate between 1.1 and 4.8% of the Carcinogenic impact in the nine case studies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Desalinating a Real Hyper-Saline Pre-Treated Produced Water via Direct-Heat Vacuum Membrane Distillation

Membrane distillation (MD) is an emerging thermal desalination technology capable of desalinating waters of any salinity. During typical MD processes, the saline feedwater is heated and acts as the thermal energy carrier; however, temperature polarization (as well as thermal energy loss) contributes to low distillate fluxes, low single-pass water recovery and poor thermal efficiency. An alternative approach is to integrate an extra thermal energy carrier as part of the membrane and/or module assembly, which can channel externally provided heat directly to the membrane-feedwater interface and/or along the feed channel length. This direct-heat delivery has been demonstrated to increase single-pass water recovery and enhance the overall thermal efficiency. We developed a bench-scale direct-heated vacuum MD (DHVMD) process to desalinate pre-treated oil and gas "produced water" with an initial total dissolved solids of 115,500 ppm at a feed temperature ranging between 24 and 32 degrees C. We evaluated both water flux and specific energy consumption (SEC) as a function of water recovery. The system achieved a 50% water recovery without significant scaling, with an average flux >6 kg m-2 hr-1 and a SEC as low as 2,530 kJ kg-1. The major species of mineral scales (i.e., NaCl, CaSO4, and SrSO4) that limited the water recovery to 68% were modeled in terms of thermodynamics and identified by scanning electron microscopy and energy-dispersive X-ray spectroscopy. In addition, we further developed and employed a physics-based process model to estimate temperature, salinity, water transport and energy flows for full-scale vacuum MD and DHVMD modules. Model results show that a direct-heat input rate of 3,600 W can increase single-pass water recovery from 2.1% to 3.1% while lowering the thermal SEC from 7,800 kJ kg-1 to 6,517 kJ kg-1 in an unoptimized module. Finally, the scaling up potential of DHVMD process is briefly discussed.

brine↗
DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗
DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗
DOE OSTI2023

Automatic detection of cataclysmic variables from SDSS images

Abstract Investigating rare and new objects have always been an important direction in astronomy. Cataclysmic variables (CVs) are ideal and natural celestial bodies for studying the accretion process of semi-detached binaries with accretion processes. However, the sample size of CVs must increase because a lager gap exists between the observational and the theoretical expanding CVs. Astronomy has entered the big data era and can provide massive images containing CV candidates. CVs as a type of faint celestial objects, are highly challenging to be identified directly from images using automatic manners. Deep learning has rapidly developed in intelligent image processing and has been widely applied in some astronomical fields with excellent detection results. YOLOX, as the latest YOLO framework, is advantageous in detecting small and dark targets. This work proposes an improved YOLOX-based framework according to the characteristics of CVs and Sloan Digital Sky Survey (SDSS) photometric images to train and verify the model to realise CV detection. We use the Convolutional Block Attention Module to increase the number of output features with the feature extraction network and adjust the feature fusion network to obtain fused features. Accordingly, the loss function is modified. Experimental results demonstrate that the improved model produces satisfactory results, with average accuracy (mean average Precision at 0.5) of 92.0%, Precision of 92.9%, Recall of 94.3%, and $F1-score$ of 93.6% on the test set. The proposed method can efficiently achieve the identification of CVs in test samples and search for CV candidates in unlabeled images. The image data vastly outnumber the spectra in the SDSS-released data. With supplementary follow-up observations or spectra, the proposed model can help astronomers in seeking and detecting CVs in a new manner to ensure that a more extensive CV catalog can be built. The proposed model may also be applied to the detection of other kinds of celestial objects.

Astronomy & Astrophysics↗
DOE OSTI2020

Improved performance of stellarator coil design optimization

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Crystal structure and synchrotron X-ray powder reference pattern for the porous pillared cyanonickelate, Ni(3-amino-4,4′-bipyridine)[Ni(CN) 4 ]

The structure of Ni(3-amino-4,4′-bipyridine)[Ni(CN) 4 ] (or known as Ni-BpyNH 2 ) in powder form was determined using synchrotron X-ray diffraction and refined using the Rietveld refinement technique (R = 8.8%). The orthorhombic (Cmca) cell parameters were determined to be a = 14.7218(3) Å, b = 22.6615(3) Å, c = 12.3833(3) Å, V = 4131.29(9) Å 3 , and Z = 8. Ni-BpyNH 2 forms a 3-D network, with a 2-D Ni(CN) 4 net connecting to each other via the BpyNH 2 ligands. Further, there are two independent Ni sites on the net. The 2-D nets are connected to each other via the bonding of the pyridine “N” atom to Ni2. The Ni2 site is of six-fold coordination to N with relatively long Ni2–N distances (average of 2.118 Å) as compared to the four-fold coordinated Ni1–C distances (average of 1.850 Å). The Ni(CN) 4 net is arranged in a wave-like fashion. The functional group, –NH 2 , is disordered and was found to be in the m-position relative to the N atom of the pyridine ring. Instead of having a unique position, N has ¼ site occupancy in each of the four m-positions. The powder reference diffraction pattern for Ni-BpyNH 2 was prepared and submitted to the Powder Diffraction File (PDF) at the International Centre of Diffraction Data (ICDD).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advanced Method Optimization with Categorical and Constrained Continuous Parameters

Traditional approaches to analytical method optimization (e.g., univariate and “guess-and-check”) can be time-consuming, costly, and often fail to identify true optima within the parameter space. Previous work defined and implemented a generalized technique for method optimization for continuous method parameters, but a knowledge gap remains for the incorporation of categorical variables into these advanced method optimization schemes. This work presents and validates a generalized optimization approach that incorporates both continuous and categorical variables while also utilizing a multivariate, multiobjective optimization scheme with Karush–Kuhn–Tucker conditions to bound the optimization space to solutions within the physical limitations of the parameter space. Method optimization from a case study using GC–MS for the analysis of 11 analytical standards with objectives to minimize peak width and maximize peak height resulted in a 3 orders of magnitude improvement in the average peak height and a 2 orders of magnitude improvement in the average peak width compared to the least optimal (but reasonable) instrumental parameters utilized in this study. This approach to optimization allows for a customizable method optimization in which users can include both continuous and categorical variables to achieve objectives specific to their analytical goals. This approach significantly reduces the labor and cost associated with traditional method development approaches and can be applied in a variety of scientific fields across a range of laboratory techniques (e.g., instrument method development, sample preparation, and extraction techniques).

Amorphous materials↗
DOE OSTI2026

Position-Specific Carbon Isotope Fingerprinting of Fluorinated Organics and Degradation Products

Fluorinated organic compounds are of growing environmental and forensic relevance due to their widespread use in pharmaceuticals, agrochemicals, and consumer products, their environmental persistence, and potential ecological and human health impacts. Elucidating their sources and transformation pathways is therefore a major focus of current research. Stable carbon isotope analysis provides a powerful approach for tracing molecular origins and linking parent compounds to degradation products. Recent isotope measurements have largely relied on mass-spectrometry techniques, which provide only an average isotope ratio across a compound. In this work, we employ a novel nuclear magnetic resonance (NMR) spectroscopy tool to determine position-specific carbon isotope ratios ( 13 C/ 12 C) in organofluorine compounds and their degradation products. This approach enables isotope measurements without combustion or extensive purification and, crucially, resolves ratios at individual carbon positions rather than bulk averages. The resulting intramolecular isotope fingerprints are unique to a molecule’s source. Applied to selected pharmaceuticals and pesticides, these fingerprints allow discrimination of chemically identical compounds. Moreover, we show that the 13 C/ 12 C signature at the fluorinated carbon persists through degradation, demonstrated for lansoprazole and fipronil. The 19 F NMR data produced for the 13 C/ 12 C analyses are also well suited for impurity profiling, providing an additional dimension for fingerprinting fluorinated organics. These findings suggest that position-specific isotope analysis can serve as part of a broader suite of tools for source characterization of organofluorine compounds and their breakdown derivatives, with potential applications in product validation, forensics, and linking these compounds to their breakdown products.

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

Rapid and Simultaneous Characterization of Drug Conjugation in Heavy and Light Chains of a Monoclonal Antibody Revealed by High-Resolution Ion Mobility Separations in SLIM

Antibody-drug conjugates (ADCs) have recently gained traction in the biomedical community due to their promise for human therapeutics and an alternative to chemotherapy for cancer. Crucial metrics for ADC efficacy, safety, and selectivity are their drug-antibody ratios (DARs). However, DAR characterization (i.e., determining the average number of conjugated drugs on the antibody) through analytical methods remains challenging due to the heterogeneity of drug conjugation as well as the numerous post-translational modifications possible in the monoclonal antibody. In this work, we report on the use of high-resolution ion mobility spectrometry separations in our structures for lossless ion manipulations coupled to mass spectrometry (SLIM IMS-MS) platform for the rapid and simultaneous characterization of the drug load profile (i.e., stoichiometric distribution of the number of conjugated drugs present on the mAb), determination of the weighted average DAR in both the heavy and light chains of a model antibody-drug conjugate, and calculation of the overall DAR of the ADC. After chemical reduction of the ADC and a subsequent 31.5 meter SLIM IMS separation the various drug bound antibody species could be well resolved for both chains. Furthermore, we also show that higher resolution separations were possible with SLIM IMS as compared to ones performed on a commercially available (~1 meter) drift tube IMS-MS platform. Interestingly, additional gas-phase substructures become evident with longer path length SLIM IMS, potentially providing insight related to how drug conjugation may alter the structural properties of a monoclonal antibody and thus affect its overall potency and/or selectivity for a given antigen. We expect high-resolution SLIM IMS separations will augment the existing toolbox for ADC characterization, particularly for structural confirmation to enable the rapid optimization of DAR for a given ADC and thus better understand its potential toxicity and potency.

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

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo–Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

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

Local Distortions and Metal–Semiconductor–Metal Transition in Quasi-One-Dimensional Nanowire Compounds AV 3 Q 3 O δ (A = K, Rb, Cs and Q = Se, Te)

Metal cluster compounds have garnered renewed interest in the search for novel superconductors and topological semimetals owing to structural instabilities of metal-cluster geometries and broken symmetries. Here we synthesized needle-like crystals of the V-cluster-based quasi-one-dimensional (Q1D) materials AV 3 Q 3 O δ (A = K, Rb, Cs, Q= Se, Te) which can also be viewed as being composed of parallel nanowires. We examine how changes in their average and local structure control their electronic properties. All compounds crystallize in the TIFe 3 Te 3 -type structure (P6 3 /m space group) with infinite (V 3 )Q 3 ) - double-walled columnar chains separated by A + cations. Our single-crystal and synchrotron powder diffraction studies indicate oxygen atoms partially occupy the center site of the V 6 octahedral metal cluster cages in KV 3 Te 3 O 0.33 , RbV 3 Te 3 O 0.32 , and CsV 3 Te 3 O 0.35 , whereas KV 3 Se 3 is structurally oxygen-free. Our synchrotron X-ray pair distribution function (PDF) analyses indicate that the oxygen-free V 6 cluster octahedra in KV 3 Se 3 are highly distorted perpendicular to the chain direction even at room temperature, reducing the symmetry of the average structure from hexagonal P6 3 /m to monoclinic P2 1 /m. Our theoretical calculation supports this P2 1 /m distortion and suggests the structure further distorts to P2 1 or P2 1 /c at lower temperatures. In contrast, the oxygen-centered V-cluster in KV 3 Te 3 O 0.33 exhibits a V-3-triangle-trimerization along the chain direction. This feature is discernible from the local PDF and is consistent with lattice dynamical calculations based on density functional theory. Resistivity measurements indicate that KV 3 Se 3 exhibits metallic behavior, whereas a dramatic metal-semiconductor-metal transition emerges in KV 3 Te 3 O 6.33 , RbV 3 Te 3 O 0.32 , and CsV 3 Te 3 O 0.35 because of oxygen disorder and changes in local structure captured from our electronic structure analyses of the Fermi surface. Our investigation of the AV 3 Q 3 O 6 family demonstrates the importance of understanding local changes in structure driven by electronic instabilities, which can guide the search for new quantum materials in other low-dimensional cluster-compound materials.

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

Influence of Pelletization and Moisture Content of Oil Palm Empty Fruit Bunches (EFBs) on Dynamic Gasification Performance

Wet-ground empty fruit bunches (EFBs), dry-ground EFB, and pelletized EFB were gasified in a laboratory-scale updraft fixed-bed reactor at differing temperatures and equivalent ratios (ERs) between 0.19 and 0.78 to evaluate the performance and quality of the syngas under autothermal operating conditions. In addition, the theoretical estimation of the key performance indicator of the ashes on the performance of the gasification process was obtained. The results show that for ERs > 0.4 the autothermal operating condition was achieved for pretreatments with pelletized and dry-ground EFBs, with reactor operating temperatures T > 700 °C. With wet EFB, this condition is achieved at any ER value. The average composition of the syngas from pelletized EFB presented the best average concentrations of the combustible species CO, H2, and CH4 when the ER values were between 0.25 and 0.40. The highest hot and cold thermal efficiency values were obtained for the gasification of EFB pellets with ER values between 0.25 and 0.40. For ground EFB, increases in the ER value led to improvements in hot thermal efficiency. Finally, during the EFB pellet gasification tests, no agglomeration and sintering problems were evidenced in comparison with wet-ground and dry-ground EFBs that showed slag formation, especially at ER < 0.46.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

LaNi x Fe 1– x O 3–δ as a Robust Redox Catalyst for CO 2 Splitting and Methane Partial Oxidation

The current study reports LaNi 0.5 Fe 0.5 O 3–δ as a robust redox catalyst for CO 2 splitting and methane partial oxidation at relatively low temperatures (~700 °C) in the context of a hybrid redox process. Specifically, perovskite-structured LaNi x Fe 1–x O 3–δ (LNFs) with nine different compositions (x = 0.05–0.5) were prepared and investigated. Among the samples evaluated, LaNi 0.4 Fe 0.6 O 3–δ and LaNi 0.5 Fe 0.5 O 3–δ showed superior redox performance, with ~90% CO 2 and methane conversions and >90% syngas selectivity. The standalone LNFs also demonstrated performance comparable to that of LNF promoted by mixed conductive Ce 0.85 Gd 0.1 Cu 0.05 O 2–δ (CGCO). Long-term testing of LaNi 0.5 Fe 0.5 O 3–δ indicated that the redox catalyst gradually loses its activity over repeated redox cycles, amounting to approximately 0.02% activity loss each cycle, averaged over 500 cycles. This gradual deactivation was found to be reversible by deep oxidation with air. Further characterizations indicated that the loss of activity resulted from a slow accumulation of iron carbide (Fe 3 C and Fe 5 C 2 ) phases, which cannot be effectively removed during the CO 2 splitting step. Reoxidation with air removed the carbide phases, increased the availability of Fe for the redox reactions via solid-state reactions with La 2 O 3 , and decreased the average crystallite size of La 2 O 3 . As a result, reactivating the redox catalyst periodically, e.g., once every 40 cycles, was shown to be highly effective, as confirmed by operating the redox catalyst over 900 cumulative cycles while maintaining satisfactory redox performance.

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

Aqueous Bromide Discharges from U.S. Coal-Fired Power Plants: Points of Origin, Concentration Ranges, and Effluent Treatment Costs

Bromide discharges from coal-fired power plants have received increased attention from regulatory bodies due to their contribution to the formation of disinfection byproducts (DBPs) in downstream drinking water treatment plants. Here, this paper characterizes the relative contributions of bromide from coal feedstocks and bromine-based mercury control processes, estimates the distribution of bromide concentrations at 85 active coal-fired power plants across the United States (U.S.) with wet flue gas desulfurization units, and estimates the cost of bromide removal from wastewater discharge using year 2020 data. Bromide discharges are estimated at the plant level using a combination of the reported coal rank and composition combusted, estimates of bromide addition in mercury control techniques under multiple halogen addition scenarios, and the estimated flue gas desulfurization (FGD) wastewater flow rate. The median, simulated plant-level estimation of total FGD wastewater flow is 18.3 gallons/min at a bromide concentration of 319 mg/L, equivalent to ~11.6 tonnes/year of bromide discharges to the environment. Next, we evaluated the expected cost of employing the best available technology (BAT) to control bromide discharges in FGD wastewater to prevent contributions to DBP formation. Treatment would need to remove more than 99.8% of bromide to reach the 0.2 mg/L voluntary incentive program (VIP) limit. The total cost of treatment depends on whether disposal is on- or off-site; the average costs for all plants combined come to an average of $110 million ($95.2/kgal) in 2021 U.S. dollars for on-site disposal, or $134 million ($115/kgal) for off-site disposal.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Geochemical and Microstructural Characterization of Shale Rock: A Workflow for Site Evaluation for Energy–Fluids Storage

Shales are a central component of petroleum systems, as source, seal, and unconventional reservoir rocks. Unlike traditional unconventional shale, a newly emerging Caney Shale play is regarded as an unconventional of unconventional shales (UUS), due to its high content of fine-grained materials, lower reservoir porosity, dominance of nanopores, and a scarcity of visible natural fractures. Other developed unconventional shales such as the Woodford and Barnett have larger average pore size, higher porosity, and extensive natural fractures at core scale that contribute to reservoir quality. In this work, the Caney Shale was cored in entirety to characterize its interbedded ductile and reservoir intervals and establish criteria for their recognition. Representative ductile and reservoir intervals known as D2 and R3, respectively, were selected for detailed analysis including variations in elemental composition, mineralogy, facies, and the presence of natural fractures at well and core scales using X-ray fluorescence, X-ray diffraction, and thin section and core description. Characteristics of microstructure and microgeochemistry at nano- and microscales are compared using field emission scanning electron microscopy with energy dispersive spectroscopy and lowpressure nitrogen adsorption isotherms with fractal dimension analysis. The ductile (detrital clay-rich) member D2 is characterized by higher concentrations of Ti and Al and lower Si, while the reservoir R3 (possibly biogenic silica-accumulated) is characterized by lower Ti and Al and higher Si. Eight mixed carbonate–siliciclastic facies are recognized, and R3 shows a higher heterogeneity of facies stacking and average fracture abundance than D2. Further, R3 shows a more microscopically heterogeneous fabric/texture of matrix and a higher microporosity than D2 that has a higher pore surface heterogeneity. A fundamental understanding of the compositional and microstructural characteristics of UUS and further ductile/reservoir intervals will allow for a better assessment of reservoir quality, more effective production of hydrocarbons, and optimized selection of safe caprocks in carbon sequestration and subsurface hydrogen storage.

58 GEOSCIENCES↗
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗