Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Exploration”

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 163 records · Page 9

Applying Corrective Machine Learning in the E3SM Atmosphere Model in C++ (EAMxx)

The Simplified Cloud-Resolving E3SM Atmosphere Model (SCREAM) is the newest addition to the family of Earth System Models capable of explicitly resolving convective systems. SCREAM is a kilometer-scale configuration of the advanced E3SM Atmosphere Model (EAMxx), designed for heterogeneous systems. While the enhanced accuracy of kilometer-scale modeling offers significant benefits, it comes with a substantial computational cost, limiting feasible simulation durations to only a few years, even on the fastest supercomputers. Machine learning presents an opportunity for scientists to achieve the high accuracy of storm-resolving models at a significantly reduced cost. Building on the previous success of applying corrective machine learning (ML) to the FV3 model, this study explores the effects of implementing corrective ML in EAMxx-SCREAM. We also address the computational challenges of integrating the corrective ML, which is written in Python, with the C++/Kokkos EAMxx driver, as well as the potential pitfalls of generalizing an approach that was effective with one atmosphere model to another.

54 ENVIRONMENTAL SCIENCES↗

Robust Interpretation of Electrochemical Impedance Spectra Using Numerical Complex Analysis

Electrochemical Impedance Spectroscopy (EIS) is a noninvasive technique widely used for understanding charge transfer and charge transport processes in electrochemical systems and devices. Standard approaches for the interpretation of EIS data involve starting with a hypothetical circuit model for the physical processes in the device based on experience/intuition and then fitting the EIS data to this circuit model. This work explores a mathematical approach for extracting key characteristic features from EIS data by relying on fundamental principles of complex analysis. These characteristic features can suggest the presence of inductors and constant phase elements (nonideal capacitors) from impedance data and enable us to answer questions about the identifiability and nonuniqueness of equivalent circuit models. In certain scenarios such as models with only resistors and capacitors, we are able to enumerate all possible families of circuit models. Finally, we apply the mathematical framework presented here to real-world electrochemical systems and highlight results using impedance measurements from a lithium-ion battery coin cell.

charge transfer↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

STITCHES: a Python package to amalgamate existing Earth system model output into new scenario realizations

Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.

97 MATHEMATICS AND COMPUTING↗

Exploring Multidimensional Spatial-Temporal Hydropower Operational Flexibilities by Modeling and Optimizing Water-Constrained Cascading Hydroelectric Systems

Because of unique characteristics such as clean and cost-competitive electricity as well as fast-ramping and storage abilities, the power industry continues to evolve its operation strategies for cascading hydroelectric (CHE) systems for providing enhanced values to the grid, especially under the deeper renewable resource integration. However, existing operation practices of CHEs predate the integration of renewables, which could prohibit the effective utilization of their inherent flexibilities in delivering maximum financial benefits and providing valuable grid services to the power system and electricity market operations. Indeed, modeling and optimizing these resource-limited while flexible CHE assets with uncertainties and imperfect information across multiple spatial-temporal dimensions present significant challenges. To facilitate CHE facility operators in effectively coordinating water usage and hydropower plant operations across multiple timescales, this project aims to fill the existing gaps by developing a suite of accurate water inflow (WI) forecast models as well as enhanced CHE modeling and optimization approaches with proper consideration of their unique characteristics, which would help explore their multidimensional spatial-temporal operational flexibility potentials. The developed approaches could better align reservoir operation strategies with variability and uncertainty of future water availability. They can also promote more effective utilization of multidimensional spatial-temporal hydropower operational flexibility potentials by designing long-term evacuation plans of reservoirs and short-term operation of CHEs, along with their coordination with other types of renewables. The project leverages various resources to facilitate the research and development activities, including actual characteristics data of CHE systems and a library of current and future cases of Portland General Electric (PGE). These realistic data enable the project team to study how to maximize the value of CHEs under current and future portfolios and evaluate opportunities to improve operation practices.

13 HYDRO ENERGY↗

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

Modeling with uncertainty quantification reveals the essentials of a non-canonical algal carbon-concentrating mechanism

The thermoacidophilic red alga Cyanidioschyzon merolae survives its challenging environment likely in part by operating a carbon-concentrating mechanism (CCM). Here, we demonstrated that C. merolae 's cellular affinity for CO 2 is stronger than the affinity of its rubisco for CO 2 . This finding provided additional evidence that C. merolae operates a CCM while lacking the structures and functions characteristic of CCMs in other organisms. To test how such a CCM could function, we created a mathematical compartmental model of a simple CCM, distinct from those we have seen previously described in detail. The results of our modeling supported the feasibility of this proposed minimal and non-canonical CCM in C. merolae . To facilitate the robust modeling of this process, we measured and incorporated physiological and enzymatic parameters into the model. Additionally, we trained a surrogate machine-learning model to emulate the mechanistic model and characterized the effects of model parameters on key outputs. This parameter exploration enabled us to identify model features that influenced whether the model met the experimentally derived criteria for functional carbon concentration and efficient energy usage. Such parameters included cytosolic pH, bicarbonate pumping cost and kinetics, cell radius, carboxylation velocity, number of thylakoid membranes, and CO 2 membrane permeability. Our exploration thus suggested that a non-canonical CCM could exist in C. merolae and illuminated the essential features generally necessary for CCMs to function.

Steensma, Anne K. [Michigan State Univ., East Lans↗

Diverse fates of ancient horizontal gene transfers in extremophilic red algae

Abstract Horizontal genetic transfer (HGT) is a common phenomenon in eukaryotic genomes. However, the mechanisms by which HGT‐derived genes persist and integrate into other pathways remain unclear. This topic is of significant interest because, over time, the stressors that initially favoured the fixation of HGT may diminish or disappear. Despite this, the foreign genes may continue to exist if they become part of a broader stress response or other pathways. The conventional model suggests that the acquisition of HGT equates to adaptation. However, this model may evolve into more complex interactions between gene products, a concept we refer to as the ‘Integrated HGT Model’ (IHM). To explore this concept further, we studied specialized HGT‐derived genes that encode heavy metal detoxification functions. The recruitment of these genes into other pathways could provide clear examples of IHM. In our study, we exposed two anciently diverged species of polyextremophilic red algae from the Galdieria genus to arsenic and mercury stress in laboratory cultures. We then analysed the transcriptome data using differential and coexpression analysis. Our findings revealed that mercury detoxification follows a ‘one gene‐one function’ model, resulting in an indivisible response. In contrast, the ars H gene in the arsenite response pathway demonstrated a complex pattern of duplication, divergence and potential neofunctionalization, consistent with the IHM. Our research sheds light on the fate and integration of ancient HGTs, providing a novel perspective on the ecology of extremophiles.

59 BASIC BIOLOGICAL SCIENCES↗

Anomalous thermal transport in Eshelby twisted van der Waals nanowires

Dislocations in van der Waals (vdW) layered nanomaterials induce strain and structural changes that substantially impact thermal transport. Understanding these effects could enable the manipulation of dislocations for improved thermoelectric and optoelectronic applications, but experimental insights remain limited. Here, in this study, we use synthetic Eshelby twisted vdW GeS nanowires (NWs) with single screw dislocations as a model system to explore the interplay between dislocation-induced structural modifications and lattice thermal conductivity. Our measurements reveal a monoclinic structure stabilized by the dislocation, leading to a substantial drop in thermal conductivity for larger-diameter NWs (70% at room temperature), supported by first-principles calculations. Interestingly, we also find an anomalous enhancement of thermal conductivity with decreasing diameter in twisted NWs, contrary to typical trends in non-twisted GeS NWs. This is attributed to increased conductivity near the NW cores due to compressive strain around the central dislocations, and aligns with a density-functional-theory-informed core–shell model. Our results highlight the critical role of dislocations in thermal conduction, providing fundamental insights for defect and strain engineering in advanced thermal applications.

36 MATERIALS SCIENCE↗

Spatial analysis of cell patterning to aid genetic and phenotypic understanding of grass stomatal density: A case study in maize

Biological processes involve complex hierarchies where composite traits result from multiple component traits. However, holistically understanding of how sets of component traits interact to underpin genotype-to-phenotype relationships is generally lacking. Stomatal density (SD) is a tractable model system for exploring how high-throughput phenotyping (HTP) data could be exploited by a new spatial analysis approach to better understand a developmentally and functionally important trait. SD is a composite trait, resulting from various components related to cell identity and size, which are themselves governed by a series of spatio-developmental processes. Data from 192 recombinant inbred lines of maize [Zea mays (L.)] were analyzed by a new stomatal patterning phenotype (SPP) to (1) describe the average spatial probability distribution of the nearest neighboring stomata; (2) derive a core set of component traits related to cell size, cell packing, and positional probabilities; (3) build a structural equation model of component traits underlying SD; and (4) identify stomatal patterning quantitative trait loci (QTL). The core set of SPP-derived traits explained 74% of the variation in SD. Analyzing SPP component traits allowed some loci previously identified as generic SD QTL to be recognized as specific to lateral versus longitudinal elements of stomatal patterning. Therefore, this study highlights how novel insights can be gained by decomposing a composite trait (e.g., SD) into a set of component traits that were present in HTP data but not previously exploited.

59 BASIC BIOLOGICAL SCIENCES↗

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗

Applying Transfer Learning for Street-Scale Nuisance Flood Forecasting in Coastal-Urban Cities

An important challenge with Machine Learning (ML) is its transferability; that is, whether a ML model trained on one set of data can be applied to a second set of data without requiring a full re-training of the model. Transfer Learning (TL) addresses this challenge by transferring knowledge learned in the source domain (the data it was trained on) to the target domain (a second set of data that is statistically different but related, which the model was not trained on). This study investigates the use of TL for street-scale nuisance flood forecasting by exploring whether a ML model trained on data collected for one set of streets can effectively forecast flooding for another set of streets in the same city using TL. The envisioned use case is a city deploying a new flood depth monitoring sensor on a street and using TL to apply a ML model, trained on sensor data from an existing flood depth sensor network, to this new street. Eventually, the new flood depth sensor will have a sufficient dataset for training its own ML model, but TL can be used to fill the gap in time while this new dataset is being generated. This method is explored using a Long Short-Term Memory (LSTM) model trained on data for the flood-prone streets of Norfolk City, Virginia. The data used for training includes environmental time series (rainfall, tide), topographic features (Digital Elevation Model (DEM), Topographic Wetness Index (TWI), Depth To Water (DTW)), and street-scale flood depth time series obtained from a high-fidelity physics-based model, acting as a synthetic street-scale stream depth sensor dataset since actual stream depth sensor data is generally unavailable for most cities. A set of 180 flood-prone streets was used to train a base model, while another set of 180 flood-prone streets was used to re-train that model using different TL strategies. The results show that full-weight re-training proved most effective and minimal re-training of only the output layer was insufficient. The advantage of TL was most pronounced when target data was limited, meaning data collected at the new water depth sensor location included generally less than 18 flood events. As target data increased beyond 18 flood events, the benefit of TL diminished relative to training a ML model directly on the local flood events. These findings can assist cities as they implement street-scale flood sensing systems to create accurate forecasts for new sensing locations that do not yet have sufficient data records to train a local ML model.

Roy, Binata [Univ. of Virginia, Charlottesville, V↗

Optimization of the BDX experiment for Light Dark Matter searches at JLab

The Light Dark Matter (LDM) hypothesis postulates the existence of a new class of sub GeV particles, neutral under Standard Model (SM) interactions. In its simplest form, LDM consists of particles ¿ with masses below 1 GeV/c2, interacting with SM particles via a new force mediated by a light, spin-1 boson A0, commonly referred to as “Dark Photon”. This framework envisions a distinct “Dark Sector” with its own particles and interactions, offering a theoretically well motivated explanation for Dark Matter, consistent with astrophysical observations and a thermal production mechanism. The Beam Dump eXperiment (BDX) is an approved experiment at Jefferson Lab designed to search for Light Dark Matter. BDX will utilize an 11 GeV electron beam impinging on a thick target to produce a forward-boosted secondary beam of Light Dark Matter particles, which will then be detected by a dedicated downstream detector. Approved in 2018, BDX is expected to be commissioned in 2026 and run in 2027-2029. My thesis focuses on the preparatory work to deploy the BDX experiment. The detector design has been optimized to balance practicality with enhanced Light Dark Matter detection capabilities. Extensive characterization of detector components has been performed to ensure a precise understanding of detector response. A custom Monte Carlo framework has been developed to simulate Light Dark Matter signal and explore various theoretical models of interest. Additionally, a comprehensive data analysis framework has been developed to maximize the experiment sensitivity to Light Dark Matter, with optimizations based on expected detector performance. The ultimate goal of this thesis is to optimize BDX as a flagship experiment in Light Dark Matter searches, enabling it to probe different Dark Matter models.

Spreafico, M. [Univ. of Genova (Italy)]↗

Functional but not taxonomic diversity increases productivity of Populus in the southeastern United States

Plant interactions like competition and facilitation impact ecosystem function and resilience. Improving our understanding of the relationships between these interactions and community productivity has important implications for managers of production systems in forestry and agriculture as well as conservation science. Populus spp. are an excellent model system for exploring how inter‐ and intraspecific interactions impact ecosystem functions, such as productivity, in forest plantations. In this study, we compared aboveground productivity of six Populus clones from three different taxa grown in monoclonal and mixed‐clonal plots. The different mixture treatments were intended to experimentally test aboveground biomass response to contrasting levels of taxonomic diversity and functional diversity based on nitrogen use characteristics of Populus clones. We hypothesize that functional diversity would be more important than taxonomic diversity in increasing aboveground productivity of mixed‐clonal plantings compared to monocultures. In addition, a subset of treatments was carried out on additional sites representing a productivity gradient in order to determine if the relationship between biodiversity and productivity in these systems diminished at more productive sites as suggested by the stress‐gradient hypothesis. We found that functionally diverse mixtures of clones had greater yield of aboveground biomass than the average of their constituent monocultures, while more taxonomically diverse mixes of clones did not differ from the average of their constituent monocultures. However, when reestablished on sites with extremely high or low productivity, the best performing clone mixture also did not differ from the average of its constituent monocultures. Our results suggest that intimate clone mixtures of Populus have the potential to significantly increase productivity, but results vary by mixture and by site. To capitalize on positive biodiversity effects on yield in production systems, targeted mixtures based on divergent functional traits linked to different use and acquisition strategies for site‐specific limiting resources are most likely to be successful.

BEF↗

Influence of thermal treatment on structure and catalytic performance of ceria-zirconia supported copper oxide (CuO x /Ce y Zr 1-y O 2 ) catalysts for CO oxidation

Copper oxide (CuO x ) supported on ceria-zirconia (Ce y Zr 1-y O 2 , y = 1.0, 0.5, 0.0) catalysts were investigated to elucidate the effects of thermal treatment on their physicochemical properties and catalytic performance in carbon monoxide (CO) oxidation. Here, the catalysts were synthesized via a one-pot chemical vapor deposition (OP-CVD) method at 700˚C and 900˚C with controlled Cu loading. Characterization techniques, including synchrotron X-ray diffraction (S-XRD), Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), inductively coupled plasma spectroscopy (ICP), and N 2 adsorption-desorption, were implemented to probe the crystalline structure, molecular and electronic structure, oxygen vacancies, specific surface area (SSA) and metal loading. CO oxidation was chosen as a model reaction to explore the structure-catalytic performance relationship. A ∼100% CO conversion was achieved at < 150˚C, particularly with the CuO x /CeO 2 catalyst calcined at 700˚C. In contrast, calcination at 900˚C caused a ∼90% decrease in SSA and a ∼24% increase in T 50 . Activity tests revealed that increasing ZrO 2 content lowered CO oxidation activity despite generating more defect sites. In-situ measurement of the 700 °C calcined samples revealed the presence of stable and unstable defects in CuO x /Ce 0.5 Zr 0.5 O 2 and CeO 2 respectively, which play a key role in the activity of the catalysts. The results highlight that catalytic performance is closely related to the SSA. Furthermore, an optimum calcination temperature favor significant oxygen vacancy formation with required CuO x -support interactions, enhancing redox properties and catalytic performance.

36 MATERIALS SCIENCE↗

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗

Hydrogenation of calcite and change in chemical bonding at high pressure: Diamond formation above 100GPa

Synchrotron X-ray diffraction (XRD) and Raman spectroscopy in laser heated diamond anvil cells and first principles molecular dynamics (FPMD) calculations have been used to investigate the reactivity of calcite and molecular hydrogen (H 2 ) at high pressures up to 120 GPa. We find that hydrogen reacts with calcite starting below 0.5 GPa at room temperature forming chemical bonds with carbon and oxygen. This results in the unit cell volume expansion; the hydrogenation level is much higher for powdered samples. Single-crystal XRD measurements at 8–24 GPa reveal the presence of previously reported III, IIIb, and VI calcite phases; some crystallites show up to 4% expansion, which is consistent with the incorporation of ≤ 1 hydrogen atom per formula unit. At 40–102 GPa XRD patterns of hydrogenated calcite demonstrate broadened features consistent with the calcite VI structure with incorporated hydrogen atoms. Above 80 GPa, the C–O stretching mode of calcite splits suggesting a change in the coordination of C–O bonds. Laser heating at 110 GPa results in the formation of C–C bonds manifested in the crystallization of diamond recorded by in situ XRD at 300 K and 110 GPa and by Raman spectroscopy on recovered samples commenced with C 13 calcite. We explored several theoretical models, which show that incorporation of atomic hydrogen results in local distortions of CO 3 groups, formation of corner-shared C–O polyhedra, and chemical bonding of H to C and O, which leads to the lattice expansion and vibrational features consistent with the experiments. In conclusion, the experimental and theoretical results support recent reports on tetrahedral C coordination in high-pressure carbonate glasses and suggest a possible source of the origin of ultradeep diamonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shared mooring system designs and cost estimates for wave energy arrays

For floating renewable energy devices to become more cost-efficient and commercially scalable, their mooring system designs must be low-cost and suited for large-scale installations. Large arrays of floating devices, such as wave energy converters (WECs), will likely be designed with an individual mooring system for each device in the array. However, new mooring technology advancements provide options to use shared mooring lines to connect adjacent devices to one another, reducing the total number of anchors in the array, thereby reducing material use and cost. Here, this paper explores the design, modeling, and cost analysis of shared mooring systems for various sizes of arrays consisting of heaving oscillating water column (OWC) WECs. Shared mooring systems for WEC arrays sized in 2 x N and N x N grid layouts are designed to meet the relevant design standards, checking their performance with a nonlinear time-domain dynamic simulation, and the costs of each are calculated and compared. Several assumptions are taken in the design process to produce efficient results, providing a preliminary optimization for guidance on design decisions rather than a full, detailed design analysis. Mooring system costs per WEC were found to decrease as the number of WECs in the array increase, up to certain array sizes. The 2 x 3 array had the lowest mooring system cost per WEC out of all arrays considered, with a 60% cost reduction relative to using individual mooring systems. The 3 x 3 and 4 x 4 arrays achieved a 50% cost per WEC reduction. In addition to these significant cost reductions, the shared mooring system designs can provide advantages through smaller mooring system footprints, lower installation times, and less seabed disturbance.

16 TIDAL AND WAVE POWER↗