Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feature impact”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System: Preprint

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2023

The Interplay Between Ferroelectricity and Electrochemical Reactivity on the Surface of Binary Ferroelectric Al x B 1–x N

Polarization dynamics and domain structure evolution in ferroelectric Al 0.93 B 0.07 N are studied using piezoresponse force microscopy and spectroscopies in ambient and controlled atmosphere environments. The application of negative unipolar and bipolar first-order reverse curve (FORC) waveforms leads to a protrusion-like feature on the Al 0.93 B 0.07 N surface and a reduction of electromechanical response due to electrochemical reactivity. A surface change is also observed on the application of fast alternating current bias. At the same time, the application of positive biases does not lead to surface changes. Comparatively in a controlled glove box atmosphere, stable polarization patterns can be observed, with minuscule changes in surface morphology. This surface morphology change is not isolated to applying biases to free surface, a similar topographical change is also observed at the electrode edges when cycling a capacitor in an ambient environment. The study suggests that surface electrochemical reactivity may have a significant impact on the functionality of this material in the ambient environment. However, even in the controlled atmosphere, the participation of the surface ions in polarization switching phenomena and ionic compensation is possible.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anatomy of the electroweak phase transition for dark sector induced baryogenesis

We investigate the electroweak phase transition patterns for a recently proposed baryogenesis model with CP violation originated in the dark sector. The model includes a complex scalar singlet-Higgs boson portal, a U(1) l gauge lepton symmetry with a Z ' gauge boson portal and a fermionic dark matter particle. We find a novel thermal history of the scalar sector, featuring a Z 2 breaking singlet vacuum in the early Universe driven by a dark Yukawa coupling, that induces a one-step strongly first order electroweak phase transition. We explore the parameter space that generates the observed matter-antimatter asymmetry and dark matter relic abundance, while being consistent with constraints from electric dipole moment, collider searches, and dark matter direct detection bounds. The complex singlet can be produced via the Higgs portal and decays into Standard Model particles after traveling a certain distance. We explore the reach for long-lived singlet scalars at the 13 TeV Large Hadron Collider with \( \mathcal{L} \) = 139 fb - 1 and show its impact on the parameter space of the model. Setting aside currently unresolved theoretical uncertainties, we estimate the gravitational wave signatures detectable at future observatories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Lattice thermal transport in two-dimensional alloys and fractal heterostructures

Abstract Engineering thermal transport in two dimensional materials, alloys and heterostructures is critical for the design of next-generation flexible optoelectronic and energy harvesting devices. Direct experimental characterization of lattice thermal conductivity in these ultra-thin systems is challenging and the impact of dopant atoms and hetero-phase interfaces, introduced unintentionally during synthesis or as part of deliberate material design, on thermal transport properties is not understood. Here, we use non-equilibrium molecular dynamics simulations to calculate lattice thermal conductivity of $${\mathrm {(Mo|W)Se_2}}$$ ( Mo | W ) Se 2 monolayer crystals including $${\mathrm {Mo}}_{1-x}{\mathrm {W}}_x{\mathrm {Se_2}}$$ Mo 1 - x W x Se 2 alloys with substitutional point defects, periodic $${\mathrm {MoSe_2}|\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures with characteristic length scales and scale-free fractal $${\mathrm {MoSe_2}}|{\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures. Each of these features has a distinct effect on phonon propagation in the crystal, which can be used to design fractal and periodic alloy structures with highly tunable thermal conductivities. This control over lattice thermal conductivity will enable applications ranging from thermal barriers to thermoelectrics.

42 ENGINEERING↗
DOE OSTI2025

Meeting Global Health Needs via Infectious Disease Forecasting: Development of a Reliable Data-Driven Framework

Infectious diseases (IDs) have a significant detrimental impact on global health. Timely and accurate ID forecasting can result in more informed implementation of control measures and prevention policies. To meet the operational decision-making needs of real-world circumstances, we aimed to build a standardized, reliable, and trustworthy ID forecasting pipeline and visualization dashboard that is generalizable across a wide range of modeling techniques, IDs, and global locations. We forecasted 6 diverse, zoonotic diseases (brucellosis, campylobacteriosis, Middle East respiratory syndrome, Q fever, tick-borne encephalitis, and tularemia) across 4 continents and 8 countries. We included a wide range of statistical, machine learning, and deep learning models (n=9) and trained them on a multitude of features (average n=2326) within the One Health landscape, including demography, landscape, climate, and socioeconomic factors. The pipeline and dashboard were created in consideration of crucial operational metrics—prediction accuracy, computational efficiency, spatiotemporal generalizability, uncertainty quantification, and interpretability—which are essential to strategic data-driven decisions. While no single best model was suitable for all disease, region, and country combinations, our ensemble technique selects the best-performing model for each given scenario to achieve the closest prediction. For new or emerging diseases in a region, the ensemble model can predict how the disease may behave in the new region using a pretrained model from a similar region with a history of that disease. The data visualization dashboard provides a clean interface of important analytical metrics, such as ID temporal patterns, forecasts, prediction uncertainties, and model feature importance across all geographic locations and disease combinations. As the need for real-time, operational ID forecasting capabilities increases, this standardized and automated platform for data collection, analysis, and reporting is a major step forward in enabling evidence-based public health decisions and policies for the prevention and mitigation of future ID outbreaks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Efficiency, Power Loss, and Power Factor Measurement of Quadrupole Magnet Power Supplies at the Spallation Neutron Source

The linear accelerator (LINAC) quadrupole magnets are powered by 42 silicon-controlled rectifier (SCR) based power supplies at the Spallation Neutron Source (SNS) facility of Oak Ridge National Laboratory. These 35V, 525A power supplies are bulky, inefficient and require both air and water cooling. The reliability of the SNS facility is impacted due to water leaks internal to power supplies and current readback issues associated with their unique control system interface, resulting in multiple downtime events. Hence, an alternate solution is necessary for the continued reliable operation of the SNS. To mitigate the above-mentioned problems, this paper proposes the use of off-the-shelf Switch Mode Power Supplies (SMPS) rated for 20V, 500A with serial control interface. These SMPS are air-cooled, more efficient and more compact owing to their switching speeds of approximately 160 kHz. The performance enhancements of the SMPS in comparison with the existing SCR power supply are discussed in detail in this paper. The features of the SMPS, along with experimental results for both power supplies, like efficiency, power losses and stability, are presented. Ongoing work is also discussed.

Harave, Sudarshan↗
DOE OSTI2024

Predicting Li-ion Battery Performance for Impurity-doped NMC Cathodes Using Deep Learning

With the electric vehicle (EV) market expansion and the energy sector's shift towards electrification, the demand for battery metals, including lithium (Li), cobalt (Co), and nickel (Ni), is set to surpass supply. A critical knowledge gap exists in the purity standards for battery precursors and the impact of impurities on battery performance. Addressing this, our study employs a Deep Machine Learning (DL) based multi-objective optimization approach to interpret the relationship between metal impurities in domestic battery resources and their effects on battery performance. We analyze experimental data from Li-ion batteries with NMC (Nickel-Manganese-Cobalt oxide) cathodes over 1000 cycles, representing approximately ~6-8 months of operation, to establish a baseline of performance without impurities. Leveraging this data, we develop a Physics-Informed Deep Learning (PIDL) framework to extend our findings to cases that include metal impurities (e.g., Fe, Cu, Al) ranging from (0.001 - 0.01) %, respectively. By incorporating physics-based features, our PIDL model can accurately estimate the performance of NMC cathodes doped with various metal impurities to provide rapid design decisions. This research paves the way for informed decisions in Li-ion battery material design and optimization, ensuring the sustainable growth of the EV market and the broader energy sector.

25 ENERGY STORAGE↗
DOE OSTI2023

Probing transport energies and defect states in organic semiconductors using energy resolved electrochemical impedance spectroscopy

Abstract Determining the relative energies of transport states in organic semiconductors is critical to understanding the properties of electronic devices and in designing device stacks. Futhermore, defect states are also highly important and can greatly impact material properties and device performance. Recently, energy‐resolved electrochemical impedance spectroscopy (ER‐EIS) is developed to probe both the ionization energy (IE) and electron affinity (EA) as well as sub‐bandgap defect states in organic semiconductors. Herein, ER‐EIS is compared to cyclic voltammetry (CV) and photoemission spectroscopies for extracting IE and EA values, and to photothermal deflection spectroscopy (PDS) for probing defect states in both polymer and molecular organic semiconductors. The results show that ER‐EIS determined IE and EA are in better agreement with photoemission spectroscopy measurements as compared to CV for both polymer and molecular materials. Furthermore, the defect states detected by ER‐EIS agree with sub‐bandgap features detected by PDS. Surprisingly, ER‐EIS measurements of regiorandom and regioregular poly(3‐hexylthiophene) (P3HT) show clear defect bands that occur at significantly different energies. In regioregular P3HT the defect band is near the edge of the occupied states while it is near the edge of the unoccupied states in regiorandom P3HT.

14 SOLAR ENERGY↗
DOE OSTI2021

The impact of different parameterizations on the interpretation of CP violation in neutrino oscillations

CP violation in the lepton mass matrix will be probed with good precision in upcoming experiments. The amount of CP violation present in oscillations can be quantified in numerous ways and is typically parameterized by the complex phase δ PDG in the standard PDG definition of the lepton mixing matrix. There are additional parameterizations of the lepton mixing matrix as well. Through various examples, we explore how, given the current data, different parameterizations can lead to different conclusions when working with parameterization dependent variables, such as δ. We demonstrate how the smallness of |U e3 | governs the scale of these results. We then demonstrate how δ can be misleading and argue that the Jarlskog is the cleanest means of presenting the amount of CP violation in the lepton sector. We also confirm that, among the different parameterizations considered, the standard PDG parameterization has a number of convenient features.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

On-demand fabrication of piezoelectric sensors for in-space structural health monitoring

Abstract Inflatable structures, promising for future deep space exploration missions, are vulnerable to damage from micrometeoroid and orbital debris impacts. Polyvinylidene fluoride-trifluoroethylene (PVDF-trFE) is a flexible, biocompatible, and chemical-resistant material capable of detecting impact forces due to its piezoelectric properties. This study used a state-of-the-art material extrusion system that has been validated for in-space manufacturing, to facilitate fast-prototyping of consistent and uniform PVDF-trFE films. By systematically investigating ink synthesis, printer settings, and post-processing conditions, this research established a comprehensive understanding of the process-structure-property relationship of printed PVDF-trFE. Consequently, this study consistently achieved the printing of PVDF-trFE films with a thickness of around 40 µ m, accompanied by an impressive piezoelectric coefficient of up to 25 pC N −1 . Additionally, an all-printed dynamic force sensor, featuring a sensitivity of 1.18 V N −1 , was produced by mix printing commercial electrically-conductive silver inks with the customized PVDF-trFE inks. This pioneering on-demand fabrication technique for PVDF-trFE films empowers future astronauts to design and manufacture piezoelectric sensors while in space, thereby significantly enhancing the affordability and sustainability of deep space exploration missions.

Instruments & Instrumentation↗
DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗
DOE OSTI2023

Experimental Characterization of OpenMP Offloading Memory Operations and Unified Shared Memory Support

The OpenMP specification recently introduced support for unified shared memory, allowing implementation to leverage underlying system software to provide a simpler GPU offloading model where explicit mapping of variables is optional. Support for this feature is becoming more available in different OpenMP implementations on several hardware platforms. A deeper understanding of the different implementation’s execution profile and performance is crucial for applications as they consider the performance portability implications of adopting a unified memory offloading programming style. This work introduces a benchmark tool to characterize unified memory support in several OepnMP compilers and runtimes, with emphasis on identifying discrepancies between different OpenMP implementations as to how they various memory allocation strategies interact with unified shared memory. The benchmark tool is used to characterize OpenMP compilers on three leading High Performance Computing platforms supporting different CPU and device architectures. The benchmark tool is used to assess the impact of enabling unified shared memory on the performance of memory-bound code, highlighting implementation differences that should be accounted for when applications consider performance portability across platforms and compilers.

Elwasif, Wael↗
DOE OSTI2021

Adoption of image-driven machine learning for microstructure characterization and materials design: A Perspective

Microstructure characterization enables the development of structure-processing-property relationships critical to several research areas within the broad field of materials science, from alloy design to the assessment of corrosion resistance, and failure analysis. Conventional approaches to material characterization have relied on either qualitative inference by the human ex-pert or software applications that can extract high-level features from images, such as boundary segmentation, average grain diameter, etc. Such approaches rely heavily on subject matter expert user intervention and knowledge of what phases or more generally, what microstructural features, are of interest. The recent surge in the adoption of machine learning techniques to address problems in materials engineering has brought with it an increased interest and application of Image Driven Machine Learning (IDML) approaches. In this work, we review the applications of IDML to the field of materials characterization. A canonical hierarchy of stages is defined, which when put sequentially together completes an IDML study: problem definition, dataset building, model selection and training, model evaluation, and integration with existing instrumentation or simulation workflow. The studies reviewed in this work are analyzed from the perspective of each of these stages. Such a review permits agranular assessment of the field, for example the impact of IDML on materials characterization at the nanoscale, the size of a typical dataset required to train a semantic segmentation model on electron microscopy images, ubiquitousness of transfer learning in the domain, etc. Finally, we discuss the importance of interpretability and explainability in the field of IDML for materials characterization, and provide an overview of two emerging techniques in the field: semantic segmentation and generative adversarial networks.

Baskaran, Arun↗
DOE OSTI2020

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Characterizing Pore-Scale Geochemical Alterations in Eagle Ford and Barnett Shale from Exposure to Hydraulic Fracturing Fluid and CO 2 /H 2 O

As demand increases for an affordable energy source that is tied to an environmental obligation to reduce greenhouse gas emissions and water usage, there is a growing consideration in shale production utilizing processes such as 1) enhancing hydrocarbon recovery via carbon dioxide (CO 2 ) flooding, 2) using CO 2 as a fracturing agent to minimize water use, and 3) storing CO 2 in depleted shale formations to mitigate emissions to the atmosphere. Understanding the geochemical reactions and alterations that occur as shale is exposed to fluids and CO 2 is necessary to develop and optimize each of these processes for field applications. While the majority of shale formations are stimulated using traditional fracturing fluid, some may be fractured using CO 2 or other non-traditional means. We examine the effect fracturing fluid has on shale and how it behaves with secondary exposure to dry CO 2 or CO 2 -saturated water using in situ Fourier Transform infrared spectroscopy (FTIR), feature relocation scanning electron microscopy (SEM), and surface area and pore size analysis using volumetric gas sorption. These techniques were performed on Eagle Ford and Barnett shale samples that were exposed to fracturing fluid and unexposed (as received). Shales that have been exposed to traditional fracturing fluid experienced two reaction fronts. The first reaction front was formed during exposure to the fracturing fluid (pH of ~1.4). A secondary reaction front was formed as a result of CO 2 -saturated fluid exposure in the form of carbonic acid (pH ~5.6). These two different reaction mechanisms drove multiple dissolution and precipitation cycles which altered petrophysical properties of the shale and could lead to a significant impact on flow pathways. FTIR showed that equilibration of carbonate dissolution and precipitation cycles could take as long as 35 days. Samples exposed to fracturing fluid showed significantly less carbonate reactivity compared to those exposed to water. Pore size analysis results indicate exposure to fracturing fluid blocked small nanopores (10-15 nm) reducing BET surface area and total pore volume. SEM results show barite precipitated heavily during exposure to fracturing fluid. It appeared that carbonic acid was able to extract sulfur from organic matter to form gypsum evaporites. The mineralogical (barite precipitation and calcite dissolution/precipitation) and pore-scale alterations observed in these samples may lead to enhancement of flow pathways for injected CO 2 or produced hydrocarbons.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Surface parameters and bedrock properties covary across a mountainous watershed: Insights from machine learning and geophysics

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

58 GEOSCIENCES↗
DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗