Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “upscaling”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗
DOE OSTI2025

Discrete fracture network model benchmarks developed and applied in a DECOVALEX-2023 repository performance assessment study

This study presents newly developed benchmarks for modeling flow and transport within discrete fracture networks (DFNs) and useful methods for analyzing the results. The new benchmarks are designed to test modeling approaches for use in probabilistic performance assessment models of deep geologic repositories in fractured rock. The benchmarks simulate flow and transport through a 1 km 3 block of fractured rock. The first simulates migration of a short pulse of tracer through a simple network of four intersecting fractures. The second adds 1089 stochastically generated fractures. The third changes the pulse to a continuous point source. Evaluation of model performance relies on moment analysis and comparison of the results of different models. The expected nondimensional first moment of the conservative tracer for each benchmark is 1. The benchmarks were simulated by teams from Canada, Czechia, Germany, Korea, Sweden, Taiwan, and the United States as part of a DECOVALEX-2023 study (decovalex.org). The teams used various approaches, including explicit DFN modeling, DFN upscaling to an equivalent continuous porous medium (ECPM), and a combination of both methods. Transport mechanisms are modeled using either the advection-dispersion equation or particle tracking. Results demonstrate strong agreement among the models in breakthrough behavior up to the 75th percentile. Significant deviations in first moments and well-clustered outputs led to the identification of inaccuracies in several models. Such findings exemplify the benefit of exercising these benchmarks and using the presented methods to test DFN flow and transport models.

Benchmark↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
DOE OSTI2022

Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model

This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗
DOE OSTI2021

Scalable fabrication of high activity nanoporous copper powders for electrochemical CO 2 reduction via ball milling and dealloying

Electrochemical CO 2 reduction (ECR) is a promising technology to close the anthropic CO 2 circle using renewable energy to achieve carbon neutrality. Future commercialization of ECR will require the development of new catalyst synthesis routes that will allow significant upscaling of catalyst production from current research-level milligram quantities to the kilogram scale and beyond while maintaining the activity and selectivity demonstrated at the research level. Here, we report on generating and testing submicron-sized nanoporous copper (npCu) particles by using a scalable approach consisting of ball milling brittle Cu-based intermetallics followed by dealloying to add nanoporosity for high surface area. The resulting npCu particles have been tested in an industry-relevant large area (25 cm 2 ) electrolyzer platform and showed Faraday efficiencies (FE) for ethylene up to 34 % at current densities of 75–100 mA/cm 2 while keeping FE for hydrogen less than 30 %. Our results demonstrate that a combination of ball milling and dealloying is a promising approach to generate large quantities of high activity and high surface area npCu particles for ECR at an industry-relevant scale.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

A pilot-process for calcium hydroxide production from iron slag by low-temperature precipitation

The production of Ca(OH) 2 via thermal decomposition of limestone is an energy intensive process resulting in significant CO 2 emissions. However, if produced in a manner that obviates the need for the thermal decomposition of limestone, Ca(OH) 2 could be a “CO 2 -negative” material. Herein, we design and demonstrate the operation of a continuous, low-temperature (< 100 °C), aqueous-phase pilot-process to produce Ca(OH) 2 using calcium extracted from alkaline industrial wastes. The three-step process encompasses unit operations including: (i) calcium leaching from basic oxygen furnace BOF slag, (ii) leachate concentration by reverse osmosis (RO), and (iii) Ca(OH) 2 precipitation by temperature swing. Here this process presents several advantages compared to the traditional route. First, it advances the recycling of waste from the steel industry. Second, it operates significantly below the temperature used in conventional lime production (900 °C) as it bypasses limestone calcination, responsible for ~65% of the CO 2 emissions from the current process. Thus, the new process has a lower CO 2 footprint and enables the use of industrial waste heat. The mass and energy balances were quantified to reveal that increasing the calcium concentration of the feed solution and the precipitation temperature, decrease the energy demands of the RO step, thereby reducing the process’s overall CO 2 footprint. The pilot system operated continuously and achieved a production rate of nearly 1 kg per day of Ca(OH) 2 with a purity greater than 95 wt%. The average particle size of the precipitates depended on the residence time in the precipitation reactor, demonstrating an ability to produce size-controlled particulates. Importantly, the process achieved full water recirculation/reutilization level indicative of a low-consumable water demand. The outcomes offer new insights and understanding relevant to developing and upscaling low-CO 2 processes for cement, lime, and portlandite.

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

Enhanced joining strength in additive-manufactured polylactic-acid structures fused by embedded heated metallic meshes

Additively manufactured thermoplastic polymers, such as polylactic acid (PLA), hold significant promise for sustainable engineering structures, including wind turbine blades. Upscaling these structures beyond the limitations of 3D printer build volumes is a challenge; fusion joining presents a potential solution. This paper introduces a displacement-controlled resistance welding process for PLA, as an alternative to the typical force-controlled methods. Here, we investigated the bonding quality of resistance-welded and adhesive-bonded PLA beams through three-point bending and measured the surface deformations using digital image correlation. Different metal meshes (30%/0.11 mm Ni—Cu, 34%/0.07 mm Ni—Cu, and 36%/0.25 mm Co—Ni) served as heating elements. The process parameters were varied for the 34%/0.07 mm Ni—Cu mesh to identify an optimum set of parameters. Results showed that this optimized displacement-controlled welding achieved 94% of the original strength of monolithic samples. This indicates that the new welding process not only ensures high-quality bonding and fine surface finishing but also promotes sustainability, recyclability, and economic efficiency in various polymer and composite structural applications.

Additive manufacturing↗
DOE OSTI2022

OLEDs on planarized light outcoupling-enhancing structures in plastic

Light extraction from OLEDs remains a challenge. While consumer demand for OLEDs in display technology continues to grow as prices decline, commercial use of OLEDs in solid state lighting (SSL) applications is lagging due to low light extraction, which results in low efficiency. Previous studies have been directed toward increasing the light extraction factor η out , but cost effective approaches that will enable high throughput and potential upscaling, as well as elucidation of the extraction process, are still lagging. One promising approach is the use of buried light extraction enhancing patterns that are planarized with a high refractive index (RI) layer. When adding substrate mode extracting means, the highest external quantum efficiencies (EQEs) were achieved using a complex approach of planarized vacuum nano holes. Here we present η out -enhancing planarized extraction structures (PES) in plastic. While plastic substrates currently present long-term stability issues and other drawbacks, as we show, they provide a tool for simple, inexpensive, and rapid generation of various enhancing structures that importantly can be transferred to the preferred rigid and flexible glass substrates. Moreover, plastic substrates are attractive for biomedical applications. The successful structures we present are easily scalable periodic corrugations, e.g., with pitch a ~ 4.25 µm and pitch/depth ratio a/h ~2.4 planarized by a layer of RI ~1.9, resulting in maximal EQEs exceeding 60% for a green OLED and 48.5% for a white OLED. Comparable increased EQEs were achieved for structures with a ~ 15 and 16 µm, and a/h ~ 2.1 and 4, respectively. Enhanced extraction was also observed with a very simply fabricated quasiperiodic nanohole array, and surprisingly, with a random shallow design. The latter is at variance with simulations based on diffraction only. Results are discussed in terms of a, h, and structural parameters together with scattering matrix simulations to assess the effect of various PES parameters and the nanoparticle-embedded planarizing layer on η out . The experiments and simulations indicate contributions of diffraction, reflection, and scattering to EQE enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stochastic multiscale modeling for quantifying statistical and model errors with application to composite materials

This paper provides a coherent and efficient computational framework for stochastic multiscale analysis of material systems in the presence of parametric uncertainties and modeling errors. Uncertainty in those model parameters that are not deduced as upscaled quantities is attributed to an uncertainty “germ”. While such parameters can appear at any scale, they are predominant at the finest analysis scale. Additional uncertainties stemming from statistical estimation, attributed to lack of data and model error, are associated with each submodel contributing to the multiscale system. Here, a robust and efficient framework based on a generalized extended polynomial chaos expansion (gEPCE) is proposed to simultaneously propagate all these uncertainties in order to provide a probabilistic representation of specific quantities of interest (QoI). We characterize the full probability distribution of the QoI and the uncertainty in the failure probability pertaining to its tails. By combining gEPCE with kernel density estimation (KDE) and directional derivatives, we construct sensitivity measures that connect these statistical metrics of QoI to the various sources of uncertainty to assess their individual and combined impacts. An illustrative problem featuring three-point bending of a composite beam is investigated to demonstrate the presented approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving estimates of sub-daily gross primary production from solar-induced chlorophyll fluorescence by accounting for light distribution within canopy

Solar-induced chlorophyll fluorescence (SIF) has long been regarded as a proxy for photosynthesis and has shown superiority in estimating gross primary production (GPP) compared to traditional vegetation indices, especially in evergreen ecosystems. However, current SIF-based GPP estimations regard the canopy as a large leaf and seldom consider the impact of interactions among light, canopy structure, and leaf physiology. In this study, we proposed GPP estimation models with different descriptions of light–structure–physiology interactions (including the layered model, the two-leaf model, and the layered two-leaf model) and compared their performances with the big-leaf model using half-hourly (or hourly) observations at evergreen needleleaf forest sites. First, we found that the big-leaf model underestimated GPP, especially at noon. All models showed higher accuracy than that of the big-leaf model. Second, we investigated the diurnal dynamics of GPP estimations in each canopy layer and found that models with a two-leaf assumption captured the diurnal variations in GPP better than that with the layered assumption. Here we also deduced that the poor performance of the big-leaf model was related to its overestimation of the overall light stress on the redox state of PSII reaction centers (qL). Finally, we noticed that the qL at the canopy scale had lower sensitivity to light change than the single-leaf qL and that the light response of canopy-scale qL was influenced by the leaf area index during seasonal cycles. Overall, this study describes methods to accurately estimate sub-daily GPP from SIF in evergreen needleleaf forests and demonstrates that the interactions among light, canopy structure, and leaf physiology regulate the SIF-GPP relationship at the canopy scale. Further, it indicates the need to consider the description of light distribution within the canopy in next-generation terrestrial biosphere models, even if they incorporate SIF to constrain their parameterization. Thus, upscaling the established leaf-scale mechanistic SIF-GPP relationship or findings to canopy-scale applications still requires much work, especially when there are significant changes in environmental conditions and their within-canopy distributions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Observed impacts of large wind farms on grassland carbon cycling

Deployment of wind energy is an essential renewable energy source that mitigates climate change and reduces air pollution. Over the last several decades, wind energy development has increased worldwide, expanding from ~20 to ~900 GW (gigawatt) during 2001-2022. Nonetheless, researchers have identified unintended consequences of wind energy on microclimate via turbine-altered surface-atmosphere exchanges of energy, momentum, mass, and trace gases. Based on multi-source observations and models, researchers also have drawn some conclusions that wind farms could warm the land surface, especially at night, at regional and continental scales. Consequently, altered microclimates at wind farms may affect vegetation productivity and carbon sequestration, two critically important ecosystem services related to carbon dynamics; however, such potential impacts and driving mechanisms remain poorly understood. Wind energy deployment is increasing globally to meet carbon neutrality goals, with upscaling of onshore wind power capacity projected to grow from 542 GW in 2018 to 1787 and 5044 GW by 2030 and 2050, respectively. Furthermore, increased demand for wind energy deployment may lead to much larger wind farms in open, expansive landscapes. In turn, a large array of geographically clustered wind turbines could collectively modify local microclimate and amplify turbine-atmosphere interactions, which, if large enough, may produce detectable impacts on ecosystem dynamics. Thus, identifying and quantifying the potential impacts of wind farms on carbon-related ecosystem services may facilitate sustainable wind energy development globally.

17 WIND ENERGY↗
DOE OSTI2020

Modeling of thermal pressurization in tight claystone using sequential THM coupling: Benchmarking and validation against in-situ heating experiments in COx claystone

We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.

58 GEOSCIENCES↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2021

Quantum Computing with Exciton Qubits in Colloidal Semiconductor Nanocrystals

The present work evaluates the feasibility of quantum computing with exciton qubits in coupled colloidal semiconductor nanocrystals (NCs). A strategy for manipulating two-qubit states of colloidal NC hetero-dimers is described. We show that a sequence of laser pulses with the same photon energy can bring excitonic states of a nanocrystal hetero-dimer into entanglement and perform arbitrary qubit rotations (quantum gates). Our simulations of a realistic two-particle assembly of CdSe/CdS core/shell NCs demonstrate that such two-qubit gate operations can be driven by optical parametric oscillators with a theoretical error of 0.1%. Furthermore, a strategy for upscaling two-qubit hetero-dimers to N-qubit exciton gates in semiconductor NC assemblies is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗