SEARCH · Search NASA
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.
One of the tasks completed under MRCSP was the analysis of CO2 flow patterns using DTS data. This task is based on DTS data acquired from the Chester 16 pinnacle reef located in Otsego County, Michigan.
Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.
This report describes the implementation and testing of domain-overlapping coupling of the Pronghorn and SAM codes. Both Pronghorn and SAM are codes developed by the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for analysis of advanced nuclear reactors. Pronghorn focuses on analyzing multi-dimensional core flow conditions, while SAM focuses on analysis of the entire system including, among other components, piping, pumps, and heat exchangers. It is desirable for many advanced reactor thermal-hydraulics simulations to couple Pronghorn and SAM to obtain self-consistent solutions in the core and the system. The domain-overlapping coupling approach provides a robust and numerically efficient candidate for this coupling. We describe the algorithm for obtaining self-consistent solutions between a SAM network of flow channels and an overlapped multi-dimensional Pronghorn domain. The algorithm is tested for a variety of simple canonical test problems and a model of the Molten Salt Reactor Experiment (MSRE). Results demonstrate that the domain overlapping coupling approach provides consistent pressure, mass flow rate, enthalpy, and passive scalar distributions between Pronghorn and SAM and converges reliably and efficiently (at most 15 iterations, usually much less than 10) for the considered test problems.
In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.
The deposition of thick B 4 C films with low residual stress by conventional direct-current magnetron sputtering is accompanied by the formation of dust particulates contaminating the target, chamber, and substrates and leading to the formation of nodular defects in films. Here, we demonstrate that the formation of particulates is greatly reduced during radio-frequency magnetron sputtering (RFMS). We systematically study properties of B 4 C films deposited by RFMS with a substrate temperature of 330 °C, a target-to-substrate distance of 10 cm, Ar working gas pressure in the range of 4.5–12.0 mTorr (0.6–1.6 Pa), and substrate tilt angles of 0°–80°. All films are x-ray amorphous. A columnar structure develops with increasing either Ar pressure or substrate tilt. For columnar films, the column tilt angle decreases with increasing Ar pressure, which we attribute to a corresponding increase in the width of the distribution of impact angles of deposition flux. In contrast to the Keller–Simmons rule, the deposition rate increases with increasing Ar pressure, which suggests a better coupling of the RF energy to the plasma processes that lead to target sputtering at higher pressures. There is a critical substrate tilt angle above which the total residual stress is close to zero. This critical substrate tilt angle is ~0° for an Ar pressure of 12 mTorr (1.6 Pa). The lower residual stress state, necessary for depositing ultrathick films, is characterized by a larger concentration of nanoscale inhomogeneities and decreased mechanical properties. Based on these results, RFMS deposition of 60-μm-thick B 4 C films is demonstrated.
Commercial scale decarbonization through carbon capture and storage may likely involve many CO 2 storage projects located in close proximity. The close proximity could raise concerns over caprock integrity associated with reservoir pressure buildup and interference among adjacent projects. Commercial-scale injection will also require large prospective CO 2 storage resource and high injectivity in the targeted storage formations. To accommodate the need for both large resource and high injectivity, project operators could consider injecting CO 2 into a stacked sequence of formations. This analysis investigates the benefits of injecting CO 2 into a vertically stacked sequence of saline formations, over injecting the same amount of CO 2 into a single saline formation, in addressing these challenges. Our analysis shows that injecting into the stacked sequence mitigates the extent of pressure buildup among the stacked formations, while still achieving the same or greater target CO 2 storage volumes. Among cases modeled, the resulting pressure buildup front is most reduced when each storage site distributes injection volumes over several wells, each of which injects a portion of the total CO 2 mass across the stacked sequence. This favorable case not only results in the smallest CO 2 aerial footprint, but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 46% compared to the single-formation storage), the result of which is crucial to maintain caprock integrity. This analysis provides insights into required decision-making when considering multi-project deployment in a shared basin.
Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.
Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.
The Integrated Waste Treatment Unit (IWTU) has been tasked with using steam reforming technology to treat 3.4 million liters (900,000 gallons) of radioactive sodium-bearing waste (SBW) at the Idaho National Laboratory (INL). The project has been hampered by a number of processing issues associated with process equipment during simulant testing, with the two main issues involving the fluidizing gas distribution system in the main reactor vessel, the denitration mineralization reformer (DMR), and high differential pressures in the main filter vessel, the process gas filter (PGF). Following the replacement of the fluidizing gas distribution system in the DMR, a successful 30-day demonstration run was performed that revealed unacceptable bridging in the filter bundles of the PGF. The bundles were replaced with new ones, and a 50-day demonstration run was performed. The main purpose of this run was to operate the IWTU over a range of operating conditions. Overall, the run was a success; however, there were multiple challenges that arose from the run. (authors)
Mercury injection capillary pressure (MICP) analysis was used to study organic-matter-rich Devonian shales varying in maturity from marginally mature to postmature to discuss controls on drainage and imbibition in shale reservoirs. Six samples come from the New Albany Shale; their total organic carbon (TOC) ranged from 3.3–8.4 % and vitrinite reflectance (R o ) from 0.54–1.42 %. One Marcellus Shale sample with R o 2.50 % was used to extend the maturity range to the advanced dry gas window. Three coal samples of varying rank were also included for comparison. Porosity of the shales ranges from 1.5–5.0 %, and MICP-derived permeability ranges from 0.17–9.9 mD. Injection capillary pressure curves have a range of entry pressures from 0.055 to0.421 MPa, and no relationship of entry pressure with maturity was revealed. Maturity, however, has an influence on pore size distribution, with the higher-maturity samples having significantly lower volumes of pores with throats less than 10 nm in diameter. Withdrawal capillary pressure curves document a wide range of withdrawal efficiencies from 55 % to almost 100 %, and pores having throats between 3–50 nm in diameter play a major role in the ability of mercury to be released from the pore system. The obtained data provide valuable information about shale properties and their behavior not only in the air-mercury system but also under reservoir conditions.
In this work, we report on defect generation in multilayer GaSe through hydrostatic pressure quenching and UV laser irradiation. The Raman line width from the UV 266 nm irradiated sample is much wider than that in pressure-quenched GaSe, corresponding to a wider defect energy distribution range in the former sample than the latter. After quenching from 11.2 GPa, three photoluminescence (PL) peaks from defect states are observed at 657, 681, and 695 nm at a low temperature of 93 K. Defect-related peaks at 649, 694, 750, and 774 nm also appear in low-temperature PL spectra after UV laser irradiation, with a nonmonotonous intensity dependence on irradiation duration. There are common features in defects produced by these two methods: the PL peaks with the lowest energy are sharp, and their PL intensities increase linearly with the excitation laser power and saturate above a certain excitation laser power. These two features are similar to those in defects for single-photon emission (SPE) in other 2D materials at even lower temperatures. Fluorescence lifetime imaging shows distinguished short (2.3 ns) and long (75.6 nm) lifetimes of the 695 nm PL line in pressure-quenched GaSe. The density functional theory predicts defect energy levels related to Se vacancy.
Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.
Design of an additively manufactured molten salt (MS) to supercritical carbon dioxide (sCO 2 ) primary heat exchanger (PHE) for solar thermal power generation is presented. The PHE is designed to handle temperatures up to 720 °C on the MS side and an internal pressure of 200 bar on the sCO 2 side. In the core, MS flows through a three-dimensional periodic lattice network, while sCO 2 flows within pin arrays. The design includes integrated sCO 2 headers located within the MS flow, allowing for a counterflow design of the PHE. The sCO 2 headers are configured to enable uniform flow distribution into each sCO 2 plate while withstanding an internal pressure of 200 bar and minimizing obstruction to the flow of MS around it. The structural integrity of the design is verified on additively manufactured (AM) 316 stainless steel sub-scale specimens. An experimentally validated, correlation-based sectional PHE core thermofluidic model is developed to study the impact of flow and geometrical parameters on the PHE performance, with varied parameters including the mass flowrate, surface roughness, and PHE dimensions. A process-based cost model is used to determine the impact of parameter variation on build cost. The model results show that a heat exchanger with a power density of 18.6 MW/m 3 (including sCO 2 header volume) and effectiveness of 0.88 can be achieved at a heat capacity rate ratio of 0.8. As a result, the impact of design and AM machine parameters on the cost of the PHE are assessed.
Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.
An X-pinch load driven by an intense current pulse (>100 kA in ~100 ns) can result in the formation of a small radius, runaway compressional micro-pinch. A micro-pinch is characterized by a hot (>1 keV), current-driven (>100 kA), high-density plasma column (near solid density) with a small neck diameter (1–10 µm), a short axial extent (<1 mm), and a short duration (≲1 ns). With material pressures often well into the multi-Mbar regime, a micro-pinch plasma often radiates an intense, sub-ns burst of sub-keV to multi-keV x rays. A low-density coronal plasma immediately surrounding the dense plasma neck could potentially shunt current away from the neck and thus reduce the magnetic drive pressure applied to the neck. Here, to study the current distribution in the coronal plasma, a Faraday rotation imaging diagnostic (1064 nm) capable of producing simultaneous high-magnification polarimetric and interferometric images has been developed for the MAIZE facility at the University of Michigan. Designed with a variable magnification (1–10×), this diagnostic achieves a spatial resolution of ~35 µm, which is useful for resolving the ~100-μm-scale coronal plasma immediately surrounding the dense core. This system has now been used on a reduced-output MAIZE (100–200 kA, 150 ns) to assess the radial distribution of drive current immediately surrounding the dense micro-pinch neck. The total current enclosed was found to increase as a function of radius, r, from a value of ≈ 50 ± 25 kA at r ≈ 140 µm (at the edge of the dense neck) to a maximal value of ≈150 ± 75 kA for r ≥ 225 µm. This corresponds to a peak magnetic drive pressure of ≈75 ± 50 kbar at r ≈ 225 µm. The limitations of these measurements are discussed in the paper.
A flow cell is a necessary measurement interface for some important optical analyzers. In our application of the flow cell, we utilize a state-of-the-art optical analyzer (Malvern Insitec) to measure particle size distribution and concentration in-situ in a sampled flow from a novel pressurized oxygen-fuel combustion process. However, since this sampling flow is a flow of moist flue gas under high temperature and pressure, and the flow contains particles and corrosive acid gases, it is an extreme challenge to obtain a flow cell with a high optical quality that does not perturb the measurement. To address this challenge, we propose a new design for an optical flow cell. By using a unique flow field in the proposed flow cell, the measurement zone can be well defined by the sampling flow, minimizing the influence of purge flow. To demonstrate this flow cell, we have built a test system, and conduct measurements utilizing polydisperse-particle standards (10-100 µm and 1-10 µm). The results reveal that the optical windows are well protected by the purge flow field, without risk of deposition from the sampling flow, and the Malvern Insitec can measure the particle size distribution by using this flow cell, without generating sample bias.
Vapor co-deposited metal nanocomposites with immiscible components have shown to have several superior properties including mechanical performance, radiation tolerance, and enhanced functionalities. The morphology and structural design of these nanocomposites are important in determining the expression of these properties. The conditions used to deposit thin film nanocomposites are the governing factors in determining the final nanostructure. In this work, we link the effects of process gas pressure, substrate temperature, and sputtering target power on resultant Cu/Nb thin film nanostructures. Here, we find that for films with equiatomic phase fractions of Cu and Nb, inhomogeneous precipitate segregation dominated the structure at lower substrate temperature depositions while higher temperatures distributed the Cu and Nb phases more uniformly. Increasing process gas pressure during deposition lead to a spatially homogeneous, nanocrystalline structure of Cu/Nb. In situ annealing of this sample shows that the phase separated morphology of the two immiscible materials depends on the initial composition. Finally, compositions of Cu and Nb that were highly disparate lead to interface segregation of the minority phase. The morphological evolution observed in the Cu/Nb system is compatible with modified structural zone models for immiscible thin film systems.