Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust optimization”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Robust Heat-Flux Sensors for Coal-Fired Boiler Extreme Environments

In this project, robust heat-flux measurement systems were developed. The heat-flux sensors utilize thermoelectric effects to directly transduce the heat-flux inputs to analog electrical voltage signals. They were constructed from dedicated materials that can withstand temperatures of at least 1000°C and maintain adequate performance at these conditions for prolonged periods of time. The proposed approaches took into account numerous considerations, including system cost, sensor head resilience, sensor footprint, data accuracy, response time, and maintenance requirements. Through modern thermoelectric materials design, methodical materials selection and rigorous testing in materials characterization labs and medium-scale fire research facilities, we have demonstrated functioning laboratory prototypes, upon which one could base industrial heat-flux sensing platforms capable of operating in the challenging high-temperature, corrosive environments of the boilers of coal-fired power plants. A distributed sensor array for heat-flux measurements throughout the furnace water-wall, the superheater area and the economizer coils can provide critical data for the power plant control systems to increase efficiency, improve safety and reduce down times. For example, the combined heat-flux sensor/control systems can contribute to the optimization of burner and boiler operations under flexible loads, the optimization of heat-exchange conditions and overall reduction of heat rate and emissions, the prediction of imminent overheating conditions, and the optimization of the soot-blowing protocols.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Prediction of gas hydrate saturation using machine learning and optimal set of well-logs

We report resistivity and acoustic logs are widely used to estimate gas hydrate saturation in various sedimentary systems using one of the two popular methods ((1) acoustic velocity and (2) electrical resistivity), but the limitations of these two methods are often overlooked, which include (i) well-specific calibration of empirical exponents in the electrical resistivity method, (ii) assumption of known pore morphology for gas hydrates in the acoustic velocity method, and (iii) presence of unknown mineralogy and bulk modulus terms in the acoustic velocity method. NMR-density porosity-derived gas hydrate saturation based on the analysis of the transverse magnetization relaxation time (T2) is considered the most precise method, but acquisition of NMR-based logs is limited at relatively recent drilled sites; additionally, its use in conventional oil and gas reservoirs is not that common due to higher cost and operational deployment limitations associated with acquiring NMR well-logs. This study proposes a new method that predicts gas hydrate saturation (S h ) for any well using porosity, bulk density, and compressional wave (P wave) velocity well-logs with neural network (or stochastic gradient descent regression) without any well-specific calibration and/or other aforementioned shortcomings of the existing methods. The method is developed by examining the underlying dependency between S h and different combinations of well-logs, chosen from 6 routine logs, with 12 different machine learning (ML) algorithms. The accuracy of the proposed method in predicting S h is ~ 84%, which is better than the accuracy of seismic and electrical resistivity methods (≤ 75%) per the results reported by three different studies. The robustness of the method in the specific case of permafrost-associated gas hydrates is demonstrated with well-log data from two wells drilled on the Alaska North Slope.

58 GEOSCIENCES↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: Cosmology from galaxy×CMB lensing and shear×CMB lensing (Final Report for Grant SC0021429 "Cosomology with Galaxy and CMB Surveys: Mitigating Systematic Effects Through Combining Datasets")

We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langle δ_g$$^κ\text{CMB}\rangle$, and source galaxy weak lensing shear with CMB lensing convergence $\langle^{γκ}$$_{\text{CMB}}\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of $\langle δ_g$$^κ\text{CMB}\rangle + \langle^{γκ}$$_{\text{CMB}}\rangle$ changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from $\langle δ_g$$^κ\text{CMB}\rangle$ and $\langle^{γκ}$$_{\text{CMB}}\rangle$ are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in$\langle δ_g$$^κ\text{CMB}\rangle$ becomes subdominant to cosmic variance, preventing $\langle δ_g$$^κ\text{CMB}\rangle$ to further improve the constraints. This implies that to maximize the cosmological constraints from future $\langle δ_g$$^κ\text{CMB}\rangle$ + $\langle^{γκ}$$_{\text{CMB}}\rangle$ analyses, we should focus more on the requirements on$^{γκ}$$_{\text{CMB}}\rangle$ instead of$\langle δ_g$$^κ\text{CMB}\rangle$. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗
DOE OSTI2026

Extracting Critical Metals from Authentic Bauxite Residue

Securing domestic sources of critical minerals (CM) is paramount to ensuring a robust and self-sustaining technology infrastructure. Utilizing untapped non-conventional sources, like acid mine drainage, coal ash, bauxite residue/red mud, and more is an emerging approach that addresses this national concern while identifying new value-added pathways originating from waste streams. Red mud is a byproduct of the Bayer process that produces alumina and contains a wealth of CM – 50 µg/g Ga, 2.9 mg/g total rare earth elements plus yttrium, 12 mg/g Mn, 58 mg/g Al, and others. About 170 MM tonnes of red mud are co-produced annually alongside the 142 MM tonnes of alumina generated, highlighting the abundancy of this feedstock. In this work, different strategies were explored for extracting CM from authentic bauxite residue that involve thermal heating, microwave heating, and sonication all with different acids and buffers. CM recovery was then explored one step further by testing the adsorption of the extracted metals onto NETL’s Multi-functional Sorbent Technology (MUST). Microwave treatment of red mud proved the most effective CM extraction, whereas sonication was less desirable due to scale-up concerns. Successful recovery of CM from this leachate with MUST supports further studies toward optimizing the overall process.

bauxite residue↗
DOE OSTI2023

Ligand Control of Oxidation and Crystallographic Disorder in the Isolation of Hexavalent Uranium Mono-Oxo Complexes

Here, the development of high-valent transuranic chemistry requires robust methodologies to access and fully characterize reactive species. We have recently demonstrated that the reducing nature of imidophosphorane ligands supports the two-electron oxidation of U 4+ to U 6+ and established the use of this ligand to evaluate the inverse-trans-influence (ITI) in actinide metal-ligand multiple bond (MLMB) complexes. To extend this methodology and analysis to transuranic complexes, new small-scale synthetic strategies and lower-symmetry ligand derivatives are necessary to improve crys-tallinity and reduce crystallographic disorder. To this end, the synthesis of two new imidophosphorane ligands, [N=P t Bu(pip) 2 ]¯ (NPC 1 ) and [N=P t Bu(pyrr) 2 ]¯ (NPC 2 ) (pip = piperidinyl; pyrr = pyrrolidinyl) is presented, which break pseudo-C 3 axes in the tetravalent complexes, U[NPC 1 ] 4 and U[NPC 2 ] 4 . The reaction of these complexes with two-electron oxygen atom transfer reagents (N 2 O, trimethylamine N-oxide (TMAO) and 2,3:5,6-dibenzo-7-azabicyclo[2.2.1]hepta-2,5-diene (dbabhNO)) yields the U 6+ mono-oxo complexes U(O)[NPC 1 ] 4 and U(O)[NPC 2 ] 4 . This methodology is optimized for direct translation to transuranic elements. Of the two ligands, the NPC 2 framework is most suitable for facilitating detailed bonding analysis and assessment of the ITI. Theoretical evaluation of the U–(NPC) bonding confirms a substantial difference between axially and equatorially-bonded N atoms, revealing mark-edly more covalent U–N ax interactions. The U 6d+5f combined contribution for U–N ax is nearly double that of U–N eq , accounting for ITI shortening and increased bond order of the axial bond. Two distinct N-atom hybridizations in the pyrrolidine/piperidine rings are noted across the complexes, with approximate sp 2 and sp 3 configurations describing the slightly shorter P–N “planar” and slightly longer P–N “pyramidal” bonds, respectively. In all complexes, the NPC 2 ligands feature more planar N atoms than NPC 1 , in accordance with a higher electron-donating capacity of the former.

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

Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays

This is the final technical report for the Virginia Tech Center for Photonics Technology (VT-CPT) research project entitled “Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays”. In the the four-year effort, VT-CPT collaborated with the Departments of Mathmatics at VT and Sentek Instrument, LLC (Sentek) to develop a fiber-optic sensing system capable of real-time simultaneous distributed measurement of multiple subsurface, drilling, and production parameters. The ultra-sensitive fiber optic distributed acoustic sensing technology, picoDAS, developed by Sentek was integrated with a novel multi-material optical sensing fiber fabricated by CPT-VT to obtain distributed acoustic and magnetic field measurements with ultrahigh sensitivity and high spatial resolution. The sensing technology successfully developed and demonstrated under this research program is truly unique has application beyond subsurface imaging to include carbon storage and electrical grid monitoring, as well as healthcare and nuclear fusion. The project, sponsored by the Advanced Technology Program, began in 2019 and had the original goal of developing the next generation of harsh environment sensing systems. Exhaustive theoretical modeling was performed to optimize the multi-material optical fiber design. VT developed the processing techniques employed to successfully fabricate single mode fibers with metal cladding wires. Sensing fibers with nickel and Metglas cladding wires were fully characterized, inscribed with sensors, and fully integrated with the picoDAS system. Laboratory scale testing was performed to demonstrate performance upon exposure to lateral and transverse magnetic fields. Preliminary field trials were performed for prototype magnetic and acoustic sensing systems upon near surface deployment at VT. The distributed acoustic and magnetic field sensing system was advanced from a (Technology Readiness Level) TRL=2 to a TRL=5 via laboratory scale system validation in relevant environments. The truly one-of-its kind distributed magnetic and acoustic sensing system is expected to find immediate applications that require subsurface monitoring, to include carbon storage site monitoring, and has the potential to be a distributive technology in healthcare and electric grid markets. The technologies developed by Virginia Tech’s Center for Photonics Technology will support the mission of the National Energy Technology Laboratory (NETL) to drive innovation and deliver solutions for an environmentally sustainable and prosperous energy future by ensuring affordable, abundant, and reliable energy that drives a robust economy and national security. Technical accomplishments during the program are summarized briefly here and described in detail in the remainder of the report.

02 PETROLEUM↗
DOE OSTI2021

Crystallization and Phase Transformations of Aluminum (Oxy)hydroxide Polymorphs in Caustic Aqueous Solution

Gibbsite, bayerite, and boehmite are important aluminum (oxy)hydroxide minerals in nature and have been widely deployed in various industrial applications. They are also major components in caustic nuclear wastes stored at various U.S. locations. Knowledge of their crystallization and phase transformation processes contributes to understanding their occurrence and could help optimize waste treatment processes. While it has been reported that partial conversion of bayerite and gibbsite to boehmite occurs in basic solutions at elevated temperatures, systematic studies of factors affecting the phase transformation as well as the underlying reaction mechanisms are non-existent, particularly in highly alkaline solutions. We explored the effects of sodium hydroxide concentrations (0.1~3 M), reaction temperature (60~100 ?) and aluminum concentrations (0.1~1 M) on the crystallization and transformation of these aluminum (oxy)hydroxides. Detailed structural and morphological characterization by X-ray diffraction (XRD), scanning electron microscopy (SEM), and nuclear magnetic resonance (NMR) spectrometry revealed that these processes depend largely on the reaction temperature and the Al/OH- ratio. When 1 = Al/OH- = 2.5, the reactions favor formation of high crystallinity precipitates, whereas at Al/OH- ratio ? 2.5 precipitation ceases unless the Al concentration is higher than 1 M. We identified pseudoboehmite, bayerite and gibbsite as intermediate phases to bayerite, gibbsite and boehmite, respectively, all of which transform via dissolution-reprecipitation. Gibbsite transforms to boehmite in both acidic and weak caustic environments at temperatures above 80 oC. However, a ‘bar-shaped’ gibbsite morphology dominates in highly caustic environments (3 M NaOH). The findings enable a robust basis for selection of various solid phases by tuning the reaction conditions.

aluminum (oxy)hydroxide polymorphs, crystallizatio↗
DOE OSTI2023

Regulation of hair cell and stomatal size by a hair cell-specific peroxidase in the grass Brachypodium distachyon

The leaf epidermis is the outermost cell layer forming the interface between plants and the atmosphere that must both provide a robust barrier against (a)biotic stressors and facilitate carbon dioxide uptake and leaf transpiration. To achieve these opposing requirements, the plant epidermis developed a wide range of specialized cell types such as stomata and hair cells. While factors forming these individual cell types are known, it is poorly understood how their number and size is coordinated. Here, in this study, we identified a role for BdPRX76/BdPOX, a class III peroxidase, in regulating hair cell and stomatal size in the model grass Brachypodium distachyon. In bdpox mutants prickle hair cells were smaller and stomata were longer. Because stomatal density remained unchanged, the negative correlation between stomatal size and density was disrupted in bdpox and resulted in higher stomatal conductance and lower intrinsic water-use efficiency. BdPOX was exclusively expressed in hair cells suggesting that BdPOX cell autonomously promotes hair cell size and indirectly restricts stomatal length. Cell wall autofluorescence and lignin stainings indicated a role for BdPOX in lignification or crosslinking of related phenolic compounds at the hair cell base. Ectopic expression of BdPOX in the stomatal lineage increased phenolic autofluorescence in guard cell walls and restricted stomatal elongation in bdpox. Together, we highlight a developmental interplay between hair cells and stomata that optimizes epidermal functionality. We propose that cell-type-specific changes disrupt this interplay and lead to compensatory developmental defects in other epidermal cell types.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

First-principles elucidation of defect-mediated Li transport in hexagonal boron nitride

Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B 3 N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.

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

Dark Energy Survey Year 3 results: A 2.7% measurement of baryon acoustic oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of baryon acoustic oscillations (BAO) in an optimized sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg 2 with 0.6 < z photo < 1.1 and a typical redshift uncertainty of 0.03 ( 1 + z ) . The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift, and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive preunblinding protocol. We used two different statistics; angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance D M at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain D M ( z eff = 0.835 ) / r d = 18.92 ± 0.51 , which is consistent with, but smaller than, the Planck prediction assuming flat Λ CDM , at the level of 2.3 σ . The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at z > 0.75 . When combined with DES 3 x 2 pt + SNIa , they lead to improvements in H 0 and Ω m constraints by ∼ 20 % .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Dark Energy Survey: A 2.1% measurement of the angular baryonic acoustic oscillation scale at redshift z eff = 0.85 from the final dataset

Here, we present the angular diameter distance measurement obtained with the baryonic acoustic oscillation (BAO) feature from galaxy clustering in the completed Dark Energy Survey, consisting of six years (Y6) of observations. We use the Y6 BAO galaxy sample, optimized for BAO science in the redshift range 0.6 < z <1.2, with an effective redshift at z eff = 0.85 and split into six tomographic bins. The sample has nearly 16 million galaxies over 4,273 square degrees. Our consensus measurement constrains the ratio of the angular distance to sound horizon scale to D M ⁡(z eff )/r d = 19.51 ± 0.41 (at 68.3% confidence interval), resulting from comparing the BAO position in our data to that predicted by planck Λ⁢CDM via the BAO shift parameter α =(D M /r d )/(D M /r d ) PLANCK . To achieve this, the BAO shift is measured with three different methods, angular correlation function (ACF), angular power spectrum (APS), and projected correlation function (PCF), obtaining α = 0.952 ± 0.023, 0.962 ± 0.022, and 0.955 ± 0.020, respectively, which we combine to α = 0.957 ± 0.020, including systematic errors. When compared with the Λ⁢CDM model that best fits planck data, this measurement is found to be 4.3% and 2.1⁢σ below the angular BAO scale predicted. To date, it represents the most precise angular BAO measurement at z > 0.75 from any survey and the most precise measurement at any redshift from photometric surveys. The analysis was performed blinded to the BAO position, and it is shown to be robust against analysis choices, data removal, redshift calibrations, and observational systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

FiXR : a framework to reconstruct fiber cross-sections from X-ray fiber diffraction experiments

Ab initio reconstruction methods have revolutionized the capabilities of small-angle X-ray scattering (SAXS), allowing the data-driven discovery of previously unknown molecular conformations, exploiting optimization heuristics and assumptions behind the composition of globular molecules. While these methods have been successful for the analysis of small particles, their impact on fibrillar assemblies has been more limited. The micrometre-range size of these assemblies and the complex interaction of their periodicities in their scattering profiles indicate that the discovery of fibril structures from SAXS measurements requires novel approaches beyond extending existing tools for molecular discovery. In this work, it is proposed to use SAXS measurements, together with diffraction theory, to infer the electron distribution of the average cross-section of a fiber. This cross-section is modeled as a discrete electron density with continuous support, allowing representations beyond binary distributions. Additional constraints, such as non-negativity or smoothness/connectedness, can also be added to the framework. The proposed approach is tested using simulated SAXS data from amyloid β fibril models and using measured data of Tobacco mosaic virus from SAXS experiments, recovering the geometry and density of the cross-sections in all cases. The approach is further tested by analyzing SAXS data from different amyloid β fibril assemblies, with results that are in agreement with previously proposed models from cryo-EM measurements. The limitations of the proposed method, together with an analysis of the robustness of the method and the combination with different experimental sources, are also discussed.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Holocene paleoclimate change in the western US: The importance of chronology in discerning patterns and drivers

We report sediment in lakes and meadows forms a powerful archive that can be used to reconstruct environmental change through time. Reconstructions of lake level, of chemical, biological, and hydrological conditions, and of surrounding vegetation provide detailed information about past climate conditions, both locally and regionally. Indeed, most of our current knowledge of centennial- to millennial-scale climate variability in the arid western United States, where information about past hydroclimate is particularly important, comes from such sediment-based reconstructions. The pressing need for robust, precise predictions of future conditions is a significant motivation for paleoclimate science, and current research questions frequently require Holocene reconstructions to be resolved at sub-centennial timescales. Increasingly, regional syntheses seek to identify synoptic-scale patterns similar to those defined from modern observations (seasonal, interannual, multi-decadal, etc.) or to compare with the output of climate model simulations. However, the age control on existing records, especially those more than about 20 years old, is often sufficient only for millennial-scale interpretation. Here we assess the age control for 84 published and unpublished records from lakes and meadows in the Great Basin, California, and desert southwest, and use Bayesian modeling to evaluate the 95% uncertainty ranges for the 42 best-dated records. In the Late Holocene, about half of the 42 records have <400-year mean uncertainty ranges; however, high-precision age control is especially critical for young records, used to develop an accurate understanding of a proxy’s response to known climate variations. In the Middle Holocene, records vary from 400 to >800-year mean uncertainty and records of the Early Holocene have 600- to >1400-year mean uncertainty ranges. We find that the largest control on modeled uncertainties is dating density, with at least 2 dates/kyr being optimal and suggest obtaining “range-finder” dates at the onset of a study to better predict the total number of dates needed for an adequate age model. Such a density avoids a commonly observed phenomenon of significant peaks in uncertainty arising in gaps between age control points. Analysis of the uncertainties associated with proxy shifts reveal that more than half are >400 years. Although such large uncertainties currently prevent sub-centennial interpretations in most cases, increased dating density, strategic use of limited funds (including budgeting for a 2 date/kyr minimum at the proposal stage), construction of age-depth models with Bayesian methods, and critical evaluation of chronological uncertainty will shed light on past climate variability at finer timescales, enhancing our understanding of global and regional drivers of western U.S. climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗
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↗
DOE OSTI2022

Scaling up carboxylic acid production from cheese whey and brewery wastewater via methane-arrested anaerobic digestion

In a circular economy, organic waste streams are valuable resources for sustainably producing chemicals and fuels. This study investigates a new methane-arrested anaerobic digestion (MAAD) process that converts high-strength cheese whey and brewery wastewater into carboxylic acids. The process was developed and optimized under various bench-scale semi-continuous (fed-batch) operating conditions (e.g., retention time, organic loading rate, pH, and feed/harvest frequency). The MAAD responses to various control-failure scenarios were also systematically investigated. The highest total acid productivity was 26g/(L liq ·d) with a substrate conversion of 0.79g COD digested /g COD fed at a hydraulic retention time (HRT) of approximately 2 d in a 14-L digester. The most stable conditions for digester operation (HRT 3 d at pH 6.0 and 40°C) were selected for process scale-up to 100gal (~380 L). Semi-continuous, pilot-scale MAAD successfully produced a total acid concentration of 40.6+/-1.1 g/L with 8.1% acetic acid, 45.1% butyric acid, and 44.5% lactic acid. The links between wastewater characteristics, operation mode, digester scale, and microbial community structure were statistically analyzed. The results show genera Sporolactobacillus and Clostridium positively correlate with butyric acid production (Pearson correlation coefficient>0.5). Moreover, four kinetic models were developed and fit to batch MAAD experimental datasets (R 2 >95%) and were successfully applied to predict the total acid production in both bench-and pilot-scale semi-continuous MAADs. In conclusion, this study shows MAAD has the potential for industrial-scale applications and is a robust platform to valorize low- or negative-value waste streams into high-value bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗