Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relative permeability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

Tensile Properties of Hydrogels and of Snake Skin

Stimulus-responsive or 'smart' gels are of potential interest as sensors and actuators, in industrial separations, and as permeable delivery systems. In most applications, a certain degree of mechanical strength and toughness will be required, yet the large-strain behavior of gels has not been widely reported. Some exceptions include work on gelatin and other food gels, some characterization of soft gels applicable for in-vitro cell growth studies, and toughness determinations on commercial contact lens materials. In general, it can be anticipated that the gel stiffness will increase with increasing degree of crosslinking, but the tensile strength may go through a maximum. Gel properties can be tailored by varying not only the degree of crosslinking, but also the polymer concentration and the nature of the polymer backbone (e.g. its stiffness or solubility). Polypeptides provide an especially interesting case, where secondary structure affects trends in moduli and conformational transitions may accompany phase changes. A few papers on the tensile properties of responsive gels have begun to appear. The responsive hydrogel chosen for the present study, crosslinked hydroxypropylcellulose, shrinks over a rather narrow temperature range near 44 C. Some vertebrate skin is also subject to substantial strain. Among reptiles, the morphologies of the skin and scales show wide variations. Bauer et al. described the mechanical properties and histology of gecko skin; longitudinal tensile properties of snake skin were examined by Jayne with reference to locomotion. The present measurements focus on adaptations related to feeding, including the response of the skin to circumferential tension. Tensile properties will be related to interspecific and regional variation in skin structure and folding.

Hinkley, Jeffrey A.↗

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

Prototyping and Testing a New Volumetric Curvature Tool for Modeling Reservoir Compartments and Leakage Pathways in the Arbuckle Saline Aquifer: Reducing Uncertainty in CO 2 Storage and Permanence

This DOE-funded project evaluates the utility of seismic volumetric curvature (VC) for predicting stratal and structural architecture diagnostic of paleokarst reservoirs. Of special interest are applications geared toward carbon capture, utilization, and storage (CCUS). VC has been championed for identifying faults (offset <¼ λ) that cannot be imaged by conventional 3-D seismic attributes such as coherence. The objective of this research was to evaluate VC-techniques for reducing uncertainties in reservoir compartmentalization studies and seal risk assessments especially for saline aquifers. A 2000-ft horizontal lateral was purposefully drilled across VC-imaged lineaments—interpreted to record a fractured and a fault-bounded doline—to physically confirm their presence. The 15-mi² study area is located in southeastern Bemis-Shutts Field, which is situated along the crest of the Central Kansas Uplift (CKU) in Ellis County, Kansas. The uppermost Arbuckle (200+ ft) has extensive paleokarst including collapsed paleocaverns and dolines related to exceedingly prolonged pre-Simpson (Sauk–Tippecanoe) and/or pre-Pennsylvanian subaerial exposure. A lateral borehole was successfully drilled across the full extent (~1100 ft) of a VC-inferred paleokarst doline. Triple combo (GR-neutron/density-resistivity), full-wave sonic, and borehole micro-imager logs were successfully run to TD on drill-pipe. Results from the formation evaluation reveal breccias (e.g., crackle, mosaic, chaotic), fractures, faults, vugs (1-6"), and unaffected host strata consistent with the pre-spud interpretation. Well-rounded pebbles were also observed on the image log. VC-inferred lineaments coincide with 20–80-ft wide intervals of high GR values (100+ API), matrix-rich breccias, and faults. To further demonstrate their utility, VC attributes are integrated into a geocellular modeling workflow: 1) to constrain the structural model; 2) to generate facies probability grids, and; 3) to collocate petrophysical models to separate-vug rock fabrics along solution-enlarged fault and fracture systems. Simulation-based studies demonstrate a potential alternative field development model for developing CO 2 storage sites that target carbonate reservoirs overprinted by paleokarst. Simulation results for this complex reservoir indicate that individual fault blocks could function as discrete containers for CO 2 storage thereby reducing the risk of plume migration outside the legally defined extent of the permitted storage site. Vertically extensive, anastomosing, solution-enlarged fault/fracture systems — infilled by clay-rich sediments — would operate as non-to-low permeability vertical "curtains" that restrict CO 2 movement beyond the confines of the CO 2 storage site. Such a location could be developed in a checker-board fashion with CO 2 injection operations occurring in one block and surveillance operations occurring in the adjacent block. Such naturally partitioned reservoirs may be ideal candidates for reducing risks associated with CO 2 plume breakthrough.

02 PETROLEUM↗

High-Reliability Waveguide Vacuum/Pressure Window

The NASA Deep Space Network (DSN) uses commercial waveguide windows on the output waveguide of Ka-band (32 GHz) low-noise amplifiers. Mechanical failure of these windows resulted in an unacceptable loss in tracking time. To address this issue, a new Ka-band WR-28 waveguide window has been designed, fabricated, and tested. The window uses a slab of low-loss, low-dielectric constant foam that is bonded into a 1/2-wave-thick waveguide/flange. The foam is a commercially available, rigid, closed-cell polymethacrylimide. It has excellent electrical properties with a dielectric constant of 1.04, and a loss tangent of 0.01. It is relatively strong with a tensile strength of 1 MPa. The material is virtually impermeable to helium. The finished window exhibits a leak rate of less than 3x10(exp -3)cu cm/s with helium. The material is also chemically resistant and can be cleaned with acetone. The window is constructed by fabricating a window body by brazing a short length of WR-28 copper waveguide into a standard rectangular flange, and machining the resulting part to a thickness of 4.6 mm. The foam is machined to a rectangular shape with a dimension of 7.06x3.53 mm. The foam is bonded into the body with a two-part epoxy. After curing, the excess glue and foam are knife-trimmed by hand. The finished window has a loss of less than 0.08 dB (2%) and a return loss of greater than 25 dB at 32 GHz. This meets the requirements for the DSN application. The window is usable for most applications over the entire 26-to-40-GHz waveguide band. The window return loss can be tuned to a required frequency by var y in g the thickness of the window slightly. Most standard waveguide windows use a thin membrane of material bonded into a recess in a waveguide flange, or sandwiched between two flanges with a polymer seal. Designs using the recessed window are prone to mechanical failure over time due to constraints on the dimensions of the recess that allow the bond to fail. Designs using the sandwich method are often permeable to helium, which prohibits the use of helium leak detection. At the time of this reporting, 40 windows have been produced. Twelve are in operation with a combined operating time of over 30,000 hours without a failure.

Britcliffe, Michael J.↗

The non-Newtonian heat and mass transport of He 2 in porous media used for vapor-liquid phase separation

This investigation of vapor-liquid phase separation (VLPS) of He 2 is related to long-term storage of cryogenic liquid. The VLPS system utilizes porous plugs in order to generate thermomechanical (thermo-osmotic) force which in turn prevents liquid from flowing out of the cryo-vessel (e.g., Infrared Astronomical Satellite). An apparatus was built and VLPS data were collected for a 2 and a 10 micrometer sintered stainless steel plug and a 5 to 15 micrometer sintered bronze plug. The VLPS data obtained at high temperature were in the nonlinear turbulent regime. At low temperature, the Stokes regime was approached. A turbulent flow model was developed, which provides a phenomenological description of the VLPS data. According to the model, most of the phase separation data are in the turbulent regime. The model is based on concepts of the Gorter-Mellink transport involving the mutual friction known from the zero net mass flow (ZNMF) studies. The latter had to be modified to obtain agreement with the present experimental VLPS evidence. In contrast to the well-known ZNMF mode, the VLPS results require a geometry dependent constant (Gorter-Mellink constant). A theoretical interpretation of the phenomenological equation for the VLPS data obtained, is based on modelling of the dynamics of quantized vortices proposed by Vinen. In extending Vinen's model to the VLPS transport of He 2 in porous media, a correlation between the K*(GM) and K(p) was obtained which permits an interpretation of the present findings. As K(p) is crucial, various methods were introduced to measure the permeability of the porous media at low temperatures. Good agreement was found between the room temperature and the low temperature K(p)-value of the plugs.

Yuan, S. W. K.↗

Effects of ingested nanocellulose on intestinal microbiota and homeostasis in Wistar Han rats

Micron scale crystalline and fibrillar cellulose materials are “generally regarded as safe” (GRAS) as binders and thickeners in food products. Yet, partly due to a lack of relevant toxicological data, nanocellulose (NC) materials, which have unique properties that can be exploited to improve food quality and safety, have yet to receive FDA approval as food ingredients. Results of in vitro and in vivo toxicological studies of ingested NC, detailed in a recent companion report, revealed minimal acute in vitro cytotoxicity, and no toxicity in a subacute rat gavage model, suggesting that NC materials are non-hazardous. However, ingested materials may also modulate gut microbial populations, or alter aspects of intestinal function not evaluated in standard totoxicity testing, which could have important health implications. Here, we report the results of studies of the effects of ingested cellulose nanofibrils (CNF) on the fecal microbiome and metabolome, intestinal (ileal) epithelial cell mRNA expression of cell junction genes, and ileal cytokine production. Feces, plasma, and ileum samples were collected from Wistar Han rats before and after five weeks of biweekly gavages of water or cream, with or without 1% w/w CNF. Analysis of fecal microbial populations revealed that CNF caused changes in both genus and species diversity, with specific effects on species that produce short chain fatty acids, and that have been associated with increased IgA production and elevation of insulin levels. Fecal metabolomic analysis revealed relatively few effects of CNF, with significant changes in the levels of only ten metabolites out of 366 measured. Exposure to CNF also caused differential regulation of mRNA expression of several genes involved in epithelial cell junctions, and increased production of cytokines that both promote and inhibit proliferation of CD8 T cells. These perturbations in intestinal homeostasis contrast somewhat with the absence of in vitro and in vivo toxicity observed previously but would appear to represent minor effects. Additionally, the gut microbiome is highly sensitive to ingested substances, and such disturbances of the intestinal microbial ecosystem do not necessarily represent or predict meaningful pathology. Further studies, including chronic feeding studies, are needed to assess the real health implications, if any, of these changes.

59 BASIC BIOLOGICAL SCIENCES↗

A Review of Tungsten Heavy Alloy Utilization in Isotope Transport Containers - 13380

A common requirement for radioisotope transport containers is that they provide both durable and efficient shielding of penetrating gamma radiation. This is the case for transport of both spent nuclear fuel as well as intentionally created radioisotopes for medical or other uses. Tungsten heavy alloy (WHA) provides a unique engineering property set for such shielding - easily surpassing more commonly used lead alloys in both strength and attenuation. This family of alloys contains typically 90-98 wt.% W in combination with transition metals such as Ni and Fe. WHA is manufactured in near net shape blanks by liquid phase sintering of compacted powder shapes to full metallurgical density parts. This powder metallurgy approach is described in its ability to provide excellent material utilization and affords efficient manufacturing of various shapes required for gamma shields or collimators. WHAs offer very high density (approaching 19 g/cc) in combination with relatively high thermal conductivity, low thermal expansion, ambient corrosion resistance, and can be provided with mechanical properties comparable to many medium carbon steels. As such, they can be machined to complex, damage resistant geometries using common metal cutting tools and methods. WHA additionally provides a lower toxicity alternative to Pb- or U-based gamma shielding. Given the specialty nature of WHA, specific metallurgical characteristics are reviewed to assist shielding designers who may otherwise encounter difficulties locating important alloy selection and fabrication details. Contained within this materials and applications overview are guidelines for WHA component design, alloy selection, and practical machining, finishing, and assembly considerations. The microstructure of WHA is that of a metal matrix composite. This factor has specific implications in the design of components for stress service as well as their protection in the presence of electrolytes. WHA is also discussed in the broader context of materials compatibility, as it is rarely used in isolated monolithic shapes. Alloy selection for new applications is often made primarily on the basis of density. An alternative strategy to this selection approach is presented which proposes that mechanical requirements for a given shielding use be the primary selection criterion over density. Standard commercial grades of WHA for radiation shielding are defined by specifications such as AMS-T-21014 and ASTM B777. These specifications define 4 density classes of as-sintered WHA. Compositional options, as well as post-sinter processing of WHA, are discussed for shielding components that must exhibit higher levels of ductility or very low magnetic permeability. In addition to the mechanical advantages over Pb-based shielding, the higher linear attenuation of energetic photons for various grades of WHA (as calculated by the NIST XCOM routine) are presented for selected photon energies of interest to illustrate the shield volume reduction generally possible through the use of tungsten-based shielding. While W provides inefficient attenuation of neutrons in a mixed radiation environment, its secondary role in shielding gamma radiation produced as a result of neutron capture is also described. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Establishing a silica gel zone in well annulus and evaluating its performance in blocking vertical water flow

Wells are often constructed for monitoring purposes with relatively long screen lengths (e.g., >10 m). Vertical water flows can occur within the artificial or natural filterpack annulus that surrounds the screened interval, bypassing packer assemblies installed inside the wellbore. Attempts to isolate discrete vertical zones during groundwater sampling are unsuccessful when annular flow occurs and lead to remedy decisions based on biased or incorrect interpretations. Blocking vertical annular water flow and contaminant transport will help obtain more accurate concentrations of contaminants from sampling in targeted depth intervals. The application of silica gels formed from the injected colloidal silica CS suspensions is a novel approach to minimize or prevent movement of vertical movement of groundwater in the surrounding filterpack annulus. In this work, we tested the feasibility of injecting CS suspensions to target locations and developed a modified CS formulation that is injectable and prevents gravity sinking. We studied the distribution and penetration of silica gel at laboratory scale in mock well annulus with surrounding formations. We evaluated the performance of the silica gel in blocking vertical water flow in the annulus and in minimizing chemical transport through the gel zone. CS suspension formulations have been defined that are ready for injection, stay in target locations, and form gel within desired time frames. Injection of CS suspensions achieved uniform distribution in a well annulus filter pack, fully occupied the annulus pore space, and penetrated the formation surrounding the filter packer with a sufficient distance to create a hydraulic annular seal when the injection was applied at a sufficient rate. The depth of penetration into the formation was dependent on the permeability contrast between the filter pack and the surrounding formation. Silica gel that formed in the annulus blocked vertical water flow and stopped the chemical transport through the gel zone. In conclusion, this research reveals that using CS suspension injection and sequential gelation (CS-GEL) is a promising technology for blocking vertical water flow and chemical transport through the filter pack in targeted zones within the annulus of long-screened well systems.

Colloidal silica suspension↗

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

Carbon Dioxide Shuttling Thermochemical Storage Using Strontium Carbonate

Phase I concludes with significant progress made towards the SunShot ELEMENTS goals of high energy density, high power density, and high temperature by virtue of a SrO/SrCO 3 based material. A detailed exploration of sintering inhibitors has been conducted and relatively stable materials supported by YSZ or SrZO 3 have been identified as the leading candidates. In 15 cycle runs using a 3 hour carbonation duration, several materials demonstrated energy densities of roughly 1500 MJ/m 3 or greater. The peak power density for the most productive materials consistently exceeded 40 MW/m 3 —an order of magnitude greater than the SOPO milestone. The team currently has a material demonstrating nearly 1000 MJ/m 3 after 100 abbreviated (1 hour carbonation) cycles. A subsequent 8 hour carbonation after the 100 cycle test exhibited over 1500 MJ/m 3 , which is evidence that the material still has capacity for high storage albeit with slower kinetics. Kinetic carbonation experiments have shown three distinct periods: induction, kinetically-controlled, and finally a diffusion-controlled period. In contrast to thermodynamic equilibrium prediction, higher carbonation temperatures lead to greater conversions over a 1 hour periods, as diffusion of CO 2 is more rapid at higher temperatures. A polynomial expression was fit to describe the temperature dependence of the linear kinetically-controlled regime, which does not obey a traditional Arrhenius relationship. Temperature and CO 2 partial pressure effects on the induction period were also investigated. The CO 2 partial pressure has a strong effect on the reaction progress at high temperatures but is insignificant at temperatures under 900°C. Tomography data for porous SrO/SrCO 3 structures at initial stage and after multiple carbonation/decomposition cycles have been obtained. Both 2D slices and 3D reconstructed representations have been obtained from the raw data. Porosity and surface area have been computed based on tomography. The two parameters, permeability and Forchheimer constants, in the Darcy’s law for linear and non-linear range have been obtained from the fluid dynamics simulation. The structure is nearly isotropic from the hydrodynamics point of view. Thermal conductivity at elevated temperature has been obtained from the 3-D pore-scale simulations. The computed thermal conductivities of the structure compare well with an existing correlation. Computational models have been developed for solving the energy and mass transport equations with chemical reaction on the surface and inside the particles. The chemical kinetics or rate law is needed on the particle scale; however it is generally measured using a lumped model based on the average of reacted as well as unreacted mass. Many simulations have been conducted on a 2-D model to determine the effects of those pore scale parameters on the heat release based on the measured data. The simulation results using 2-D model show that average volumetric heat flow rate of 4 MW/m 3 in the first 5 minutes can be obtained. The insight gained from the 2-D simulation will be useful for the next stage simulation based on realistic 3-D structures. A probabilistic analysis has been conducted on the baseline system of recovering heat from the hot CO 2 evolved from the endothermic reaction and compressing the cooled CO 2 into a storage tank at 20-21 bar. Using a Monte Carlo simulation using 2,000,000 samples and both uniform and triangular probability distributions of four variables (heat exchanger pinch points, compressor efficiency, and power block efficiency), exergetic efficiency of 90% appears to be very likely.

14 SOLAR ENERGY↗

Properties and processing characteristics of low density carbon cloth phenolic composites

Ply-lift and pocketing are two critical anomalies of carbon cloth phenolic composites (CCPC) in rocket nozzle applications. Ply lift occurs at low temperatures when the A/P and in-plane permeabilities of the composite materials are still very low and in-plane porous paths are blocked. Pocketing occurs at elevated temperatures when in-plane permeability is reduced by the A/P compressive stress. The thermostructural response of CCPC in a rapid heating environment involves simultaneous heat, mass, and momentum transfer along with the degradation of phenolic resin in a multiphase system with temperature- and time-dependent material properties as well as dynamic processing conditions. Three temperature regions represent the consequent chemical reactions, material transformations, and property transitions, and provide a quick qualitative method for characterizing the thermostructural behavior of a CCPC. In order to optimize the FM5939 LDCCP (low density carbon cloth phenolic) for the nozzle performance required in the Advanced Solid Rocket Motor (ASRM) program, a fundamental study on LDCCP materials was conducted. The cured composite has a density of 1.0 +/- 0.5 gm/cc which includes 10 to 25 percent void volume. The weight percent of carbon microballoon is low (7-15 percent). However, they account for approximately one third of the volume and historically their percentages have not been controlled very tightly. In addition, the composite properties show no correlation with microballoon weight percent or fiber properties (e.g. fiber density or fiber moisture adsorption capacity). Test results concerning the ply-lift anomaly in the MNASA motor firings were: (1) Steeper ply angle (shorter path lenght) designs minimized/eliminated by lifting, (2) material with higher void volume ply lifted less frequently, (3) materials with high (greater than 9 percent) microballoon content had a higher rate of ply lifting, and (4) LDCCP materials failed at microballoon-resin interfaces. The objectives of this project are: (1) to investigate the effects of carbon microballoon and cabosil fillers as well as fiber heat treatment on plylift-related mechanical properties, (2) to develop a science-based thermostructural process model for the carbon phenolics. The model can be used in the future for the selection of the improved ASRM materials, (3) to develop the micro-failure mechanisms for the ply-lift initiation and propagation processes during the thermoelastic region of phenolic degradation, i.e. postcuring and devolatilization.

Wang, C. Jeff↗

Modeling earthquake effects on groundwater levels: evidences from the 2012 Emilia earthquake (Italy)

Changes in water level are commonly reported in regions struck by a seismic event. The sign and amplitude of such changes depend on the relative position of measuring points with respect to the hypocenter, and on the poroelastic properties of the rock. We apply a porous media flow model (TOUGH2) to describe groundwater flow and water-level changes associated with the first M L 5.9 mainshock of the 2012 seismic sequence in Emilia (Italy). We represent the earthquake as an instantaneous pressure step, whose amplitude was inferred from the properties of the seismic source inverted from geodetic data. The results are consistent with the evolution recorded in both deep and shallow water wells in the area and suggest that our description of the seismic event is suitable to capture both timing and magnitude of water-level changes. Finally, we draw some conclusions about the influence of material heterogeneity on the pore pressure evolution, and we show that to reproduce the observed maximum amplitude it is necessary to take into account compaction in the shallow layer.

2012 Emilia earthquake↗

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of the transport properties of thermally rearranged (TR) polymers and polymers of intrinsic microporosity (PIM) relative to upper bound performance

A critical analysis comparing the diffusion selectivity, solubility selectivity, diffusivity, and solubility coefficients of thermally rearranged (TR) polymers and polymers of intrinsic microporosity (PIM), two types of polymers which consistently perform at or beyond the polymer upper bound for certain gas pairs (O 2 /N 2 and CO 2 /CH 4 ), is reported in this paper. Although most polymers in these two classes exhibit transport properties in the range of typical glassy polymers, several variants offer exceptional performance. The diffusivity selectivity for O 2 /N 2 and CO 2 /CH 4 for TR polymers and PIM lies outside the range of the glassy polymer database. The solubility coefficients for O 2 are at the upper end of the range of both TR and PIM polymers, as is also the case for CO 2 solubility coefficients for PIM polymers. The O 2 /N 2 and CO 2 /CH 4 solubility selectivities for both PIM and TR polymers are in the range of typical glassy polymers. Thus, unique separation values for these polymers are a manifestation of maximizing diffusivity selectivity combined with very high gas solubility. Additionally, the previously determined diffusion gas kinetic diameters, d g , were found to correlate best with diffusion coefficients of PIM and TR polymers and serve to better analyze transport properties of upper bound materials.

36 MATERIALS SCIENCE↗