Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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 685 records · Page 38

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SECARB-USA: Needs Assessment Framework for Storage Complexes (Task 2.1.b)

A team of SECARB storage experts examined 63 formations at 39 sites that were selected to represent the range of diversity of newly assessed, as well as well-advanced, storage prospects in the SECARB region. We inventoried the data needs triggered by the requirements to obtain a Class VI UIC storage permit, the data needs that results from the requirement to create geocellular fluid flow models to support that permit, and by pragmatic and best practice inputs such as public acceptance, regulatory readiness, and pore space leasing. We anchored both the needs inventory and the processes for and cost of meeting the needs with data from 9 sites in in the SECARB area that have advanced far in characterization. Results show that the total cost of characterization prior to obtaining a permit is convergent, because the permit and modeling requirements drive projects to obtain the same types of data for all cases. The high cost data that control cost are 1) drilling, coring, core-testing, logging, sampling and testing a characterization well and 2) collection of a 3-D seismic survey to map reservoir and confining system properties over the area of the plume or the area of elevated pressure. In 5 of our case study sites we determined that one or both of these costs could be avoided because the needs are met by available data.

54 ENVIRONMENTAL SCIENCES↗

Advanced single crystal for SSME turbopumps

The objective of this program was to evaluate the influence of high thermal gradient casting, hot isostatic pressing (HIP) and alternate heat treatments on the microstructure and mechanical properties of a single crystal nickel base superalloy. The alloy chosen for the study was PWA 1480, a well characterized, commercial alloy which had previously been chosen as a candidate for the Space Shuttle Main Engine high pressure turbopump turbine blades. Microstructural characterization evaluated the influence of casting thermal gradient on dendrite arm spacing, casting porosity distribution and alloy homogeneity. Hot isostatic pressing was evaluated as a means of eliminating porosity as a preferred fatigue crack initiation site. The alternate heat treatment was chosen to improve hydrogen environment embrittlement resistance and for potential fatigue life improvement. Mechanical property evaluation was aimed primarily at determining improvements in low cycle and high cycle fatigue life due to the advanced processing methods. Statistically significant numbers of tests were conducted to quantitatively demonstrate life differences. High thermal gradient casting improves as-cast homogeneity, which facilitates solution heat treatment of PWA 1480 and provides a decrease in internal pore size, leading to increases in low cycle and high cycle fatigue lives.

Fritzemeier, L. G.↗

Single-Wall Carbon Nanotube Production by the Arc Process: A Parametric Study

Single wall carbon nanotubes are produced using the arc discharge process. Graphite anodes are filled with a mixture of nickel and yttrium metallic powders, then vaporized by creating a high current arc. By varying the current, gap distance, and ambient pressure it is shown that the best yield of single wall carbon nanotubes is obtained within a narrow range of conditions. The relative yield and purity of the product are indicated semi-quantitatively from scanning electric microscopy (SEM) and thermogravimetric analysis (TGA). Two types of anodes have been investigated. The first is hollow and filled with a powder mixture of graphite, nickel and yttrium. The second is filled with a paste made of a mixture of metal nitrates, graphite powder and carbon adhesive, then reduced in an argon atmosphere at high temperature. Product purity and yield will be compared for the two types of anodes. The graphite in the anodes may have hydrogen attached in the pores. To remove this impurity anodes have been baked up to 1400 - 1500 C. The effect of baking the anodes on impurities in the product will be given.

Scott, Carl D.↗

In situ investigation of phosphonate retarder interaction in oil well cements at elevated temperature and pressure conditions

In this paper, the effect of a high-performance retarding additive in oil well cements was investigated under elevated temperature (165°C) and pressure (1000 psi) conditions via in situ synchrotron-based X-ray diffraction (XRD) and quasielastic neutron scattering (QENS) techniques. Under these temperature and pressure conditions, crystalline calcium silicate hydrates (C–S–H) are formed through the cement hydration process. From in situ XRD experiments, the retardation effect was observed by a change in the rate of the appearance of 11 Å tobermorites as well as a change in the rate of the α-C 2 SH generation and depletion. QENS analysis revealed that the retardation effect was related to the non-conversion of free water to chemical and constrained water components. A high presence of free water components was attributed to a decrease in 11Å tobermorites along with slower consumption of the quartz and portlandite phases. Furthermore, QENS results infer that the water molecules experienced confinement in the restricted pore spaces. The retarder inhibited this initial water confinement by slowing the bulk diffusion of free water in the confined region.

36 MATERIALS SCIENCE↗

3D Printed Carbon Aerogels for Polymer-Electrolyte Fuel Cells

Carbon aerogels (CA) are known for ultra-low density, very high surface area, fine porous structure, and excellent electrical conductivity. Recent advances in high resolution 3D printing at LLNL enable the manufacturing of complex structured CA electrodes which have been demonstrated in electrolysis, flow battery, and supercapacitor applications. Polymer-electrolyte fuel cells (PEFC) may be an excellent candidate for these materials/methods as designers seek to optimize membrane-electrode assemblies to achieve higher efficiencies. A major challenge to PEFC design is two-phase flow due to water production which causes pore flooding and leads to cell instability and degradation. To mitigate these effects, groups have investigated the creation of liquid specific pathways in the porous layers to alleviate liquid buildup and enhance gas diffusion. While some progress has been made, traditionally manufactured GDLs, which consist of randomly distributed carbon fibers, are limiting. 3D printed CAs are a promising solution as complex macroporous electrode structures can be manufactured out of micro and nano porous filaments with a high level of control. Additionally, the process may allow for the consolidation of multiple layers (MPL, GDL, & flow channels) into a single printed architecture. In this work various CA lattice structures were printed using a direct ink write (DIW) method and liquid breakthrough pressures were measured to demonstrate the ability to tune the water pathways to achieve specific pressure thresholds. Furthermore, a discussion of the prospects for printed CA electrodes in PEFC will focus on material properties, synthesis techniques, and current challenges.

36 MATERIALS SCIENCE↗

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗

Gas-phase Radiolysis in Plutonium Dioxide Powder

A sample of PuO 2 powder undergoes self radiolysis, causing water present in the solid to generate hydrogen, oxygen, and minor species (e.g. H 2 O 2 ). This is a matter of concern for stored PuO 2 , since the gases could react exothermically and rupture the container. In experiments at LANL, most PuO 2 /salt samples generate mostly H 2 , but a few samples generate a mixture of H 2 and O 2 together. In these special cases, the mixture of H 2 and O 2 in the pore volumes of the powder will then be subject to radiolysis, causing the reverse reaction H 2 + O 2 $\Rightarrow$ H 2 O. In order to understand the steady-state quantities of these gases, it is necessary to estimate the amount of energy from the alpha irradiation that drives this reaction. This report outlines one possible approach for this process, yielding a value for the fractional amount of radiation deposited in the gas phase of 294 ppm. This corresponds to the case of a powder with porosity of 78%. We then use a comprehensive chemical kinetics model based on elementary gas phase reactions to calculate the rate of conversion back to water. We estimate a reduction in the H 2 pressure of about 12 kPa/month, which applies only when O 2 is simultaneously present in the gas phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

Consolidation and Permeability of the B1 and D1 Gas Hydrate Bearing Sands and Associated Seal Sediments of the Extended-Duration Gas Production Test Site on the Alaska North Slope

Gas hydrate, a solid combination of gas (mostly methane in nature) and water molecules stable at low temperatures and elevated pressures, occurs naturally in marine and permafrost-associated environments. Gas hydrate reservoirs, such as those in the Alaska North Slope, have been considered potential energy resources for gas production. To understand the petrophysical and geo-mechanical characteristics of the reservoir, core samples retrieved from the site of the JOGMEC-DOE-USGS collaborative gas hydrate R&D project have been analyzed in the laboratory for their hydraulic and mechanical properties. This paper focuses on both seal and reservoir samples associated with the B1 and D1 sands, which are evaluated for index properties (including porosity, grain size distribution, liquid and plastic limits, specific surface area, and specific gravity), consolidation, permeability, and water retention. Furthermore, the reservoir core samples were tested with pore-filling, laboratory-grown tetrahydrofuran hydrate, in order to assess reservoir behavior during gas production from hydrates. Under simulated in situ stress conditions, the seal and hydrate-free reservoir cores had a permeability anisotropy ratio of k h /k v = 3.0−5.0, and k h /k v = 2.4−3.0 for the reservoir tetrahydrofuran hydrate-bearing cores. The data suggest that depressurizing the reservoir to induce hydrate dissociation alters the reservoir effective permeability in three ways: permeabilities decrease due to porosity lost (e.g., the initial reservoir thickness can decrease by up to 5% upon 7 MPa depressurization), permeability increases due to the loss of solid hydrate in the pore space, and permeability anisotropy k h /k v decreases in response to the evolving pore-space geometry. We show that given the simulated in situ gas hydrate saturations (i.e., S h = 32% in core 7P-2E and S h = 21% in core 20P-4), gas production from the dissociation of tetrahydrofuran hydrate in the two tested cores results in a net increase in effective permeability and a decrease in k h /k v . This study highlights the importance of investigating seal and reservoir sediments and the impacts of depressurization on the porosity and permeability responses during production.

Geological materials↗

The Magnetic Properties Of Aggregate Polycrystalline Diamond: Implications For Carbonado Petrogenesis

Carbonados are a type of diamond, which are made up of many aggregrates of small crystalline diamonds or microdiamonds. The term "carbonado" comes from the Portuguese word carbonated. They are only found in sedimentary deposits in the Central African Republic (CAR) and the Bahia Province of Brazil. They were once the source of the world's supply of industrial diamonds. Their origin is uncertain but several mutually exclusive hypotheses have been proposed. This theories are: (1) extraterrestrial, that is they formed from the dust cloud of original solar nebulae; (2) produced by the high temperatures and pressures of the Earth's mantle; (3) or as the result of an extra-terrestrial impact into a carbon rich layer of sediment. Our study was done to further the understanding of their origin. We measured the magnetic properties on some twenty samples from the CAR. An earlier study was done on whole samples of carbonados and the "common" or kimberlitic diamond. Our work differed in that we started at the surface and subsequently removed the surface layers (by days of acid immersion) into the interior; measuring the magnetic properties at each interval. This procedure permits us to monitor the distribution of magnetic substances within the samples. Our results showed that the magnetic carriers are distributed on the surface including the open pores and that the carbonado interior is essentially non-magnetic. This result suggests that the initial formation environment was deficient in magnetic particles. Such a situation could indicate that their formation was the result of an extra-terrestrial body impacting carbon-rich sediment. Obviously, more work will be required on isotopic and chemical analyses before a more detailed ori-in can be determined.

Kleteschka, Gunther↗

Ionic liquids for carbon capture

This article covers cutting-edge research in ionic liquid-based materials and their impact on carbon capture. Ionic liquids boast high thermal stability, low vapor pressure, and great tunability. They can be deployed as either sorbents or supported liquid membranes for carbon capture and other gas separations. Absorption of CO 2 in ionic liquids can be achieved via physisorption or chemisorption, depending on the strength of the interaction. For physisorption, both CO 2 -anion and cation–anion interactions are important factors in dictating CO 2 solubility. For chemisorption, the reactivity with CO 2 can be tuned by the basicity of the anion while the stability of the ionic liquid can be enhanced by controlling the cation. Marrying ionicity of ionic liquids with porosity leads to novel porous ionic liquids that feature unique properties for facilitating gas transport, while interfacing ionic liquids with microporous membranes allow the ions to gate the pores for selective gas transport. Further, composite materials of ionic liquids with crystalline materials such as metal–organic frameworks offer scalable heterogeneous interfaces to promote selectivity. Future opportunities include controlling the interaction with CO 2 and its transport at the charged or electrified interface with the ionic-liquid electric double layer, as well as combining CO 2 capture and conversion with ionic liquids.

36 MATERIALS SCIENCE↗

Heteroatom-Modified and Compacted Zeolite-Templated Carbons for Gas Storage

Methane, the primary component of natural gas, is an energy-rich fuel by weight but exhibits lower energy density by volume than diesel and gasoline; efforts to increase its “volumetric” energy density have been focused on the design and synthesis of a porous scaffold that can bind the methane favorably. Novel approaches typically invoke the highly designable class of materials known as metal-organic frameworks (MOFs) owing to the relative difficulty of controlling carbon structure and composition precisely, without the use of a metal-based node. In this project, we instead explore the concept of carbon templating to achieve a porous scaffold with atomistically thin walls separated at a distance of ~1.2 nm, ideal for the formation of precisely two layers of methane per pore. This bottom-up templating approach permits the use of any small molecular precursor to build up the carbon scaffold, begging the question: which other elements (besides carbon) can be used to tune the surface for optimal methane binding without compromising the structure of the overall material? Boron and nitrogen were both explored and nitrogen was identified as the ideal agent for increasing methane binding to the optimal strength for use at near room temperature. Densification of the resulting materials by mechanical compaction was then explored to produce free-standing pellets with the highest methane “delivery” capacities of any known material. Here we define the amount “delivered” as the amount stored minus the amount remaining in the tank at 73 psi, the lowest useful pressure for vehicular transportation applications. The key features of the resulting methane storage system are rapid discharging and refueling speeds, extremely long cycle life (infinite cycle life for high-purity methane), a near-constant heat of refueling (making the heat generated upon refueling easier to manage), and a low maximum pressure of operation (<1500 psi), enabling the use of all-metal Type I cylinders and significantly reducing the costs of conversion of existing vehicles to cleaner burning, abundantly available natural gas). The technical effectiveness of this approach was further improved by exploring several cost-reduction strategies and investigating the reliability of such storage systems in the presence of “real” natural gas mixtures which contain a mixture of gaseous species.

03 NATURAL GAS↗

Improved Understanding of Zircaloy Hydrogen Pickup Mechanisms in BWRs. NSUF Report (Final Report)

Using an instrumented autoclave operating at 288°C and 1100 psi, electrical impedance spectroscopy (EIS) measurements were taken from three neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). The three high priority fuel channel samples were tested at 288°C using a simulated PWR water with the AC frequency range extending down through 10 -4 Hz, with corresponding DC measurements to validate the resistance measurements. EIS measurements were also taken at room temperature using a simulated PWR water from four neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). and four water rod samples (Zircaloy-2). EIS measurements were also performed on unirradiated Zircaloy-2 tubes that were previously oxidized in autoclaves at various temperatures for different time durations to obtain oxide layers of different thicknesses. These experiments proved to be more challenging than anticipated because of the small size of the samples, which required development of fixtures to allow valid testing of the materials, some of which are radioactive, under pressure and temperature. Equivalent circuits were used to model the data. Some samples resulted in Nyquist plots where the arcs were too close to each other resulting in the inability to separate the different oxide morphological thicknesses. Microstructural characterization in the form of optical microscopy, scanning electron microscopy and (scanning) transmission electron microscopy was performed to document the microstructural features of the oxidized Zr-2 and Zr-4 samples (oxide and metal/oxide interface) and to interpret the EIS results. At the start of this work, it was assumed that despite the laminar cracks and porosity visible in both optical and SEM, the oxide scale is dense near the metal, effectively preventing the water from contacting the metal directly. While it is generally accepted that the large laminar cracks pores allow the water to percolate into the oxide film, it is unknown just how deep it can go, and how interconnected is the porosity, which will dictate how much the porosity is connected electrically to the water within the oxide scale. The EIS results demonstrated that there is in fact a distribution of oxide thicknesses over the various samples, and within each material condition, that are communicating to the liquid. This is derived from the fact that Nyquist plots have an arc radius whose center point is often below the Z’ axis, implying that instead of a single oxide thickness, these oxides have a distribution of thicknesses. This distribution of thicknesses produces arcs that are not clearly separable in most of the samples.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗