Characterization of Mason Gully (H5): The second recovered fall from the Desert Fireball Network
Abstract not provided
SEARCH · Engineering Papers
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.
Abstract not provided
Pelletized magnesium manganese oxide shows promise for high temperature thermochemical energy storage. It can be thermally reduced in the temperature range between 1250 °C and 1500 °C and re-oxidized with air at typical gas-turbine inlet pressures (1–25 bar) in the temperature range between 600 °C and 1500 °C. The combined thermal and chemical volumetric energy density is approximately 2300 MJ/m3. The rate at which a thermochemical storage module can be charged is limited by heat transfer inside the solid packed bed. Hence, the effective thermal conductivity of packed beds of magnesium-manganese oxide pellets is a crucial parameter for engineering Mg-Mn-O redox storage devices. We have measured the effective thermal conductivity of a packed bed of 3.66 ± 0.516 mm sized magnesium manganese oxide (Mn to Mg molar ratio of 1:1) pellets in the temperature range of 300–1400 °C. Since the material is electrically conductive at temperatures above 600 °C, the sheathed transient hot wire method is used for measurements. Raw data is analyzed using the Blackwell solution to extract the bed thermal conductivity. The effective thermal conductivity standard deviation is less than 10% for a minimum of three repeat measurements at each temperature. Experimental results show an increase in the effective thermal conductivity with temperature from 0.50 W/m °C around 300 °C to 1.81 W/m °C close to 1400 °C. We propose a dual porosity model to express the effective thermal conductivity as a function of temperature. This model also considers the effect of radiation within the bed, as this is the dominant heat transfer mode at high temperatures. The proposed model accounts for microscale pellet porosity, macroscale bed porosity, pellet size, solid thermal conductivity (phonon transport), and radiation (photon transport). The coefficient of determination between the proposed model and the experimental results is greater than 0.90.
Calendering serves as a multifunctional step in dry electrode processing that not only densifies the electrode but also induces polytetrafluoroethylene (PTFE) fibrillation and reorganizes the microstructure. These coupled effects are essential for achieving electrical connectivity and sufficient cohesion, yet they also introduce trade-offs, such as active material particle fracture, pore collapse, and excessive porosity loss, that can hinder ionic transport. This research systematically maps the calendering parameter space for LiNi0.6Mn0.2Co0.2O2 (NMC622) dry cathodes with a target thickness of ∼100 µm and porosity of ∼30% by varying roll gaps, roll temperature, roll speed, and the number of passes. A practical processing window for this formulation and electrode architecture is identified that achieves sufficient PTFE fibrillation and strong interfacial contact while minimizing particle fracture and preserving the porosity required for efficient ionic transport. In particular, gradual-gap calendering with moderate per-pass compression mitigates fracture and pore collapse while still reaching the target thickness with reasonable throughput, and slower roll speeds with modest roll temperatures further reduce mechanical damage. These results provide actionable guidance for scaling NMC622-based thick dry-processed cathodes.
Abstract The design and fabrication of Inconel 718 open-pore lattice structures via Laser Powder Bed Fusion (LPBF) has been investigated in this research, focusing on applications such as catalyst supports in jet fuel production. The study explores the impact of laser power and scanning speed on the geometrical resolution of these structures aiming to achieve high porosity (porosity > 60%) and specific pore sizes ranging from 500–1000 μm, intending to serve as catalyst supports, replacing conventionally manufactured foams to reduce costs. Results demonstrate the significant influence of processing parameters on the geometrical aspects of printed lattice structures, with laser power having a more pronounced effect on geometrical accuracy than scanning speed. Additionally, the mechanical properties of the printed lattice structures showed a correlation with the lattice strut sizes, as lattices with less porosity and thicker struts resulted in higher maximum shear stress.
Porosity defects are currently a major factor that hinders the widespread adoption of laser-based metal additive manufacturing technologies. One common porosity occurs when an unstable vapor depression zone (keyhole) forms because of excess laser energy input. With simultaneous high-speed synchrotron x-ray imaging and thermal imaging, coupled with multiphysics simulations, we discovered two types of keyhole oscillation in laser powder bed fusion of Ti-6Al-4V. Amplifying this understanding with machine learning, we developed an approach for detecting the stochastic keyhole porosity generation events with submillisecond temporal resolution and near-perfect prediction rate. Finally, the highly accurate data labeling enabled by operando x-ray imaging allowed us to demonstrate a facile and practical way to adopt our approach in commercial systems.
This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.
Miscanthus is a leading bioenergy crop with enormous lignocellulose production potential for biofuels and chemicals. However, lignocellulose recalcitrance leads to biomass process difficulty for an efficient bioethanol production. Hence, it becomes essential to identify the integrative impact of lignocellulose recalcitrant factors on cellulose accessibility for biomass enzymatic hydrolysis. In this study, we analyzed four typical pairs of Miscanthus accessions that showed distinct cell wall compositions and sorted out three major factors that affected biomass saccharification for maximum bioethanol production. Among the three optimal (i.e., liquid hot water, H 2 SO 4 and NaOH) pretreatments performed, mild alkali pretreatment (4% NaOH at 50 °C) led to almost complete biomass saccharification when 1% Tween-80 was co-supplied into enzymatic hydrolysis in the desirable Miscanthus accessions. Consequently, the highest bioethanol yields were obtained at 19% (% dry matter) from yeast fermentation, with much higher sugar–ethanol conversion rates by 94–98%, compared to the other Miscanthus species subjected to stronger pretreatments as reported in previous studies. By comparison, three optimized pretreatments distinctively extracted wall polymers and specifically altered polymer features and inter-linkage styles, but the alkali pretreatment caused much increased biomass porosity than that of the other pretreatments. Based on integrative analyses, excellent equations were generated to precisely estimate hexoses and ethanol yields under various pretreatments and a hypothetical model was proposed to outline an integrative impact on biomass saccharification and bioethanol production subjective to a predominate factor (CR stain) of biomass porosity and four additional minor factors (DY stain, cellulose DP, hemicellulose X/A, lignin G-monomer). Using four pairs of Miscanthus samples with distinct cell wall composition and varied biomass saccharification, this study has determined three main factors of lignocellulose recalcitrance that could be significantly reduced for much-increased biomass porosity upon optimal pretreatments. It has also established a novel standard that should be applicable to judge any types of biomass process technology for high biofuel production in distinct lignocellulose substrates. Hence, this study provides a potential strategy for precise genetic modification of lignocellulose in all bioenergy crops.
Shallow permafrost distribution and characteristics are important for predicting ecosystem feedbacks to a changing climate over decadal to century timescales because they can drive active layer deepening and land surface deformation, which in turn can significantly affect hydrologic and biogeochemical responses, including greenhouse gas dynamics. As part of the U.S. Department of Energy Next-Generation Ecosystem Experiments- Arctic, we have investigated shallow Arctic permafrost characteristics at a site in Barrow, Alaska, with the objective of improving our understanding of the spatial distribution of shallow permafrost, its associated properties, and its links with landscape microtopography. To meet this objective, we have acquired and integrated a variety of information, including electric resistance tomography data, frequency-domain electromagnetic induction data, laboratory core analysis, petrophysical studies, high-resolution digital surface models, and color mosaics inferred from kite-based landscape imaging. The results of our study provide a comprehensive and high-resolution examination of the distribution and nature of shallow permafrost in the Arctic tundra, including the estimation of ice content, porosity, and salinity. Among other results, porosity in the top 2 m varied between 85% (besides ice wedges) and 40%, and was negatively correlated with fluid salinity. Salinity directly influenced ice and unfrozen water content and indirectly influenced the soil organic matter content. A relatively continuous but depth-variable increase in salinity led to a partially unfrozen saline layer (cryopeg) located below the top of the permafrost. The cryopeg environment could lead to year-round microbial production of greenhouse gases. Results also indicated a covariability between topography and permafrost characteristics including icewedge and salinity distribution. In addition to providing insight about the Arctic ecosystem, through integration of lab-based petrophysical results with field data, this study also quantified the key controls on electric resistivity at this Arctic permafrost site, including salinity, porosity, water content, ice content, soil organic matter content, and lithologic properties.
Geologic carbon sequestration involves the injection of captured carbon dioxide ([Formula: see text]) into subsurface formations for long-term storage. The movement and fate of the injected [Formula: see text] plume is of great concern to regulators because monitoring helps to identify potential leakage zones and determines the possibility of safe long-term storage. To address this concern, we design a deep-learning framework for [Formula: see text] saturation monitoring to determine the geologic controls on the storage of the injected [Formula: see text]. We use different combinations of porosities and permeabilities for a given reservoir to generate saturation and velocity models. We train the deep-learning model with a few time-lapse seismic images and their corresponding changes in saturation values for a particular [Formula: see text] injection site. The deep-learning model learns the mapping from the change in the time-lapse seismic response to the change in [Formula: see text] saturation during the training phase. We then apply the trained model to data sets comprising different time-lapse seismic image slices (corresponding to different time instances) generated using different porosity and permeability distributions that are not part of the training to estimate the [Formula: see text] saturation values along with the plume extent. Our algorithm provides a deep-learning assisted framework for the direct estimation of [Formula: see text] saturation values and plume migration in heterogeneous formations using the time-lapse seismic data. Our method improves the efficiency of time-lapse inversion by streamlining the large number of intermediate steps in the conventional time-lapse inversion workflow. This method also helps to incorporate the geologic uncertainty for a given reservoir by accounting for the statistical distribution of porosity and permeability during the training phase. Tests on different examples verify the effectiveness of our approach.
Calendering is a technique used to maximize the volumetric energy density of battery electrodes. However, higher amounts of calendering result in increased tortuosity and particle cracking. We propose a novel packing structure of electrode particles to maximize calendering benefits while minimizing particle fracture. Cobalt-free layered oxide cathode LiNi0.92Mn0.04Al0.04O2 (NMA) particles are pulverized through ball-milling and coated with lithium phosphate. Pulverized and pristine NMA are fabricated into “bimodal” electrodes, whereas “unimodal” electrodes consist of only pristine NMA. Each electrode type was made into 30% porosity, 40% porosity, and uncalendered coin cell samples. X-ray diffraction suggests that the unimodal samples suffer from more particle fracture than the bimodal samples when calendered to the same porosity of 30%. Electrochemical impedance spectroscopy suggests that the bimodal electrodes exhibit lower surface film resistance. This is supported by enhanced capacity retention for the bimodal samples after 100 cycles.
The study of propped fracture conductivity began in earnest with the development of the Cooke cell which later became part of the initial API standard. Subsequent developments included a patented multicell design to conduct 4 tests in a press at the same time. Other modifications have been used by various investigators. Recent studies by the Stim-Lab proppant consortium have indicated that the flow field across a Cooke proppant conductivity testing cell may not be uniform as initially believed which resulted is significantly different conductivity results. Post test analysis of low temperature metal alloy injections at the termination of proppant testing prior to the release of the applied stress suggest that higher flow is to be expected along the sides and top of the proppant pack than compared to the middle of the pack. To evaluate these experimental findings, a physics-based two-dimensional (2-D) discrete element model (DEM) was developed and applied to simulate proppant rearrangement during stress loading in the Cooke conductivity cell and the resulting porosity field. Analysis of these simulations are critical to understanding the impact of modification to the testing cell as well as understanding key proppant conductivity issues such as how these effects are manifested in proppant concentration testing results. The 2-D DEM model was constructed to represent a realistic cross section of the Cooke cell with a distribution of four material properties, three that represented the Cooke cell (steel, sandstone,square rings), and one representing the proppant. In principle, Cooke cell materials can be approximated as assemblies of independent discrete elements (particles) of various sizes and material properties that interact via cohesive interactions, repulsive forces, and frictional forces. The macroscopic behavior can then be modeled as the collective behavior of many interacting discrete elements. This DEM model is particularly suitable for modeling proppant mechanical interactions subjected to an applied stress, where the experimental cell is represented as a cohesive body composed of a large number of discrete elements, and proppants can be modeled as the individual discrete particles with various sizes (following the proppant size distribution-density function used in the test) that exhibit no cohesive strength between the particles. Initial 2-D DEM modeling results suggest that proppant rearrangement and non-uniform stress distribution across the proppant pack results in significant non-uniform porosity distribution across the Cooke cell. Larger porosities develop along the edge of the proppant pack beneath the square ring seal and would result in a disproportionate higher flow field along these edges as compared to the middle of the proppant pack. These results suggest that reported conductivity values determined by the Cooke cell may be biased to overestimate the actual conductivity of the proppant at high stresses and that modifications to the standard Cooke cell will affect the magnitude of this bias.
Understanding proppant-transport and -deposition patterns in a hydraulic fracture is vital for effective and economical production of petroleum hydrocarbons. In this research, a numerical-modeling approach, combining the discrete-element method (DEM) with single-/multiphase lattice Boltzmann (LB) simulation, was adopted to advance understanding of the interaction between reservoir depletion, proppant-particle compaction, and single-/multiphase flows in a hydraulic fracture. DEM was used to simulate effective-stress increase and the resultant proppant-particle movement and rearrangement during the process of reservoir depletion caused by hydrocarbon production. Simulated pore structure of the proppant pack was extracted and used as boundary conditions in the LB simulation to calculate the time-dependent permeability of the proppant pack. We first validated the DEM-LB numerical work flow, and the simulated proppant-pack permeabilities as functions of effective stress were in good agreement with laboratory measurements. Furthermore, three proppant packs with the same average particle diameter but different diameter distributions were generated to study the role of proppant-size heterogeneity (variation in particle diameter). Specifically, we used the coefficient of variation (COV) of diameter, defined as the ratio of standard deviation of diameter to mean diameter, to characterize the heterogeneity of particle size. We obtained proppant-pack porosity, permeability, and fracture-width reduction (closure distance) as functions of effective stress. Under the same effective stress, a proppant pack with a higher diameter COV had lower porosity and permeability and larger fracture-width reduction. This was because the high diameter COV gave rise to a wider diameter distribution of proppant particles; smaller particles were compressed into the pore space between larger particles with increasing stress, leading to larger closure distance and lower porosity and permeability. With multiphase LB simulation, relative permeability curves were obtained, which are critical for larger-scale reservoir simulations under various effective stresses. The relative permeability of the oil phase was more sensitive to changes in geometry and stress, compared with the wetting phase. This was because the oil phase occupied larger pores; compaction of the proppant pack affected the structure of the pores, because the pores were farther from the grain contacts and thus more susceptible to collapse. It is also interesting to note that when effective stress increased continuously, oil relative permeability increased first and then decreased. This nonmonotonic behavior was the result of the nonmonotonic development of pore structure and oil connectivity under increasing stress.
Summary This study continues the work of presenting a novel approach for making petrophysical assessments of tight core samples. This method, the full-immersion pressure-pulse decay (Hannon 2019), involves applying a rapid increase in pressure in a chamber surrounding the entire outer surface area of a cylindrical sample, shutting the system in, and monitoring the pressure decay in the chamber as it reaches a new equilibrium. A precursor article covered the numerical simulator designed to model flow through the sample, demonstrating its performance and accuracy in addition to providing a first-order comparison between the speed and shape of the pressure-decay responses of the full-immersion method with those of other similar transient methods. This study covers the parameter-estimation procedure and experimental verification through a proof-of-concept laboratory investigation. The investigations provided here demonstrate that under appropriate, achievable experimental conditions, the pressure data can be analyzed in such a way that returns an estimate of the porosity and apparent permeabilities both parallel and perpendicular to bedding from a single test performed on a single cylindrical sample. After determining these experimental conditions (the uniqueness window), this report outlines a data-inversion strategy to estimate the petrophysical properties (porosity, horizontal permeability, and vertical permeability) from each test. This strategy is put to the test through comparisons with measurements performed by a commercial core laboratory. A common set of samples recovered from an outcrop of a tight-gas sandstone formation were investigated using the full-immersion method, and their results are compared with those from conventional steady-state measurement procedures performed by the commercial laboratory. Comparisons between petrophysical characterizations of these samples, which had permeabilities between 25 nd and 2.3 μd, demonstrated close agreement in most cases. However, whereas steady-state measurements performed at the professional laboratory required 4 to 5 hours of testing time per measurement of a single permeability, similar assessments using the full-immersion technique, requiring approximately 5 to 10 minutes to complete, returned estimates of the horizontal and vertical permeability simultaneously. Additional analyses are provided to determine principal reasons of discrepancies in instances where agreement was not as strong. Based on lessons learned from these experiences, the report closes with suggestions on areas of improvement in the experimental approach. Once complete, these developments should propel this technology to fill a critical need to determine petrophysical properties (porosity and permeability) of tight rocks in a time-efficient manner and in a way that does not compromise their accuracy.
Summary Carbon dioxide (CO2) injection into a deep saline aquifer can dissolve into formation brine and generate carbonic acid. The resulting acid can drive fluid-rock geochemical reactions. The impact of these fluid-rock geochemical reactions on porosity, permeability, and multiphase flow responses is relevant to the determination of CO2 storage capacity of deep saline aquifers. In this research, carbonic acid flooding experiments were performed on core samples consisting of poorly sorted, quartz-rich sand with laminated bedding from a possible CO2 storage target in northwest Wyoming. Complementary pre- and post-injection porosity and permeability, thin-section, Brunauer-Emmett-Teller (BET) surface area, mercury intrusion capillary pressure (MICP), and time-domain nuclear magnetic resonance (TD-NMR) measurements were conducted. Overall, both core porosity and permeability increased after a 7-day carbonic acid injection, from 6.2 to 8.4% and 1.6 to 3.7 md, respectively. We attributed these changes to carbonate mineral dissolution, which was evidenced by the effluent brine geochemistry, pore-throat size distribution (PTSD), and BET surface area. To be more specific, within the more-permeable section of core samples containing larger pore size, the permeability increment is apparent due to dolomite mineral grains and cements dissolution. However, for the lower-permeability section corresponding to the smaller pore size, mineral precipitation possibly lessened dissolution effects, leading to insignificant petrophysical properties changes. Consequently, the observed heterogeneous carbonic acid-rock interactions resulted in alterations of CO2/brine relative permeability (i.e., the initial CO2 saturation decreased and the CO2 flow capacity was enhanced). This research provides a fundamental understanding regarding effects of fluid-rock reactions on changes in static and multiphase flow properties of eolian sandstones, which lays the foundation for more accurate prediction/simulation of CO2 injection into deep saline aquifers.
Core Ideas A field‐scale experiment and numerical simulations confirm salt backfill behavior. Simulations closely match temperature around and under the piled salt backfill. Results indicate limited dissolution–precipitation reactions around the heat source. Alteration of backfill is unlikely if the drift is allowed to dry before emplacement. Investigations relating to in‐drift disposal of heat‐generating nuclear waste in salt have raised questions about heat–brine interactions in the unsaturated run‐of‐mine (RoM) salt pile used as backfill. These interactions have the potential to change the structure of the RoM salt surrounding the canister, possibly altering long‐term containment of the source. An experiment is in progress at the Waste Isolation Pilot Plant (WIPP), New Mexico, in which a heated canister was placed on the floor of an open drift, covered in a pile of RoM salt, and energized with 1000 W. Temperature in the RoM salt pile had stabilized after about 15 d, allowing evaluation of the heat‐up period of the ongoing experiment. Using a multiphase porous flow simulator that has been modified to handle salt‐specific coupled processes, we examined coupled thermal–hydrological–chemical behavior in the RoM salt pile. Our simulations suggest that for the relatively dry cases examined, porosity changes within RoM salt in a generic salt repository are likely to be minor in the period between waste emplacement and plastic closure of the drift. The primary sensitivity for porosity change is to the early moisture content of the RoM salt used to cover the canister. Secondary influences include moisture availability from the disturbed rock zone (DRZ) surrounding the drift and the capillary pressure ratio between the DRZ and the RoM salt. Early changes in porosity and permeability may be affected by moisture content, but this was not observed in the test. Such changes would be most likely to occur when using damp RoM salt or if waste is emplaced in a drift immediately following opening of the drift before evaporative dewatering of the drift walls occurs.
This article's purpose is to review some of the recent research in which geochemists have examined precipitation of solid phases in porous media, particularly in pores a few nanometers in diameter (nanopores). While this is a “review,” it is actually more forward-looking in that the list of things about this phenomenon that we do not know or cannot control at this time is likely longer than what we do know and can control. For example, there are three directly contradictory theories on how to predict how precipitation proceeds in a medium of varying pore size, as will be discussed below. The confusion on this subject likely stems from the complexity of the phenomenon itself: One can easily clog a porous medium by inducing a rapid, homogeneous precipitation directly from solution, or have limited precipitation occur that does not affect permeability or even porosity substantially. It is more difficult to engineer mineral precipitation in order to obtain a specific outcome, such as filling all available pore space over a targeted area for the purposes of contaminant sequestration. However, breakthrough discoveries could occur in the next five to ten years that enhance our ability to predict robustly and finely control precipitation in porous media by understanding how porosity and permeability evolve in response to system perturbations. These discoveries will likely stem (at least in part) from advances in our ability to 1) perform and interpret X-ray/neutron scattering experiments that reveal the extent of precipitation and its locales within porous media (Anovitz and Cole 2015, this volume), and 2) utilize increasingly powerful simulations to test concepts and models about the evolution of porosity and permeability as precipitation occurs (Steefel et al. 2015, this volume). A further important technique to isolate specific phenomena and understand reactivity is also microfluidics cell experiments that allow specific control of flow paths and fluid velocities (Yoon et al. 2012). An improved ability to synthesize idealized porous media will allow for tailored control of pore distributions, mineralogy and will allow more reproducible results. This in turn may allow us to isolate specific processes without the competing and obfuscatory effects that hinder generalization of observations when working with solely natural samples. It is likely that no one single experiment, or simulation technique will provide the key discoveries: to make substantive progress will require a collaborative effort to understand the interplay between fluid transport and geochemistry. Finally, where rock fracturing and elevated pressures are of concern, an understanding and capability to model geomechanical properties are necessary (Scherer 1999).
An integrated geologic and geophysical study of the Bakken Petroleum System, in the Williston basin of North Dakota and Montana indicates that: (1) dolomite is needed for good reservoir performance in the Middle Bakken; (2) regional and local fractures play a significant role in enhancing permeability and well production, and it is important to recognize both because local fractures will dominate in on-structure locations; and (3) the organic-rich Bakken shale serves as both a source and reservoir rock. Results from the lithofacies, mineral, and fracture analyses of this study were used to construct a dual porosity Petrel geo-model for a portion of the Elm Coulee Field. In this field, dolomitization enhances reservoir porosity and permeability. First year cumulative production helps locate areas of high well productivity and in deriving fracture swarm distribution. A fracture model was developed based on high productivity well distribution, and regional fracture distribution, and was combined with favorable matrix properties to build a dual porosity geo-model.
Modeling of reacting flows in porous media has become particularly important with the increased interest in hydrogen solid-storage beds. An advanced type of storage bed has been proposed that utilizes oxidation of uranium hydride to heat and decompose the hydride, releasing the hydrogen. To reduce the cost and time required to develop these systems experimentally, a valid computational model is required that simulates the reaction of uranium hydride and oxygen gas in a hydrogen storage bed using multiphysics finite element modeling. This SAND report discusses the advancements made in FY12 (since our last SAND report SAND2011-6939) to the model developed as a part of an ASC-P&EM project to address the shortcomings of the previous model. The model considers chemical reactions, heat transport, and mass transport within a hydride bed. Previously, the time-varying permeability and porosity were considered uniform. This led to discrepancies between the simulated results and experimental measurements. In this work, the effects of non-uniform changes in permeability and porosity due to phase and thermal expansion are accounted for. These expansions result in mechanical stresses that lead to bed deformation. To describe this, a simplified solid mechanics model for the local variation of permeability and porosity as a function of the local bed deformation is developed. By using this solid mechanics model, the agreement between our reacting bed model and the experimental data is improved. Additionally, more accurate uranium hydride oxidation kinetics parameters are obtained by fitting the experimental results from a pure uranium hydride oxidation measurement to the ones obtained from the coupled transport-solid mechanics model. Finally, the coupled transport-solid mechanics model governing equations and boundary conditions are summarized and recommendations are made for further development of ARIA and other Sandia codes in order for them to sufficiently implement the model.